Detecting adaptive trait introgression between Iris fulva and I. brevicaulis in highly selective field conditions

The idea that natural hybridization has served as an important force in evolutionary and adaptive diversification has gained considerable momentum in recent years. By combining genome analyses with a highly selective field experiment, we provide evidence for adaptive trait introgression between two naturally hybridizing Louisiana Iris species, flood-tolerant Iris fulva and dry-adapted I. brevicaulis. We planted reciprocal backcross (BC1) hybrids along with pure-species plants into natural settings that, due to a flooding event, favored I. fulva. As expected, I. fulva plants survived at much higher rates than I. brevicaulis plants. Backcross hybrids toward I. fulva (BCIF) also survived at significantly higher rates than the reciprocal backcross toward I. brevicaulis (BCIB). Survivorship of BCIB hybrids was strongly influenced by the presence of a number of introgressed I. fulva alleles located throughout the genome, while survivorship in the reciprocal BCIF hybrids was heavily influenced by two epistatically acting QTL of opposite effects. These results demonstrate the potential for adaptive trait introgression between these two species and may help to explain patterns of genetic variation observed in naturally occurring hybrid zones.     

The genetic architecture of reproductive isolation in Louisiana irises: pollination syndromes and pollinator preferences

In animal-pollinated plants, pollinator preferences for divergent floral forms can lead to partial reproductive isolation. We describe regions of plant genomes that affect pollinator preferences for two species of Louisiana Irises, Iris brevicaulis and Iris fulva, and their artificial hybrids. Iris brevicaulis and I. fulva possess bee and bird-pollination syndromes, respectively. Hummingbirds preferred I. fulva and under-visited both I. brevicaulis and backcrosses toward this species. Lepidopterans preferred I. fulva and backcrosses toward I. fulva, but also under-visited I. brevicaulis and I. brevicaulis backcrosses. Bumblebees preferred I. brevicaulis and F1 hybrids and rarely visited I. fulva. Although all three pollen vectors preferred one or the other species, these preferences did not prevent visitation to other hybrid/parental classes. Quantitative trait locus (QTL) mapping, in reciprocal BC1 mapping populations, defined the genetic architecture of loci that affected pollinator behavior. We detected six and nine QTLs that affected pollinator visitation rates in the BCIb and BCIf mapping populations, respectively, with as many as three QTLs detected for each trait. Overall, this study reflects the possible role of quantitative genetic factors in determining (1) reproductive isolation, (2) the pattern of pollinator-mediated genetic exchange, and thus (3) hybrid zone evolution.     

Loci affecting long-term hybrid survivorship in Louisiana irises: implications for reproductive isolation and introgression

Iris fulva and I. brevicaulis are long-lived plant species known to hybridize where they coexist in nature. Year-to-year survival contributes significantly to overall fitness for both species and their hybrid derivatives, and differences in hybrid survivability may have important consequences to interspecific gene flow in nature. We examined the genetic architecture of long-term survivorship of reciprocal backcross I. fulva x I. brevicaulis hybrids in a common-garden, greenhouse environment. Differences in mortality were found between the two backcross (BC1) hybrid classes, with hybrids crossed toward I. fulva (BCIF) revealing twice the mortality of those hybrids backcrossed toward I. brevicaulis (BCIB). Using genomic scans on two separate genetic linkage maps derived from the reciprocal hybrid populations, we found that hybrid survivorship is influenced by several genetic regions. Multiple interval mapping (MIM) revealed four quantitative trait loci (QTLs) in BCIF hybrids that were significantly associated with survivorship. Introgressed I. brevicaulis DNA increased survivorship at three of the four QTLs. For the fourth QTL, introgressed I. brevicaulis DNA was associated with decreased survivorship. No QTLs were detected in BCIB hybrids; however, single-marker analysis revealed five unlinked loci that were significantly associated with survivorship. At all five markers, survivorship was positively associated with introgressed I. fulva DNA. The present findings have important implications for the evolutionary dynamics of naturally occurring hybrid zones. Regions of the genome that increase survivorship when in a heterozygous (i.e., hybrid) state should have an increased likelihood of passing across species boundaries, whereas those that decrease survivorship will be less likely to introgress.     

