The cycling and aging mouse female reproductive tract at single-cell resolution.

The female reproductive tract (FRT) undergoes extensive remodeling during reproductive cycling. This recurrent remodeling and how it shapes organ-specific aging remains poorly explored. Using single-cell and spatial transcriptomics, we systematically characterized morphological and gene expression changes occurring in ovary, oviduct, uterus, cervix, and vagina at each phase of the mouse estrous cycle, during decidualization, and into aging. These analyses reveal that fibroblasts play central-and highly organ-specific-roles in FRT remodeling by orchestrating extracellular matrix (ECM) reorganization and inflammation. Our results suggest a model wherein recurrent FRT remodeling over reproductive lifespan drives the gradual, age-related development of fibrosis and chronic inflammation. This hypothesis was directly tested using chemical ablation of cycling, which reduced fibrotic accumulation during aging. Our atlas provides extensive detail into how estrus, pregnancy, and aging shape the organs of the female reproductive tract and reveals the unexpected cost of the recurrent remodeling required for reproduction.

Unwinding Limb Development.

The molecular mechanisms underpinning vertebrate body plan evolution are beginning to be unravelled. In this issue of Cell, Kvon et al. spectacularly demonstrate how transplanting snake-specific genetic changes found uniquely in serpent enhancers leads to limb loss in mice.

Enhancer evolution across 20 mammalian species.

The mammalian radiation has corresponded with rapid changes in noncoding regions of the genome, but we lack a comprehensive understanding of regulatory evolution in mammals. Here, we track the evolution of promoters and enhancers active in liver across 20 mammalian species from six diverse orders by profiling genomic enrichment of H3K27 acetylation and H3K4 trimethylation. We report that rapid evolution of enhancers is a universal feature of mammalian genomes. Most of the recently evolved enhancers arise from ancestral DNA exaptation, rather than lineage-specific expansions of repeat elements. In contrast, almost all liver promoters are partially or fully conserved across these species. Our data further reveal that recently evolved enhancers can be associated with genes under positive selection, demonstrating the power of this approach for annotating regulatory adaptations in genomic sequences. These results provide important insight into the functional genetics underpinning mammalian regulatory evolution.

Strand-resolved mutagenicity of DNA damage and repair.

DNA base damage is a major source of oncogenic mutations1. Such damage can produce strand-phased mutation patterns and multiallelic variation through the process of lesion segregation2. Here we exploited these properties to reveal how strand-asymmetric processes, such as replication and transcription, shape DNA damage and repair. Despite distinct mechanisms of leading and lagging strand replication3,4, we observe identical fidelity and damage tolerance for both strands. For small alkylation adducts of DNA, our results support a model in which the same translesion polymerase is recruited on-the-fly to both replication strands, starkly contrasting the strand asymmetric tolerance of bulky UV-induced adducts5. The accumulation of multiple distinct mutations at the site of persistent lesions provides the means to quantify the relative efficiency of repair processes genome wide and at single-base resolution. At multiple scales, we show DNA damage-induced mutations are largely shaped by the influence of DNA accessibility on repair efficiency, rather than gradients of DNA damage. Finally, we reveal specific genomic conditions that can actively drive oncogenic mutagenesis by corrupting the fidelity of nucleotide excision repair. These results provide insight into how strand-asymmetric mechanisms underlie the formation, tolerance and repair of DNA damage, thereby shaping cancer genome evolution.

CTCF and cohesin: linking gene regulatory elements with their targets.

Current epigenomics approaches have facilitated the genome-wide identification of regulatory elements based on chromatin features and transcriptional regulator binding and have begun to map long-range interactions between regulatory elements and their targets. Here, we focus on the emerging roles of CTCF and the cohesin in coordinating long-range interactions between regulatory elements. We discuss how species-specific transposable elements may influence such interactions by remodeling the CTCF binding repertoire and suggest that cohesin's association with enhancers, promoters, and sites defined by CTCF binding has the potential to form developmentally regulated networks of long-range interactions that reflect and promote cell-type-specific transcriptional programs.

Cooperativity and rapid evolution of cobound transcription factors in closely related mammals.

