PoPoolation DB: a user-friendly web-based database for the retrieval of natural polymorphisms in Drosophila

Sox21b underlies the rapid diversification of a novel male genital structure between Drosophila species

The emergence and diversification of morphological novelties is a major feature of animal evolution.1,2,3,4,5,6,7,8,9 However, relatively little is known about the genetic basis of the evolution of novel structures and the mechanisms underlying their diversification. The epandrial posterior lobes of male genitalia are a novelty of particular Drosophila species.10,11,12,13 The lobes grasp the female ovipositor and insert between her abdominal tergites and, therefore, are important for copulation and species recognition.10,11,12,14,15,16,17 The posterior lobes likely evolved from co-option of a Hox-regulated gene network from the posterior spiracles10 and have since diversified in morphology in the D. simulans clade, in particular, over the last 240,000 years, driven by sexual selection.18,19,20,21 The genetic basis of this diversification is polygenic but, to the best of our knowledge, none of the causative genes have been identified.22,23,24,25,26,27,28,29,30 Identifying the genes underlying the diversification of these secondary sexual structures is essential to understanding the evolutionary impact on copulation and species recognition. Here, we show that Sox21b negatively regulates posterior lobe size. This is consistent with expanded Sox21b expression in D. mauritiana, which develops smaller posterior lobes than D. simulans. We tested this by generating reciprocal hemizygotes and confirmed that changes in Sox21b underlie posterior lobe evolution between these species. Furthermore, we found that posterior lobe size differences caused by the species-specific allele of Sox21b significantly affect copulation duration. Taken together, our study reveals the genetic basis for the sexual-selection-driven diversification of a novel morphological structure and its functional impact on copulatory behavior.

Genetic analyses uncover pleiotropic compensatory roles for Drosophila Nucleobindin-1 in inositol trisphosphate-mediated intracellular calcium homeostasis

Nucleobindin-1 is an EF-hand calcium-binding protein with a distinctive profile, predominantly localized to the Golgi in insect and wide-ranging vertebrate cell types, alike. Its putative involvements in intracellular calcium (Ca²⁺) homeostasis have never been phenotypically characterized in any model organism. We have analyzed an adult-viable mutant that completely disrupts the G protein α-subunit binding and activating (GBA) motif of Drosophila Nucleobindin-1 (dmNUCB1). Such disruption does not manifest any obvious fitness-related, morphological/developmental, or behavioral abnormalities. A single copy of this mutation or the knockdown of dmnucb1 in restricted sets of cells variously rescues pleiotropic mutant phenotypes arising from impaired inositol 1,4,5-trisphosphate receptor (IP3R) activity (in turn depleting cytoplasmic Ca²⁺ levels across diverse tissue types). Additionally, altered dmNUCB1 expression or function considerably reverses lifespan and mobility improvements effected by IP3R mutants, in a Drosophila model of amyotrophic lateral sclerosis. Homology modeling-based analyses further predict a high degree of conformational conservation in Drosophila, of biochemically validated structural determinants in the GBA motif that specify in vertebrates, the unconventional Ca²⁺-regulated interaction of NUCB1 with Gαi subunits. The broad implications of our findings are hypothetically discussed, regarding potential roles for NUCB1 in GBA-mediated, Golgi-associated Ca²⁺ signaling, in health and disease.

tartan underlies the evolution of Drosophila male genital morphology

The morphology of male genitalia evolves rapidly, probably driven by sexual selection. However, little is known about the genes underlying genitalia differences between species. Identifying these genes is key to understanding how sexual selection acts to produce rapid phenotypic change. We have found that the gene tartan underlies differences between male Drosophila mauritiana and Drosophila simulans in the size and bristle number of the claspers—genital projections that grasp the female during copulation. Moreover, since tartan encodes a protein that is involved in cell interactions, this may represent an alternative developmental mechanism for morphological change. Therefore, our study provides insights into the genetic and developmental bases for the rapid evolution of male genitalia and organ size more generally.

