Do Heliconius butterfly species exchange mimicry alleles?

Hybridization has the potential to transfer beneficial alleles across species boundaries, and there are a growing number of examples in which this has apparently occurred. Recent studies suggest that Heliconius butterflies have transferred wing pattern mimicry alleles between species via hybridization, but ancestral polymorphism could also produce a signature of shared ancestry around mimicry genes. To distinguish between these alternative hypotheses, we measured DNA sequence divergence around putatively introgressed mimicry loci and compared this with the rest of the genome. Our results reveal that putatively introgressed regions show strongly reduced sequence divergence between co-mimetic species, suggesting that their divergence times are younger than the rest of the genome. This is consistent with introgression and not ancestral variation. We further show that this signature of introgression occurs at sites throughout the genome, not just around mimicry genes.     

Gene flow-dependent genomic divergence between Anopheles gambiae M and S forms

Anopheles gambiae sensu stricto exists as two often-sympatric races termed the M and S molecular forms, characterized by fixed differences at an X-linked marker. Extreme divergence between M and S forms at pericentromeric "genomic islands" suggested that selection on variants therein could be driving interform divergence in the presence of ongoing gene flow, but recent work has detected much more widespread genomic differentiation. Whether such genomic islands are important in reproductive isolation or represent ancestral differentiation preserved by low recombination is currently unclear. A critical test of these competing hypotheses could be provided by comparing genomic divergence when rates of recent introgression vary. We genotyped 871 single nucleotide polymorphisms (SNPs) in A. gambiae sensu stricto from locations of M and S sympatry and allopatry, encompassing the full range of observed hybridization rates (0-25%). M and S forms were readily partitioned based on genomewide SNP variation in spite of evidence for ongoing introgression that qualitatively reflects hybridization rates. Yet both the level and the heterogeneity of genomic divergence varied markedly in line with levels of introgression. A few genomic regions of differentiation between M and S were common to each sampling location, the most pronounced being two centromere-proximal speciation islands identified previously but with at least one additional region outside of areas expected to exhibit reduced recombination. Our results demonstrate that extreme divergence at genomic islands does not simply represent segregating ancestral polymorphism in regions of low recombination and can be resilient to substantial gene flow. This highlights the potential for islands comprising a relatively small fraction of the genome to play an important role in early-stage speciation when reproductive isolation is limited.     

A genetic map in the Mimulus guttatus species complex reveals transmission ratio distortion due to heterospecific interactions.

As part of a study of the genetics of floral adaptation and speciation in the Mimulus guttatus species complex, we constructed a genetic linkage map of an interspecific cross between M. guttatus and M. nasutus. We genotyped an F(2) mapping population (N = 526) at 255 AFLP, microsatellite, and gene-based markers and derived a framework map through repeated rounds of ordering and marker elimination. The final framework map consists of 174 marker loci on 14 linkage groups with a total map length of 1780 cM Kosambi. Genome length estimates (2011-2096 cM) indicate that this map provides thorough coverage of the hybrid genome, an important consideration for QTL mapping. Nearly half of the markers in the full data set (49%) and on the framework map (48%) exhibited significant transmission ratio distortion (alpha = 0.05). We localized a minimum of 11 transmission ratio distorting loci (TRDLs) throughout the genome, 9 of which generate an excess of M. guttatus alleles and a deficit of M. nasutus alleles. This pattern indicates that the transmission ratio distortion results from particular interactions between the heterospecific genomes and suggests that substantial genetic divergence has occurred between these Mimulus species. We discuss possible causes of the unequal representation of parental genomes in the F(2) generation.     

Transmission ratio distortion in intraspecific hybrids of Mimulus guttatus: implications for genomic divergence

