Mutational bias in spermatogonia impacts the anatomy of regulatory sites in the human genome.

Mutation in the germline is the ultimate source of genetic variation, but little is known about the influence of germline chromatin structure on mutational processes. Using ATAC-seq, we profile the open chromatin landscape of human spermatogonia, the most proliferative cell type of the germline, identifying transcription factor binding sites (TFBSs) and PRDM9 binding sites, a subset of which will initiate meiotic recombination. We observe an increase in rare structural variant (SV) breakpoints at PRDM9-bound sites, implicating meiotic recombination in the generation of structural variation. Many germline TFBSs, such as NRF1, are also associated with increased rates of SV breakpoints, apparently independent of recombination. Singleton short insertions (≥5 bp) are highly enriched at TFBSs, particularly at sites bound by testis active TFs, and their rates correlate with those of structural variant breakpoints. Short insertions often duplicate the TFBS motif, leading to clustering of motif sites near regulatory regions in this male-driven evolutionary process. Increased mutation loads at germline TFBSs disproportionately affect neural enhancers with activity in spermatogonia, potentially altering neurodevelopmental regulatory architecture. Local chromatin structure in spermatogonia is thus pervasive in shaping both evolution and disease.

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.

Chromatin loop anchors are associated with genome instability in cancer and recombination hotspots in the germline.

BACKGROUND: Chromatin loops form a basic unit of interphase nuclear organization, with chromatin loop anchor points providing contacts between regulatory regions and promoters. However, the mutational landscape at these anchor points remains under-studied. Here, we describe the unusual patterns of somatic mutations and germline variation associated with loop anchor points and explore the underlying features influencing these patterns. RESULTS: Analyses of whole genome sequencing datasets reveal that anchor points are strongly depleted for single nucleotide variants (SNVs) in tumours. Despite low SNV rates in their genomic neighbourhood, anchor points emerge as sites of evolutionary innovation, showing enrichment for structural variant (SV) breakpoints and a peak of SNVs at focal CTCF sites within the anchor points. Both CTCF-bound and non-CTCF anchor points harbour an excess of SV breakpoints in multiple tumour types and are prone to double-strand breaks in cell lines. Common fragile sites, which are hotspots for genome instability, also show elevated numbers of intersecting loop anchor points. Recurrently disrupted anchor points are enriched for genes with functions in cell cycle transitions and regions associated with predisposition to cancer. We also discover a novel class of CTCF-bound anchor points which overlap meiotic recombination hotspots and are enriched for the core PRDM9 binding motif, suggesting that the anchor points have been foci for diversity generated during recent human evolution. CONCLUSIONS: We suggest that the unusual chromatin environment at loop anchor points underlies the elevated rates of variation observed, marking them as sites of regulatory importance but also genomic fragility.

When TADs go bad: chromatin structure and nuclear organisation in human disease.

Chromatin in the interphase nucleus is organised as a hierarchical series of structural domains, including self-interacting domains called topologically associating domains (TADs). This arrangement is thought to bring enhancers into closer physical proximity with their target genes, which often are located hundreds of kilobases away in linear genomic distance. TADs are demarcated by boundary regions bound by architectural proteins, such as CTCF and cohesin, although much remains to be discovered about the structure and function of these domains. Recent studies of TAD boundaries disrupted in engineered mouse models show that boundary mutations can recapitulate human developmental disorders as a result of aberrant promoter-enhancer interactions in the affected TADs. Similar boundary disruptions in certain cancers can result in oncogene overexpression, and CTCF binding sites at boundaries appear to be hyper-mutated across cancers. Further insights into chromatin organisation, in parallel with accumulating whole genome sequence data for disease cohorts, are likely to yield additional valuable insights into the roles of noncoding sequence variation in human disease.

Mutational Biases Drive Elevated Rates of Substitution at Regulatory Sites across Cancer Types.

