Anthony J. Geneva

Genome sequencing reveals complex speciation in the Drosophila simulans clade

The three species of the Drosophila simulans clade--the cosmopolitan species, D. simulans, and the two island endemic species, D. mauritiana and D. sechellia--are important models in speciation genetics, but some details of their phylogenetic and speciation history remain unresolved. The order and timing of speciation are disputed, and the existence, magnitude, and timing of gene flow among the three species remain unclear. Here we report on the analysis of a whole-genome four-species sequence alignment that includes all three D. simulans clade species as well as the D. melanogaster reference sequence. The alignment comprises novel, paired short-read sequence data from a single highly inbred line each from D. simulans, D. mauritiana, and D. sechellia. We are unable to reject a species phylogeny with a basal polytomy; the estimated age of the polytomy is 242,000 yr before the present. However, we also find that up to 4.6% of autosomal and 2.2% of X-linked regions have evolutionary histories consistent with recent gene flow between the mainland species (D. simulans) and the two island endemic species (D. mauritiana and D. sechellia). Our findings thus show that gene flow has occurred throughout the genomes of the D. simulans clade species despite considerable geographic, ecological, and intrinsic reproductive isolation. Last, our analysis of lineage-specific changes confirms that the D. sechellia genome has experienced a significant excess of slightly deleterious changes and a dearth of presumed favorable changes. The relatively reduced efficacy of natural selection in D. sechellia is consistent with its derived, persistently reduced historical effective population size.

ANOLIS SEX CHROMOSOMES ARE DERIVED FROM A SINGLE ANCESTRAL PAIR

To explain the frequency and distribution of heteromorphic sex chromosomes in the lizard genus Anolis, we compared the relative roles of sex chromosome conservation versus turnover of sex‐determining mechanisms. We used model‐based comparative methods to reconstruct karyotype evolution and the presence of heteromorphic sex chromosomes onto a newly generated Anolis phylogeny. We found that heteromorphic sex chromosomes evolved multiple times in the genus. Fluorescent in situ hybridization (FISH) of repetitive DNA showed variable rates of Y chromosome degeneration among Anolis species and identified previously undetected, homomorphic sex chromosomes in two species. We confirmed homology of sex chromosomes in the genus by performing FISH of an X‐linked bacterial artificial chromosome (BAC) and quantitative PCR of X‐linked genes in multiple Anolis species sampled across the phylogeny. Taken together, these results are consistent with long‐term conservation of sex chromosomes in the group. Our results pave the way to address additional questions related to Anolis sex chromosome evolution and describe a conceptual framework that can be used to evaluate the origins and evolution of heteromorphic sex chromosomes in other clades.

Molecular phylogenetic analysis of tropical freshwater mussels (Mollusca: Bivalvia: Unionoida) resolves the position of Coelatura and supports a monophyletic Unionidae

In previous molecular phylogenetic analyses of the freshwater mussel family Unionidae (Bivalvia: Unionoida), the Afrotropical genus Coelatura had been recovered in various positions, generally indicating a paraphyletic Unionidae. However that result was typically poorly supported and in conflict with morphology-based analyses. We set out to test the phylogenetic position of Coelatura by sampling tropical lineages omitted from previous studies. Forty-one partial 28S nuclear rDNA and partial COI mtDNA sequences (1130 total aligned nucleotides) were analyzed separately and in combination under both maximum parsimony and likelihood, as well as Bayesian inference. There was significant phylogenetic incongruence between the character sets (partition homogeneity test, p<0.01), but a novel heuristic for comparing bootstrap values among character sets analyzed separately and in combination illustrated that the observed conflict was due to homoplasy rather than separate gene histories. Phylogenetic analyses robustly supported a monophyletic Unionidae, with Coelatura recovered as part of a well-supported Africa-India clade (=Parreysiinae). The implications of this result are discussed in the context of Afrotropical freshwater mussel evolution and the classification of the family Unionidae.

Genome Diversity and Divergence in Drosophila mauritiana: Multiple Signatures of Faster X Evolution

