Stephen W. Schaeffer

Chromosomal Rearrangement Inferred From Comparisons of 12 Drosophila Genomes

The availability of 12 complete genomes of various species of genus Drosophila provides a unique opportunity to analyze genome-scale chromosomal rearrangements among a group of closely related species. This article reports on the comparison of gene order between these 12 species and on the fixed rearrangement events that disrupt gene order. Three major themes are addressed: the conservation of syntenic blocks across species, the disruption of syntenic blocks (via chromosomal inversion events) and its relationship to the phylogenetic distribution of these species, and the rate of rearrangement events over evolutionary time. Comparison of syntenic blocks across this large genomic data set confirms that genetic elements are largely (95%) localized to the same Muller element across genus Drosophila species and paracentric inversions serve as the dominant mechanism for shuffling the order of genes along a chromosome. Gene-order scrambling between species is in accordance with the estimated evolutionary distances between them and we find it to approximate a linear process over time (linear to exponential with alternate divergence time estimates). We find the distribution of synteny segment sizes to be biased by a large number of small segments with comparatively fewer large segments. Our results provide estimated chromosomal evolution rates across this set of species on the basis of whole-genome synteny analysis, which are found to be higher than those previously reported. Identification of conserved syntenic blocks across these genomes suggests a large number of conserved blocks with varying levels of embryonic expression correlation in Drosophila melanogaster. On the other hand, an analysis of the disruption of syntenic blocks between species allowed the identification of fixed inversion breakpoints and estimates of breakpoint reuse and lineage-specific breakpoint event segregation.

Polytene Chromosomal Maps of 11 Drosophila Species: The Order of Genomic Scaffolds Inferred From Genetic and Physical Maps

The sequencing of the 12 genomes of members of the genus Drosophila was taken as an opportunity to reevaluate the genetic and physical maps for 11 of the species, in part to aid in the mapping of assembled scaffolds. Here, we present an overview of the importance of cytogenetic maps to Drosophila biology and to the concepts of chromosomal evolution. Physical and genetic markers were used to anchor the genome assembly scaffolds to the polytene chromosomal maps for each species. In addition, a computational approach was used to anchor smaller scaffolds on the basis of the analysis of syntenic blocks. We present the chromosomal map data from each of the 11 sequenced non-Drosophila melanogaster species as a series of sections. Each section reviews the history of the polytene chromosome maps for each species, presents the new polytene chromosome maps, and anchors the genomic scaffolds to the cytological maps using genetic and physical markers. The mapping data agree with Mullers idea that the majority of Drosophila genes are syntenic. Despite the conservation of genes within homologous chromosome arms across species, the karyotypes of these species have changed through the fusion of chromosomal arms followed by subsequent rearrangement events.

Evolutionary genomics of inversions in Drosophila pseudoobscura: Evidence for epistasis

Drosophila pseudoobscura harbors a rich polymorphism for paracentric inversions on the third chromosome, and the clines in the inversion frequencies across the southwestern United States indicate that strong natural selection operates on them. Isogenic inversion strains were made from isofemale lines collected from four localities, and eight molecular markers were mapped on the third chromosome. Nucleotide diversity was measured for these loci and formed the basis of an evolutionary genomic analysis. The loci were differentiated among inversions. The inversions did not show significant differences among populations, however, likely the result of extensive gene flow among populations. Some loci had significant reductions in nucleotide diversity within inversions compared with interspecies divergence, suggesting that these loci are near inversion breakpoints or are near targets of directional selection. Linkage disequilibrium (LD) levels tended to decrease with distance between loci, indicating that some genetic exchange occurs among gene arrangements despite the presence of inversions. In some cases, however, adjacent genes had low levels of interlocus LD and loosely linked genes had high levels of interlocus LD, suggesting strong epistatic selection. Our results support the hypothesis that the inversions of D. pseudoobscura have emerged as suppressors of recombination to maintain positive epistatic relationships among loci within gene arrangements that developed as the species adapted to a heterogeneous environment.

Divergence Between the Drosophila pseudoobscura and D. persimilis Genome Sequences in Relation to Chromosomal Inversions

