Genomic divergence in sympatry indicates strong reproductive barriers and cryptic species within Eucalyptus salubris

Genetic studies are increasingly detecting cryptic taxa that likely represent a significant component of global biodiversity. However, cryptic taxa are often criticized because they are typically detected serendipitously and may not receive the follow‐up study required to verify their geographic or evolutionary limits. Here, we follow‐up a study of Eucalyptus salubris that unexpectedly detected two divergent lineages but was not sampled sufficiently to make clear interpretations. We undertook comprehensive sampling for an independent genomic analysis (3,605 SNPs) to investigate whether the two purported lineages remain discrete genetic entities or if they intergrade throughout the species’ range. We also assessed morphological and ecological traits, and sequenced chloroplast DNA. SNP results showed strong genome‐wide divergence (F _ST = 0.252) between two discrete lineages: one dominated the north and one the southern regions of the species’ range. Within lineages, gene flow was high, with low differentiation (mean F _ST = 0.056) spanning hundreds of kilometers. In the central region, the lineages were interspersed but maintained their genomic distinctiveness: an indirect demonstration of reproductive isolation. Populations of the southern lineage exhibited significantly lower specific leaf area and occurred on soils with lower phosphorus relative to the northern lineage. Finally, two major chloroplast haplotypes were associated with each lineage but were shared between lineages in the central distribution. Together, these results suggest that these lineages have non‐contemporary origins and that ecotypic adaptive processes strengthened their divergence more recently. We conclude that these lineages warrant taxonomic recognition as separate species and provide fascinating insight into eucalypt speciation.     

Shennongjia–Wushan Mountains—One cryptic glacial refugium introduced by the phylogeographical study of the Geometridae moth Ourapteryx szechuana Wehrli

The origin and evolution of biodiversity in the Shennongjia and Wushan Mountains, located in central China, are little known. In this study, we used Ourapteryx szechuana, which is widely distributed in China and northern Nepal, to explore whether these mountains acted as glacial refugia during climate oscillations of the Quaternary. In total, 192 samples of O. szechuana were collected throughout much of the distribution range. Phylogenetic analysis, molecular dating, demographic history reconstructions, and MAXENT were used to investigate the evolutionary history and differentiation mechanisms and predict the potential species distributions during four different periods. The phylogenetic tree and the star‐like median‐joining network strongly supported two reciprocally monophyletic and allopatric lineages. Lineage I was restricted to the Shennongjia and Wushan Mountains. The divergence time of O. szechuana from its sister species O. thibetaria was approximately 1.94 Ma. The differentiation processes of the two intraspecific lineages occurred at approximately 0.47 Ma. The demographic history reconstruction and the ecological niche model suggested that Lineage II experienced an expansion after the LGM (Last Glacial Maximum), whereas Lineage I did not experience any expansion. Our results suggested the Naynayxungla glaciation promoted the divergence of the two lineages by restricting them to different refugia. The valleys of the Shennongjia–Wushan Mountains may have kept stable and warm (thus ice‐free) environments during Quaternary glaciations, allowing this region to act as a glacial refugia. Our studies show that the Shennongjia and Wushan Mountains are likely to be important but little studied glacial refugia for the insect and thus worthy of more attention.

The Naynayxungla Glaciation promoted the divergence of the two lineages by restricting them to different refugia. The Shennongjia and Wushan Mountains are one important and easily neglected refugia for the insect groups and are worthy of more attention.     

DNA metabarcoding quantifies the relative biomass of arthropod taxa in songbird diets: Validation with camera‐recorded diets

