Tissue-specific effects of genetic and epigenetic variation on gene regulation and splicing

Expression of Truncated Products at the 5'-Terminal Region of RIPK2 and Evolutive Aspects that Support Their Biological Importance

Alternative splicing is the process of generating different mRNAs from the same primary transcript, which contributes to increase the transcriptome and proteome diversity. Abnormal splicing has been associated with the development of several diseases including cancer. Given that mutations and abnormal levels of the RIPK2 transcript and RIP-2 protein are frequent in tumors, and that RIP-2 modulates immune and inflammatory responses, we investigated alternative splicing events that result in partial deletions of the kinase domain at the N-terminus of RIP-2. We also investigated the structure and expression of the RIPK2 truncated variants and isoforms in different environments. In addition, we searched data throughout Supraprimates evolution that could support the biological importance of RIPK2 alternatively spliced products. We observed that human variants and isoforms were differentially regulated following temperature stress, and that the truncated transcript was more expressed than the long transcript in tumor samples. The inverse was found for the longer protein isoform. The truncated variant was also detected in chimpanzee, gorilla, hare, pika, mouse, rat, and tree shrew. The fact that the same variant has been preserved in mammals with divergence times up to 70 million years raises the hypothesis that it may have a functional significance.

History of tuberculosis disease is associated with genetic regulatory variation in Peruvians

A quarter of humanity is estimated to have been exposed to Mycobacterium tuberculosis (Mtb) with a 5-10% risk of developing tuberculosis (TB) disease. Variability in responses to Mtb infection could be due to host or pathogen heterogeneity. Here, we focused on host genetic variation in a Peruvian population and its associations with gene regulation in monocyte-derived macrophages and dendritic cells (DCs). We recruited former household contacts of TB patients who previously progressed to TB (cases, n = 63) or did not progress to TB (controls, n = 63). Transcriptomic profiling of monocyte-derived DCs and macrophages measured the impact of genetic variants on gene expression by identifying expression quantitative trait loci (eQTL). We identified 330 and 257 eQTL genes in DCs and macrophages (False Discovery Rate (FDR) < 0.05), respectively. Four genes in DCs showed interaction between eQTL variants and TB progression status. The top eQTL interaction for a protein-coding gene was with FAH, the gene encoding fumarylacetoacetate hydrolase, which mediates the last step in mammalian tyrosine catabolism. FAH expression was associated with genetic regulatory variation in cases but not controls. Using public transcriptomic and epigenomic data of Mtb-infected monocyte-derived dendritic cells, we found that Mtb infection results in FAH downregulation and DNA methylation changes in the locus. Overall, this study demonstrates effects of genetic variation on gene expression levels that are dependent on history of infectious disease and highlights a candidate pathogenic mechanism through pathogen-response genes. Furthermore, our results point to tyrosine metabolism and related candidate TB progression pathways for further investigation.

Alterations in Blood Methylome as Potential Epigenetic Biomarker in Sporadic Parkinson's Disease

OBJECTIVE

To characterize DNA methylation (DNAm) differences between sporadic Parkinson's disease (PD) and healthy control (HC) individuals enrolled in the Parkinson's Progression Markers Initiative (PPMI).

METHODS

Using whole blood, we characterized longitudinal differences in DNAm between sporadic PD patients (n = 196) and HCs (n = 86) enrolled in PPMI. RNA sequencing (RNAseq) was used to conduct gene expression analyses for genes mapped to differentially methylated cytosine-guanine sites (CpGs).

RESULTS

At the time of patient enrollment, 5,178 CpGs were differentially methylated (2,683 hypermethylated and 2,495 hypomethylated) in PD compared to HC. Of these, 579 CpGs underwent significant methylation changes over 3 years. Several differentially methylated CpGs were found near the cytochrome P450 family 2 subfamily E member 1 (CYP2E1) gene. Additionally, multiple hypermethylated CpGs were associated with the N-myc downregulated gene family member 4 (NDRG4) gene. RNA-Seq analyses showed 75 differentially expressed genes in PD patients compared to controls. An integrative analysis of both differentially methylated sites and differentially expressed genes revealed 20 genes that exhibited hypomethylation concomitant with overexpression. Additionally, 1 gene, cathepsin H (CTSH), displayed hypermethylation that was associated with its decreased expression.

INTERPRETATION

We provide initial evidence of alterations in DNAm in blood of PD patients that may serve as potential epigenetic biomarker of disease. To evaluate the significance of these changes throughout the progression of PD, additional profiling at longer intervals and during the prodromal stages of disease will be necessary. ANN NEUROL 2024;95:1162-1172.

Functional genomics implicates natural killer cells in the pathogenesis of ankylosing spondylitis

Objective

Multiple lines of evidence indicate that ankylosing spondylitis (AS) is a lymphocyte-driven disease. However, which lymphocyte populations are critical in AS pathogenesis is not known. In this study, we aimed to identify the key cell types mediating the genetic risk in AS using an unbiased functional genomics approach.

Methods

We integrated genome-wide association study (GWAS) data with epigenomic and transcriptomic datasets of human immune cells. To quantify enrichment of cell type-specific open chromatin or gene expression in AS risk loci, we used three published methods that have successfully identified relevant cell types in other diseases. We performed co-localization analyses between GWAS risk loci and genetic variants associated with gene expression (eQTL) to find putative target genes.

Results

Natural killer (NK) cell-specific open chromatin regions are significantly enriched in heritability for AS, compared to other immune cell types such as T cells, B cells, and monocytes. This finding was consistent between two AS GWAS. Using RNA-seq data, we validated that genes in AS risk loci are enriched in NK cell-specific gene expression. Using the human Space-Time Gut Cell Atlas, we also found significant upregulation of AS-associated genes predominantly in NK cells. Co-localization analysis revealed four AS risk loci affecting regulation of candidate target genes in NK cells: two known loci, ERAP1 and TNFRSF1A, and two under-studied loci, ENTR1 (aka SDCCAG3) and B3GNT2.

Conclusion

Our findings suggest that NK cells may play a crucial role in AS development and highlight four putative target genes for functional follow-up in NK cells.

Former smoking associated with epigenetic modifications in human granulosa cells among women undergoing assisted reproduction

Smoking exposure during adulthood can disrupt oocyte development in women, contributing to infertility and possibly adverse birth outcomes. Some of these effects may be reflected in epigenome profiles in granulosa cells (GCs) in human follicular fluid. We compared the epigenetic modifications throughout the genome in GCs from women who were former (N = 15) versus never smokers (N = 44) undergoing assisted reproductive technologies (ART). This study included 59 women undergoing ART. Smoking history including time since quitting was determined by questionnaire. GCs were collected during oocyte retrieval and DNA methylation (DNAm) levels were profiled using the Infinium MethylationEPIC BeadChip. We performed an epigenome-wide association study with robust linear models, regressing DNAm level at individual loci on smoking status, adjusting for age, ovarian stimulation protocol, and three surrogate variables. We performed differentially methylated regions (DMRs) analysis and over-representation analysis of the identified CpGs and corresponding gene set. 81 CpGs were differentially methylated among former smokers compared to never smokers (FDR < 0.05). We identified 2 significant DMRs (KCNQ1 and RHBDD2). The former smoking-associated genes were enriched in oxytocin signaling, adrenergic signaling in cardiomyocytes, platelet activation, axon guidance, and chemokine signaling pathway. These epigenetic variations have been associated with inflammatory responses, reproductive outcomes, cancer development, neurodevelopmental disorder, and cardiometabolic health. Secondarily, we examined the relationships between time since quitting and DNAm at significant CpGs. We observed three CpGs in negative associations with the length of quitting smoking (p < 0.05), which were cg04254052 (KCNIP1), cg22875371 (OGDHL), and cg27289628 (LOC148145), while one in positive association, which was cg13487862 (PLXNB1). As a pilot study, we demonstrated epigenetic modifications associated with former smoking in GCs. The study is informative to potential biological pathways underlying the documented association between smoking and female infertility and biomarker discovery for smoking-associated reproductive outcomes.

Identification of region-specific gene isoforms in the human brain using long-read transcriptome sequencing

