Genetic and epigenetic variation among inbred mouse littermates: identification of inter-individual differentially methylated regions

Epigenetic Modulation of GPER Expression in Gastric and Colonic Smooth Muscle of Male and Female Non-Obese Diabetic (NOD) Mice: Insights into H3K4me3 and H3K27ac Modifications

Type 1 diabetes (T1D) affects gastrointestinal (GI) motility, favoring gastroparesis, constipation, and fecal incontinence, which are more prevalent in women. The mechanisms are unknown. Given the G-protein-coupled estrogen receptor's (GPER) role in GI motility, we investigated sex-related diabetes-induced epigenetic changes in GPER. We assessed GPER mRNA and protein expression levels using qPCR and Western blot analyses, and quantified the changes in nuclear DNA methyltransferases and histone modifications (H3K4me3, H3Ac, and H3K27Ac) by ELISA kits. Targeted bisulfite and chromatin immunoprecipitation assays were used to evaluate DNA methylation and histone modifications around the GPER promoter by chromatin immunoprecipitation assays in gastric and colonic smooth muscle tissues of male and female control (CTR) and non-obese diabetic (NOD) mice. GPER expression was downregulated in NOD, with sex-dependent variations. In the gastric smooth muscle, not in colonic smooth muscle, downregulation coincided with differences in methylation ratios between regions 1 and 2 of the GPER promoter of NOD. DNA methylation was higher in NOD male colonic smooth muscle than in NOD females. H3K4me3 and H3ac enrichment decreased in NOD gastric smooth muscle. H3K4me3 levels diminished in the colonic smooth muscle of NOD. H3K27ac levels were unaffected, but enrichment decreased in NOD male gastric smooth muscle; however, it increased in the NOD male colonic smooth muscle and decreased in the female NOD colonic smooth muscle. Male NOD colonic smooth muscle exhibited decreased H3K27ac levels, not female, whereas female NOD colonic smooth muscle demonstrated diminished enrichment of H3ac at the GPER promoter, contrary to male NOD. Sex-specific epigenetic mechanisms contribute to T1D-mediated suppression of GPER expression in the GI tract. These insights advance our understanding of T1D complications and suggest promising avenues for targeted therapeutic interventions.

Metastable epialleles in humans

First identified in isogenic mice, metastable epialleles (MEs) are loci where the extent of DNA methylation (DNAm) is variable between individuals but correlates across tissues derived from different germ layers within a given individual. This property, termed systemic interindividual variation (SIV), is attributed to stochastic methylation establishment before germ layer differentiation. Evidence suggests that some putative human MEs are sensitive to environmental exposures in early development. In this review we introduce key concepts pertaining to human MEs, describe methods used to identify MEs in humans, and review their genomic features. We also highlight studies linking DNAm at putative human MEs to early environmental exposures and postnatal (including disease) phenotypes.

Between therapy effect and false-positive result in animal experimentation

Despite the animal models' complexity, researchers tend to reduce the number of animals in experiments for expenses and ethical concerns. This tendency makes the risk of false-positive results, as statistical significance, the primary criterion to validate findings, often fails if testing small samples. This study aims to highlight such risks using an example from experimental regenerative therapy and propose a machine-learning solution to validate treatment effects. The example analysed was the pharmacological treatment of ear pinna punch wound healing in mice. Wound closure data analysed included eight groups treated with an epigenetic inhibitor, zebularine, and eight control groups receiving vehicle alone, of six mice each. We confirmed the zebularine healing effect for all 64 pairwise comparisons between treatment and control groups but also determined minor yet statistically significant differences between control groups in five of 28 possible comparisons. The occurrences of significant differences between the control groups, regardless of standardised experimental conditions, indicate a risk of statistically significant effects in the case a compound lacking the desired biological activity is tested. Since the criterion of statistical significance itself can be confusing, we demonstrate a machine-learning algorithm trained on datasets representing treatment and control experiments as a helpful tool for validating treatment outcomes. We tested two machine-learning approaches, Naïve Bayes and Support Vector Machine classifiers. In contrast to the Mann-Whitney U-test, indicating enhanced healing effects for some control groups receiving saline alone, both machine-learning algorithms faultlessly assigned all animal groups receiving saline to the controls.

Epigenomic reprogramming in iAs-mediated carcinogenesis

Arsenic is a naturally occurring metal carcinogen found in the Earth's crust. Millions of people worldwide are chronically exposed to arsenic through drinking water and food. Exposure to inorganic arsenic has been implicated in many diseases ranging from acute toxicities to malignant transformations. Despite the well-known deleterious health effects of arsenic exposure, the molecular mechanisms in arsenic-mediated carcinogenesis are not fully understood. Since arsenic is non-mutagenic, the mechanism by which arsenic causes carcinogenesis is via alterations in epigenetic-regulated gene expression. There are two possible ways by which arsenic may modify the epigenome-indirectly through an arsenic-induced generation of reactive oxygen species which then impacts chromatin remodelers, or directly through interaction and modulation of chromatin remodelers. Whether directly or indirectly, arsenic modulates epigenetic gene regulation and our understanding of the direct effect of this modulation on chromatin structure is limited. In this chapter we will discuss the various ways by which inorganic arsenic affects the epigenome with consequences in health and disease.

Interfering With Contextual Fear Memories by Post-reactivation Administration of Propranolol in Mice: A Series of Null Findings

Post-reactivation amnesia of contextual fear memories by blockade of noradrenergic signaling has been shown to have limited replicability in rodents. This is usually attributed to several boundary conditions that gate the destabilization of memory during its retrieval. How these boundary conditions can be overcome, and what neural mechanisms underlie post-reactivation changes in contextual fear memories remain largely unknown. Here, we report a series of experiments in a contextual fear-conditioning paradigm in mice, that were aimed at solving these issues. We first attempted to obtain a training paradigm that would consistently result in contextual fear memory that could be destabilized upon reactivation, enabling post-retrieval amnesia by the administration of propranolol. Unexpectedly, our attempts were unsuccessful to this end. Specifically, over a series of experiments in which we varied different parameters of the fear acquisition procedure, at best small and inconsistent effects were observed. Additionally, we found that propranolol did not alter retrieval-induced neural activity, as measured by the number of c-Fos⁺ cells in the hippocampal dentate gyrus. To determine whether propranolol was perhaps ineffective in interfering with reactivated contextual fear memories, we also included anisomycin (i.e., a potent and well-known amnesic drug) in several experiments, and measures of synaptic glutamate receptor subunit GluA2 (i.e., a marker of memory destabilization). No post-retrieval amnesia by anisomycin and no altered GluA2 expression by reactivation was observed, suggesting that the memories did not undergo destabilization. The null findings are surprising, given that the training paradigms we implemented were previously shown to result in memories that could be modified upon reactivation. Together, our observations illustrate the elusive nature of reactivation-dependent changes in non-human fear memory.

