Sergey V. Nuzhdin

Genomic Islands of Speciation in Anopheles gambiae

The African malaria mosquito, Anopheles gambiae sensu stricto (A. gambiae), provides a unique opportunity to study the evolution of reproductive isolation because it is divided into two sympatric, partially isolated subtaxa known as M form and S form. With the annotated genome of this species now available, high-throughput techniques can be applied to locate and characterize the genomic regions contributing to reproductive isolation. In order to quantify patterns of differentiation within A. gambiae, we hybridized population samples of genomic DNA from each form to Affymetrix GeneChip microarrays. We found that three regions, together encompassing less than 2.8 Mb, are the only locations where the M and S forms are significantly differentiated. Two of these regions are adjacent to centromeres, on Chromosomes 2L and X, and contain 50 and 12 predicted genes, respectively. Sequenced loci in these regions contain fixed differences between forms and no shared polymorphisms, while no fixed differences were found at nearby control loci. The third region, on Chromosome 2R, contains only five predicted genes; fixed differences in this region were also verified by direct sequencing. These “speciation islands” remain differentiated despite considerable gene flow, and are therefore expected to contain the genes responsible for reproductive isolation. Much effort has recently been applied to locating the genes and genetic changes responsible for reproductive isolation between species. Though much can be inferred about speciation by studying taxa that have diverged for millions of years, studying differentiation between taxa that are in the early stages of isolation will lead to a clearer view of the number and size of regions involved in the genetics of speciation. Despite appreciable levels of gene flow between the M and S forms of A. gambiae, we were able to isolate three small regions of differentiation where genes responsible for ecological and behavioral isolation are likely to be located. We expect reproductive isolation to be due to changes at a small number of loci, as these regions together contain only 67 predicted genes. Concentrating future mapping experiments on these regions should reveal the genes responsible for reproductive isolation between forms.

RNA-seq: technical variability and sampling.

BackgroundRNA-seq is revolutionizing the way we study transcriptomes. mRNA can be surveyed without prior knowledge of gene transcripts. Alternative splicing of transcript isoforms and the identification of previously unknown exons are being reported. Initial reports of differences in exon usage, and splicing between samples as well as quantitative differences among samples are beginning to surface. Biological variation has been reported to be larger than technical variation. In addition, technical variation has been reported to be in line with expectations due to random sampling. However, strategies for dealing with technical variation will differ depending on the magnitude. The size of technical variance, and the role of sampling are examined in this manuscript.ResultsIn this study three independent Solexa/Illumina experiments containing technical replicates are analyzed. When coverage is low, large disagreements between technical replicates are apparent. Exon detection between technical replicates is highly variable when the coverage is less than 5 reads per nucleotide and estimates of gene expression are more likely to disagree when coverage is low. Although large disagreements in the estimates of expression are observed at all levels of coverage.ConclusionsTechnical variability is too high to ignore. Technical variability results in inconsistent detection of exons at low levels of coverage. Further, the estimate of the relative abundance of a transcript can substantially disagree, even when coverage levels are high. This may be due to the low sampling fraction and if so, it will persist as an issue needing to be addressed in experimental design even as the next wave of technology produces larger numbers of reads. We provide practical recommendations for dealing with the technical variability, without dramatic cost increases.

Sure facts, speculations, and open questions about the evolution of transposable element copy number

Transposable elements (TEs) are sequences capable of multiplying in their hosts genome. They survive by increasing copy numbers due to transpositions, and natural selection washes them out because hosts with heavier loads of TEs have lower fitness. The available phylogenetic evidence supports the view that TEs have existed in living organisms for hundreds of millions of years. A fundamental question facing the field is how can an equilibrium be attained between transposition and selection which allows these parasitic genetic elements to persist for such a long time period? To answer this question, it is necessary to understand how the rate of TE transposition is controlled and to describe the mechanisms with which natural selection opposes TE accumulation. Perhaps the best models for such a study are copia and gypsy retrotransposons in Drosophila. Their average rate of transposition in nature is between 10−5 − 10−4 transpositions per copy per generation. Unlike nature, transposition rates vary widely, from zero to 10−2, between laboratory lines. This variability in transposition rate is controlled by host genes. It is probable that in nature TE site heterogeneity is caused by frequent transpositions in rare flies with permissive alleles, and no transpositions happen in the rest of flies. The average rate of TE transposition in nature may thus depend on the frequency of permissive alleles, which is a function of the rate of mutation from restrictive to permissive alleles, the mechanism and the strength of selection opposing TE multiplication, and population size. Thus, evolution of the frequency of permissive alleles of genes controlling transposition must be accounted for to understand evolution of TE copy numbers.

Natural Genetic Variation in Cuticular Hydrocarbon Expression in Male and Female Drosophila melanogaster