In neurological and neuropsychiatric diseases, different brain regions are affected, and differences in gene expression patterns could potentially explain this mechanism. However, limited studies have precisely explored gene expression in different regions of the human brain. In this study, we performed long-read RNA sequencing on three different brain regions of the same individuals: the cerebellum, hypothalamus, and temporal cortex. Despite stringent filtering criteria excluding isoforms predicted to be artifacts, over half of the isoforms expressed in multiple samples across multiple regions were found to be unregistered in the GENCODE reference. We then especially focused on genes with different major isoforms in each brain region, even with similar overall expression levels, and identified that many of such genes including GAS7 might have distinct roles in dendritic spine and neuronal formation in each region. We also found that DNA methylation might, in part, drive different isoform expressions in different regions. These findings highlight the significance of analyzing isoforms expressed in disease-relevant sites.

DNA methylome analysis provides insights into gene regulatory mechanism for better performance of rice under fluctuating environmental conditions: epigenomics of adaptive plasticity

MAIN CONCLUSION

Roots play an important role in adaptive plasticity of rice under dry/direct-sown conditions. However, hypomethylation of genes in leaves (resulting in up-regulated expression) complements the adaptive plasticity of Nagina-22 under DSR conditions. Rice is generally cultivated by transplanting which requires plenty of water for irrigation. Such a practice makes rice cultivation a challenging task under global climate change and reducing water availability. However, dry-seeded/direct-sown rice (DSR) has emerged as a resource-saving alternative to transplanted rice (TPR). Though some of the well-adapted local cultivars are used for DSR, only limited success has been achieved in developing DSR varieties mainly because of a limited knowledge of adaptability of rice under fluctuating environmental conditions. Based on better morpho-physiological and agronomic performance of Nagina-22 (N-22) under DSR conditions, N-22 and IR-64 were grown by transplanting and direct-sowing and used for whole genome methylome analysis to unravel the epigenetic basis of adaptive plasticity of rice. Comparative methylome and transcriptome analyses indicated a large number (4078) of genes regulated through DNA methylation/demethylation in N-22 under DSR conditions. Gene × environment interactions play important roles in adaptive plasticity of rice under direct-sown conditions. While genes for pectinesterase, LRK10, C2H2 zinc-finger protein, splicing factor, transposable elements, and some of the unannotated proteins were hypermethylated, the genes for regulation of transcription, protein phosphorylation, etc. were hypomethylated in CG context in the root of N-22, which played important roles in providing adaptive plasticity to N-22 under DSR conditions. Hypomethylation leading to up-regulation of gene expression in the leaf complements the adaptive plasticity of N-22 under DSR conditions. Moreover, differential post-translational modification of proteins and chromatin assembly/disassembly through DNA methylation in CHG context modulate adaptive plasticity of N-22. These findings would help developing DSR cultivars for increased water-productivity and ecological efficiency.

Exploring multisite heterogeneity of human basal cell carcinoma proteome and transcriptome

Basal cell carcinoma (BCC) is the most common type of skin cancer. Due to multiple, potential underlying molecular tumor aberrations, clinical treatment protocols are not well-defined. This study presents multisite molecular heterogeneity profiles of human BCC based on RNA and proteome profiling. Three areas from lesions excised from 9 patients were analyzed. The focus was gene expression profiles based on proteome and RNA measurements of intra-tumor heterogeneity from the same patient and inter-tumor heterogeneity in nodular, infiltrative, and superficial BCC tumor subtypes from different patients. We observed significant overlap in intra- and inter-tumor variability of proteome and RNA expression profiles, showing significant multisite heterogeneity of protein expression in the BCC tumors. Inter-subtype analysis has also identified unique proteins for each BCC subtype. This profiling leads to a deeper understanding of BCC molecular heterogeneity and potentially contributes to developing new sampling tools for personalized diagnostics therapeutic approaches to BCC.

Identification of a regulatory pathway governing TRAF1 via an arthritis-associated non-coding variant

TRAF1/C5 was among the first loci shown to confer risk for inflammatory arthritis in the absence of an associated coding variant, but its genetic mechanism remains undefined. Using Immunochip data from 3,939 patients with juvenile idiopathic arthritis (JIA) and 14,412 control individuals, we identified 132 plausible common non-coding variants, reduced serially by single-nucleotide polymorphism sequencing (SNP-seq), electrophoretic mobility shift, and luciferase studies to the single variant rs7034653 in the third intron of TRAF1. Genetically manipulated experimental cells and primary monocytes from genotyped donors establish that the risk G allele reduces binding of Fos-related antigen 2 (FRA2), encoded by FOSL2, resulting in reduced TRAF1 expression and enhanced tumor necrosis factor (TNF) production. Conditioning on this JIA variant eliminated attributable risk for rheumatoid arthritis, implicating a mechanism shared across the arthritis spectrum. These findings reveal that rs7034653, FRA2, and TRAF1 mediate a pathway through which a non-coding functional variant drives risk of inflammatory arthritis in children and adults.

Estrogen administration and withdrawal in a model of hormone-simulated pregnancy lead to alterations in behavior and gene expression but do not induce depression-like phenotypes in mice

Pregnancy and the post-partum period are associated with substantial fluctuations in hormone levels and are frequently associated with significant stress. Many individuals also experience affective disturbances during the peri‑partum period, including anxiety, the 'baby blues,' and post-partum depression. However, the extent to which these affective changes result from rapidly altering hormone levels, increased stress, or the combination of both remains largely unknown. The current study sought to evaluate the consequences of pregnancy-like hormonal changes on behavior and gene expression in c57BL/6 mice in the absence of stress using a hormone-simulated pregnancy model. Our results reveal that animals receiving hormone injections to simulate the high levels of estrogen observed in late pregnancy and animals withdrawn from estrogen to mimic the rapid decline in this hormone following parturition both exhibit increased anxiety-like behavior compared to ovariectomized controls in the novel open field test. However, no other significant anxiety- or depression-like alterations were observed in either hormone-treated group compared to ovariectomized controls. Both hormone administration and estrogen withdrawal were shown to induce several significant alterations in gene expression in the bed nucleus of the stria terminalis and the paraventricular nucleus of the hypothalamus. In contrast to the estrogen withdrawal hypothesis of post-partum depression, our results suggest that this method estrogen withdrawal following hormone-simulated pregnancy in the absence of stress does not induce phenotypes consistent with post-partum depression in c57BL/6 mice. However, given that estrogen withdrawal does lead to significant gene expression changes in two stress-sensitive brain regions, it remains possible that estrogen withdrawal could still contribute to affective dysregulation in the peri-partum period by influencing susceptibility to stress. Future research is required to evaluate this possibility.

History of tuberculosis disease is associated with genetic regulatory variation in Peruvians

A quarter of humanity is estimated to be latently infected with Mycobacterium tuberculosis (Mtb) with a 5-10% risk of developing tuberculosis (TB) disease. Variability in responses to Mtb infection could be due to host or pathogen heterogeneity. Here, we focused on host genetic variation in a Peruvian population and its associations with gene regulation in monocyte-derived macrophages and dendritic cells (DCs). We recruited former household contacts of TB patients who previously progressed to TB (cases, n=63) or did not progress to TB (controls, n=63). Transcriptomic profiling of monocyte-derived dendritic cells (DCs) and macrophages measured the impact of genetic variants on gene expression by identifying expression quantitative trait loci (eQTL). We identified 330 and 257 eQTL genes in DCs and macrophages (False Discovery Rate (FDR) < 0.05), respectively. Five genes in DCs showed interaction between eQTL variants and TB progression status. The top eQTL interaction for a protein-coding gene was with FAH, the gene encoding fumarylacetoacetate hydrolase, which mediates the last step in mammalian tyrosine catabolism. FAH expression was associated with genetic regulatory variation in cases but not controls. Using public transcriptomic and epigenomic data of Mtb-infected monocyte-derived dendritic cells, we found that Mtb infection results in FAH downregulation and DNA methylation changes in the locus. Overall, this study demonstrates effects of genetic variation on gene expression levels that are dependent on history of infectious disease and highlights a candidate pathogenic mechanism through pathogen-response genes. Furthermore, our results point to tyrosine metabolism and related candidate TB progression pathways for further investigation.

Knockdown of sirtuin6 positively regulates acetylation of DNMT1 to inhibit NOTCH signaling pathway in non-small cell lung cancer cell lines

BACKGROUND AND AIM

