Rob J. Kulathinal

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.

The genomics of speciation in Drosophila: diversity, divergence, and introgression estimated using low-coverage genome sequencing.

In nature, closely related species may hybridize while still retaining their distinctive identities. Chromosomal regions that experience reduced recombination in hybrids, such as within inversions, have been hypothesized to contribute to the maintenance of species integrity. Here, we examine genomic sequences from closely related fruit fly taxa of the Drosophila pseudoobscura subgroup to reconstruct their evolutionary histories and past patterns of genic exchange. Partial genomic assemblies were generated from two subspecies of Drosophila pseudoobscura (D. ps.) and an outgroup species, D. miranda. These new assemblies were compared to available assemblies of D. ps. pseudoobscura and D. persimilis, two species with overlapping ranges in western North America. Within inverted regions, nucleotide divergence among each pair of the three species is comparable, whereas divergence between D. ps. pseudoobscura and D. persimilis in non-inverted regions is much lower and closer to levels of intraspecific variation. Using molecular markers flanking each of the major chromosomal inversions, we identify strong crossover suppression in F1 hybrids extending over 2 megabase pairs (Mbp) beyond the inversion breakpoints. These regions of crossover suppression also exhibit the high nucleotide divergence associated with inverted regions. Finally, by comparison to a geographically isolated subspecies, D. ps. bogotana, our results suggest that autosomal gene exchange between the North American species, D. ps. pseudoobscura and D. persimilis, occurred since the split of the subspecies, likely within the last 200,000 years. We conclude that chromosomal rearrangements have been vital to the ongoing persistence of these species despite recent hybridization. Our study serves as a proof-of-principle on how whole genome sequencing can be applied to formulate and test hypotheses about species formation in lesser-known non-model systems.

Bayesian Analysis Suggests that Most Amino Acid Replacements in Drosophila Are Driven by Positive Selection

One of the principal goals of population genetics is to understand the processes by which genetic variation within species (polymorphism) becomes converted into genetic differences between species (divergence). In this transformation, selective neutrality, near neutrality, and positive selection may each play a role, differing from one gene to the next. Synonymous nucleotide sites are often used as a uniform standard of comparison across genes on the grounds that synonymous sites are subject to relatively weak selective constraints and so may, to a first approximation, be regarded as neutral. Synonymous sites are also interdigitated with nonsynonymous sites and so are affected equally by genomic context and demographic factors. Hence a comparison of levels of polymorphism and divergence between synonymous sites and amino acid replacement sites in a gene is potentially informative about the magnitude of selective forces associated with amino acid replacements. We have analyzed 56 genes in which polymorphism data from D. simulans are compared with divergence from a reference strain of D. melanogaster. The framework of the analysis is Bayesian and assumes that the distribution of selective effects (Malthusian fitnesses) is Gaussian with a mean that differs for each gene. In such a model, the average scaled selection intensity (γ =Nes) of amino acid replacements eligible to become polymorphic or fixed is −7.31, and the standard deviation of selective effects within each locus is 6.79 (assuming homoscedasticity across loci). For newly arising mutations of this type that occur in autosomal or X-linked genes, the average proportion of beneficial mutations is 19.7%. Among the amino acid polymorphisms in the sample, the expected average proportion of beneficial mutations is 47.7%, and among amino acid replacements that become fixed the average proportion of beneficial mutations is 94.3%. The average scaled selection intensity of fixed mutations is +5.1. The presence of positive selection is pervasive with the single exception of kl-5, a Y-linked fertility gene. We find no evidence that a significant fraction of fixed amino acid replacements is neutral or nearly neutral or that positive selection drives amino acid replacements at only a subset of the loci. These results are model dependent and we discuss possible modifications of the model that might allow more neutral and nearly neutral amino acid replacements to be fixed.

Fine-scale mapping of recombination rate in Drosophila refines its correlation to diversity and divergence

