Intergenerational Inheritance of Hepatic Steatosis in a Mouse Model of Childhood Obesity: Potential Involvement of Germ-Line microRNAs.

Childhood obesity increases the risk of developing metabolic syndrome later in life. Moreover, metabolic dysfunction may be inherited into the following generation through non-genomic mechanisms, with epigenetics as a plausible candidate. The pathways involved in the development of metabolic dysfunction across generations in the context of childhood obesity remain largely unexplored. We have developed a mouse model of early adiposity by reducing litter size at birth (small litter group, SL: 4 pups/dam; control group, C: 8 pups/dam). Mice raised in small litters (SL) developed obesity, insulin resistance and hepatic steatosis with aging. Strikingly, the offspring of SL males (SL-F1) also developed hepatic steatosis. Paternal transmission of an environmentally induced phenotype strongly suggests epigenetic inheritance. We analyzed the hepatic transcriptome in C-F1 and SL-F1 mice to identify pathways involved in the development of hepatic steatosis. We found that the circadian rhythm and lipid metabolic process were the ontologies with highest significance in the liver of SL-F1 mice. We explored whether DNA methylation and small non-coding RNAs might be involved in mediating intergenerational effects. Sperm DNA methylation was largely altered in SL mice. However, these changes did not correlate with the hepatic transcriptome. Next, we analyzed small non-coding RNA content in the testes of mice from the parental generation. Two miRNAs (miR-457 and miR-201) appeared differentially expressed in the testes of SL-F0 mice. They are known to be expressed in mature spermatozoa, but not in oocytes nor early embryos, and they may regulate the transcription of lipogenic genes, but not clock genes, in hepatocytes. Hence, they are strong candidates to mediate the inheritance of adult hepatic steatosis in our murine model. In conclusion, litter size reduction leads to intergenerational effects through non-genomic mechanisms. In our model, DNA methylation does not seem to play a role on the circadian rhythm nor lipid genes. However, at least two paternal miRNAs might influence the expression of a few lipid-related genes in the first-generation offspring, F1.

Mechanisms of epigenetic inheritance of variable traits through the germline.

During the past half century, evidence for inheritance of variable traits has accumulated from experiments in plants and animals and epidemiological studies in humans. Here, we summarize some of the reported cases of epigenetic inheritance and the proposed mechanisms involved in the transmission of non-genetic information between generations in plants, nematodes, flies and mammals. It has long been accepted that information is epigenetically inherited in plants. Although many questions regarding the underlying mechanisms remain to be answered, it is now evident that epigenetic mechanisms are also responsible for the transmission of phenotypes in animals. We highlight similarities and differences between models and species.

Human sperm displays rapid responses to diet.

The global rise in obesity and steady decline in sperm quality are two alarming trends that have emerged during recent decades. In parallel, evidence from model organisms shows that paternal diet can affect offspring metabolic health in a process involving sperm tRNA-derived small RNA (tsRNA). Here, we report that human sperm are acutely sensitive to nutrient flux, both in terms of sperm motility and changes in sperm tsRNA. Over the course of a 2-week diet intervention, in which we first introduced a healthy diet followed by a diet rich in sugar, sperm motility increased and stabilized at high levels. Small RNA-seq on repeatedly sampled sperm from the same individuals revealed that tsRNAs were up-regulated by eating a high-sugar diet for just 1 week. Unsupervised clustering identified two independent pathways for the biogenesis of these tsRNAs: one involving a novel class of fragments with specific cleavage in the T-loop of mature nuclear tRNAs and the other exclusively involving mitochondrial tsRNAs. Mitochondrial involvement was further supported by a similar up-regulation of mitochondrial rRNA-derived small RNA (rsRNA). Notably, the changes in sugar-sensitive tsRNA were positively associated with simultaneous changes in sperm motility and negatively associated with obesity in an independent clinical cohort. This rapid response to a dietary intervention on tsRNA in human sperm is attuned with the paternal intergenerational metabolic responses found in model organisms. More importantly, our findings suggest shared diet-sensitive mechanisms between sperm motility and the biogenesis of tsRNA, which provide novel insights about the interplay between nutrition and male reproductive health.

