Exome sequencing from nanogram amounts of starting DNA: comparing three approaches

Hybridization-based target enrichment protocols require relatively large starting amounts of genomic DNA, which is not always available. Here, we tested three approaches to pre-capture library preparation starting from 10 ng of genomic DNA: (i and ii) whole-genome amplification of DNA samples with REPLI-g (Qiagen) and GenomePlex (Sigma) kits followed by standard library preparation, and (iii) library construction with a low input oriented ThruPLEX kit (Rubicon Genomics). Exome capture with Agilent SureSelectXT2 Human AllExon v4+UTRs capture probes, and HiSeq2000 sequencing were performed for test libraries along with the control library prepared from 1 µg of starting DNA. Tested protocols were characterized in terms of mapping efficiency, enrichment ratio, coverage of the target region, and reliability of SNP genotyping. REPLI-g- and ThruPLEX-FD-based protocols seem to be adequate solutions for exome sequencing of low input samples

Haplotype analysis of APOE intragenic SNPs

BACKGROUND: APOE epsilon4 allele is most common genetic risk factor for Alzheimer\u27s disease (AD) and cognitive decline. However, it remains poorly understood why only some carriers of APOE epsilon4 develop AD and how ethnic variabilities in APOE locus contribute to AD risk. Here, to address the role of APOE haplotypes, we reassessed the diversity of APOE locus in major ethnic groups and in Alzheimer\u27s Disease Neuroimaging Initiative (ADNI) dataset on patients with AD, and subjects with mild cognitive impairment (MCI), and control non-demented individuals. RESULTS: We performed APOE gene haplotype analysis for a short block of five SNPs across the gene using the ADNI whole genome sequencing dataset. The compilation of ADNI data with 1000 Genomes identified the APOE epsilon4 linked haplotypes, which appeared to be distant for the Asian, African and European populations. The common European epsilon4-bearing haplotype is associated with AD but not with MCI, and the Africans lack this haplotype. Haplotypic inference revealed alleles that may confer protection against AD. By assessing the DNA methylation profile of the APOE haplotypes, we found that the AD-associated haplotype features elevated APOE CpG content, implying that this locus can also be regulated by genetic-epigenetic interactions. CONCLUSIONS: We showed that SNP frequency profiles within APOE locus are highly skewed to population-specific haplotypes, suggesting that the ancestral background within different sites at APOE gene may shape the disease phenotype. We propose that our results can be utilized for more specific risk assessment based on population descent of the individuals and on higher specificity of five site haplotypes associated with AD

Changes in the transcriptome of the prefrontal cortex of OXYS rats as the signs of Alzheimer’s disease development

Alzheimer’s disease (AD) is the most prevalent neuro­degenerative disease. It produces atrophic changes in the brain, which cause dementia. The incidence of AD is increasing with increasing life expectancy and gradual aging of the population in developed countries. There are no effective prophylactic inter­ventions because of insufficient understanding of the AD pathogenesis and the absence of adequate experimental models. Recently, we showed that senescence-accelerated OXYS rats represent a promis­ing model of AD; in these rats, accelerated aging of the brain is accompanied by the typical signs of AD: degenerative alterations and death of neurons, a de­crease in synaptic density, mitochondrial dysfunction, hyperphosphorylation of the tau protein, an increased level of amyloid β (Aβ1–42), and the formation of amyloid plaques. To elucidate how these signs develop, we used a nextgeneration RNA sequencing technique (RNA-Seq) to study the prefron­tal-cortex transcriptome of OXYS rats during the manifestation of AD signs (at an age of 5 months) and during their active progres­sion (at an age of 18 months), using age-matched Wistar rats (parental strain) as controls. At the age of 5 months, there were significant differences between OXYS and Wistar rats (p < 0.01) in the mRNA expression of more than 900 genes (> 2000 genes at the age of 18 months) in the prefrontal cortex. Most of these genes were related to neuronal plasticity, protein phosphorylation, Са2+ homeostasis, hypoxia, immune processes, and apoptosis. Between the ages of 5 and 18 months, there were changes in the expression of 499 genes in Wistar rats and changes in the expres­sion of 5500 genes in OXYS rats. Only 333 genes were common between these sets. This finding points to differences in the mechanisms and rates of age-related changes in the brain between normal aging and the period of development of AD-specific neuro­degene­rative processes

Quantifying human genome parameters in aging

