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

BackgroundThe 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.ResultsWe 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.ConclusionsNonsense 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

Molecular evolution under low recombination

Analyzing regions in the genome with low levels of recombination helps understand the prevalence of sexual reproduction. Here, I show that variability in regions of reduced recombination in Drosophila can be explained by interference among strongly deleterious mutations; selection becomes progressively less effective in influencing the behaviour of neighbouring sites as the number of closely linked sites on a chromosome increases. I also show that the accumulation of loss-of-function mutations on the neo-Y chromosome of Drosophila miranda is compatible with a model of selection against such mutations alone, without the need to invoke the action of selective sweeps. I describe the discovery of two new sex-linked genes in the plant Silene latifolia, SlCyt and SlX9/SlY9. SlCyt has been recently translocated from an autosome to the X and shows signs of a selective sweep. Its possible role in having caused recombination arrest between the evolving X and Y chromosome is discussed. SlX9 still has an intact Y-linked copy that is presumably functional. Nucleotide diversity at SlY9 is very low, whereas SlX9 has an unusually high diversity and shows signs of introgression from S. dioica into S. latifolia, but the effect of this seems very localized

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

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

Exploring Large-scale Structure with Billions of Galaxies

We consider cosmological applications of galaxy number density correlations to be inferred from future deep and wide multi-band optical surveys. We mostly focus on very large scales as a probe of possible features in the primordial power spectrum. We find the proposed survey of the Large Synoptic Survey Telescope may be competitive with future all-sky CMB experiments over a broad range of scales. On very large scales the inferred power spectrum is robust to photometric redshift errors, and, given a sufficient number density of galaxies, to angular variations in dust extinction and photometric calibration errors. We also consider other applications, such as constraining dark energy with the two CMB-calibrated standard rulers in the matter power spectrum, and controlling the effect of photometric redshift errors to facilitate the interpretation of cosmic shear data. We find that deep photometric surveys over wide area can provide constraints that are competitive with spectroscopic surveys in small volumes.Comment: 11 pages, 7 figures, ApJ accepted, references added, expanded discussion in Sec. 3.

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

The cellular and synaptic architecture of the mechanosensory dorsal horn

The deep dorsal horn is a poorly characterized spinal cord region implicated in processing low-threshold mechanoreceptor (LTMR) information. We report an array of mouse genetic tools for defining neuronal components and functions of the dorsal horn LTMR-recipient zone (LTMR-RZ), a role for LTMR-RZ processing in tactile perception, and the basic logic of LTMR-RZ organization. We found an unexpectedly high degree of neuronal diversity in the LTMR-RZ: seven excitatory and four inhibitory subtypes of interneurons exhibiting unique morphological, physiological, and synaptic properties. Remarkably, LTMRs form synapses on between four and 11 LTMR-RZ interneuron subtypes, while each LTMR-RZ interneuron subtype samples inputs from at least one to three LTMR classes, as well as spinal cord interneurons and corticospinal neurons. Thus, the LTMR-RZ is a somatosensory processing region endowed with a neuronal complexity that rivals the retina and functions to pattern the activity of ascending touch pathways that underlie tactile perception

Pervasive lesion segregation shapes cancer genome evolution

Cancers arise through the acquisition of oncogenic mutations and grow through clonal expansion. Here we reveal that most mutagenic DNA lesions are not resolved as mutations 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 characterise 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 multi-allelic and combinatorial genetic diversity. The phasing of lesions enables the 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.This work was supported by: Cancer Research UK (20412, 22398), the European Research Council (615584, 682398), the Wellcome Trust (WT108749/Z/15/Z, WT106563/Z/14/A, WT202878/B/16/Z), the European Molecular Biology Laboratory, the MRC Human Genetics Unit core funding programme grants (MC_UU_00007/11, MC_UU_00007/16), and the ERDF/Spanish Ministry of Science, Innovation and Universities-Spanish State Research Agency/DamReMap Project (RTI2018-094095-B-I00)

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

Abstract 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

The state of the Martian climate

60°N was +2.0°C, relative to the 1981–2010 average value (Fig. 5.1). This marks a new high for the record. The average annual surface air temperature (SAT) anomaly for 2016 for land stations north of starting in 1900, and is a significant increase over the previous highest value of +1.2°C, which was observed in 2007, 2011, and 2015. Average global annual temperatures also showed record values in 2015 and 2016. Currently, the Arctic is warming at more than twice the rate of lower latitudes