To mechanistically characterize the microevolutionary processes active in altering transcription factor (TF) binding among closely related mammals, we compared the genome-wide binding of three tissue-specific TFs that control liver gene expression in six rodents. Despite an overall fast turnover of TF binding locations between species, we identified thousands of TF regions of highly constrained TF binding intensity. Although individual mutations in bound sequence motifs can influence TF binding, most binding differences occur in the absence of nearby sequence variations. Instead, combinatorial binding was found to be significant for genetic and evolutionary stability; cobound TFs tend to disappear in concert and were sensitive to genetic knockout of partner TFs. The large, qualitative differences in genomic regions bound between closely related mammals, when contrasted with the smaller, quantitative TF binding differences among Drosophila species, illustrate how genome structure and population genetics together shape regulatory evolution.

Waves of retrotransposon expansion remodel genome organization and CTCF binding in multiple mammalian lineages.

CTCF-binding locations represent regulatory sequences that are highly constrained over the course of evolution. To gain insight into how these DNA elements are conserved and spread through the genome, we defined the full spectrum of CTCF-binding sites, including a 33/34-mer motif, and identified over five thousand highly conserved, robust, and tissue-independent CTCF-binding locations by comparing ChIP-seq data from six mammals. Our data indicate that activation of retroelements has produced species-specific expansions of CTCF binding in rodents, dogs, and opossum, which often functionally serve as chromatin and transcriptional insulators. We discovered fossilized repeat elements flanking deeply conserved CTCF-binding regions, indicating that similar retrotransposon expansions occurred hundreds of millions of years ago. Repeat-driven dispersal of CTCF binding is a fundamental, ancient, and still highly active mechanism of genome evolution in mammalian lineages.

DNA lesion bypass and the stochastic dynamics of transcription-coupled repair.

DNA base damage is a major source of oncogenic mutations and disruption to gene expression. The stalling of RNA polymerase II (RNAP) at sites of DNA damage and the subsequent triggering of repair processes have major roles in shaping the genome-wide distribution of mutations, clearing barriers to transcription, and minimizing the production of miscoded gene products. Despite its importance for genetic integrity, key mechanistic features of this transcription-coupled repair (TCR) process are controversial or unknown. Here, we exploited a well-powered in vivo mammalian model system to explore the mechanistic properties and parameters of TCR for alkylation damage at fine spatial resolution and with discrimination of the damaged DNA strand. For rigorous interpretation, a generalizable mathematical model of DNA damage and TCR was developed. Fitting experimental data to the model and simulation revealed that RNA polymerases frequently bypass lesions without triggering repair, indicating that small alkylation adducts are unlikely to be an efficient barrier to gene expression. Following a burst of damage, the efficiency of transcription-coupled repair gradually decays through gene bodies with implications for the occurrence and accurate inference of driver mutations in cancer. The reinitation of transcription from the repair site is not a general feature of transcription-coupled repair, and the observed data is consistent with reinitiation never taking place. Collectively, these results reveal how the directional but stochastic activity of TCR shapes the distribution of mutations following DNA damage.

Pervasive lesion segregation shapes cancer genome evolution.

Cancers arise through the acquisition of oncogenic mutations and grow by clonal expansion1,2. Here we reveal that most mutagenic DNA lesions are not resolved into a mutated DNA base pair within a single cell cycle. Instead, DNA lesions segregate, unrepaired, into daughter cells for multiple cell generations, resulting in the chromosome-scale phasing of subsequent mutations. We characterize this process in mutagen-induced mouse liver tumours and show that DNA replication across persisting lesions can produce multiple alternative alleles in successive cell divisions, thereby generating both multiallelic and combinatorial genetic diversity. The phasing of lesions enables accurate measurement of strand-biased repair processes, quantification of oncogenic selection and fine mapping of sister-chromatid-exchange events. Finally, we demonstrate that lesion segregation is a unifying property of exogenous mutagens, including UV light and chemotherapy agents in human cells and tumours, which has profound implications for the evolution and adaptation of cancer genomes.

Polyploidisation pleiotropically buffers ageing in hepatocytes.

