CpG Islands Undermethylation in Human Genomic Regions under Selective Pressure

DNA methylation at CpG islands (CGIs) is one of the most intensively studied epigenetic mechanisms. It is fundamental for cellular differentiation and control of transcriptional potential. DNA methylation is involved also in several processes that are central to evolutionary biology, including phenotypic plasticity and evolvability. In this study, we explored the relationship between CpG islands methylation and signatures of selective pressure in Homo Sapiens, using a computational biology approach. By analyzing methylation data of 25 cell lines from the Encyclopedia of DNA Elements (ENCODE) Consortium, we compared the DNA methylation of CpG islands in genomic regions under selective pressure with the methylation of CpG islands in the remaining part of the genome. To define genomic regions under selective pressure, we used three different methods, each oriented to provide distinct information about selective events. Independently of the method and of the cell type used, we found evidences of undermethylation of CGIs in human genomic regions under selective pressure. Additionally, by analyzing SNP frequency in CpG islands, we demonstrated that CpG islands in regions under selective pressure show lower genetic variation. Our findings suggest that the CpG islands in regions under selective pressure seem to be somehow more “protected” from methylation when compared with other regions of the genome

Promoter Methylation in Head and Neck Squamous Cell Carcinoma Cell Lines Is Significantly Different than Methylation in Primary Tumors and Xenografts

Studies designed to identify novel methylation events related to cancer often employ cancer cell lines in the discovery phase of the experiments and have a relatively low rate of discovery of cancer-related methylation events. An alternative algorithm for discovery of novel methylation in cancer uses primary tumor-derived xenografts instead of cell lines as the primary source of nucleic acid for evaluation. We evaluated DNA extracted from primary head and neck squamous cell carcinomas (HNSCC), xenografts grown from these primary tumors in nude mice, HNSCC-derived cell lines, normal oral mucosal samples, and minimally transformed oral keratinocyte-derived cell lines using Illumina Infinum Humanmethylation 27 genome-wide methylation microarrays. We found >2,200 statistically significant methylation differences between cancer cell lines and primary tumors and when comparing normal oral mucosa to keratinocyte cell lines. We found no statistically significant promoter methylation differences between primary tumor xenografts and primary tumors. This study demonstrates that tumor-derived xenografts are highly accurate representations of promoter methylation in primary tumors and that cancer derived cell lines have significant drawbacks for discovery of promoter methylation alterations in primary tumors. These findings also support use of primary tumor xenografts for the study of methylation in cancer, drug discovery, and the development of personalized cancer treatments

Genetic Control of mRNA Splicing as a Potential Mechanism for Incomplete Penetrance of Rare Coding Variants

Exonic variants present some of the strongest links between genotype and phenotype. However, these variants can have significant inter-individual pathogenicity differences, known as variable penetrance. In this study, we propose a model where genetically controlled mRNA splicing modulates the pathogenicity of exonic variants. By first cataloging exonic inclusion from RNA-sequencing data in GTEx V8, we find that pathogenic alleles are depleted on highly included exons. Using a large-scale phased whole genome sequencing data from the TOPMed consortium, we observe that this effect may be driven by common splice-regulatory genetic variants, and that natural selection acts on haplotype configurations that reduce the transcript inclusion of putatively pathogenic variants, especially when limiting to haploinsufficient genes. Finally, we test if this effect may be relevant for autism risk using families from the Simons Simplex Collection, but find that splicing of pathogenic alleles has a penetrance reducing effect here as well. Overall, our results indicate that common splice-regulatory variants may play a role in reducing the damaging effects of rare exonic variants

Interaction between genetic and epigenetic variation defines gene expression patterns at the asthma-associated locus 17q12-q21 in lymphoblastoid cell lines

