Centromeric transposable elements and epigenetic status drive karyotypic variation in the eastern hoolock gibbon.

Great apes have maintained a stable karyotype with few large-scale rearrangements; in contrast, gibbons have undergone a high rate of chromosomal rearrangements coincident with rapid centromere turnover. Here we characterize assembled centromeres in the Eastern hoolock gibbon, Hoolock leuconedys (HLE), finding a diverse group of transposable elements (TEs) that differ from the canonical alpha satellites found across centromeres of other apes. We find that HLE centromeres contain a CpG methylation centromere dip region, providing evidence this epigenetic feature is conserved in the absence of satellite arrays; nevertheless, we report a variety of atypical centromeric features, including protein-coding genes and mismatched replication timing. Further, large structural variations define HLE centromeres and distinguish them from other gibbons. Combined with differentially methylated TEs, topologically associated domain boundaries, and segmental duplications at chromosomal breakpoints, we propose that a "perfect storm" of multiple genomic attributes with propensities for chromosome instability shaped gibbon centromere evolution.

Complete sequencing of ape genomes.

We present haplotype-resolved reference genomes and comparative analyses of six ape species, namely: chimpanzee, bonobo, gorilla, Bornean orangutan, Sumatran orangutan, and siamang. We achieve chromosome-level contiguity with unparalleled sequence accuracy (<1 error in 500,000 base pairs), completely sequencing 215 gapless chromosomes telomere-to-telomere. We resolve challenging regions, such as the major histocompatibility complex and immunoglobulin loci, providing more in-depth evolutionary insights. Comparative analyses, including human, allow us to investigate the evolution and diversity of regions previously uncharacterized or incompletely studied without bias from mapping to the human reference. This includes newly minted gene families within lineage-specific segmental duplications, centromeric DNA, acrocentric chromosomes, and subterminal heterochromatin. This resource should serve as a definitive baseline for all future evolutionary studies of humans and our closest living ape relatives.

Astrocyte extracellular matrix modulates neuronal dendritic development.

Major developmental events occurring in the hippocampus during the third trimester of human gestation and neonatally in altricial rodents include rapid and synchronized dendritic arborization and astrocyte proliferation and maturation. We tested the hypothesis that signals sent by developing astrocytes to developing neurons modulate dendritic development in vivo. We altered neuronal development by neonatal (third trimester-equivalent) ethanol exposure in mice; this treatment increased dendritic arborization in hippocampal pyramidal neurons. We next assessed concurrent changes in the mouse astrocyte translatome by translating ribosomal affinity purification (TRAP)-seq. We followed up on ethanol-inhibition of astrocyte Chpf2 and Chsy1 gene translation because these genes encode for biosynthetic enzymes of chondroitin sulfate glycosaminoglycan (CS-GAG) chains (extracellular matrix components that inhibit neuronal development and plasticity) and have not been explored before for their roles in dendritic arborization. We report that Chpf2 and Chsy1 are enriched in astrocytes and their translation is inhibited by ethanol, which also reduces the levels of CS-GAGs measured by Liquid Chromatography/Mass Spectrometry. Finally, astrocyte-conditioned medium derived from Chfp2-silenced astrocytes increased neurite branching of hippocampal neurons in vitro. These results demonstrate that CS-GAG biosynthetic enzymes in astrocytes regulates dendritic arborization in developing neurons.

Deletion of an evolutionarily conserved TAD boundary compromises spermatogenesis in mice.

Spermatogenesis is a complex process that can be disrupted by genetic and epigenetic changes, potentially leading to male infertility. Recent research has rapidly increased the number of protein coding mutations causally linked to impaired spermatogenesis in humans and mice. However, the role of non-coding mutations remains largely unexplored. As a case study to evaluate the effects of non-coding mutations on spermatogenesis, we first identified an evolutionarily conserved topologically associated domain (TAD) boundary near two genes with important roles in mammalian testis function: Dmrtb1 and Lrp8 . We then used CRISPR-Cas9 to generate a mouse line where 26kb of the boundary was removed including a strong and evolutionarily conserved CTCF binding site. ChIP-seq and Hi-C experiments confirmed the removal of the CTCF site and a resulting increase in the DNA-DNA interactions across the domain boundary. Mutant mice displayed significant changes in testis gene expression, abnormal testis histology, a 35% drop in the estimated efficiency of spermatogenesis and a 28% decrease in daily sperm production compared to littermate controls. Despite these quantitative changes in testis function, mutant mice show no significant changes in fertility. This suggests that non-coding deletions affecting testis gene regulation may have smaller effects on fertility compared to coding mutations of the same genes. Our results demonstrate that disruption of a TAD boundary can have a negative impact on sperm production and highlight the importance of considering non-coding mutations in the analysis of patients with male infertility.