BACKGROUND & AIMS: Polyploidy in hepatocytes has been proposed as a genetic mechanism to buffer against transcriptional dysregulation. Here, we aim to demonstrate the role of polyploidy in modulating gene regulatory networks in hepatocytes during ageing. METHODS: We performed single-nucleus RNA sequencing in hepatocyte nuclei of different ploidy levels isolated from young and old wild-type mice. Changes in the gene expression and regulatory network were compared to three independent strains that were haploinsufficient for HNF4A, CEBPA or CTCF, representing non-deleterious perturbations. Phenotypic characteristics of the liver section were additionally evaluated histologically, whereas the genomic allele composition of hepatocytes was analysed by BaseScope. RESULTS: We observed that ageing in wild-type mice results in nuclei polyploidy and a marked increase in steatosis. Haploinsufficiency of liver-specific master regulators (HFN4A or CEBPA) results in the enrichment of hepatocytes with tetraploid nuclei at a young age, affecting the genomic regulatory network, and dramatically suppressing ageing-related steatosis tissue wide. Notably, these phenotypes are not the result of subtle disruption to liver-specific transcriptional networks, since haploinsufficiency in the CTCF insulator protein resulted in the same phenotype. Further quantification of genotypes of tetraploid hepatocytes in young and old HFN4A-haploinsufficient mice revealed that during ageing, tetraploid hepatocytes lead to the selection of wild-type alleles, restoring non-deleterious genetic perturbations. CONCLUSIONS: Our results suggest a model whereby polyploidisation leads to fundamentally different cell states. Polyploid conversion enables pleiotropic buffering against age-related decline via non-random allelic segregation to restore a wild-type genome. IMPACT AND IMPLICATIONS: The functional role of hepatocyte polyploidisation during ageing is poorly understood. Using single-nucleus RNA sequencing and BaseScope approaches, we have studied ploidy dynamics during ageing in murine livers with non-deleterious genetic perturbations. We have identified that hepatocytes present different cellular states and the ability to buffer ageing-associated dysfunctions. Tetraploid nuclei exhibit robust transcriptional networks and are better adapted to genomically overcome perturbations. Novel therapeutic interventions aimed at attenuating age-related changes in tissue function could be exploited by manipulation of ploidy dynamics during chronic liver conditions.

Chromosome assembly of large and complex genomes using multiple references.

Despite the rapid development of sequencing technologies, the assembly of mammalian-scale genomes into complete chromosomes remains one of the most challenging problems in bioinformatics. To help address this difficulty, we developed Ragout 2, a reference-assisted assembly tool that works for large and complex genomes. By taking one or more target assemblies (generated from an NGS assembler) and one or multiple related reference genomes, Ragout 2 infers the evolutionary relationships between the genomes and builds the final assemblies using a genome rearrangement approach. By using Ragout 2, we transformed NGS assemblies of 16 laboratory mouse strains into sets of complete chromosomes, leaving <5% of sequence unlocalized per set. Various benchmarks, including PCR testing and realigning of long Pacific Biosciences (PacBio) reads, suggest only a small number of structural errors in the final assemblies, comparable with direct assembly approaches. We applied Ragout 2 to the Mus caroli and Mus pahari genomes, which exhibit karyotype-scale variations compared with other genomes from the Muridae family. Chromosome painting maps confirmed most large-scale rearrangements that Ragout 2 detected. We applied Ragout 2 to improve draft sequences of three ape genomes that have recently been published. Ragout 2 transformed three sets of contigs (generated using PacBio reads only) into chromosome-scale assemblies with accuracy comparable to chromosome assemblies generated in the original study using BioNano maps, Hi-C, BAC clones, and FISH.

Repeat associated mechanisms of genome evolution and function revealed by the Mus caroli and Mus pahari genomes.

Understanding the mechanisms driving lineage-specific evolution in both primates and rodents has been hindered by the lack of sister clades with a similar phylogenetic structure having high-quality genome assemblies. Here, we have created chromosome-level assemblies of the Mus caroli and Mus pahari genomes. Together with the Mus musculus and Rattus norvegicus genomes, this set of rodent genomes is similar in divergence times to the Hominidae (human-chimpanzee-gorilla-orangutan). By comparing the evolutionary dynamics between the Muridae and Hominidae, we identified punctate events of chromosome reshuffling that shaped the ancestral karyotype of Mus musculus and Mus caroli between 3 and 6 million yr ago, but that are absent in the Hominidae. Hominidae show between four- and sevenfold lower rates of nucleotide change and feature turnover in both neutral and functional sequences, suggesting an underlying coherence to the Muridae acceleration. Our system of matched, high-quality genome assemblies revealed how specific classes of repeats can play lineage-specific roles in related species. Recent LINE activity has remodeled protein-coding loci to a greater extent across the Muridae than the Hominidae, with functional consequences at the species level such as reproductive isolation. Furthermore, we charted a Muridae-specific retrotransposon expansion at unprecedented resolution, revealing how a single nucleotide mutation transformed a specific SINE element into an active CTCF binding site carrier specifically in Mus caroli, which resulted in thousands of novel, species-specific CTCF binding sites. Our results show that the comparison of matched phylogenetic sets of genomes will be an increasingly powerful strategy for understanding mammalian biology.