Healthy human longevity is a global goal of the world health system. Determining the causes and processes influencing human longevity is the primary fundamental goal facing the scientific community. Currently, the main efforts of the scientific community are aimed at identifying the qualitative characteristics of the genome that determine the trait. At the same time, when evaluating qualitative characteristics, there are many challenges that make it difficult to establish associations. Quantitative traits are burdened with such problems to a lesser extent, but they are largely overlooked in current genomic studies of aging and longevity. Although there is a wide repertoire of quantitative trait analyses based on genomic data, most opportunities are ignored by authors, which, along with the inaccessibility of published data, leads to the loss of this important information. This review focuses on describing quantitative traits important for understanding aging and necessary for analysis in further genomic studies, and recommends the inclusion of the described traits in the analysis. The review considers the relationship between quantitative characteristics of the mitochondrial genome and aging, longevity, and age-related neurodegenerative diseases, such as the frequency of extensive mitochondrial DNA (mtDNA) deletions, mtDNA half-life, the frequency of A>G replacements in the mtDNA heavy chain, the number of mtDNA copies; special attention is paid to the mtDNA methylation sign. A separate section of this review is devoted to the correlation of telomere length parameters with age, as well as the association of telomere length with the amount of mitochondrial DNA. In addition, we consider such a quantitative feature as the rate of accumulation of somatic mutations with aging in relation to the lifespan of living organisms. In general, it may be noted that there are quite serious reasons to suppose that various quantitative characteristics of the genome may be directly or indirectly associated with certain aspects of aging and longevity. At the same time, the available data are clearly insufficient for definitive conclusions and the determination of causal relationships

Mosaic loss of the Y chromosome in human neurodegenerative and oncological diseases

The development of new biomarkers for prediction and early detection of human diseases, as well as for monitoring the response to therapy is one of the most relevant areas of modern human genetics and genomics. Until recently, it was believed that the function of human Y chromosome genes was limited to determining sex and controlling spermatogenesis. Thanks to occurance of large databases of the genome-wide association study (GWAS), there has been a transition to the use of large samples for analyzing genetic changes in both normal and pathological conditions. This has made it possible to assess the association of mosaic aneuploidy of the Y chromosome in somatic cells with a shorter lifespan in men compared to women. Based on data from the UK Biobank, an association was found between mosaic loss of the Y chromosome (mLOY) in peripheral blood leukocytes and the age of men over 70, as well as a number of oncological, cardiac, metabolic, neurodegenerative, and psychiatric diseases. As a result, mLOY in peripheral blood cells has been considered a potential marker of biological age in men and as a marker of certain age-related diseases. Currently, numerous associations have been identified between mLOY and genes based on GWAS and transcriptomes in affected tissues. However, the exact cause of mLOY and the impact and consequences of this phenomenon at the whole organism level have not been established. In particular, it is unclear whether aneuploidy of the Y chromosome in blood cells may affect the development of pathologies that manifest in other organs, such as the brain in Alzheimer’s disease, or whether it is a neutral biomarker of general genomic instability. This review examines the main pathologies and genetic factors associated with mLOY, as well as the hypotheses regarding their interplay. Special attention is given to recent studies on mLOY in brain cells in Alzheimer’s disease

Chromatin profiling of cortical neurons identifies individual epigenetic signatures in schizophrenia

Both heritability and environment contribute to risk for schizophrenia. However, the molecular mechanisms of interactions between genetic and non-genetic factors remain unclear. Epigenetic regulation of neuronal genome may be a presumable mechanism in pathogenesis of schizophrenia. Here, we performed analysis of open chromatin landscape of gene promoters in prefrontal cortical (PFC) neurons from schizophrenic patients. We cataloged cell-type-based epigenetic signals of transcriptional start sites (TSS) marked by histone H3-K4 trimethylation (H3K4me3) across the genome in PFC from multiple schizophrenia subjects and age-matched control individuals. One of the top-ranked chromatin alterations was found in the major histocompatibility (MHC) locus on chromosome 6 highlighting the overlap between genetic and epigenetic risk factors in schizophrenia. The chromosome conformation capture (3C) analysis in human brain cells revealed the architecture of multipoint chromatin interactions between the schizophrenia-associated genetic and epigenetic polymorphic sites and distantly located HLA-DRB5 and BTNL2 genes. In addition, schizophrenia-specific chromatin modifications in neurons were particularly prominent for non-coding RNA genes, including an uncharacterized LINC01115 gene and recently identified BNRNA_052780. Notably, protein-coding genes with altered epigenetic state in schizophrenia are enriched for oxidative stress and cell motility pathways. Our results imply the rare individual epigenetic alterations in brain neurons are involved in the pathogenesis of schizophrenia

ANALYSIS OF THE INTERACTIONS OF NEURONAL APOPTOSIS GENES IN THE ASSOCIATIVE GENE NETWORK OF PARKINSON’S DISEASE