Disruption of gene regulation is known to play major roles in carcinogenesis and tumour progression. Here, we comprehensively characterize the mutational profiles of diverse transcription factor binding sites (TFBSs) across 1,574 completely sequenced cancer genomes encompassing 11 tumour types. We assess the relative rates and impact of the mutational burden at the binding sites of 81 transcription factors (TFs), by comparing the abundance and patterns of single base substitutions within putatively functional binding sites to control sites with matched sequence composition. There is a strong (1.43-fold) and significant excess of mutations at functional binding sites across TFs, and the mutations that accumulate in cancers are typically more disruptive than variants tolerated in extant human populations at the same sites. CTCF binding sites suffer an exceptionally high mutational load in cancer (3.31-fold excess) relative to control sites, and we demonstrate for the first time that this effect is seen in essentially all cancer types with sufficient data. The sub-set of CTCF sites involved in higher order chromatin structures has the highest mutational burden, suggesting a widespread breakdown of chromatin organization. However, we find no evidence for selection driving these distinctive patterns of mutation. The mutational load at CTCF-binding sites is substantially determined by replication timing and the mutational signature of the tumor in question, suggesting that selectively neutral processes underlie the unusual mutation patterns. Pervasive hyper-mutation within transcription factor binding sites rewires the regulatory landscape of the cancer genome, but it is dominated by mutational processes rather than selection.

Homozygous loss-of-function variants in European cosmopolitan and isolate populations.

Homozygous loss of function (HLOF) variants provide a valuable window on gene function in humans, as well as an inventory of the human genes that are not essential for survival and reproduction. All humans carry at least a few HLOF variants, but the exact number of inactivated genes that can be tolerated is currently unknown­as are the phenotypic effects of losing function for most human genes. Here, we make use of 1432 whole exome sequences from five European populations to expand the catalogue of known human HLOF mutations; after stringent filtering of variants in our dataset, we identify a total of 173 HLOF mutations, 76 (44%) of which have not been observed previously. We find that population isolates are particularly well suited to surveys of novel HLOF genes because individuals in such populations carry extensive runs of homozygosity, which we show are enriched for novel, rare HLOF variants. Further, we make use of extensive phenotypic data to show that most HLOFs, ascertained in population-based samples, appear to have little detectable effect on the phenotype. On the contrary, we document several genes directly implicated in disease that seem to tolerate HLOF variants. Overall HLOF genes are enriched for olfactory receptor function and are expressed in testes more often than expected, consistent with reduced purifying selection and incipient pseudogenisation.

De novo transcriptome assembly reveals sex-specific selection acting on evolving neo-sex chromosomes in Drosophila miranda.

BACKGROUND: The Drosophila miranda neo-sex chromosome system is a useful resource for studying recently evolved sex chromosomes. However, the neo-Y genomic assembly is fragmented due to the accumulation of repetitive sequence. Furthermore, the separate assembly of the neo-X and neo-Y chromosomes into genomic scaffolds has proven to be difficult, due to their low level of sequence divergence, which in coding regions is about 1.5%. Here, we de novo assemble the transcriptome of D. miranda using RNA-seq data from several male and female tissues, and develop a bioinformatic pipeline to separately reconstruct neo-X and neo-Y transcripts. RESULTS: We obtain 2,141 transcripts from the neo-X and 1,863 from the neo-Y. Neo-Y transcripts are generally shorter than their homologous neo-X transcripts (N50 of 2,048-bp vs. 2,775-bp) and expressed at lower levels. We find that 24% of expressed neo-Y transcripts harbor nonsense mutation within their open reading frames, yet most non-functional neo-Y genes are expressed throughout all of their length. We find evidence of gene loss of male-specific genes on the neo-X chromosome, and transcriptional silencing of testis-specific genes from the neo-X. CONCLUSIONS: Nonsense mediated decay (NMD) has been implicated to degrade transcripts containing pre-mature termination codons (PTC) in Drosophila, but rampant description of neo-Y genes with pre-mature stop codons suggests that it does not play a major role in down-regulating transcripts from the neo-Y. Loss or transcriptional down-regulation of genes from the neo-X with male-biased function provides evidence for beginning demasculinization of the neo-X. Thus, evolving sex chromosomes can rapidly shift their gene content or patterns of gene expression in response to their sex-biased transmission, supporting the idea that sex-specific or sexually antagonistic selection plays a major role in the evolution of heteromorphic sex chromosomes.

The epigenome of evolving Drosophila neo-sex chromosomes: dosage compensation and heterochromatin formation.