IL-7-dependent compositional changes within the γδ T cell pool in lymph nodes during ageing lead to an unbalanced anti-tumour response.

How the age-associated decline of immune function leads to increased cancer incidence is poorly understood. Here, we have characterised the cellular composition of the γδ T-cell pool in peripheral lymph nodes (pLNs) upon ageing. We find that ageing has minimal cell-intrinsic effects on function and global gene expression of γδ T cells, and γδTCR diversity remains stable. However, ageing alters TCRδ chain usage and clonal structure of γδ T-cell subsets. Importantly, IL-17-producing γδ17 T cells dominate the γδ T-cell pool of aged mice-mainly due to the selective expansion of Vγ6+ γδ17 T cells and augmented γδ17 polarisation of Vγ4+ T cells. Expansion of the γδ17 T-cell compartment is mediated by increased IL-7 expression in the T-cell zone of old mice. In a Lewis lung cancer model, pro-tumourigenic Vγ6+ γδ17 T cells are exclusively activated in the tumour-draining LN and their infiltration into the tumour correlates with increased tumour size in aged mice. Thus, upon ageing, substantial compositional changes in γδ T-cell pool in the pLN lead to an unbalanced γδ T-cell response in the tumour that is associated with accelerated tumour growth.

Latent regulatory potential of human-specific repetitive elements.

At least half of the human genome is derived from repetitive elements, which are often lineage specific and silenced by a variety of genetic and epigenetic mechanisms. Using a transchromosomic mouse strain that transmits an almost complete single copy of human chromosome 21 via the female germline, we show that a heterologous regulatory environment can transcriptionally activate transposon-derived human regulatory regions. In the mouse nucleus, hundreds of locations on human chromosome 21 newly associate with activating histone modifications in both somatic and germline tissues, and influence the gene expression of nearby transcripts. These regions are enriched with primate and human lineage-specific transposable elements, and their activation corresponds to changes in DNA methylation at CpG dinucleotides. This study reveals the latent regulatory potential of the repetitive human genome and illustrates the species specificity of mechanisms that control it.

Functional signatures of evolutionarily young CTCF binding sites.

BACKGROUND: The introduction of novel CTCF binding sites in gene regulatory regions in the rodent lineage is partly the effect of transposable element expansion, particularly in the murine lineage. The exact mechanism and functional impact of evolutionarily novel CTCF binding sites are not yet fully understood. We investigated the impact of novel subspecies-specific CTCF binding sites in two Mus genus subspecies, Mus musculus domesticus and Mus musculus castaneus, that diverged 0.5 million years ago. RESULTS: CTCF binding site evolution is influenced by the action of the B2-B4 family of transposable elements independently in both lineages, leading to the proliferation of novel CTCF binding sites. A subset of evolutionarily young sites may harbour transcriptional functionality as evidenced by the stability of their binding across multiple tissues in M. musculus domesticus (BL6), while overall the distance of subspecies-specific CTCF binding to the nearest transcription start sites and/or topologically associated domains (TADs) is largely similar to musculus-common CTCF sites. Remarkably, we discovered a recurrent regulatory architecture consisting of a CTCF binding site and an interferon gene that appears to have been tandemly duplicated to create a 15-gene cluster on chromosome 4, thus forming a novel BL6 specific immune locus in which CTCF may play a regulatory role. CONCLUSIONS: Our results demonstrate that thousands of CTCF binding sites show multiple functional signatures rapidly after incorporation into the genome.

Crisflash: open-source software to generate CRISPR guide RNAs against genomes annotated with individual variation.

SUMMARY: CRISPR/Cas9 system requires short guide RNAs (sgRNAs) to direct genome modification. Most currently available tools for sgRNA design operate only with standard reference genomes, and are best suited for small-scale projects. To address these limitations, we developed Crisflash, a software tool for fast sgRNA design and potential off-target discovery, built for performance and flexibility. Crisflash can rapidly design CRISPR guides against any sequenced genome or genome sequences, and can optimize guide accuracy by incorporating user-supplied variant data. Crisflash is over an order of magnitude faster than comparable tools, even using a single CPU core, and efficiently and robustly scores the potential off-targeting of all possible candidate CRISPR guide oligonucleotides. AVAILABILITY AND IMPLEMENTATION: https://github.com/crisflash. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Evolution of transcription factor binding in metazoans - mechanisms and functional implications.