Assessment of kinship detection using RNA-seq data.

Analysis of RNA sequencing (RNA-seq) data from related individuals is widely used in clinical and molecular genetics studies. Prediction of kinship from RNA-seq data would be useful for confirming the expected relationships in family based studies and for highlighting samples from related individuals in case-control or population based studies. Currently, reconstruction of pedigrees is largely based on SNPs or microsatellites, obtained from genotyping arrays, whole genome sequencing and whole exome sequencing. Potential problems with using RNA-seq data for kinship detection are the low proportion of the genome that it covers, the highly skewed coverage of exons of different genes depending on expression level and allele-specific expression. In this study we assess the use of RNA-seq data to detect kinship between individuals, through pairwise identity by descent (IBD) estimates. First, we obtained high quality SNPs after successive filters to minimize the effects due to allelic imbalance as well as errors in sequencing, mapping and genotyping. Then, we used these SNPs to calculate pairwise IBD estimates. By analysing both real and simulated RNA-seq data we show that it is possible to identify up to second degree relationships using RNA-seq data of even low to moderate sequencing depth.

The Polycomb-Dependent Epigenome Controls ß Cell Dysfunction, Dedifferentiation, and Diabetes.

To date, it remains largely unclear to what extent chromatin machinery contributes to the susceptibility and progression of complex diseases. Here, we combine deep epigenome mapping with single-cell transcriptomics to mine for evidence of chromatin dysregulation in type 2 diabetes. We find two chromatin-state signatures that track ß cell dysfunction in mice and humans: ectopic activation of bivalent Polycomb-silenced domains and loss of expression at an epigenomically unique class of lineage-defining genes. ß cell-specific Polycomb (Eed/PRC2) loss of function in mice triggers diabetes-mimicking transcriptional signatures and highly penetrant, hyperglycemia-independent dedifferentiation, indicating that PRC2 dysregulation contributes to disease. The work provides novel resources for exploring ß cell transcriptional regulation and identifies PRC2 as necessary for long-term maintenance of ß cell identity. Importantly, the data suggest a two-hit (chromatin and hyperglycemia) model for loss of ß cell identity in diabetes.

Double deficiency of Trex2 and DNase1L2 nucleases leads to accumulation of DNA in lingual cornifying keratinocytes without activating inflammatory responses.

The cornification of keratinocytes on the surface of skin and oral epithelia is associated with the degradation of nuclear DNA. The endonuclease DNase1L2 and the exonuclease Trex2 are expressed specifically in cornifying keratinocytes. Deletion of DNase1L2 causes retention of nuclear DNA in the tongue epithelium but not in the skin. Here we report that lack of Trex2 results in the accumulation of DNA fragments in the cytoplasm of cornifying lingual keratinocytes and co-deletion of DNase1L2 and Trex2 causes massive accumulation of DNA fragments throughout the cornified layers of the tongue epithelium. By contrast, cornification-associated DNA breakdown was not compromised in the epidermis. Aberrant retention of DNA in the tongue epithelium was associated neither with enhanced expression of DNA-driven response genes, such as Ifnb, Irf7 and Cxcl10, nor with inflammation. Of note, the expression of Tlr9, Aim2 and Tmem173, key DNA sensor genes, was markedly lower in keratinocytes and keratinocyte-built tissues than in macrophages and immune tissues, and DNA-driven response genes were not induced by introduction of DNA in keratinocytes. Altogether, our results indicate that DNase1L2 and Trex2 cooperate in the breakdown and degradation of DNA during cornification of lingual keratinocytes and aberrant DNA retention is tolerated in the oral epithelium.

Impaired DNA replication derepresses chromatin and generates a transgenerationally inherited epigenetic memory.