As whole-genome sequence assemblies accumulate, a challenge is to determine how these can be used to address fundamental evolutionary questions, such as inferring the process of speciation. Here, we use the sequence assemblies of Drosophila pseudoobscura and D. persimilis to test hypotheses regarding divergence with gene flow. We observe low differentiation between the two genome sequences in pericentromeric and peritelomeric regions. We interpret this result as primarily a remnant of the correlation between levels of variation and local recombination rate observed within populations. However, we also observe lower differentiation far from the fixed chromosomal inversions distinguishing these species and greater differentiation within and near these inversions. This finding is consistent with models suggesting that chromosomal inversions facilitate species divergence despite interspecies gene flow. We also document heterogeneity among the inverted regions in their degree of differentiation, suggesting temporal differences in the origin of each inverted region consistent with the inversions arising during a process of divergence with gene flow. While this study provides insights into the speciation process using two single-genome sequences, it was informed by lower throughput but more rigorous examinations of polymorphism and divergence. This reliance highlights the need for complementary genomic and population genetic approaches for tackling fundamental evolutionary questions such as speciation.

S -allele sequence diversity in natural populations of Solanum carolinense (Horsenettle)

S-allele diversity in Solanum carolinense was surveyed in two natural populations, located in Tennessee and North Carolina, with a molecular assay to determine the genotype of individual plants. A total of 13 different S-alleles were identified and sequenced. There is high overlap between the two populations sampled, with 10 alleles shared in common, one allele found only in Tennessee, and two found only in North Carolina. The number of alleles in this species appears to be extremely low compared with other species with gametophytic self-incompatibility. Sequence comparisons show that most alleles are extremely different one from another in their primary sequence and a phylogenetic analysis indicates extensive trans-specific evolution of S-lineages. In addition, some alleles appear to be derived much more recently. The implications of these observations are discussed in the light of recent theoretical results on S-allele population diversity and persistence.

Tubeworm succession at hydrothermal vents: use of biogenic cues to reduce habitat selection error?

Abstract Species colonizing new deep-sea hydrothermal vents along the East Pacific Rise show a distinct successional sequence: pioneer assemblages dominated by the vestimentiferan tubeworm Tevnia jerichonana being subsequently invaded by another vestimentiferan Riftia pachyptila, and eventually the mussel Bathymodiolus thermophilus. Using a manipulative approach modified from shallow-water ecological studies, we test three alternative hypotheses to explain the initial colonization by T. jerichonana and its subsequent replacement by R. pachyptila. We show that R. pachyptila and another vestimentiferan, Oasisia alvinae, colonized new surfaces only if the surfaces also were colonized by T. jerichonana. This pattern does not appear to be due to restricted habitat tolerances or inferior dispersal capabilities of R. pachyptila and O. alvinae, and we argue the alternative explanation that T. jerichonana facilitates the settlement of the other two species and is eventually outcompeted by R. pachyptila. Unlike the classic model of community succession, in which facilitating species promote their own demise by modifying the environment to make it more hospitable for competitors, we suggest that T. jerichonana may produce a chemical substance that induces settlement of these competitors. This process of selecting habitat based on biogenic cues may be especially adaptive and widespread among later-successional species that occupy a physically variable and unpredictable environment. In these cases, the presence of weedy species implies some integrated period of environmental suitability, whereas an instantaneous assessment of physical habitat conditions, such as water temperature for vent tubeworms, provides a poorer predictor of long-term habitat suitability.

Molecular population genetics of sequence length diversity in the Adh region of Drosophila pseudoobscura

Positive and negative selection on indel variation may explain the correlation between intron length and recombination levels in natural populations of Drosophila. A nucleotide sequence analysis of the 3.5 kilobase sequence of the alcohol dehydrogenase (Adh) region from 139 Drosophila pseudoobscura strains and one D. miranda strain was used to determine whether positive or negative selection acts on indel variation in a gene that experiences high levels of recombination. A total of 30 deletion and 36 insertion polymorphisms were segregating within D. pseudoobscura populations and no indels were fixed between D. pseudoobscura and its two sibling species D. miranda and D. persimilis. The ratio of Tajimas D to its theoretical minimum value (D(min)) was proposed as a metric to assess the heterogeneity in D among D. pseudoobscura loci when the number of segregating sites differs among loci. The magnitude of the D/D(min) ratio was found to increase as the rate of population expansion increases, allowing one to assess which loci have an excess of rare variants due to population expansion versus purifying selection. D. pseudoobscura populations appear to have had modest increases in size accounting for some of the observed excess of rare variants. The D/D(min) ratio rejected a neutral model for deletion polymorphisms. Linkage disequilibrium among pairs of indels was greater than between pairs of segregating nucleotides. These results suggest that purifying selection removes deletion variation from intron sequences, but not insertion polymorphisms. Genome rearrangement and size-dependent intron evolution are proposed as mechanisms that limit runaway intron expansion.