Differences in transcription factor binding can contribute to organismal evolution by altering downstream gene expression programmes. Genome-wide studies in Drosophila melanogaster and mammals have revealed common quantitative and combinatorial properties of in vivo DNA binding, as well as marked differences in the rate and mechanisms of evolution of transcription factor binding in metazoans. Here, we review the recently discovered rapid 're-wiring' of in vivo transcription factor binding between related metazoan species and summarize general principles underlying the observed patterns of evolution. We then consider what might explain the differences in genome evolution between metazoan phyla and outline the conceptual and technological challenges facing this research field.

Independence of repressive histone marks and chromatin compaction during senescent heterochromatic layer formation.

The expansion of repressive epigenetic marks has been implicated in heterochromatin formation during embryonic development, but the general applicability of this mechanism is unclear. Here we show that nuclear rearrangement of repressive histone marks H3K9me3 and H3K27me3 into nonoverlapping structural layers characterizes senescence-associated heterochromatic foci (SAHF) formation in human fibroblasts. However, the global landscape of these repressive marks remains unchanged upon SAHF formation, suggesting that in somatic cells, heterochromatin can be formed through the spatial repositioning of pre-existing repressively marked histones. This model is reinforced by the correlation of presenescent replication timing with both the subsequent layered structure of SAHFs and the global landscape of the repressive marks, allowing us to integrate microscopic and genomic information. Furthermore, modulation of SAHF structure does not affect the occupancy of these repressive marks, nor vice versa. These experiments reveal that high-order heterochromatin formation and epigenetic remodeling of the genome can be discrete events.

Specificity of RNAi, LNA and CRISPRi as loss-of-function methods in transcriptional analysis.

Loss-of-function (LOF) methods such as RNA interference (RNAi), antisense oligonucleotides or CRISPR-based genome editing provide unparalleled power for studying the biological function of genes of interest. However, a major concern is non-specific targeting, which involves depletion of transcripts other than those intended. Little work has been performed to characterize the off-target effects of these common LOF methods at the whole-transcriptome level. Here, we experimentally compared the non-specific activity of RNAi, antisense oligonucleotides and CRISPR interference (CRISPRi). All three methods yielded non-negligible off-target effects in gene expression, with CRISPRi also exhibiting strong clonal effects. As an illustrative example, we evaluated the performance of each method for determining the role of an uncharacterized long noncoding RNA (lncRNA). Several LOF methods successfully depleted the candidate lncRNA but yielded different sets of differentially expressed genes as well as a different cellular phenotype upon depletion. Similar discrepancies between methods were observed with a protein-coding gene (Ch-TOG/CKAP5) and another lncRNA (MALAT1). We suggest that the differences between methods arise due to method-specific off-target effects and provide guidelines for mitigating such effects in functional studies. Our recommendations provide a framework with which off-target effects can be managed to improve functional characterization of genes of interest.

Codon-Driven Translational Efficiency Is Stable across Diverse Mammalian Cell States.

Whether codon usage fine-tunes mRNA translation in mammals remains controversial, with recent papers suggesting that production of proteins in specific Gene Ontological (GO) pathways can be regulated by actively modifying the codon and anticodon pools in different cellular conditions. In this work, we compared the sequence content of genes in specific GO categories with the exonic genome background. Although a substantial fraction of variability in codon usage could be explained by random sampling, almost half of GO sets showed more variability in codon usage than expected by chance. Nevertheless, by quantifying translational efficiency in healthy and cancerous tissues in human and mouse, we demonstrated that a given tRNA pool can equally well translate many different sets of mRNAs, irrespective of their cell-type specificity. This disconnect between variations in codon usage and the stability of translational efficiency is best explained by differences in GC content between gene sets. GC variation across the mammalian genome is most likely a result of the interplay between genome repair and gene duplication mechanisms, rather than selective pressures caused by codon-driven translational rates. Consequently, codon usage differences in mammalian transcriptomes are most easily explained by well-understood mutational biases acting on the underlying genome.