CAGEE: Computational Analysis of Gene Expression Evolution.

Despite the increasing abundance of whole transcriptome data, few methods are available to analyze global gene expression across phylogenies. Here, we present a new software package (Computational Analysis of Gene Expression Evolution [CAGEE]) for inferring patterns of increases and decreases in gene expression across a phylogenetic tree, as well as the rate at which these changes occur. In contrast to previous methods that treat each gene independently, CAGEE can calculate genome-wide rates of gene expression, along with ancestral states for each gene. The statistical approach developed here makes it possible to infer lineage-specific shifts in rates of evolution across the genome, in addition to possible differences in rates among multiple tissues sampled from the same species. We demonstrate the accuracy and robustness of our method on simulated data and apply it to a data set of ovule gene expression collected from multiple self-compatible and self-incompatible species in the genus Solanum to test hypotheses about the evolutionary forces acting during mating system shifts. These comparisons allow us to highlight the power of CAGEE, demonstrating its utility for use in any empirical system and for the analysis of most morphological traits. Our software is available at https://github.com/hahnlab/CAGEE/.

Inferring the Genetic Basis of Sex Determination from the Genome of a Dioecious Nightshade.

Dissecting the genetic mechanisms underlying dioecy (i.e., separate female and male individuals) is critical for understanding the evolution of this pervasive reproductive strategy. Nonetheless, the genetic basis of sex determination remains unclear in many cases, especially in systems where dioecy has arisen recently. Within the economically important plant genus Solanum (∼2,000 species), dioecy is thought to have evolved independently at least 4 times across roughly 20 species. Here, we generate the first genome sequence of a dioecious Solanum and use it to ascertain the genetic basis of sex determination in this species. We de novo assembled and annotated the genome of Solanum appendiculatum (assembly size: ∼750 Mb scaffold N50: 0.92 Mb; ∼35,000 genes), identified sex-specific sequences and their locations in the genome, and inferred that males in this species are the heterogametic sex. We also analyzed gene expression patterns in floral tissues of males and females, finding approximately 100 genes that are differentially expressed between the sexes. These analyses, together with observed patterns of gene-family evolution specific to S. appendiculatum, consistently implicate a suite of genes from the regulatory network controlling pectin degradation and modification in the expression of sex. Furthermore, the genome of a species with a relatively young sex-determination system provides the foundational resources for future studies on the independent evolution of dioecy in this clade.

Hybrid incompatibility between Drosophila virilis and D. lummei is stronger in the presence of transposable elements.

Mismatches between parental genomes in selfish elements are frequently hypothesized to underlie hybrid dysfunction and drive speciation. However, because the genetic basis of most hybrid incompatibilities is unknown, testing the contribution of selfish elements to reproductive isolation is difficult. Here, we evaluated the role of transposable elements (TEs) in hybrid incompatibilities between Drosophila virilis and D. lummei by experimentally comparing hybrid incompatibility in a cross where active TEs are present in D. virilis (TE+) and absent in D. lummei, to a cross where these TEs are absent from both D. virilis (TE-) and D. lummei genotypes. Using genomic data, we confirmed copy number differences in TEs between the D. virilis (TE+) strain and both the D. virilis (TE-) strain and D. lummei. We observed F1 postzygotic reproductive isolation exclusively in the interspecific cross involving TE+ D. virilis but not in crosses involving TE- D. virilis. This mirrors intraspecies dysgenesis where atrophied testes only occur when TE+ D. virilis is the paternal parent. A series of backcross experiments, that accounted for alternative models of hybrid incompatibility, showed that both F1 hybrid incompatibility and intrastrain dysgenesis are consistent with the action of TEs rather than genic interactions. Thus, our data suggest that this TE mechanism manifests as two different incompatibility phenotypes. A further Y-autosome interaction contributes to additional, sex-specific, inviability in one direction of this cross-combination. These experiments demonstrate that TEs that cause intraspecies dysgenesis can increase reproductive isolation between closely related lineages, thereby adding to the processes that consolidate speciation.

Testing potential mechanisms of conspecific sperm precedence in Drosophila pseudoobscura.

