Phylogeography and population genetics of a widespread cold-adapted ant, Prenolepis imparis

Historical climate fluctuations have left genetic signatures on species and populations across North America. Here, we used phylogenetic and population genetic analyses from 1,402 orthologous sequences of 75 individuals obtained through sequencing of Ultraconserved Elements (UCEs) to identify population genetic structure and historical demographic patterns across the range of a widespread, cold-adapted ant, the winter ant, Prenolepis imparis. We relate the genomic patterns to those expected as a result of in situ diversification, maintained connectivity, or recent migration. We recovered five well-supported, genetically isolated clades across the distribution: 1) a basal lineage located in Florida, 2) populations across the southern United States, 3) populations that span the midwestern and northeastern United States, 4) populations from the western United States, and 5) populations in southwestern Arizona and Mexico. Using Bayesian clustering analysis in STRUCTURE and k-means clustering in ADEGENET, we investigated gene flow between these major genetic clades and did not find evidence of gene flow between clades. We did find evidence of localized structure with migration in the western United States clade. High support for five major geographic lineages and lack of evidence of contemporary gene flow indicate in situ diversification across the species’ range, probably influenced by glacial cycles of the late Quaternary.

Evidence for extensive hybridisation and past introgression events in feather grasses using genome-wide SNP genotyping

Background The proper identification of feather grasses in nature is often limited due to phenotypic variability and high morphological similarity between many species. Among plausible factors influencing this issue are hybridisation and introgression recently detected in the genus. Nonetheless, to date, only a bounded set of taxa have been investigated using integrative taxonomy combining morphological and molecular data. Here, we report the first large-scale study on five feather grass species across several hybrid zones in Russia and Central Asia. In total, 302 specimens were sampled in the field and classified based on the current descriptions of these taxa. They were then genotyped with high density genome-wide markers and measured based on a set of morphological characters to delimitate species and assess levels of hybridisation and introgression. Moreover, we tested species for past introgression and estimated divergence times between them. Results Our findings demonstrated that 250 specimens represent five distinct species: S. baicalensis, S. capillata, S. glareosa, S. grandis and S. krylovii. The remaining 52 individuals provided evidence for extensive hybridisation between S. capillata and S. baicalensis, S. capillata and S. krylovii, S. baicalensis and S. krylovii, as well as to a lesser extent between S. grandis and S. krylovii, S. grandis and S. baicalensis. We detected past reticulation events between S. baicalensis, S. krylovii, S. grandis and inferred that diversification within species S. capillata, S. baicalensis, S. krylovii and S. grandis started ca. 130–96 kya. In addition, the assessment of genetic population structure revealed signs of contemporary gene flow between populations across species from the section Leiostipa, despite significant geographical distances between some of them. Lastly, we concluded that only 5 out of 52 hybrid taxa were properly identified solely based on morphology. Conclusions Our results support the hypothesis that hybridisation is an important mechanism driving evolution in Stipa. As an outcome, this phenomenon complicates identification of hybrid taxa in the field using morphological characters alone. Thus, integrative taxonomy seems to be the only reliable way to properly resolve the phylogenetic issue of Stipa. Moreover, we believe that feather grasses may be a suitable genus to study hybridisation and introgression events in nature.

Changes in selection pressure can facilitate hybridization during biological invasion in a Cuban lizard

Hybridization is among the evolutionary mechanisms most frequently hypothesized to drive the success of invasive species, in part because hybrids are common in invasive populations. One explanation for this pattern is that biological invasions coincide with a change in selection pressures that limit hybridization in the native range. To investigate this possibility, we studied the introduction of the brown anole (Anolis sagrei) in the southeastern United States. We find that native populations are highly genetically structured. In contrast, all invasive populations show evidence of hybridization among native-range lineages. Temporal sampling in the invasive range spanning 15 y showed that invasive genetic structure has stabilized, indicating that large-scale contemporary gene flow is limited among invasive populations and that hybrid ancestry is maintained. Additionally, our results are consistent with hybrid persistence in invasive populations resulting from changes in natural selection that occurred during invasion. Specifically, we identify a large-effect X chromosome locus associated with variation in limb length, a well-known adaptive trait in anoles, and show that this locus is often under selection in the native range, but rarely so in the invasive range. Moreover, we find that the effect size of alleles at this locus on limb length is much reduced in hybrids among divergent lineages, consistent with epistatic interactions. Thus, in the native range, epistasis manifested in hybrids can strengthen extrinsic postmating isolation. Together, our findings show how a change in natural selection can contribute to an increase in hybridization in invasive populations.

Genetic Diversity and Population Structure of the Antarctic Toothfish, Dissostichus mawsoni, Using Mitochondrial and Microsatellite DNA Markers