Faster skin wound healing predicts survival after myocardial infarction

Both skin wound healing and the cardiac response to myocardial infarction (MI) progress through similar pathways involving inflammation, resolution, tissue repair, and scar formation. Due to the similarities, we hypothesized that the healing response to skin wounding would predict future response to MI. Mice were given a 3-mm skin wound using a disposable biopsy punch and the skin wound was imaged daily until closure. The same set of animals was given MI by permanent coronary artery ligation 28 days later and followed for 7 days. Cardiac physiology was measured by echocardiography at baseline and MI days 3 and 7. Animals that survived until day 7 were grouped as survivors, and animals that died from MI were grouped as nonsurvivors. Survivors had faster skin wound healing than nonsurvivors. Faster skin wound healing predicted MI survival better than commonly used cardiac functional variables (e.g., infarct size, fractional shortening, and end diastolic dimension). N-glycoproteome profiling of MI day 3 plasma revealed α2-macroglobulin and ELL-associated factor 1 as strong predictors of future MI death and progression to heart failure. A second cohort of MI mice validated these findings. To investigate the clinical relevance of α2-macroglobulin, we mapped the plasma glycoproteome in patients with MI 48 h after admission and in healthy controls. In patients, α2-macroglobulin was increased 48 h after MI. Apolipoprotein D, another plasma glycoprotein, detrimentally regulated both skin and cardiac wound healing in male but not female mice by promoting inflammation. Our results reveal that the skin is a mirror to the heart and common pathways link wound healing across organs.NEW & NOTEWORTHY Faster skin wound healers had more efficient cardiac healing after myocardial infarction (MI). Two plasma proteins at D3 MI, EAF1 and A2M, predicted MI death in 66% of cases. ApoD regulated both skin and cardiac wound healing in male mice by promoting inflammation. The skin was a mirror to the heart and common pathways linked wound healing across organs.

Transmission of trained immunity and heterologous resistance to infections across generations

Intergenerational inheritance of immune traits linked to epigenetic modifications has been demonstrated in plants and invertebrates. Here we provide evidence for transmission of trained immunity across generations to murine progeny that survived a sublethal systemic infection with Candida albicans or a zymosan challenge. The progeny of trained mice exhibited cellular, developmental, transcriptional and epigenetic changes associated with the bone marrow-resident myeloid effector and progenitor cell compartment. Moreover, the progeny of trained mice showed enhanced responsiveness to endotoxin challenge, alongside improved protection against systemic heterologous Escherichia coli and Listeria monocytogenes infections. Sperm DNA of parental male mice intravenously infected with the fungus C. albicans showed DNA methylation differences linked to immune gene loci. These results provide evidence for inheritance of trained immunity in mammals, enhancing protection against infections.

Sequence features of retrotransposons allow for epigenetic variability

Transposable elements (TEs) are mobile genetic elements that make up a large fraction of mammalian genomes. While select TEs have been co-opted in host genomes to have function, the majority of these elements are epigenetically silenced by DNA methylation in somatic cells. However, some TEs in mice, including the Intracisternal A-particle (IAP) subfamily of retrotransposons, have been shown to display interindividual variation in DNA methylation. Recent work has revealed that IAP sequence differences and strain-specific KRAB zinc finger proteins (KZFPs) may influence the methylation state of these IAPs. However, the mechanisms underlying the establishment and maintenance of interindividual variability in DNA methylation still remain unclear. Here, we report that sequence content and genomic context influence the likelihood that IAPs become variably methylated. IAPs that differ from consensus IAP sequences have altered KZFP recruitment that can lead to decreased KAP1 recruitment when in proximity of constitutively expressed genes. These variably methylated loci have a high CpG density, similar to CpG islands, and can be bound by ZF-CxxC proteins, providing a potential mechanism to maintain this permissive chromatin environment and protect from DNA methylation. These observations indicate that variably methylated IAPs escape silencing through both attenuation of KZFP binding and recognition by ZF-CxxC proteins to maintain a hypomethylated state.

Sex Differences in the Relation Between Frailty and Endothelial Dysfunction in Old Mice

Vascular endothelial function declines with age on average, but there is high variability in the magnitude of this decline within populations. Measurements of frailty, known as frailty index (FI), can be used as surrogates for biological age, but it is unknown if frailty relates to the age-related decline in vascular function. To examine this relation, we studied young (4-9 months) and old (23-32 months) C57BL6 mice of both sexes. We found that FI was greater in old compared with young mice, but did not differ between old male and female mice. Middle cerebral artery (MCA) and mesenteric artery endothelium-dependent dilation (EDD) also did not differ between old male and female mice; however, there were sex differences in the relations between FI and EDD. For the MCA, FI was inversely related to EDD among old female mice, but not old male mice. In contrast, for the mesenteric artery, FI was inversely related to EDD among old male mice, but not old female mice. A higher FI was related to a greater improvement in EDD with the superoxide scavenger TEMPOL in the MCAs for old female mice and in the mesenteric arteries for old male mice. FI related to mesenteric artery gene expression negatively for extracellular superoxide dismutase (Sod3) and positively for interleukin-1β (Il1b). In summary, we found that the relation between frailty and endothelial function is dependent on sex and the artery examined. Arterial oxidative stress and pro-inflammatory signaling are potential mediators of the relations of frailty and endothelial function.

Activation of a hippocampal CREB-pCREB-miRNA-MEF2 axis modulates individual variation of spatial learning and memory capability

Phenotypic variation is a fundamental prerequisite for cell and organism evolution by natural selection. Whereas the role of stochastic gene expression in phenotypic diversity of genetically identical cells is well studied, not much is known regarding the relationship between stochastic gene expression and individual behavioral variation in animals. We demonstrate that a specific miRNA (miR-466f-3p) is upregulated in the hippocampus of a portion of individual inbred mice upon a Morris water maze task. Significantly, miR-466f-3p positively regulates the neuron morphology, function and spatial learning, and memory capability of mice. Mechanistically, miR-466f-3p represses translation of MEF2A, a negative regulator of learning/memory. Finally, we show that varied upregulation of hippocampal miR-466f-3p results from randomized phosphorylation of hippocampal cyclic AMP (cAMP)-response element binding (CREB) in individuals. This finding of modulation of spatial learning and memory via a randomized hippocampal signaling axis upon neuronal stimulation represents a demonstration of how variation in tissue gene expression lead to varied animal behavior.

Defective folate metabolism causes germline epigenetic instability and distinguishes Hira as a phenotype inheritance biomarker