Sirtuin proteins (1-7) are nicotinamide adenine dinucleotide (NAD)-dependent deacetylases and ADP-ribosyl transferases (class III histone deacetylase enzymes (HDAC)) mainly involved in the removal of the acetyl group from histone proteins. SIRT6, one of the sirtuins, plays a major role in cancer progression in many types of cancer conditions. We recently reported that SIRT6 acts as an oncogene in NSCLC; thus, silencing of SIRT6 inhibits cell proliferation and induces apoptosis in NSCLC cell lines. NOTCH signaling has been reported to be involved in cell survival and regulates cell proliferation and differentiation. However, recent studies from different groups have converged on the notion that NOTCH1 may be an important oncogene in NSCLC. The abnormal expression of NOTCH signaling pathway members is a relatively frequent event in patients with NSCLC. SIRT6 and the NOTCH signaling pathway might play a critical role in tumorigenesis since they are highly expressed in NSCLC. This study has been performed to explore the exact mechanism by which SIRT6 inhibits cell proliferation and induces the apoptosis of NSCLC cell lines and its correlation with NOTCH signaling.

MAIN METHODS

In vitro experiments with human NSCLC cells have been performed. Immunocytochemistry study was used to analyze the expression of NOTCH1 and DNMT1 in A549 and NCI-H460 cell lines. RT-qPCR, Western Blot, Methylated DNA specific PCR, and Co-Immunoprecipitation were performed to explore the key events in the regulation of NOTCH signaling by silencing SIRT6 in NSCLC cell lines.

KEY FINDINGS

The findings of this study suggest that silencing of SIRT6 significantly promotes the acetylation status of DNMT1 and stabilizes it. Consequently, acetylated DNMT1 translocates into the nucleus and methylates the NOTCH1 promoter region, resulting in the hindering of NOTCH1-mediated NOTCH signaling.

Distinct genome-wide DNA methylation and gene expression signatures in classical monocytes from African American patients with systemic sclerosis

BACKGROUND

Systemic sclerosis (SSc) is a multisystem autoimmune disorder that has an unclear etiology and disproportionately affects women and African Americans. Despite this, African Americans are dramatically underrepresented in SSc research. Additionally, monocytes show heightened activation in SSc and in African Americans relative to European Americans. In this study, we sought to investigate DNA methylation and gene expression patterns in classical monocytes in a health disparity population.

METHODS

Classical monocytes (CD14+ + CD16-) were FACS-isolated from 34 self-reported African American women. Samples from 12 SSc patients and 12 healthy controls were hybridized on MethylationEPIC BeadChip array, while RNA-seq was performed on 16 SSc patients and 18 healthy controls. Analyses were computed to identify differentially methylated CpGs (DMCs), differentially expressed genes (DEGs), and CpGs associated with changes in gene expression (eQTM analysis).

RESULTS

We observed modest DNA methylation and gene expression differences between cases and controls. The genes harboring the top DMCs, the top DEGs, as well as the top eQTM loci were enriched for metabolic processes. Genes involved in immune processes and pathways showed a weak upregulation in the transcriptomic analysis. While many genes were newly identified, several other have been previously reported as differentially methylated or expressed in different blood cells from patients with SSc, supporting for their potential dysregulation in SSc.

CONCLUSIONS

While contrasting with results found in other blood cell types in largely European-descent groups, the results of this study support that variation in DNA methylation and gene expression exists among different cell types and individuals of different genetic, clinical, social, and environmental backgrounds. This finding supports the importance of including diverse, well-characterized patients to understand the different roles of DNA methylation and gene expression variability in the dysregulation of classical monocytes in diverse populations, which might help explaining the health disparities.

DNA methylation QTL mapping across diverse human tissues provides molecular links between genetic variation and complex traits

Studies of DNA methylation (DNAm) in solid human tissues are relatively scarce; tissue-specific characterization of DNAm is needed to understand its role in gene regulation and its relevance to complex traits. We generated array-based DNAm profiles for 987 human samples from the Genotype-Tissue Expression (GTEx) project, representing 9 tissue types and 424 subjects. We characterized methylome and transcriptome correlations (eQTMs), genetic regulation in cis (mQTLs and eQTLs) across tissues and e/mQTLs links to complex traits. We identified mQTLs for 286,152 CpG sites, many of which (>5%) show tissue specificity, and mQTL colocalizations with 2,254 distinct GWAS hits across 83 traits. For 91% of these loci, a candidate gene link was identified by integration of functional maps, including eQTMs, and/or eQTL colocalization, but only 33% of loci involved an eQTL and mQTL present in the same tissue type. With this DNAm-focused integrative analysis, we contribute to the understanding of molecular regulatory mechanisms in human tissues and their impact on complex traits.

Maternal and Placental DNA Methylation Changes Associated with the Pathogenesis of Gestational Diabetes Mellitus

Gestational diabetes mellitus (GDM) is an important metabolic complication of pregnancy, which affects the future health of both the mother and the newborn. The pathogenesis of GDM is not completely clear, but what is clear is that with the development and growth of the placenta, GDM onset and blood glucose is difficult to control, while gestational diabetes patients' blood glucose drops and reaches normal after placenta delivery. This may be associated with placental secretion of insulin-like growth factor, adipokines, tumor necrosis factor-α, cytokines and insulin resistance. Therefore, endocrine secretion of placenta plays a key role in the pathogenesis of GDM. The influence of DNA methylation of these molecules and pathway-related genes on gene expression is also closely related to the pathogenesis of GDM. Here, this review attempts to clarify the pathogenesis of GDM and the related maternal and placental DNA methylation changes and how they affect metabolic pathways.

Spatiotemporal specificity of correlated DNA methylation and gene expression pairs across different human tissues and stages of brain development

DNA methylation (DNAm) that occurs on promoter regions is primarily considered to repress gene expression. Previous studies indicated that DNAm could also show positive correlations with gene expression. Both DNAm and gene expression profiles are known to be tissue- and development-specific. This study aims to investigate how DNAm and gene expression are coordinated across different human tissues and developmental stages, as well as the biological significance of such correlations. By analyzing 2,239 samples with both DNAm and gene expression data in the same human subjects obtained from six published datasets, we evaluated the correlations between gene and CpG pairs (GCPs) at cis-regions and compared significantly correlated GCPs (cGCPs) across different tissues and brains at different age groups. A total of 37,363 cGCPs was identified in the six datasets; approximately 38% of the cGCPs were positively correlated. The majority (>90%) of cGCPs was tissue- or development-specific. We also observed that the correlation direction can be opposite in different tissues and ages. Further analysis highlights the importance of cGCPs for their cellular functions and potential roles in complex traits and human diseases. For instance, the early developmental brain possessed a highly unique set of cGCPs that were associated with neurogenesis and psychiatric disorders. By assessing the epigenetic factors involved in cGCPs, we discovered novel regulatory mechanisms of positive cGCPs distinct from negative cGCPs, which were related to multiple factors, such as H3K27me3, CTCF, and JARD2. The catalogue of cGCPs compiled can be used to guide functional interpretation of genetic and epigenetic studies.

Epigenetic signatures relating to disease-associated genotypic burden in familial risk of bipolar disorder

Environmental factors contribute to risk of bipolar disorder (BD), but how environmental factors impact the development of psychopathology within the context of elevated genetic risk is unknown. We herein sought to identify epigenetic signatures operating in the context of polygenic risk for BD in young people at high familial risk (HR) of BD. Peripheral blood-derived DNA was assayed using Illumina PsychArray, and Methylation-450K or -EPIC BeadChips. Polygenic risk scores (PRS) were calculated using summary statistics from recent genome-wide association studies for BD, major depressive disorder (MDD) and cross-disorder (meta-analysis of eight psychiatric disorders). Unrelated HR participants of European ancestry (n = 103) were stratified based on their BD-PRS score within the HR-population distribution, and the top two quintiles (High-BD-PRS; n = 41) compared against the bottom two quintiles (Low-BD-PRS; n = 41). The High-BD-PRS stratum also had higher mean cross-disorder-PRS and MDD-PRS (ANCOVA p = 0.035 and p = 0.024, respectively). We evaluated DNA methylation differences between High-BD-PRS and Low-BD-PRS strata using linear models. One differentially methylated probe (DMP) (cg00933603; p = 3.54 × 10^-7) in VARS2, a mitochondrial aminoacyl-tRNA synthetase, remained significantly hypomethylated after multiple-testing correction. Overall, BD-PRS appeared to broadly impact epigenetic processes, with 1,183 genes mapped to nominal DMPs (p < 0.05); these displayed convergence with genes previously associated with BD, schizophrenia, chronotype, and risk taking. We tested poly-methylomic epigenetic profiles derived from nominal DMPs in two independent samples (n = 54 and n = 82, respectively), and conducted an exploratory evaluation of the effects of family environment, indexing cohesion and flexibility. This study highlights an important interplay between heritable risk and epigenetic factors, which warrant further exploration.