Ecological research is often hampered by the inability to quantify animal diets. Diet composition can be tracked through DNA metabarcoding of fecal samples, but whether (complex) diets can be quantitatively determined with metabarcoding is still debated and needs validation using free‐living animals. This study validates that DNA metabarcoding of feces can retrieve actual ingested taxa, and most importantly, that read numbers retrieved from sequencing can also be used to quantify the relative biomass of dietary taxa. Validation was done with the hole‐nesting insectivorous Pied Flycatcher whose diet was quantified using camera footage. Size‐adjusted counts of food items delivered to nestlings were used as a proxy for provided biomass of prey orders and families, and subsequently, nestling feces were assessed through DNA metabarcoding. To explore potential effects of digestion, gizzard and lower intestine samples of freshly collected birds were subjected to DNA metabarcoding. For metabarcoding with Cytochrome Oxidase subunit I (COI), we modified published invertebrate COI primers LCO1490 and HCO1777, which reduced host reads to 0.03%, and amplified Arachnida DNA without significant changing the recovery of other arthropod taxa. DNA metabarcoding retrieved all commonly camera‐recorded taxa. Overall, and in each replicate year (N = 3), the relative scaled biomass of prey taxa and COI read numbers correlated at R = .85 (95CI:0.68–0.94) at order level and at R = .75 (CI:0.67–0.82) at family level. Similarity in arthropod community composition between gizzard and intestines suggested limited digestive bias. This DNA metabarcoding validation demonstrates that quantitative analyses of arthropod diet is possible. We discuss the ecological applications for insectivorous birds.     

Assessing diversity of King Crab Lithodes spp. in the south‐eastern pacific using phylogeny and molecular species delimitation methods

The purpose of this study was to test the hypothesis that the genetic diversity of commercially significant species of King Crabs (Lithodes spp.) along the south‐eastern Pacific (SEP) comprises different independent evolutionary units (IEUs) with spatially isolated distribution. Nine localities from inner and open waters along the SEP Chilean coast (39°S‐55°S) were sampled. We analyzed sequences from 173 individuals for the mitochondrial gene Cytochrome oxidase I (COX‐I), 151 individuals for the Internal Transcribed Spacer 1 (ITS) and 135 for the structural ribosomal RNA (28S). Genetic delimitation was performed through three analytical methods: ABGD, GMYC, and its Bayesian implementation, bGMYC. Bayesian phylogenetic analyses and haplotype networks were also performed. Divergence time between clades was assessed for the COX‐I marker and estimated from known evolutionary rates for this marker in other crustacean species and fossil calibration from other Anomuran species. Delimitation analyses, phylogenetic analyses, and mitochondrial haplotype networks suggested the presence of two deeply divergent mitochondrial lineages of Lithodes in the SEP, referred to as Clade1 and Clade 2. Nuclear markers showed low phylogenetic resolution and therefore were unsuitable for molecular species delimitation. Divergence time analysis of the mitochondrial lineages suggests a separation between Clades of approximately 2.3 Mya. The divergence time obtained suggested that Pliocene glaciations and deglaciations cycles could be involved in hybridization events between Lithodes IEUs at southern tip of South American coasts. The different frequencies of Lithodes haplotypes in inner and open water environments along SEP coasts could be explained by events such as the last glacial maximum or by differences in the adaptation of each clade to different environments. These findings support the necessity of evaluating the taxonomic status of Lithodes individuals found along SEP coasts under an integrative taxonomy approach or through markers with other evolution rates than those already used.     

Pleistocene dynamics of the Eurasian steppe as a driving force of evolution: Phylogenetic history of the genus Capsella (Brassicaceae)

Capsella is a model plant genus of the Brassicaceae closely related to Arabidopsis. To disentangle its biogeographical history and intrageneric phylogenetic relationships, 282 individuals of all five currently recognized Capsella species were genotyped using a restriction digest‐based next‐generation sequencing method. Our analysis retrieved two main lineages within Capsella that split c. one million years ago, with western C. grandiflora and C. rubella forming a sister lineage to the eastern lineage consisting of C. orientalis. The split was attributed to continuous latitudinal displacements of the Eurasian steppe belt to the south during Early Pleistocene glacial cycles. During the interglacial cycles of the Late Pleistocene, hybridization of the two lineages took place in the southwestern East European Plain, leading to the allotetraploid C. bursa‐pastoris. Extant genetic variation within C. orientalis postdated any extensive glacial events. Ecological niche modeling showed that suitable habitat for C. orientalis existed during the Last Glacial Maximum around the north coast of the Black Sea and in southern Kazakhstan. Such a scenario is also supported by population genomic data that uncovered the highest genetic diversity in the south Kazakhstan cluster, suggesting that C. orientalis originated in continental Asia and migrated north‐ and possibly eastwards after the last ice age. Post‐glacial hybridization events between C. bursa‐pastoris and C. grandiflora/rubella in the southwestern East European Plain and the Mediterranean gave rise to C. thracica. Introgression of C. grandiflora/rubella into C. bursa‐pastoris resulted in a new Mediterranean cluster within the already existing Eurasian C. bursa‐pastoris cluster. This study shows that the continuous displacement and disruption of the Eurasian steppe belt during the Pleistocene was the driving force in the evolution of Capsella.     