The mechanism behind transgenerational epigenetic inheritance is unclear, particularly through the maternal grandparental line. We previously showed that disruption of folate metabolism in mice by the Mtrr hypomorphic mutation results in transgenerational epigenetic inheritance of congenital malformations. Either maternal grandparent can initiate this phenomenon, which persists for at least four wildtype generations. Here, we use genome-wide approaches to reveal genetic stability in the Mtrr model and genome-wide differential DNA methylation in the germline of Mtrr mutant maternal grandfathers. We observe that, while epigenetic reprogramming occurs, wildtype grandprogeny and great grandprogeny exhibit transcriptional changes that correlate with germline methylation defects. One region encompasses the Hira gene, which is misexpressed in embryos for at least three wildtype generations in a manner that distinguishes Hira transcript expression as a biomarker of maternal phenotypic inheritance.

CRISPRi screens reveal a DNA methylation-mediated 3D genome dependent causal mechanism in prostate cancer

Prostate cancer (PCa) risk-associated SNPs are enriched in noncoding cis-regulatory elements (rCREs), yet their modi operandi and clinical impact remain elusive. Here, we perform CRISPRi screens of 260 rCREs in PCa cell lines. We find that rCREs harboring high risk SNPs are more essential for cell proliferation and H3K27ac occupancy is a strong indicator of essentiality. We also show that cell-line-specific essential rCREs are enriched in the 8q24.21 region, with the rs11986220-containing rCRE regulating MYC and PVT1 expression, cell proliferation and tumorigenesis in a cell-line-specific manner, depending on DNA methylation-orchestrated occupancy of a CTCF binding site in between this rCRE and the MYC promoter. We demonstrate that CTCF deposition at this site as measured by DNA methylation level is highly variable in prostate specimens, and observe the MYC eQTL in the 8q24.21 locus in individuals with low CTCF binding. Together our findings highlight a causal mechanism synergistically driven by a risk SNP and DNA methylation-mediated 3D genome architecture, advocating for the integration of genetics and epigenetics in assessing risks conferred by genetic predispositions.

PHACTR1 genetic variability is not critical in small vessel ischemic disease patients and PcomA recruitment in C57BL/6J mice

Recently, several genome-wide association studies identified PHACTR1 as key locus for five diverse vascular disorders: coronary artery disease, migraine, fibromuscular dysplasia, cervical artery dissection and hypertension. Although these represent significant risk factors or comorbidities for ischemic stroke, PHACTR1 role in brain small vessel ischemic disease and ischemic stroke most important survival mechanism, such as the recruitment of brain collateral arteries like posterior communicating arteries (PcomAs), remains unknown. Therefore, we applied exome and genome sequencing in a multi-ethnic cohort of 180 early-onset independent familial and apparently sporadic brain small vessel ischemic disease and CADASIL-like Caucasian patients from US, Portugal, Finland, Serbia and Turkey and in 2 C57BL/6J stroke mouse models (bilateral common carotid artery stenosis [BCCAS] and middle cerebral artery occlusion [MCAO]), characterized by different degrees of PcomAs patency. We report 3 very rare coding variants in the small vessel ischemic disease-CADASIL-like cohort (p.Glu198Gln, p.Arg204Gly, p.Val251Leu) and a stop-gain mutation (p.Gln273*) in one MCAO mouse. These coding variants do not cluster in PHACTR1 known pathogenic domains and are not likely to play a critical role in small vessel ischemic disease or brain collateral circulation. We also exclude the possibility that copy number variants (CNVs) or a variant enrichment in Phactr1 may be associated with PcomA recruitment in BCCAS mice or linked to diverse vascular traits (cerebral blood flow pre-surgery, PcomA size, leptomeningeal microcollateral length and junction density during brain hypoperfusion) in C57BL/6J mice, respectively. Genetic variability in PHACTR1 is not likely to be a common susceptibility factor influencing small vessel ischemic disease in patients and PcomA recruitment in C57BL/6J mice. Nonetheless, rare variants in PHACTR1 RPEL domains may influence the stroke outcome and are worth investigating in a larger cohort of small vessel ischemic disease patients, different ischemic stroke subtypes and with functional studies.

Genomic properties of variably methylated retrotransposons in mouse

BACKGROUND

Transposable elements (TEs) are enriched in cytosine methylation, preventing their mobility within the genome. We previously identified a genome-wide repertoire of candidate intracisternal A particle (IAP) TEs in mice that exhibit inter-individual variability in this methylation (VM-IAPs) with implications for genome function.

RESULTS

Here we validate these metastable epialleles and discover a novel class that exhibit tissue specificity (tsVM-IAPs) in addition to those with uniform methylation in all tissues (constitutive- or cVM-IAPs); both types have the potential to regulate genes in cis. Screening for variable methylation at other TEs shows that this phenomenon is largely limited to IAPs, which are amongst the youngest and most active endogenous retroviruses. We identify sequences enriched within cVM-IAPs, but determine that these are not sufficient to confer epigenetic variability. CTCF is enriched at VM-IAPs with binding inversely correlated with DNA methylation. We uncover dynamic physical interactions between cVM-IAPs with low methylation ranges and other genomic loci, suggesting that VM-IAPs have the potential for long-range regulation.

CONCLUSION

Our findings indicate that a recently evolved interplay between genetic sequence, CTCF binding, and DNA methylation at young TEs can result in inter-individual variability in transcriptional outcomes with implications for phenotypic variation.

DNA methylation differs extensively between strains of the same geographical origin and changes with age in Daphnia magna

BACKGROUND

Patterns of methylation influence lifespan, but methylation and lifespan may also depend on diet, or differ between genotypes. Prior to this study, interactions between diet and genotype have not been explored together to determine their influence on methylation. The invertebrate Daphnia magna is an excellent choice for testing the epigenetic response to the environment: parthenogenetic offspring are identical to their siblings (making for powerful genetic comparisons), they are relatively short lived and have well-characterised inter-strain life-history trait differences. We performed a survival analysis in response to caloric restriction and then undertook a 47-replicate experiment testing the DNA methylation response to ageing and caloric restriction of two strains of D. magna.

RESULTS

Methylated cytosines (CpGs) were most prevalent in exons two to five of gene bodies. One strain exhibited a significantly increased lifespan in response to caloric restriction, but there was no effect of food-level CpG methylation status. Inter-strain differences dominated the methylation experiment with over 15,000 differently methylated CpGs. One gene, Me31b, was hypermethylated extensively in one strain and is a key regulator of embryonic expression. Sixty-one CpGs were differentially methylated between young and old individuals, including multiple CpGs within the histone H3 gene, which were hypermethylated in old individuals. Across all age-related CpGs, we identified a set that are highly correlated with chronological age.

CONCLUSIONS