The complete sequence and comparative analysis of ape sex chromosomes.

Apes possess two sex chromosomes-the male-specific Y chromosome and the X chromosome, which is present in both males and females. The Y chromosome is crucial for male reproduction, with deletions being linked to infertility1. The X chromosome is vital for reproduction and cognition2. Variation in mating patterns and brain function among apes suggests corresponding differences in their sex chromosomes. However, owing to their repetitive nature and incomplete reference assemblies, ape sex chromosomes have been challenging to study. Here, using the methodology developed for the telomere-to-telomere (T2T) human genome, we produced gapless assemblies of the X and Y chromosomes for five great apes (bonobo (Pan paniscus), chimpanzee (Pan troglodytes), western lowland gorilla (Gorilla gorilla gorilla), Bornean orangutan (Pongo pygmaeus) and Sumatran orangutan (Pongo abelii)) and a lesser ape (the siamang gibbon (Symphalangus syndactylus)), and untangled the intricacies of their evolution. Compared with the X chromosomes, the ape Y chromosomes vary greatly in size and have low alignability and high levels of structural rearrangements-owing to the accumulation of lineage-specific ampliconic regions, palindromes, transposable elements and satellites. Many Y chromosome genes expand in multi-copy families and some evolve under purifying selection. Thus, the Y chromosome exhibits dynamic evolution, whereas the X chromosome is more stable. Mapping short-read sequencing data to these assemblies revealed diversity and selection patterns on sex chromosomes of more than 100 individual great apes. These reference assemblies are expected to inform human evolution and conservation genetics of non-human apes, all of which are endangered species.

Massively parallel jumping assay decodes Alu retrotransposition activity.

The human genome contains millions of retrotransposons, several of which could become active due to somatic mutations having phenotypic consequences, including disease. However, it is not thoroughly understood how nucleotide changes in retrotransposons affect their jumping activity. Here, we developed a novel massively parallel jumping assay (MPJA) that can test the jumping potential of thousands of transposons en masse. We generated nucleotide variant library of selected four Alu retrotransposons containing 165,087 different haplotypes and tested them for their jumping ability using MPJA. We found 66,821 unique jumping haplotypes, allowing us to pinpoint domains and variants vital for transposition. Mapping these variants to the Alu-RNA secondary structure revealed stem-loop features that contribute to jumping potential. Combined, our work provides a novel high-throughput assay that assesses the ability of retrotransposons to jump and identifies nucleotide changes that have the potential to reactivate them in the human genome.

Structurally divergent and recurrently mutated regions of primate genomes.

We sequenced and assembled using multiple long-read sequencing technologies the genomes of chimpanzee, bonobo, gorilla, orangutan, gibbon, macaque, owl monkey, and marmoset. We identified 1,338,997 lineage-specific fixed structural variants (SVs) disrupting 1,561 protein-coding genes and 136,932 regulatory elements, including the most complete set of human-specific fixed differences. We estimate that 819.47 Mbp or ∼27% of the genome has been affected by SVs across primate evolution. We identify 1,607 structurally divergent regions wherein recurrent structural variation contributes to creating SV hotspots where genes are recurrently lost (e.g., CARD, C4, and OLAH gene families) and additional lineage-specific genes are generated (e.g., CKAP2, VPS36, ACBD7, and NEK5 paralogs), becoming targets of rapid chromosomal diversification and positive selection (e.g., RGPD gene family). High-fidelity long-read sequencing has made these dynamic regions of the genome accessible for sequence-level analyses within and between primate species.