Sex chromosomes originated from autosomes but have evolved a highly specialized chromatin structure. Drosophila Y chromosomes are composed entirely of silent heterochromatin, while male X chromosomes have highly accessible chromatin and are hypertranscribed as a result of dosage compensation. Here, we dissect the molecular mechanisms and functional pressures driving heterochromatin formation and dosage compensation of the recently formed neo-sex chromosomes of Drosophila miranda. We show that the onset of heterochromatin formation on the neo-Y is triggered by an accumulation of repetitive DNA. The neo-X has evolved partial dosage compensation and we find that diverse mutational paths have been utilized to establish several dozen novel binding consensus motifs for the dosage compensation complex on the neo-X, including simple point mutations at pre-binding sites, insertion and deletion mutations, microsatellite expansions, or tandem amplification of weak binding sites. Spreading of these silencing or activating chromatin modifications to adjacent regions results in massive mis-expression of neo-sex linked genes, and little correspondence between functionality of genes and their silencing on the neo-Y or dosage compensation on the neo-X. Intriguingly, the genomic regions being targeted by the dosage compensation complex on the neo-X and those becoming heterochromatic on the neo-Y show little overlap, possibly reflecting different propensities along the ancestral chromosome that formed the sex chromosome to adopt active or repressive chromatin configurations. Our findings have broad implications for current models of sex chromosome evolution, and demonstrate how mechanistic constraints can limit evolutionary adaptations. Our study also highlights how evolution can follow predictable genetic trajectories, by repeatedly acquiring the same 21-bp consensus motif for recruitment of the dosage compensation complex, yet utilizing a diverse array of random mutational changes to attain the same phenotypic outcome.

Conservation and de novo acquisition of dosage compensation on newly evolved sex chromosomes in Drosophila.

Dosage compensation has arisen in response to the evolution of distinct male (XY) and female (XX) karyotypes. In Drosophila melanogaster, the MSL complex increases male X transcription approximately twofold. X-specific targeting is thought to occur through sequence-dependent binding to chromatin entry sites (CESs), followed by spreading in cis to active genes. We tested this model by asking how newly evolving sex chromosome arms in Drosophila miranda acquired dosage compensation. We found evidence for the creation of new CESs, with the analogous sequence and spacing as in D. melanogaster, providing strong support for the spreading model in the establishment of dosage compensation.

Sex-biased gene expression at homomorphic sex chromosomes in emus and its implication for sex chromosome evolution.

Sex chromosomes originate from autosomes. The accumulation of sexually antagonistic mutations on protosex chromosomes selects for a loss of recombination and sets in motion the evolutionary processes generating heteromorphic sex chromosomes. Recombination suppression and differentiation are generally viewed as the default path of sex chromosome evolution, and the occurrence of old, homomorphic sex chromosomes, such as those of ratite birds, has remained a mystery. Here, we analyze the genome and transcriptome of emu (Dromaius novaehollandiae) and confirm that most genes on the sex chromosome are shared between the Z and W. Surprisingly, however, levels of gene expression are generally sex-biased for all sex-linked genes relative to autosomes, including those in the pseudoautosomal region, and the male-bias increases after gonad formation. This expression bias suggests that the emu sex chromosomes have become masculinized, even in the absence of ZW differentiation. Thus, birds may have taken different evolutionary solutions to minimize the deleterious effects imposed by sexually antagonistic mutations: some lineages eliminate recombination along the protosex chromosomes to physically restrict sexually antagonistic alleles to one sex, whereas ratites evolved sex-biased expression to confine the product of a sexually antagonistic allele to the sex it benefits. This difference in conflict resolution may explain the preservation of recombining, homomorphic sex chromosomes in other lineages and illustrates the importance of sexually antagonistic mutations driving the evolution of sex chromosomes.

Nonrandom gene loss from the Drosophila miranda neo-Y chromosome.

A lack of recombination leads to the degeneration of an evolving Y chromosome. However, it is not known whether gene loss is largely a random process and primarily driven by the order in which mutations occur or whether certain categories of genes are lost less quickly than others; the latter would imply that selection counteracts the degeneration of Y chromosomes to some extent. In this study, we investigate the relationship between putative ancestral expression levels of neo-Y-linked genes in Drosophila miranda and their rates of degeneration. We use RNA-Seq data from its close relative Drosophila pseudoobscura to show that genes that have become nonfunctional on the D. miranda neo-Y had, on average, lower ancestral transcript levels and were expressed in fewer tissues compared with genes with intact reading frames. We also show that genes with male-biased expression are retained for longer on the neo-Y compared with female-biased genes. Our results imply that gene loss on the neo-Y is not a purely random, mutation-driven process. Instead, selection is--at least to some extent--preserving the function of genes that are more costly to lose, despite the strongly reduced efficacy of selection on the neo-Y chromosome.

