Variable B chromosomes frequencies between males and females of two species of pufferfishes (Tetraodontiformes)

Cytogenetic analysis in two sympatric species of pufferfish from Paraná Coast, Brazil, Sphoeroides spengleri (Tetraodontidae) and Chilomycterus spinosus (Diodontidae) revealed the presence of 2n = 46 and 2n = 50 chromosomes, respectively. In S. spengleri, the double FISH technique with 5S and 18S rDNA probes showed no syntenical organization, with both multigene families occupying a terminal position on the chromosomes. In S. spengleri the poor pattern of heterochromatin distribution is discussed about its relation with the compactness of their genome. On the other hand C. spinosus presented conspicuous blocks occupying the short arms of several chromosomes, possibly consequence of an accumulative process. Heterochromatic supernumerary chromosomes were detected in both species, however in S. spengleri they were presented as microchromosomes varying of 0–3 and being more frequent in females than in males, possibly as result of a preferential segregation. Conversely in C. spinosus an extra chromosome was restricted to the males, which to behave as univalent during meiosis. In view of these peculiarities we hypothesized a possible relationship between B chromosome and sex determination, although others aspects on accumulation, maintenance, and frequency of these extra elements were also considered, focusing on the parasitic model. In spite of the karyotypic diversity and others genetic peculiarities in Tetraodontiformes, this corresponds the third report of the occurrence of chromosomes B in the group.     

Genomic footprints of an old affair: Single nucleotide polymorphism data reveal historical hybridization and the subsequent evolution of reproductive barriers in two recently diverged grasshoppers with partly overlapping distributions

Secondary contact in close relatives can result in hybridization and the admixture of previously isolated gene pools. However, after an initial period of hybridization, reproductive isolation can evolve through different processes and lead to the interruption of gene flow and the completion of the speciation process. Omocestus minutissimus and O. uhagonii are two closely related grasshoppers with partially overlapping distributions in the Central System mountains of the Iberian Peninsula. To analyse spatial patterns of historical and/or contemporary hybridization between these two taxa and understand how species boundaries are maintained in the region of secondary contact, we sampled sympatric and allopatric populations of the two species and obtained genome‐wide single nucleotide polymorphism data using a restriction site‐associated DNA sequencing approach. We used Bayesian clustering analyses to test the hypothesis of contemporary hybridization in sympatric populations and employed a suite of phylogenomic approaches and a coalescent‐based simulation framework to evaluate alternative hypothetical scenarios of interspecific gene flow. Our analyses rejected the hypothesis of contemporary hybridization but revealed past introgression in the area where the distributions of the two species overlap. Overall, these results point to a scenario of historical gene flow after secondary contact followed by the evolution of reproductive isolation that currently prevents hybridization among sympatric populations.     

Maintenance of species boundaries in a Neotropical radiation of Begonia

A major goal of evolutionary biology is to determine the mechanisms generating biodiversity. In Begonia, one of the largest plant genera (1900+ species), it has been postulated that the high number of endemic species is a by‐product of low gene flow among populations, which predisposes the group to speciation. However, this model of divergence requires that reproductive barriers accumulate rapidly among diverging species that overlap in their geographic ranges, otherwise speciation will be opposed by homogenizing gene flow in zones of secondary contact. Here, we test the outcomes of secondary contact in Begonia by genotyping multiple sympatric sites with 12 nuclear and seven plastid loci. We show that three sites of secondary contact between B. heracleifolia and B. nelumbiifolia are highly structured, mostly containing parental genotypes, with few F1 hybrids. A sympatric site between B. heracleifolia and B. sericoneura contains a higher proportion of F1s, but little evidence of introgression. The lack of later‐generation hybrids contrasts with that documented in many other plant taxa, where introgression is extensive. Our results, in conjunction with previous genetic work, show that Begonia demonstrate properties making them exceptionally prone to speciation, at multiple stages along the divergence continuum. Not only are populations weakly connected by gene flow, promoting allopatric speciation, but species often show strong reproductive barriers in secondary contact. Whether similar mechanisms contribute to diversification in other large genera remains to be tested.     