Cuticular hydrocarbons (CHCs) act as contact pheromones in Drosophila melanogaster and are an important component of several ecological traits. Segregating genetic variation in the expression of CHCs at the population level in D. melanogaster is likely to be important for mate choice and climatic adaptation; however, this variation has never been characterized. Using a panel of recombinant inbred lines (RILs) derived from a natural population, we found significant between-line variation for nearly all CHCs in both sexes. We identified 25 QTL in females and 15 QTL in males that pleiotropically influence CHC expression. There was no evidence of colocalization of QTL for homologous traits across the sexes, indicating that sexual dimorphism and low intersex genetic correlations between homologous CHCs are a consequence of largely independent genetic control. This is consistent with a pattern of divergent sexual and natural selection between the sexes.

Regulatory Divergence in Drosophila melanogaster and D. simulans, a Genomewide Analysis of Allele-Specific Expression

Species-specific regulation of gene expression contributes to the development and maintenance of reproductive isolation and to species differences in ecologically important traits. A better understanding of the evolutionary forces that shape regulatory variation and divergence can be developed by comparing expression differences among species and interspecific hybrids. Once expression differences are identified, the underlying genetics of regulatory variation or divergence can be explored. With the goal of associating cis and/or trans components of regulatory divergence with differences in gene expression, overall and allele-specific expression levels were assayed genomewide in female adult heads of Drosophila melanogaster, D. simulans, and their F1 hybrids. A greater proportion of cis differences than trans differences were identified for genes expressed in heads and, in accordance with previous studies, cis differences also explained a larger number of species differences in overall expression level. Regulatory divergence was found to be prevalent among genes associated with defense, olfaction, and among genes downstream of the Drosophila sex determination hierarchy. In addition, two genes, with critical roles in sex determination and micro RNA processing, Sxl and loqs, were identified as misexpressed in hybrid female heads, potentially contributing to hybrid incompatibility.

GENOMIC BASIS OF AGING AND LIFE HISTORY EVOLUTION IN DROSOPHILA MELANOGASTER

Natural diversity in aging and other life‐history patterns is a hallmark of organismal variation. Related species, populations, and individuals within populations show genetically based variation in life span and other aspects of age‐related performance. Population differences are especially informative because these differences can be large relative to within‐population variation and because they occur in organisms with otherwise similar genomes. We used experimental evolution to produce populations divergent for life span and late‐age fertility and then used deep genome sequencing to detect sequence variants with nucleotide‐level resolution. Several genes and genome regions showed strong signatures of selection, and the same regions were implicated in independent comparisons, suggesting that the same alleles were selected in replicate lines. Genes related to oogenesis, immunity, and protein degradation were implicated as important modifiers of late‐life performance. Expression profiling and functional annotation narrowed the list of strong candidate genes to 38, most of which are novel candidates for regulating aging. Life span and early age fecundity were negatively correlated among populations; therefore, the alleles we identified also are candidate regulators of a major life‐history trade‐off. More generally, we argue that hitchhiking mapping can be a powerful tool for uncovering the molecular bases of quantitative genetic variation.

A Versatile Method for Cell-Specific Profiling of Translated mRNAs in Drosophila

In Drosophila melanogaster few methods exist to perform rapid cell-type or tissue-specific expression profiling. A translating ribosome affinity purification (TRAP) method to profile actively translated mRNAs has been developed for use in a number of multicellular organisms although it has only been implemented to examine limited sets of cell- or tissue-types in these organisms. We have adapted the TRAP method for use in the versatile GAL4/UAS system of Drosophila allowing profiling of almost any tissue/cell-type with a single genetic cross. We created transgenic strains expressing a GFP-tagged ribosomal protein, RpL10A, under the control of the UAS promoter to perform cell-type specific translatome profiling. The GFP::RpL10A fusion protein incorporates efficiently into ribosomes and polysomes. Polysome affinity purification strongly enriches mRNAs from expected genes in the targeted tissues with sufficient sensitivity to analyze expression in small cell populations. This method can be used to determine the unique translatome profiles in different cell-types under varied physiological, pharmacological and pathological conditions.

Somatic sex-specific transcriptome differences in Drosophila revealed by whole transcriptome sequencing