Evolution of sex chromosomes in insects.

Sex chromosomes have many unusual features relative to autosomes. Y (or W) chromosomes lack genetic recombination, are male- (female-) limited, and show an abundance of genetically inert heterochromatic DNA but contain few functional genes. X (or Z) chromosomes also show sex-biased transmission (i.e., X chromosomes show female-biased and Z-chromosomes show male-biased inheritance) and are hemizygous in the heterogametic sex. Their unusual ploidy level and pattern of inheritance imply that sex chromosomes play a unique role in many biological processes and phenomena, including sex determination, epigenetic chromosome-wide regulation of gene expression, the distribution of genes in the genome, genomic conflict, local adaptation, and speciation. The vast diversity of sex chromosome systems in insects--ranging from the classical male heterogametic XY system in Drosophila to ZW systems in Lepidoptera or mobile genes determining sex as found in house flies--implies that insects can serve as unique model systems to study various functional and evolutionary aspects of these different processes.

Gene loss on Y chromosomes: a consequence of purifying selection?

Recombination restriction between evolving sex chromosomes leads to the degeneration of the chromosome that is present only in the heterogametic sex (the Y chromosome in XY species). The evolutionary forces driving Y chromosome degeneration, however, are still under debate and include positive and negative selection models. In a recent study, we showed that the rate of accumulation of loss-of-function mutations on the neo-Y chromosome of Drosophila miranda is compatible with the process of Muller's ratchet, the stochastic loss of the best mutational class of individuals from a small asexual population. Purifying selection at amino acid sites can accelerate the ratchet, and the speed of degeneration depends on the number of genes still present on the evolving Y chromosome. Our study shows that Y chromosome degeneration does not require the action of selective sweeps at linked sites, and can take place under realistic parameters of purifying selection only.

Nucleotide diversity in Silene latifolia autosomal and sex-linked genes.

The plant Silene latifolia has separate sexes and sex chromosomes, and is of interest for studying the early stages of sex chromosome evolution, especially the evolution of non-recombining regions on the Y chromosome. Hitch-hiking processes associated with ongoing genetic degeneration of the non-recombining Y chromosome are predicted to reduce Y-linked genes' effective population sizes, and S. latifolia Y-linked genes indeed have lower diversity than X-linked ones. We tested whether this represents a true diversity reduction on the Y, versus the alternative possibility, elevated diversity at X-linked genes, by collecting new data on nucleotide diversity for autosomal genes, which had previously been little studied. We find clear evidence that Y-linked genes have reduced diversity. However, another alternative explanation for a low Y effective size is a high variance in male reproductive success. Autosomal genes should then also have lower diversity than expected, relative to the X, but this is not found in our loci. Taking into account the higher mutation rate of Y-linked genes, their low sequence diversity indicates a strong effect of within-population hitch-hiking on the Y chromosome.

Muller's ratchet and the degeneration of the Drosophila miranda neo-Y chromosome.

Since its formation about 1.75 million years ago, the Drosophila miranda neo-Y chromosome has undergone a rapid process of degeneration, having lost approximately half of the genes that it originally contained. Using estimates of mutation rates and selection coefficients for loss-of-function mutations, we show that the high rate of accumulation of these mutations can largely be explained by Muller's ratchet, the process of stochastic loss of the least-loaded mutational class from a finite, nonrecombining population. We show that selection at nonsynonymous coding sites can accelerate the process of gene loss and that this effect varies with the number of genes still present on the degenerating neo-Y chromosome.

Slcyt, a newly identified sex-linked gene, has recently moved onto the X chromosome in Silene latifolia (Caryophyllaceae).

