Cis-regulatory evolution in prokaryotes revealed by interspecific archaeal hybrids.

The study of allele-specific expression (ASE) in interspecific hybrids has played a central role in our understanding of a wide range of phenomena, including genomic imprinting, X-chromosome inactivation, and cis-regulatory evolution. However across the hundreds of studies of hybrid ASE, all have been restricted to sexually reproducing eukaryotes, leaving a major gap in our understanding of the genomic patterns of cis-regulatory evolution in prokaryotes. Here we introduce a method to generate stable hybrids between two species of halophilic archaea, and measure genome-wide ASE in these hybrids with RNA-seq. We found that over half of all genes have significant ASE, and that genes encoding kinases show evidence of lineage-specific selection on their cis-regulation. This pattern of polygenic selection suggested species-specific adaptation to low phosphate conditions, which we confirmed with growth experiments. Altogether, our work extends the study of ASE to archaea, and suggests that cis-regulation can evolve under polygenic lineage-specific selection in prokaryotes.

Noninvasive prenatal screening at low fetal fraction: comparing whole-genome sequencing and single-nucleotide polymorphism methods.

OBJECTIVE: Performance of noninvasive prenatal screening (NIPS) methodologies when applied to low fetal fraction samples is not well established. The single-nucleotide polymorphism (SNP) method fails samples below a predetermined fetal fraction threshold, whereas some laboratories employing the whole-genome sequencing (WGS) method report aneuploidy calls for all samples. Here, the performance of the two methods was compared to determine which approach actually detects more fetal aneuploidies. METHODS: Computational models were parameterized with up-to-date published data and used to compare the performance of the two methods at calling common fetal trisomies (T21, T18, T13) at low fetal fractions. Furthermore, clinical experience data were reviewed to determine aneuploidy detection rates based on compliance with recent invasive screening recommendations. RESULTS: The SNP method's performance is dependent on the origin of the trisomy, and is lowest for the most common trisomies (maternal M1 nondisjunction). Consequently, the SNP method cannot maintain acceptable performance at fetal fractions below ~3%. In contrast, the WGS method maintains high specificity independent of fetal fraction and has >80% sensitivity for trisomies in low fetal fraction samples. CONCLUSION: The WGS method will detect more aneuploidies below the fetal fraction threshold at which many labs issue a no-call result, avoiding unnecessary invasive procedures. © 2017 Counsyl Inc. Prenatal Diagnosis published by John Wiley & Sons, Ltd.

Dissecting the Genetic Basis of a Complex cis-Regulatory Adaptation.

Although single genes underlying several evolutionary adaptations have been identified, the genetic basis of complex, polygenic adaptations has been far more challenging to pinpoint. Here we report that the budding yeast Saccharomyces paradoxus has recently evolved resistance to citrinin, a naturally occurring mycotoxin. Applying a genome-wide test for selection on cis-regulation, we identified five genes involved in the citrinin response that are constitutively up-regulated in S. paradoxus. Four of these genes are necessary for resistance, and are also sufficient to increase the resistance of a sensitive strain when over-expressed. Moreover, cis-regulatory divergence in the promoters of these genes contributes to resistance, while exacting a cost in the absence of citrinin. Our results demonstrate how the subtle effects of individual regulatory elements can be combined, via natural selection, into a complex adaptation. Our approach can be applied to dissect the genetic basis of polygenic adaptations in a wide range of species.

Accounting for biases in riboprofiling data indicates a major role for proline in stalling translation.

The recent advent of ribosome profiling-sequencing of short ribosome-bound fragments of mRNA-has offered an unprecedented opportunity to interrogate the sequence features responsible for modulating translational rates. Nevertheless, numerous analyses of the first riboprofiling data set have produced equivocal and often incompatible results. Here we analyze three independent yeast riboprofiling data sets, including two with much higher coverage than previously available, and find that all three show substantial technical sequence biases that confound interpretations of ribosomal occupancy. After accounting for these biases, we find no effect of previously implicated factors on ribosomal pausing. Rather, we find that incorporation of proline, whose unique side-chain stalls peptide synthesis in vitro, also slows the ribosome in vivo. We also reanalyze a method that implicated positively charged amino acids as the major determinant of ribosomal stalling and demonstrate that it produces false signals of stalling in low-coverage data. Our results suggest that any analysis of riboprofiling data should account for sequencing biases and sparse coverage. To this end, we establish a robust methodology that enables analysis of ribosome profiling data without prior assumptions regarding which positions spanned by the ribosome cause stalling.