Demography and speciation history of the homoploid hybrid pine Pinus densata on the Tibetan Plateau

Pinus densata is an ecologically successful homoploid hybrid that inhabits vast areas of heterogeneous terrain on the south‐eastern Tibetan Plateau as a result of multiple waves of colonization. Its region of origin, route of colonization onto the plateau and the directions of introgression with its parental species have previously been defined, but little is known about the isolation and divergence history of its populations. In this study, we surveyed nucleotide polymorphism over eight nuclear loci in 19 representative populations of P. densata and its parental species. Using this information and coalescence simulations, we assessed the historical changes in its population size, gene flow and divergence in time and space. The results indicate a late Miocene origin for P. densata associated with the recent uplift of south‐eastern Tibet. The subsequent differentiation between geographical regions of this species began in the late Pliocene and was induced by regional topographical changes and Pleistocene glaciations. The ancestral P. densata population had a large effective population size but the central and western populations were established by limited founders, suggesting that there were severe bottlenecks during the westward migration out of the ancestral hybrid zone. After separating from their ancestral populations, population expansion occurred in all geographical regions especially in the western range. Gene flow in P. densata was restricted to geographically neighbouring populations, resulting in significant differentiation between regional groups. The new information on the divergence and demographic history of P. densata reported herein enhances our understanding of its speciation process on the Tibetan Plateau.     

Hybridisation and introgression between Brassica napus and B. rapa in the Netherlands

We used flow cytometry, chromosome counting and AFLP markers to investigate gene flow from the crop plant oilseed rape, Brassica napus (AACC) to wild B. rapa (AA) in the Netherlands. From 89 B. napus source populations investigated, all near cropping fields or at transhipment sites, only 19 contained a B. rapa population within a 2.5‐km radius. During our survey we found only three populations with F1 hybrids (AAC), as recognized by their nine extra chromosomes and by flow cytometry. These hybrids were all collected in mixed populations where the two species grew in close proximity. Populations with F1 hybrids were not close to crops, but instead were located on road verges with highly disturbed soils, in which both species were probably recruited from the soil seed bank. Many plants in the F2, BC1 or higher backcrosses are expected to carry one to eight C chromosomes. However, these plants were not observed among the hybrids. We further investigated introgression with molecular markers (AFLP) and compared sympatric B. rapa populations (near populations of B. napus) with control populations of B. rapa (no B. napus within at least 7 km). We found no difference between sympatric and control populations in the number of C markers in B. rapa, nor did we find that these sympatric populations closely resembled B. napus. Our data show that hybrids occur but also suggest no recent introgression of alleles from the crop plant B. napus into wild B. rapa in the Dutch populations studied.     

Maintaining their genetic distance: Little evidence for introgression between widely hybridizing species of Geum with contrasting mating systems

Within the plant kingdom, many genera contain sister lineages with contrasting outcrossing and inbreeding mating systems that are known to hybridize. The evolutionary fate of these sister lineages is likely to be influenced by the extent to which they exchange genes. We measured gene flow between outcrossing Geum rivale and selfing Geum urbanum, sister species that hybridize in contemporary populations. We generated and used a draft genome of G. urbanum to develop dd‐RAD data scorable in both species. Coalescent analysis of RAD data from allopatric populations indicated that the species diverged 2–3 Mya, and that historical gene flow between them was extremely low (1 migrant every 25 generations). Comparison of genetic divergence between species in sympatry and allopatry, together with an analysis of allele frequencies in potential parental and hybrid populations, provided no evidence of contemporary introgression in sympatric populations. Cluster‐ and species‐specific marker analyses revealed that, apart from four early‐generation hybrids, individuals in sympatric populations fell into two genetically distinct groups that corresponded exactly to their morphological species classification with maximum individual admixture estimates of only 1–3%. However, we did observe joint segregation of four putatively introgressed SNPs across two scaffolds in the G. urbanum population that was associated with significant morphological variation, interpreted as tentative evidence for rare, recent interspecific gene flow. Overall, our results indicate that despite the presence of hybrids in contemporary populations, genetic exchange between G. rivale and G. urbanum has been extremely limited throughout their evolutionary history.     