THE GENETIC ARCHITECTURE OF REPRODUCTIVE ISOLATION IN LOUISIANA IRISES: HYBRID FITNESS IN NATURE

Negative epistasis in hybrid genomes commonly results in postzygotic isolation between divergent lineages. However, some genomic regions may be selectively neutral or adaptive in hybrids and thus may potentially cross species barriers. We examined postzygotic isolation between ecologically similar species of Louisiana Iris: Iris brevicaulis and I. fulva to determine the potential for adaptive introgression in nature. Line-cross analyses allowed us a general overview of the gene action responsible for fitness-related traits. We then used a QTL mapping approach to detect genomic regions responsible for variation in these traits. Although hybrid classes suffered reduced fitness for many traits, hybrid means were equivalent to at least one of the parental species in overall estimates of maternal and paternal fitness during the two years of the field study. The genetic architecture underlying the fitness-related traits varied across field site and year of the study, thus emphasizing the importance of the environment in determining the degree of postzygotic isolation and potential for introgression across natural hybrid zones.     

QTL mapping reveals the genetic architecture of loci affecting pre- and post-zygotic isolating barriers in Louisiana Iris

ABSTRACT: BACKGROUND: Hybridization among Louisiana Irises has been well established and the genetic architecture of reproductive isolation is known to affect the potential for and the directionality of introgression between taxa. Here we use co-dominant markers to identify regions where QTL are located both within and between backcross maps to compare the genetic architecture of reproductive isolation and fitness traits across treatments and years. RESULTS: QTL mapping was used to elucidate the genetic architecture of reproductive isolation between Iris fulva and Iris brevicaulis. Homologous co-dominant EST-SSR markers scored in two backcross populations between I. fulva and I. brevicaulis were used to generate genetic linkage maps. These were used as the framework for mapping QTL associated with variation in 11 phenotypic traits likely responsible for reproductive isolation and fitness. QTL were dispersed throughout the genome, with the exception of one region of a single linkage group (LG) where QTL for flowering time, sterility, and fruit production clustered. In most cases, homologous QTL were not identified in both backcross populations, however, homologous QTL for flowering time, number of growth points per rhizome, number of nodes per inflorescence, and number of flowers per node were identified on several linkage groups. CONCLUSION: Two different traits affecting reproductive isolation, flowering time and sterility, exhibit different genetic architectures, with numerous QTL across the Iris genome controlling flowering time and fewer, less distributed QTL affecting sterility. QTL for traits affecting fitness are largely distributed across the genome with occasional overlap, especially on LG 4, where several QTL increasing fitness and decreasing sterility cluster. Given the distribution and effect direction of QTL affecting reproductive isolation and fitness, we have predicted genomic regions where introgression may be more likely to occur (those regions associated with an increase in fitness and unlinked to loci controlling reproductive isolation) and those that are less likely to exhibit introgression (those regions linked to traits decreasing fitness and reproductive isolation).     

Nectar Production in Louisiana Iris Hybrids

Nectar is an important attractant for pollinators, and a plant's success in sexual reproduction can be influenced by the amount and concentration of nectar produced by its flowers. We studied nectar production over flower lifetime in Iris fulva, Iris brevicaulis, and four classes of hybrids-reciprocal F1's and backcrosses-between these species. Iris fulva produced less concentrated nectar than did I. brevicaulis, whereas I. brevicaulis flowers had a shorter life span. Hybrids were not intermediate, but they had the high nectar concentration of I. brevicaulis combined with the long life span of I. fulva flowers. Nectar production and concentration declined after the first day in all classes, but flowers continued to produce nectar until they were completely wilted. Backcrosses did not show a shift in mean or increased variation for the characters that distinguished the parental species; backcrosses toward I. fulva retained the high nectar concentration of I. brevicaulis, and backcrosses toward I. brevicaulis did not have a reduced flower life span. Overall, F1 hybrid flowers produced the highest amounts of nectar and nectar sugar over their life spans. These results, together with previously obtained data on pollinator choice in mixed arrays of the same flower classes, show that F1 hybrids between these species do not suffer from reduced attractiveness to pollinators. F1 individuals produced more nectar and nectar sugar than did their parents, and thus, they are possibly even more attractive to pollinators that forage for nectar.     