We constructed a genetic linkage map between two divergent populations of Mimulus guttatus. We genotyped an F(2) mapping population (N = 539) at 154 AFLP, microsatellite, and gene-based markers. A framework map was constructed consisting of 112 marker loci on 14 linkage groups with a total map length of 1518 cM Kosambi. Nearly half of all markers (48%) exhibited significant transmission ratio distortion (alpha = 0.05). By using a Bayesian multipoint mapping method and visual inspection of significantly distorted markers, we detected 12 transmission ratio distorting loci (TRDL) throughout the genome. The high degree of segregation distortion detected in this intraspecific map indicates substantial genomic divergence that perhaps suggests genomic incompatibilities between these two populations. We compare the pattern of transmission ratio distortion in this map to an interspecific map constructed between M. guttatus and M. nasutus. A similar level of segregation distortion is detected in both maps. Collinear regions between maps are compared to determine if there are shared genetic patterns of non-Mendelian segregation distortion within and among Mimulus species.     

Do highly divergent loci reside in genomic regions affecting reproductive isolation? A test using next-generation sequence data in Timema stick insects

ABSTRACT: BACKGROUND: Genetic divergence during speciation with gene flow is heterogeneous across the genome, with some regions exhibiting stronger differentiation than others. Exceptionally differentiated regions are often assumed to experience reduced introgression, i.e., reduced flow of alleles from one population into another because such regions are affected by divergent selection or cause reproductive isolation. In contrast, the remainder of the genome can be homogenized by high introgression. Although many studies have documented variation across the genome in genetic differentiation, there are few tests of this hypothesis that explicitly quantify introgression. Here, we provide such a test using 38,304 SNPs in populations of Timema cristinae stick insects. We quantify whether loci that are highly divergent between geographically separated ('allopatric') populations exhibit unusual patterns of introgression in admixed populations. To the extent this is true, highly divergent loci between allopatric populations contribute to reproductive isolation in admixed populations. RESULTS: As predicted, we find a substantial association between locus-specific divergence between allopatric populations and locus-specific introgression in admixed populations. However, many loci depart from this relationship, sometimes strongly so. We also report evidence for selection against foreign alleles due to local adaptation. CONCLUSIONS: Loci that are strongly differentiated between allopatric populations sometimes contribute to reproductive isolation in admixed populations. However, geographic variation in selection and local adaptation, in aspects of genetic architecture (such as organization of genes, recombination rate variation, number and effect size of variants contributing to adaptation, etc.), and in stochastic evolutionary processes such as drift can cause strong differentiation of loci that do not always contribute to reproductive isolation. The results have implications for the theory of 'genomic islands of speciation'.     

Insight into the roles of selection in speciation from genomic patterns of divergence and introgression in secondary contact in venomous rattlesnakes

Investigating secondary contact of historically isolated lineages can provide insight into how selection and drift influence genomic divergence and admixture. Here, we studied the genomic landscape of divergence and introgression following secondary contact between lineages of the Western Diamondback Rattlesnake (Crotalus atrox) to determine whether genomic regions under selection in allopatry also contribute to reproductive isolation during introgression. We used thousands of nuclear loci to study genomic differentiation between two lineages that have experienced recent secondary contact following isolation, and incorporated sampling from a zone of secondary contact to identify loci that are resistant to gene flow in hybrids. Comparisons of patterns of divergence and introgression revealed a positive relationship between allelic differentiation and resistance to introgression across the genome, and greater‐than‐expected overlap between genes linked to lineage‐specific divergence and loci that resist introgression. Genes linked to putatively selected markers were related to prominent aspects of rattlesnake biology that differ between populations of Western Diamondback rattlesnakes (i.e., venom and reproductive phenotypes). We also found evidence for selection against introgression of genes that may contribute to cytonuclear incompatibility, consistent with previously observed biased patterns of nuclear and mitochondrial alleles suggestive of partial reproductive isolation due to cytonuclear incompatibilities. Our results provide a genome‐scale perspective on the relationships between divergence and introgression in secondary contact that is relevant for understanding the roles of selection in maintaining partial isolation of lineages, causing admixing lineages to not completely homogenize.     

Gene-Flow in a Mosaic Hybrid Zone: Is Local Introgression Adaptive?

