Search

Genome-wide association studies (GWAS) have identified many genomic loci associated with risk for schizophrenia, but unambiguous identification of the relationship between disease-associated variants and specific genes, and in particular their effect on risk conferring transcripts, has proven difficult. To better understand the specific molecular mechanism(s) at the schizophrenia locus in 11q25, we undertook cis expression quantitative trait loci (cis-eQTL) mapping for this 2 megabase genomic region using postmortem human brain samples. To comprehensively assess the effects of genetic risk upon local expression, we evaluated multiple transcript features: genes, exons, and exon-exon junctions in multiple brain regions-dorsolateral prefrontal cortex (DLPFC), hippocampus, and caudate. Genetic risk variants strongly associated with expression of SNX19 transcript features that tag multiple rare classes of SNX19 transcripts, whereas they only weakly affected expression of an exon-exon junction that tags the majority of abundant transcripts. The most prominent class of SNX19 risk-associated transcripts is predicted to be overexpressed, defined by an exon-exon splice junction between exons 8 and 10 (junc8.10) and that is predicted to encode proteins that lack the characteristic nexin C terminal domain. Risk alleles were also associated with either increased or decreased expression of multiple additional classes of transcripts. With RACE, molecular cloning, and long read sequencing, we found a number of novel SNX19 transcripts that further define the set of potential etiological transcripts. We explored epigenetic regulation of SNX19 expression and found that DNA methylation at CpG sites near the primary transcription start site and within exon 2 partially mediate the effects of risk variants on risk-associated expression. ATAC sequencing revealed that some of the most strongly risk-associated SNPs are located within a region of open chromatin, suggesting a nearby regulatory element is involved. These findings indicate a potentially complex molecular etiology, in which risk alleles for schizophrenia generate epigenetic alterations and dysregulation of multiple classes of SNX19 transcripts.

RNA sequencing (RNA-seq) is a powerful approach for measuring gene expression levels in cells and tissues, but it relies on high-quality RNA. We demonstrate here that statistical adjustment using existing quality measures largely fails to remove the effects of RNA degradation when RNA quality associates with the outcome of interest. Using RNA-seq data from molecular degradation experiments of human primary tissues, we introduce a method—quality surrogate variable analysis (qSVA)—as a framework for estimating and removing the confounding effect of RNA quality in differential expression analysis. We show that this approach results in greatly improved replication rates (>3×) across two large independent postmortem human brain studies of schizophrenia and also removes potential RNA quality biases in earlier published work that compared expression levels of different brain regions and other diagnostic groups. Our approach can therefore improve the interpretation of differential expression analysis of transcriptomic data from human tissue.

Histone modification including H3 lysine 79 methylation (H3K79me) plays a key role during gene transcription and DNA damage repair. DOT1L, the sole methyltransferase for three states of H3K79me, is implicated in leukemia, co-lorectal cancer, and dilated cardiomyopathy. However, understanding of DOT1L and H3K79me in these pathways and disease pathogenesis has been limited due to the difficulty of working with DOT1L protein. For instance, locus-specific or genome-wide binding sites of DOT1L revealed by chromatin immunoprecipitation (ChIP)-based methods are necessary for inferring its functions, but high-quality ChIP-grade antibodies are currently not available. Herein we have developed a knock-in approach to tag endogenous DOT1L with 3 × Flag at its C-terminal domain to follow functional analyses. The knock-in was facilitated by using TALENs to induce a targeted double-strand break at the endogenous DOTIL to stimulate local homologous recombination at that site. The single cell colonies with successful knock-in were isolated and verified by different methods. We also demonstrated that tagged DOT1L maintains its normal function in terms of methylation and that the engineered cells would be very useful for further studies.

Genetic disorders have been a major concern for public health in China, especially in the rural regions. However, there is little information available about prevalence of many common single-gene disorders in Guizhou Province in the south western part of China. In the present study, we performed a molecular newborn screening for four genetic disorders, including beta-thalassemia (β-thal), glucose-6-phosphate dehydrogenase (G6PD) deficiency, phenylketonuria (PKU), and non-syndromic hearing loss and deafness (NSHL) in this region. A total of 515 newborns were genotyped using matrix-assisted laser desorption/ionization time of flight mass spectrometry (MALDI-TOF MS) developed for screening the mutations causing these four disorders, and then confirmed by Sanger sequencing. The results showed that 48 out of 515 newborns were carriers of mutations related to these four diseases, with a frequency of 1 in 11 (9.32%). The carrier frequencies for each disease are: β-thal 2.72%; G6PD deficiency 1.94%; PKU 0.78% and NSHL 4.47%. The genotyping results by MALDI-TOF MS were concordant with Sanger sequencing results within 30 randomly selected samples. This is the first study that reveals carrier frequencies of these four diseases in Guizhou Province. These data provide valuable information for the genetic counseling and disease prevention in Guizhou and southwest China.