Drosophila pseudoobscura females that co-occur with sister species D. persimilis show elevated fertilization by conspecific sperm when they mate with both a heterospecific and a conspecific male. This phenomenon, known as conspecific sperm precedence (CSP), has evolved as a mechanism to avoid maladaptive hybridization with D. persimilis. In this study, we assessed pericopulatory (during mating) and postcopulatory (after mating) traits in crosses with sympatric or allopatric D. pseudoobscura females and conspecific or heterospecific males to evaluate potential mechanisms of CSP in this system. We observed shorter copulation duration in crosses with sympatric females, but found no difference in quantity of sperm transferred or female reproductive tract toxicity between sympatry and allopatry. Our data show some support for the hypothesis that parasperm, a short, sterile sperm morph, can protect fertile eusperm from the D. pseudoobscura female reproductive tract, though it is unclear how this might affect patterns of sperm use in sympatry vs. allopatry. Overall, these results suggest that copulation duration could potentially contribute to the elevated CSP observed in sympatry.

Regional differences in the abiotic environment contribute to genomic divergence within a wild tomato species.

The wild currant tomato Solanum pimpinellifolium inhabits a wide range of abiotic habitats across its native range of Ecuador and Peru. Although it has served as a key genetic resource for the improvement of domestic cultivars, little is known about the genetic basis of traits underlying local adaptation in this species, nor what abiotic variables are most important for driving differentiation. Here we use redundancy analysis (RDA) and other multivariate statistical methods (structural equation modelling [SEM] and generalized dissimilarity modelling [GDM]) to quantify the relationship of genomic variation (6,830 single nucleotide polymorphisms [SNPs]) with climate and geography, among 140 wild accessions. RDA, SEM and GDM each identified environment as explaining more genomic variation than geography, suggesting that local adaptation to heterogeneous abiotic habitats may be an important source of genetic diversity in this species. Environmental factors describing temporal variation in precipitation and evaporative demand explained the most SNP variation among accessions, indicating that these forces may represent key selective agents. Lastly, by studying how SNP-environment associations vary throughout the genome (44,064 SNPs), we mapped the location and investigated the functions of loci putatively contributing to climatic adaptations. Together, our findings indicate an important role for selection imposed by the abiotic environment in driving genomic differentiation between populations.

Intraspecific Genetic Variation Underlying Postmating Reproductive Barriers between Species in the Wild Tomato Clade (Solanum sect. Lycopersicon).

A goal of speciation genetics is to understand how the genetic components underlying interspecific reproductive barriers originate within species. Unilateral incompatibility (UI) is a postmating prezygotic barrier in which pollen rejection in the female reproductive tract (style) occurs in only one direction of an interspecific cross. Natural variation in the strength of UI has been observed among populations within species in the wild tomato clade. In some cases, molecular loci underlying self-incompatibility (SI) are associated with this variation in UI, but the mechanistic connection between these intra- and inter-specific pollen rejection behaviors is poorly understood in most instances. We generated an F2 population between SI and SC genotypes of a single species, Solanum pennellii, to examine the genetic basis of intraspecific variation in UI against other species, and to determine whether loci underlying SI are genetically associated with this variation. We found that F2 individuals vary in the rate at which UI rejection occurs. One large effect QTL detected for this trait co-localized with the SI-determining S-locus. Moreover, individuals that expressed S-RNase-the S-locus protein involved in SI pollen rejection-in their styles had much more rapid UI responses compared with those without S-RNase protein. Our analysis shows that intraspecific variation at mate choice loci-in this case at loci that prevent self-fertilization-can contribute to variation in the expression of interspecific isolation, including postmating prezygotic barriers. Understanding the nature of such intraspecific variation can provide insight into the accumulation of these barriers between diverging lineages.

Pervasive antagonistic interactions among hybrid incompatibility loci.