The genetic architecture of reproductive isolation in Louisiana irises: flowering phenology

Despite the potential importance of divergent reproductive phenologies as a barrier to gene flow, we know less about the genetics of this factor than we do about any other isolating barrier. Here, we report on the genetic architecture of divergent flowering phenologies that result in substantial reproductive isolation between the naturally hybridizing plant species Iris fulva and I. brevicaulis. I. fulva initiates and terminates flowering significantly earlier than I. brevicaulis. We examined line crosses of reciprocal F1 and backcross (BC1) hybrids and determined that flowering time was polygenic in nature. We further defined quantitative trait loci (QTL) that affect the initiation of flowering in each of these species. QTL analyses were performed separately for two different growing seasons in the greenhouse, as well as in two field plots where experimental plants were placed into nature. For BCIF hybrids (BC1 toward I. fulva), 14 of 17 detected QTL caused flowering to occur later in the season when I. brevicaulis alleles were present, while the remaining 3 caused flowering to occur earlier. In BCIB hybrids (BC1 toward I. brevicaulis), 11 of 15 detected QTL caused flowering to occur earlier in the season when introgressed I. fulva alleles were present, while the remaining 4 caused flowering to occur later. These ratios are consistent with expectations of selection (as opposed to drift) promoting flowering divergence in the evolutionary history of these species. Furthermore, epistatic interactions among the QTL also reflected the same trends, with the majority of epistatic effects causing later flowering than expected in BCIF hybrids and earlier flowering in BCIB hybrids. Overlapping QTL that influenced flowering time across all four habitat/treatment types were not detected, indicating that increasing the sample size of genotyped plants would likely increase the number of significant QTL found in this study.     

Parapatric genetic introgression and phenotypic assimilation: testing conditions for introgression between Hercules beetles (Dynastes,Dynastinae)

The prevalence and consequences of genetic introgression between species have been intensively debated. I used Hercules beetles as examples to test for conditions that may be associated with the occurrence of introgression. RADseq data were used to reconstruct the species tree and history of introgression between Hercules beetles. Image data from museum specimens were used to investigate the phenotypic similarity of two adaptive traits between species from two distinct climatic realms (Nearctic vs. Neotropical). Genetic introgression was identified between Hercules beetles living in geographic proximity (parapatric). Phylogenetic relatedness and phenotypic similarity did not predict nor preclude genetic introgression between species. Phenotypic assimilation in body coloration was evident between distantly related Hercules beetles codistributed in Central America, where directional introgression was also statistically supported from the putative donor to receiver lineages. The number of introgressed loci was significantly higher between species with than without phenotypic similarity. I discuss the implications of recent studies on adaptive genetic introgression by providing supporting evidence from the Hercules beetle system.     

Intimately linked or hardly speaking? The relationship between genotype and environmental gradients in a Louisiana Iris hybrid population

Several models of hybrid zone evolution predict the same spatial patterns of genotypic distribution whether or not structuring is due to environment-dependent or -independent selection. In this study, we tested for evidence of environment-dependent selection in an Iris fulva x Iris brevicaulis hybrid population by examining the distribution of genotypes in relation to environmental gradients. We selected 201 Louisiana Iris plants from within a known hybrid population (80 m x 80 m) and placed them in four different genotypic classes (I. fulva, I. fulva-like hybrid, I. brevicaulis-like hybrid and I. brevicaulis) based on seven species-specific random amplified polymorphic DNA (RAPD) markers and two chloroplast DNA haplotypes. Environmental variables were then measured. These variables included percentage cover by tree canopy, elevation from the high water mark, soil pH and percentage soil organic matter. Each variable was sampled for all 201 plants. Canonical discriminant analysis (CDA) was used to infer the environmental factors most strongly associated with the different genotypic groups. Slight differences in elevation (-0.5 m to +0.4 m) were important for distinguishing habitat distributions described by CDA, even though there were no statistical differences between mean elevations alone. I. brevicaulis occurred in a broad range of habitats, while I. fulva had a narrower distribution. Of all the possible combinations, I. fulva-like hybrids and I. brevicaulis-like hybrids occurred in the most distinct habitat types relative to one another. Each hybrid class was not significantly different from its closest parent with regard to habitat occupied, but was statistically unique from its more distant parental species. Within the hybrid genotypes, most, but not all, RAPD loci were individually correlated with environmental variables. This study suggests that, at a very fine spatial scale, environment-dependent selection contributed to the genetic structuring of this hybrid zone.     