The sex chromosomes of the plant species Silene latifolia (white campion) are very young (only 5-10 My old), and all 11 X-linked genes so far described have Y-linked homologues. Theory predicts that X chromosomes should accumulate a nonrandom set of genes. However, little is known about the importance of gene movements between the X and the autosomes in plants, or in any very young sex chromosome system. Here, we isolate from cDNA a new gene, Slcyt, on the S. latifolia X, which encodes a cytochrome B protein. We genetically mapped SlCyt and found that it is located approximately 1 cM from the pseudoautosomal region. Genes in this region of the X chromosome have low divergence values from their homologous Y-linked genes, indicating that the X only recently stopped recombining with the Y. Genetic mapping in Silene vulgaris suggests that Slcyt originally belonged to a different linkage group from that of the other S. latifolia X-linked genes. Silene latifolia has no Y-linked homologue of Slcyt, and also no autosomal paralogues seem to exist. Slcyt moved from an autosome to the X very recently, as the Cyt gene is also X linked in Silene dioica, the sister species to S. latifolia, but is probably autosomal in Silene diclinis, implying that a translocation to the X probably occurred after the split between S. diclinis and S. latifolia/S. dioica. Diversity at Slcyt is extremely low (pi(syn) = 0.16%), and we find an excess of high frequency-derived variants and a negative Tajima's D, suggesting that the translocation was driven by selection.

The effects of deleterious mutations on evolution in non-recombining genomes.

Analyzing regions of the Drosophila genome that have low levels of genetic recombination helps us understand the prevalence of sexual reproduction. Here, we show that genetic variability in these regions can be explained by interference among strongly deleterious mutations and that selection becomes progressively less effective in influencing the behaviour of neighbouring sites as the number of closely linked sites on a chromosome increases.

Insertion events of CR1 retrotransposable elements elucidate the phylogenetic branching order in galliform birds.

Using standard phylogenetic methods, it can be hard to resolve the order in which speciation events took place when new lineages evolved in the distant past and within a short time frame. As an example, phylogenies of galliform birds (including well-known species such as chicken, turkey, and quail) usually show low bootstrap support values at short internal branches, reflecting the rapid diversification of these birds in the Eocene. However, given the key role of chicken and related poultry species in agricultural, evolutionary, general biological and disease studies, it is important to know their internal relationships. Recently, insertion patterns of transposable elements such as long and short interspersed nuclear element markers have proved powerful in revealing branching orders of difficult phylogenies. Here we decipher the order of speciation events in a group of 27 galliform species based on insertion events of chicken repeat 1 (CR1) transposable elements. Forty-four CR1 marker loci were identified from the draft sequence of the chicken genome, and from turkey BAC clone sequence, and the presence or absence of markers across species was investigated via electrophoretic size separation of amplification products and subsequent confirmation by DNA sequencing. Thirty markers proved possible to type with electrophoresis of which 20 were phylogenetically informative. The distribution of these repeat elements supported a single homoplasy-free cladogram, which confirmed that megapodes, cracids, New World quail, and guinea fowl form outgroups to Phasianidae and that quails, pheasants, and partridges are each polyphyletic groups. Importantly, we show that chicken is an outgroup to turkey and quail, an observation which does not have significant support from previous DNA sequence- and DNA-DNA hybridization-based trees and has important implications for evolutionary studies based on sequence or karyotype data from galliforms. We discuss the potential and limitations of using a genome-based retrotransposon approach in resolving problematic phylogenies among birds.

Nonrandom distribution of genes with sex-biased expression in the chicken genome.

Evolutionary theory predicts that sexually antagonistic genes should show a nonrandom genomic distribution with sex chromosomes usually being enriched for such genes. However, empirical observations from model organisms (Drosophila melanogaster, Caenorhabditis elegans, mammals) on the genomic location of genes with sex-biased expression have provided conflicting data and are not easily explained by a unified framework based on standard models of the evolution of sexually antagonistic genes. Previous studies have been confined to organisms with male heterogamety, meaning that effects related to homo- or heterozygosity of sex chromosomes cannot be separated from effects related to sex-specific characteristics. We therefore studied the genomic distribution of genes with sex-biased expression in the chicken, that is, in an organism with female heterogamety (males ZZ, females ZW). From the abundance of transcripts in expressed sequence tag libraries, we found an underrepresentation of female-specific genes (germ line and somatic tissue) and an overrepresentation of male-specific genes (somatic) on the Z chromosome. This is consistent with theoretical predictions only if mutations beneficial to one sex generally tend to be at least partly dominant (h > 0.5). We also note that sexual selection for a male-biased trait is facilitated by Z-linkage, because sons in organisms with female heterogamety will always inherit a Z chromosome from their fathers.