Species barriers, expressed as hybrid inviability and sterility, are often due to epistatic interactions between divergent loci from two lineages. Theoretical models indicate that the strength, direction, and complexity of these genetic interactions can strongly affect the expression of interspecific reproductive isolation and the rates at which new species evolve. Nonetheless, empirical analyses have not quantified the frequency with which loci are involved in interactions affecting hybrid fitness, and whether these loci predominantly interact synergistically or antagonistically, or preferentially involve loci that have strong individual effects on hybrid fitness. We systematically examined the prevalence of interactions between pairs of short chromosomal regions from one species (Solanum habrochaites) co-introgressed into a heterospecific genetic background (Solanum lycopersicum), using lines containing pairwise combinations of 15 chromosomal segments from S. habrochaites in the background of S. lycopersicum (i.e., 95 double introgression lines). We compared the strength of hybrid incompatibility (either pollen sterility or seed sterility) expressed in each double introgression line to the expected additive effect of its two component single introgressions. We found that epistasis was common among co-introgressed regions. Interactions for hybrid dysfunction were substantially more prevalent in pollen fertility compared to seed fertility phenotypes, and were overwhelmingly antagonistic (i.e., double hybrids were less unfit than expected from additive single introgression effects). This pervasive antagonism is expected to attenuate the rate at which hybrid infertility accumulates among lineages over time (i.e., giving diminishing returns as more reproductive isolation loci accumulate), as well as decouple patterns of accumulation of sterility loci and hybrid incompatibility phenotypes. This decoupling effect might explain observed differences between pollen and seed fertility in their fit to theoretical predictions of the accumulation of isolation loci, including the 'snowball' effect.

Desiccation resistance and pigmentation variation reflects bioclimatic differences in the Drosophila americana species complex.

BACKGROUND: Disentangling the selective factors shaping adaptive trait variation is an important but challenging task. Many studies-especially in Drosophila-have documented trait variation along latitudinal or altitudinal clines, but frequently lack resolution about specific environmental gradients that could be causal selective agents, and often do not investigate covariation between traits simultaneously. Here we examined variation in multiple macroecological factors across geographic space and their associations with variation in three physiological traits (desiccation resistance, UV resistance, and pigmentation) at both population and species scales, to address the role of abiotic environment in shaping trait variation. RESULTS: Using environmental data from collection locations of three North American Drosophila species-D. americana americana, D. americana texana and D. novamexicana-we identified two primary axes of macroecological variation; these differentiated species habitats and were strongly loaded for precipitation and moisture variables. In nine focal populations (three per species) assayed for each trait, we detected significant species-level variation for both desiccation resistance and pigmentation, but not for UV resistance. Species-level trait variation was consistent with differential natural selection imposed by variation in habitat water availability, although patterns of variation differed between desiccation resistance and pigmentation, and we found little evidence for pleiotropy between traits. CONCLUSIONS: Our multi-faceted approach enabled us to identify potential agents of natural selection and examine how they might influence the evolution of multiple traits at different evolutionary scales. Our findings highlight that environmental factors influence functional trait variation in ways that can be complex, and point to the importance of studies that examine these relationships at both population- and species-levels.

Conspecific sperm precedence is reinforced, but postcopulatory sexual selection weakened, in sympatric populations of Drosophila.

Sexual selection can accelerate speciation by driving the evolution of reproductive isolation, but forces driving speciation could also reciprocally feedback on sexual selection. This might be particularly important in the context of 'reinforcement', where selection acts directly to increase prezygotic barriers to reduce the cost of heterospecific matings. Using assays of sperm competition within and between two sister species, we show a signature of reinforcement where these species interact: populations of Drosophila pseudoobscura that co-occur with sister species D. persimilis have an elevated ability to outcompete heterospecific sperm, consistent with selection for increased postcopulatory isolation. We also find these D. pseudoobscura populations have decreased sperm competitive ability against conspecifics, reducing the opportunity for sexual selection within these populations. Our findings demonstrate that direct selection to increase reproductive isolation against other species can compromise the efficacy of sexual selection within species, a collateral effect of reproductive traits responding to heterospecific interactions.

A simple genetic architecture and low constraint allow rapid floral evolution in a diverse and recently radiating plant genus.

