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.

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.

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.

De novo assembly and annotation of the singing mouse genome.

BACKGROUND: Developing genomic resources for a diverse range of species is an important step towards understanding the mechanisms underlying complex traits. Specifically, organisms that exhibit unique and accessible phenotypes-of-interest allow researchers to address questions that may be ill-suited to traditional model organisms. We sequenced the genome and transcriptome of Alston's singing mouse (Scotinomys teguina), an emerging model for social cognition and vocal communication. In addition to producing advertisement songs used for mate attraction and male-male competition, these rodents are diurnal, live at high-altitudes, and are obligate insectivores, providing opportunities to explore diverse physiological, ecological, and evolutionary questions. RESULTS: Using PromethION, Illumina, and PacBio sequencing, we produced an annotated genome and transcriptome, which were validated using gene expression and functional enrichment analyses. To assess the usefulness of our assemblies, we performed single nuclei sequencing on cells of the orofacial motor cortex, a brain region implicated in song coordination, identifying 12 cell types. CONCLUSIONS: These resources will provide the opportunity to identify the molecular basis of complex traits in singing mice as well as to contribute data that can be used for large-scale comparative analyses.

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.

Comparative Analyses of Gibbon Centromeres Reveal Dynamic Genus-Specific Shifts in Repeat Composition.

Centromeres are functionally conserved chromosomal loci essential for proper chromosome segregation during cell division, yet they show high sequence diversity across species. Despite their variation, a near universal feature of centromeres is the presence of repetitive sequences, such as DNA satellites and transposable elements (TEs). Because of their rapidly evolving karyotypes, gibbons represent a compelling model to investigate divergence of functional centromere sequences across short evolutionary timescales. In this study, we use ChIP-seq, RNA-seq, and fluorescence in situ hybridization to comprehensively investigate the centromeric repeat content of the four extant gibbon genera (Hoolock, Hylobates, Nomascus, and Siamang). In all gibbon genera, we find that CENP-A nucleosomes and the DNA-proteins that interface with the inner kinetochore preferentially bind retroelements of broad classes rather than satellite DNA. A previously identified gibbon-specific composite retrotransposon, LAVA, known to be expanded within the centromere regions of one gibbon genus (Hoolock), displays centromere- and species-specific sequence differences, potentially as a result of its co-option to a centromeric function. When dissecting centromere satellite composition, we discovered the presence of the retroelement-derived macrosatellite SST1 in multiple centromeres of Hoolock, whereas alpha-satellites represent the predominate satellite in the other genera, further suggesting an independent evolutionary trajectory for Hoolock centromeres. Finally, using de novo assembly of centromere sequences, we determined that transcripts originating from gibbon centromeres recapitulate the species-specific TE composition. Combined, our data reveal dynamic shifts in the repeat content that define gibbon centromeres and coincide with the extensive karyotypic diversity within this lineage.

Integrated analysis of behavioral, epigenetic, and gut microbiome analyses in AppNL-G-F, AppNL-F, and wild type mice.

Epigenetic mechanisms occurring in the brain as well as alterations in the gut microbiome composition might contribute to Alzheimer's disease (AD). Human amyloid precursor protein knock-in (KI) mice contain the Swedish and Iberian mutations (AppNL-F) or those two and also the Arctic mutation (AppNL-G-F). In this study, we assessed whether behavioral and cognitive performance in 6-month-old AppNL-F, AppNL-G-F, and C57BL/6J wild-type (WT) mice was associated with the gut microbiome, and whether the genotype modulates this association. The genotype effects observed in behavioral tests were test-dependent. The biodiversity and composition of the gut microbiome linked to various aspects of mouse behavioral and cognitive performance but differences in genotype modulated these relationships. These genotype-dependent associations include members of the Lachnospiraceae and Ruminococcaceae families. In a subset of female mice, we assessed DNA methylation in the hippocampus and investigated whether alterations in hippocampal DNA methylation were associated with the gut microbiome. Among other differentially methylated regions, we identified a 1 Kb region that overlapped ing 3'UTR of the Tomm40 gene and the promoter region of the Apoe gene that and was significantly more methylated in the hippocampus of AppNL-G-F than WT mice. The integrated gut microbiome hippocampal DNA methylation analysis revealed a positive relationship between amplicon sequence variants (ASVs) within the Lachnospiraceae family and methylation at the Apoe gene. Hence, these microbes may elicit an impact on AD-relevant behavioral and cognitive performance via epigenetic changes in AD-susceptibility genes in neural tissue or that such changes in the epigenome can elicit alterations in intestinal physiology that affect the growth of these taxa in the gut microbiome.

Dynamic changes in chromatin accessibility, altered adipogenic gene expression, and total versus de novo fatty acid synthesis in subcutaneous adipose stem cells of normal-weight polycystic ovary syndrome (PCOS) women during adipogenesis: evidence of cellular programming.