Regional rates of recombination often correlate with levels of nucleotide diversity, and either selective or neutral hypotheses can explain this relationship. Regional recombination rates also correlate with nucleotide differences between human and chimpanzee, consistent with models where recombination is mutagenic; however, a lack of correlation is observed in the Drosophila melanogaster group, consistent with models invoking natural selection. Here, we revisit the relationship among recombination, diversity, and interspecies difference by generating empirical estimates of these parameters in Drosophila pseudoobscura. To measure recombination rate, we genotyped 1,294 backcross hybrids at 50 markers across the largest assembled linkage group in this species. Genome-wide diversity was estimated by sequencing a second isolate of D. pseudoobscura at shallow coverage. Alignment to the sequenced genome of the closely related species, Drosophila persimilis, provided nucleotide site orthology. Our findings demonstrate that scale is critical in determining correlates to recombination rate: fine-scale cross-over rate estimates are far stronger predictors of both diversity and interspecies difference than broad-scale estimates. The correlation of fine-scale recombination rate to diversity and interspecies difference appears to be genome-wide, evidenced by examination of an X-linked region in greater detail. Because we observe a strong correlation of cross-over rate with interspecies difference, even after correcting for segregating ancestral variation, we suggest that both mutagenic and selective forces generate these correlations, the latter in regions of low crossing over. We propose that it is not cross-overs per se that are mutagenic, but rather repair of DNA double-strand break precursors via crossing over and gene conversion.

Population Genomic Inferences from Sparse High-Throughput Sequencing of Two Populations of Drosophila melanogaster

Short-read sequencing techniques provide the opportunity to capture genome-wide sequence data in a single experiment. A current challenge is to identify questions that shallow-depth genomic data can address successfully and to develop corresponding analytical methods that are statistically sound. Here, we apply the Roche/454 platform to survey natural variation in strains of Drosophila melanogaster from an African (n = 3) and a North American (n = 6) population. Reads were aligned to the reference D. melanogaster genomic assembly, single nucleotide polymorphisms were identified, and nucleotide variation was quantified genome wide. Simulations and empirical results suggest that nucleotide diversity can be accurately estimated from sparse data with as little as 0.2× coverage per line. The unbiased genomic sampling provided by random short-read sequencing also allows insight into distributions of transposable elements and copy number polymorphisms found within populations and demonstrates that short-read sequencing methods provide an efficient means to quantify variation in genome organization and content. Continued development of methods for statistical inference of shallow-depth genome-wide sequencing data will allow such sparse, partial data sets to become the norm in the emerging field of population genomics.

CYTOLOGICAL CHARACTERIZATION OF PREMEIOTIC VERSUS POSTMEIOTIC DEFECTS PRODUCING HYBRID MALE STERILITY AMONG SIBLING SPECIES OF THE DROSOPHILA MELANOGASTER COMPLEX

In accordance with Haldanes rule, hybridizations between species of the Drosophila simulans clade produce fertile females but sterile males. In this study, a comprehensive characterization was undertaken on the six types of F1 males that were the result of the crosses between D. simulans, D. sechellia, and D. mauritiana. With the use of light and electron microscopy, it was shown that while each particular hybrid genotype exhibited a specific sterility phenotype, these phenotypes fell into two distinct classes. The two hybrid genotypes that possessed D. mauritiana X‐chromosomes contained spermatogenic defects that caused arrests in premeiotic spermatogenic stages. The other four F1 hybrids possessed postmeiotic spermatogenic defects. Nonsynchronous cell divisions, underdeveloped mitochondrial derivative‐axonemal associations, and microtubule abnormalities were common to all of these hybrids. Each particular postmeiotically defective hybrid genotype demonstrated characteristically distinct profiles in sperm bundle number in addition to characteristic spermiogenic arrests in the furthest developed spermatids. These results in species hybrids contrast with the absence of significant differences in spermatogenic characters between species of this clade. In addition, by utilizing an attached‐X cross, we investigated the influence of maternal effects and cytoplasmic factors on the sterility of D. simulans F1 hybrids and found none. However, we discovered a strain of D. simulans (2119) that caused a large shift in sterility from postmeiotic to premeiotic when crossed to D. sechellia. This suggests that D. simulans is polymorphic for genes involving premeiotic and postmeiotic sterility and that the two types of sterilities between species may have a simple genetic basis.

The molecular basis of speciation: from patterns to processes, rules to mechanisms

The empirical study of speciation has brought us closer to unlocking the origins of life’s vast diversity. By examining recently formed species, a number of general patterns, or rules, become apparent. Among fixed differences between species, sexual genes and traits are one of the most rapidly evolving and novel functional classes, and premating isolation often develops earlier than postmating isolation. Among interspecific hybrids, sterility evolves faster than inviability, the X-chromosome has a greater effect on incompatibilities than autosomes, and hybrid dysfunction affects the heterogametic sex more frequently than the homogametic sex (Haldane’s rule). Haldane’s rule, in particular, has played a major role in reviving interest in the genetics of speciation. However, the large genetic and reproductive differences between taxa and the multi-factorial nature of each rule have made it difficult to ascribe general mechanisms. Here, we review the extensive progress made since Darwin on understanding the origin of species. We revisit the rules of speciation, regarding them as landmarks as species evolve through time. We contrast these ‘rules’ of speciation to ‘mechanisms’ of speciation representing primary causal factors ranging across various levels of organization—from genic to chromosomal to organismal. To explain the rules, we propose a new ‘hierarchical faster-sex’ theory: the rapid evolution of sex and reproduction-related (SRR) genes (faster-SRR evolution), in combination with the preferential involvement of the X-chromosome (hemizygous X-effects) and sexually selected male traits (faster-male evolution). This unified theory explains a comprehensive set of speciation rules at both the prezyotic and postzygotic levels and also serves as a cohesive alternative to dominance, composite, and recent genomic conflict interpretations of Haldane’s rule.