The Antarctic toothfish, Dissostichus mawsoni, serves as a valuable fishery resource around the Antarctic Continent since 1997, managed by the Commission for the Conservation of Antarctic Marine Living Resources (CCAMLR). Although delineating genetic or stock structure of populations is crucial for improving fishery management of this species, its number of genetic populations and genetic diversity levels remain ambiguous. In the present study, we assessed the population genetic and phylogeographic structure of the Antarctic toothfish across 20 geographic localities spanning from Subareas 88 (88.1, 88.2, and 88.3) to Subareas 58 (58.4 and 58.5) by using mitochondrial DNA (mtDNA) cytochrome oxidase I (COI) and 16S rRNA (16S) sequences and seven nuclear microsatellite loci. MtDNA revealed a low level of polymorphism (h = 0.571, π = 0.0006) with 40 haplotypes in 392 individuals, connected only by 1–5 mutational steps, which is indicative of shallow evolutionary history. Microsatellites showed a range of allelic richness (AR) from 6.328 (88.3 RB3) to 7.274 (88.3 RB6) within populations. Overall genetic diversity was generally higher in Subareas 58 than in Subareas 88, suggesting that effective population size (NE) is larger in Subareas 58. The results of population analyses using microsatellites suggest that the sampled populations are likely to comprise a well-admixed single gene pool (i.e., one genetic stock), perhaps due to high contemporary gene flow occurring during the prolonged larval phase of this fish. However, given weak, but significant microsatellite differentiation found in six population-pairs, the possibility of existence of multiple genetic populations could not be completely excluded. The mtDNA AMOVA suggests a genetic break between the Subareas 88 and 58 groups (FCT = 0.011, P = 0.004). Moreover, mtDNA genetic distances (FST) between populations were proportionally greater as geographic distances increase. The patterns of isolation by distance (IBD) shown only in mtDNA, but not in microsatellites might suggest that population differentiation or divergence processes underwent faster in mtDNA than microsatellites, due to its NE being only one-quarter of nuclear DNA. Temporal stability in the genetic structure of D. mawsoni is also indicated by the results of no genetic differentiation between juveniles and adults. The findings of this study will help to design effective stock management strategies for this valuable fishery resource. We suggest that a long-term genetic monitoring is needed to understand the population structure and dynamics of toothfish in response to ongoing climate changes.

Genomic evidence for the Chinese mountain cat as a wildcat conspecific (Felis silvestris bieti) and its introgression to domestic cats

The Qinghai-Tibet Plateau endemic Chinese mountain cat has a controversial taxonomic status, whether it is a true species or a wildcat (Felis silvestris) subspecies and whether it has contributed to cat (F. s. catus) domestication in East Asia. Here, we sampled F. silvestris lineages across China and sequenced 51 nuclear genomes, 55 mitogenomes, and multilocus regions from 270 modern or museum specimens. Genome-wide analyses classified the Chinese mountain cat as a wildcat conspecific F. s. bieti, which was not involved in cat domestication of China, thus supporting a single domestication origin arising from the African wildcat (F. s. lybica). A complex hybridization scenario including ancient introgression from the Asiatic wildcat (F. s. ornata) to F. s. bieti, and contemporary gene flow between F. s. bieti and sympatric domestic cats that are likely recent Plateau arrivals, raises the prospect of disrupted wildcat genetic integrity, an issue with profound conservation implications.

Depicting the mating system and patterns of contemporary pollen flow in trees of the genus Anadenanthera (Fabaceae)

Anadenanthera (Fabaceae) is endemic to the Neotropics and consists of two tree species: A. colubrina (Vell.) Brenan and A. peregrina (L.) Speg. This study examined the mating system and contemporary gene flow of A. colubrina (Acol) and A. peregrina (Aper) in a highly fragmented area of the Atlantic Forest to provide valuable information that informs conservation strategies. Reproductive adults from forest remnants [nA. colubrina = 30 (2.7 ha), nA. peregrina = 55 (4.0 ha)] and progeny-arrays (nA. colubrina = 322, nA. peregrina = 300) were genotyped for seven nuclear microsatellite markers. Mating system analyses revealed that A. colubrina is a mixed mating species (tm = 0.619) while A. peregrina is a predominantly outcrossing species (tm = 0.905). For both Anadenanthera species, high indices of biparental inbreeding were observed (Acol = 0.159, Aper = 0.216), resulting in low effective pollination neighborhood sizes. Categorical paternity analysis revealed different scales of pollen dispersal distance: the majority of crossings occurring locally (i.e., between nearby trees within the same population), with moderate pollen dispersal coming from outside the forest fragments boundaries (Acolmp = 30%, Apermp = 35%). Nevertheless, pollen immigration from trees outside the populations for both species suggests that the populations are not reproductively isolated. This study highlights the importance of evaluating both mating system and contemporary gene flow for a better understanding of the biology of Anadenanthera species. This information should be considered to ensure the effective conservation and management practices of these plant species.

Genetic Consequences of Forest Fragmentation in a Widespread Forest Bat (Natalus mexicanus, Chiroptera: Natalidae)