Methylated cytosines are concentrated in early exons of gene sequences indicative of a directed, non-random, process despite the low overall DNA methylation percentage in this species. We identify no effect of caloric restriction on DNA methylation, contrary to our previous results, and established impacts of caloric restriction on phenotype and gene expression. We propose our approach here is more robust in invertebrates given genome-wide CpG distributions. For both strain and ageing, a single gene emerges as differentially methylated that for each factor could have widespread phenotypic effects. Our data showed the potential for an epigenetic clock at a subset of age positions, which is exciting but requires confirmation.

Atypical Genotypes for Canine Agouti Signaling Protein Suggest Novel Chromosomal Rearrangement

Canine coat color is a readily observed phenotype of great interest to dog enthusiasts; it is also an excellent avenue to explore the mechanisms of genetics and inheritance. As such, multiple commercial testing laboratories include basic color alleles in their popular screening panels, allowing for the creation of genotyped datasets at a scale not before appreciated in canine genetic research. These vast datasets have revealed rare genotype anomalies that encourage further exploration of color and pattern inheritance. We previously reported the simultaneous presence of greater than two allele variants at the Agouti Signaling Protein (ASIP) locus in a commercial genotype cohort of 11,790 canids. Here we present additional data to characterize the occurrence of anomalous ASIP genotypes. We document the detection of combinations of three or four ASIP allele variants in 17 dog breeds and Dingoes, at within-breed frequencies of 1.32–63.34%. We analyze the potential impact on phenotype that these allele combinations present, and propose mechanisms that could account for the findings, including: gene recombination, duplication, and incorrect causal variant identification. These findings speak to the accuracy of industry-wide protocols for commercial ASIP genotyping and imply that ASIP should be analyzed via haplotype, rather than using only the existing allele hierarchy, in the future.

DNA methylation at the crossroads of gene and environment interactions

DNA methylation is an epigenetic mark involved in regulating genome function and is critical for normal development in mammals. It has been observed that the developmental environment can lead to permanent changes in gene expression and DNA methylation, at least at 'metastable epialleles'. These are defined as regions of the genome that show a variable epigenetic state that is established early in development and maintained through subsequent cell divisions. However, the majority of the known genome does not behave in this manner. Here, we use the developmental origins of adult disease hypothesis to understand environmental epigenomics. Some challenges to studying how DNA methylation is influenced by the environment include identifying DNA methylation changes associated with an environmental exposure in tissues with a complex cellular composition and at genomic regions for which DNA methylation is dynamically regulated in a cell-type specific manner. We also offer a perspective of how emerging technologies may be useful for dissecting the functional contribution of exposure-associated epigenetic changes and highlight recent evidence that suggests that genomic regions that are absent from genome assemblies may be unappreciated hotspots for environmental modulation of the epigenetic state.

Metastable epialleles and their contribution to epigenetic inheritance in mammals

Many epigenetic differences between individuals are driven by genetic variation. Mammalian metastable epialleles are unusual in that they show variable DNA methylation states between genetically identical individuals. The occurrence of such states across generations has resulted in their consideration by many as strong evidence for epigenetic inheritance in mammals, with the classic A^vy and Axin^Fu mouse models - each products of repeat element insertions - being the most widely accepted examples. Equally, there has been interest in exploring their use as epigenetic biosensors given their susceptibility to environmental compromise. Here we review the classic murine metastable epialleles as well as more recently identified candidates, with the aim of providing a more holistic understanding of their biology. We consider the extent to which epigenetic inheritance occurs at metastable epialleles and explore the limited mechanistic insights into the establishment of their variable epigenetic states. We discuss their environmental modulation and their potential relevance in genome regulation. In light of recent whole-genome screens for novel metastable epialleles, we point out the need to reassess their biological relevance in multi-generational studies and we highlight their value as a model to study repeat element silencing as well as the mechanisms and consequences of mammalian epigenetic stochasticity.

Impact of Historic Migrations and Evolutionary Processes on Human Immunity

The evolution of mankind has constantly been influenced by the pathogens encountered. The various populations of modern humans that ventured out of Africa adapted to different environments and faced a large variety of infectious agents, resulting in local adaptations of the immune system for these populations. The functional variation of immune genes as a result of evolution is relevant in the responses against infection, as well as in the emergence of autoimmune and inflammatory diseases observed in modern populations. Understanding how host-pathogen interactions have influenced the human immune system from an evolutionary perspective might contribute to unveiling the causes behind different immune-mediated disorders and promote the development of new strategies to detect and control such diseases.

Whole Genome Sequencing of the Mutamouse Model Reveals Strain- and Colony-Level Variation, and Genomic Features of the Transgene Integration Site

The MutaMouse transgenic rodent model is widely used for assessing in vivo mutagenicity. Here, we report the characterization of MutaMouse's whole genome sequence and its genetic variants compared to the C57BL/6 reference genome. High coverage (>50X) next-generation sequencing (NGS) of whole genomes from multiple MutaMouse animals from the Health Canada (HC) colony showed ~5 million SNVs per genome, ~20% of which are putatively novel. Sequencing of two animals from a geographically separated colony at Covance indicated that, over the course of 23 years, each colony accumulated 47,847 (HC) and 17,677 (Covance) non-parental homozygous single nucleotide variants. We found no novel nonsense or missense mutations that impair the MutaMouse response to genotoxic agents. Pairing sequencing data with array comparative genomic hybridization (aCGH) improved the accuracy and resolution of copy number variants (CNVs) calls and identified 300 genomic regions with CNVs. We also used long-read sequence technology (PacBio) to show that the transgene integration site involved a large deletion event with multiple inversions and rearrangements near a retrotransposon. The MutaMouse genome gives important genetic context to studies using this model, offers insight on the mechanisms of structural variant formation, and contributes a framework to analyze aCGH results alongside NGS data.

Epigenome modulated xenobiotic detoxification pathways control DMBA-induced breast cancer in agouti Avy/a mice

Environmental xenobiotics with genotoxic activity are carcinogenic. However, individual differences in the susceptibility to xenobiotic-induced breast cancer remain unclear. Since epigenetic modifications could control the expression of metabolic enzymes, our goal was to determine whether epigenome modulated metabolic networks determine susceptibility to xenobiotic-induced breast cancer. The effect of epigenetic background on predisposition to carcinogen 7,12-dimethylbenz(a)anthracene (DMBA)-induced breast cancer development and progression was assessed using the A^vy/a mouse model. In a randomized block design, 22 isogenic A^vy/a (8 yellow, 7 slightly mottled, 7 pseudoagouti) and 8 wild type non-agouti (a/a black) age matched female mice were subjected to DMBA (30 mg/kg per mouse weight) once a week for 6 weeks to induce breast cancer. Compared to pseudoagouti littermates, a significant decrease in tumour latency with increased tumour burden was observed in slightly mottled and yellow littermates (p ≤ 0.05). However, tumour latency and tumour burden were similar in non-agouti a/a mice and A^vy/a cohorts. Network analysis of differentially expressed liver genes identified altered metabolic gene networks among agouti phenotypes. Consequently, in HPLC analyses, DMBA metabolites were significantly increased in A^vy/a pseudoagouti mice (p ≤ 0.05). Relative to A^vy/a slightly mottled, A^vy/a yellow and non-agouti a/a black mice, DMBA metabolites increased nine-, eight-, and four-fold, respectively, in A^vy/a pseudoagouti mice. In agreement with this, seven phase 2 xenobiotic detoxification genes were significantly upregulated in A^vy/a pseudoagouti mice (p ≤ 0.05). The Results from this study suggest that epigenome modulation of xenobiotic detoxification pathways may control xenobiotic-induced breast cancer susceptibility in A^vy/a mice.