Genome-wide scans of genetic differentiation between hybridizing taxa can identify genome regions with unusual rates of introgression. Regions of high differentiation might represent barriers to gene flow, while regions of low differentiation might indicate adaptive introgression—the spread of selectively beneficial alleles between reproductively isolated genetic backgrounds. Here we conduct a scan for unusual patterns of differentiation in a mosaic hybrid zone between two mussel species, Mytilus edulis and M. galloprovincialis. One outlying locus, mac-1, showed a characteristic footprint of local introgression, with abnormally high frequency of edulis-derived alleles in a patch of M. galloprovincialis enclosed within the mosaic zone, but low frequencies outside of the zone. Further analysis of DNA sequences showed that almost all of the edulis allelic diversity had introgressed into the M. galloprovincialis background in this patch. We then used a variety of approaches to test the hypothesis that there had been adaptive introgression at mac-1. Simulations and model fitting with maximum-likelihood and approximate Bayesian computation approaches suggested that adaptive introgression could generate a “soft sweep,” which was qualitatively consistent with our data. Although the migration rate required was high, it was compatible with the functioning of an effective barrier to gene flow as revealed by demographic inferences. As such, adaptive introgression could explain both the reduced intraspecific differentiation around mac-1 and the high diversity of introgressed alleles, although a localized change in barrier strength may also be invoked. Together, our results emphasize the need to account for the complex history of secondary contacts in interpreting outlier loci.     

Characterization of six microsatellite loci in Myrica faya (Myricaceae) and cross amplification in the endangered endemic M. rivas-martinezii in Canary Islands, Spain

Six novel polymorphic microsatellite markers were isolated from enriched libraries in Myrica faya Ait., recently renamed Morella faya, (fayatree, firetree, or firebush) in order to examine the genetic diversity in natural populations. Also, test cross-specific amplification and genetic diversity in Myrica rivas-martinezii, which is endemic on the Canary islands. Microsatellite loci were screened in 225 individuals of both species from different islands of the Canarian archipelago. All markers were successfully amplified from both Myrica species, with an average number of 6.5 and 9.3 alleles per locus in M. rivas-martinezii and M. faya, respectively. There was no evidence for linkage disequilibrium between loci, and the probability of null alleles ranged from 0.01 to 0.17.     

A library of Solanum lycopersicoides introgression lines in cultivated tomato.

A set of introgression lines (ILs), containing individual chromosome segments from the wild nightshade Solanum lycopersicoides bred into the genetic background of cultivated tomato (Lycopersicon esculentum), has been developed. A primary group of 56 lines was selected for maximum representation of the S. lycopersicoides genome (approximately 96% of the total map units), homozygosity, and a minimum number of introgressed segments per line. A secondary set of 34 lines provides increased map resolution in certain regions. Approximately 34% of the lines were sterile in the homozygous condition, but could be maintained by heterozygotes. To facilitate identification of segregating ILs, restriction fragment length polymorphism probes were converted to higher throughput cleaved amplified polymorphic sequence markers, which supplement allozyme and morphological loci. Strong segregation distortion was observed in F2 progeny of heterozygous ILs, with an excess of L. esculentum alleles in most regions. For introgressions on distal chromosome 1L, a preferential transmission of S. lycopersicoides alleles was observed in the male germ line. Homozygous ILs generally yielded less seed from self pollination than corresponding heterozygotes, indicating that sterility effects were recessive. This IL library provides a novel resource for genetic studies of traits found in S. lycopersicoides.     

THE VARIABLE GENOMIC ARCHITECTURE OF ISOLATION BETWEEN HYBRIDIZING SPECIES OF HOUSE MICE

Studies of the genetics of hybrid zones can provide insight into the genomic architecture of species boundaries. By examining patterns of introgression of multiple loci across a hybrid zone, it may be possible to identify regions of the genome that have experienced selection. Here, we present a comparison of introgression in two replicate transects through the house mouse hybrid zone through central Europe, using data from 41 single nucleotide markers. Using both genomic and geographic clines, we found many differences in patterns of introgression between the two transects, as well as some similarities. We found that many loci may have experienced the effects of selection at linked sites, including selection against hybrid genotypes, as well as positive selection in the form of genotypes introgressed into a foreign genetic background. We also found many positive associations of conspecific alleles among unlinked markers, which could be caused by epistatic interactions. Different patterns of introgression in the two transects highlight the challenge of using hybrid zones to identify genes underlying isolation and raise the possibility that the genetic basis of isolation between these species may be dependent on the local population genetic make-up or the local ecological setting.     