Phenotypic variation results from variation in gene expression, which is modulated by genetic and/or epigenetic factors. To understand the molecular basis of human disease, interaction between genetic and epigenetic factors needs to be taken into account. The asthma-associated region 17q12-q21 harbors three genes, the zona pellucida binding protein 2 (ZPBP2), gasdermin B (GSDMB) and ORM1-like 3 (ORMDL3), that show allele-specific differences in expression levels in lymphoblastoid cell lines (LCLs) and CD4+ T cells. Here, we report a molecular dissection of allele-specific transcriptional regulation of the genes within the chromosomal region 17q12-q21 combining in vitro transfection, formaldehyde-assisted isolation of regulatory elements, chromatin immunoprecipitation and DNA methylation assays in LCLs. We found that a single nucleotide polymorphism rs4795397 influences the activity of ZPBP2 promoter in vitro in an allele-dependent fashion, and also leads to nucleosome repositioning on the asthma-associated allele. However, variable methylation of exon 1 of ZPBP2 masks the strong genetic effect on ZPBP2 promoter activity in LCLs. In contrast, the ORMDL3 promoter is fully unmethylated, which allows detection of genetic effects on its transcription. We conclude that the cis-regulatory effects on 17q12-q21 gene expression result from interaction between several regulatory polymorphisms and epigenetic factors within the cis-regulatory haplotype region

Variation in Array Size, Monomer Composition and Expression of the Macrosatellite DXZ4

Macrosatellites are some of the most polymorphic regions of the human genome, yet many remain uncharacterized despite the association of some arrays with disease susceptibility. This study sought to explore the polymorphic nature of the X-linked macrosatellite DXZ4. Four aspects of DXZ4 were explored in detail, including tandem repeat copy number variation, array instability, monomer sequence polymorphism and array expression. DXZ4 arrays contained between 12 and 100 3.0 kb repeat units with an average array containing 57. Monomers were confirmed to be arranged in uninterrupted tandem arrays by restriction digest analysis and extended fiber FISH, and therefore DXZ4 encompasses 36–288 kb of Xq23. Transmission of DXZ4 through three generations in three families displayed a high degree of meiotic instability (8.3%), consistent with other macrosatellite arrays, further highlighting the unstable nature of these sequences in the human genome. Subcloning and sequencing of complete DXZ4 monomers identified numerous single nucleotide polymorphisms and alleles for the three microsatellite repeats located within each monomer. Pairwise comparisons of DXZ4 monomer sequences revealed that repeat units from an array are more similar to one another than those originating from different arrays. RNA fluorescence in situ hybridization revealed significant variation in DXZ4 expression both within and between cell lines. DXZ4 transcripts could be detected originiating from both the active and inactive X chromosome. Expression levels of DXZ4 varied significantly between males, but did not relate to the size of the array, nor did inheritance of the same array result in similar expression levels. Collectively, these studies provide considerable insight into the polymorphic nature of DXZ4, further highlighting the instability and variation potential of macrosatellites in the human genome

Human JAK1 gain of function causes dysregulated myelopoeisis and severe allergic inflammation

Primary atopic disorders are a group of inborn errors of immunity that skew the immune system toward severe allergic disease. Defining the biology underlying these extreme monogenic phenotypes reveals shared mechanisms underlying common polygenic allergic disease and identifies potential drug targets. Germline gain-of-function (GOF) variants in JAK1 are a cause of severe atopy and eosinophilia. Modeling the JAK1GOF (p.A634D) variant in both zebrafish and human induced pluripotent stem cells (iPSCs) revealed enhanced myelopoiesis. RNA-Seq of JAK1GOF human whole blood, iPSCs, and transgenic zebrafish revealed a shared core set of dysregulated genes involved in IL-4, IL-13, and IFN signaling. Immunophenotypic and transcriptomic analysis of patients carrying a JAK1GOF variant revealed marked Th cell skewing. Moreover, long-term ruxolitinib treatment of 2 children carrying the JAK1GOF (p.A634D) variant remarkably improved their growth, eosinophilia, and clinical features of allergic inflammation. This work highlights the role of JAK1 signaling in atopic immune dysregulation and the clinical impact of JAK1/2 inhibition in treating eosinophilic and allergic disease.Transplantation and immunomodulatio