Impaired DNA replication is a hallmark of cancer and a cause of genomic instability. We report that, in addition to causing genetic change, impaired DNA replication during embryonic development can have major epigenetic consequences for a genome. In a genome-wide screen, we identified impaired DNA replication as a cause of increased expression from a repressed transgene in Caenorhabditis elegans. The acquired expression state behaved as an "epiallele," being inherited for multiple generations before fully resetting. Derepression was not restricted to the transgene but was caused by a global reduction in heterochromatin-associated histone modifications due to the impaired retention of modified histones on DNA during replication in the early embryo. Impaired DNA replication during development can therefore globally derepress chromatin, creating new intergenerationally inherited epigenetic expression states.

Transgenerational transmission of environmental information in C. elegans.

The environment experienced by an animal can sometimes influence gene expression for one or a few subsequent generations. Here, we report the observation that a temperature-induced change in expression from a Caenorhabditis elegans heterochromatic gene array can endure for at least 14 generations. Inheritance is primarily in cis with the locus, occurs through both oocytes and sperm, and is associated with altered trimethylation of histone H3 lysine 9 (H3K9me3) before the onset of zygotic transcription. Expression profiling reveals that temperature-induced expression from endogenous repressed repeats can also be inherited for multiple generations. Long-lasting epigenetic memory of environmental change is therefore possible in this animal.

The epigenetic landscape of Alu repeats delineates the structural and functional genomic architecture of colon cancer cells.

Cancer cells exhibit multiple epigenetic changes with prominent local DNA hypermethylation and widespread hypomethylation affecting large chromosomal domains. Epigenome studies often disregard the study of repeat elements owing to technical complexity and their undefined role in genome regulation. We have developed NSUMA (Next-generation Sequencing of UnMethylated Alu), a cost-effective approach allowing the unambiguous interrogation of DNA methylation in more than 130,000 individual Alu elements, the most abundant retrotransposon in the human genome. DNA methylation profiles of Alu repeats have been analyzed in colon cancers and normal tissues using NSUMA and whole-genome bisulfite sequencing. Normal cells show a low proportion of unmethylated Alu (1%-4%) that may increase up to 10-fold in cancer cells. In normal cells, unmethylated Alu elements tend to locate in the vicinity of functionally rich regions and display epigenetic features consistent with a direct impact on genome regulation. In cancer cells, Alu repeats are more resistant to hypomethylation than other retroelements. Genome segmentation based on high/low rates of Alu hypomethylation allows the identification of genomic compartments with differential genetic, epigenetic, and transcriptomic features. Alu hypomethylated regions show low transcriptional activity, late DNA replication, and its extent is associated with higher chromosomal instability. Our analysis demonstrates that Alu retroelements contribute to define the epigenetic landscape of normal and cancer cells and provides a unique resource on the epigenetic dynamics of a principal, but largely unexplored, component of the primate genome.

The Exonuclease Trex2 Shapes Psoriatic Phenotype.

Trex2 is a keratinocyte-specific 3'-deoxyribonuclease that participates in the maintenance of skin homeostasis after DNA damage. Here, we show that this exonuclease is strongly upregulated in human psoriasis, a hyperproliferative and inflammatory skin disease. Similarly, the imiquimod (IMQ)- and Il23-induced mouse psoriasis was associated with a substantial upregulation of Trex2, which was recruited into fragmented chromatin in keratinocytes that were undergoing impaired proliferation, differentiation, and cell death, indicating an important role in DNA processing. Using Trex2 knockout mice, we have found that Trex2 deficiency attenuated IMQ-induced psoriasis-like skin inflammation and enhanced IMQ-induced parakeratosis. Also, Il23-induced ear swelling was diminished in Trex2 knockout mice in comparison with wild-type (wt) mice. Transcriptome analysis identified multiple genes that were deregulated by Trex2 loss after treatment with IMQ. Specifically, immune response genes and pathways normally associated with inflammation were downregulated, whereas those related to skin differentiation and chromatin biology showed increased expression. Interestingly, Trex2 deficiency led to decreased IMQ-induced keratinocyte death via both cell autonomous and noncell autonomous mechanisms. Hence, our data indicate that Trex2 acts as a critical factor in the pathogenesis of psoriasis by promoting keratinocyte apoptosis and enucleation and thereby influencing skin immune responses.