Environment-dependent performance and fitness of Iris brevicaulis, I. fulva (Iridaceae), and hybrids

We tested the relative fitness of two Louisiana Iris species (Iris brevicaulis and I. fulva) and their first-generation backcross hybrids in three experimental watering treatments: dry, field capacity, and flooded. Leaf area expansion rate, gas exchange (A(max), g(s), c(i)), and biomass at final harvest were measured for each species and hybrid class in all three environmental treatments. Fitness (based on total biomass) of the four genotypic classes differed significantly with environment. All genotypic classes performed most poorly in the dry treatment. The fitness ranking of genotypic class also changed across environments (significant genotypic class by treatment interaction) with hybrid genotype fitness shifting relative to parental genotypes. Integrating over all treatments, backcrosses to I. fulva showed the lowest fitness, whereas backcrosses to I. brevicaulis outperformed I. fulva. The differences in fitness were apparently achieved by a combination of differences in photosynthesis and allocation. In this system, hybrids are not necessarily less fit than their parents, and the relationship between hybrid and parental fitness is influenced by environmental conditions, lending support to the Hybrid Novelty model of hybrid zone evolution.     

Conservation genetics of bush mango from central/west Africa: implications from random amplified polymorphic DNA analysis

Genetic variation was assessed in the two bush mango species, Irvingia gabonensis and I. wombolu, valuable multipurpose fruit trees from central and west Africa that are currently undergoing domestication. A total of 130 individuals sampled from Cameroon, Nigeria and Gabon were analysed using 74 random amplified polymorphic DNAs (RAPDs). Significant genetic integrity was found in the two morphologically similar species (among-species analysis of molecular variance [AMOVA] variance component 25.8%, P < 0.001), with no evidence of hybridization, even between individuals from areas of sympatry where hybridization was considered probable. Results suggest that large-scale transplantation of either species into new habitats will probably not lead to genetic introgression from or into the other species. Therefore, subsequent cultivation of the two species should not be hindered by this consideration, although further studies on the potential for hybridization/introgression between these species would be prudent. Significant genetic differentiation of both species (among-countries within species, nested AMOVA variance component 9.8%, P < 0.001) was observed over the sampled regions, and genetic similarity of samples decreased significantly with increasing geographical distance, according to number of alleles in common (NAC) analysis. 'Hot spots' of genetic diversity were found clustered in southern Nigeria and southern Cameroon for I. wombolu, and in southern Nigeria, southern Cameroon and central Gabon for I. gabonensis. The possible reasons for this distribution of genetic variation are discussed, but it may reflect evolutionary history, as these populations occur in areas of postulated Pleistocene refugia. The application of these results to domestication programmes and, in the light of extensive deforestation in the region, conservation approaches, is discussed.     

NATURAL HYBRIDIZATION BETWEEN IRIS FULVA AND IRIS HEXAGONA: PATTERN OF RIBOSOMAL DNA VARIATION