Transcript length mediates developmental timing of gene expression across Drosophila.

The time required to transcribe genes with long primary transcripts may limit their ability to be expressed in cells with short mitotic cycles, a phenomenon termed intron delay. As such short cycles are a hallmark of the earliest stages of insect development, we tested the impact of intron delay on the Drosophila developmental transcriptome. We find that long zygotically expressed genes show substantial delay in expression relative to their shorter counterparts, which is not observed for maternally deposited transcripts. Patterns of RNA-seq coverage along transcripts show that this delay is consistent with their inability to completely transcribe long transcripts, but not with transcriptional initiation-based regulatory control. We further show that highly expressed zygotic genes maintain compact transcribed regions across the Drosophila phylogeny, allowing conservation of embryonic expression patterns. We propose that the physical constraints of intron delay affect patterns of expression and the evolution of gene structure of a substantial portion of the Drosophila transcriptome.

Comparative validation of the D. melanogaster modENCODE transcriptome annotation.

Accurate gene model annotation of reference genomes is critical for making them useful. The modENCODE project has improved the D. melanogaster genome annotation by using deep and diverse high-throughput data. Since transcriptional activity that has been evolutionarily conserved is likely to have an advantageous function, we have performed large-scale interspecific comparisons to increase confidence in predicted annotations. To support comparative genomics, we filled in divergence gaps in the Drosophila phylogeny by generating draft genomes for eight new species. For comparative transcriptome analysis, we generated mRNA expression profiles on 81 samples from multiple tissues and developmental stages of 15 Drosophila species, and we performed cap analysis of gene expression in D. melanogaster and D. pseudoobscura. We also describe conservation of four distinct core promoter structures composed of combinations of elements at three positions. Overall, each type of genomic feature shows a characteristic divergence rate relative to neutral models, highlighting the value of multispecies alignment in annotating a target genome that should prove useful in the annotation of other high priority genomes, especially human and other mammalian genomes that are rich in noncoding sequences. We report that the vast majority of elements in the annotation are evolutionarily conserved, indicating that the annotation will be an important springboard for functional genetic testing by the Drosophila community.

Evolution at two levels of gene expression in yeast.

Despite the greater functional importance of protein levels, our knowledge of gene expression evolution is based almost entirely on studies of mRNA levels. In contrast, our understanding of how translational regulation evolves has lagged far behind. Here we have applied ribosome profiling--which measures both global mRNA levels and their translation rates--to two species of Saccharomyces yeast and their interspecific hybrid in order to assess the relative contributions of changes in mRNA abundance and translation to regulatory evolution. We report that both cis- and trans-acting regulatory divergence in translation are abundant, affecting at least 35% of genes. The majority of translational divergence acts to buffer changes in mRNA abundance, suggesting a widespread role for stabilizing selection acting across regulatory levels. Nevertheless, we observe evidence of lineage-specific selection acting on several yeast functional modules, including instances of reinforcing selection acting at both levels of regulation. Finally, we also uncover multiple instances of stop-codon readthrough that are conserved between species. Our analysis reveals the underappreciated complexity of post-transcriptional regulatory divergence and indicates that partitioning the search for the locus of selection into the binary categories of "coding" versus "regulatory" may overlook a significant source of selection, acting at multiple regulatory levels along the path from genotype to phenotype.

Mediation of Drosophila autosomal dosage effects and compensation by network interactions.

