Paternal aging impacts expression and epigenetic markers as early as the first embryonic tissue lineage differentiation

BACKGROUND

Advanced paternal age (APA) is associated with adverse outcomes to offspring health, including increased risk for neurodevelopmental disorders. The aim of this study was to investigate the methylome and transcriptome of the first two early embryonic tissue lineages, the inner cell mass (ICM) and the trophectoderm (TE), from human blastocysts in association with paternal age and disease risk. High quality human blastocysts were donated with patient consent from donor oocyte IVF cycles from either APA (≥ 50 years) or young fathers. Blastocysts were mechanically separated into ICM and TE lineage samples for both methylome and transcriptome analyses.

RESULTS

Significant differential methylation and transcription was observed concurrently in ICM and TE lineages of APA-derived blastocysts compared to those from young fathers. The methylome revealed significant enrichment for neuronal signaling pathways, as well as an association with neurodevelopmental disorders and imprinted genes, largely overlapping within both the ICM and TE lineages. Significant enrichment of neurodevelopmental signaling pathways was also observed for differentially expressed genes, but only in the ICM. In stark contrast, no significant signaling pathways or gene ontology terms were identified in the trophectoderm. Despite normal semen parameters in aged fathers, these significant molecular alterations can adversely contribute to downstream impacts on offspring health, in particular neurodevelopmental disorders like autism spectrum disorder and schizophrenia.

CONCLUSIONS

An increased risk for neurodevelopmental disorders is well described in children conceived by aged fathers. Using blastocysts derived from donor oocyte IVF cycles to strategically control for maternal age, our data reveals evidence of methylation dysregulation in both tissue lineages, as well as transcription dysregulation in neurodevelopmental signaling pathways associated with APA fathers. This data also reveals that embryos derived from APA fathers do not appear to be compromised for initial implantation potential with no significant pathway signaling disruption in trophectoderm transcription. Collectively, our work provides insights into the complex molecular mechanisms that occur upon paternal aging during the first lineage differentiation in the preimplantation embryo. Early expression and epigenetic markers of APA-derived preimplantation embryos highlight the susceptibility of the future fetus to adverse health outcomes.

Identifying risk variants for embryo aneuploidy using ultra-low coverage whole-genome sequencing from preimplantation genetic testing

Aneuploidy frequently arises during human meiosis and is the primary cause of early miscarriage and in vitro fertilization (IVF) failure. Individuals undergoing IVF exhibit significant variability in aneuploidy rates, although the exact genetic causes of the variability in aneuploid egg production remain unclear. Preimplantation genetic testing for aneuploidy (PGT-A) using next-generation sequencing is a standard test for identifying and selecting IVF-derived euploid embryos. The wealth of embryo aneuploidy data and ultra-low coverage whole-genome sequencing (ulc-WGS) data from PGT-A have the potential to discover variants in parental genomes that are associated with aneuploidy risk in their embryos. Using ulc-WGS data from ∼10,000 PGT-A biopsies, we imputed genotype likelihoods of genetic variants in embryo genomes. We then used the imputed variants and embryo aneuploidy calls to perform a genome-wide association study of aneuploidy incidence. Finally, we carried out functional evaluation of the identified candidate gene in a mouse oocyte system. We identified one locus on chromosome 3 that is significantly associated with meiotic aneuploidy risk. One candidate gene, CCDC66, encompassed by this locus, is involved in chromosome segregation during meiosis. Using mouse oocytes, we showed that CCDC66 regulates meiotic progression and chromosome segregation fidelity, especially in older mice. Our work extended the research utility of PGT-A ulc-WGS data by allowing robust association testing and improved the understanding of the genetic contribution to maternal meiotic aneuploidy risk. Importantly, we introduce a generalizable method that has potential to be leveraged for similar association studies that use ulc-WGS data.

Identifying risk genes for embryo aneuploidy using ultra-low coverage whole-genome sequencing

Background