Genetic differentiation between two varieties of Oreocharis benthamii (Gesneriaceae) in sympatric and allopatric regions

The pattern of genetic differentiation between diverging species receives much attention as one of the key observable features of speciation. It has often been suggested that introgression between closely related species occurs commonly where their distributions overlap, leading to their becoming more morphologically and genetically similar, but there are a few opposite results. However, most of these studies have been carried out with animals and separate species; few have looked at intraspecific cases, especially in plants. Here, we conduct a comparative study on patterns of genetic differentiation among populations of two varieties of Oreocharis benthamii in allopatry and sympatry based on ISSR data for 754 individuals from 26 populations, in order to understand the processes leading to speciation. Contrary to expectations, the facultative xenogamy (mixed mating) species O. benthamii has a relatively low genetic diversity within populations (H = 0.1014, I = 0.1528) and high genetic differentiation among populations (G_ST = 0.5867, Ф_ST = 0.659), as is typically found for selfing species. Genetic variance between the two varieties in sympatric populations (44%, Ф_ST = 0.444) is significantly more than that in allopatric populations (14%, Ф_ST = 0.138). Consistent with the taxonomical delimitation of the two varieties, all sampled individuals of O. benthamii clustered into two genetic groups. Moreover, the genetic structures of populations of both varieties are correlated with their different geographical origins. Our studies show that significant divergence between sympatric populations of the two varieties could be attributed primarily to reinforcement by genetic divergent selection in sympatry where secondary contact had occurred. The major proportion of the genetic variation in outcrossing and mixed mating plants may exist among populations when the populations are distributed in fragmented habitats, due to the paucity of suitable habitat combined with inefficient seed dispersal mechanism and limited pollinator foraging area that may limit the gene flow.     

The role of hybridization and introgression in maintaining species integrity and cohesion in naturally isolated inselberg bromeliad populations

• Hybridization is a widespread phenomenon present in numerous lineages across the tree of life. Its evolutionary consequences range from effects on the origin and maintenance, to the loss of biodiversity.
• We studied genetic diversity and intra‐ and interspecific gene flow between two sympatric populations of closely‐related species, Pitcairnia flammea and P. corcovadensis (Bromeliaceae), which are adapted to naturally fragmented Neotropical inselbergs, based on nuclear and plastidial DNA.
• Our main results indicate a strong reproductive isolation barrier, although low levels of interspecific gene flow were observed in both sympatric populations. The low rates of intraspecific gene flow observed for both P. corcovadensis and P. flammea populations corroborate the increasing body of evidence that inselberg bromeliad species are maintained as discrete evolutionary units despite the presence of low genetic connectivity. Nuclear patterns of genetic diversity and gene flow revealed that hybridization and introgression might not cause species extinction via genetic assimilation of the rare P. corcovadensis.
• In the face of reduced intraspecific gene exchange, hybridization and introgression may be important aspects of the Pitcairnia diversification process, with a positive evolutionary impact at the bromeliad community level, and thus contribute to increasing and maintaining genetic diversity in local isolated inselberg populations.

     

Rapid sexual and genomic isolation in sympatric Drosophila without reproductive character displacement