BACKGROUND: Gene dosage change is a mild perturbation that is a valuable tool for pathway reconstruction in Drosophila. While it is often assumed that reducing gene dose by half leads to two-fold less expression, there is partial autosomal dosage compensation in Drosophila, which may be mediated by feedback or buffering in expression networks. RESULTS: We profiled expression in engineered flies where gene dose was reduced from two to one. While expression of most one-dose genes was reduced, the gene-specific dose responses were heterogeneous. Expression of two-dose genes that are first-degree neighbors of one-dose genes in novel network models also changed, and the directionality of change depended on the response of one-dose genes. CONCLUSIONS: Our data indicate that expression perturbation propagates in network space. Autosomal compensation, or the lack thereof, is a gene-specific response, largely mediated by interactions with the rest of the transcriptome.

The developmental transcriptome of Drosophila melanogaster.

Drosophila melanogaster is one of the most well studied genetic model organisms; nonetheless, its genome still contains unannotated coding and non-coding genes, transcripts, exons and RNA editing sites. Full discovery and annotation are pre-requisites for understanding how the regulation of transcription, splicing and RNA editing directs the development of this complex organism. Here we used RNA-Seq, tiling microarrays and cDNA sequencing to explore the transcriptome in 30 distinct developmental stages. We identified 111,195 new elements, including thousands of genes, coding and non-coding transcripts, exons, splicing and editing events, and inferred protein isoforms that previously eluded discovery using established experimental, prediction and conservation-based approaches. These data substantially expand the number of known transcribed elements in the Drosophila genome and provide a high-resolution view of transcriptome dynamics throughout development.

Identification of functional elements and regulatory circuits by Drosophila modENCODE.

To gain insight into how genomic information is translated into cellular and developmental programs, the Drosophila model organism Encyclopedia of DNA Elements (modENCODE) project is comprehensively mapping transcripts, histone modifications, chromosomal proteins, transcription factors, replication proteins and intermediates, and nucleosome properties across a developmental time course and in multiple cell lines. We have generated more than 700 data sets and discovered protein-coding, noncoding, RNA regulatory, replication, and chromatin elements, more than tripling the annotated portion of the Drosophila genome. Correlated activity patterns of these elements reveal a functional regulatory network, which predicts putative new functions for genes, reveals stage- and tissue-specific regulators, and enables gene-expression prediction. Our results provide a foundation for directed experimental and computational studies in Drosophila and related species and also a model for systematic data integration toward comprehensive genomic and functional annotation.

Molecular evidence for increased regulatory conservation during metamorphosis, and against deleterious cascading effects of hybrid breakdown in Drosophila.

BACKGROUND: Speculation regarding the importance of changes in gene regulation in determining major phylogenetic patterns continues to accrue, despite a lack of broad-scale comparative studies examining how patterns of gene expression vary during development. Comparative transcriptional profiling of adult interspecific hybrids and their parental species has uncovered widespread divergence of the mechanisms controlling gene regulation, revealing incompatibilities that are masked in comparisons between the pure species. However, this has prompted the suggestion that misexpression in adult hybrids results from the downstream cascading effects of a subset of genes improperly regulated in early development. RESULTS: We sought to determine how gene expression diverges over development, as well as test the cascade hypothesis, by profiling expression in males of Drosophila melanogaster, D. sechellia, and D. simulans, as well as the D. simulans (female) x D. sechellia (male) male F1 hybrids, at four different developmental time points (3rd instar larval, early pupal, late pupal, and newly-emerged adult). Contrary to the cascade model of misexpression, we find that there is considerable stage-specific autonomy of regulatory breakdown in hybrids, with the larval and adult stages showing significantly more hybrid misexpression as compared to the pupal stage. However, comparisons between pure species indicate that genes expressed during earlier stages of development tend to be more conserved in terms of their level of expression than those expressed during later stages, suggesting that while Von Baer's famous law applies at both the level of nucleotide sequence and expression, it may not apply necessarily to the underlying overall regulatory network, which appears to diverge over the course of ontogeny and which can only be ascertained by combining divergent genomes in species hybrids. CONCLUSION: Our results suggest that complex integration of regulatory circuits during morphogenesis may lead to it being more refractory to divergence of underlying gene regulatory mechanisms--more than that suggested by the conservation of gene expression levels between species during earlier stages. This provides support for a 'developmental hourglass' model of divergence of gene expression in Drosophila resulting in a highly conserved pupal stage.