Aneuploidy, the state of a cell containing extra or missing chromosomes, frequently arises during human meiosis and is the primary cause of early miscarriage and maternal age-related in vitro fertilization (IVF) failure. IVF patients exhibit significant variability in aneuploidy rates, although the exact genetic causes of the variability in aneuploid egg production remain unclear. Preimplantation genetic testing for aneuploidy (PGT-A) using ultra-low coverage whole-genome sequencing (ulc-WGS) is a standard test for identifying and selecting IVF-derived embryos with a normal chromosome complement. The wealth of embryo aneuploidy data and ulc-WGS data from PGT-A has potential for discovering variants in paternal genomes that are associated with aneuploidy risk in their embryos.

Methods

Using ulc-WGS data from ∼10,000 PGT-A biopsies, we imputed genotype likelihoods of genetic variants in parental genomes. We then used the imputed variants and aneuploidy calls from the embryos to perform a genome-wide association study of aneuploidy incidence. Finally, we carried out functional evaluation of the identified candidate gene in a mouse oocyte system.

Results

We identified one locus on chromosome 3 that is significantly associated with maternal meiotic aneuploidy risk. One candidate gene, CCDC66, encompassed by this locus, is involved in chromosome segregation during meiosis. Using mouse oocytes, we showed that CCDC66 regulates meiotic progression and chromosome segregation fidelity, especially in older mice.

Conclusions

Our work extended the research utility of PGT-A ulc-WGS data by allowing robust association testing and improved the understanding of the genetic contribution to maternal meiotic aneuploidy risk. Importantly, we introduce a generalizable method that can be leveraged for similar association studies using ulc-WGS data.

Segmental aneuploid hotspots identified across the genome concordant on reanalysis

The aim of this study was to characterize a large set of full segmental aneuploidies identified in trophectoderm (TE) biopsies and evaluate concordance in human blastocysts. Full segmental aneuploid errors were identified in TE biopsies (n = 2766) from preimplantation genetic testing for aneuploid (PGT-A) cycles. Full segmental deletions (n = 1872; 66.1%) presented twice as many times as duplications (n = 939; 33.9%), mapped more often to the q-arm (n = 1696; 61.3%) than the p-arm (n = 847; 31.0%) or both arms (n = 223; 8.1%; P < 0.05), and were eight times more likely to include the distal end of a chromosome than not (P < 0.05). Additionally, 37 recurring coordinates (each ≥ 10 events) were discovered across 17 different chromosomes, which were also significantly enriched for distal regions (P = 4.1 × 10-56). Blinded concordance analysis of 162 dissected blastocysts validated the original TE PGT-A full segmental result for a concordance of 96.3% (n = 156); remaining dissected blastocysts were identified as mosaic (n = 6; 3.7%). Origin of aneuploid analysis revealed full segmental aneuploid errors were mostly paternally derived (67%) in contrast to whole chromosome aneuploid errors (5.8% paternally derived). Errors from both parental gametes were observed in 6.5% of aneuploid embryos when multiple whole chromosomes were affected. The average number of recombination events was significantly less in paternally derived (1.81) compared to maternally derived (3.81) segmental aneuploidies (P < 0.0001). In summary, full segmental aneuploidies were identified at hotspots across the genome and were highly concordant upon blinded analysis. Nevertheless, future studies assessing the reproductive potential of full (non-mosaic) segmental aneuploid embryos are critical to rule out potential harmful reproductive risks.

Spatially resolved transcriptomic profiling of ovarian aging in mice

Ovarian aging precedes that of any other mammalian organ and is the primary cause of female age-related infertility. The biological mechanisms responsible for ovarian aging remain unclear. Previous studies have been limited by their use of bulk RNA-sequencing, which masks the dynamic and heterogeneous nature of the ovary. In this study, we spatially resolved the transcriptomic landscape of ovaries from young and aged outbred mice. In total, we defined eight main ovarian cell populations, all of which were characterized by significant transcriptomic changes between young and aged samples. Further sub-cluster analysis revealed separate transcriptomes for distinct granulosa cell populations found in young versus aged mice, in addition to an oocyte sub-cluster population completely absent from aged mouse ovaries. This study provides a new perspective on mammalian ovarian aging using spatial transcriptomics to achieve deeper understanding of the localization and cell-population-specific mechanisms underlying age-related fertility decline.

•

Spatial transcriptomic analysis of murine ovarian aging characterizes cell populations

•

Distinct granulosa cell populations exist in aged mice compared to young

•

These are associated with the disruption of inhibins, activins, and gap junction activity

•