The rapid evolution of sexual isolation in sympatry has long been associated with reinforcement (i.e., selection to avoid maladaptive hybridization). However, there are many species pairs in sympatry that have evolved rapid sexual isolation without known costs to hybridization. A major unresolved question is what evolutionary processes are involved in driving rapid speciation in such cases. Here, we focus on one such system; the Drosophila athabasca species complex, which is composed of three partially sympatric and interfertile semispecies: WN, EA, and EB. To study speciation in this species complex, we assayed sexual and genomic isolation within and between these semispecies in both sympatric and allopatric populations. First, we found no evidence of reproductive character displacement (RCD) in sympatric zones compared to distant allopatry. Instead, semispecies were virtually completely sexually isolated from each other across their entire ranges. Moreover, using spatial approaches and coalescent demographic simulations, we detected either zero or only weak heterospecific gene flow in sympatry. In contrast, within each semispecies we found only random mating and little population genetic structure, except between highly geographically distant populations. Finally, we determined that speciation in this system is at least an order of magnitude older than previously assumed, with WN diverging first, around 200K years ago, and EA and EB diverging 100K years ago. In total, these results suggest that these semispecies should be given full species status and we adopt new nomenclature: WN—D. athabasca, EA—D. mahican, and EB—D. lenape. While the lack of RCD in sympatry and interfertility do not support reinforcement, we discuss what additional evidence is needed to further decipher the mechanisms that caused rapid speciation in this species complex.     

Molecular population genetics of a host‐associated sibling species complex of phytophagous ladybird beetles (Coleoptera: Coccinellidae: Epilachninae)

To clarify evolutionary relationships and assess the intensity and persistence of reproductive isolation because of host fidelity in three closely related host‐associated phytophagous ladybird beetles (Henosepilachna pustulosa, Henosepilachna niponica and Henosepilachna yasutomii), we determined sequences of part of the nuclear DNA internal transcribed spacer‐II (ITS‐2) region (388 bp long) for 70 individuals from seven populations, and part of the mitochondrial cytochrome c oxidase subunit I (COI) gene (985 bp long) for 280 individuals from 28 populations across the three species. We failed to detect any polymorphic sites in the ITS‐2 region. We detected 59 COI haplotypes, among which two were shared by different species. While a haplotype network and a phylogenetic gene tree of the COI haplotypes did not support monophyly for any of the three species, the results of AMOVA did not invalidate the classification of these beetles into three species. The isolation‐with‐migration analytic (IMa) test suggested that gene flow between allopatric species (i.e., H. pustulosa versus H. niponica; H. pustulosa versus H. yasutomii) was well prevented, probably because of their limited dispersal power. On the other hand, the IMa test detected a low level of unidirectional gene flow between the sympatric species H. niponica and H. yasutomii, from the former to the latter. This result was consistent with a survival rate of F1 hybrids between the two species that is higher on the host plant of H. yasutomii than on the host plant of H. niponica. High F_ST values among two sympatric and three nearly sympatric population pairs of H. niponica and H. yasutomii, however, indicated that gene flow between the two species has been quite restricted even if it has occurred. Because no reproductive barriers other than the difference in host plants is known for sympatric H. niponica and H. yasutomii, our results suggest that host differentiation and associated ecological divergence function as an effective reproductive barrier between the two species.     

Genetic differentiation between the ant Myrmica rubra and its microgynous social parasite

Hymenopteran inquiline species have been proposed to originate by sympatric speciation through intraspecific social parasitism. One such parasite, Myrmica microrubra, was recently synonymized with its Myrmica rubra host, because comparisons across Europe indicated insufficient genetic differentiation. Here, we use microsatellite markers to study genetic differentiation more precisely in a sample of Finnish M. rubra and its inquilines collected at two localities, supplemented with mitochondrial DNA sequences. The parasite had much lower genetic variation than the host at three of the four loci studied. Genetic differentiation between the host populations was moderate (F _ST = 0.089), whereas the parasite populations were more strongly subdivided (F _ST = 0.440). The host and parasite were highly genetically differentiated both across populations (F _ST = 0.346) and in strict sympatry (0.327, 0.364), a result that remained robust both in a haplotype network and in PCA ordination. Individual assignments of genotypes indicated that gene flow between sympatric host and inquiline populations is reduced by about an order of magnitude relative to the gene flow within the morphs. Our results suggest that the parasitic morph of M. rubra may be an incipient species, but it remains unclear to what extent the observed genetic differentiation between host and inquiline is due to possible assortative mating and selection against hybrids or to recurrent bottlenecking and genetic drift. We conclude that an explicitly functional species concept would be unambiguous in treating this inquiline as a full species, as it begets its own kind and maintains its integrity in spite of occasional interbreeding with the host.     