BackgroundUnderstanding animal development and physiology at a molecular-biological level has been advanced by the ability to determine at high resolution the repertoire of mRNA molecules by whole transcriptome resequencing. This includes the ability to detect and quantify rare abundance transcripts and isoform-specific mRNA variants produced from a gene.The sex hierarchy consists of a pre-mRNA splicing cascade that directs the production of sex-specific transcription factors that specify nearly all sexual dimorphism. We have used deep RNA sequencing to gain insight into how the Drosophila sex hierarchy generates somatic sex differences, by examining gene and transcript isoform expression differences between the sexes in adult head tissues.ResultsHere we find 1,381 genes that differ in overall expression levels and 1,370 isoform-specific transcripts that differ between males and females. Additionally, we find 512 genes not regulated downstream of transformer that are significantly more highly expressed in males than females. These 512 genes are enriched on the × chromosome and reside adjacent to dosage compensation complex entry sites, which taken together suggests that their residence on the × chromosome might be sufficient to confer male-biased expression. There are no transcription unit structural features, from a set of features, that are robustly significantly different in the genes with significant sex differences in the ratio of isoform-specific transcripts, as compared to random isoform-specific transcripts, suggesting that there is no single molecular mechanism that generates isoform-specific transcript differences between the sexes, even though the sex hierarchy is known to include three pre-mRNA splicing factors.ConclusionsWe identify thousands of genes that show sex-specific differences in overall gene expression levels, and identify hundreds of additional genes that have differences in the abundance of isoform-specific transcripts. No transcription unit structural feature was robustly enriched in the sex-differentially expressed transcript isoforms. Additionally, we found that many genes with male-biased expression were enriched on the × chromosome and reside adjacent to dosage compensation entry sites, suggesting that differences in sex chromosome composition contributes to dimorphism in gene expression. Taken together, this study provides new insight into the molecular underpinnings of sexual differentiation.

Sex-Specific Splicing in Drosophila: Widespread Occurrence, Tissue Specificity and Evolutionary Conservation

Many genes in eukaryotic genomes produce multiple transcripts through a variety of molecular mechanisms including alternative splicing. Alternatively spliced transcripts often encode functionally distinct proteins, indicating that gene regulation at this level makes an important contribution to organismal complexity. The multilevel splicing cascade that regulates sex determination and sex-specific development in Drosophila is a classical example of the role of alternative splicing in cell differentiation. Recent evidence suggests that a large proportion of genes in the Drosophila genome may be spliced in a sex-biased fashion, raising the possibility that alternative splicing may play a more general role in sexually dimorphic development and physiology. However, the prevalence of sex-specific splicing and the extent to which it is shared among genotypes are not fully understood. Genetic variation in the splicing of key components of the sex determination pathway is known to influence the expression of downstream target genes, suggesting that alternative splicing at other loci may also vary in functionally important ways. In this study, we used exon-specific microarrays to examine 417 multitranscript genes for evidence of sex-specific and genotype-specific splicing in 80 different genotypes of Drosophila melanogaster. Most of these loci showed sex-biased splicing, whereas genotype-specific splicing was rare. One hundred thirty-five genes showed different alternative transcript use in males vs. females. Real-time PCR analysis of 6 genes chosen to represent a broad range of biological functions showed that most sex-biased splicing occurs in the gonads. However, somatic tissues, particularly adult heads, also show evidence of sex-specific splicing. Comparison of splicing patterns at orthologous loci in seven Drosophila species shows that sexual biases in alternative exon representation are highly conserved, indicating that sex-specific splicing is an ancient feature of Drosophila biology. To investigate potential mechanisms of sex-biased splicing, we used real-time PCR to examine the expression of six known regulators of alternative splicing in males vs. females. We found that all six loci are themselves spliced sex specifically in gonads and heads, suggesting that regulatory hierarchies based on alternative splicing may be an important feature of sexual differentiation.

De Novo Assembly of the Manila Clam Ruditapes philippinarum Transcriptome Provides New Insights into Expression Bias, Mitochondrial Doubly Uniparental Inheritance and Sex Determination

Males and females share the same genome, thus, phenotypic divergence requires differential gene expression and sex-specific regulation. Accordingly, the analysis of expression patterns is pivotal to the understanding of sex determination mechanisms. Many bivalves are stable gonochoric species, but the mechanism of gonad sexualization and the genes involved are still unknown. Moreover, during the period of sexual rest, a gonad is not present and sex cannot be determined. A mechanism associated with germ line differentiation in some bivalves, including the Manila clam Ruditapes philippinarum, is the doubly uniparental inheritance (DUI) of mitochondria, a variation of strict maternal inheritance. Two mitochondrial lineages are present, one transmitted through eggs and the other through sperm, as well as a mother-dependent sex bias of the progeny. We produced a de novo annotation of 17,186 transcripts from R. philippinarum and compared the transcriptomes of males and females and identified 1,575 genes with strong sex-specific expression and 166 sex-specific single nucleotide polymorphisms, obtaining preliminary information about genes that could be involved in sex determination. Then we compared the transcriptomes between a family producing predominantly females and a family producing predominantly males to identify candidate genes involved in regulation of sex-specific aspects of DUI system, finding a relationship between sex bias and differential expression of several ubiquitination genes. In mammalian embryos, sperm mitochondria are degraded by ubiquitination. A modification of this mechanism is hypothesized to be responsible for the retention of sperm mitochondria in male embryos of DUI species. Ubiquitination can additionally regulate gene expression, playing a role in sex determination of several animals. These data enable us to develop a model that incorporates both the DUI literature and our new findings.