Trim28 Haploinsufficiency Triggers Bi-stable Epigenetic Obesity.

More than one-half billion people are obese, and despite progress in genetic research, much of the heritability of obesity remains enigmatic. Here, we identify a Trim28-dependent network capable of triggering obesity in a non-Mendelian, "on/off" manner. Trim28(+/D9) mutant mice exhibit a bi-modal body-weight distribution, with isogenic animals randomly emerging as either normal or obese and few intermediates. We find that the obese-"on" state is characterized by reduced expression of an imprinted gene network including Nnat, Peg3, Cdkn1c, and Plagl1 and that independent targeting of these alleles recapitulates the stochastic bi-stable disease phenotype. Adipose tissue transcriptome analyses in children indicate that humans too cluster into distinct sub-populations, stratifying according to Trim28 expression, transcriptome organization, and obesity-associated imprinted gene dysregulation. These data provide evidence of discrete polyphenism in mouse and man and thus carry important implications for complex trait genetics, evolution, and medicine.

The small RNA content of human sperm reveals pseudogene-derived piRNAs complementary to protein-coding genes.

At the end of mammalian sperm development, sperm cells expel most of their cytoplasm and dispose of the majority of their RNA. Yet, hundreds of RNA molecules remain in mature sperm. The biological significance of the vast majority of these molecules is unclear. To better understand the processes that generate sperm small RNAs and what roles they may have, we sequenced and characterized the small RNA content of sperm samples from two human fertile individuals. We detected 182 microRNAs, some of which are highly abundant. The most abundant microRNA in sperm is miR-1246 with predicted targets among sperm-specific genes. The most abundant class of small noncoding RNAs in sperm are PIWI-interacting RNAs (piRNAs). Surprisingly, we found that human sperm cells contain piRNAs processed from pseudogenes. Clusters of piRNAs from human testes contain pseudogenes transcribed in the antisense strand and processed into small RNAs. Several human protein-coding genes contain antisense predicted targets of pseudogene-derived piRNAs in the male germline and these piRNAs are still found in mature sperm. Our study provides the most extensive data set and annotation of human sperm small RNAs to date and is a resource for further functional studies on the roles of sperm small RNAs. In addition, we propose that some of the pseudogene-derived human piRNAs may regulate expression of their parent gene in the male germline.

Paternal diet defines offspring chromatin state and intergenerational obesity.

The global rise in obesity has revitalized a search for genetic and epigenetic factors underlying the disease. We present a Drosophila model of paternal-diet-induced intergenerational metabolic reprogramming (IGMR) and identify genes required for its encoding in offspring. Intriguingly, we find that as little as 2 days of dietary intervention in fathers elicits obesity in offspring. Paternal sugar acts as a physiological suppressor of variegation, desilencing chromatin-state-defined domains in both mature sperm and in offspring embryos. We identify requirements for H3K9/K27me3-dependent reprogramming of metabolic genes in two distinct germline and zygotic windows. Critically, we find evidence that a similar system may regulate obesity susceptibility and phenotype variation in mice and humans. The findings provide insight into the mechanisms underlying intergenerational metabolic reprogramming and carry profound implications for our understanding of phenotypic variation and evolution.

Sperm epigenomics: challenges and opportunities.

Sperm is a highly differentiated cell type whose function is to deliver a haploid genome to the oocyte. The sperm "epigenomes" were traditionally considered to be insignificant - the sperm is transcriptionally inactive, its genome is packaged in sperm-specific protamine toroids instead of nucleosomes, and its DNA methylation profile is erased immediately post-fertilization. Yet, in recent years there has been an increase in the number of reported cases of apparent epigenetic inheritance through the male germline, suggesting that the sperm epigenome may transmit information between generations. At the same time, technical advances have made the genome-wide profiling of different layers of the sperm epigenome feasible. As a result, a large number of datasets have been recently generated and analyzed with the aim to better understand what non-genetic material is contained within the sperm and whether it has any function post-fertilization. Here, we provide an overview of the current knowledge of the sperm epigenomes as well as the challenges in analysing them and the opportunities in understanding the potential non-genetic carriers of information in sperm.