物种形成过程中的分化基因组岛及其形成机制

理解物种形成机制是生态和进化领域的重要任务。得益于测序技术的快速发展, 越来越多研究发现分化种群(亚种、物种)间的基因组常呈现异质性分化景观, 存在分化基因组岛, 这被认为是基因流存在下的歧化选择引起的, 支持基因流存在下的成种假说。然而, 基因渐渗、祖先多态性的差异分选、连锁选择等其他进化过程也可导致分化基因组岛的形成。现有实证研究在解析分化基因组岛的形成机制时, 往往忽略了上述其他进化过程的作用。为此, 本文在辨析分化基因组岛相关概念的基础上, 总结了利用种群基因组数据鉴定分化基因组岛的方法, 对比了不同进化过程形成分化基因组岛的特征, 指出在区分不同机制时联用基因渐渗程度、绝对分化指数(d_XY)、相对节点深度(RND)、重组率等多个指标的必要性, 归纳了物种形成过程中分化基因组岛形成机制解析的研究思路, 并对未来在生殖隔离机制上的深入探索以及实证研究的整合分析等方面进行了展望。     

Evidence for bimodal hybrid zones between two species of char (Pisces: Salvelinus) in northwestern North America

Dolly Varden (Salvelinus malma, Pisces: Salmonidae) and bull trout (Salvelinus confluentus) have widely overlapping, but largely parapatric ranges in watersheds in northwestern North America from Washington State to northern British Columbia. Genetic analysis of natural populations using diagnostic molecular markers revealed widespread local sympatry and hybridization with hybrids comprising 0-25% of the local samples. In a detailed analysis of hybridization using four nuclear DNA markers and mitochondrial DNA within the Thutade Lake watershed, northcentral British Columbia, hybrid genotypes constituted up to 9% of the population of juvenile char. There were significant deviations from Hardy-Weinberg, gametic, and cytonuclear equilibria, and local samples showed bimodal frequency distributions of genotypes. Pure parental and inferred backcross genotypes were most common, and F1 and F(n) hybrids were comparatively rare. Interspecific hybridization was asymmetrical, with most F1 hybrids (five of six) bearing S. confluentus mtDNA. The introgression of nuclear and mitochondrial alleles was asymmetrical, with S. confluentus mtDNA and Growth Hormone 2 introgressing into S. malma significantly more than either introgression of the three other nuclear loci, or introgression of S. malma alleles into S. confluentus. Substantial prezygotic isolation between the species likely depends on the large body size difference between them in sympatry: S. malma have small bodies and a stream resident life history (12-21 cm adult fork length at maturity), while S. confluentus are larger and adfluvial, i.e., they migrate to Thutade Lake where they grow to maturity before returning to tributary streams to spawn (40-90 cm at maturity). These traits may limit interspecific pairings because of size assortative pairing and size-dependent reproductive habitat use.     

Genetics of postzygotic isolation in Eucalyptus: whole-genome analysis of barriers to introgression in a wide interspecific cross of Eucalyptus grandis and E. globulus

The genetic architecture of hybrid fitness characters can provide valuable insights into the nature and evolution of postzygotic reproductive barriers in diverged species. We determined the genome-wide distribution of barriers to introgression in an F(1) hybrid of two Eucalyptus tree species, Eucalyptus grandis (W. Hill ex Maiden.) and E. globulus (Labill.). Two interspecific backcross families (N = 186) were used to construct comparative, single-tree, genetic linkage maps of an F(1) hybrid individual and two backcross parents. A total of 1354 testcross AFLP marker loci were evaluated in the three parental maps and a substantial proportion (27.7% average) exhibited transmission ratio distortion (alpha = 0.05). The distorted markers were located in distinct regions of the parental maps and marker alleles within each region were all biased toward either of the two parental species. We used a Bayesian approach to estimate the position and effect of transmission ratio distorting loci (TRDLs) in the distorted regions of each parental linkage map. The relative viability of TRDL alleles ranged from 0.20 to 0.72. Contrary to expectation, heterospecific (donor) alleles of TRDLs were favored as often as recurrent alleles in both backcrosses, suggesting that positive and negative heterospecific interactions affect introgression rates in this wide interspecific pedigree.     