Drosophila mauritiana is an Indian Ocean island endemic species that diverged from its two sister species, Drosophila simulans and Drosophila sechellia, approximately 240,000 years ago. Multiple forms of incomplete reproductive isolation have evolved among these species, including sexual, gametic, ecological, and intrinsic postzygotic barriers, with crosses among all three species conforming to Haldane’s rule: F1 hybrid males are sterile and F1 hybrid females are fertile. Extensive genetic resources and the fertility of hybrid females have made D. mauritiana, in particular, an important model for speciation genetics. Analyses between D. mauritiana and both of its siblings have shown that the X chromosome makes a disproportionate contribution to hybrid male sterility. But why the X plays a special role in the evolution of hybrid sterility in these, and other, species remains an unsolved problem. To complement functional genetic analyses, we have investigated the population genomics of D. mauritiana, giving special attention to differences between the X and the autosomes. We present a de novo genome assembly of D. mauritiana annotated with RNAseq data and a whole-genome analysis of polymorphism and divergence from ten individuals. Our analyses show that, relative to the autosomes, the X chromosome has reduced nucleotide diversity but elevated nucleotide divergence; an excess of recurrent adaptive evolution at its protein-coding genes; an excess of recent, strong selective sweeps; and a large excess of satellite DNA. Interestingly, one of two centimorgan-scale selective sweeps on the D. mauritiana X chromosome spans a region containing two sex-ratio meiotic drive elements and a high concentration of satellite DNA. Furthermore, genes with roles in reproduction and chromosome biology are enriched among genes that have histories of recurrent adaptive protein evolution. Together, these genome-wide analyses suggest that genetic conflict and frequent positive natural selection on the X chromosome have shaped the molecular evolutionary history of D. mauritiana, refining our understanding of the possible causes of the large X-effect in speciation.

RWTY (R We There Yet): An R Package for Examining Convergence of Bayesian Phylogenetic Analyses

Bayesian inference using Markov chain Monte Carlo (MCMC) has become one of the primary methods used to infer phylogenies from sequence data. Assessing convergence is a crucial component of these analyses, as it establishes the reliability of the posterior distribution estimates of the tree topology and model parameters sampled from the MCMC. Numerous tests and visualizations have been developed for this purpose, but many of the most popular methods are implemented in ways that make them inconvenient to use for large data sets. RWTY is an R package that implements established and new methods for diagnosing phylogenetic MCMC convergence in a single convenient interface.

Molecular phylogenetic analysis supports a Gondwanan origin of the Hyriidae (Mollusca: Bivalvia: Unionida) and the paraphyly of Australasian taxa

The freshwater mussel family Hyriidae (Mollusca: Bivalvia: Unionida) has a disjunct trans-Pacific distribution in Australasia and South America. Previous phylogenetic analyses have estimated the evolutionary relationships of the family and the major infra-familial taxa (Velesunioninae and Hyriinae: Hyridellini in Australia; Hyriinae: Hyriini, Castaliini, and Rhipidodontini in South America), but taxon and character sampling have been too incomplete to support a predictive classification or allow testing of biogeographical hypotheses. We sampled 30 freshwater mussel individuals representing the aforementioned hyriid taxa, as well as outgroup species representing the five other freshwater mussel families and their marine sister group (order Trigoniida). Our ingroup included representatives of all Australian genera. Phylogenetic relationships were estimated from three gene fragments (nuclear 28S, COI and 16S mtDNA) using maximum parsimony, maximum likelihood, and Bayesian inference, and we applied a Bayesian relaxed clock model calibrated with fossil dates to estimate node ages. Our analyses found good support for monophyly of the Hyriidae and the subfamilies and tribes, as well as the paraphyly of the Australasian taxa (Velesunioninae, (Hyridellini, (Rhipidodontini, (Castaliini, Hyriini)))). The Hyriidae was recovered as sister to a clade comprised of all other Recent freshwater mussel families. Our molecular date estimation supported Cretaceous origins of the major hyriid clades, pre-dating the Tertiary isolation of South America from Antarctica/Australia. We hypothesize that early diversification of the Hyriidae was driven by terrestrial barriers on Gondwana rather than marine barriers following disintegration of the super-continent.

Using genomic data to revisit an early example of reproductive character displacement in Haitian Anolis lizards

The pattern of reproductive character displacement (RCD)—in which traits associated with reproductive isolation are more different where two species occur together than where they occur in isolation—is frequently attributed to reinforcement, a process during which natural selection acting against maladaptive mating events leads to enhanced prezygotic isolation between species or incipient species. One of the first studies of RCD to include molecular genetic data was described 40 years ago in a complex of Haitian trunk anole lizards using a small number of allozyme loci. In this example, Anolis caudalis appears to experience divergence in the color and pattern of an extensible throat fan, or dewlap, in areas of contact with closely related species at the northern and southern limits of its range. However, this case study has been largely overlooked for decades; meanwhile, explanations for geographic variation in dewlap color and pattern have focused primarily on adaptation to local signalling environments. We reinvestigate this example using amplified fragment length polymorphism (AFLP) genome scans, mtDNA sequence data, information on dewlap phenotypes and GIS data on environmental variation to test the hypothesis of RCD generated by reinforcement in Haitian trunk anoles. Together, our phenotypic and genetic results are consistent with RCD at the southern and northern limits of the range of A. caudalis. We evaluate the evidence for reinforcement as the explanation for RCD in Haitian trunk anoles, consider alternative explanations and provide suggestions for future work on the relationship between dewlap variation and speciation in Haitian trunk anoles.