Genetic variation associated with the natural hybridization of Iris fulva and I. hexagona was investigated to test for the occurrence of introgression. These species have been viewed as a classic example of the process of introgressive hybridization (Anderson, 1949). However, more recent studies have concluded that there has not been an exchange of genetic material between I. fulva and I. hexagona (Randolph et al., 1967). Our analysis has involved the examination of both allopatric and parapatric populations of I. fulva and I. hexagona with reference to diagnostic ribosomal DNA markers. The pattern of variation in the parapatric population indicates the presence of the repeated backcrossing necessary to the process of introgressive hybridization. Indeed, in the region of parapatry, we suggest that localized introgression of ribosomal sequences has occurred into both I. fulva and I. hexagona. Significantly, we have also detected the presence of the diagnostic ribosomal markers from each species in allopatric populations of the alternate species. Our findings suggest that not only is introgressive hybridization presently occurring in parapatry between I. fulva and I. hexagona, but that past hybridization between these species has resulted in introgression into areas of allopatry.     

Genetic diversity enhanced by ancient introgression and secondary contact in East Pacific black mangroves

Regional distribution of genetic diversity in widespread species may be influenced by hybridization with locally restricted, closely related species. Previous studies have shown that Central American East Pacific populations of the wide-ranged Avicennia germinans , the black mangrove, harbour higher genetic diversity than the rest of its range. Genetic diversity in this region might be enhanced by introgression with the locally restricted Avicennia bicolor . We tested the hypotheses of ancient hybridization using phylogenetic analysis of the internal transcribed spacer region (ITS) of the nuclear ribosomal DNA and intergenic chloroplast DNA; we also tested for current hybridization by population level analysis of nuclear microsatellites. Our results unveiled ancient ITS introgression between a northern Pacific Central American A. germinans lineage and A. bicolor . However, microsatellite data revealed contemporary isolation between the two species. Polymorphic ITS sequences from Costa Rica and Panama are consistent with a zone of admixture between the introgressant ITS A. germinans lineage and a southern Central American lineage of A. germinans . Interspecific introgression influenced lineage diversity and divergence at the nuclear ribosomal DNA; intraspecific population differentiation and secondary contact are more likely to have enhanced regional genetic diversity in Pacific Central American populations of the widespread A. germinans .     

A method for rapid and simultaneous mapping of genetic Loci and introgression sizes in nematode species

Caenorhabditis briggsae is emerging as an attractive model organism not only in studying comparative biology against C. elegans, but also in developing novel experimentation avenues. In particular, recent identification of a new Caenorhabditis species, C. sp.9 with which it can mate and produce viable progeny provides an opportunity for studying the genetics of hybrid incompatibilities (HI) between the two. Mapping of a specific HI locus demands repeated backcrossing to get hold of the specific genomic region underlying an observed phenotype. To facilitate mapping of HI loci between C. briggsae and C. sp.9, an efficient mapping method and a genetic map ideally consisting of dominant markers are required for systematic introgression of genomic fragments between the two species. We developed a fast and cost-effective method for high throughput mapping of dominant loci with resolution up to 1 million bps in C. briggsae. The method takes advantage of the introgression between C. briggsae and C. sp.9 followed by PCR genotyping using C. briggsae specific primers. Importantly, the mapping results can not only serve as an effective way for estimating the chromosomal position of a genetic locus in C. briggsae, but also provides size information for the introgression fragment in an otherwise C. sp.9 background. In addition, it also helps generate introgression line as a side-product that is invaluable for the subsequent mapping of HI loci. The method will greatly facilitate the construction of a genetic map consisting of dominant markers and pave the way for systematic isolation of HI loci between C. briggsae and C. sp.9 which has so far not been attempted between nematode species. The method is designed for mapping of a dominant allele, but can be easily adapted for mapping of any other type of alleles in any other species if introgression between a sister species pair is feasible.     

The genetic basis of behavioral isolation between Drosophila mauritiana and D. sechellia