Phylogeography of sugar kelp: Northern ice‐age refugia in the Gulf of Alaska

Many Northeast (NE) Pacific fishes and invertebrates survived Pleistocene glaciations in northern refugia, but the extent that kelps survived in northern areas is uncertain. Here, we test the hypothesis that populations of sugar kelp (Saccharina latissima) persisted in the Gulf of Alaska during ice‐age maxima when the western margin of the Cordilleran ice sheet covered coastal areas around the NE Pacific Ocean. We estimated genetic diversities within and phylogeographical relationships among 14 populations along 2,800 km in the NE Pacific and Bering Sea with partial sequences of mitochondrial DNA 5′‐cytochrome oxidase subunit I (COI, bp = 624, n = 543), chloroplast DNA ribulose‐1,5‐bisphosphate carboxylase large subunit‐3′ (rbcL, bp = 735, n = 514), and 11 microsatellite loci. Concatenated sequences of rbcL and COI showed moderate levels of within‐population genetic diversity (mean h = 0.200) but substantial differences among populations (Φ_ST = 0.834, p < .0001). Microsatellites showed moderate levels of heterozygosity within populations (mean H _E = 0.391). Kelps in the same organellar lineage tended to cluster together, regardless of geographic origins, as indicated in a principal coordinate analysis (PCoA) of microsatellite genotypes. The PCoA also showed evidence of nuclear hybridizations between co‐occurring organellar lineages. Individual admixture plots with population clusters of K = 2, 6, and 9 showed increasing complexity with considerable historical admixture between some clusters. A time‐calibrated phylogeny placed divergences between rbcL‐COI lineages at 1.4 million years at most. The time frames of mutation in the rbcL‐COI lineages and microsatellite population clusters differed among locations. The existence of ancient lineages in the Gulf of Alaska, moderate levels of genetic diversity, and the absence of departures from neutrality are consistent with northern refugia during multiple Croll‐Milankovitch climate cycles in the Pleistocene Epoch.     

Phylogeography and genetic diversity of the widespread katydid Ducetia japonica (Thunberg, 1815) across China

Habitat fragmentation can lower migration rates and genetic connectivity among remaining populations of native species. Ducetia japonica is one of the most widespread katydids in China, but little is known about its genetic structure and phylogeographic distribution. We combined the five‐prime region of cytochrome c oxidase subunit I (COI‐5P), 11 newly developed microsatellite loci coupled with an ecological niche model (ENM) to examine the genetic diversity and population structure of D. japonica in China and beyond to Laos and Singapore. Both Bayesian inference (BI) and haplotype network methods revealed six mitochondrial COI‐5P lineages. The distribution of COI‐5P haplotypes may not demonstrate significant phylogeographic structure (N _ST > G _ST, p > .05). The STRUCTURE analysis based on microsatellite data also revealed six genetic clusters, but discordant with those obtained from COI‐5P haplotypes. For both COI‐5P and microsatellite data, Mantel tests revealed a significant positive correlation between geographic and genetic distances in mainland China. Bayesian skyline plot (BSP) analyses indicated that the population size of D. japonica''s three major mitochondrial COI‐5P lineages were seemingly not affected by last glacial maximum (LGM, 0.015–0.025 Mya). The ecological niche models showed that the current distribution of D. japonica was similar to the species’ distribution during the LGM period and only slightly extended in northern China. Further phylogeographic studies based on more extensive sampling are needed to identify specific locations of glacial refugia in northern China.     

A multi‐year dormancy strategy in a cabbage beetle population in southeastern China