Male sex drive and the maintenance of sex: evidence from Drosophila.

The resolution of the paradoxes surrounding the evolutionary origins and maintenance of sexual reproduction has been a major focus in biology. The operation of sexual selection-which is very common among multicellular organisms-has been proposed as an important factor in the maintenance of sex, though in order for this hypothesis to hold, the strength of sexual selection must be stronger in males than in females. Sexual selection poses its own series of evolutionary questions, including how genetic variability is maintained in the face of sustained directional selection (known as the "paradox of the lek"). In this short review, we present evidence obtained from recent comparative genomics projects arguing that 1) the genomic consequences of sexual selection clearly show that its effect is stronger in males and 2) this sustained selection over evolutionary timescales also has an effect of capturing de novo genes and expression patterns influencing male fitness, thus providing a mechanism via which new genetic variation can be input into to male traits. Furthermore, we argue that this latter process of genomic "masculinization" has an additional effect of making males difficult to purge from populations, as evidence from Drosophila indicates that, for example, many male sexually selected seminal fluid factors are required to ensure maximally efficient reproduction. Newly arising parthenogenic mutations would suffer an immediate reproductive rate disadvantage were these proteins lost. We show that recent studies confirm that genomic masculinization, as a result of "male sex drive," has important consequences for the evolution of sexually dimorphic species.

Demystifying phenotypes: The comparative genomics of evo-devo.

Developmental geneticists have spearheaded the synthesis of evolutionary and developmental biology, a.k.a 'evo-devo', leading to a wealth of recent insights about how morphological diversity has evolved. However, there exists a gap between these disciplines, and evo-devo has not benefited from an integration of the principles derived from population genetics and molecular evolution. In order to contribute to the remediation of this deficiency, we recently performed a study investigating how genes diverge among closely related species of Drosophila as a function of when they are expressed during development. We found that patterns of genetic divergence parallels morphology: interspecific divergence accumulates as development progresses. We also sought to test whether this positive gradient of divergence over ontogeny is best explained by purifying selection constraining the divergence of early-expressed genes or positive selection driving the evolution of those expressed later. Interestingly, we found evidence that both processes occur simultaneously. We argue that comparative genomics approaches, by juxtaposing gene- and phenotypelevel divergence, particularly among closely related species, have much to contribute to the ongoing evo-devo synthesis, complementing traditional genetics-based techniques with largescale screening analyses uncovering the mechanisms underlying developmental change.

Ontogeny and phylogeny: molecular signatures of selection, constraint, and temporal pleiotropy in the development of Drosophila.

BACKGROUND: Karl Ernst Von Baer noted that species tend to show greater morphological divergence in later stages of development when compared to earlier stages. Darwin originally interpreted these observations via a selectionist framework, suggesting that divergence should be greatest during ontogenic stages in which organisms experienced varying 'conditions of existence' and opportunity for differential selection. Modern hypotheses have focused on the notion that genes and structures involved in early development will be under stronger purifying selection due to the deleterious pleiotropic effects of mutations propagating over the course of ontogeny, also known as the developmental constraint hypothesis. RESULTS: Using developmental stage-specific expressed sequence tag (EST) libraries, we tested the 2 hypotheses by comparing the rates of evolution of 7,180 genes obtained from 6 species of the Drosophila melanogaster group with respect to ontogeny, and sex and reproduction-related functions in gonadal tissues. Supporting morphological observations, we found evidence of a pattern of increasing mean evolutionary rate in genes that are expressed in subsequent stages of development. Furthermore, supporting expectations that early expressed genes are constrained in divergence, we found that embryo stage genes are involved in a higher mean number of interactions as compared to later stages. We noted that the accelerated divergence of genes in the adult stage is explained by those expressed specifically in the male gonads, whose divergence is driven by positive selection. In addition, accelerated gonadal gene divergence occurs only in the adult stage, suggesting that the effects of selection are observed primarily at the stages during which they are expected occur. Finally, we also found a significant correlation between temporal specificity of gene expression and evolutionary rate, supporting expectations that genes with ubiquitous expression are under stronger constraint. CONCLUSION: Taken together, these results support both the developmental constraint hypothesis limiting the divergence of early expressed developmentally important genes, leading to a gradient of divergence rates over ontogeny (embryonic < larval/pupal < adult), as well as Darwin's 'selection opportunity' hypothesis leading to increased divergence in adults, particularly in the case of reproductive tissues. We suggest that a constraint early/opportunity late model best explains divergence over ontogeny.