Identification, Characterization, and Heritability of Murine Metastable Epialleles: Implications for Non-genetic Inheritance

Generally repressed by epigenetic mechanisms, retrotransposons represent around 40% of the murine genome. At the Agouti viable yellow (Avy) locus, an endogenous retrovirus (ERV) of the intracisternal A particle (IAP) class retrotransposed upstream of the agouti coat-color locus, providing an alternative promoter that is variably DNA methylated in genetically identical individuals. This results in variable expressivity of coat color that is inherited transgenerationally. Here, a systematic genome-wide screen identifies multiple C57BL/6J murine IAPs with Avy epigenetic properties. Each exhibits a stable methylation state within an individual but varies between individuals. Only in rare instances do they act as promoters controlling adjacent gene expression. Their methylation state is locus-specific within an individual, and their flanking regions are enriched for CTCF. Variably methylated IAPs are reprogrammed after fertilization and re-established as variable loci in the next generation, indicating reconstruction of metastable epigenetic states and challenging the generalizability of non-genetic inheritance at these regions.

A survey of inter-individual variation in DNA methylation identifies environmentally responsive co-regulated networks of epigenetic variation in the human genome

While population studies have resulted in detailed maps of genetic variation in humans, to date there are few robust maps of epigenetic variation. We identified sites containing clusters of CpGs with high inter-individual epigenetic variation, termed Variably Methylated Regions (VMRs) in five purified cell types. We observed that VMRs occur preferentially at enhancers and 3' UTRs. While the majority of VMRs have high heritability, a subset of VMRs within the genome show highly correlated variation in trans, forming co-regulated networks that have low heritability, differ between cell types and are enriched for specific transcription factor binding sites and biological pathways of functional relevance to each tissue. For example, in T cells we defined a network of 95 co-regulated VMRs enriched for genes with roles in T-cell activation; in fibroblasts a network of 34 co-regulated VMRs comprising all four HOX gene clusters enriched for control of tissue growth; and in neurons a network of 18 VMRs enriched for roles in synaptic signaling. By culturing genetically-identical fibroblasts under varying environmental conditions, we experimentally demonstrated that some VMR networks are responsive to the environment, with methylation levels at these loci changing in a coordinated fashion in trans dependent on cellular growth. Intriguingly these environmentally-responsive VMRs showed a strong enrichment for imprinted loci (p<10-80), suggesting that these are particularly sensitive to environmental conditions. Our study provides a detailed map of common epigenetic variation in the human genome, showing that both genetic and environmental causes underlie this variation.

Extensive sex differences at the initiation of genetic recombination

Meiotic recombination differs between males and females; however, when and how these differences are established is unknown. Here we identify extensive sex differences at recombination initiation by mapping hotspots of meiotic DNA double-strand breaks in male and female mice. Contrary to past findings in humans, few hotspots are used uniquely in either sex. Instead, grossly different recombination landscapes result from up to 15-fold differences in hotspot usage between males and females. Indeed, most recombination occurs at sex-biased hotspots. Sex-biased hotspots appear to be partly determined by chromosome structure, and DNA methylation, absent in females at the onset of meiosis, plays a substantial role. Sex differences are also evident later in meiosis as the repair frequency of distal meiotic breaks as crossovers diverges in males and females. Suppression of distal crossovers may help to minimize age-related aneuploidy that arises due to cohesion loss during dictyate arrest in females.

Establishment of environmentally sensitive DNA methylation states in the very early human embryo

The molecular mechanisms responsible for the developmental origins of later disease are currently unknown. We previously demonstrated that women's periconceptional nutrition predicts their offspring's DNA methylation at metastable epialleles (MEs). We present a genome-wide screen yielding 687 MEs and track their trajectories across nine developmental stages in human in vitro fertilization embryos. MEs exhibit highly unusual methylation dynamics across the implantation-gastrulation transition, producing a large excess of intermediate methylation states, suggesting the potential for differential programming in response to external signals. Using a natural experiment in rural Gambia, we show that genomic regions sensitive to season of conception are highly enriched for MEs and show similar atypical methylation patterns. MEs are enriched for proximal enhancers and transcription start sites and are influenced by genotype. Together, these observations position MEs as distinctive epigenomic features programmed in the early embryo, sensitive to genetic and periconceptional environment, and with the potential to influence phenotype.

Cell-type-specific brain methylomes profiled via ultralow-input microfluidics

Methylomic analyses typically require substantial amounts of DNA, thus hindering studies involving scarce samples. Here, we show that microfluidic diffusion-based reduced representative bisulfite sequencing (MID-RRBS) permits high-quality methylomic profiling with nanogram-to-single-cell quantities of starting DNA. We used the microfluidic device, which allows for efficient bisulfite conversion with high DNA recovery, to analyse genome-wide DNA methylation in cell nuclei isolated from mouse brains and sorted into NeuN+ (primarily neuronal) and NeuN- (primarily glial) fractions, and to establish cell-type-specific methylomes. Genome-wide methylation and methylation in low-CpG-density promoter regions showed distinct patterns for NeuN+ and NeuN- fractions from the mouse cerebellum. The identification of substantial variations in the methylomic landscapes of the NeuN+ fraction of the frontal cortex of mice chronically treated with an atypical antipsychotic drug suggests that this technology can be broadly used for cell-type-specific drug profiling and for the study of drug-methylome interactions.

Two predominant MUPs, OBP3 and MUP13, are male pheromones in rats

BACKGROUND

In rats, urine-borne male pheromones comprise organic volatile compounds and major urinary proteins (MUPs). A number of volatile pheromones have been reported, but no MUP pheromones have been identified in rat urine.

RESULTS