The life cycle of the cabbage beetle Colaphellus bowringi in southeastern China is complex due to four options for adult development: summer diapause, winter diapause, prolonged diapsuse, and nondiapause. However, detailed information on the multi‐year emergence patterns of diapausing individuals in this beetle has not been documented. In this study, we monitored the adult emergence patterns of diapausing individuals and estimated the influence of the diapause‐inducing temperature and photoperiod on the incidence of prolonged diapause under seminatural conditions for several years. The duration of diapause for adults collected from the vegetable fields in different years varied from several months to 5 years. Approximately 25.9%–29.2% of individuals showed prolonged diapause (emergence more than 1 year after entering diapause) over the 5 years of observation. Furthermore, regardless of insect age, the emergence of diapausing adults from the soil always occurred between mid‐February and March in spring and between late August and mid‐October in autumn, when the host plants were available. The influence of diapause‐inducing temperatures (22, 25, and 28°C) combined with different photoperiods (L:D 12:12 h and L:D 14:10 h) on diapause duration was tested under seminatural conditions. Pairwise comparisons of diapause duration performed by the log‐rank test revealed that the low temperature of 22°C combined with the long photoperiod of L:D 14:10 h induced the longest diapause duration, whereas the low temperature of 22°C combined with the short photoperiod of L:D 12:12 h induced the highest proportion of prolonged diapause. This study indicates that C. bowringi adopts a multi‐year dormancy strategy to survive local environmental conditions and unpredictable risks.     

Mito‐nuclear discordance at a mimicry color transition zone in bumble bee Bombus melanopygus

As hybrid zones exhibit selective patterns of gene flow between otherwise distinct lineages, they can be especially valuable for informing processes of microevolution and speciation. The bumble bee, Bombus melanopygus, displays two distinct color forms generated by Müllerian mimicry: a northern “Rocky Mountain''’ color form with ferruginous mid‐abdominal segments (B. m. melanopygus) and a southern “Pacific''’ form with black mid‐abdominal segments (B. m. edwardsii). These morphs meet in a mimetic transition zone in northern California and southern Oregon that is more narrow and transitions further west than comimetic bumble bee species. To understand the historical formation of this mimicry zone, we assessed color distribution data for B. melanopygus from the last 100 years. We then examined gene flow among the color forms in the transition zone by comparing sequences from mitochondrial COI barcode sequences, color‐controlling loci, and the rest of the nuclear genome. These data support two geographically distinct mitochondrial haplogroups aligned to the ancestrally ferruginous and black forms that meet within the color transition zone. This clustering is also supported by the nuclear genome, which, while showing strong admixture across individuals, distinguishes individuals most by their mitochondrial haplotype, followed by geography. These data suggest the two lineages most likely were historically isolated, acquired fixed color differences, and then came into secondary contact with ongoing gene flow. The transition zone, however, exhibits asymmetries: mitochondrial haplotypes transition further south than color pattern, and both transition over shorter distances in the south. This system thus demonstrates alternative patterns of gene flow that occur in contact zones, presenting another example of mito‐nuclear discordance. Discordant gene flow is inferred to most likely be driven by a combination of mimetic selection, dominance effects, and assortative mating.     

Effect of winter cold duration on spring phenology of the orange tip butterfly,Anthocharis cardamines

The effect of spring temperature on spring phenology is well understood in a wide range of taxa. However, studies on how winter conditions may affect spring phenology are underrepresented. Previous work on Anthocharis cardamines (orange tip butterfly) has shown population‐specific reaction norms of spring development in relation to spring temperature and a speeding up of post‐winter development with longer winter durations. In this experiment, we examined the effects of a greater and ecologically relevant range of winter durations on post‐winter pupal development of A. cardamines of two populations from the United Kingdom and two from Sweden. By analyzing pupal weight loss and metabolic rate, we were able to separate the overall post‐winter pupal development into diapause duration and post‐diapause development. We found differences in the duration of cold needed to break diapause among populations, with the southern UK population requiring a shorter duration than the other populations. We also found that the overall post‐winter pupal development time, following removal from winter cold, was negatively related to cold duration, through a combined effect of cold duration on diapause duration and on post‐diapause development time. Longer cold durations also lead to higher population synchrony in hatching. For current winter durations in the field, the A. cardamines population of southern UK could have a reduced development rate and lower synchrony in emergence because of short winters. With future climate change, this might become an issue also for other populations. Differences in winter conditions in the field among these four populations are large enough to have driven local adaptation of characteristics controlling spring phenology in response to winter duration. The observed phenology of these populations depends on a combination of winter and spring temperatures; thus, both must be taken into account for accurate predictions of phenology.     