Challenges in Analyzing Functional Epigenetic Data in Perspective of Adolescent Psychiatric Health

The formative period of adolescence plays a crucial role in the development of skills and abilities for adulthood. Adolescents who are affected by mental health conditions are at risk of suicide and social and academic impairments. Gene–environment complementary contributions to the molecular mechanisms involved in psychiatric disorders have emphasized the need to analyze epigenetic marks such as DNA methylation (DNAm) and non-coding RNAs. However, the large and diverse bioinformatic and statistical methods, referring to the confounders of the statistical models, application of multiple-testing adjustment methods, questions regarding the correlation of DNAm across tissues, and sex-dependent differences in results, have raised challenges regarding the interpretation of the results. Based on the example of generalized anxiety disorder (GAD) and depressive disorder (MDD), we shed light on the current knowledge and usage of methodological tools in analyzing epigenetics. Statistical robustness is an essential prerequisite for a better understanding and interpretation of epigenetic modifications and helps to find novel targets for personalized therapeutics in psychiatric diseases.

Functional Characterization of Genetic Variant Effects on Expression

Thousands of common genetic variants in the human population have been associated with disease risk and phenotypic variation by genome-wide association studies (GWAS). However, the majority of GWAS variants fall into noncoding regions of the genome, complicating our understanding of their regulatory functions, and few molecular mechanisms of GWAS variant effects have been clearly elucidated. Here, we set out to review genetic variant effects, focusing on expression quantitative trait loci (eQTLs), including their utility in interpreting GWAS variant mechanisms. We discuss the interrelated challenges and opportunities for eQTL analysis, covering determining causal variants, elucidating molecular mechanisms of action, and understanding context variability. Addressing these questions can enable better functional characterization of disease-associated loci and provide insights into fundamental biological questions of the noncoding genetic regulatory code and its control of gene expression. Expected final online publication date for the Annual Review of Biomedical Data Science, Volume 5 is August 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Epigenetics in Ecology, Evolution, and Conservation

Epigenetic variation is often characterized by modifications to DNA that do not alter the underlying nucleotide sequence, but can influence behavior, morphology, and physiological phenotypes by affecting gene expression and protein synthesis. In this review, we consider how the emerging field of ecological epigenetics (eco-epi) aims to use epigenetic variation to explain ecologically relevant phenotypic variation and predict evolutionary trajectories that are important in conservation. Here, we focus on how epigenetic data have contributed to our understanding of wild populations, including plants, animals, and fungi. First, we identified published eco-epi literature and found that there was limited taxonomic and ecosystem coverage and that, by necessity of available technology, these studies have most often focused on the summarized epigenome rather than locus- or nucleotide-level epigenome characteristics. We also found that while many studies focused on adaptation and heritability of the epigenome, the field has thematically expanded into topics such as disease ecology and epigenome-based ageing of individuals. In the second part of our synthesis, we discuss key insights that have emerged from the epigenetic field broadly and use these to preview the path toward integration of epigenetics into ecology. Specifically, we suggest moving focus to nucleotide-level differences in the epigenome rather than whole-epigenome data and that we incorporate several facets of epigenome characterization (e.g., methylation, chromatin structure). Finally, we also suggest that incorporation of behavior and stress data will be critical to the process of fully integrating eco-epi data into ecology, conservation, and evolutionary biology.

Identification of autosomal cis expression quantitative trait methylation (cis eQTMs) in children's blood

Background

The identification of expression quantitative trait methylation (eQTMs), defined as associations between DNA methylation levels and gene expression, might help the biological interpretation of epigenome-wide association studies (EWAS). We aimed to identify autosomal cis eQTMs in children's blood, using data from 832 children of the Human Early Life Exposome (HELIX) project.

Methods

Blood DNA methylation and gene expression were measured with the Illumina 450K and the Affymetrix HTA v2 arrays, respectively. The relationship between methylation levels and expression of nearby genes (1 Mb window centered at the transcription start site, TSS) was assessed by fitting 13.6 M linear regressions adjusting for sex, age, cohort, and blood cell composition.

Results

We identified 39,749 blood autosomal cis eQTMs, representing 21,966 unique CpGs (eCpGs, 5.7% of total CpGs) and 8,886 unique transcript clusters (eGenes, 15.3% of total transcript clusters, equivalent to genes). In 87.9% of these cis eQTMs, the eCpG was located at <250 kb from eGene's TSS; and 58.8% of all eQTMs showed an inverse relationship between the methylation and expression levels. Only around half of the autosomal cis-eQTMs eGenes could be captured through annotation of the eCpG to the closest gene. eCpGs had less measurement error and were enriched for active blood regulatory regions and for CpGs reported to be associated with environmental exposures or phenotypic traits. In 40.4% of the eQTMs, the CpG and the eGene were both associated with at least one genetic variant. The overlap of autosomal cis eQTMs in children's blood with those described in adults was small (13.8%), and age-shared cis eQTMs tended to be proximal to the TSS and enriched for genetic variants.

Conclusions

This catalogue of autosomal cis eQTMs in children's blood can help the biological interpretation of EWAS findings and is publicly available at https://helixomics.isglobal.org/ and at Dryad (doi:10.5061/dryad.fxpnvx0t0).

Funding