Sympatric diversification vs. immigration: deciphering host‐plant specialization in a polyphagous insect,the stolbur phytoplasma vector Hyalesthes obsoletus (Cixiidae)

The epidemiology of vector transmitted plant diseases is highly influenced by dispersal and the host‐plant range of the vector. Widening the vector's host range may increase transmission potential, whereas specialization may induce specific disease cycles. The process leading to a vector's host shift and its epidemiological outcome is therefore embedded in the frameworks of sympatric evolution vs. immigration of preadapted populations. In this study, we analyse whether a host shift of the stolbur phytoplasma vector, Hyalesthes obsoletus from field bindweed to stinging nettle in its northern distribution range evolved sympatrically or by immigration. The exploitation of stinging nettle has led to outbreaks of the grapevine disease bois noir caused by a stinging nettle‐specific phytoplasma strain. Microsatellite data from populations from northern and ancestral ranges provide strong evidence for sympatric host‐race evolution in the northern range: Host‐plant associated populations were significantly differentiated among syntopic sites (0.054 < F_HT < 0.098) and constant over 5 years. While gene flow was asymmetric from the old into the predicted new host race, which had significantly reduced genetic diversity, the genetic identity between syntopic host‐race populations in the northern range was higher than between these populations and syntopic populations in ancestral ranges, where there was no evidence for genetic host races. Although immigration was detected in the northern field bindweed population, it cannot explain host‐race diversification but suggests the introduction of a stinging nettle‐specific phytoplasma strain by plant‐unspecific vectors. The evolution of host races in the northern range has led to specific vector‐based bois noir disease cycles.     

Hybridization and introgression across different ploidy levels in the Neotropical orchids Epidendrum fulgens and E. puniceoluteum (Orchidaceae)

The hypothesis of gene flow between species with large differences in chromosome numbers has rarely been tested in the wild, mainly because species of different ploidy are commonly assumed to be reproductively isolated from each other because of instantaneous and strong postzygotic barriers. In this study, a broad‐scale survey of molecular variation was carried out between two orchid species with different ploidy levels: Epidendrum fulgens (2n = 2x = 24 chromosomes) and Epidendrum puniceoluteum (2n = 4x = 52 chromosomes). To test the strength of their reproductive barriers, we investigated the distribution of genetic variation in sympatric and allopatric populations of these two species and conducted crossing experiments. Nuclear and plastid microsatellite loci were used to genotype 463 individuals from eight populations across the geographical range of both species along the Brazilian coastal plain. All six sympatric populations analysed presented hybrid zones, indicating that hybridization between E. fulgens and E. puniceoluteum is a common phenomenon. Bayesian assignment analysis detected the presence of F₁ and F₂ individuals and also signs of introgression, demonstrating a high potential for interspecific gene flow. Introgression occurs preferentially from E. fulgens to E. puniceoluteum. Pure parental individuals of both species display strong genotype–habitat associations, indicating that environment‐dependent selection could be acting in all hybrid zones. This study suggests that hybridization and introgression are evolutionary processes playing a role in the diversification of Epidendrum and indicates the importance of investigations of hybrid zones in understanding reproductive barriers and speciation processes in Neotropical orchid species.     

Spatial and host related genomic variation in partially sympatric cactophagous moth species