Mechanisms of Genetic Exchange within the Chromosomal Inversions of Drosophila Pseudoobscura

We have used the inversion system of Drosophila pseudoobscura to investigate how genetic flux occurs among the gene arrangements. The patterns of nucleotide polymorphism at seven loci were used to infer gene conversion events between pairs of different gene arrangements. We estimate that the average gene conversion tract length is 205 bp and that the average conversion rate is 3.4 × 10−6, which is 2 orders of magnitude greater than the mutation rate. We did not detect gene conversion events between all combinations of gene arrangements even though there was sufficient nucleotide variation for detection and sufficient opportunity for exchanges to occur. Genetic flux across the inverted chromosome resulted in higher levels of differentiation within 0.1 Mb of inversion breakpoints, but a slightly lower level of differentiation in central inverted regions. No gene conversion events were detected within 17 kb of an inversion breakpoint suggesting that the formation of double-strand breaks is reduced near rearrangement breakpoints in heterozygotes. At least one case where selection rather than proximity to an inversion breakpoint is responsible for reduction in polymorphism was identified.

Selection in Heterogeneous Environments Maintains the Gene Arrangement Polymorphism of Drosophila pseudoobscura

Abstract Chromosomal rearrangements may play an important role in how populations adapt to a local environment. The gene arrangement polymorphism on the third chromosome of Drosophila pseudoobscura is a model system to help determine the role that inversions play in the evolution of this species. The gene arrangements are the likely target of strong selection because they form classical clines across diverse geographic habitats, they cycle in frequency over seasons, and they form stable equilibria in population cages. A numerical approach was developed to estimate the fitness sets for 15 gene arrangement karyotypes in six niches based on a model of selection–migration balance. Gene arrangement frequencies in the six different niches were able to reach a stable meta-population equilibrium that matched the observed gene arrangement frequencies when recursions used the estimated fitnesses with a variety of initial inversion frequencies. These analyses show that a complex pattern of selection is operating in the six niches to maintain the D. pseudoobscura gene arrangement polymorphism. Models of local adaptation predict that the new inversion mutations were able to invade populations because they held combinations of two to 13 local adaptation loci together.

Restriction to large‐scale gene flow vs. regional panmixia among cold seep Escarpia spp. (Polychaeta, Siboglinidae)

The history of colonization and dispersal in fauna distributed among deep‐sea chemosynthetic ecosystems remains enigmatic and poorly understood because of an inability to mark and track individuals. A combination of molecular, morphological and environmental data improves understanding of spatial and temporal scales at which panmixia, disruption of gene flow or even speciation may occur. Vestimentiferan tubeworms of the genus Escarpia are important components of deep ‐sea cold seep ecosystems, as they provide long‐term habitat for many other taxa. Three species of Escarpia, Escarpia spicata [Gulf of California (GoC)], Escarpia laminata [Gulf of Mexico (GoM)] and Escarpia southwardae (West African Cold Seeps), have been described based on morphology, but are not discriminated through the use of mitochondrial markers (cytochrome oxidase subunit 1; large ribosomal subunit rDNA, 16S; cytochrome b). Here, we also sequenced the exon‐primed intron‐crossing Haemoglobin subunit B2 intron and genotyped 28 microsatellites to (i) determine the level of genetic differentiation, if any, among the three geographically separated entities and (ii) identify possible population structure at the regional scale within the GoM and West Africa. Results at the global scale support the occurrence of three genetically distinct groups. At the regional scale among eight sampling sites of E. laminata (n = 129) and among three sampling sites of E. southwardae (n = 80), no population structure was detected. These findings suggest that despite the patchiness and isolation of seep habitats, connectivity is high on regional scales.