We used sodium dodecyl sulphate-polyacrylamide gel electrophoresis (SDS-PAGE), isoelectric focusing electrophoresis (IEF), nano-liquid chromatography-tandem mass spectrometry (nLC-MS/MS) after in gel digestion of the proteins and quantitative real-time PCR (qRT-PCR) and showed that the levels of two MUPs, odorant-binding protein 3 (OBP3) (i.e. PGCL4) and MUP13 (i.e. PGCL1), in urine and their mRNAs in liver were higher in males than in females and were suppressed by orchidectomy and restored by testosterone treatment (T treatment). We then generated recombinant MUPs (rMUPs) and found that the sexual attractiveness of urine from castrated males to females significantly increased after the addition of either recombinant OBP3 (rOBP3) or recombinant MUP13 (rMUP13). Using c-Fos immunohistochemistry, we further examined neuronal activation in the brains of female rats after they sniffed rOBP3 or rMUP13. Both rOBP3 and rMUP13 activated the accessory olfactory bulb (AOB), medial preoptic area (MPA), bed nucleus of the stria terminalis (BST), medial amygdala (MeA), posteromedial cortical amygdala (PMCo) and ventromedial nucleus of the hypothalamus (VMH), which participate in the neural circuits responsible for pheromone-induced sexual behaviours. In particular, more c-Fos-immunopositive (c-Fos-ir) cells were observed in the posterior AOB than in the anterior AOB.

CONCLUSIONS

The expression of OBP3 and MUP13 was male-biased and androgen-dependent. They attracted females and activated brain areas related to sexual behaviours in female rats, suggesting that both OBP3 and MUP13 are male pheromones in rats. Particularly, an OBP excreted into urine was exemplified to be a chemical signal.

Adiposity QTL Adip20 decomposes into at least four loci when dissected using congenic strains

An average mouse in midlife weighs between 25 and 30 g, with about a gram of tissue in the largest adipose depot (gonadal), and the weight of this depot differs between inbred strains. Specifically, C57BL/6ByJ mice have heavier gonadal depots on average than do 129P3/J mice. To understand the genetic contributions to this trait, we mapped several quantitative trait loci (QTLs) for gonadal depot weight in an F2 intercross population. Our goal here was to fine-map one of these QTLs, Adip20 (formerly Adip5), on mouse chromosome 9. To that end, we analyzed the weight of the gonadal adipose depot from newly created congenic strains. Results from the sequential comparison method indicated at least four rather than one QTL; two of the QTLs were less than 0.5 Mb apart, with opposing directions of allelic effect. Different types of evidence (missense and regulatory genetic variation, human adiposity/body mass index orthologues, and differential gene expression) implicated numerous candidate genes from the four QTL regions. These results highlight the value of mouse congenic strains and the value of this sequential method to dissect challenging genetic architecture.

Genetic and phenotypic characterization of the novel mouse substrain C57BL/6N Korl with increased body weight

In inbred mouse lines, there is generally little genetic difference between individuals. This small genetic variability facilitates carrying out research on minute changes of various traits and the gene pool. Also, characterizing the diversity and detecting selective genetic and phenotypic signatures are crucial to understanding the genomic basis of a population and to identify specific patterns of evolutionary change. In this study, we investigated the underlying genetic profiles of a newly developed mouse strain, C57BL/6NKorl (Korl), established through sibling mating over 30 generations. To analyse the distinctive genomic features of Korl mice, we used whole-genome sequencing from six samples, which were compared to those of other C57BL/6N-based mouse strains. Korl strain-specific polymorphisms were identified and signatures of a selective sweep were detected. In particular, the candidate genes related to the increased body weight of the Korl strain were identified. Establishment of the genetic profile of Korl mice can provide insight into the inbreeding-induced changes to the gene pool, and help to establish this strain as a useful model for practical and targeted research purposes.

Inhibition of microRNA-214 promotes epithelial-mesenchymal transition process and induces interstitial cystitis in postmenopausal women by upregulating Mfn2

Our study aims to investigate the roles that microRNA-214 (miR-214) plays in the epithelial mesenchymal transition (EMT) process and the development of interstitial cystitis (IC) in postmenopausal women by targeting Mitofusin 2 (Mfn2). IC bladder tissues and adjacent normal bladder tissues were collected from postmenopausal women. Immunohistochemistry (IHC) staining was conducted. The target relationship between miR-214 and Mfn2 was determined by a dual luciferase reporter gene assay. Adipose-derived mesenchymal stem cells (ADMSCs) were extracted from postmenopausal rats and assigned to the blank, mimics, miR-214 inhibitors, mimics negative control (NC), inhibitors NC, Mfn2 siRNA, miR-214 inhibitors and Mfn2 siRNA groups. Exosomes secreted by transfected ADMSCs were instilled into the bladders of postmenopausal rats. The expression of miR-214 and Mfn2 mRNA and EMT markers was assessed by qRT-PCR and western blotting. It was confirmed that Mfn2 was the target gene of miR-214 in IC. Compared with the normal bladder tissues, miR-214 decreased, but Mfn2 increased in IC bladder tissues. Compared with the blank group, the expression of miR-214 and the expression levels of N-cadherin, Fibronectin, Twist1, Snail and Vimentin mRNA and protein increased, whereas the expression levels of Mfn2, E-cadherin and ZO-1 mRNA and protein decreased in the miR-214 mimics and Mfn2 groups. The expression of MiR-214 and the expression levels of N-cadherin, Fibronectin, Twist1, Snail and Vimentin mRNA and protein decreased, whereas the expression levels of Mfn2, E-cadherin and ZO-1 mRNA and protein increased in the miR-214 inhibitors group. Our findings indicate that the inhibition of miR-214 promotes the EMT process and contributes to bladder wall fibrosis by up-regulating Mfn2, thus leading to the occurrence of IC in postmenopausal women.

Dynamic silencing of somatic L1 retrotransposon insertions reflects the developmental and cellular contexts of their genomic integration

Background

The ongoing mobilization of mammalian transposable elements (TEs) contributes to natural genetic variation. To survey the epigenetic control and expression of reporter genes inserted by L1 retrotransposition in diverse cellular and genomic contexts, we engineered highly sensitive, real-time L1 retrotransposon reporter constructs.

Results

Here we describe different patterns of expression and epigenetic controls of newly inserted sequences retrotransposed by L1 in various somatic cells and tissues including cultured human cancer cells, mouse embryonic stem cells, and tissues of pseudofounder transgenic mice and their progeny. In cancer cell lines, the newly inserted sequences typically underwent rapid transcriptional gene silencing, but they lacked cytosine methylation even after many cell divisions. L1 reporter expression was reversible and oscillated frequently. Silenced or variegated reporter expression was strongly and uniformly reactivated by treatment with inhibitors of histone deacetylation, revealing the mechanism for their silencing. By contrast, de novo integrants retrotransposed by L1 in pluripotent mouse embryonic stem (ES) cells underwent rapid silencing by dense cytosine methylation. Similarly, de novo cytosine methylation also was identified at new integrants when studied in several distinct somatic tissues of adult founder mice. Pre-existing L1 elements in cultured human cancer cells were stably silenced by dense cytosine methylation, whereas their transcription modestly increased when cytosine methylation was experimentally reduced in cells lacking DNA methyltransferases DNMT1 and DNMT3b. As a control, reporter genes mobilized by piggyBac (PB), a DNA transposon, revealed relatively stable and robust expression without apparent silencing in both cultured cancer cells and ES cells.