Presence of two mitochondrial genomes in the mytilid Perumytilus purpuratus: Phylogenetic evidence for doubly uniparental inheritance

This study presents evidence, using sequences of ribosomal 16S and COI mtDNA, for the presence of two mitochondrial genomes in Perumytilus purpuratus. This may be considered evidence of doubly uniparental mtDNA inheritance. The presence of the two types of mitochondrial genomes differentiates females from males. The F genome was found in the somatic and gonadal tissues of females and in the somatic tissues of males; the M genome was found in the gonads and mantle of males only. For the mitochondrial 16S region, ten haplotypes were found for the F genome (nucleotide diversity 0.004), and 7 haplotypes for the M genome (nucleotide diversity 0.001), with a distance Dxy of 0.125 and divergence Kxy of 60.33%. For the COI gene 17 haplotypes were found for the F genome (nucleotide diversity 0.009), and 10 haplotypes for the M genome (nucleotide diversity 0.010), with a genetic distance Dxy of 0.184 and divergence Kxy of 99.97%. Our results report the presence of two well-differentiated, sex-specific types of mitochondrial genome (one present in the male gonad, the other in the female gonad), implying the presence of DUI in P. purpuratus. These results indicate that care must be taken in phylogenetic comparisons using mtDNA sequences of P. purpuratus without considering the sex of the individuals.     

Tracking the progression of speciation: variable patterns of introgression across the genome provide insights on the species delimitation between progenitor–derivative spruces (Picea mariana × P. rubens)

The genic species concept implies that while most of the genome can be exchanged somewhat freely between species through introgression, some genomic regions remain impermeable to interspecific gene flow. Hence, interspecific differences can be maintained despite ongoing gene exchange within contact zones. This study assessed the heterogeneous patterns of introgression at gene loci across the hybrid zone of an incipient progenitor–derivative species pair, Picea mariana (black spruce) and Picea rubens (red spruce). The spruce taxa likely diverged in geographic isolation during the Pleistocene and came into secondary contact during late Holocene. A total of 300 SNPs distributed across the 12 linkage groups (LG) of black spruce were genotyped for 385 individual trees from 33 populations distributed across the allopatric zone of each species and within the zone of sympatry. An integrative framework combining three population genomic approaches was used to scan the genomes, revealing heterogeneous patterns of introgression. A total of 23 SNPs scattered over 10 LG were considered impermeable to introgression and putatively under diverging selection. These loci revealed the existence of impermeable genomic regions forming the species boundary and are thus indicative of ongoing speciation between these two genetic lineages. Another 238 SNPs reflected selectively neutral diffusion across the porous species barrier. Finally, 39 highly permeable SNPs suggested ancestral polymorphism along with balancing selection. The heterogeneous patterns of introgression across the genome indicated that the speciation process between black spruce and red spruce is young and incomplete, albeit some interspecific differences are maintained, allowing ongoing species divergence even in sympatry. The approach developed in this study can be used to track the progression of ongoing speciation processes.     

COMPARATIVE GENOMIC AND POPULATION GENETIC ANALYSES INDICATE HIGHLY POROUS GENOMES AND HIGH LEVELS OF GENE FLOW BETWEEN DIVERGENT HELIANTHUS SPECIES

