Genomic imprinting variations in the mouse type 3 deiodinase gene between tissues and brain regions.

The Dio3 gene, which encodes for the type 3 deiodinase (D3), controls thyroid hormone (TH) availability. The lack of D3 in mice results in tissue overexposure to TH and a broad neuroendocrine phenotype. Dio3 is an imprinted gene, preferentially expressed from the paternally inherited allele in the mouse fetus. However, heterozygous mice with paternal inheritance of the inactivating Dio3 mutation exhibit an attenuated phenotype when compared with that of Dio3 null mice. To investigate this milder phenotype, the allelic expression of Dio3 was evaluated in different mouse tissues. Preferential allelic expression of Dio3 from the paternal allele was observed in fetal tissues and neonatal brain regions, whereas the biallelic Dio3 expression occurred in the developing eye, testes, and cerebellum and in the postnatal brain neocortex, which expresses a larger Dio3 mRNA transcript. The newborn hypothalamus manifests the highest degree of Dio3 expression from the paternal allele, compared with other brain regions, and preferential allelic expression of Dio3 in the brain relaxed in late neonatal life. A methylation analysis of two regulatory regions of the Dio3 imprinted domain revealed modest but significant differences between tissues, but these did not consistently correlate with the observed patterns of Dio3 allelic expression. Deletion of the Dio3 gene and promoter did not result in significant changes in the tissue-specific patterns of Dio3 allelic expression. These results suggest the existence of unidentified epigenetic determinants of tissue-specific Dio3 imprinting. The resulting variation in the Dio3 allelic expression between tissues likely explains the phenotypic variation that results from paternal Dio3 haploinsufficiency.This is the final version of the article. It is available from the Endocrine Society in Molecular Endocrinology here: http://press.endocrine.org/doi/pdf/10.1210/me.2014-1210

Blood Gene Expression and Immune Cell Subtypes Associated with Chronic Obstructive Pulmonary Disease Exacerbations.

Rationale: Acute exacerbations of chronic obstructive pulmonary disease (AE-COPDs) are associated with a significant disease burden. Blood immune phenotyping may improve our understanding of a COPD endotype at increased risk of exacerbations. Objective: To determine the relationship between the transcriptome of circulating leukocytes and COPD exacerbations. Methods: Blood RNA sequencing data (n = 3,618) from the COPDGene (Genetic Epidemiology of COPD) study were analyzed. Blood microarray data (n = 646) from the ECLIPSE (Evaluation of COPD Longitudinally to Identify Predictive Surrogate Endpoints) study were used for validation. We tested the association between blood gene expression and AE-COPDs. We imputed the abundance of leukocyte subtypes and tested their association with prospective AE-COPDs. Flow cytometry was performed on blood in SPIROMICS (Subpopulations and Intermediate Outcomes in COPD Study) (n = 127), and activation markers for T cells were tested for association with prospective AE-COPDs. Measurements and Main Results: Exacerbations were reported 4,030 and 2,368 times during follow-up in COPDGene (5.3 ± 1.7 yr) and ECLIPSE (3 yr), respectively. We identified 890, 675, and 3,217 genes associated with a history of AE-COPDs, persistent exacerbations (at least one exacerbation per year), and prospective exacerbation rate, respectively. In COPDGene, the number of prospective exacerbations in patients with COPD (Global Initiative for Chronic Obstructive Lung Disease stage ⩾2) was negatively associated with circulating CD8+ T cells, CD4+ T cells, and resting natural killer cells. The negative association with naive CD4+ T cells was replicated in ECLIPSE. In the flow-cytometry study, an increase in CTLA4 on CD4+ T cells was positively associated with AE-COPDs. Conclusions: Individuals with COPD with lower circulating lymphocyte counts, particularly decreased CD4+ T cells, are more susceptible to AE-COPDs, including persistent exacerbations

Analysis of genetically driven alternative splicing identifies FBXO38 as a novel COPD susceptibility gene

While many disease-associated single nucleotide polymorphisms (SNPs) are associated with gene expression (expression quantitative trait loci, eQTLs), a large proportion of complex disease genome-wide association study (GWAS) variants are of unknown function. Some of these SNPs may contribute to disease by regulating gene splicing. Here, we investigate whether SNPs that are associated with alternative splicing (splice QTL or sQTL) can identify novel functions for existing GWAS variants or suggest new associated variants in chronic obstructive pulmonary disease (COPD). RNA sequencing was performed on whole blood from 376 subjects from the COPDGene Study. Using linear models, we identified 561,060 unique sQTL SNPs associated with 30,333 splice sites corresponding to 6,419 unique genes. Similarly, 708,928 unique eQTL SNPs involving 15,913 genes were detected at 10% FDR. While there is overlap between sQTLs and eQTLs, 55.3% of sQTLs are not eQTLs. Co-localization analysis revealed that 7 out of 21 loci associated with COPD (p</p

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

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

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

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