High-throughput methylome analysis reveals differential methylation for early and late onset preeclampsia for mothers and their children.

Preeclampsia is a hypertensive disorder of pregnancy that affects ∼2%-5% of all pregnancies, contributes to 4 of the top 10 causes of pregnancy-related deaths, and remains a long-term risk factor for cardiometabolic diseases. Yet, little is still known about the molecular mechanisms that lead to this disease. There is evidence that some cases have a genetic cause. However, it is well appreciated that harmful factors in the environment, such as poor nutrition, stress, and toxins, may lead to epigenetics changes that can contribute to this disease. DNA methylation is one of the epigenetic modifications known to be fairly stable and impact gene expression. Using DNA from buccal swabs, we analyzed global DNA methylation among three groups of individuals: mothers who experienced 1) early-stage preeclampsia (<32 wk), 2) late-stage preeclampsia (>37 wk), or 3) no complications during their pregnancies, as well as the children from these three groups. We found significant differentially methylated regions (DMRs) between mothers who experienced preeclampsia compared with those with no complications adjacent or within genes that are important for placentation, embryonic development, cell adhesion, and inflammation (e.g., the cadherin pathway). A significant portion of DMR genes showed expression in tissues relevant to preeclampsia (i.e., the brain, heart, kidney, uterus, ovaries, and placenta). As this study was performed on DNA extracted from cheek swabs, this opens the way to future studies in different tissues, aimed at identifying possible biomarkers of risk and early detection, developing targeted interventions, and reducing the progression of this life-threatening disease.NEW & NOTEWORTHY Preeclampsia is a life-threatening hypertensive disorder, affecting 2%-5% of pregnancies, that remains poorly understood. This study analyzed DNA methylation from buccal swabs from mothers who experienced early and late-stage preeclampsia and those with uncomplicated pregnancies, along with their children. Differentially methylated regions were found near and within genes crucial for placental function, embryonic development, and inflammation. Many of these genes are expressed in preeclampsia-related tissues, offering hope for future biomarker development for this condition.

The Complete Sequence and Comparative Analysis of Ape Sex Chromosomes.

Apes possess two sex chromosomes-the male-specific Y and the X shared by males and females. The Y chromosome is crucial for male reproduction, with deletions linked to infertility. The X chromosome carries genes vital for reproduction and cognition. Variation in mating patterns and brain function among great apes suggests corresponding differences in their sex chromosome structure and evolution. However, due to their highly repetitive nature and incomplete reference assemblies, ape sex chromosomes have been challenging to study. Here, using the state-of-the-art experimental and computational methods developed for the telomere-to-telomere (T2T) human genome, we produced gapless, complete assemblies of the X and Y chromosomes for five great apes (chimpanzee, bonobo, gorilla, Bornean and Sumatran orangutans) and a lesser ape, the siamang gibbon. These assemblies completely resolved ampliconic, palindromic, and satellite sequences, including the entire centromeres, allowing us to untangle the intricacies of ape sex chromosome evolution. We found that, compared to the X, ape Y chromosomes vary greatly in size and have low alignability and high levels of structural rearrangements. This divergence on the Y arises from the accumulation of lineage-specific ampliconic regions and palindromes (which are shared more broadly among species on the X) and from the abundance of transposable elements and satellites (which have a lower representation on the X). Our analysis of Y chromosome genes revealed lineage-specific expansions of multi-copy gene families and signatures of purifying selection. In summary, the Y exhibits dynamic evolution, while the X is more stable. Finally, mapping short-read sequencing data from >100 great ape individuals revealed the patterns of diversity and selection on their sex chromosomes, demonstrating the utility of these reference assemblies for studies of great ape evolution. These complete sex chromosome assemblies are expected to further inform conservation genetics of nonhuman apes, all of which are endangered species.

TAD evolutionary and functional characterization reveals diversity in mammalian TAD boundary properties and function.