BACKGROUND: Normal-weight polycystic ovary syndrome (PCOS) women exhibit adipose resistance in vivo accompanied by enhanced subcutaneous (SC) abdominal adipose stem cell (ASC) development to adipocytes with accelerated lipid accumulation per cell in vitro. The present study examines chromatin accessibility, RNA expression and fatty acid (FA) synthesis during SC abdominal ASC differentiation into adipocytes in vitro of normal-weight PCOS versus age- and body mass index-matched normoandrogenic ovulatory (control) women to study epigenetic/genetic characteristics as well as functional alterations of PCOS and control ASCs during adipogenesis. RESULTS: SC abdominal ASCs from PCOS women versus controls exhibited dynamic chromatin accessibility during adipogenesis, from significantly less chromatin accessibility at day 0 to greater chromatin accessibility by day 12, with enrichment of binding motifs for transcription factors (TFs) of the AP-1 subfamily at days 0, 3, and 12. In PCOS versus control cells, expression of genes governing adipocyte differentiation (PPARγ, CEBPα, AGPAT2) and function (ADIPOQ, FABP4, LPL, PLIN1, SLC2A4) was increased two-sixfold at days 3, 7, and 12, while that involving Wnt signaling (FZD1, SFRP1, and WNT10B) was decreased. Differential gene expression in PCOS cells at these time points involved triacylglycerol synthesis, lipid oxidation, free fatty acid beta-oxidation, and oxidative phosphorylation of the TCA cycle, with TGFB1 as a significant upstream regulator. There was a broad correspondence between increased chromatin accessibility and increased RNA expression of those 12 genes involved in adipocyte differentiation and function, Wnt signaling, as well as genes involved in the triacylglycerol synthesis functional group at day 12 of adipogenesis. Total content and de novo synthesis of myristic (C14:0), palmitic (C16:0), palmitoleic (C16:1), and oleic (C18:1) acid increased from day 7 to day 12 in all cells, with total content and de novo synthesis of FAs significantly greater in PCOS than controls cells at day 12. CONCLUSIONS: In normal-weight PCOS women, dynamic chromatin remodeling of SC abdominal ASCs during adipogenesis may enhance adipogenic gene expression as a programmed mechanism to promote greater fat storage.

Co-option of the lineage-specific LAVA retrotransposon in the gibbon genome.

Co-option of transposable elements (TEs) to become part of existing or new enhancers is an important mechanism for evolution of gene regulation. However, contributions of lineage-specific TE insertions to recent regulatory adaptations remain poorly understood. Gibbons present a suitable model to study these contributions as they have evolved a lineage-specific TE called LAVA (LINE-AluSz-VNTR-AluLIKE), which is still active in the gibbon genome. The LAVA retrotransposon is thought to have played a role in the emergence of the highly rearranged structure of the gibbon genome by disrupting transcription of cell cycle genes. In this study, we investigated whether LAVA may have also contributed to the evolution of gene regulation by adopting enhancer function. We characterized fixed and polymorphic LAVA insertions across multiple gibbons and found 96 LAVA elements overlapping enhancer chromatin states. Moreover, LAVA was enriched in multiple transcription factor binding motifs, was bound by an important transcription factor (PU.1), and was associated with higher levels of gene expression in cis We found gibbon-specific signatures of purifying/positive selection at 27 LAVA insertions. Two of these insertions were fixed in the gibbon lineage and overlapped with enhancer chromatin states, representing putative co-opted LAVA enhancers. These putative enhancers were located within genes encoding SETD2 and RAD9A, two proteins that facilitate accurate repair of DNA double-strand breaks and prevent chromosomal rearrangement mutations. Co-option of LAVA in these genes may have influenced regulation of processes that preserve genome integrity. Our findings highlight the importance of considering lineage-specific TEs in studying evolution of gene regulatory elements.

Myeloid lineage enhancers drive oncogene synergy in CEBPA/CSF3R mutant acute myeloid leukemia.

Acute Myeloid Leukemia (AML) develops due to the acquisition of mutations from multiple functional classes. Here, we demonstrate that activating mutations in the granulocyte colony stimulating factor receptor (CSF3R), cooperate with loss of function mutations in the transcription factor CEBPA to promote acute leukemia development. The interaction between these distinct classes of mutations occurs at the level of myeloid lineage enhancers where mutant CEBPA prevents activation of a subset of differentiation associated enhancers. To confirm this enhancer-dependent mechanism, we demonstrate that CEBPA mutations must occur as the initial event in AML initiation. This improved mechanistic understanding will facilitate therapeutic development targeting the intersection of oncogene cooperativity.

Epigenetic maintenance of topological domains in the highly rearranged gibbon genome.