We have characterized the landscape of DNA methylation (DNAm) across the first two decades of human neocortical development in NeuN+ neurons using whole-genome bisulfite sequencing and compared them to non-neurons (primarily glia) and prenatal homogenate cortex. We show that DNAm changes more dramatically during the first five years of postnatal life than during the entire remaining period. We further refined global patterns of increasingly divergent neuronal CpG and CpH methylation (mCpG and mCpH) into six developmental trajectories and found that in contrast to genome-wide patterns, neighboring mCpG and mCpH levels within these regions were highly correlated. We then integrated paired RNA-seq data and identified direct regulation of hundreds of transcripts and their splicing events exclusively by mCpH levels, independently from mCpG levels, across this period. We finally explored the relationship between DNAm patterns and development of brain-related phenotypes and found enriched heritability for many phenotypes within DNAm features we identify.

Differentiating pluripotent cells from fibroblast progenitors is a potentially transformative tool in personalized medicine. We previously identified relatively greater success culturing dura-derived fibroblasts than scalp-derived fibroblasts from postmortem tissue. We hypothesized that these differences in culture success were related to epigenetic differences between the cultured fibroblasts by sampling location, and therefore generated genome-wide DNA methylation and transcriptome data on 11 intrinsically matched pairs of dural and scalp fibroblasts from donors across the lifespan (infant to 85 years). While these cultured fibroblasts were several generations removed from the primary tissue and morphologically indistinguishable, we found widespread epigenetic differences by sampling location at the single CpG (N = 101,989), region (N = 697), “block” (N = 243), and global spatial scales suggesting a strong epigenetic memory of original fibroblast location. Furthermore, many of these epigenetic differences manifested in the transcriptome, particularly at the region-level. We further identified 7,265 CpGs and 11 regions showing significant epigenetic memory related to the age of the donor, as well as an overall increased epigenetic variability, preferentially in scalp-derived fibroblasts—83% of loci were more variable in scalp, hypothesized to result from cumulative exposure to environmental stimuli in the primary tissue. By integrating publicly available DNA methylation datasets on individual cell populations in blood and brain, we identified significantly increased inter-individual variability in our scalp- and other skin-derived fibroblasts on a similar scale as epigenetic differences between different lineages of blood cells. Lastly, these epigenetic differences did not appear to be driven by somatic mutation—while we identified 64 probable de-novo variants across the 11 subjects, there was no association between mutation burden and age of the donor (p = 0.71). These results depict a strong component of epigenetic memory in cell culture from primary tissue, even after several generations of daughter cells, related to cell state and donor age., Author Summary Regenerative medicine specialists have been using a type of cell commonly found in the skin called the fibroblast because it is easily obtained from skin samples, grows well in culture, and can be manipulated in the laboratory to de-differentiate into a primordial state known as the induced pluripotent stem cell. These primitive stem cells can then be transformed into mature tissues, such as liver or pancreas cells. Here we show that fibroblasts, coming from different locations in the same individual, vary significantly in epigenetic marks called DNA methylation, which are involved in the regulation of gene expression. In addition to location-specific patterns of DNA methylation, we also find that fibroblasts from different anatomical locations respond differently in epigenetic patterns related to aging. As the field of regenerative medicine advances, our study demonstrates that deciding upon the source of fibroblasts from an individual to generate new tissues and organs may be an important consideration.

The purpose of this study was to evaluate the effects of butyrate infusion on rumen epithelial transcriptome. Next-generation sequencing (NGS) and bioinformatics are used to accelerate our understanding of regulation in rumen epithelial transcriptome of cattle in the dry period induced by butyrate infusion at the level of the whole transcriptome. Butyrate, as an essential element of nutrients, is a histone deacetylase (HDAC) inhibitor that can alter histone acetylation and methylation, and plays a prominent role in regulating genomic activities influencing rumen nutrition utilization and function. Ruminal infusion of butyrate was following 0-hour sampling (baseline controls) and continued for 168 hours at a rate of 5.0 L/day of a 2.5 M solution as a continuous infusion. Following the 168-hour infusion, the infusion was stopped, and cows were maintained on the basal lactation ration for an additional 168 hours for sampling. Rumen epithelial samples were serially collected via biopsy through rumen fistulae at 0-, 24-, 72-, and 168-hour (D1, D3, D7) and 168-hour post-infusion (D14). In comparison with pre-infusion at 0 hours, a total of 3513 genes were identified to be impacted in the rumen epithelium by butyrate infusion at least once at different sampling time points at a stringent cutoff of false discovery rate (FDR)