Recent historical and anthropogenic changes in the landscape causing habitat fragmentation can disrupt the connectivity of wild populations and pose a threat to the genetic diversity of multiple species. This study investigated the effect of habitat fragmentation on the structure and genetic diversity of the Mexican greater funnel-eared bat (Natalus mexicanus) throughout its distribution range in Mexico, whose natural habitat has decreased dramatically in recent years. Genetic structure and diversity were measured using the HVII hypervariable domain of the mitochondrial control region and ten nuclear microsatellite loci, to analyze historical and contemporary information, respectively. The mitochondrial and nuclear results pointed to a differential genetic structuring, derived mainly from philopatry in females. Our results also showed that genetic diversity was historically high and currently moderate; additionally, the contemporary gene flow between the groups observed was null. These findings confirm that the effects of habitat fragmentation have started to be expressed in populations and that forest loss is already building barriers to contemporary gene flow. The concern is that gene flow is a process essential to ensure that the genetic diversity of N. mexicanus populations (and probably of many other forest species) distributed in Mexico is preserved or increased in the long term by maintaining forest connectivity between locations.

Whole-genome analysis of multiple wood ant population pairs supports similar speciation histories, but different degrees of gene flow, across their European range

The application of demographic history modeling and inference to the study of divergence between species is becoming a cornerstone of speciation genomics. The demographic history is usually reconstructed by analysing a single population from each species, assuming that the divergence history inferred between these populations represents the actual speciation history. However, this assumption is rarely explicitly tested, and it may not be met when species diverge with gene flow. For instance, secondary contact between two species after a range expansion may be confined into a specific geographic region. In this study, we tested to what extent the divergence history inferred from two heterospecific populations would vary depending on their geographic locations, using mound-building red wood ants. The wood ant species Formica polyctena and F. aquilonia have contrasting distributions in Europe and naturally hybridize in Finland. We first performed whole-genome resequencing of 20 individuals sampled in multiple populations across both species ranges. We then reconstructed the divergence histories of distinct heterospecific population pairs using a coalescent-based approach. We found that the analysis of these different population pairs always supported a scenario of divergence with gene flow, suggesting that species divergence started in the Pleistocene (ca. 500 kya) and occurred with continuous asymmetrical gene flow from F. aquilonia to F. polyctena until a recent time, when migration stopped (2-19 kya, depending on the population pair considered). However, we found support for contemporary gene flow in the sympatric population pair from Finland, where hybrids have been described. Overall, our results suggest that divergence histories reconstructed from a few individuals may be reliable and applicable at the species level. Nonetheless, the geographical context of populations chosen to represent their species should be taken into account, as it may affect estimates of migration rates between species when gene flow is heterogeneous across their geographical ranges.

Ocean Currents Drove Genetic Structure of Seven Dominant Mangrove Species Along the Coastlines of Southern China

Mangrove forest ecosystems, which provide important ecological services for marine environments and human activities, are being destroyed worldwide at an alarming rate. The objective of our study was to use molecular data and analytical techniques to separate the effects of historical and contemporary processes on the distribution of mangroves and patterns of population genetic differentiation. Seven mangrove species (Acanthus ilicifolius, Aegiceras corniculatum, Avicennia marina, Bruguiera gymnorrhiza, Kandelia obovata, Lumnitzera racemosa, and Rhizophora stylosa), which are predominant along the coastlines of South China, were genotyped at nuclear (nSSR) and chloroplast (cpSSR) microsatellite markers. We estimated historical and contemporary gene flow, the genetic diversity and population structure of seven mangrove species in China. All of these seven species exhibited few haplotypes, low levels of genetic diversity (HE = 0.160–0.361, with the exception of K. obovata) and high levels of inbreeding (FIS = 0.104–0.637), which may be due to their marginal geographical distribution, human-driven and natural stressors on habitat loss and fragmentation. The distribution patterns of haplotypes and population genetic structures of seven mangrove species in China suggest historical connectivity between populations over a large geographic area. In contrast, significant genetic differentiation [FST = 0.165–0.629 (nSSR); GST = 0.173–0.923 (cpSSR)] indicates that populations of mangroves are isolated from one another with low levels of contemporary gene flow among populations. Our results suggest that populations of mangroves were historically more widely inter-connected and have recently been isolated, likely through a combination of ocean currents and human activities. In addition, genetic admixture in Beibu Gulf populations and populations surrounding Hainan Island and southern mainland China were attributed to asymmetric gene flow along prevailing oceanic currents in China in historical times. Even ocean currents promote genetic exchanges among mangrove populations, which are still unable to offset the effects of natural and anthropogenic fragmentation. The recent isolation and lack of gene flow among populations of mangroves may affect their long-term survival along the coastlines of South China. Our study enhances the understanding of oceanic currents contributing to population connectivity, and the effects of anthropogenic and natural habitat fragmentation on mangroves, thereby informing future conservation efforts and seascape genetics toward mangroves.