Conclusions

We hypothesize that the de novo methylation marks at newly inserted sequences retrotransposed by L1 in early pre-implantation development are maintained or re-established in adult somatic tissues. By contrast, histone deacetylation reversibly silences L1 reporter insertions that had mobilized at later timepoints in somatic development and differentiation, e.g., in cancer cell lines. We conclude that the cellular contexts of L1 retrotransposition can determine expression or silencing of newly integrated sequences. We propose a model whereby reporter expression from somatic TE insertions reflects the timing, molecular mechanism, epigenetic controls and the genomic, cellular and developmental contexts of their integration.

Electronic supplementary material

The online version of this article (doi:10.1186/s13100-017-0091-2) contains supplementary material, which is available to authorized users.

Effects of G6pc2 deletion on body weight and cholesterol in mice

Genome-wide association study (GWAS) data have linked the G6PC2 gene to variations in fasting blood glucose (FBG). G6PC2 encodes an islet-specific glucose-6-phosphatase catalytic subunit that forms a substrate cycle with the beta cell glucose sensor glucokinase. This cycle modulates the glucose sensitivity of insulin secretion and hence FBG. GWAS data have not linked G6PC2 to variations in body weight but we previously reported that female C57BL/6J G6pc2-knockout (KO) mice were lighter than wild-type littermates on both a chow and high-fat diet. The purpose of this study was to compare the effects of G6pc2 deletion on FBG and body weight in both chow-fed and high-fat-fed mice on two other genetic backgrounds. FBG was reduced in G6pc2 KO mice largely independent of gender, genetic background or diet. In contrast, the effect of G6pc2 deletion on body weight was markedly influenced by these variables. Deletion of G6pc2 conferred a marked protection against diet-induced obesity in male mixed genetic background mice, whereas in 129SvEv mice deletion of G6pc2 had no effect on body weight. G6pc2 deletion also reduced plasma cholesterol levels in a manner dependent on gender, genetic background and diet. An association between G6PC2 and plasma cholesterol was also observed in humans through electronic health record-derived phenotype analyses. These observations suggest that the action of G6PC2 on FBG is largely independent of the influences of environment, modifier genes or epigenetic events, whereas the action of G6PC2 on body weight and cholesterol are influenced by unknown variables.

Inter-individual variation in DNA methylation is largely restricted to tissue-specific differentially methylated regions in maize

BACKGROUND

Variation in DNA methylation across distinct genetic populations, or in response to specific biotic or abiotic stimuli, has typically been studied in leaf DNA from pooled individuals using either reduced representation bisulfite sequencing, whole genome bisulfite sequencing (WGBS) or methylation sensitive amplified polymorphism (MSAP). The latter represents a useful alterative when sample size is large, or when analysing methylation changes in genomes that have yet to be sequenced. In this study we compared variation in methylation across ten individual leaf and endosperm samples from maize hybrid and inbred lines using MSAP. We also addressed the methodological implications of analysing methylation variation using pooled versus individual DNA samples, in addition to the validity of MSAP compared to WGBS. Finally, we analysed a subset of variable and non-variable fragments with respect to genomic location, vicinity to repetitive elements and expression patterns across leaf and endosperm tissues.

RESULTS

On average, 30% of individuals showed inter-individual methylation variation, mostly of leaf and endosperm-specific differentially methylated DNA regions. With the exception of low frequency demethylation events, the bulk of inter-individual methylation variation (84 and 80% in leaf and endosperm, respectively) was effectively captured in DNA from pooled individuals. Furthermore, available genome-wide methylation data largely confirmed MSAP leaf methylation profiles. Most variable methylation that mapped within genes was associated with CG methylation, and many of such genes showed tissue-specific expression profiles. Finally, we found that the hAT DNA transposon was the most common class II transposable element found in close proximity to variable DNA regions.

CONCLUSIONS

The relevance of our results with respect to future studies of methylation variation is the following: firstly, the finding that inter-individual methylation variation is largely restricted to tissue-specific differentially methylated DNA regions, underlines the importance of tissue-type when analysing the methylation response to a defined stimulus. Secondly, we show that pooled sample-based MSAP studies are methodologically appropriate to study methylation variation. Thirdly, we confirm that MSAP is a powerful tool when WGBS is not required or feasible, for example in plant species that have yet to be sequenced.

Erasure and reestablishment of random allelic expression imbalance after epigenetic reprogramming

Clonal level random allelic expression imbalance and random monoallelic expression provides cellular heterogeneity within tissues by modulating allelic dosage. Although such expression patterns have been observed in multiple cell types, little is known about when in development these stochastic allelic choices are made. We examine allelic expression patterns in human neural progenitor cells before and after epigenetic reprogramming to induced pluripotency, observing that loci previously characterized by random allelic expression imbalance (0.63% of expressed genes) are generally reset to a biallelic state in induced pluripotent stem cells (iPSCs). We subsequently neuralized the iPSCs and profiled isolated clonal neural stem cells, observing that significant random allelic expression imbalance is reestablished at 0.65% of expressed genes, including novel loci not found to show allelic expression imbalance in the original parental neural progenitor cells. Allelic expression imbalance was associated with altered DNA methylation across promoter regulatory regions, with clones characterized by skewed allelic expression being hypermethylated compared to their biallelic sister clones. Our results suggest that random allelic expression imbalance is established during lineage commitment and is associated with increased DNA methylation at the gene promoter.

Preclinical identification of vaccine induced protective correlates in human leukocyte antigen expressing transgenic mice infected with Coccidioides posadasii