Comparative analysis of function and interaction of transcription factors in nematodes: extensive conservation of orthology coupled to rapid sequence evolution.

BACKGROUND: Much of the morphological diversity in eukaryotes results from differential regulation of gene expression in which transcription factors (TFs) play a central role. The nematode Caenorhabditis elegans is an established model organism for the study of the roles of TFs in controlling the spatiotemporal pattern of gene expression. Using the fully sequenced genomes of three Caenorhabditid nematode species as well as genome information from additional more distantly related organisms (fruit fly, mouse, and human) we sought to identify orthologous TFs and characterized their patterns of evolution. RESULTS: We identified 988 TF genes in C. elegans, and inferred corresponding sets in C. briggsae and C. remanei, containing 995 and 1093 TF genes, respectively. Analysis of the three gene sets revealed 652 3-way reciprocal 'best hit' orthologs (nematode TF set), approximately half of which are zinc finger (ZF-C2H2 and ZF-C4/NHR types) and HOX family members. Examination of the TF genes in C. elegans and C. briggsae identified the presence of significant tandem clustering on chromosome V, the majority of which belong to ZF-C4/NHR family. We also found evidence for lineage-specific duplications and rapid evolution of many of the TF genes in the two species. A search of the TFs conserved among nematodes in Drosophila melanogaster, Mus musculus and Homo sapiens revealed 150 reciprocal orthologs, many of which are associated with important biological processes and human diseases. Finally, a comparison of the sequence, gene interactions and function indicates that nematode TFs conserved across phyla exhibit significantly more interactions and are enriched in genes with annotated mutant phenotypes compared to those that lack orthologs in other species. CONCLUSION: Our study represents the first comprehensive genome-wide analysis of TFs across three nematode species and other organisms. The findings indicate substantial conservation of transcription factors even across distant evolutionary lineages and form the basis for future experiments to examine TF gene function in nematodes and other divergent phyla.

Sexual selection and maintenance of sex: evidence from comparisons of rates of genomic accumulation of mutations and divergence of sex-related genes in sexual and hermaphroditic species of Caenorhabditis.

Several hypotheses have been proposed to explain the persistence of dioecy despite the reproductive advantages conferred to hermaphrodites, including greater efficiency at purging deleterious mutations in the former. Dioecy can benefit from both mutation purging and accelerated evolution by bringing together beneficial mutations in the same individual via recombination and shuffling of genotypes. In addition, mathematical treatment has shown that sexual selection is also capable of mitigating the cost of maintaining separate sexes by increasing the overall fitness of sexual populations, and genomic comparisons have shown that sexual selection can lead to accelerated evolution. Here, we examine the advantages of dioecy versus hermaphroditism by comparing the rate of evolution in sex-related genes and the rate of accumulation of deleterious mutations using a large number of orthologs (11,493) in the dioecious Caenorhabditis remanei and the hermaphroditic Caenorhabditis briggsae. We have used this data set to estimate the deleterious mutation rate per generation, U, in both species and find that although it is significantly higher in hermaphrodites, both species are at least 2 orders of magnitude lower than the value required to explain the persistence of sex by efficiency at purging deleterious mutations alone. We also find that genes expressed in sperm are evolving rapidly in both species; however, they show a greater increase in their rate of evolution relative to genes expressed in other tissues in C. remanei, suggesting stronger sexual selection pressure acting on these genes in dioecious species. Interestingly, the persistence of a signal of rapid evolution of sperm genes in C. briggsae suggests a recent evolutionary origin of hermaphrodism in this lineage. Our results provide empirical evidence of increased sexual selection pressure in dioecious animals, supporting the possibility that sexual selection may play an important role in the maintenance of sexual reproduction.