The study has received funding from the European Community's Seventh Framework Programme (FP7/2007-206) under grant agreement no 308333 (HELIX project); the H2020-EU.3.1.2. - Preventing Disease Programme under grant agreement no 874583 (ATHLETE project); from the European Union's Horizon 2020 research and innovation programme under grant agreement no 733206 (LIFECYCLE project), and from the European Joint Programming Initiative "A Healthy Diet for a Healthy Life" (JPI HDHL and Instituto de Salud Carlos III) under the grant agreement no AC18/00006 (NutriPROGRAM project). The genotyping was supported by the projects PI17/01225 and PI17/01935, funded by the Instituto de Salud Carlos III and co-funded by European Union (ERDF, "A way to make Europe") and the Centro Nacional de Genotipado-CEGEN (PRB2-ISCIII). BiB received core infrastructure funding from the Wellcome Trust (WT101597MA) and a joint grant from the UK Medical Research Council (MRC) and Economic and Social Science Research Council (ESRC) (MR/N024397/1). INMA data collections were supported by grants from the Instituto de Salud Carlos III, CIBERESP, and the Generalitat de Catalunya-CIRIT. KANC was funded by the grant of the Lithuanian Agency for Science Innovation and Technology (6-04-2014_31V-66). The Norwegian Mother, Father and Child Cohort Study is supported by the Norwegian Ministry of Health and Care Services and the Ministry of Education and Research. The Rhea project was financially supported by European projects (EU FP6-2003-Food-3-NewGeneris, EU FP6. STREP Hiwate, EU FP7 ENV.2007.1.2.2.2. Project No 211250 Escape, EU FP7-2008-ENV-1.2.1.4 Envirogenomarkers, EU FP7-HEALTH-2009- single stage CHICOS, EU FP7 ENV.2008.1.2.1.6. Proposal No 226285 ENRIECO, EU- FP7- HEALTH-2012 Proposal No 308333 HELIX), and the Greek Ministry of Health (Program of Prevention of obesity and neurodevelopmental disorders in preschool children, in Heraklion district, Crete, Greece: 2011-2014; "Rhea Plus": Primary Prevention Program of Environmental Risk Factors for Reproductive Health, and Child Health: 2012-15). We acknowledge support from the Spanish Ministry of Science and Innovation through the "Centro de Excelencia Severo Ochoa 2019-2023" Program (CEX2018-000806-S), and support from the Generalitat de Catalunya through the CERCA Program. MV-U and CR-A were supported by a FI fellowship from the Catalan Government (FI-DGR 2015 and #016FI_B 00272). MC received funding from Instituto Carlos III (Ministry of Economy and Competitiveness) (CD12/00563 and MS16/00128).

Genome-Wide Analysis for the Regulation of Gene Alternative Splicing by DNA Methylation Level in Glioma and its Prognostic Implications

Glioma is a primary high malignant intracranial tumor with poorly understood molecular mechanisms. Previous studies found that both DNA methylation modification and gene alternative splicing (AS) play a key role in tumorigenesis of glioma, and there is an obvious regulatory relationship between them. However, to date, no comprehensive study has been performed to analyze the influence of DNA methylation level on gene AS in glioma on a genome-wide scale. Here, we performed this study by integrating DNA methylation, gene expression, AS, disease risk methylation at position, and clinical data from 537 low-grade glioma (LGG) and glioblastoma (GBM) individuals. We first conducted a differential analysis of AS events and DNA methylation positions between LGG and GBM subjects, respectively. Then, we evaluated the influence of differential methylation positions on differential AS events. Further, Fisher’s exact test was used to verify our findings and identify potential key genes in glioma. Finally, we performed a series of analyses to investigate influence of these genes on the clinical prognosis of glioma. In total, we identified 130 glioma-related genes whose AS significantly affected by DNA methylation level. Eleven of them play an important role in glioma prognosis. In short, these results will help to better understand the pathogenesis of glioma.

Pup1 QTL Regulates Gene Expression Through Epigenetic Modification of DNA Under Phosphate Starvation Stress in Rice

Cytosine methylation, epigenetic DNA modification, is well known to regulate gene expression. Among the epigenetic modifications, 5-methylcytosine (5-mC) has been one of the extensively studied epigenetic changes responsible for regulating gene expression in animals and plants. Though a dramatic change in 5-mC content is observed at the genome level, the variation in gene expression is generally less than that it is expected. Only less is understood about the significance of 5-mC in gene regulation under P-starvation stress in plants. Using whole-genome bisulfite sequencing of a pair of rice [Pusa-44 and its near-isogenic line (NIL)-23 harboring Pup1 QTL] genotypes, we could decipher the role of Pup1 on DNA (de)methylation-mediated regulation of gene expression under P-starvation stress. We observed 13-15% of total cytosines to be methylated in the rice genome, which increased significantly under the stress. The number of differentially methylated regions (DMRs) for hypomethylation (6,068) was higher than those (5,279) for hypermethylated DMRs under the stress, particularly in root of NIL-23. Hypomethylation in CHH context caused upregulated expression of 489 genes in shoot and 382 genes in root of NIL-23 under the stress, wherein 387 genes in shoot and 240 genes in root were upregulated exclusively in NIL-23. Many of the genes for DNA methylation, a few for DNA demethylation, and RNA-directed DNA methylation were upregulated in root of NIL-23 under the stress. Methylation or demethylation of DNA in genic regions differentially affected gene expression. Correlation analysis for the distribution of DMRs and gene expression indicated the regulation of gene mainly through (de)methylation of promoter. Many of the P-responsive genes were hypomethylated or upregulated in roots of NIL-23 under the stress. Hypermethylation of gene body in CG, CHG, and CHH contexts caused up- or downregulated expression of transcription factors (TFs), P transporters, phosphoesterases, retrotransposon proteins, and other proteins. Our integrated transcriptome and methylome analyses revealed an important role of the Pup1 QTL in epigenetic regulation of the genes for transporters, TFs, phosphatases, carbohydrate metabolism, hormone-signaling, and chromatin architecture or epigenetic modifications in P-starvation tolerance. This provides insights into the molecular function of Pup1 in modulating gene expression through DNA (de)methylation, which might be useful in improving P-use efficiency or productivity of rice in P-deficient soil.

Epigenome-wide association studies: current knowledge, strategies and recommendations

The aetiology and pathophysiology of complex diseases are driven by the interaction between genetic and environmental factors. The variability in risk and outcomes in these diseases are incompletely explained by genetics or environmental risk factors individually. Therefore, researchers are now exploring the epigenome, a biological interface at which genetics and the environment can interact. There is a growing body of evidence supporting the role of epigenetic mechanisms in complex disease pathophysiology. Epigenome-wide association studies (EWASes) investigate the association between a phenotype and epigenetic variants, most commonly DNA methylation. The decreasing cost of measuring epigenome-wide methylation and the increasing accessibility of bioinformatic pipelines have contributed to the rise in EWASes published in recent years. Here, we review the current literature on these EWASes and provide further recommendations and strategies for successfully conducting them. We have constrained our review to studies using methylation data as this is the most studied epigenetic mechanism; microarray-based data as whole-genome bisulphite sequencing remains prohibitively expensive for most laboratories; and blood-based studies due to the non-invasiveness of peripheral blood collection and availability of archived DNA, as well as the accessibility of publicly available blood-cell-based methylation data. Further, we address multiple novel areas of EWAS analysis that have not been covered in previous reviews: (1) longitudinal study designs, (2) the chip analysis methylation pipeline (ChAMP), (3) differentially methylated region (DMR) identification paradigms, (4) methylation quantitative trait loci (methQTL) analysis, (5) methylation age analysis and (6) identifying cell-specific differential methylation from mixed cell data using statistical deconvolution.

Gaining insight into metabolic diseases from human genetic discoveries

Human large-scale genetic association studies have identified sequence variations at thousands of genetic risk loci that are more common in patients with diverse metabolic disease compared with healthy controls. While these genetic associations have been replicated in multiple large cohorts and sometimes can explain up to 50% of heritability, the molecular and cellular mechanisms affected by common genetic variation associated with metabolic disease remains mostly unknown. A variety of new genome-wide data types, in conjunction with novel biostatistical and computational analytical methodologies and foundational experimental technologies, are paving the way for a principled approach to systematic variant-to-function (V2F) studies for metabolic diseases, turning associated regions into causal variants, cell types and states of action, effector genes, and cellular and physiological mechanisms. Identification of new target genes and cellular programs for metabolic risk loci will improve mechanistic understanding of disease biology and identification of novel therapeutic strategies.

GWAS loci associated with Chagas cardiomyopathy influences DNA methylation levels

A recent genome-wide association study (GWAS) identified a locus in chromosome 11 associated with the chronic cardiac form of Chagas disease. Here we aimed to elucidate the potential functional mechanism underlying this genetic association by analyzing the correlation among single nucleotide polymorphisms (SNPs) and DNA methylation (DNAm) levels as cis methylation quantitative trait loci (cis-mQTL) within this region. A total of 2,611 SNPs were tested against 2,647 DNAm sites, in a subset of 37 chronic Chagas cardiomyopathy patients and 20 asymptomatic individuals from the GWAS. We identified 6,958 significant cis-mQTLs (False Discovery Rate [FDR]<0.05) at 1 Mb each side of the GWAS leading variant, where six of them potentially modulate the expression of the SAC3D1 gene, the reported gene in the previous GWAS. In addition, a total of 268 cis-mQTLs showed differential methylation between chronic Chagas cardiomyopathy patients and asymptomatic individuals. The most significant cis-mQTLs mapped in the gene bodies of POLA2 (FDR = 1.04x10^-11), PLAAT3 (FDR = 7.22x10^-03), and CCDC88B (FDR = 1.89x10^-02) that have been associated with cardiovascular and hematological traits in previous studies. One of the most relevant interactions correlated with hypermethylation of CCDC88B. This gene is involved in the inflammatory response, and its methylation and expression levels have been previously reported in Chagas cardiomyopathy. Our findings support the functional relevance of the previously associated genomic region, highlighting the regulation of novel genes that could play a role in the chronic cardiac form of the disease.

Genome-wide association studies (GWAS) have provided extensive information regarding the genetic component of complex traits, including parasitic diseases such as Chagas disease. However, these associations mapped in regulatory regions of the genome and assigning them a functional consequence have been cumbersome. In this study we aimed to evaluate the functional mechanism underlying the previously reported genomic association with chronic Chagas cardiomyopathy, by assessing the correlation between methylation changes and the underlying genetic variations within the region. These methylation quantitative trait loci (mQTLs) may be involved in gene expression regulation. We identified mQTLs in three genes that have been associated with cardiovascular diseases in previous studies. Interestingly, one of these genes was previously identified as differentially methylated and expressed in heart biopsies of chronic Chagas cardiomyopathy patients. Our results suggest novel genes that could play a role in the chronic Chagas cardiomyopathy, evidencing the functional relevance of the previously associated loci.

SLC25A24 gene methylation and gray matter volume in females with and without conduct disorder: an exploratory epigenetic neuroimaging study

Conduct disorder (CD), a psychiatric disorder characterized by a repetitive pattern of antisocial behaviors, results from a complex interplay between genetic and environmental factors. The clinical presentation of CD varies both according to the individual's sex and level of callous-unemotional (CU) traits, but it remains unclear how genetic and environmental factors interact at the molecular level to produce these differences. Emerging evidence in males implicates methylation of genes associated with socio-affective processes. Here, we combined an epigenome-wide association study with structural neuroimaging in 51 females with CD and 59 typically developing (TD) females to examine DNA methylation in relation to CD, CU traits, and gray matter volume (GMV). We demonstrate an inverse pattern of correlation between CU traits and methylation of a chromosome 1 region in CD females (positive) as compared to TD females (negative). The identified region spans exon 1 of the SLC25A24 gene, central to energy metabolism due to its role in mitochondrial function. Increased SLC25A24 methylation was also related to lower GMV in multiple brain regions in the overall cohort. These included the superior frontal gyrus, prefrontal cortex, and supramarginal gyrus, secondary visual cortex and ventral posterior cingulate cortex, which are regions that have previously been implicated in CD and CU traits. While our findings are preliminary and need to be replicated in larger samples, they provide novel evidence that CU traits in females are associated with methylation levels in a fundamentally different way in CD and TD, which in turn may relate to observable variations in GMV across the brain.

Heterogeneous information network and its application to human health and disease

The molecular components with the functional interdependencies in human cell form complicated biological network. Diseases are mostly caused by the perturbations of the composite of the interaction multi-biomolecules, rather than an abnormality of a single biomolecule. Furthermore, new biological functions and processes could be revealed by discovering novel biological entity relationships. Hence, more and more biologists focus on studying the complex biological system instead of the individual biological components. The emergence of heterogeneous information network (HIN) offers a promising way to systematically explore complicated and heterogeneous relationships between various molecules for apparently distinct phenotypes. In this review, we first present the basic definition of HIN and the biological system considered as a complex HIN. Then, we discuss the topological properties of HIN and how these can be applied to detect network motif and functional module. Afterwards, methodologies of discovering relationships between disease and biomolecule are presented. Useful insights on how HIN aids in drug development and explores human interactome are provided. Finally, we analyze the challenges and opportunities for uncovering combinatorial patterns among pharmacogenomics and cell-type detection based on single-cell genomic data.

Identification of tissue-specific and common methylation quantitative trait loci in healthy individuals using MAGAR

Background

Understanding the influence of genetic variants on DNA methylation is fundamental for the interpretation of epigenomic data in the context of disease. There is a need for systematic approaches not only for determining methylation quantitative trait loci (methQTL), but also for discriminating general from cell type-specific effects.

Results

Here, we present a two-step computational framework MAGAR (https://bioconductor.org/packages/MAGAR), which fully supports the identification of methQTLs from matched genotyping and DNA methylation data, and additionally allows for illuminating cell type-specific methQTL effects. In a pilot analysis, we apply MAGAR on data in four tissues (ileum, rectum, T cells, B cells) from healthy individuals and demonstrate the discrimination of common from cell type-specific methQTLs. We experimentally validate both types of methQTLs in an independent data set comprising additional cell types and tissues. Finally, we validate selected methQTLs located in the PON1, ZNF155, and NRG2 genes by ultra-deep local sequencing. In line with previous reports, we find cell type-specific methQTLs to be preferentially located in enhancer elements.

Conclusions

Our analysis demonstrates that a systematic analysis of methQTLs provides important new insights on the influences of genetic variants to cell type-specific epigenomic variation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13072-021-00415-6.

Genome-wide sequencing-based identification of methylation quantitative trait loci and their role in schizophrenia risk

DNA methylation (DNAm) is an epigenetic regulator of gene expression and a hallmark of gene-environment interaction. Using whole-genome bisulfite sequencing, we have surveyed DNAm in 344 samples of human postmortem brain tissue from neurotypical subjects and individuals with schizophrenia. We identify genetic influence on local methylation levels throughout the genome, both at CpG sites and CpH sites, with 86% of SNPs and 55% of CpGs being part of methylation quantitative trait loci (meQTLs). These associations can further be clustered into regions that are differentially methylated by a given SNP, highlighting the genes and regions with which these loci are epigenetically associated. These findings can be used to better characterize schizophrenia GWAS-identified variants as epigenetic risk variants. Regions differentially methylated by schizophrenia risk-SNPs explain much of the heritability associated with risk loci, despite covering only a fraction of the genomic space. We provide a comprehensive, single base resolution view of association between genetic variation and genomic methylation, and implicate schizophrenia GWAS-associated variants as influencing the epigenetic plasticity of the brain.

The authors provide a comprehensive, single base resolution view of association between genetic variation and DNA methylation in human brain. They also show that heritability attributed to schizophrenia GWAS-associated variants reflects the epigenetic plasticity of the brain.

Epigenetic Alterations Related to Gestational Diabetes Mellitus

Gestational diabetes mellitus (GDM) is the most common metabolic complication in pregnancy, which affects the future health of both the mother and the newborn. Its pathophysiology involves nutritional, hormonal, immunological, genetic and epigenetic factors. Among the latter, it has been observed that alterations in DNA (deoxyribonucleic acid) methylation patterns and in the levels of certain micro RNAs, whether in placenta or adipose tissue, are related to well-known characteristics of the disease, such as hyperglycemia, insulin resistance, inflammation and excessive placental growth. Furthermore, epigenetic alterations of gestational diabetes mellitus are observable in maternal blood, although their pathophysiological roles are completely unknown. Despite this, it has not been possible to determine the causes of the epigenetic characteristics of GDM, highlighting the need for integral and longitudinal studies. Based on this, this article summarizes the most relevant and recent studies on epigenetic alterations in placenta, adipose tissue and maternal blood associated with GDM in order to provide the reader with a general overview of the subject and indicate future research topics.

MOSTWAS: Multi-Omic Strategies for Transcriptome-Wide Association Studies

Traditional predictive models for transcriptome-wide association studies (TWAS) consider only single nucleotide polymorphisms (SNPs) local to genes of interest and perform parameter shrinkage with a regularization process. These approaches ignore the effect of distal-SNPs or other molecular effects underlying the SNP-gene association. Here, we outline multi-omics strategies for transcriptome imputation from germline genetics to allow more powerful testing of gene-trait associations by prioritizing distal-SNPs to the gene of interest. In one extension, we identify mediating biomarkers (CpG sites, microRNAs, and transcription factors) highly associated with gene expression and train predictive models for these mediators using their local SNPs. Imputed values for mediators are then incorporated into the final predictive model of gene expression, along with local SNPs. In the second extension, we assess distal-eQTLs (SNPs associated with genes not in a local window around it) for their mediation effect through mediating biomarkers local to these distal-eSNPs. Distal-eSNPs with large indirect mediation effects are then included in the transcriptomic prediction model with the local SNPs around the gene of interest. Using simulations and real data from ROS/MAP brain tissue and TCGA breast tumors, we show considerable gains of percent variance explained (1-2% additive increase) of gene expression and TWAS power to detect gene-trait associations. This integrative approach to transcriptome-wide imputation and association studies aids in identifying the complex interactions underlying genetic regulation within a tissue and important risk genes for various traits and disorders.

A gene-level methylome-wide association analysis identifies novel Alzheimer's disease genes

MOTIVATION

Transcriptome-wide association studies (TWAS) have successfully facilitated the discovery of novel genetic risk loci for many complex traits, including late-onset Alzheimer's disease (AD). However, most existing TWAS methods rely only on gene expression and ignore epigenetic modification (i.e., DNA methylation) and functional regulatory information (i.e., enhancer-promoter interactions), both of which contribute significantly to the genetic basis of AD.

RESULTS

We develop a novel gene-level association testing method that integrates genetically regulated DNA methylation and enhancer-target gene pairs with genome-wide association study (GWAS) summary results. Through simulations, we show that our approach, referred to as the CMO (cross methylome omnibus) test, yielded well controlled type I error rates and achieved much higher statistical power than competing methods under a wide range of scenarios. Furthermore, compared with TWAS, CMO identified an average of 124% more associations when analyzing several brain imaging-related GWAS results. By analyzing to date the largest AD GWAS of 71,880 cases and 383,378 controls, CMO identified six novel loci for AD, which have been ignored by competing methods.

AVAILABILITY

Software: https://github.com/ChongWuLab/CMO.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Differential DNA Methylation Landscape in Skin Fibroblasts from African Americans with Systemic Sclerosis

The etiology and reasons underlying the ethnic disparities in systemic sclerosis (SSc) remain unknown. African Americans are disproportionally affected by SSc and yet are underrepresented in research. The aim of this study was to comprehensively investigate the association of DNA methylation levels with SSc in dermal fibroblasts from patients of African ancestry. Reduced representation bisulfite sequencing (RRBS) was performed on primary dermal fibroblasts from 15 SSc patients and 15 controls of African ancestry, and over 3.8 million CpG sites were tested for differential methylation patterns between cases and controls. The dermal fibroblasts from African American patients exhibited widespread reduced DNA methylation. Differentially methylated CpG sites were most enriched in introns and intergenic regions while depleted in 5' UTR, promoters, and CpG islands. Seventeen genes and eleven promoters showed significant differential methylation, mostly in non-coding RNA genes and pseudogenes. Gene set enrichment analysis (GSEA) and gene ontology (GO) analyses revealed an enrichment of pathways related to interferon signaling and mesenchymal differentiation. The hypomethylation of DLX5 and TMEM140 was accompanied by these genes' overexpression in patients but underexpression for lncRNA MGC12916. These data show that differential methylation occurs in dermal fibroblasts from African American patients with SSc and identifies novel coding and non-coding genes.

Allele-specific expression: applications in cancer and technical considerations

Allele-specific gene expression can influence disease traits. Non-coding germline genetic variants that alter regulatory elements can cause allele-specific gene expression and contribute to cancer susceptibility. In tumors, both somatic copy number alterations and somatic single nucleotide variants have been shown to lead to allele-specific expression of genes, many of which are considered drivers of tumor growth. Here, we review recent studies revealing the pervasive presence of this phenomenon in cancer susceptibility and progression. Furthermore, we underscore the importance of careful experimental design and computational analysis for accurate allelic expression quantification and avoidance of false positives. Finally, we discuss additional methodological challenges encountered in cancer studies and in the burgeoning field of single-cell transcriptomics.

Machine learning-based integrative analysis of methylome and transcriptome identifies novel prognostic DNA methylation signature in uveal melanoma

Uveal melanoma (UVM) is the most common primary intraocular human malignancy with a high mortality rate. Aberrant DNA methylation has rapidly emerged as a diagnostic and prognostic signature in many cancers. However, such DNA methylation signature available in UVM remains limited. In this study, we performed a genome-wide integrative analysis of methylome and transcriptome and identified 40 methylation-driven prognostic genes (MDPGs) associated with the tumorigenesis and progression of UVM. Then, we proposed a machine-learning-based discovery and validation strategy to identify a DNA methylation-driven signature (10MeSig) composing of 10 MDPGs (AZGP1, BAI1, CCDC74A, FUT3, PLCD1, S100A4, SCN8A, SEMA3B, SLC25A38 and SLC44A3), which stratified 80 patients of the discovery cohort into two risk subtypes with significantly different overall survival (HR = 29, 95% CI: 6.7-126, P < 0.001). The 10MeSig was validated subsequently in an independent cohort with 57 patients and yielded a similar prognostic value (HR = 2.1, 95% CI: 1.2-3.7, P = 0.006). Multivariable Cox regression analysis showed that the 10MeSig is an independent predictive factor for the survival of patients with UVM. With a prospective validation study, this 10MeSig will improve clinical decisions and provide new insights into the pathogenesis of UVM.

Profiling Novel Alternative Splicing within Multiple Tissues Provides Useful Insights into Porcine Genome Annotation

Alternative splicing (AS) is a process during gene expression that results in a single gene coding for different protein variants. AS contributes to transcriptome and proteome diversity. In order to characterize AS in pigs, genome-wide transcripts and AS events were detected using RNA sequencing of 34 different tissues in Duroc pigs. In total, 138,403 AS events and 29,270 expressed genes were identified. An alternative donor site was the most common AS form and accounted for 44% of the total AS events. The percentage of the other three AS forms (exon skipping, alternative acceptor site, and intron retention) was approximately 19%. The results showed that the most common AS events involving alternative donor sites could produce different transcripts or proteins that affect the biological processes. The expression of genes with tissue-specific AS events showed that gene functions were consistent with tissue functions. AS increased proteome diversity and resulted in novel proteins that gained or lost important functional domains. In summary, these findings extend porcine genome annotation and highlight roles that AS could play in determining tissue identity.

Involvement of lncRNAs in celiac disease pathogenesis

Celiac disease (CD) is an immune-mediated disease that develops in genetically susceptible individuals upon gluten exposure. Human Leukocyte Antigen (HLA) genes in the Major Histocompatibility Complex (MHC) have been described to represent the 40% of the genetic risk to develop CD. Aiming to gain understanding of the genetic involvement in CD, high throughput studies have been performed, revealing that many CD-associated variants are located in non-coding regions, hindering the study of the functional implications of these single nucleotide polymorphisms (SNPs). In the last decade, long non-coding RNAs (lncRNAs) have been described to be influenced by disease-associated SNPs and to drive many important mechanisms involved in the development of inflammatory diseases. Here we describe the lncRNAs identified and characterized in the context of celiac disease and highlight the importance of the study of these molecules in inflammatory and autoimmune disorders.

Association of CNVs with methylation variation

Germline copy number variants (CNVs) and single-nucleotide polymorphisms (SNPs) form the basis of inter-individual genetic variation. Although the phenotypic effects of SNPs have been extensively investigated, the effects of CNVs is relatively less understood. To better characterize mechanisms by which CNVs affect cellular phenotype, we tested their association with variable CpG methylation in a genome-wide manner. Using paired CNV and methylation data from the 1000 genomes and HapMap projects, we identified genome-wide associations by methylation quantitative trait locus (mQTL) analysis. We found individual CNVs being associated with methylation of multiple CpGs and vice versa. CNV-associated methylation changes were correlated with gene expression. CNV-mQTLs were enriched for regulatory regions, transcription factor-binding sites (TFBSs), and were involved in long-range physical interactions with associated CpGs. Some CNV-mQTLs were associated with methylation of imprinted genes. Several CNV-mQTLs and/or associated genes were among those previously reported by genome-wide association studies (GWASs). We demonstrate that germline CNVs in the genome are associated with CpG methylation. Our findings suggest that structural variation together with methylation may affect cellular phenotype.

Alternative splicing of MR1 regulates antigen presentation to MAIT cells

Mucosal Associated Invariant T (MAIT) cells can sense intracellular infection by a broad array of pathogens. These cells are activated upon encountering microbial antigen(s) displayed by MR1 on the surface of an infected cell. Human MR1 undergoes alternative splicing. The full-length isoform, MR1A, can activate MAIT cells, while the function of the isoforms, MR1B and MR1C, are incompletely understood. In this report, we sought to characterize the expression and function of these splice variants. Using a transcriptomic analysis in conjunction with qPCR, we find that that MR1A and MR1B transcripts are widely expressed. However only MR1A can present mycobacterial antigen to MAIT cells. Coexpression of MR1B with MR1A decreases MAIT cell activation following bacterial infection. Additionally, expression of MR1B prior to MR1A lowers total MR1A abundance, suggesting competition between MR1A and MR1B for either ligands or chaperones required for folding and/or trafficking. Finally, we evaluated CD4/CD8 double positive thymocytes expressing surface MR1. Here, we find that relative expression of MR1A/MR1B transcript is associated with the prevalence of MR1 + CD4/CD8 cells in the thymus. Our results suggest alternative splicing of MR1 represents a means of regulating MAIT activation in response to microbial ligand(s).

A case-control study on the SNP309T → G and 40-bp Del1518 of the MDM2 gene and a systematic review for MDM2 polymorphisms in the patients with breast cancer

Objective

The current research was conducted to study the association between the SNP309 and del1518 polymorphisms with the breast cancer in the patients with the Kurdish ethnic background from western Iran. Also, a systematic review of the relevant case‐control studies on the MDM2 polymorphisms in the patients with breast cancer was performed.

Methodology

Two mL of peripheral blood was taken from 100 patients with breast cancer and 100 healthy individuals. The frequencies of MDM2 SNP309 and del1518 genotypes and alleles were determined using the PCR‐RFLP and PCR methods, respectively.

Results

The frequency of the TT, TG, and GG of MDM2‐SNP309 genotypes in the patients was obtained as 23%, 52%, and 25%, and they were equal to 22%, 40%, and 38% in the control group, respectively. Also, considering the MDM2‐del1518 polymorphism, the frequencies of ins/ins, ins/del, and del/del genotypes were equal to 52%, 41%, and 7% in the breast cancer group and they were equal to 62, 30, and 8% in the control group, respectively. Analysis of the results indicated that the GG genotype plays a protective role for the breast cancer in the recessive model (GG vs TT + TG) of SNP309 (χ² = 3.916, P = .048, and OR = 0.54).

Conclusion

Our findings revealed that the GG genotype of MDM2‐SNP309 can play a protective role in the breast cancer disease. Also, our systematic review indicated that the SNP309, SNP285, and del1518 of MDM2 gene in different populations mostly did not have a significant association with the risk of breast cancer.

sn-spMF: matrix factorization informs tissue-specific genetic regulation of gene expression

Genetic regulation of gene expression, revealed by expression quantitative trait loci (eQTLs), exhibits complex patterns of tissue-specific effects. Characterization of these patterns may allow us to better understand mechanisms of gene regulation and disease etiology. We develop a constrained matrix factorization model, sn-spMF, to learn patterns of tissue-sharing and apply it to 49 human tissues from the Genotype-Tissue Expression (GTEx) project. The learned factors reflect tissues with known biological similarity and identify transcription factors that may mediate tissue-specific effects. sn-spMF, available at https://github.com/heyuan7676/ts_eQTLs , can be applied to learn biologically interpretable patterns of eQTL tissue-specificity and generate testable mechanistic hypotheses.

Expression Quantitative Trait Methylation Analysis Reveals Methylomic Associations With Gene Expression in Childhood Asthma

BACKGROUND

Nasal (airway) epithelial methylation profiles have been associated with asthma, but the effects of such profiles on expression of distant cis-genes are largely unknown.

RESEARCH QUESTION

To identify genes whose expression is associated with proximal and distal CpG probes (within 1 Mb), and to assess whether and how such genes are differentially expressed in atopic asthma.

STUDY DESIGN AND METHODS

Genome-wide expression quantitative trait methylation (eQTM) analysis in nasal epithelium from Puerto Rican subjects (aged 9-20 years) with (n = 219) and without (n = 236) asthma. After the eQTM analysis, a Gene Ontology Enrichment analysis was conducted for the top 500 eQTM genes, and mediation analyses were performed to identify paths from DNA methylation to atopic asthma through gene expression. Asthma was defined as physician-diagnosed asthma and wheeze in the previous year, and atopy was defined as at least one positive IgE to allergens. Atopic asthma was defined as the presence of both atopy and asthma.

RESULTS

We identified 16,867 significant methylation-gene expression pairs (false-discovery rate-adjusted P < .01) in nasal epithelium from study participants. Most eQTM methylation probes were distant (average distance, ∼378 kb) from their target genes, and also more likely to be located in enhancer regions of their target genes in lung tissue than control probes. The top 500 eQTM genes were enriched in pathways for immune processes and epithelial integrity and were more likely to have been previously identified as differentially expressed in atopic asthma. In a mediation analysis, we identified 5,934 paths through which methylation markers could affect atopic asthma through gene expression in nasal epithelium.

INTERPRETATION

Previous epigenome-wide association studies of asthma have estimated the effects of DNA methylation markers on expression of nearby genes in airway epithelium. Our findings suggest that distant epigenetic regulation of gene expression in airway epithelium plays a role in atopic asthma.

Integration of methylation QTL and enhancer-target gene maps with schizophrenia GWAS summary results identifies novel genes

MOTIVATION

Most trait-associated genetic variants identified in genome-wide association studies (GWASs) are located in non-coding regions of the genome and thought to act through their regulatory roles.

RESULTS

To account for enriched association signals in DNA regulatory elements, we propose a novel and general gene-based association testing strategy that integrates enhancer-target gene pairs and methylation quantitative trait locus data with GWAS summary results; it aims to both boost statistical power for new discoveries and enhance mechanistic interpretability of any new discovery. By reanalyzing two large-scale schizophrenia GWAS summary datasets, we demonstrate that the proposed method could identify some significant and novel genes (containing no genome-wide significant SNPs nearby) that would have been missed by other competing approaches, including the standard and some integrative gene-based association methods, such as one incorporating enhancer-target gene pairs and one integrating expression quantitative trait loci.

AVAILABILITY AND IMPLEMENTATION

Software: wuchong.org/egmethyl.html.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

The therapeutic potential of exercise for neuropsychiatric diseases: A review

Exercise is known to have a myriad of health benefits. There is much to be learned from the effects of exercise and its potential for prevention, attenuation and treatment of multiple neuropsychiatric diseases and behavioral disorders. Furthermore, recent data and research on exercise benefits with respect to major health crises, such as, that of opioid and general substance use disorders, make it very important to better understand and review the mechanisms of exercise and how it could be utilized for effective treatments or adjunct treatments for these diseases. In addition, mechanisms, epigenetics and sex differences are examined and discussed in terms of future research implications.

Collective effects of long-range DNA methylations predict gene expressions and estimate phenotypes in cancer

DNA methylation of various genomic regions has been found to be associated with gene expression in diverse biological contexts. However, most genome-wide studies have focused on the effect of (1) methylation in cis, not in trans and (2) a single CpG, not the collective effects of multiple CpGs, on gene expression. In this study, we developed a statistical machine learning model, geneEXPLORE (gene expression prediction by long-range epigenetics), that quantifies the collective effects of both cis- and trans- methylations on gene expression. By applying geneEXPLORE to The Cancer Genome Atlas (TCGA) breast and 10 other types of cancer data, we found that most genes are associated with methylations of as much as 10 Mb from the promoters or more, and the long-range methylation explains 50% of the variation in gene expression on average, far greater than cis-methylation. geneEXPLORE outperforms competing methods such as BioMethyl and MethylXcan. Further, the predicted gene expressions could predict clinical phenotypes such as breast tumor status and estrogen receptor status (AUC = 0.999, 0.94 respectively) as accurately as the measured gene expression levels. These results suggest that geneEXPLORE provides a means for accurate imputation of gene expression, which can be further used to predict clinical phenotypes.

IGREX for quantifying the impact of genetically regulated expression on phenotypes

By leveraging existing GWAS and eQTL resources, transcriptome-wide association studies (TWAS) have achieved many successes in identifying trait-associations of genetically regulated expression (GREX) levels. TWAS analysis relies on the shared GREX variation across GWAS and the reference eQTL data, which depends on the cellular conditions of the eQTL data. Considering the increasing availability of eQTL data from different conditions and the often unknown trait-relevant cell/tissue-types, we propose a method and tool, IGREX, for precisely quantifying the proportion of phenotypic variation attributed to the GREX component. IGREX takes as input a reference eQTL panel and individual-level or summary-level GWAS data. Using eQTL data of 48 tissue types from the GTEx project as a reference panel, we evaluated the tissue-specific IGREX impact on a wide spectrum of phenotypes. We observed strong GREX effects on immune-related protein biomarkers. By incorporating trans-eQTLs and analyzing genetically regulated alternative splicing events, we evaluated new potential directions for TWAS analysis.

Large epigenome-wide association study of childhood ADHD identifies peripheral DNA methylation associated with disease and polygenic risk burden

Epigenetic variation in peripheral tissues is being widely studied as a molecular biomarker of complex disease and disease-related exposures. To date, few studies have examined differences in DNA methylation associated with attention-deficit hyperactivity disorder (ADHD). In this study, we profiled genetic and methylomic variation across the genome in saliva samples from children (age 7-12 years) with clinically established ADHD (N = 391) and nonpsychiatric controls (N = 213). We tested for differentially methylated positions (DMPs) associated with both ADHD diagnosis and ADHD polygenic risk score, by using linear regression models including smoking, medication effects, and other potential confounders in our statistical models. Our results support previously reported associations between ADHD and DNA methylation levels at sites annotated to VIPR2, and identify several novel disease-associated DMPs (p < 1e-5), although none of them were genome-wide significant. The two top-ranked, ADHD-associated DMPs (cg17478313 annotated to SLC7A8 and cg21609804 annotated to MARK2) are also significantly associated with nearby SNPs (p = 1.2e-46 and p = 2.07e-59), providing evidence that disease-associated DMPs are under genetic control. We also report a genome-wide significant association between ADHD polygenic risk and variable DNA methylation at a site annotated to the promoter of GART and SON (p = 6.71E-8). Finally, we show that ADHD-associated SNPs colocalize with SNPs associated with methylation levels in saliva. This is the first large-scale study of DNA methylation in children with ADHD. Our results represent novel epigenetic biomarkers for ADHD that may be useful for patient stratification, reinforce the importance of genetic effects on DNA methylation, and provide plausible molecular mechanisms for ADHD risk variants.