Genomic and proteomic dissection and characterization of the human sperm chromatin.

The mammalian spermatozoon has a unique chromatin structure where the majority of DNA is packaged by protamines, while a small fraction (∼8%) remains associated with nucleosomes. However, the chromatin affinity and repertoire of the additional proteins constituting the different sperm chromatin fractions have not yet been explored. To address this we have carried out a genomic and proteomic characterization of human sperm samples subjected to chromatin fractionation using either 0.65 M NaCl extraction followed by EcoRI/BamHI DNA restriction enzyme digestion, or micrococcal nuclease digestion. DNA fractions corresponding to the nucleosome-packaged DNA were sequenced, confirming an appropriate dissection of the sperm chromatin. In addition we detected and sequenced a subnucleosomal particle. Although both fractions were highly enriched at gene promoters, some sequences were found to be exclusively associated with one of those. The results of the proteomic analyses demonstrate that there are two distinct sets of sperm proteins which differ in chromatin affinity. Histone variants, transcription factors, chromatin-associated and modifying proteins involved in regulatory roles were identified as weakly attached to the sperm DNA, whereas proteins with structural roles were identified in the condensed fraction. Many factors, such as the histone lysine demethylase PHF8 identified for the first time in the human sperm cell in this study, were identified exclusively in soluble fraction. Our results provide additional support to the possibility that all of these factors may constitute additional layers of sperm epigenetic information or have structural or regulatory roles transmitted by the sperm cell to the oocyte at fertilization.

Parallel evolution of chordate cis-regulatory code for development.

Urochordates are the closest relatives of vertebrates and at the larval stage, possess a characteristic bilateral chordate body plan. In vertebrates, the genes that orchestrate embryonic patterning are in part regulated by highly conserved non-coding elements (CNEs), yet these elements have not been identified in urochordate genomes. Consequently the evolution of the cis-regulatory code for urochordate development remains largely uncharacterised. Here, we use genome-wide comparisons between C. intestinalis and C. savignyi to identify putative urochordate cis-regulatory sequences. Ciona conserved non-coding elements (ciCNEs) are associated with largely the same key regulatory genes as vertebrate CNEs. Furthermore, some of the tested ciCNEs are able to activate reporter gene expression in both zebrafish and Ciona embryos, in a pattern that at least partially overlaps that of the gene they associate with, despite the absence of sequence identity. We also show that the ability of a ciCNE to up-regulate gene expression in vertebrate embryos can in some cases be localised to short sub-sequences, suggesting that functional cross-talk may be defined by small regions of ancestral regulatory logic, although functional sub-sequences may also be dispersed across the whole element. We conclude that the structure and organisation of cis-regulatory modules is very different between vertebrates and urochordates, reflecting their separate evolutionary histories. However, functional cross-talk still exists because the same repertoire of transcription factors has likely guided their parallel evolution, exploiting similar sets of binding sites but in different combinations.

Human genes with CpG island promoters have a distinct transcription-associated chromatin organization.

BACKGROUND: More than 50% of human genes initiate transcription from CpG dinucleotide-rich regions referred to as CpG islands. These genes show differences in their patterns of transcription initiation, and have been reported to have higher levels of some activation-associated chromatin modifications. RESULTS: Here we report that genes with CpG island promoters have a characteristic transcription-associated chromatin organization. This signature includes high levels of the transcription elongation-associated histone modifications H4K20me1, H2BK5me1 and H3K79me1/2/3 in the 5' end of the gene, depletion of the activation marks H2AK5ac, H3K14ac and H3K23ac immediately downstream of the transcription start site (TSS), and characteristic epigenetic asymmetries around the TSS. The chromosome organization factor CTCF may be bound upstream of RNA polymerase in most active CpG island promoters, and an unstable nucleosome at the TSS may be specifically marked by H4K20me3, the first example of such a modification. H3K36 monomethylation is only detected as enriched in the bodies of active genes that have CpG island promoters. Finally, as expression levels increase, peak modification levels of the histone methylations H3K9me1, H3K4me1, H3K4me2 and H3K27me1 shift further away from the TSS into the gene body. CONCLUSIONS: These results suggest that active genes with CpG island promoters have a distinct step-like series of modified nucleosomes after the TSS. The identity, positioning, shape and relative ordering of transcription-associated histone modifications differ between genes with and without CpG island promoters. This supports a model where chromatin organization reflects not only transcription activity but also the type of promoter in which transcription initiates.

