Abstract
Genetic variants can alter the profile of heritable molecules such as small RNAs in sperm and oocytes, and in this manner ancestral genetic variants can have a significant effect on offspring phenotypes even if they are not themselves inherited. Here we show that wild type female mice descended from ancestors with a mutation in the mammalian germ cell gene Khdc3 have hepatic metabolic defects that persist over multiple generations. We find that genetically wild type females descended from Khdc3 mutants have transcriptional dysregulation of critical hepatic metabolic genes, which persist over multiple generations and pass through both female and male lineages. This was associated with dysregulation of hepatically-metabolized molecules in the blood of these wild type mice with mutational ancestry. The oocytes of Khdc3-null females, as well as their wild type descendants, had dysregulation of multiple small RNAs, suggesting that these epigenetic changes in the gametes transmit the phenotype between generations. Our results demonstrate that ancestral mutation in Khdc3 can produce transgenerational inherited phenotypes, potentially indefinitely.
Introduction
It is becoming increasingly apparent that non-DNA molecules inherited from germ cells contribute significantly to the transmission of traits and diseases across generations. This has been demonstrated in response to a variety of exposures such as diet or stress in model organisms (1, 2) and has been confirmed in human epidemiological studies (3, 4). The mechanisms driving this process remain poorly defined and are complicated by the erasure of acquired epigenetic molecules such as DNA methylation early in embryonic development. Expression of germ cell small RNAs such as microRNA (miRNA) and tRNA fragments (tsRNA) are responsive to exposures and inherited at fertilization, making them potential candidates for passage of phenotypic information across generations (5). Microinjection of these small RNAs from the germ cells of an exposed mouse into an unexposed embryo is sufficient to recapitulate the descendant’s phenotypes, demonstrating causality (1). However, the absence of an RNA-dependent RNA polymerase in mammals implies that small RNAs by themselves cannot propagate phenotypes across multiple generations.
Non-inherited ancestral genetic variants can also affect descendants’ risk of disease, as has been demonstrated in cancer (6, 7), thyroid hormone metabolism (8), body weight (9), anxiety (10), type I diabetes (11), and folate metabolism (12). This phenomenon is often referred to as “genetic nurture” based on the effect of non-inherited DNA variants nurturing phenotypes in the next generation (13). In all described cases, the observed phenotypes disappear after a few generations. The molecular mechanisms driving these cross-generational traits, in the absence of inheritance of the causal variant, remains largely undescribed, with the exception of paramutation effects that are driven by inherited RNA molecules (14, 15). Some of the genes associated with this phenomenon are thought to directly alter the molecular profile of the gamete, for example with testicular germ cell tumors and RNA-modifying genes that are expressed in sperm (6). Other reported examples involve genes that do not have any known function in the germ cells, and likely alter descendants’ phenotypes indirectly, by changing underlying physiology in a manner that subsequently alters germ cell heritable molecules, not dissimilar from the mechanisms of exposure-based changes to heritable germ cell molecules.
Khdc3 is a mammalian gene expressed in the male and female germ cells that encodes a protein containing an RNA-binding KH domain that localizes to a multi-protein complex at the oocyte periphery known as the subcortical maternal complex (16, 17). Khdc3-null female mice have decreased fertility caused by defects in maintaining euploidy during embryogenesis (18). Human females with homozygous mutations in the ortholog KHDC3L are infertile, which has been associated with abnormal DNA methylation in oocytes (16). KHDC3L (and Khdc3) does not localize to the nucleus and does not contain a DNA methyltransferase domain, suggesting that the DNA methylation defects are secondary outcome and not the main function of this gene.
We demonstrate that wild type female descendants of Khdc3-null mice have dysfunctional expression of hepatic metabolic genes that persists over multiple generations, and that this effect is transmitted from both male and female mutant ancestors. This corresponds with abnormal levels of hepatic metabolites in the serum of these wild type mice with mutant ancestors. The persistence of these abnormalities in genetically wild type mice suggests that altered epigenetic information in the germ cells is transmitting the inherited metabolic phenotypes. Accordingly, we observed that the oocytes of Khdc3-null females and their wild type descendants have multiple dysregulated miRNAs and tsRNAs, suggesting a mechanism of inheritance.
Results
Wild type females descended from Khdc3-null ancestors display hepatic transcriptional dysregulation
Global transcriptome analysis of Khdc3-null oocytes revealed significant dysregulation of genes important in metabolic processes regulated by the liver, especially the metabolism of lipids (atherosclerosis, PPAR signaling, and pantothenate and CoA metabolism) and of carbohydrates (glycolysis/gluconeogenesis, and fructose and mannose metabolism) (Fig. 1A) (19–22). Despite the abnormal expression of metabolic genes in the Khdc3-null oocyte, this gene is only expressed in female reproductive tissue, with no detectable expression in the liver and minimal expression in other metabolic tissues such as adipose or pancreas (Fig. 1B).
Expression of Cyp17a1, a gene central to lipid metabolism (23, 24), was increased in the livers of Khdc3-null (knockout, or “KO”) females, despite a lack of expression or any known function of Khdc3 in wild type livers. Elevated expression of Cyp17a1 was observed in KO females generated from KO parents (referred to as “KOKO”) as well as KO females generated by mating Khdc3 heterozygote parents (referred to as “KO*”) (Fig. 1C). Cyp17a1 expression was also significantly increased, to the same extent as KOKO and KO* females, in genetically wild type females descended from heterozygous parents (referred to as “WT*” to indicate genetically wild type mice that have descended from Khdc3 heterozygous mutant parents and a combination of wild type and Khdc3-null grandparents). Thus, having an ancestor that carried a Khdc3 mutation was associated with elevated Cyp17a1 expression, and the Khdc3 genotype had no effect on Cyp17a1 expression, revealing a DNA-independent form of inheritance (Fig. 1C). Increased Cyp17a1 mRNA expression occurred in WT* females generated from all possible combinations of male and female KO grandparents (Fig. 1D). In addition to conventional genotyping, the wild type genotype of these mice was confirmed with RNA-Seq of ovaries (Supplemental Figure 1).
RNA-Seq of female WT* livers revealed a pattern of global transcriptional dysregulation that overlapped significantly with KO* and KOKO females (Fig. 1E-F), demonstrating that shared ancestry rather than genotype accounts for the transcriptional abnormalities in WT* mice. This is supported by the observation that most significantly dysregulated genes in WT* livers were similarly dysregulated in KO* livers (Fig. 1G). Gene ontology of the common dysregulated genes in livers of WT* and KO* mice revealed enrichment in metabolic processes in which the liver plays a central role, especially lipid and glucose metabolism (22–25) (Fig. 1H).
WT* defects persist over multiple generations
When WT* male and female mice were mated with each other to generate the next generation (referred to as WT**, signifying the 2nd generation of genetically wild type mice), hepatic transcriptional dysregulation of metabolic genes persisted. There was a similar pattern of dysregulation in WT** females derived from either male or female Khdc3-null ancestors (denoted WT**(P) and WT**(M), respectively, denoting Paternal or Maternal ancestral mutant history) (Fig. 2A). Comparison of dysregulated genes between WT**(P) and WT**(M) females showed significant overlap (Fig. 2C) revealing that the same pattern of abnormal gene expression is observed in wild type females descended from both male and female mutant ancestors.
WT(P) male and female mice were mated with each other in order to examine passage these defects over successive generations. Descendants of male mutants were examined to avoid potential confounding from non-germ cell mechanisms of inheritance such as mitochondrial inheritance, fetal-maternal communication across the placenta, or transmission of molecules via breastmilk. Expression of Cyp17a1 and 2610507I01Rik, two genes dysregulated in WT*, KO*, KOKO, and WT**(P) female mice, demonstrated persistent and worsening dysregulated expression in the 2nd through 6th generation of females (Fig. 2D). In a minority of mice, expression reverted back to WT levels (Fig. 2D). Importantly, those mice that reverted to WT levels of Cyp17a1 maintained elevated expression of 2610507I01Rik, and vice-versa, suggesting that the abnormal expression of these genes is driven by independent units of inheritance.
Khdc3-null females were outcrossed with a true wild type male without any ancestral history of Khdc3 mutation, in order to examine the persistence of the observed defects through only the maternal line without any contribution from males that descended from mutants. 100% of the F2 outcrossed females, descended from a maternal grandmother mutant, had persistent transcriptional dysregulation of metabolism-related genes in the liver (Fig. 2E). Many of these dysregulated genes overlapped with WT** females descended from both male (WT**(P)) and female (WT**(M)) homozygous-null ancestors (Fig 2F). In this mating scheme, persistence of hepatic transcriptional dysregulation in all 7 examined females is not consistent with mechanisms of inheritance involving in-cis transmission, including an unaccounted DNA polymorphism/mutation, DNA methylation, or histone protein modification, in which any chromosome from the Khdc3-null grandmother would have been inherited by ∼50% of the F2 outcross generation.
WT* females have abnormal levels of serum hepatic metabolites
Global serum metabolomic analysis was performed to examine the metabolic consequences of the observed transcriptional dysregulation. Both WT**(P) and WT****(P) female mice had abnormal levels of multiple bile acids, which are synthesized in the liver (26–29)(Fig. 3A). The magnitude of dysregulation was similar to that observed in KOKO females (Fig. 3A), suggesting that mutational ancestry was more important than the organism’s genotype. These females also had abnormal concentrations of various xenobiotics (Fig. 3B) and metabolic cofactors (Fig. 3C), both of which are associated with hepatic metabolism (30–33).
WT****(P) mice were metabolically challenged with 8 weeks of a high fat diet (HFD). When compared with true wild type mice exposed to a HFD, WT****(P) females exposed to a HFD showed significant increase in multiple lipid molecules that were not dysregulated in WT****(P) females that consumed a conventional diet, or in WT mice consuming a HFD (Fig. 3D). Thus, there were latent metabolic abnormalities that could not be detected unless stressed with a high fat diet. The HFD WT****(P) serum also had decreased levels of multiple other hepatic metabolites (34–36)(Fig. 3E), some of which were also decreased in WT HFD females although at a much lesser magnitude (Fig. 3F). These findings demonstrate that the hepatic transcriptional dysregulation observed in these mice has significant effects on multiple metabolites found in the blood.
Small RNA dysregulation
Livers of WT**(P) females had dysfunctional expression of multiple small RNA-processing genes, especially genes that regulate tRNA (Trmt9b (37), Ang (38), Nsun6 (39), Elac1 (40)) (Fig. 3A) and miRNA processing (Mettl1 (41, 42), Ago2 (43), Exosc10, Syncrip (44)) (Fig. 4A). Based on this observation, in combination with the fact that gamete small RNAs can transmit non-genetically inherited phenotypes in mice, we hypothesized that the observed phenotypes were driven by the defective inheritance of small RNAs from the germ cells of Khdc3 mutant ancestors. Small RNA-Seq of KOKO oocytes revealed dysregulation of multiple miRNAs and tRNA fragments, with minimal piRNA dysregulation (Fig. 4B). In WT**(P) oocytes, there was also abnormal expression of tRNA fragments and miRNA, most of which were different small RNAs than the ones dysregulated in KOKO oocytes (Fig. 4C). Of note, the WT**(P) oocytes had downregulation of the tRNA fragment Gly-GCC, for which expression in sperm has been associated with inherited metabolic and hepatic gluconeogenesis phenotypes in offspring (1, 45, 46). There were 3 tRNA fragments and 18 miRNAs that were commonly dysregulated in both KOKO and WT**(P) oocytes (Fig. 3D), suggesting that their abnormal expression is established in KOKO oocytes and is not normalized with reintroduction of a wild type Khdc3 allele. The function of most of these commonly dysregulated tsRNAs and miRNAs remains undescribed, however miR-107, which was upregulated in both KOKO and WT**(P) oocytes, is involved in hepatic lipid metabolism (47).
Discussion
The presence of abnormal phenotypes in wild type organisms descended from mutant-carrying ancestors has been described in mice and humans, however all previously observed defects disappear after one or two generations. This can occur from mutations that affect epigenetic phenomenon in the germ cell (6, 7, 48), but in other cases the mutation affects a process that has no clear connection with heritable molecules in the germ cell, such as with thyroid metabolism or type I diabetes (8, 11).
We find that loss-of-function mutation in Khdc3 alters the hepatic metabolism of female wild type descendants over at least 6 generations in a manner that cannot be rescued with re-introduction of the Khdc3 allele. This is paralleled by abnormal levels of multiple metabolites in the serum of these mice. Because there is no detectable Khdc3 expression in the liver, we suspect that the defects in WT* mice are caused by inherited molecules that affect hepatic metabolism independent of the inheritance of a functional Khdc3 gene. We expect that other metabolic tissues such as the pancreas and adipose will also demonstrate evidence of metabolic dysregulation.
We have demonstrated that Khdc3-null oocytes have abnormal expression of multiple miRNAs and tsRNAs, some of which persist over multiple generations in WT* mice, providing a potential mechanism of inheritance. Microinjection experiments using sperm RNA have demonstrated that small RNAs are sufficient to drive metabolic derangements in offspring (1, 49), however the maintenance of defects to the successive generation remains unexplained. The amount of RNA needed for these experiments currently prevents using oocyte RNA in such an experimental paradigm, because of the low total RNA yield obtained from the oocytes of a single mouse.
We propose that the metabolic phenotype is at least partially driven by abnormal oocyte small RNAs, however the persistence of a phenotype over multiple generations requires maintenance of small RNAs via a RNA-dependent RNA polymerase, which mammals do not possess. It is possible that small RNAs inherited from the oocyte drive abnormal metabolic profiles in offspring, and oocyte small RNA defects are re-established in the germ cells of the following generation in response to the abnormal metabolism. A recent study demonstrated that transfer of serum from a stress-exposed mouse into a non-exposed mouse can recapitulate the effects observed in the exposed’s offspring, suggesting serum factors can modify the composition of heritable molecules in the germ cells and affect offspring traits (50).
We utilized an outcrossing scheme to demonstrate that the observed defects persisted when passed only through the maternal line. In this mating scheme, 100% of the observed F2 outcrossed female mice had evidence of hepatic metabolic dysregulation. This pattern of inheritance is not consistent with in-cis mechanisms such as an unaccounted DNA variant, DNA methylation, or histone protein modification, because in these scenarios the alleles in the Khdc3-null maternal grandmother would be present in ∼50% of the F2 outcrossed females.
We chose to examine females because of the reported fertility phenotypes in KHDC3 human (16, 17) and Khdc3 mouse female mutants (18). Human KHDC3-null females are infertile, while mouse Khdc3-null females are subfertile. The fertility defects are associated with oocyte DNA methylation abnormalities that we suspect are sequelae of a more primary defect in mutant organisms. There are no reports of abnormal phenotypes in Khdc3-null males, however our detected metabolic abnormalities could affect males as well, which will be an important future investigation. Furthermore, it will be important to examine how Khdc3-null mice transmit exposure-based information across generations to affect phenotypes in non-exposed offspring, as has been observed in numerous studies, most thoroughly with a high fat diet or stress.
One important conclusion from studies on genetic nurture is that wild type mice descended from mutant ancestors should be used thoughtfully and with caution, because phenotypes can be obscured that would be apparent with the use of wild type mice without any ancestral history of mutation. In sum, the effect of non-inherited genetic variants on offspring phenotypes remains a relatively unexplored and potentially significant contributor to traits and disease risk.
Materials and Methods Mice
FVB/N WT (Jax 001800) mice were purchased from The Jackson Laboratory and bred in-house. Frozen Khdc3-null sperm was obtained from the Mutant Mouse Resource and Research Center (MMRRC, North Carolina) and the Khdc3 KO mouse was rederived. Female mice were used at 8 weeks of age. In individual experiments, all animals were age-matched. All mice were maintained under specific pathogen-free (SPF) conditions on a 12-hour light/dark cycle, and provided food and water ad libitum. All mouse experiments were approved by, and performed in accordance with, the Institutional Animal Care and Use Committee guidelines at Weill Cornell Medicine.
Mouse Genotyping
As described in this study, WT mice were generated from WT parents. WT* and KO* mice were generated from the mating of Khdc3 heterozygote parents that were themselves generated from mating of WT and KOKO mice. KOKO mice were generated from Khdc3 KO parents. DNA was extracted from tail biopsies by incubating the tails biopsies in an alkaline lysis buffer (24mM NaOH and 0.2mM disodium EDTA) for 30 minutes at 95°C and then on ice for 10 minutes. A neutralization buffer (40 mM Tris-HCl) was added to the samples and diluted 1:10 with water. PCR to detect Khdc3 deletion was performed with 50–200 ng of DNA using the Phire Green Hot Start II DNA Polymerase Master Mix (Thermo Scientific) in the Bio-Rad T100 Thermal Cycler. PCR for Khdc3 deletion was performed using a forward primer for the wild type allele F1, a separate forward primer that incorporates the deletion allele F2, and a common reverse primer R1. The primers were as follows: P1: 5’-TGCCTGGGCAGGTTATTTAG-3’, P2: 5’-CGAGCGTCTGAAACCTCTTC-3’ and P3: 5’-AGCTAGCTTGGCTGGACGTA-3’. P1 and P2 amplified the wild type allele and P1 and P3 amplified the Khdc3 mutant allele (KO). The PCR products were separated on a 1% agarose gel with SYBR safe DNA gel stain. Because of the different sizes of the PCR products, the genotypes can be easily determined from the band patterns on DNA gels.
Liver RNA Extraction
Total RNA extraction was performed on the liver of 8 week old mice using the Qiagen RNeasy Lipid Tissue Mini Kit (Qiagen) according to the manufacturer’s instructions. Samples were eluted in 30μl nuclease-free water. Nucleic acid concentration, A260/280 and A260/230 ratios was determined via NanodropD-1000 (Thermo Scientific). Extracted samples were aliquoted and stored at −80°C.
RT-PCR and quantitative RT-PCR
2000ng RNA per sample was used to generate cDNA using SuperScript III First-Strand Synthesis System (Invitrogen) following the manufacturer’s instructions. cDNA was diluted to 1:4 in water before performing qPCR using SYBR Green PCR Master Mix (Applied Biosystems). qPCR reactions were performed on a QuantStudio 6 Flex Real Time PCR Instrument. Cycling conditions were as follows: Initial denaturation 95°C for 3 minutes, 40 cycles of denaturation at 95°C for 15 seconds followed by annealing/extension at 60°C for 60 seconds. Relative expression was calculated using the ΔΔCt method, using YWHAZ as the reference gene. Primers used are listed in Supplemental Table 1.
Oocyte RNA Isolation
The oocytes were dissected following a published protocol with some modifications (51). Briefly, the ovaries were first dissected from unstimulated 8-week-old WT and Khdc3-null mice and placed in PBS. The ovary was dissected from surrounding para-ovarian fat and subsequently placed in a 35-mm culture dish with 2mL PBS and 20μl collagenase and 20μl DNase in a 37°C incubator for 20 minutes. Using a dissecting microscope, the oocytes were separated from granulosa cells, theca cells, and stromal cells, and transferred to a petri dish with PBS using a mouth-controlled micropipette. RNA was extracted from the oocytes using the PicoPure RNA Isolation kit (Thermo Fisher) following the manufacturer’s instructions and eluted in 20 µL of the provided elution buffer.
Liver RNA-seq Library Preparation
Sequencing libraries were prepared using the Illumina TruSeq Stranded Total RNA kit according to the manufacturer’s protocol and sequenced to a depth of 40 million reads per sample. The paired-end (PE) libraries were sequenced on Novaseq platform.
Oocyte RNA-seq Library Preparation
Sequencing libraries were prepared using the SMART-Seq v4 Ultra Low Input RNA plus Nextera XT DNA Sample Preparation according to the manufacturer’s protocol and sequenced to a depth of 40 million reads per sample. The paired-end (PE) libraries were sequenced on Novaseq platform.
RNA-seq Bioinformatics
All sequenced reads were assessed for quality using FastQC (52). Adapter trimming and filtering of low-quality bases (<20) were performed using cutadapt (53). The trimmed reads were then aligned with Hisat2 (54) against the mouse reference (GRCm38) genome. Per gene counts were computed using htseq-count (55) (genocode M15). DESeq2 (56) bioconductor package was used for normalization and differential gene expression analysis.
Oocyte small RNA isolation
Eight week old female mice were sacrificed and 30 oocytes were dissected under stereomicroscopy from the ovaries of each mouse, as described above. Four separate WT, Khdc3-null, and WT(P)** RNA samples were generated from separate mice. Total RNA was isolated using a magnet-based method (ChargeSwitch Total RNA Cell Kit).
Oocyte small RNA-Seq Library Preparation
Small RNA quality and concentration was assessed using an Agilent Bioanalyzer and small RNA Chip to ensure a clear small RNA population of more than 1ng. Subsequently, small RNA libraries were generated and sequenced to a depth of 10 million reads per sample. Small RNA libraries were generated using the NEXTFLEX Small RNA-seq Kit, which were sequenced on the Novaseq platform.
Small RNA-seq Bioinformatics
Read quality was assessed using FastQC (v0.11.8), and adapter sequences were trimmed using Trimmomatic (v0.39). After adapter trimming, reads were mapped sequentially to rRNA mapping reads, miRbase, murine tRNAs, pachytene piRNA clusters (57), repeatmasker and Refseq using Bowtie 2 alignment algorithm (v2.3.5) and totaled using Feature Counts on the Via Foundry (v1.6.4) platform (58). To assess for differentially expressed small RNAs, data was loaded into R Statistical Software and analyzed using the DESeq2 package (56). Differentially abundant small RNAs were determined as those with a log fold change > 0.58 and P-value < 0.05.
Serum preparation
Blood was collected from the submandibular (facial) vein of mice. Approximately 400-500μl of whole blood was collected from each mouse and placed in an SST amber microtainer blood collection tube (BD). Blood collection tubes were placed in a 37°C incubator for 30 minutes, followed by 10 minutes at 4°C. Samples were then centrifuged at 5500 RPM for 10 minutes at 4°C. The serum was collected from the top of the tube and stored at −80°C until needed.
Metabolomic profiling
Metabolomics profiling was conducted using ultra-high-performance liquid chromatography-tandem mass-spectrometry by the metabolomics provider Metabolon Inc. (Morrisville, USA) on mouse serum samples from WT, WT**(P), WT****(P) and KOKO mice, as well as these mice after consuming a high fat diet. The metabolomic dataset measured by Metabolon includes known metabolites containing the following broad categories – amino-acids, peptides, carbohydrates, energy intermediates, lipids, nucleotides, cofactors and vitamins, and xenobiotics. Statistical differences were determined using unpaired t-test.
High fat Diet
Female WT and WT****(P) mice were divided into two diet groups, one group receiving a high-fat diet (HFD, D12492i; Research Diets Inc., New Brunswick, NJ) and the other group received a normal diet with 10% fat (D12450Ji; Research Diets Inc., New Brunswick, NJ). All mice were given access to food and water ad libitum and were maintained on a 12:12-h light-dark artificial lighting cycle. After 8 weeks of each diet, serum was collected and stored in the −80°C.
Statistical Analysis and data visualization
Gene ontology performed with WebGestalt (https://www.webgestalt.org/) using KEGG analysis. Volcano plots generated with ggVolcanoR (https://ggvolcanor.erc.monash.edu). Venn diagrams generated with InteractiVenn (http://www.interactivenn.net/). Scatter plots were generated with Prism 9.2.0. The manuscript was written using Microsoft Word v16.66.1. Figures were assembled on Microsoft PowerPoint v16.66.1
Data Availability
All RNA-Seq and small RNA-Seq data are available in the NCBI GEO database.
Acknowledgements
This work was supported by the Young Investigator Award (LS) at Weill Cornell Medicine, Pediatric Scientist Development Program under the National Institute of Child Health under award number K12HD000850 (SC), Friedman Family Clinical Scholar Award (MSR), and National Institute of Child Health and Human Development (NICHD) of the National Institutes of Health under award number P50HD104454 (MSR).
Additional Information
Supplementary Information is available for this paper.
Supplemental Material
References
- 1.Sperm tsRNAs contribute to intergenerational inheritance of an acquired metabolic disorderScience 351:397–400
- 2.Paternal Contributions to Offspring Health: Role of Sperm Small RNAs in Intergenerational Transmission of Epigenetic InformationFront Cell Dev Biol 7
- 3.Longevity determined by paternal ancestors’ nutrition during their slow growth periodActa Biotheor 49:53–59
- 4.Cardiovascular and diabetes mortality determined by nutrition during parents’ and grandparents’ slow growth periodEur J Hum Genet 10:682–688
- 5.Small RNAs Gained during Epididymal Transit of Sperm Are Essential for Embryonic Development in MiceDev Cell 46:470–480
- 6.Parent-of-origin effects of A1CF and AGO2 on testicular germ-cell tumors, testicular abnormalities, and fertilization biasProc Natl Acad Sci U S A 113:E5425–5433
- 7.Trans-generational epistasis between Dnd1Ter and other modifier genes controls susceptibility to testicular germ cell tumorsHum Mol Genet 16:2233–2240
- 8.Reduced Sensitivity to Thyroid Hormone as a Transgenerational Epigenetic Marker Transmitted Along the Human Male LineThyroid 29:778–782
- 9.Ancestral paternal genotype controls body weight and food intake for multiple generationsHum Mol Genet 19:4134–4144
- 10.Transgenerational genetic effects of the paternal Y chromosome on daughters’ phenotypesEpigenomics 2:513–521
- 11.Insulin VNTR allele-specific effect in type 1 diabetes depends on identity of untransmitted paternal allele. The IMDIAB GroupNat Genet 17:350–352
- 12.Mutation in folate metabolism causes epigenetic instability and transgenerational effects on developmentCell 155:81–93
- 13.The nature of nurture: Effects of parental genotypesScience 359:424–428
- 14.From paramutation to paradigmPLoS Genet 9
- 15.Breeding scheme and maternal small RNAs affect the efficiency of transgenerational inheritance of a paramutation in miceSci Rep 5
- 16.A KHDC3L mutation resulting in recurrent hydatidiform mole causes genome-wide DNA methylation loss in oocytes and persistent imprinting defects post-fertilisationGenome Med 11
- 17.The subcortical maternal complex: multiple functions for one biological structure?J Assist Reprod Genet 33:1431–1438
- 18.Role of Filia, a maternal effect gene, in maintaining euploidy during cleavage-stage mouse embryogenesisProc Natl Acad Sci U S A 106:7473–7478
- 19.PPARs as Metabolic Regulators in the Liver: Lessons from Liver-Specific PPAR-Null MiceInt J Mol Sci 21
- 20.It’s reticulated: the liver at the heart of atherosclerosisJ Endocrinol 238:R1–R11
- 21.The Pathophysiological Role of CoAInt J Mol Sci 21
- 22.Regulation of glucose metabolism from a liver-centric perspectiveExp Mol Med 48
- 23.Steroidogenic control of liver metabolism through a nuclear receptor-networkMol Metab 30:221–229
- 24.Mutated CYP17A1 promotes atherosclerosis and early-onset coronary artery diseaseCell Commun Signal 21
- 25.Liver lipid metabolismJ Anim Physiol Anim Nutr (Berl) 92:272–283
- 26.Bile acid metabolism and signalingCompr Physiol 3:1191–1212
- 27.Xanthosine, a purine glycoside mediates hepatic glucose homeostasis through inhibition of gluconeogenesis and activation of glycogenesis via regulating the AMPK/ FoxO1/AKT/GSK3β signaling cascadeChem Biol Interact 371
- 28.Tauroursodeoxycholate-Bile Acid with Chaperoning Activity: Molecular and Cellular Effects and Therapeutic PerspectivesCells 8
- 29.The Natural Metabolite 4-Cresol Improves Glucose Homeostasis and Enhances β-Cell FunctionCell Rep 30:2306–2320
- 30.Role of xenobiotics in the induction and progression of fatty liver diseaseToxicol Res (Camb) 7:664–680
- 31.Increased serum levels of quinolinic acid indicate enhanced severity of hepatic dysfunction in patients with liver cirrhosisHum Immunol 74:60–66
- 32.Endogenous nicotinamide riboside metabolism protects against diet-induced liver damageNat Commun 10
- 33.Nicotinamide N-methyltransferase and liver diseasesGenes Dis 10:1883–1893
- 34.Daidzein supplementation prevents non-alcoholic fatty liver disease through alternation of hepatic gene expression profiles and adipocyte metabolismInt J Obes (Lond) 35:1019–1030
- 35.Influence of Genistein on Hepatic Lipid Metabolism in an In Vitro Model of Hepatic SteatosisMolecules 26
- 36.Ferulic acid protects HepG2 cells and mouse liver from iron-induced damageBiochem Biophys Rep 35
- 37.Dysfunctional tRNA reprogramming and codon-biased translation in cancerTrends Mol Med 28:964–978
- 38.Angiogenin-induced tRNA fragments inhibit translation initiationMol Cell 43:613–623
- 39.NSUN6 is a human RNA methyltransferase that catalyzes formation of m5C72 in specific tRNAsRNA 21:1532–1543
- 40.ELAC1 Repairs tRNAs Cleaved during Ribosome-Associated Quality ControlCell Rep 30:2106–2114
- 41.METTL1 promotes tumorigenesis through tRNA-derived fragment biogenesis in prostate cancerMol Cancer 22
- 42.METTL1 Promotes let-7 MicroRNA Processing via m7G MethylationMol Cell 74:1278–1290
- 43.Overview of MicroRNA BiogenesisMechanisms of Actions, and Circulation. Front Endocrinol (Lausanne) 9
- 44.The RNA-Binding Protein SYNCRIP Is a Component of the Hepatocyte Exosomal Machinery Controlling MicroRNA SortingCell Rep 17:799–808
- 45.Paternal High-Fat Diet Altered Sperm 5’tsRNA-Gly-GCC Is Associated With Enhanced Gluconeogenesis in the OffspringFront Mol Biosci 9
- 46.Biogenesis and function of tRNA fragments during sperm maturation and fertilization in mammalsScience 351:391–396
- 47.miR-107 orchestrates ER stress induction and lipid accumulation by post-transcriptional regulation of fatty acid synthase in hepatocytesBiochim Biophys Acta 1839:334–343
- 48.Intergenerational epigenetic inheritance of cancer susceptibility in mammalsElife 8
- 49.Epigenetic inheritance of acquired traits through sperm RNAs and sperm RNA modificationsNat Rev Genet 17:733–743
- 50.Involvement of circulating factors in the transmission of paternal experiences through the germlineEMBO J 39
- 51.Mammalian Oocytes Locally Remodel Follicular Architecture to Provide the Foundation for Germline-Soma CommunicationCurr Biol 28:1124–1131
- 52.S. Andrews (2010) FastQC: a quality control tool for high throughput sequence data. P Available online at: http://www.bioinformatics.babraham.ac.uk/projects/fastqc.
- 53.Cutadapt removes adapter sequences from high-throughput sequencing readsEMBnet J
- 54.Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotypeNat Biotechnol 37:907–915
- 55.HTSeq--a Python framework to work with high-throughput sequencing dataBioinformatics 31:166–169
- 56.Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2Genome Biol 15
- 57.An ancient transcription factor initiates the burst of piRNA production during early meiosis in mouse testesMol Cell 50:67–81
- 58.DolphinNext: a distributed data processing platform for high throughput genomicsBMC Genomics 21
Article and author information
Author information
Version history
- Sent for peer review:
- Preprint posted:
- Reviewed Preprint version 1:
Copyright
© 2024, Senaldi et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
- views
- 93
- download
- 1
- citations
- 0
Views, downloads and citations are aggregated across all versions of this paper published by eLife.