Surveys of patterns of genetic variation in natural sympatric and allopatric populations of recently diverged species are necessary to understand the processes driving intra- and interspecific diversification. The South American moths Cactoblastis cactorum, Cactoblastis doddi and Cactoblastis bucyrus are specialized in the use of cacti as host plants. These species have partially different geographic ranges and differ in patterns of host plant use. However, there are areas that overlap, particularly, in northwestern Argentina, where they are sympatric. Using a combination of genome-wide SNPs and mitochondrial data we assessed intra and interspecific genetic variation and investigated the relative roles of geography and host plants on genetic divergence. We also searched for genetic footprints of hybridization between species. We identified three well delimited species and detected signs of hybridization in the area of sympatry. Our results supported a hypothetical scenario of allopatric speciation in the generalist C. cactorum and genetic interchange during secondary geographic contact with the pair of specialists C. bucyrus and C. doddi that probably speciated sympatrically. In both cases, adaptation to new host plants probably played an important role in speciation. The results also suggested the interplay of geography and host plant use as drivers of divergence and limiting gene flow at intra and interspecific levels.     

Forest refugia and riverine barriers promote diversification in the West African pygmy shrew (Crocidura obscurior complex,Soricomorpha)

The Crocidura obscurior or West African pygmy shrew complex is endemic to West African forests from south‐eastern Guinea, eastern Liberia, southern Côte d'Ivoire and south‐western Ghana. We explore the genetic and morphometric diversity of 239 individuals of the C. obscurior complex from 17 localities across its geographical range. Using genetic data from three mitochondrial (16S, cytochrome b and COI) and four nuclear markers (BRCA1, STAT5A, HDAC2 and RIOK3) and skull geometric morphometrics, we show that this complex is composed of two cryptic and sympatric species, C. obscurior and C. eburnea. We then test several hypotheses to infer their evolutionary history. The observed phylogeographical pattern based on cytochrome b and COI sequences fits the forest refuge theory: during arid phases of the Plio‐Pleistocene, around 3.5, 2.1, 1 and 0.5 Mya, a small number of populations survived in isolated forest patches and diverged allopatrically. During wetter climatic periods, forests expanded, leading to secondary contacts between previously isolated populations. Our results also suggest the possible contribution of episodes of isolation in subrefuges. Historical variation of the West African hydrographic network could also have contributed to the observed patterns of genetic differentiation. Rivers such as the Volta and Sassandra may act as past and/or current barriers to gene flow. Although these two species have sympatric distributions, their phylogeographical histories are somewhat dissimilar due to small differences in their dispersal abilities and ecological requirements.     

Whole‐genome sequencing reveals absence of recent gene flow and separate demographic histories for Anopheles punctulatus mosquitoes in Papua New Guinea

Anopheles mosquitoes are the vectors of several human diseases including malaria. In many malaria endemic areas, several species of Anopheles coexist, sometimes in the form of related sibling species that are morphologically indistinguishable. Determining the size and organization of Anopheles populations, and possible ongoing gene flow among them is important for malaria control and, in particular, for monitoring the spread of insecticide resistance alleles. However, these parameters have been difficult to evaluate in most Anopheles species due to the paucity of genetic data available. Here, we assess the extent of contemporary gene flow and historical variations in population size by sequencing and de novo assembling the genomes of wild‐caught mosquitoes from four species of the Anopheles punctulatus group of Papua New Guinea. Our analysis of more than 50 Mb of orthologous DNA sequences revealed no evidence of contemporary gene flow among these mosquitoes. In addition, investigation of the demography of two of the An. punctulatus species revealed distinct population histories. Overall, our analyses suggest that, despite their similarities in morphology, behaviour and ecology, contemporary sympatric populations of An. punctulatus are evolving independently.     

Phylogeny and biogeography of Muusoctopus (Cephalopoda: Enteroctopodidae)