Chromatin organization in sperm may be the major functional consequence of base composition variation in the human genome.

Chromatin in sperm is different from that in other cells, with most of the genome packaged by protamines not nucleosomes. Nucleosomes are, however, retained at some genomic sites, where they have the potential to transmit paternal epigenetic information. It is not understood how this retention is specified. Here we show that base composition is the major determinant of nucleosome retention in human sperm, predicting retention very well in both genic and non-genic regions of the genome. The retention of nucleosomes at GC-rich sequences with high intrinsic nucleosome affinity accounts for the previously reported retention at transcription start sites and at genes that regulate development. It also means that nucleosomes are retained at the start sites of most housekeeping genes. We also report a striking link between the retention of nucleosomes in sperm and the establishment of DNA methylation-free regions in the early embryo. Taken together, this suggests that paternal nucleosome transmission may facilitate robust gene regulation in the early embryo. We propose that chromatin organization in the male germline, rather than in somatic cells, is the major functional consequence of fine-scale base composition variation in the human genome. The selective pressure driving base composition evolution in mammals could, therefore, be the need to transmit paternal epigenetic information to the zygote.

Predicting genetic modifier loci using functional gene networks.

Most phenotypes are genetically complex, with contributions from mutations in many different genes. Mutations in more than one gene can combine synergistically to cause phenotypic change, and systematic studies in model organisms show that these genetic interactions are pervasive. However, in human association studies such nonadditive genetic interactions are very difficult to identify because of a lack of statistical power--simply put, the number of potential interactions is too vast. One approach to resolve this is to predict candidate modifier interactions between loci, and then to specifically test these for associations with the phenotype. Here, we describe a general method for predicting genetic interactions based on the use of integrated functional gene networks. We show that in both Saccharomyces cerevisiae and Caenorhabditis elegans a single high-coverage, high-quality functional network can successfully predict genetic modifiers for the majority of genes. For C. elegans we also describe the construction of a new, improved, and expanded functional network, WormNet 2. Using this network we demonstrate how it is possible to rapidly expand the number of modifier loci known for a gene, predicting and validating new genetic interactions for each of three signal transduction genes. We propose that this approach, termed network-guided modifier screening, provides a general strategy for predicting genetic interactions. This work thus suggests that a high-quality integrated human gene network will provide a powerful resource for modifier locus discovery in many different diseases.

Intrinsic protein disorder and interaction promiscuity are widely associated with dosage sensitivity.

Why are genes harmful when they are overexpressed? By testing possible causes of overexpression phenotypes in yeast, we identify intrinsic protein disorder as an important determinant of dosage sensitivity. Disordered regions are prone to make promiscuous molecular interactions when their concentration is increased, and we demonstrate that this is the likely cause of pathology when genes are overexpressed. We validate our findings in two animals, Drosophila melanogaster and Caenorhabditis elegans. In mice and humans the same properties are strongly associated with dosage-sensitive oncogenes, such that mass-action-driven molecular interactions may be a frequent cause of cancer. Dosage-sensitive genes are tightly regulated at the transcriptional, RNA, and protein levels, which may serve to prevent harmful increases in protein concentration under physiological conditions. Mass-action-driven interaction promiscuity is a single theoretical framework that can be used to understand, predict, and possibly treat the effects of increased gene expression in evolution and disease.