Genetic correlations among different components of phenotypes, especially those resulting from pleiotropy, can constrain or facilitate trait evolution. These factors could especially influence the evolution of traits that are functionally integrated, such as those comprising the flower. Indeed, pleiotropy is proposed as a main driver of repeated convergent trait transitions, including the evolution of phenotypically similar pollinator syndromes. We assessed the role of pleiotropy in the differentiation of floral and other reproductive traits between two species - Jaltomata sinuosa and J. umbellata (Solanaceae) - that have divergent suites of floral traits consistent with bee and hummingbird pollination, respectively. To do so, we generated a hybrid population and examined the genetic architecture (trait segregation and quantitative trait locus (QTL) distribution) underlying 25 floral and fertility traits. We found that most floral traits had a relatively simple genetic basis (few, predominantly additive, QTLs of moderate to large effect), as well as little evidence of antagonistic pleiotropy (few trait correlations and QTL colocalization, particularly between traits of different classes). However, we did detect a potential case of adaptive pleiotropy among floral size and nectar traits. These mechanisms may have facilitated the rapid floral trait evolution observed within Jaltomata, and may be a common component of rapid phenotypic change more broadly.

Phylogenomics Reveals Three Sources of Adaptive Variation during a Rapid Radiation.

Speciation events often occur in rapid bursts of diversification, but the ecological and genetic factors that promote these radiations are still much debated. Using whole transcriptomes from all 13 species in the ecologically and reproductively diverse wild tomato clade (Solanum sect. Lycopersicon), we infer the species phylogeny and patterns of genetic diversity in this group. Despite widespread phylogenetic discordance due to the sorting of ancestral variation, we date the origin of this radiation to approximately 2.5 million years ago and find evidence for at least three sources of adaptive genetic variation that fuel diversification. First, we detect introgression both historically between early-branching lineages and recently between individual populations, at specific loci whose functions indicate likely adaptive benefits. Second, we find evidence of lineage-specific de novo evolution for many genes, including loci involved in the production of red fruit color. Finally, using a "PhyloGWAS" approach, we detect environment-specific sorting of ancestral variation among populations that come from different species but share common environmental conditions. Estimated across the whole clade, small but substantial and approximately equal fractions of the euchromatic portion of the genome are inferred to contribute to each of these three sources of adaptive genetic variation. These results indicate that multiple genetic sources can promote rapid diversification and speciation in response to new ecological opportunity, in agreement with our emerging phylogenomic understanding of the complexity of both ancient and recent species radiations.

Dissecting the basis of novel trait evolution in a radiation with widespread phylogenetic discordance.

Phylogenetic analyses of trait evolution can provide insight into the evolutionary processes that initiate and drive phenotypic diversification. However, recent phylogenomic studies have revealed extensive gene tree-species tree discordance, which can lead to incorrect inferences of trait evolution if only a single species tree is used for analysis. This phenomenon-dubbed "hemiplasy"-is particularly important to consider during analyses of character evolution in rapidly radiating groups, where discordance is widespread. Here, we generate whole-transcriptome data for a phylogenetic analysis of 14 species in the plant genus Jaltomata (the sister clade to Solanum), which has experienced rapid, recent trait evolution, including in fruit and nectar colour, and flower size and shape. Consistent with other radiations, we find evidence for rampant gene tree discordance due to incomplete lineage sorting (ILS) and to introgression events among the well-supported subclades. As both ILS and introgression increase the probability of hemiplasy, we perform several analyses that take discordance into account while identifying genes that might contribute to phenotypic evolution. Despite discordance, the history of fruit colour evolution in Jaltomata can be inferred with high confidence, and we find evidence of de novo adaptive evolution at individual genes associated with fruit colour variation. In contrast, hemiplasy appears to strongly affect inferences about floral character transitions in Jaltomata, and we identify candidate loci that could arise either from multiple lineage-specific substitutions or standing ancestral polymorphisms. Our analysis provides a generalizable example of how to manage discordance when identifying loci associated with trait evolution in a radiating lineage.

Interspecific tests of allelism reveal the evolutionary timing and pattern of accumulation of reproductive isolation mutations.