Molecular phylogeny of thorny catfishes (Siluriformes: Doradidae).

Doradidae is a monophyletic catfish family endemic to continental South America, and composed of 93 valid species here placed in 31 genera. Existing phylogenetic hypotheses for Doradidae are derived from comprehensive analyses of morphological data, and a single molecular-based study on a limited subset of taxa. To provide a robust molecular phylogeny commensurate with those based on morphology, we gathered original and published sequence data for 86 species-level taxa (at least 70 valid species plus 16 new or questionably nominal species) and all genera of Doradidae, as well as 10 species (nine genera) of Auchenipteridae and three species and genera of Aspredinidae as outgroups. 3011 base pairs were aligned for two mitochondrial genes (cytochrome c oxidase subunit 1, and 16S ribosomal RNA) and one nuclear gene (recombination activating gene 1), and analyzed for a total of 143 specimens (130 doradids, 10 auchenipterids and three aspredinids). Tree topologies generated by Maximum Parsimony, Maximum Likelihood, and Bayesian analyses were largely congruent, and are compared to existing phylogenies based on morphology and molecules. Although many of the relationships supported by our molecular analyses corroborated those based on morphology, others are newly hypothesized or remain in conflict. The monotypic Wertheimeria, Franciscodoras and Kalyptodoras, for example, form a newly proposed clade, and the subfamily Astrodoradinae is placed at the base of the doradid tree. The monotypic Doraops and Centrochir, endemic to Caribbean drainages north and west of the Andes, are sister to Pterodoras and Platydoras, respectively, two genera that are widely distributed in Atlantic drainages. Additional biogeographic implications are discussed for hypothesized relationships among doradids. Molecular evidence strongly supports synonymization of monotypic Merodoras with Amblydoras, and transfer of Amblydoras bolivarensis to genus Scorpiodoras. Furthermore, we consider Opsodoras ternetzi to be more properly placed in the genus Nemadoras. The genus Opsodoras may warrant synonymization with Hemidoras, and the monophyly of genus Ossancora is not supported; however, we refrain from taxonomic decisions regarding those taxa until a broader spectrum of doradids can be submitted to further morphological and molecular phylogenetic analyses.

A new method to scan genomes for introgression in a secondary contact model.

Secondary contact between divergent populations or incipient species may result in the exchange and introgression of genomic material. We develop a simple DNA sequence measure, called G min, which is designed to identify genomic regions experiencing introgression in a secondary contact model. G min is defined as the ratio of the minimum between-population number of nucleotide differences in a genomic window to the average number of between-population differences. Although it is conceptually simple, one advantage of G min is that it is computationally inexpensive relative to model-based methods for detecting gene flow and it scales easily to the level of whole-genome analysis. We compare the sensitivity and specificity of G min to those of the widely used index of population differentiation, F ST, and suggest a simple statistical test for identifying genomic outliers. Extensive computer simulations demonstrate that G min has both greater sensitivity and specificity for detecting recent introgression than does F ST. Furthermore, we find that the sensitivity of G min is robust with respect to both the population mutation and recombination rates. Finally, a scan of G min across the X chromosome of Drosophila melanogaster identifies candidate regions of introgression between sub-Saharan African and cosmopolitan populations that were previously missed by other methods. These results show that G min is a biologically straightforward, yet powerful, alternative to F ST, as well as to more computationally intensive model-based methods for detecting gene flow.

Multilocus phylogenetic analyses of Hispaniolan and Bahamian trunk anoles (distichus species group)

The distichus species group includes six species and 21 subspecies of trunk ecomorph anoles distributed across Hispaniola and its satellite islands as well as the northern Bahamas. Although this group has long served as a model system for studies of reproductive character displacement, adaptation, behavior and speciation, it has never been the subject of a comprehensive phylogenetic analysis. Our goal here is to generate a multilocus phylogenetic dataset (one mitochondrial and seven nuclear loci) and to use this dataset to infer phylogenetic relationships among the majority of the taxa assigned to the distichus species group. We use these phylogenetic trees to address three topics about the groups evolution. First, we consider longstanding taxonomic controversies about the status of several species and subspecies assigned to the distichus species group. Second, we investigate the biogeographic history of the group and specifically test the hypotheses that historical division of Hispaniola into two paleo-islands contributed to the groups diversification and that Bahamian and Hispaniolan satellite island populations are derived from colonists from the main Hispaniolan landmass. Finally, third, we use comparative phylogenetic analyses to test the hypothesis that divergence between pale yellow and darkly pigmented orange or red dewlap coloration has occurred repeatedly across the distichus species group.