Parkinson’s disease (PD) affects an estimated 7–10 million people worldwide and 210000 people in Russia. PD is accompanied by degeneration of dopaminergic neurons and because of that neuronal apoptosis is an important factor in this disease. Analysis of gene networks is one of the key approaches in systems biology. We previously developed the ANDSystem tool, designed to automatically extract knowledge from scientific publications and reconstruct on this basis associative gene networks describing the molecular genetic mechanisms of biological processes. The aim of this work was prioritization of neuronal apoptosis genes by their involvement in PD pathogenesis, taking into account the structure of the PD associative gene network using ANDSystem. Analysis of the centrality of neuronal apoptosis genes, associated with PD, revealed that mean values of degree, closeness and betweenness centralities statistically significantly exceed such values of all nodes of the PD network. The APOE, CASP3 and GAPDH genes involved in neuronal apoptosis were among the most central genes. Prioritization of neuronal apoptosis genes for which there was no data in ANDSystem on their associations with PD was performed using standard methods (Endeavor and ToppGene) and the criteria of centrality and specificity of genes interactions with the PD gene network. Analysis revealed that genes involved in such processes as positive and negative regulation of neu ronal apoptosis, MAPK and ephrin receptor signaling pathways, are mainly represented among candidate genes with the highest priority (top 50, 70, 100 genes were considered). In particular, TP53, JUN, BCL2, PIK3CA and APP were among candidate genes with the highest priority

Whole exome sequencing links dental tumor to an autosomal-dominant mutation in ANO5 gene associated with gnathodiaphyseal dysplasia and muscle dystrophies

Tumors of the jaws may represent different human disorders and frequently associate with pathologic bone fractures. In this report, we analyzed two affected siblings from a family of Russian origin, with a history of dental tumors of the jaws, in correspondence to original clinical diagnosis of cementoma consistent with gigantiform cementoma (GC, OMIM: 137575). Whole exome sequencing revealed the heterozygous missense mutation c.1067G \u3e A (p.Cys356Tyr) in ANO5 gene in these patients. To date, autosomal-dominant mutations have been described in the ANO5 gene for gnathodiaphyseal dysplasia (GDD, OMIM: 166260), and multiple recessive mutations have been described in the gene for muscle dystrophies (OMIM: 613319, 611307); the same amino acid (Cys) at the position 356 is mutated in GDD. These genetic data and similar clinical phenotypes demonstrate that the GC and GDD likely represent the same type of bone pathology. Our data illustrate the significance of mutations in single amino-acid position for particular bone tissue pathology. Modifying role of genetic variations in another gene on the severity of the monogenic trait pathology is also suggested. Finally, we propose the model explaining the tissue-specific manifestation of clinically distant bone and muscle diseases linked to mutations in one gene

The influence of the LED lighting on structural-functional parameters of lettuce plants

Comparative studies were conducted on morphofunctional parameters of lettuce (Lactuca sativa) plants cultivated with the use of supplementary lighting from different light sources - LED lamps, high-pressure mercury arc lamp (HPM) and high-pressure sodium arc lamp (HPS). The work showed that morphofunctional parameters of lettuce were significantly higher with LED lighting, in comparison to plants cultivated with HPM illumination, and main parameters of leaves grown under LED were barely discernible from the samples grown under HPS lighting. We observed accelerated transfer into reproductive period during HPM lighting, what can lead to quality degradation of the product. Considering the approximate values of leaf character range in experiments with HPS and LED lighting and greater energy efficiency of LED, light emitting diodes are obviously more promising for supplementary illumination in protected ground conditions

Role of apoptosis genes in aggression revealed using combined analysis of ANDSystem gene networks, expression and genomic data in grey rats with aggressive behavior

Aggressive behavior in animals plays an important role in protecting the territory, offspring, establishing social hierarchical relations, etc. Increased aggression is observed in a number of diseases ( schizophrenia, bipolar disorder, brain degenerative disorders). Neuronal apoptosis is crucial in the maintenance of developmental processes during neurogenesis. Alterations in neuronal apoptosis are observed in aging and neuropathologies accompanied by changes in psycho­emo­ tional state (epilepsy, Alzheimer’s disease, neurotrauma). The expression of key neuronal apoptosis genes (Casp3, Bax and Bcl-xl) in the brain of highly aggressive rats is significantly altered. The aim of this work was to analyze associative networks that describe genetic interactions between genes/proteins involved in neuronal apoptosis, differentially expressed genes and genes with polymorphisms in grey rats with aggressive behavior. Analysis revealed 819 differentially expressed genes in the hypothalamus, ventral tegmental region and periaqueductus Sylvii grey matter in grey rats with aggressive and tame behavior. The Stx1a, Mbp and Th genes have the highest index of betweenness centrality in the associative network of differentially expressed genes. Genome analysis revealed 137 polymorphic genes. Three of them (Lig4, Parp1 and Pigt) were involved in neuronal apoptosis. It was shown that polymorphic and differentially expressed genes were statistically significantly overrepresented among ge nes interacting with neuronal apoptosis genes (p value < 0.01). Three molecular­genetic chains describing connections between polymorphic and neuronal apoptosis genes mediated by differentially expressed genes were reconstructed. Chains included the polymorphic genes Tsc1, Adamts4 and Lgals3, differentially expressed genes Ezr, Acan, Th and 19 neuronal apoptosis genes. It was shown that neuronal apoptosis is closely related to aggressive behavior in animals