Imprinted genes are vulnerable to environmental influences during early embryonic development, thereby contributing to the onset of disease in adulthood. Monoallelic methylation at several germline imprints has been reported as DNMT1-dependent. However, which of these two epigenetic attributes, DNMT1-dependence or allelic methylation, renders imprinted genes susceptible to environmental stressors has not been determined. Herein, we developed a new approach, referred to as NORED, to identify 2468 DNMT1-dependent DNA methylation patterns in the mouse genome. We further developed an algorithm based on a genetic variation-independent approach (referred to as MethylMosaic) to detect 2487 regions with bimodal methylation patterns. Two approaches identified 207 regions, including known imprinted germline allele-specific methylation patterns (ASMs), that were both NORED and MethylMosaic regions. Examination of methylation in four independent mouse embryonic stem cell lines shows that two regions identified by both NORED and MethylMosaic (Hcn2 and Park7) did not display parent-of-origin-dependent allelic methylation. In these four F1 hybrid cell lines, genetic variation in Cast allele at Hcn2 locus introduces a transcription factor binding site for MTF-1 that may predispose Cast allelic hypomethylation in a reciprocal cross with either C57 or 129 strains. In contrast, each allele of Hcn2 ASM in J1 inbred cell line and Park7 ASM in four F1 hybrid cell lines seems to exhibit similar propensity to be either hypo- or hypermethylated, suggesting a ‘random, switchable’ ASM. Together with published results, our data on ASMs prompted us to propose a hypothesis of regional ‘autosomal chromosome inactivation (ACI)’ that may control a subset of autosomal genes. Therefore, our results open a new avenue to understand monoallelic methylation and provide a rich resource of candidate genes to examine in environmental and nutritional exposure models.

Background High-throughput DNA methylation arrays are likely to accelerate the pace of methylation biomarker discovery for a wide variety of diseases. A potential problem with a standard set of probes measuring the methylation status of CpG sites across the whole genome is that many sites may not show inter-individual methylation variation among the biosamples for the disease outcome being studied. Inclusion of these so-called "non-variable sites" will increase the risk of false discoveries and reduce statistical power to detect biologically relevant methylation markers. Results We propose a method to estimate the proportion of non-variable CpG sites and eliminate those sites from further analyses. Our method is illustrated using data obtained by hybridizing DNA extracted from the peripheral blood mononuclear cells of 311 samples to an array assaying 1505 CpG sites. Results showed that a large proportion of the CpG sites did not show inter-individual variation in methylation. Conclusions Our method resulted in a substantial improvement in association signals between methylation sites and outcome variables while controlling the false discovery rate at the same level.

We have entirely sequenced YCR59, which is a 10·1 kb segment of the right arm of chromosome III, and is a part of the clone E5F from the Newlon collection. The segment contains two long open reading frames (ORFs): YCR591 which starts in the adjacent fragment H9G (situated towards CRY1 and the centromere), and continues with 1833 codons in YCR59. The second ORF YCR592 is 1226 codons long and encoded entirely within YCR59. The two ORFs represent 91% of the total length of the segment. Excellent agreement in both location and length is found between the ORFs YCR591 and YCR592 and the transcripts 86 and 87 respectively in the Yoshikawa and Isono (1990) map of chromosome III. The two ORFs correspond to new genes and show no significant similarity with any known genes.

We have adapted a LacZ promoter trap screen developed by Burns et al.(1994) to search for genes whose expression is dependent on Rtg2p, a protein with an N-terminal hsp70/actin/sugar kinase ATP binding domain. Rtg2p acts upstream of the basic helix–loop–helix/leucine zipper transcription factors, Rtg1p and Rtg3p. All three proteins are known to be required for the expression of the CIT2 gene, which encodes a peroxisomal isoform of citrate synthase whose expression is also dependent on the functional state of mitochondria. Using this screen, we have identified a previously uncharacterized gene, YEL071w, predicted to encode a protein of 496 amino acids that shares 80% homology and 60% sequence identity with actin interacting protein 2, encoded by the AIP2 gene; both proteins also share sequence similarity to aD-lactate dehydrogenase encoded by the DLD1 gene. Expression of YEL071w is dependent on the functional state of mitochondria and on all three of the Rtg proteins, whereas AIP2 expression is independent of the Rtg proteins and the functional state of mitochondria. Like CIT2, the 5′ flanking region of YEL071w contains two R box binding sites for the Rtg1p/Rtg3p heterodimeric transcription complex. Both R boxes are necessary for full YEL071w expression. We show that YEL071w and AIP2 encode proteins withD-lactate dehydrogenase activity, the former located in the cytoplasm and the latter in the mitochondrial matrix. Our data thus provide gene assignments for two previously unrecognizedD-lactate dehydrogenase activities in yeast. Copyright © 1999 John Wiley & Sons, Ltd.

The Saccharomyces Genome Database (SGD) provides Internet access to the complete Saccharomyces cerevisiae genomic sequence, its genes and their products, the phenotypes of its mutants, and the literature supporting these data. The amount of information and the number of features provided by SGD have increased greatly following the release of the S.cerevisiae genomic sequence, which is currently the only complete sequence of a eukaryotic genome. SGD aids researchers by providing not only basic information, but also tools such as sequence similarity searching that lead to detailed information about features of the genome and relationships between genes. SGD presents information using a variety of user-friendly, dynamically created graphical displays illustrating physical, genetic and sequence feature maps. SGD can be accessed via the World Wide Web at http://genome-www.stanford.edu/Saccharomyces/