The relationship between evolutionary genome remodeling and the three-dimensional structure of the genome remain largely unexplored. Here, we use the heavily rearranged gibbon genome to examine how evolutionary chromosomal rearrangements impact genome-wide chromatin interactions, topologically associating domains (TADs), and their epigenetic landscape. We use high-resolution maps of gibbon-human breaks of synteny (BOS), apply Hi-C in gibbon, measure an array of epigenetic features, and perform cross-species comparisons. We find that gibbon rearrangements occur at TAD boundaries, independent of the parameters used to identify TADs. This overlap is supported by a remarkable genetic and epigenetic similarity between BOS and TAD boundaries, namely presence of CpG islands and SINE elements, and enrichment in CTCF and H3K4me3 binding. Cross-species comparisons reveal that regions orthologous to BOS also correspond with boundaries of large (400-600 kb) TADs in human and other mammalian species. The colocalization of rearrangement breakpoints and TAD boundaries may be due to higher chromatin fragility at these locations and/or increased selective pressure against rearrangements that disrupt TAD integrity. We also examine the small portion of BOS that did not overlap with TAD boundaries and gave rise to novel TADs in the gibbon genome. We postulate that these new TADs generally lack deleterious consequences. Last, we show that limited epigenetic homogenization occurs across breakpoints, irrespective of their time of occurrence in the gibbon lineage. Overall, our findings demonstrate remarkable conservation of chromatin interactions and epigenetic landscape in gibbons, in spite of extensive genomic shuffling.

Complex selection on a regulator of social cognition: Evidence of balancing selection, regulatory interactions and population differentiation in the prairie vole Avpr1a locus.

Adaptive variation in social behaviour depends upon standing genetic variation, but we know little about how evolutionary forces shape genetic diversity relevant to brain and behaviour. In prairie voles (Microtus ochrogaster), variants at the Avpr1a locus predict expression of the vasopressin 1a receptor in the retrosplenial cortex (RSC), a brain region that mediates spatial and contextual memory; cortical V1aR abundance in turn predicts diversity in space use and sexual fidelity in the field. To examine the potential contributions of adaptive and neutral forces to variation at the Avpr1a locus, we explore sequence diversity at the Avpr1a locus and throughout the genome in two populations of wild prairie voles. First, we refine results demonstrating balancing selection at the locus by comparing the frequency spectrum of variants at the locus to a random sample of the genome. Next, we find that the four single nucleotide polymorphisms that predict high V1aR expression in the RSC are in stronger linkage disequilibrium than expected by chance despite high recombination among intervening variants, suggesting that epistatic selection maintains their association. Analysis of population structure and a haplotype network for two populations revealed that this excessive LD was unlikely to be due to admixture alone. Furthermore, the two populations differed considerably in the region shown to be a regulator of V1aR expression despite the extremely low levels of genomewide genetic differentiation. Together, our data suggest that complex selection on Avpr1a locus favours specific combinations of regulatory polymorphisms, maintains the resulting alleles at population-specific frequencies, and may contribute to unique patterns of spatial cognition and sexual fidelity among populations.

Individual Differences in Social Behavior and Cortical Vasopressin Receptor: Genetics, Epigenetics, and Evolution.

Social behavior is among the most complex and variable of traits. Despite its diversity, we know little about how genetic and developmental factors interact to shape natural variation in social behavior. This review surveys recent work on individual differences in the expression of the vasopressin 1a receptor (V1aR), a major regulator of social behavior, in the neocortex of the socially monogamous prairie vole. V1aR exhibits profound variation in the retrosplenial cortex (RSC), a region critical to spatial and contextual memory. RSC-V1aR abundance is associated with patterns of male space-use and sexual fidelity in the field: males with high RSC-V1aR show high spatial and sexual fidelity to partners, while low RSC-V1aR males are significantly more likely to mate outside the pair-bond. Individual differences in RSC-V1aR are predicted by a set of linked single nucleotide polymorphisms within the avpr1a locus. These alternative alleles have been actively maintained by selection, suggesting that the brain differences represent a balanced polymorphism. Lastly, the alleles occur within regulatory sequences, and result in differential sensitivity to environmental perturbation. Together the data provide insight into how genetic, epigenetic and evolutionary forces interact to shape the social brain.

Sexual fidelity trade-offs promote regulatory variation in the prairie vole brain.

Individual variation in social behavior seems ubiquitous, but we know little about how it relates to brain diversity. Among monogamous prairie voles, levels of vasopressin receptor (encoded by the gene avpr1a) in brain regions related to spatial memory predict male space use and sexual fidelity in the field. We find that trade-offs between the benefits of male fidelity and infidelity are reflected in patterns of territorial intrusion, offspring paternity, avpr1a expression, and the evolutionary fitness of alternative avpr1a alleles. DNA variation at the avpr1a locus includes polymorphisms that reliably predict the epigenetic status and neural expression of avpr1a, and patterns of DNA diversity demonstrate that avpr1a regulatory variation has been favored by selection. In prairie voles, trade-offs in the fitness consequences of social behaviors seem to promote neuronal and molecular diversity.