While speciation can be found in the presence of gene flow, it is not clear what impact this gene flow has on genome- and range-wide patterns of differentiation. Here we examine gene flow across the entire range of the common sunflower, H. annuus , its historically allopatric sister species H. argophyllus and a more distantly related, sympatric relative H. petiolaris . Analysis of genotypes at 26 microsatellite loci in 1015 individuals from across the range of the three species showed substantial introgression between geographically proximal populations of H. annuus and H. petiolaris , limited introgression between H. annuus and H. argophyllus , and essentially no gene flow between the allopatric pair, H. argophyllus and H. petiolaris. Analysis of sequence divergence levels among the three species in 1420 orthologs identified from EST databases identified a subset of loci showing extremely low divergence between H. annuus and H. petiolaris and extremely high divergence between the sister species H. annuus and H. argophyllus , consistent with introgression between H. annuus and H. petiolaris at these loci. Thus, at many loci, the allopatric sister species are more genetically divergent than the more distantly related sympatric species, which have exchanged genes across much of the genome while remaining morphologically and ecologically distinct.     

Fitness landscapes support the dominance theory of post-zygotic isolation in the mussels Mytilus edulis and M. galloprovincialis

We studied the genetic basis of post-zygotic isolation in the marine mussels Mytilus edulis and Mytilus galloprovincialis. Evidence was obtained for a high number of recessive Dobzhansky-Muller substitutions in the genome of these two mussel taxa. We analysed the segregation of unlinked diagnostic markers in the progeny of two backcrosses and an F2 cross, 36 h and 200 days after fertilization. Directional selection favouring M. galloprovincialis genotypes was observed in both kinds of cross. In the F2, epistatic interactions between each pair of chromosome fragments mapped by the markers were identified in addition. Our results imply that homozygous-homozygous interactions are required for breakdown of coadaptation, in accordance with the dominance theory of post-zygotic isolation. Endogenous post-zygotic selection distributed over many loci throughout the genome provides the missing factor explaining the astonishing persistence and strength of barriers to neutral introgression in such a dispersive taxon as Mytilus.     

Genetic association of physically unlinked islands of genomic divergence in incipient species of Anopheles gambiae

Previous efforts to uncover the genetic underpinnings of ongoing ecological speciation of the M and S forms of the African malaria vector Anopheles gambiae revealed two centromere‐proximal islands of genetic divergence on X and chromosome 2. Under the assumption of considerable ongoing gene flow between M and S, these persistently divergent genomic islands were widely considered to be ‘speciation islands’. In the course of microarray‐based divergence mapping, we discovered a third centromere‐associated island of divergence on chromosome 3, which was validated by targeted re‐sequencing. To test for genetic association between the divergence islands on all three chromosomes, SNP‐based assays were applied in four natural populations of M and S spanning West, Central and East Africa. Genotyping of 517 female M and S mosquitoes revealed nearly complete linkage disequilibrium between the centromeres of the three independently assorting chromosomes. These results suggest that despite the potential for inter‐form gene flow through hybridization, actual (realized) gene flow between M and S may be substantially less than commonly assumed and may not explain most shared variation. Moreover, the possibility of very low gene flow calls into question whether diverged pericentromeric regions—characterized by reduced levels of variation and recombination—are in fact instrumental rather than merely incidental to the speciation process.     

Detection of the introgression of genome elements of the Aegilops cylindrica host. into the Triticum aestivum L. genome by ISSR and SSR analysis

To reveal sites of the donor genome in wheat crossed with Aegilops cylindrica, which acquired conferred resistance to fungal diseases, a comparative analysis of introgressive and parental forms was conducted. Two systems of PCR analysis, ISSR and SSR-PCR, were employed. Upon use of 7 ISSR primers in genotypes of 30 individual plants BC1 F9 belonging to lines 5/55-91 and 5/20-91, 19 ISSR loci were revealed and assigned to introgressive fragments of Aegilops cylindrica genome in Triticum aestivum. The 40 pairs of SSR primers allowed the detection of seven introgressive alleles; three of these alleles were located on common wheat chromosomes in the B genome, while four alleles, in the D genome. Based on data of microsatellite analysis, it was assumed that the telomeric region of the long arm of common wheat chromosome 6A also changed. ISSR and SSR methods were shown to be effective for detecting variability caused by introgression of foreign genetic material into the genome of common wheat.