Topological associating domains (TADs) are self-interacting genomic units crucial for shaping gene regulation patterns. Despite their importance, the extent of their evolutionary conservation and its functional implications remain largely unknown. In this study, we generate Hi-C and ChIP-seq data and compare TAD organization across four primate and four rodent species and characterize the genetic and epigenetic properties of TAD boundaries in correspondence to their evolutionary conservation. We find 14% of all human TAD boundaries to be shared among all eight species (ultraconserved), while 15% are human-specific. Ultraconserved TAD boundaries have stronger insulation strength, CTCF binding, and enrichment of older retrotransposons compared to species-specific boundaries. CRISPR-Cas9 knockouts of an ultraconserved boundary in a mouse model lead to tissue-specific gene expression changes and morphological phenotypes. Deletion of a human-specific boundary near the autism-related AUTS2 gene results in the upregulation of this gene in neurons. Overall, our study provides pertinent TAD boundary evolutionary conservation annotations and showcases the functional importance of TAD evolution.

Grandmaternal allergen sensitization reprograms epigenetic and airway responses to allergen in second-generation offspring.

Asthma susceptibility is influenced by environmental, genetic, and epigenetic factors. DNA methylation is one form of epigenetic modification that regulates gene expression and is both inherited and modified by environmental exposures throughout life. Prenatal development is a particularly vulnerable time period during which exposure to maternal asthma increases asthma risk in offspring. How maternal asthma affects DNA methylation in offspring and what the consequences of differential methylation are in subsequent generations are not fully known. In this study, we tested the effects of grandmaternal house dust mite (HDM) allergen sensitization during pregnancy on airway physiology and inflammation in HDM-sensitized and challenged second-generation mice. We also tested the effects of grandmaternal HDM sensitization on tissue-specific DNA methylation in allergen-naïve and -sensitized second-generation mice. Descendants of both allergen- and vehicle-exposed grandmaternal founders exhibited airway hyperreactivity after HDM sensitization. However, grandmaternal allergen sensitization significantly potentiated airway hyperreactivity and altered the epigenomic trajectory in second-generation offspring after HDM sensitization compared with HDM-sensitized offspring from vehicle-exposed founders. As a result, biological processes and signaling pathways associated with epigenetic modifications were distinct between lineages. A targeted analysis of pathway-associated gene expression found that Smad3 was significantly dysregulated as a result of grandmaternal allergen sensitization. These data show that grandmaternal allergen exposure during pregnancy establishes a unique epigenetic trajectory that reprograms allergen responses in second-generation offspring and may contribute to asthma risk.NEW & NOTEWORTHY Asthma susceptibility is influenced by environmental, genetic, and epigenetic factors. This study shows that maternal allergen exposure during pregnancy promotes unique epigenetic trajectories in second-generation offspring at baseline and in response to allergen sensitization, which is associated with the potentiation of airway hyperreactivity. These effects are one mechanism by which maternal asthma may influence the inheritance of asthma risk.

Structurally divergent and recurrently mutated regions of primate genomes.

To better understand the pattern of primate genome structural variation, we sequenced and assembled using multiple long-read sequencing technologies the genomes of eight nonhuman primate species, including New World monkeys (owl monkey and marmoset), Old World monkey (macaque), Asian apes (orangutan and gibbon), and African ape lineages (gorilla, bonobo, and chimpanzee). Compared to the human genome, we identified 1,338,997 lineage-specific fixed structural variants (SVs) disrupting 1,561 protein-coding genes and 136,932 regulatory elements, including the most complete set of human-specific fixed differences. Across 50 million years of primate evolution, we estimate that 819.47 Mbp or ~27% of the genome has been affected by SVs based on analysis of these primate lineages. We identify 1,607 structurally divergent regions (SDRs) wherein recurrent structural variation contributes to creating SV hotspots where genes are recurrently lost (CARDs, ABCD7, OLAH) and new lineage-specific genes are generated (e.g., CKAP2, NEK5) and have become targets of rapid chromosomal diversification and positive selection (e.g., RGPDs). High-fidelity long-read sequencing has made these dynamic regions of the genome accessible for sequence-level analyses within and between primate species for the first time.

TAD Evolutionary and functional characterization reveals diversity in mammalian TAD boundary properties and function.