Evolution of genes and genomes on the Drosophila phylogeny.

Comparative analysis of multiple genomes in a phylogenetic framework dramatically improves the precision and sensitivity of evolutionary inference, producing more robust results than single-genome analyses can provide. The genomes of 12 Drosophila species, ten of which are presented here for the first time (sechellia, simulans, yakuba, erecta, ananassae, persimilis, willistoni, mojavensis, virilis and grimshawi), illustrate how rates and patterns of sequence divergence across taxa can illuminate evolutionary processes on a genomic scale. These genome sequences augment the formidable genetic tools that have made Drosophila melanogaster a pre-eminent model for animal genetics, and will further catalyse fundamental research on mechanisms of development, cell biology, genetics, disease, neurobiology, behaviour, physiology and evolution. Despite remarkable similarities among these Drosophila species, we identified many putatively non-neutral changes in protein-coding genes, non-coding RNA genes, and cis-regulatory regions. These may prove to underlie differences in the ecology and behaviour of these diverse species.

Evolution in the fast lane: rapidly evolving sex-related genes in Drosophila.

A large portion of the annotated genes in Drosophila melanogaster show sex-biased expression, indicating that sex and reproduction-related genes (SRR genes) represent an appreciable component of the genome. Previous studies, in which subsets of genes were compared among few Drosophila species, have found that SRR genes exhibit unusual evolutionary patterns. Here, we have used the newly released genome sequences from 12 Drosophila species, coupled to a larger set of SRR genes, to comprehensively test the generality of these patterns. Among 2505 SRR genes examined, including ESTs with biased expression in reproductive tissues and genes characterized as involved in gametogenesis, we find that a relatively high proportion of SRR genes have experienced accelerated divergence throughout the genus Drosophila. Several testis-specific genes, male seminal fluid proteins (SFPs), and spermatogenesis genes show lineage-specific bursts of accelerated evolution and positive selection. SFP genes also show evidence of lineage-specific gene loss and/or gain. These results bring us closer to understanding the details of the evolutionary dynamics of SRR genes with respect to species divergence.

Association between levels of coding sequence divergence and gene misregulation in Drosophila male hybrids.

Previous studies have shown widespread conservation of gene expression levels between species of the Drosophila melanogaster subgroup as well as a positive correlation between coding sequence divergence and expression level divergence between species. Meanwhile, large-scale misregulation of gene expression level has been described in interspecific sterile hybrids between D. melanogaster, D. simulans, D. mauritiana, and D. sechellia. Using data from gene expression analysis involving D. simulans, D. melanogaster, and their hybrids, we observed a significant positive correlation between protein sequence divergence and gene expression differences between hybrids and their parental species. Furthermore, we demonstrate that underexpressed misregulated genes in hybrids are evolving more rapidly at the protein sequence level than nonmisregulated genes or overexpressed misregulated genes, highlighting the possible effects of sexual and natural selection as male-biased genes and nonessential genes are the principal gene categories affected by interspecific hybrid misregulation.

A physical map of the genome of Atlantic salmon, Salmo salar.

A physical map of the Atlantic salmon (Salmo salar) genome was generated based on HindIII fingerprints of a publicly available BAC (bacterial artificial chromosome) library constructed from DNA isolated from a Norwegian male. Approximately 11.5 haploid genome equivalents (185,938 clones) were successfully fingerprinted. Contigs were first assembled via FPC using high-stringency (1e-16), and then end-to-end joins yielded 4354 contigs and 37,285 singletons. The accuracy of the contig assembly was verified by hybridization and PCR analysis using genetic markers. A subset of the BACs in the library contained few or no HindIII recognition sites in their insert DNA. BglI digestion fragment patterns of these BACs allowed us to identify three classes: (1) BACs containing histone genes, (2) BACs containing rDNA-repeating units, and (3) those that do not have BglI recognition sites. End-sequence analysis of selected BACs representing these three classes confirmed the identification of the first two classes and suggested that the third class contained highly repetitive DNA corresponding to tRNAs and related sequences.