Despite extensive theory, little is known about the empirical accumulation and evolutionary timing of mutations that contribute to speciation. Here we combined QTL (Quantitative Trait Loci) analyses of reproductive isolation, with information on species evolutionary relationships, to reconstruct the order and timing of mutations contributing to reproductive isolation between three plant (Solanum) species. To evaluate whether reproductive isolation QTL that appear to coincide in more than one species pair are homologous, we used cross-specific tests of allelism and found evidence for both homologous and lineage-specific (non-homologous) alleles at these co-localized loci. These data, along with isolation QTL unique to single species pairs, indicate that >85% of isolation-causing mutations arose later in the history of divergence between species. Phylogenetically explicit analyses of these data support non-linear models of accumulation of hybrid incompatibility, although the specific best-fit model differs between seed (pairwise interactions) and pollen (multi-locus interactions) sterility traits. Our findings corroborate theory that predicts an acceleration ('snowballing') in the accumulation of isolation loci as lineages progressively diverge, and suggest different underlying genetic bases for pollen versus seed sterility. Pollen sterility in particular appears to be due to complex genetic interactions, and we show this is consistent with a snowball model where later arising mutations are more likely to be involved in pairwise or multi-locus interactions that specifically involve ancestral alleles, compared to earlier arising mutations.

Molecular mechanisms of postmating prezygotic reproductive isolation uncovered by transcriptome analysis.

Little is known about the physiological responses and genetic mutations associated with reproductive isolation between species, especially for postmating prezygotic isolating barriers. Here, we examine changes in gene expression that accompany the expression of 'unilateral incompatibility' (UI)-a postmating prezygotic barrier in which fertilization is prevented by gamete rejection in the reproductive tract [in this case of pollen tubes (male gametophytes)] in one direction of a species cross, but is successful in the reciprocal crossing direction. We use whole-transcriptome sequencing of multiple developmental stages of male and female tissues in two Solanum species that exhibit UI to: (i) identify transcript differences between UI-competent and UI noncompetent tissues; (ii) characterize transcriptional changes specifically associated with the phenotypic expression of UI; and (iii) using these comparisons, evaluate the behaviour of a priori candidate loci for UI and identify new candidates for future manipulative work. In addition to describing transcriptome-wide changes in gene expression that accompany this isolating barrier, we identify at least five strong candidates for involvement in postmating prezygotic incompatibility between species. These include three novel candidates and two candidates that are strongly supported by prior developmental, functional, and quantitative trait locus mapping studies. These latter genes are known molecular players in the intraspecific expression of mate choice via genetic self-incompatibility, and our study supports prior evidence that these inter- and intraspecific postmating prezygotic reproductive behaviours share specific genetic and molecular mechanisms.

Intraspecific sperm competition genes enforce post-mating species barriers in Drosophila.

Sexual selection and sexual conflict are considered important drivers of speciation, based on both theoretical models and empirical correlations between sexually selected traits and diversification. However, whether reproductive isolation between species evolves directly as a consequence of intrapopulation sexual dynamics remains empirically unresolved, in part because knowledge of the genetic mechanisms (if any) connecting these processes is limited. Here, we provide evidence of a direct mechanistic link between intraspecies sexual selection and reproductive isolation. We examined genes with known roles in intraspecific sperm competition (ISC) in D. melanogaster and assayed their impact on conspecific sperm precedence (CSP). We found that two such genes (Acp36DE and CG9997) contribute to both offensive sperm competition and CSP; null/knockdown lines both had lower competitive ability against D. melanogaster conspecifics and were no longer able to displace heterospecific D. simulans sperm in competitive matings. In comparison, Sex Peptide (Acp70A)-another locus essential for ISC-does not contribute to CSP. These data indicate that two loci important for sperm competitive interactions have an additional role in similar interactions that enforce post-mating reproductive isolation between species, and show that sexual selection and sexual isolation can act on the same molecular targets in a gene-specific manner.

Morphological and anatomical determinants of mesophyll conductance in wild relatives of tomato (Solanum sect. Lycopersicon, sect. Lycopersicoides; Solanaceae).