Topological associating domains (TADs) are self-interacting genomic units crucial for shaping gene regulation patterns. Despite their importance, the extent of their evolutionary conservation and its functional implications remain largely unknown. In this study, we generate Hi-C and ChIP-seq data and compare TAD organization across four primate and four rodent species, and characterize the genetic and epigenetic properties of TAD boundaries in correspondence to their evolutionary conservation. We find that only 14% of all human TAD boundaries are shared among all eight species (ultraconserved), while 15% are human-specific. Ultraconserved TAD boundaries have stronger insulation strength, CTCF binding, and enrichment of older retrotransposons, compared to species-specific boundaries. CRISPR-Cas9 knockouts of two ultraconserved boundaries in mouse models leads to tissue-specific gene expression changes and morphological phenotypes. Deletion of a human-specific boundary near the autism-related AUTS2 gene results in upregulation of this gene in neurons. Overall, our study provides pertinent TAD boundary evolutionary conservation annotations, and showcase the functional importance of TAD evolution.

Whole-genome sequence and assembly of the Javan gibbon (Hylobates moloch).

The Javan gibbon, Hylobates moloch, is an endangered gibbon species restricted to the forest remnants of western and central Java, Indonesia, and one of the rarest of the Hylobatidae family. Hylobatids consist of 4 genera (Holoock, Hylobates, Symphalangus, and Nomascus) that are characterized by different numbers of chromosomes, ranging from 38 to 52. The underlying cause of this karyotype plasticity is not entirely understood, at least in part, due to the limited availability of genomic data. Here we present the first scaffold-level assembly for H. moloch using a combination of whole-genome Illumina short reads, 10X Chromium linked reads, PacBio, and Oxford Nanopore long reads and proximity-ligation data. This Hylobates genome represents a valuable new resource for comparative genomics studies in primates.

GoPeaks: histone modification peak calling for CUT&Tag.

Genome-wide mapping of histone modifications is critical to understanding transcriptional regulation. CUT&Tag is a new method for profiling histone modifications, offering improved sensitivity and decreased cost compared with ChIP-seq. Here, we present GoPeaks, a peak calling method specifically designed for histone modification CUT&Tag data. We compare the performance of GoPeaks against commonly used peak calling algorithms to detect histone modifications that display a range of peak profiles and are frequently used in epigenetic studies. We find that GoPeaks robustly detects genome-wide histone modifications and, notably, identifies a substantial number of H3K27ac peaks with improved sensitivity compared to other standard algorithms.

DNA Methylation Analysis of Turner Syndrome BAV.

Turner Syndrome (TS) is a rare cytogenetic disorder caused by the complete loss or structural variation of the second sex chromosome. The most common cause of early mortality in TS results from a high incidence of left-sided congenital heart defects, including bicuspid aortic valve (BAV), which occurs in about 30% of individuals with TS. BAV is also the most common congenital heart defect in the general population with a prevalence of 0.5-2%, with males being three-times more likely to have a BAV than females. TS is associated with genome-wide hypomethylation when compared to karyotypically normal males and females. Alterations in DNA methylation in primary aortic tissue are associated with BAV in euploid individuals. Here we show significant differences in DNA methylation patterns associated with BAV in TS found in peripheral blood by comparing TS BAV (n = 12), TS TAV (n = 13), and non-syndromic BAV (n = 6). When comparing TS with BAV to TS with no heart defects we identified a differentially methylated region encompassing the BAV-associated gene MYRF, and enrichment for binding sites of two known transcription factor contributors to BAV. When comparing TS with BAV to euploid women with BAV, we found significant overlapping enrichment for ChIP-seq transcription factor targets including genes in the NOTCH1 pathway, known for involvement in the etiology of non-syndromic BAV, and other genes that are essential regulators of heart valve development. Overall, these findings suggest that altered DNA methylation affecting key aortic valve development genes contributes to the greatly increased risk for BAV in TS.

Somatic retrotransposition in the developing rhesus macaque brain.