Killer whale (Orcinus orca) population dynamics in response to a period of rapid ecosystem change in the eastern North Atlantic

This study investigates survival and abundance of killer whales (Orcinus orca) in Norway in 1988–2019 using capture–recapture models of photo‐identification data. We merged two datasets collected in a restricted fjord system in 1988–2008 (Period 1) with a third, collected after their preferred herring prey shifted its wintering grounds to more exposed coastal waters in 2012–2019 (Period 2), and investigated any differences between these two periods. The resulting dataset, spanning 32 years, comprised 3284 captures of 1236 whales, including 148 individuals seen in both periods. The best‐supported models of survival included the effects of sex and time period, and the presence of transients (whales seen only once). Period 2 had a much larger percentage of transients compared to Period 1 (mean = 30% vs. 5%) and the identification of two groups of whales with different residency patterns revealed heterogeneity in recapture probabilities. This caused estimates of survival rates to be biased downward (females: 0.955 ± 0.027 SE, males: 0.864 ± 0.038 SE) compared to Period 1 (females: 0.998 ± 0.002 SE, males: 0.985 ± 0.009 SE). Accounting for this heterogeneity resulted in estimates of apparent survival close to unity for regularly seen whales in Period 2. A robust design model for Period 2 further supported random temporary emigration at an estimated annual probability of 0.148 (± 0.095 SE). This same model estimated a peak in annual abundance in 2015 at 1061 individuals (95% CI 999–1127), compared to a maximum of 731 (95% CI 505–1059) previously estimated in Period 1, and dropped to 513 (95% CI 488–540) in 2018. Our results indicate variations in the proportion of killer whales present of an undefined population (or populations) in a larger geographical region. Killer whales have adjusted their distribution to shifts in key prey resources, indicating potential to adapt to rapidly changing marine ecosystems.     

Home‐site fidelity and homing behavior of the big‐headed turtle Platysternon megacephalum

Site fidelity refers to the restriction of dispersal distance of an animal and its tendency to return to a stationary site. To our knowledge, the homing ability of freshwater turtles and their fidelity is reportedly very low in Asia. We examined mark–recapture data spanning a 4‐year period in Diaoluoshan National Nature Reserve, Hainan Province, China, to investigate the site fidelity and homing behavior of big‐headed turtles Platysternon megacephalum. A total of 11 big‐headed turtles were captured, and all individuals were used in this mark–recapture study. The site fidelity results showed that the adult big‐headed turtles (n = 4) had a 71.43% recapture rate in the original site after their release at the same site, whereas the juveniles (n = 1) showed lower recapture rates (0%). Moreover, the homing behavior results showed that the adults (n = 5) had an 83.33% homing rate after displacement. Adult big‐headed turtles were able to return to their initial capture sites (home) from 150 to 2,400 m away and precisely to their home sites from either upstream or downstream of their capture sites or even from other streams. However, none of the juveniles (n = 4) returned home, despite only being displaced 25–150 m away. These results indicated that the adult big‐headed turtles showed high fidelity to their home site and strong homing ability. In contrast, the juvenile turtles may show an opposite trend but further research is needed.     

Signature of mid‐Pleistocene lineages in the European silver fir (Abies alba Mill.) at its geographic distribution margin

In a conservation and sustainable management perspective, we identify the ecological, climatic, and demographic factors responsible for the genetic diversity patterns of the European silver fir (Abies alba Mill.) at its southwestern range margin (Pyrenees Mountains, France, Europe). We sampled 45 populations throughout the French Pyrenees and eight neighboring reference populations in the Massif Central, Alps, and Corsica. We genotyped 1,620 individuals at three chloroplast and ten nuclear microsatellite loci. We analyzed within‐ and among‐population genetic diversity using phylogeographic reconstructions, tests of isolation‐by‐distance, Bayesian population structure inference, modeling of demographic scenarios, and regression analyses of genetic variables with current and past environmental variables. Genetic diversity decreased from east to west suggesting isolation‐by‐distance from the Alps to the Pyrenees and from the Eastern to the Western Pyrenees. We identified two Pyrenean lineages that diverged from a third Alpine–Corsica–Massif Central lineage 0.8 to 1.1 M years ago and subsequently formed a secondary contact zone in the Central Pyrenees. Population sizes underwent contrasted changes, with a contraction in the west and an expansion in the east. Glacial climate affected the genetic composition of the populations, with the western genetic cluster only observed in locations corresponding to the coldest past climate and highest elevations. The eastern cluster was observed over a larger range of temperatures and elevations. All demographic events shaping the current spatial structure of genetic diversity took place during the Mid‐Pleistocene Transition, long before the onset of the Holocene. The Western Pyrenees lineage may require additional conservation efforts, whereas the eastern lineage is well protected in in situ gene conservation units. Due to past climate oscillations and the likely emergence of independent refugia, east–west oriented mountain ranges may be important reservoir of genetic diversity in a context of past and ongoing climate change in Europe.     