There is an emerging interest to develop human vaccines against medically-important fungal pathogens and a need for a preclinical animal model to assess vaccine efficacies and protective correlates. HLA-DR4 (DRB1∗0401 allele) transgenic mice express a human major histocompatibility complex class II (MHC II) receptor in such a way that CD4⁺ T-cell response is solely restricted by this human molecule. In this study HLA-DR4 transgenic mice were immunized with a live-attenuated vaccine (ΔT) and challenged by the intranasal route with 50-70 Coccidioides posadasii spores, a potentially lethal dose. The same vaccination regimen offers 100% survival for C57BL/6 mice. Conversely, ΔT-vaccinated HLA-DR4 mice displayed 3 distinct manifestations of Coccidioides infection including 40% fatal acute (FAD), 30% disseminated (DD) and 30% pulmonary disease (PD). The latter 2 groups of mice had reduced loss of body weight and survived to at least 50days postchallenge (dpc). These results suggest that ΔT vaccinated HLA-DR4 mice activated heterogeneous immunity against pulmonary Coccidioides infection. Vaccinated HLA-DR4 mice displayed early expansion of Th1 and Th17 cells and recruitment of inflammatory innate cells into Coccidioides-infected lungs during the first 9dpc. While contraction rates of Th cells and the inflammatory response during 14-35dpc significantly differed among the 3 groups of vaccinated HLA-DR4 mice. The FAD group displayed a sharply reduced Th1 and Th17 response, while overwhelmingly recruiting neutrophils into lungs during 9-14days. The FAD group approached moribund by 14dpc. In contrast, vaccinated HLA-DR4 survivors gradually contracted Th cells and inflammatory response with the greatest rate in the PD group. While vaccinated HLA-DR4 mice are susceptible to Coccidioides infection, they are useful for evaluation of vaccine efficacy and identification of immunological correlates against this mycosis.

Epigenetic studies in Developmental Origins of Health and Disease: pitfalls and key considerations for study design and interpretation

The field of Developmental Origins of Health and Disease (DOHaD) seeks to understand the relationships between early-life environmental exposures and long-term health and disease. Until recently, the molecular mechanisms underlying these phenomena were poorly understood; however, epigenetics has been proposed to bridge the gap between the environment and phenotype. Epigenetics involves the study of heritable changes in gene expression, which occur without changes to the underlying DNA sequence. Different types of epigenetic modifications include DNA methylation, post-translational histone modifications and non-coding RNAs. Increasingly, changes to the epigenome have been associated with early-life exposures in both humans and animal models, offering both an explanation for how the environment may programme long-term health, as well as molecular changes that could be developed as biomarkers of exposure and/or future disease. As such, epigenetic studies in DOHaD hold much promise; however, there are a number of factors which should be considered when designing and interpreting such studies. These include the impact of the genome on the epigenome, the tissue-specificity of epigenetic marks, the stability (or lack thereof) of epigenetic changes over time and the importance of associating epigenetic changes with changes in transcription or translation to demonstrate functional consequences. In this review, we discuss each of these key concepts and provide practical strategies to mitigate some common pitfalls with the aim of providing a useful guide for future epigenetic studies in DOHaD.

Chromatin variation associated with liver metabolism is mediated by transposable elements

Background

Functional regulatory regions in eukaryotic genomes are characterized by the disruption of nucleosomes leading to accessible chromatin. The modulation of chromatin accessibility is one of the key mediators of transcriptional regulation, and variation in chromatin accessibility across individuals has been linked to complex traits and disease susceptibility. While mechanisms responsible for chromatin variation across individuals have been investigated, the overwhelming majority of chromatin variation remains unexplained. Furthermore, the processes through which the variation of chromatin accessibility contributes to phenotypic diversity remain poorly understood.

Results

We profiled chromatin accessibility in liver from seven strains of mice with phenotypic diversity in response to a high-fat/high-sucrose (HF/HS) diet and identified reproducible chromatin variation across the individuals. We found that sites of variable chromatin accessibility were more likely to coincide with particular classes of transposable elements (TEs) than sites with common chromatin signatures. Evolutionarily younger long interspersed nuclear elements (LINEs) are particularly likely to harbor variable chromatin sites. These younger LINEs are enriched for binding sites of immune-associated transcription factors, whereas older LINEs are enriched for liver-specific transcription factors. Genomic region enrichment analysis indicates that variable chromatin sites at TEs may function to regulate liver metabolic pathways. CRISPR-Cas9 deletion of a number of variable chromatin sites at TEs altered expression of nearby metabolic genes. Finally, we show that polymorphism of TEs and differential DNA methylation at TEs can both influence chromatin variation.

Conclusions

Our results demonstrate that specific classes of TEs show variable chromatin accessibility across strains of mice that display phenotypic diversity in response to a HF/HS diet. These results indicate that chromatin variation at TEs is an important contributor to phenotypic variation among populations.

Electronic supplementary material

The online version of this article (doi:10.1186/s13072-016-0078-0) contains supplementary material, which is available to authorized users.

Individuality and Transgenerational Inheritance of Social Dominance and Sex Pheromones in Isogenic Male Mice

Phenotypic variation and its epigenetic regulations within the inbred isogenic mice have long intrigued biologists. Here, we used inbred C57BL/6 mice to examine the individual differences and the inheritance of social dominance and male pheromones, expecting to create a model for studying the underlying epigenetic mechanisms for the evolution of these traits. We used a repeated male-male contest paradigm to form stable dominance-submission relationships between paired males and make superior or inferior quality manifest. Females showed olfactory preferences for the urine of dominant males to that of subordinate opponents. Gas chromatography-mass spectrometer analysis revealed that dominance-related or superior quality related pheromones were actually exaggerated male pheromone components (e.g., E-β-farnesene, hexadecanol, and 1-hexadecanol acetate) of preputial gland origin. Although the socially naïve sons of both dominant and subordinate males elicited the same female attraction when reaching adulthood, the former could dominated over the latter during undergoing the male-male competition and then gained more attraction of females. Our results demonstrated that social dominance or superior quality and the related pheromones were heritable and could be expressed through the interaction between aggression-related epigenotypes and male-male contests. It suggested that the evolution of sexually selected traits could be epigenetically determined and promoted through female mate choice. The epigenetic mechanisms driving the individual differences in behavior and male pheromones deserve further studies.

Epigenetic Inheritance and Its Role in Evolutionary Biology: Re-Evaluation and New Perspectives

Epigenetics increasingly occupies a pivotal position in our understanding of inheritance, natural selection and, perhaps, even evolution. A survey of the PubMed database, however, reveals that the great majority (>93%) of epigenetic papers have an intra-, rather than an inter-generational focus, primarily on mechanisms and disease. Approximately ~1% of epigenetic papers even mention the nexus of epigenetics, natural selection and evolution. Yet, when environments are dynamic (e.g., climate change effects), there may be an “epigenetic advantage” to phenotypic switching by epigenetic inheritance, rather than by gene mutation. An epigenetically-inherited trait can arise simultaneously in many individuals, as opposed to a single individual with a gene mutation. Moreover, a transient epigenetically-modified phenotype can be quickly “sunsetted”, with individuals reverting to the original phenotype. Thus, epigenetic phenotype switching is dynamic and temporary and can help bridge periods of environmental stress. Epigenetic inheritance likely contributes to evolution both directly and indirectly. While there is as yet incomplete evidence of direct permanent incorporation of a complex epigenetic phenotype into the genome, doubtlessly, the presence of epigenetic markers and the phenotypes they create (which may sort quite separately from the genotype within a population) will influence natural selection and, so, drive the collective genotype of a population.