Natural selection on photosynthetic performance is a primary factor determining leaf phenotypes. The complex CO2 diffusion path from substomatal cavities to the chloroplasts - the mesophyll conductance (g(m)) - limits photosynthetic rate in many species and hence shapes variation in leaf morphology and anatomy. Among sclerophyllous and succulent taxa, structural investment in leaves, measured as the leaf dry mass per area (LMA), has been implicated in decreased gm . However, in herbaceous taxa with high g(m), it is less certain how LMA impacts CO2 diffusion and whether it significantly affects photosynthetic performance. We addressed these questions in the context of understanding the ecophysiological significance of leaf trait variation in wild tomatoes, a closely related group of herbaceous perennials. Although g(m) was high in wild tomatoes, variation in g(m) significantly affected photosynthesis. Even in these tender-leaved herbaceous species, greater LMA led to reduced g(m). This relationship between g(m) and LMA is partially mediated by cell packing and leaf thickness, although amphistomy (equal distribution of stomata on both sides of the leaf) mitigates the effect of leaf thickness. Understanding the costs of increased LMA will inform future work on the adaptive significance of leaf trait variation across ecological gradients in wild tomatoes and other systems.

No evidence for phylogenetic constraint on natural defense evolution among wild tomatoes.

Plant defense traits can be shaped by evolutionary and physiological constraints, as well as local ecological selection. We assessed the relative importance of these factors in shaping defense trait variation across the wild tomato clade (a group of 13 closely related species) using an herbivore bioassay (Manduca sexta). With phylogenetic comparative methods, we evaluated patterns of constitutive and induced defense variation, and the extent of coupling between alternative defense strategies. We detected substantial variation among species and found no evidence for phylogenetic conservatism among defensive traits, unlike for two other ecologically relevant (reproductive) traits. In addition, constitutive and induced defense syndromes were unassociated. These data indicate that, in this group, there is no evidence for either phylogenetic conservatism of shared consumer guilds that shape defense traits, or for constraints on defense trait evolution, including mechanistic trade-offs between defense strategies. Our data suggest that defense trait variation in this clade instead results from rapid responses to local ecological conditions.

Merging ecology and genomics to dissect diversity in wild tomatoes and their relatives.

To understand the origin, history, and function, of natural biological variation, from nucleotide to community levels, is a fundamental promise of ecological genomics. The most fruitful systems for this work are those that possess both ecological and genomic resources. Such systems provide an opportunity to precisely dissect genetic and developmental mechanisms, and to connect genotypes to phenotypes, as well as to directly demonstrate the ecological and evolutionary relevance of this phenotypic variation. Here we synthesize findings emerging from our efforts to understand two fundamental evolutionary processes - speciation and adaptation - using ecological genomics approaches. Many of these studies have been in the wild tomato clade (Solanum section Lycopersicon), a group that has both exceptional diversity and genomic tools. We also highlight the expanding taxonomic reach of this work, especially in two genera - Capsicum and Jaltomata - that are closely related to Solanum. Parallel approaches in these ecologically and reproductively diverse clades enable us to examine novel questions and traits that are not captured within Solanum, while leveraging the power of comparative studies to understand shared ecological and evolutionary patterns. By synthesizing findings from phenotypic, ecophysiological, genetic, and comparative perspectives, our ultimate goal is to understand the complex mechanistic and evolutionary contributions to the formation of new traits and species diversity.

Long-sperm precedence and other cryptic female choices in Drosophila melanogaster.

Females that mate multiply make postmating choices about which sperm fertilize their eggs (cryptic female choice); however, the male characteristics they use to make such choices remain unclear. In this study, we sought to understand female sperm use patterns by evaluating whether Drosophila melanogaster females adjust sperm use (second male paternity) in response to four main factors: male genotype, male courtship effort, male pheromone alteration, and male postmating reproductive morphology. Our experiment was replicated across four different D. melanogaster lines, in a full factorial design, including a pheromone manipulation in which second males were perfumed to resemble heterospecific (D. yakuba) males. We found that females prefer longer sperm-regardless of mating order-in almost all contexts; this observed pattern of 'long-sperm precedence' is consistent with female postcopulatory choice of high-fitness male traits. Nonetheless, we also found that this general preference can be plastically altered by females in response to effects including perfuming treatment; this differential female sperm use is between otherwise identical males, and therefore solely female-mediated. Furthermore, our finding that females exercise choice using diverse criteria suggests a possible mechanism for the maintenance of variation in sexually selected male traits.