Elevated Evolutionary Rates among Functionally Diverged Reproductive Genes across Deep Vertebrate Lineages

Among closely related taxa, proteins involved in reproduction generally evolve more rapidly than other proteins. Here, we apply a functional and comparative genomics approach to compare functional divergence across a deep phylogenetic array of egg-laying and live-bearing vertebrate taxa. We aligned and annotated a set of 4,986 1 : 1 : 1 : 1 : 1 orthologs in Anolis carolinensis (green lizard), Danio rerio (zebrafish), Xenopus tropicalis (frog), Gallus gallus (chicken), and Mus musculus (mouse) according to function using ESTs from available reproductive (including testis and ovary) and non-reproductive tissues as well as Gene Ontology. For each species lineage, genes were further classified as tissue-specific (found in a single tissue) or tissue-expressed (found in multiple tissues). Within independent vertebrate lineages, we generally find that gonadal-specific genes evolve at a faster rate than gonadal-expressed genes and significantly faster than non-reproductive genes. Among the gonadal set, testis genes are generally more diverged than ovary genes. Surprisingly, an opposite but nonsignificant pattern is found among the subset of orthologs that remained functionally conserved across all five lineages. These contrasting evolutionary patterns found between functionally diverged and functionally conserved reproductive orthologs provide evidence for pervasive and potentially cryptic lineage-specific selective processes on ancestral reproductive systems in vertebrates.

A BIOGEOGRAPHIC GENETIC APPROACH FOR TESTING THE ROLE OF REINFORCEMENT: THE CASE OF DROSOPHILA PSEUDOOBSCURA AND D. PERSIMILIS

Abstract.— The role of reinforcement in speciation can be explained by two distinct models. In model I, two diverged populations hybridize and produce fertile hybrids that successfully backcross (hybridization with gene flow). In model II, two populations hybridize but succeeding backcrosses are unproductive (hybridization without gene flow). Using Drosophila persimilis and D. pseudoobscura, we have tested model I by comparing the extent of heterospecific introgression in sympatric versus allopatric populations. We show that certain expectations of this particular model of reinforcement, which is based on hybridization and gene flow between divergent populations after secondary contact, are not realized in these two species. The evidence consists of the similarity of genetic distances as well as proportions of unique/rare alleles between sympatric and allopatric heterospecific populations and a negative correlation between genetic distance and geographical distance between heterospecific populations, which suggests ecological differentiation. This approach in quantifying differential gene flow has important consequences to studies that compare sympatric and allopatric isolation using genetic distance. Following model I, one would expect a pattern of higher prezygotic isolation in sympatric species compared to allopatric species of the same genetic distance simply as a result of an underestimation of genetic distance due to introgression between sympatric populations. We suggest more parsimonious explanations such as reinforcement without genetic exchange (model II) and ecological differentiation, which require high levels of preexisting reproductive isolation between populations.

The Nature of Genetic Variation in Sex and Reproduction-related Genes Among Sibling Species of the Drosophila melanogaster Complex

Much is known about the biology of Drosophila melanogaster. As a model organism, a comprehensive understanding of its development, physiology and reproduction has been acquired. As a result, a broad variety of transferable genetic tools and information has allowed sibling species of the D. melanogaster complex to emerge as an important speciation model system. By comparing D. melanogaster with its close relative, Drosophila simulans, as well as its other sibling species, we are beginning to understand the nature of genetic changes during the early stages of speciation. In general, we find that genes and traits involved in sex and reproduction are more variable. A large assortment of genes and traits that are involved in various aspects of mating and fertility reveal diagnostic differences between these sibling species. Sex and reproduction-related (SRR) genes are, on average, more diverged than genes with no apparent reproductive function. Furthermore, SRR genes appear more permissive at opting in novel function. These results follow a general trend observed in other taxa and demonstrate the preferential involvement of SRR genes in reproductive isolation and species formation.