Male genital structures are among the most rapidly evolving morphological traits and are often the only features that can distinguish closely related species. This process is thought to be driven by sexual selection and may reinforce species separation. However, while the genetic bases of many phenotypic differences have been identified, we still lack knowledge about the genes underlying evolutionary differences in male genital organs and organ size more generally. The claspers (surstyli) are periphallic structures that play an important role in copulation in insects. Here, we show that divergence in clasper size and bristle number between Drosophila mauritiana and Drosophila simulans is caused by evolutionary changes in tartan (trn), which encodes a transmembrane leucine-rich repeat domain protein that mediates cell–cell interactions and affinity. There are no fixed amino acid differences in trn between D. mauritiana and D. simulans, but differences in the expression of this gene in developing genitalia suggest that cis-regulatory changes in trn underlie the evolution of clasper morphology in these species. Finally, analyses of reciprocal hemizygotes that are genetically identical, except for the species from which the functional allele of trn originates, determined that the trn allele of D. mauritiana specifies larger claspers with more bristles than the allele of D. simulans. Therefore, we have identified a gene underlying evolutionary change in the size of a male genital organ, which will help to better understand not only the rapid diversification of these structures, but also the regulation and evolution of organ size more broadly.

Rapid genome-wide evolution in Brassica rapa populations following drought revealed by sequencing of ancestral and descendant gene pools