Extensive sharing of chloroplast haplotypes among East Asian Cerris oaks: The imprints of shared ancestral polymorphism and introgression

Shared ancestral polymorphism and introgression are two main causes of chloroplast DNA (cpDNA) haplotype sharing among closely related angiosperms. In this study, we explored the roles of these two processes in shaping the phylogeographic patterns of East Asian Cerris oaks by examining the geographic distributions of randomly and locally distributed shared haplotypes, which coincide with the expectations of shared ancestry and introgression, respectively. We sequenced 1340 bp of non‐coding cpDNA from Quercus acutissima (n = 418) and Q. chenii (n = 183) and compiled previously published sequence data of Q. variabilis (n = 439). The phylogenetic relationships among haplotypes were examined using a median‐joining network. The geographic patterns of interspecifically shared haplotypes were assessed to test whether nearby populations have a higher degree of interspecific cpDNA sharing than distant ones. We identified a total of 27 haplotypes that were grouped into three non‐species‐specific lineages with overlapping distributions. Ancestral haplotypes were extensively shared and randomly distributed across populations of the three species. Some young haplotypes were locally shared in mountainous areas that may have been shared refugia. The local exchange of cpDNA resulted in an excess of similar haplotypes between nearby populations. Our study demonstrated that the haplotype sharing pattern among East Asian Cerris oaks reflected the imprints of both shared ancestral polymorphism and introgression. This pattern was also associated with the relatively stable climates and complex landscapes in East Asia, which not only allowed the long‐term persistence of ancestral lineages but also connected the survived populations across refugia.     

Revised Timeline and Distribution of the Earliest Diverged Human Maternal Lineages in Southern Africa

The oldest extant human maternal lineages include mitochondrial haplogroups L0d and L0k found in the southern African click-speaking forager peoples broadly classified as Khoesan. Profiling these early mitochondrial lineages allows for better understanding of modern human evolution. In this study, we profile 77 new early-diverged complete mitochondrial genomes and sub-classify another 105 L0d/L0k individuals from southern Africa. We use this data to refine basal phylogenetic divergence, coalescence times and Khoesan prehistory. Our results confirm L0d as the earliest diverged lineage (∼172 kya, 95%CI: 149–199 kya), followed by L0k (∼159 kya, 95%CI: 136–183 kya) and a new lineage we name L0g (∼94 kya, 95%CI: 72–116 kya). We identify two new L0d1 subclades we name L0d1d and L0d1c4/L0d1e, and estimate L0d2 and L0d1 divergence at ∼93 kya (95%CI:76–112 kya). We concur the earliest emerging L0d1’2 sublineage L0d1b (∼49 kya, 95%CI:37–58 kya) is widely distributed across southern Africa. Concomitantly, we find the most recent sublineage L0d2a (∼17 kya, 95%CI:10–27 kya) to be equally common. While we agree that lineages L0d1c and L0k1a are restricted to contemporary inland Khoesan populations, our observed predominance of L0d2a and L0d1a in non-Khoesan populations suggests a once independent coastal Khoesan prehistory. The distribution of early-diverged human maternal lineages within contemporary southern Africans suggests a rich history of human existence prior to any archaeological evidence of migration into the region. For the first time, we provide a genetic-based evidence for significant modern human evolution in southern Africa at the time of the Last Glacial Maximum at between ∼21–17 kya, coinciding with the emergence of major lineages L0d1a, L0d2b, L0d2d and L0d2a.     