Deep‐sea octopuses of the genus Muusoctopus are thought to have originated in the Pacific Northern Hemisphere and then diversified throughout the Pacific and into the rest of the World Ocean. However, this hypothesis was inferred only from molecular divergence times. Here, the ancestral distribution and dispersal routes are estimated by Bayesian analysis based on a new phylogeny including 38 specimens from the south‐eastern Pacific Ocean. Morphological data and molecular sequences of three mitochondrial genes (16S rRNA, COI and COIII) are presented. The morphological data confirm that specimens newly acquired from off the coast of Chile comprise two species: Muusoctopus longibrachus and the poorly described species, Muusoctopus eicomar. The latter is here redescribed and is clearly distinguished from M. longibrachus and other closely related species in the region. A gene tree was built using Bayesian analysis to infer the phylogenetic position of these species within the species group, revealing that a large genetic distance separates the two sympatric Chilean species. M. longibrachus is confirmed as the sister species of Muusooctopus eureka from the Falkland Islands; while M. eicomar is a sister species of Muusoctopus yaquinae from the North Pacific, most closely related to the amphi‐Atlantic species Muusoctopus januarii. Molecular divergence times and ancestral distribution analyses suggest that genus Muusoctopus may have originated in the North Atlantic: one lineage dispersed directly southward to the Magellan region and another dispersed southward along the Eastern Pacific to the Southern Ocean and Antarctica. The Muusoctopus species in the Southern Hemisphere have different phylogenetic origins and represent independent invasions of this region.     

Testing an hypothesis of hybrid zone movement for toads in France

Hybrid zone movement may result in substantial unidirectional introgression of selectively neutral material from the local to the advancing species, leaving a genetic footprint. This genetic footprint is represented by a trail of asymmetric tails and displaced cline centres in the wake of the moving hybrid zone. A peak of admixture linkage disequilibrium is predicted to exist ahead of the centre of the moving hybrid zone. We test these predictions of the movement hypothesis in a hybrid zone between common (Bufo bufo) and spined toads (B. spinosus), using 31 nuclear and one mtDNA SNPs along a transect in the northwest of France. Average effective selection in Bufo hybrids is low and clines vary in shape and centre. A weak pattern of asymmetric introgression is inferred from cline discordance of seven nuclear markers. The dominant direction of gene flow is from B. spinosus to B. bufo and is in support of southward movement of the hybrid zone. Conversely, a peak of admixture linkage disequilibrium north of the hybrid zone suggests northward movement. These contrasting results can be explained by reproductive isolation of the B. spinosus and B. bufo gene pools at the southern (B. spinosus) side of the hybrid zone. The joint occurrence of asymmetric introgression and admixture linkage disequilibrium can also be explained by the combination of low dispersal and random genetic drift due to low effective population sizes.     

Population genetic structure in six sympatric and widespread aquatic plants inhabiting diverse lake environments in China

Many aquatic plant species are distributed over large areas and diverse environments with populations interconnected by abiotic and biotic mediators. Here, we examined differences and similarities in the population genetic structure of six sympatric and widespread aquatic plant species. We sampled the aquatic species from six Chinese lakes found on plateaus, plains, and different river systems and analyzed them using inter‐simple sequence repeat (ISSR) markers. Samples originating from each lake tended to cluster together. Of the six species, only Nymphoides peltata and Myriophyllum spicatum could be divided into plateau and plain groups, once Taihu Lake individuals were excluded. Genetic similarities between populations connected by the Yangtze River were not consistently higher than unconnected populations. Populations from Taihu Lake and/or Weishanhu Lake were distant from other lake populations for all species except Potamogeton lucens. The Taihu and Weishanhu populations clustered for Ceratophyllum demersum and Typha latifolia. Hydrophilous C. demersum had the lowest gene flow (Nm = 0.913), whereas the entomophilous Hydrocharis dubia (Nm = 2.084) and N. peltata (Nm = 2.204) had the highest gene flow. The genetic relationships among distant populations of aquatic plants reflect the comprehensive effects of environmental selection pressure and biotic and abiotic connectivity. Differences in environmental factors between plateau and plain lakes and long distance hydrochory have limited importance on aquatic plant genetic structures. Among multiple evolutionary forces, gene flow mediated by birds may play the most important role in the formation of genetic patterns in the six species examined. For example, the close genetic relationship between Taihu Lake and Weishanhu Lake populations, each in different river systems and with different climates, may be related to the migration routes of birds. Differences in gene flow among the six aquatic plants may be attributable to different bird‐transport and the fruit traits of each species.