Oocyte transcriptomes captured in vivo and within follicles of varying maturity

The inherited methylome landscape is directly altered with paternal aging and associated with offspring neurodevelopmental disorders

Paternal aging and the prevalence of neurodevelopmental disorders in offspring are well documented. Yet, the underlying mechanism and the mode of inheritance have not been conclusively established. Advancing paternal age is a subtle and varying phenotype. As such, it is likely that a threshold for cumulative risk may exist that, if surpassed, culminates in a predisposition to disease and ultimately an observed phenotype in offspring. Epigenetic regulation provides a plausible explanation for the nongenetic paternal transmission of disease susceptibility. With the use of whole‐genome methylation sequencing, the data described herein substantiate an increasingly compromised DNA methylation profile as sperm ages and, for the first time, also demonstrate a generational correlation in sperm and blastocyst of an altered methylome associated with advanced paternal age. Methylation alterations are not randomly distributed across the genome, but appear clustered at certain chromosomal locations, and significantly colocalize with regions of nucleosome retention. Genes associated with autism spectrum disorder, schizophrenia, and bipolar disorder are significantly enriched with causative methylation aberrations in both sperm and embryos from aged fathers. The long‐term health burden and societal economic impact of these conditions are substantial and will continue with increasingly prevalent diagnosis. This work provides a mechanistic link between the paternal age effect and offspring neurodevelopmental disorders leading to a better understanding of causation and investigation into potential future therapy.

Genetic Risk of Arrhythmic Phenotypes in Patients With Dilated Cardiomyopathy

BACKGROUND

Genotype-phenotype correlations in dilated cardiomyopathy (DCM) and, in particular, the effects of gene variants on clinical outcomes remain poorly understood.

OBJECTIVES

The purpose of this study was to investigate the prognostic role of genetic variant carrier status in a large cohort of DCM patients.

METHODS

A total of 487 DCM patients were analyzed by next-generation sequencing and categorized the disease genes into functional gene groups. The following composite outcome measures were assessed: 1) all-cause mortality; 2) heart failure-related death, heart transplantation, or destination left ventricular assist device implantation (DHF/HTx/VAD); and 3) sudden cardiac death/sustained ventricular tachycardia/ventricular fibrillation (SCD/VT/VF).

RESULTS

A total of 183 pathogenic/likely pathogenic variants were found in 178 patients (37%): 54 (11%) Titin; 19 (4%) Lamin A/C (LMNA); 24 (5%) structural cytoskeleton-Z disk genes; 16 (3.5%) desmosomal genes; 46 (9.5%) sarcomeric genes; 8 (1.6%) ion channel genes; and 11 (2.5%) other genes. All-cause mortality was no different between variant carriers and noncarriers (p = 0.99). A trend toward worse SCD/VT/VF (p = 0.062) and DHF/HTx/VAD (p = 0.061) was found in carriers. Carriers of desmosomal and LMNA variants experienced the highest rate of SCD/VT/VF, which was independent of the left ventricular ejection fraction.

CONCLUSIONS

Desmosomal and LMNA gene variants identify the subset of DCM patients who are at greatest risk for SCD and life-threatening ventricular arrhythmias, regardless of the left ventricular ejection fraction.

Transcriptome signature of ventricular arrhythmia in dilated cardiomyopathy reveals increased fibrosis and activated TP53

AIMS

One-third of DCM patients experience ventricular tachycardia (VT), but a clear biological basis for this has not been established. The purpose of this study was to identify transcriptome signatures and enriched pathways in the hearts of dilated cardiomyopathy (DCM) patients with VT.

METHODS AND RESULTS

We used RNA-sequencing in explanted heart tissue from 49 samples: 19 DCM patients with VT, 16 DCM patients without VT, and 14 non-failing controls. We compared each DCM cohort to the controls and identified the genes that were differentially expressed in DCM patients with VT but not without VT. Differentially expressed genes were evaluated using pathway analysis, and pathways of interest were investigated by qRT-PCR validation, Western blot, and microscopy. There were 590 genes differentially expressed in DCM patients with VT that are not differentially expressed in patients without VT. These genes were enriched for genes in the TGFß1 and TP53 signaling pathways. Increased fibrosis and activated TP53 signaling was demonstrated in heart tissue of DCM patients with VT.