Breeding system of diploid sexuals within the Ranunculus auricomus complex and its role in a geographical parthenogenesis scenario

The larger distribution area of asexuals compared with their sexual relatives in geographical parthenogenesis (GP) scenarios has been widely attributed to the advantages of uniparental reproduction and polyploidy. However, potential disadvantages of sexuals due to their breeding system have received little attention so far. Here, we study the breeding system of five narrowly distributed sexual lineages of Ranunculus notabilis s.l. (R. auricomus complex) and its effects on outcrossing, inbreeding, female fitness, and heterozygosity. We performed selfing and intra‐ and interlineage crossings by bagging 481 flowers (59 garden individuals) followed by germination experiments. We compared seed set and germination rates, and related them to genetic distance and genome‐wide heterozygosity (thousands of RADseq loci). Selfings (2.5%) unveiled a significantly lower seed set compared with intra‐ (69.0%) and interlineage crossings (69.5%). Seed set of intra‐ (65%) compared to interpopulation crossings (78%) was significantly lower. In contrast, all treatments showed comparable germination rates (32%–43%). Generalized linear regressions between seed set and genetic distance revealed positive relationships in general and between lineages, and a negative one within lineages. Seed set was the main decisive factor for female fitness. Germination rates were not related to genetic distance at any level, but were positively associated with heterozygosity in interlineage crossings. Experiments confirmed full crossability and predominant outcrossing among sexual R. notabilis s.l. lineages. However, up to 5% (outliers 15%–31%) of seeds were formed by selfing, probably due to semi‐self‐compatibility in a multi‐locus gametophytic SI system. Less seed set in intrapopulation crossings, and higher seed set and germination rates from crossings of genetically more distant and heterozygous lineages (interlineage) indicate negative inbreeding and positive outbreeding effects. In GP scenarios, sexual species with small and/or isolated populations can suffer from decreased female fitness due to their breeding system. This factor, among others, probably limits range expansion of sexuals.     

Toward a global DNA barcode reference library of the intolerant nonbiting midge genus Rheocricotopus Brundin, 1956

Environmental DNA metabarcoding is becoming a predominant tool in biodiversity assessment, as this time‐ and cost‐efficient tactics have the ability to increase monitoring accuracy. As a worldwide distributed genus, Rheocricotopus Brundin, 1956 still does not possess a complete and comprehensive global DNA barcode reference library for biodiversity monitoring. In the present study, we compiled a cytochrome c oxidase subunit 1 (COI) DNA barcode library of Rheocricotopus with 434 barcodes around the world, including 121 newly generated DNA barcodes of 32 morphospecies and 313 public barcodes. Automatic Barcode Gap Discovery (ABGD) was applied on the 434 COI barcodes to provide a comparison between the operational taxonomic units (OTU) number calculated from the Barcode Index Number (BIN) with the “Barcode Gap Analysis” and neighbor‐joining (NJ) tree analysis. Consequently, these 434 COI barcodes were clustered into 78 BINs, including 42 new BINs. ABGD yielded 51 OTUs with a prior intraspecific divergence of Pmax = 7.17%, while NJ tree revealed 52 well‐separated clades. Conservatively, 14 unknown species and one potential synonym were uncovered with reference to COI DNA barcodes. Besides, based on our ecological analysis, we discovered that annual mean temperature and annual precipitation could be considered as key factors associated with distribution of certain members from this genus. Our global DNA barcode reference library of Rheocricotopus provides one fundamental database for accurate species delimitation in Chironomidae taxonomy and facilitates the biodiversity monitoring of aquatic biota.     

Mitochondrial genomes reveal the extinct Hippidion as an outgroup to all living equids