The retrotransposon LINE-1 (L1) is central to the recent evolutionary history of the human genome and continues to drive genetic diversity and germline pathogenesis. However, the spatiotemporal extent and biological significance of somatic L1 activity are poorly defined and are virtually unexplored in other primates. From a single L1 lineage active at the divergence of apes and Old World monkeys, successive L1 subfamilies have emerged in each descendant primate germline. As revealed by case studies, the presently active human L1 subfamily can also mobilize during embryonic and brain development in vivo. It is unknown whether nonhuman primate L1s can similarly generate somatic insertions in the brain. Here we applied approximately 40× single-cell whole-genome sequencing (scWGS), as well as retrotransposon capture sequencing (RC-seq), to 20 hippocampal neurons from two rhesus macaques (Macaca mulatta). In one animal, we detected and PCR-validated a somatic L1 insertion that generated target site duplications, carried a short 5' transduction, and was present in ∼7% of hippocampal neurons but absent from cerebellum and nonbrain tissues. The corresponding donor L1 allele was exceptionally mobile in vitro and was embedded in PRDM4, a gene expressed throughout development and in neural stem cells. Nanopore long-read methylome and RNA-seq transcriptome analyses indicated young retrotransposon subfamily activation in the early embryo, followed by repression in adult tissues. These data highlight endogenous macaque L1 retrotransposition potential, provide prototypical evidence of L1-mediated somatic mosaicism in a nonhuman primate, and allude to L1 mobility in the brain over the past 30 million years of human evolution.

Molecular contribution to embryonic aneuploidy and karyotypic complexity in initial cleavage divisions of mammalian development.

Embryonic aneuploidy is highly complex, often leading to developmental arrest, implantation failure or spontaneous miscarriage in both natural and assisted reproduction. Despite our knowledge of mitotic mis-segregation in somatic cells, the molecular pathways regulating chromosome fidelity during the error-prone cleavage-stage of mammalian embryogenesis remain largely undefined. Using bovine embryos and live-cell fluorescent imaging, we observed frequent micro-/multi-nucleation of mis-segregated chromosomes in initial mitotic divisions that underwent unilateral inheritance, re-fused with the primary nucleus or formed a chromatin bridge with neighboring cells. A correlation between a lack of syngamy, multipolar divisions and asymmetric genome partitioning was also revealed, and single-cell DNA-seq showed propagation of primarily non-reciprocal mitotic errors. Depletion of the mitotic checkpoint protein BUB1B (also known as BUBR1) resulted in similarly abnormal nuclear structures and cell divisions, as well as chaotic aneuploidy and dysregulation of the kinase-substrate network that mediates mitotic progression, all before zygotic genome activation. This demonstrates that embryonic micronuclei sustain multiple fates, provides an explanation for blastomeres with uniparental origins, and substantiates defective checkpoints and likely other maternally derived factors as major contributors to the karyotypic complexity afflicting mammalian preimplantation development.

Gene-Specific DNA Methylation Linked to Postoperative Cognitive Dysfunction in Apolipoprotein E3 and E4 Mice.

BACKGROUND: Geriatric surgical patients are at higher risk of developing postoperative neurocognitive disorders (NCD) than younger patients. The specific mechanisms underlying postoperative NCD remain unknown, but they have been linked to genetic risk factors, such as the presence of APOE4, compared to APOE3, and epigenetic modifications caused by exposure to anesthesia and surgery. OBJECTIVE: To test the hypothesis that compared to E3 mice, E4 mice exhibit a more pronounced postoperative cognitive impairment associated with differential DNA methylation in brain regions linked to learning and memory. METHODS: 16-month-old humanized apolipoprotein-E targeted replacement mice bearing E3 or E4 were subjected to surgery (laparotomy) under general isoflurane anesthesia or sham. Postoperative behavioral testing and genome-wide DNA methylation were performed. RESULTS: Exposure to surgery and anesthesia impaired cognition in aged E3, but not E4 mice, likely due to the already lower cognitive performance of E4 prior to surgery. Cognitive impairment in E3 mice was associated with hypermethylation of specific genes, including genes in the Ephrin pathway implicated in synaptic plasticity and learning in adults and has been linked to Alzheimer's disease. Other genes, such as the Scratch Family Transcriptional Repressor 2, were altered after surgery and anesthesia in both the E3 and E4 mice. CONCLUSION: Our findings suggest that the neurocognitive and behavioral effects of surgery and anesthesia depend on baseline neurocognitive status and are associated with APOE isoform-dependent epigenetic modifications of specific genes and pathways involved in memory and learning.