CONCLUSIONS

Our study supports that distinct biological mechanisms distinguish ventricular arrhythmia in DCM patients.

The Giant Protein Titin's Role in Cardiomyopathy: Genetic, Transcriptional, and Post-translational Modifications of TTN and Their Contribution to Cardiac Disease

Dilated cardiomyopathy (DCM) is a leading cause of heart failure, sudden cardiac death and heart transplant. DCM is inherited in approximately 50% of cases, in which the most frequent genetic defects are truncation variants of the titin gene (TTNtv). TTN encodes titin, which is the largest protein in the body and is an essential component of the sarcomere. Titin serves as a biological spring, spanning half of the sarcomere and connecting the Z-disk to the M-line, with scaffold and signaling functions. Truncations of titin are believed to lead to either haploinsufficiency and loss-of-function, or to a “poison peptide” effect. However, other titin mechanisms are postulated to influence cardiac function including post-translational modifications, in particular changes in titin phosphorylation that alters the stiffness of the protein, and diversity of alternative splicing that generates different titin isoforms. In this article, we review the role of TTN mutations in development of DCM, how differential expression of titin isoforms relate to DCM pathophysiology, and discuss how post-translational modifications of titin can affect cardiomyocyte function. Current research efforts aim to elucidate the contribution of titin to myofibril assembly, stability, and signal transduction, and how mutant titin leads to cardiac dysfunction and human disease. Future research will need to translate this knowledge toward novel therapeutic approaches that can modulate titin transcriptional and post-translational defects to treat DCM and heart failure.

Autoantigen glycoprotein 70 expression is regulated by a single locus, which acts as a checkpoint for pathogenic anti-glycoprotein 70 autoantibody production and hence for the corresponding development of severe nephritis, in lupus-prone PXSB mice

Retroviral envelope glycoprotein gp70 is present in the sera of immunologically normal and autoimmune-prone strains of mice. However, only lupus-prone mice spontaneously develop gp70-anti-gp70 immune complexes (gp70IC), and these have been implicated in the development of nephritis. We investigated the genetic factors that affect the production of both free serum gp70 and gp70IC in the lupus-prone BXSB mouse strain by analyzing (BXSB x (C57BL/10 x BXSB)F(1))- and (C57BL/10 x (C57BL/10 x BXSB)F(1))-backcrossed male mice. Production of gp70 mapped to a single major locus located on chromosome 13 (Bxs6) with a maximum log likelihood of the odds of 36.7 (p = 1.6 x 10(-38)). The level of gp70IC was highly dependent on Bxs6-related gp70 production, and high titer autoantibody production only occurred when serum gp70 levels were greater than a threshold value of approximately 4.0 microg/ml. The subdivision of the (BXSB x (C57BL/10 x BXSB)F(1))-backcrossed mice into those homozygous or heterozygous for Bxs6 enabled a remarkable association to be observed between high levels of gp70IC and severe nephritis in the Bxs6 homozygote population. A further mapping study in these two subgroups identified a previously unrecognized interval associated with the production of autoantibodies.

Identification of intervals on chromosomes 1, 3, and 13 linked to the development of lupus in BXSB mice

OBJECTIVE

To identify intervals containing systemic lupus erythematosus (SLE) susceptibility alleles in the BXSB strain of mice.

METHODS

We analyzed 286 (B10 x [B10 x BXSB]F1) backcross mice for a range of phenotypic traits associated with the development of SLE in BXSB mice. The mice were genotyped using 93 microsatellite markers, and the linkage of these markers to disease was studied by extreme-phenotype and quantitative trait locus analysis.

RESULTS

The disease phenotype in these backcross mice was less severe than that in BXSB mice. However, antinuclear antibody production was increased compared with the parental strain. We identified 4 areas of genetic linkage to disease on chromosome 1 (Bxs1-4), 1 on chromosome 3 (Bxs5), and another interval on chromosome 13 which were associated with various aspects of the phenotype. Bxs4 and Bxs5 are located in regions not previously linked to disease in other models of SLE.

CONCLUSION

SLE in the BXSB mouse model has a complex genetic basis and involves at least 5 distinct intervals located on chromosomes 1 and 3. There is evidence that different intervals affect particular aspects of the SLE phenotype.