There is increasing evidence that evolution can occur rapidly in response to selection. Recent advances in sequencing suggest the possibility of documenting genetic changes as they occur in populations, thus uncovering the genetic basis of evolution, particularly if samples are available from both before and after selection. Here, we had a unique opportunity to directly assess genetic changes in natural populations following an evolutionary response to a fluctuation in climate. We analysed genome‐wide differences between ancestors and descendants of natural populations of Brassica rapa plants from two locations that rapidly evolved changes in multiple phenotypic traits, including flowering time, following a multiyear late‐season drought in California. These ancestor‐descendant comparisons revealed evolutionary shifts in allele frequencies in many genes. Some genes showing evolutionary shifts have functions related to drought stress and flowering time, consistent with an adaptive response to selection. Loci differentiated between ancestors and descendants (F_ST outliers) were generally different from those showing signatures of selection based on site frequency spectrum analysis (Tajima's D), indicating that the loci that evolved in response to the recent drought and those under historical selection were generally distinct. Very few genes showed similar evolutionary responses between two geographically distinct populations, suggesting independent genetic trajectories of evolution yielding parallel phenotypic changes. The results show that selection can result in rapid genome‐wide evolutionary shifts in allele frequencies in natural populations, and highlight the usefulness of combining resurrection experiments in natural populations with genomics for studying the genetic basis of adaptive evolution.

see also the Perspective by Hancock

Whole-genome sequencing of two North American Drosophila melanogaster populations reveals genetic differentiation and positive selection

The prevailing demographic model for Drosophila melanogaster suggests that the colonization of North America occurred very recently from a subset of European flies that rapidly expanded across the continent. This model implies a sudden population growth and range expansion consistent with very low or no population subdivision. As flies adapt to new environments, local adaptation events may be expected. To describe demographic and selective events during North American colonization, we have generated a data set of 35 individual whole-genome sequences from inbred lines of D. melanogaster from a west coast US population (Winters, California, USA) and compared them with a public genome data set from Raleigh (Raleigh, North Carolina, USA). We analysed nuclear and mitochondrial genomes and described levels of variation and divergence within and between these two North American D. melanogaster populations. Both populations exhibit negative values of Tajima's D across the genome, a common signature of demographic expansion. We also detected a low but significant level of genome-wide differentiation between the two populations, as well as multiple allele surfing events, which can be the result of gene drift in local subpopulations on the edge of an expansion wave. In contrast to this genome-wide pattern, we uncovered a 50-kilobase segment in chromosome arm 3L that showed all the hallmarks of a soft selective sweep in both populations. A comparison of allele frequencies within this divergent region among six populations from three continents allowed us to cluster these populations in two differentiated groups, providing evidence for the action of natural selection on a global scale.

Sequencing of pooled DNA samples (Pool-Seq) uncovers complex dynamics of transposable element insertions in Drosophila melanogaster

Transposable elements (TEs) are mobile genetic elements that parasitize genomes by semi-autonomously increasing their own copy number within the host genome. While TEs are important for genome evolution, appropriate methods for performing unbiased genome-wide surveys of TE variation in natural populations have been lacking. Here, we describe a novel and cost-effective approach for estimating population frequencies of TE insertions using paired-end Illumina reads from a pooled population sample. Importantly, the method treats insertions present in and absent from the reference genome identically, allowing unbiased TE population frequency estimates. We apply this method to data from a natural Drosophila melanogaster population from Portugal. Consistent with previous reports, we show that low recombining genomic regions harbor more TE insertions and maintain insertions at higher frequencies than do high recombining regions. We conservatively estimate that there are almost twice as many “novel” TE insertion sites as sites known from the reference sequence in our population sample (6,824 novel versus 3,639 reference sites, with on average a 31-fold coverage per insertion site). Different families of transposable elements show large differences in their insertion densities and population frequencies. Our analyses suggest that the history of TE activity significantly contributes to this pattern, with recently active families segregating at lower frequencies than those active in the more distant past. Finally, using our high-resolution TE abundance measurements, we identified 13 candidate positively selected TE insertions based on their high population frequencies and on low Tajima's D values in their neighborhoods.

Transposable elements (TE's) are parasitic genetic elements that spread by replicating themselves within a host genome. Most organisms are burdened with transposable elements; in fact, up to 80% of some genomes can consist of TE–derived DNA. Here, we use new sequencing technology to examine variation in genomic TE composition within a population at a finer scale and in a more unbiased fashion than has been possible before. We study a Portuguese population of D. melanogaster and find a large number of TE insertions, most of which occur in few individuals. Our analysis confirms that TE insertions are subject to purifying selection that counteracts their spread, and it suggests that the genome records waves of past TE invasions, with recently active elements occurring at low population frequency. We also find indications that TE insertions may sometimes have beneficial effects.

PoPoolation2: identifying differentiation between populations using sequencing of pooled DNA samples (Pool-Seq)

Summary: Sequencing pooled DNA samples (Pool-Seq) is the most cost-effective approach for the genome-wide comparison of population samples. Here, we introduce PoPoolation2, the first software tool specifically designed for the comparison of populations with Pool-Seq data. PoPoolation2 implements a range of commonly used measures of differentiation (F_ST, Fisher's exact test and Cochran-Mantel-Haenszel test) that can be applied on different scales (windows, genes, exons, SNPs). The result may be visualized with the widely used Integrated Genomics Viewer.

Availability and Implementation: PoPoolation2 is implemented in Perl and R. It is freely available on http://code.google.com/p/popoolation2/

Contact:christian.schloetterer@vetmeduni.ac.at

Supplementary Information: Manual: http://code.google.com/p/popoolation2/wiki/Manual Test data and tutorial: http://code.google.com/p/popoolation2/wiki/Tutorial Validation: http://code.google.com/p/popoolation2/wiki/Validation

Genome-wide genetic marker discovery and genotyping using next-generation sequencing

The advent of next-generation sequencing (NGS) has revolutionized genomic and transcriptomic approaches to biology. These new sequencing tools are also valuable for the discovery, validation and assessment of genetic markers in populations. Here we review and discuss best practices for several NGS methods for genome-wide genetic marker development and genotyping that use restriction enzyme digestion of target genomes to reduce the complexity of the target. These new methods -- which include reduced-representation sequencing using reduced-representation libraries (RRLs) or complexity reduction of polymorphic sequences (CRoPS), restriction-site-associated DNA sequencing (RAD-seq) and low coverage genotyping -- are applicable to both model organisms with high-quality reference genome sequences and, excitingly, to non-model species with no existing genomic data.