Understanding how species form is a fundamental question in evolutionary biology. Identifying the genetic bases of barriers that prevent gene flow between species provides insight into how speciation occurs. Here, I analyze a poorly understood reproductive isolating barrier, prezygotic reproductive isolation. I perform a genetic analysis of prezygotic isolation between two closely related species of Drosophila, D. mauritiana and D. sechellia. I first confirm the existence of strong behavioral isolation between D. mauritiana females and D. sechellia males. Next, I examine the genetic basis of behavioral isolation by (1) scanning an existing set of introgression lines for chromosomal regions that have a large effect on isolation; and (2) mapping quantitative trait loci (QTL) that underlie behavioral isolation via backcross analysis. In particular, I map QTL that determine whether a hybrid backcross female and a D. sechellia male will mate. I identify a single significant QTL, on the X chromosome, suggesting that few major-effect loci contribute to behavioral isolation between these species. In further work, I refine the map position of the QTL to a small region of the X chromosome.     

Integrating Bayesian genomic cline analyses and association mapping of morphological and ecological traits to dissect reproductive isolation and introgression in a Louisiana Iris hybrid zone

Hybrid zones provide unique opportunities to examine reproductive isolation and introgression in nature. We utilized 45,384 single nucleotide polymorphism (SNP) loci to perform association mapping of 14 floral, vegetative and ecological traits that differ between Iris hexagona and Iris fulva, and to investigate, using a Bayesian genomic cline (BGC) framework, patterns of genomic introgression in a large and phenotypically diverse hybrid zone in southern Louisiana. Many loci of small effect size were consistently found to be associated with phenotypic variation across all traits, and several individual loci were revealed to influence phenotypic variation across multiple traits. Patterns of genomic introgression were quite heterogeneous throughout the Louisiana Iris genome, with I. hexagona alleles tending to be favoured over those of I. fulva. Loci that were found to have exceptional patterns of introgression were also found to be significantly associated with phenotypic variation in a small number of morphological traits. However, this was the exception rather than the rule, as most loci that were associated with morphological trait variation were not significantly associated with excess ancestry. These findings provide insights into the complexity of the genomic architecture of phenotypic differences and are a first step towards identifying loci that are associated with both trait variation and reproductive isolation in nature.     

Genome sequencing and comparison of two nonhuman primate animal models, the cynomolgus and Chinese rhesus macaques

The nonhuman primates most commonly used in medical research are from the genus Macaca. To better understand the genetic differences between these animal models, we present high-quality draft genome sequences from two macaque species, the cynomolgus/crab-eating macaque and the Chinese rhesus macaque. Comparison with the previously sequenced Indian rhesus macaque reveals that all three macaques maintain abundant genetic heterogeneity, including millions of single-nucleotide substitutions and many insertions, deletions and gross chromosomal rearrangements. By assessing genetic regions with reduced variability, we identify genes in each macaque species that may have experienced positive selection. Genetic divergence patterns suggest that the cynomolgus macaque genome has been shaped by introgression after hybridization with the Chinese rhesus macaque. Macaque genes display a high degree of sequence similarity with human disease gene orthologs and drug targets. However, we identify several putatively dysfunctional genetic differences between the three macaque species, which may explain functional differences between them previously observed in clinical studies.     

A comprehensive evaluation of cattle introgression into US federal bison herds

Genetic introgression, especially from interspecies hybridization, is a significant threat to species conservation worldwide. In this study, 11 US federal bison populations were comprehensively examined for evidence of both mitochondrial and nuclear domestic cattle (Bos taurus) introgression. Mitochondrial introgression was examined using established polymerase chain reaction methods and confirmed through analysis of D-loop sequences. Nuclear introgression was assessed in 14 chromosomal regions through examination of microsatellite electromorph and sequence differences between bison and domestic cattle. Only one population was identified with domestic cattle mitochondrial DNA introgression. In contrast, evidence of nuclear introgression was found in 7 (63.6%) of the examined populations. Historic accounts of bison transfers among populations were corroborated with evidence of introgressed DNA transmission. While neither nuclear nor mitochondrial domestic cattle introgression was detected in bison from Grand Teton National Park, Sully's Hill National Game Preserve, Wind Cave National Park, or Yellowstone National Park, adequate sample sizes were available only from the last 2 populations to allow for statistical confidence (>90%) in nuclear introgression detection limits. The identification of genetically unique and undisturbed populations is critical to species conservation efforts, and this study serves as a model for the genetic evaluation of interspecies introgression.     

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.