Hippidions were equids with very distinctive anatomical features. They lived in South America 2.5 million years ago (Ma) until their extinction approximately 10 000 years ago. The evolutionary origin of the three known Hippidion morphospecies is still disputed. Based on palaeontological data, Hippidion could have diverged from the lineage leading to modern equids before 10 Ma. In contrast, a much later divergence date, with Hippidion nesting within modern equids, was indicated by partial ancient mitochondrial DNA sequences. Here, we characterized eight Hippidion complete mitochondrial genomes at 3.4–386.3-fold coverage using target-enrichment capture and next-generation sequencing. Our dataset reveals that the two morphospecies sequenced (H. saldiasi and H. principale) formed a monophyletic clade, basal to extant and extinct Equus lineages. This contrasts with previous genetic analyses and supports Hippidion as a distinct genus, in agreement with palaeontological models. We date the Hippidion split from Equus at 5.6–6.5 Ma, suggesting an early divergence in North America prior to the colonization of South America, after the formation of the Panamanian Isthmus 3.5 Ma and the Great American Biotic Interchange.     

Genetic diversity and population dynamic of Ziziphus jujuba var. spinosa (Bunge) Hu ex H. F. Chow in Central China

Phylogeographic research concerning Central China has been rarely conducted. Population genetic and phylogeography of Ziziphus jujuba var. spinosa (also called sour jujube) were investigated to improve our understanding of plant phylogeographic patterns in Central China. Single‐copy nuclear gene markers and complete chloroplast genome data were applied to 328 individuals collected from 21 natural populations of sour jujube in China. Nucleotide variation of sour jujube was relatively high (π = 0.00720, θ _w = 0.00925), which resulted from the mating system and complex population dynamics. Analysis of molecular variation analysis revealed that most of the total variation was attributed to variation within populations, and a high level of genetic differentiation among populations was detected (F _st = 0.197). Relatively low long‐distance dispersal capability and vitality of pollen contributed to high genetic differentiation among populations. Differences in the environmental conditions and long distance among populations further restricted gene flow. Structure clustering analysis uncovered intraspecific divergence between central and marginal populations. Migrate analysis found a high level of gene flow between these two intraspecific groups. Bayesian skyline plot detected population expansion of these two intraspecific groups. Network and phylogeny analysis of chloroplast haplotypes also found intraspecific divergence, and the divergence time was estimated to occur at about 55.86 Ma. Haplotype native to the Loess Plateau was more ancient, and multiple glacial refugia of sour jujube were found to locate at the Loess Plateau, areas adjacent to the Qinling Mountains and Tianmu Mountains. Species distribution model analysis found a typical contraction‐expansion model corresponding to the Quaternary climatic oscillations. In the future, the distribution of sour jujube may shift to high‐latitude areas. This study provides new insights for phylogeographic research of temperate plant species distributed in Central China and sets a solid foundation for the application of the scientific management strategy of Z. jujuba var. spinosa.     

Molecular and iridescent feather reflectance data reveal recent genetic diversification and phenotypic differentiation in a cloud forest hummingbird

The present day distribution and spatial genetic diversity of Mesoamerican biota reflects a long history of responses to habitat change. The hummingbird Lampornis amethystinus is distributed in northern Mesoamerica, with geographically disjunct populations. Based on sampling across the species range using mitochondrial DNA (mtDNA) sequences and nuclear microsatellites jointly analysed with phenotypic and climatic data, we (1) test whether the fragmented distribution is correlated with main evolutionary lineages, (2) assess body size and plumage color differentiation of populations in geographic isolation, and (3) evaluate a set of divergence scenarios and demographic patterns of the hummingbird populations. Analysis of genetic variation revealed four main groups: blue‐throated populations (Sierra Madre del Sur); two groups of amethyst‐throated populations (Trans‐Mexican Volcanic Belt and Sierra Madre Oriental); and populations east of the Isthmus of Tehuantepec (IT) with males showing an amethyst throat. The most basal split is estimated to have originated in the Pleistocene, 2.39–0.57 million years ago (MYA), and corresponded to groups of populations separated by the IT. However, the estimated recent divergence time between blue‐ and amethyst‐throated populations does not correspond to the 2‐MY needed to be in isolation for substantial plumage divergence, likely because structurally iridescent colors are more malleable than others. Results of species distribution modeling and Approximate Bayesian Computation analysis fit a model of lineage divergence west of the Isthmus after the Last Glacial Maximum (LGM), and that the species’ suitable habitat was disjunct during past and current conditions. These results challenge the generality of the contraction/expansion glacial model to cloud forest‐interior species and urges management of cloud forest, a highly vulnerable ecosystem to climate change and currently facing destruction, to prevent further loss of genetic diversity or extinction.