Intraspecific independent evolution of floral spur length in response to local flower visitor size in Japanese Aquilegia in different mountain regions

Geographic differences in floral traits may reflect geographic differences in effective pollinator assemblages. Independent local adaptation to pollinator assemblages in multiple regions would be expected to cause parallel floral trait evolution, although sufficient evidence for this is still lacking. Knowing the intraspecific evolutionary history of floral traits will reveal events that occur in the early stages of trait diversification. In this study, we investigated the relationship between flower spur length and pollinator size in 16 populations of Aquilegia buergeriana var. buergeriana distributed in four mountain regions in the Japanese Alps. We also examined the genetic relationship between yellow‐ and red‐flowered individuals, to see if color differences caused genetic differentiation by pollinator isolation. Genetic relationships among 16 populations were analyzed based on genome‐wide single‐nucleotide polymorphisms. Even among populations within the same mountain region, pollinator size varied widely, and the average spur length of A. buergeriana var. buergeriana in each population was strongly related to the average visitor size of that population. Genetic relatedness between populations was not related to the similarity of spur length between populations; rather, it was related to the geographic proximity of populations in each mountain region. Our results indicate that spur length in each population evolved independently of the population genetic structure but in parallel in response to local flower visitor size in different mountain regions. Further, yellow‐ and red‐flowered individuals of A. buergeriana var. buergeriana were not genetically differentiated. Unlike other Aquilegia species in Europe and America visited by hummingbirds and hawkmoths, the Japanese Aquilegia species is consistently visited by bumblebees. As a result, genetic isolation by flower color may not have occurred.     

Benefits and limitations of a new genome‐based PCR‐RFLP genotyping assay (GB‐RFLP): A SNP‐based detection method for identification of species in extremely young adaptive radiations

High‐throughput DNA sequencing technologies make it possible now to sequence entire genomes relatively easily. Complete genomic information obtained by whole‐genome resequencing (WGS) can aid in identifying and delineating species even if they are extremely young, cryptic, or morphologically difficult to discern and closely related. Yet, for taxonomic or conservation biology purposes, WGS can remain cost‐prohibitive, too time‐consuming, and often constitute a “data overkill.” Rapid and reliable identification of species (and populations) that is also cost‐effective is made possible by species‐specific markers that can be discovered by WGS. Based on WGS data, we designed a PCR restriction fragment length polymorphism (PCR‐RFLP) assay for 19 Neotropical Midas cichlid populations (Amphilophus cf. citrinellus), that includes all 13 described species of this species complex. Our work illustrates that identification of species and populations (i.e., fish from different lakes) can be greatly improved by designing genetic markers using available “high resolution” genomic information. Yet, our work also shows that even in the best‐case scenario, when whole‐genome resequencing information is available, unequivocal assignments remain challenging when species or populations diverged very recently, or gene flow persists. In summary, we provide a comprehensive workflow on how to design RFPL markers based on genome resequencing data, how to test and evaluate their reliability, and discuss the benefits and pitfalls of our approach.     

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.     

RAD sequencing of common whelk,Buccinum undatum,reveals fine‐scale population structuring in Europe and cryptic speciation within the North Atlantic

Buccinum undatum is a subtidal gastropod that exhibits clear spatial variation in several phenotypic shell traits (color, shape, and thickness) across its North Atlantic distribution. Studies of spatial phenotypic variation exist for the species; however, population genetic studies have thus far relied on a limited set of mitochondrial and microsatellite markers. Here, we greatly expand on previous work by characterizing population genetic structure in B. undatum across the North Atlantic from SNP variation obtained by RAD sequencing. There was a high degree of genetic differentiation between Canadian and European populations (Iceland, Faroe Islands, and England) consistent with the divergence of populations in allopatry (F _ST > 0.57 for all pairwise comparisons). In addition, B. undatum populations within Iceland, the Faroe Islands, and England are typified by weak but significant genetic structuring following an isolation‐by‐distance model. Finally, we established a significant correlation between genetic structuring in Iceland and two phenotypic traits: shell shape and color frequency. The works detailed here enhance our understanding of genetic structuring in B. undatum and establish the species as an intriguing model for future genome‐wide association studies.     

Genetic differentiation and phylogeography of rotifer Polyarthra dolichoptera and P. vulgaris populations between Southeastern China and eastern North America: High intercontinental differences

Genetic differentiations and phylogeographical patterns of small organisms may be shaped by spatial isolation, environmental gradients, and gene flow. However, knowledge about genetic differentiation of rotifers at the intercontinental scale is still limited. Polyarthra dolichoptera and P. vulgaris are cosmopolitan rotifers that are tolerant to environmental changes, offering an excellent model to address the research gap. Here, we investigated the populations in Southeastern China and eastern North America and evaluated the phylogeographical patterns from their geographical range sizes, geographic–genetic distance relationships and their responses to spatial‐environmental factors. Using the mitochondrial cytochrome c oxidase subunit I gene as the DNA marker, we analyzed a total of 170 individuals. Our results showed that some putative cryptic species, also known as entities were widely distributed, but most of them were limited to single areas. The divergence of P. dolichoptera and P. vulgaris indicated that gene flow between continents was limited while that within each continent was stronger. Oceanographic barriers do affect the phylogeographic pattern of rotifers in continental waters and serve to maintain genetic diversity in nature. The genetic distance of P. dolichoptera and P. vulgaris populations showed significant positive correlation with geographic distance. This might be due to the combined effects of habitat heterogeneity, long‐distance colonization, and oceanographic barriers. Furthermore, at the intercontinental scale, spatial distance had a stronger influence than environmental variables on the genetic differentiations of both populations. Wind‐ and animal‐mediated transport and even historical events of continental plate tectonics are potential factors for phylogeography of cosmopolitan rotifers.     

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.     

Depicting the phenotypic space of the annual plant Diplotaxis acris in hyperarid deserts

The phenotypic space encompasses the assemblage of trait combinations yielding well‐suited integrated phenotypes. At the population level, understanding the phenotypic space structure requires the quantification of among‐ and within‐population variations in traits and the correlation pattern among them. Here, we studied the phenotypic space of the annual plant Diplotaxis acris occurring in hyperarid deserts. Given the advance of warming and aridity in vast regions occupied by drylands, D. acris can indicate the successful evolutionary trajectory that many other annual plant species may follow in expanding drylands. To this end, we conducted a greenhouse experiment with 176 D. acris individuals from five Saudi populations to quantify the genetic component of variation in architectural and life history traits. We found low among‐population divergence but high among‐individual variation in all traits. In addition, all traits showed a high degree of genetic determination in our study experimental conditions. We did not find significant effects of recruitment and fecundity on fitness. Finally, all architectural traits exhibited a strong correlation pattern among them, whereas for life history traits, only higher seed germination implied earlier flowering. Seed weight appeared to be an important trait in D. acris as individuals with heavier seeds tended to advance flowering and have a more vigorous branching pattern, which led to higher fecundity. Population divergence in D. acris might be constrained by the severity of the hyperarid environment, but populations maintain high among‐individual genetic variation in all traits. Furthermore, D. acris showed phenotypic integration for architectural traits and, to a lesser extent, for life history traits. Overall, we hypothesize that D. acris may be fine‐tuned to its demanding extreme environments. Evolutionary speaking, annual plants facing increasing warming, aridity, and environmental seasonality might modify their phenotypic spaces toward new phenotypic configurations strongly dominated by correlated architectural traits enhancing fecundity and seed‐related traits advancing flowering time.     

ddRAD‐Seq reveals evolutionary insights into population differentiation and the cryptic phylogeography of Hyporhamphus intermedius in Mainland China

Species differentiation and local adaptation in heterogeneous environments have attracted much attention, although little is known about the mechanisms involved. Hyporhamphus intermedius is an anadromous, brackish‐water halfbeak that is widely distributed in coastal areas and hyperdiverse freshwater systems in China, making it an interesting model for research on phylogeography and local adaptation. Here, 156 individuals were sampled at eight sites from heterogeneous aquatic habitats to examine environmental and genetic contributions to phenotypic divergence. Using double‐digest restriction‐site‐associated DNA sequencing (ddRAD‐Seq) in the specimens from the different watersheds, 5498 single nucleotide polymorphisms (SNPs) were found among populations, with obvious population differentiation. We find that present‐day Mainland China populations are structured into distinct genetic clusters stretching from southern and northern ancestries, mirroring geography. Following a transplant event in Plateau Lakes, there were virtually no variations of genetic diversity occurred in two populations, despite the fact two main splits were unveiled in the demographic history. Additionally, dorsal, and anal fin traits varied widely between the southern group and the others, which highlighted previously unrecognized lineages. We then explore genotype–phenotype‐environment associations and predict candidate loci. Subgroup ranges appeared to correspond to geographic regions with heterogeneous hydrological factors, indicating that these features are likely important drivers of diversification. Accordingly, we conclude that genetic and phenotypic polymorphism and a moderate amount of genetic differentiation occurred, which might be ascribed to population subdivision, and the impact of abiotic factors.     

Multigenerational experimental simulation of climate change on an economically important insect pest

Long‐term multigenerational experimental simulations of climate change on insect pests of economically and socially important crops are crucial to anticipate challenges for feeding humanity in the not‐so‐far future. Mexican bean weevil Zabrotes subfasciatus, is a worldwide pest that attacks the common bean Phaseolus vulgaris seeds, in crops and storage. We designed a long term (i.e., over 10 generations), experimental simulation of climate change by increasing temperature and CO₂ air concentration in controlled conditions according to model predictions for 2100. Higher temperature and CO₂ concentrations favored pest''s egg‐to‐adult development survival, even at high female fecundity. It also induced a reduction of fat storage and increase of protein content but did not alter body size. After 10 generations of simulation, genetic adaptation was detected for total lipid content only, however, other traits showed signs of such process. Future experimental designs and methods similar to ours, are key for studying long‐term effects of climate change through multigenerational experimental designs.     

Genomic dynamics of brown trout populations released to a novel environment

Population translocations occur for a variety of reasons, from displacement due to climate change to human‐induced transfers. Such actions have adverse effects on genetic variation and understanding their microevolutionary consequences requires monitoring. Here, we return to an experimental release of brown trout (Salmo trutta) in order to monitor the genomic effects of population translocations. In 1979, fish from each of two genetically (F _ST = 0.16) and ecologically separate populations were simultaneously released, at one point in time, to a lake system previously void of brown trout. Here, whole‐genome sequencing of pooled DNA (Pool‐seq) is used to characterize diversity within and divergence between the introduced populations and fish inhabiting two lakes downstream of the release sites, sampled 30 years later (c. 5 generations). Present results suggest that while extensive hybridization has occurred, the two introduced populations are unequally represented in the lakes downstream of the release sites. One population, which is ecologically resident in its original habitat, mainly contributes to the lake closest to the release site. The other population, migratory in its natal habitat, is genetically more represented in the lake further downstream. Genomic regions putatively under directional selection in the new habitat are identified, where allele frequencies in both established populations are more similar to the introduced population stemming from a resident population than the migratory one. Results suggest that the microevolutionary consequences of population translocations, for example, hybridization and adaptation, can be rapid and that Pool‐seq can be used as an initial tool to monitor genome‐wide effects.     

Multiple introductions and overwintering shape the progressive invasion of Aedes albopictus beyond the Alps

Aedes albopictus originates from Southeast Asia and is considered one of the most invasive species globally. This mosquito is a nuisance and a disease vector of significant public health relevance. In Europe, Ae. albopictus is firmly established and widespread south of the Alps, a mountain range that forms a formidable biogeographic barrier to many organisms. Recent reports of Ae. albopictus north of the Alps raise questions of (1) the origins of its recent invasion, and (2) if this mosquito has established overwintering populations north of the Alps. To answer these questions, we analyzed population genomic data from >4000 genome‐wide SNPs obtained through double‐digest restriction site‐associated DNA sequencing. We collected SNP data from specimens from six sites in Switzerland, north and south of the Alps, and analyzed them together with specimens from other 33 European sites, five from the Americas, and five from its Asian native range. At a global level, we detected four genetic clusters with specimens from Indonesia, Brazil, and Japan as the most differentiated, whereas specimens from Europe, Hong Kong, and USA largely overlapped. Across the Alps, we detected a weak genetic structure and high levels of genetic admixture, supporting a scenario of rapid and human‐aided dispersal along transportation routes. While the genetic pattern suggests frequent re‐introductions into Switzerland from Italian sources, the recovery of a pair of full siblings in two consecutive years in Strasbourg, France, suggests the presence of an overwintering population north of the Alps. The suggestion of overwintering populations of Ae. albopictus north of the Alps and the expansion patterns identified points to an increased risk of further northward expansion and the need for increased surveillance of mosquito populations in Northern Europe.     

Mitochondrial DNA D‐loop sequence analysis reveals high variation and multiple maternal origins of indigenous Tanzanian goat populations

The Small East African (SEA) goat are widely distributed in different agro‐ecological zones of Tanzania. We report the genetic diversity, maternal origin, and phylogenetic relationship among the 12 Tanzanian indigenous goat populations, namely Fipa, Songwe, Tanga, Pwani, Iringa, Newala, Lindi, Gogo, Pare, Maasai, Sukuma, and Ujiji, based on the mitochondrial DNA (mtDNA) D‐loop. High haplotype (H _d = 0.9619–0.9945) and nucleotide (π = 0.0120–0.0162) diversities were observed from a total of 389 haplotypes. The majority of the haplotypes (n = 334) belonged to Haplogroup A which was consistent with the global scenario on the genetic pattern of maternal origin of all goat breeds in the world. Haplogroup G comprised of 45 haplotypes drawn from all populations except the Ujiji goat population while Haplogroup B with 10 haplotypes was dominated by Ujiji goats (41%). Tanzanian goats shared four haplotypes with the Kenyan goats and two with goats from South Africa, Namibia, and Mozambique. There was no sharing of haplotypes observed between individuals from Tanzanian goat populations with individuals from North or West Africa. The indigenous goats in Tanzania have high genetic diversity defined by 389 haplotypes and multiple maternal origins of haplogroup A, B, and G. There is a lot of intermixing and high genetic variation within populations which represent an abundant resource for selective breeding in the different agro‐ecological regions of the country.     

Genetic differentiation and connectivity of morphological types of the broadcast‐spawning coral Galaxea fascicularis in the Nansei Islands,Japan

Population connectivity resulting from larval dispersal is essential for the maintenance or recovery of populations in marine ecosystems, including coral reefs. Studies of species diversity and genetic connectivity within species are essential for the conservation of corals and coral reef ecosystems. We analyzed mitochondrial DNA sequence types and microsatellite genotypes of the broadcast‐spawning coral, Galaxea fascicularis, from four regions in the subtropical Nansei Islands in the northwestern Pacific Ocean. Two types (soft and hard types) of nematocyst morphology are known in G. fascicularis and are significantly correlated with the length of a mitochondrial DNA noncoding sequence (soft type: mt‐L; hard type: mt‐S type). Using microsatellites, significant genetic differentiation was detected between the mitochondrial DNA sequence types in all regions. We also found a third genetic cluster (mt‐L+), and this unexpected type may be a cryptic species of Galaxea. High clonal diversity was detected in both mt‐L and mt‐S types. Significant genetic differentiation, which was found among regions within a given type (F _ST = 0.009–0.024, all Ps ≤ 0.005 in mt‐L; 0.009–0.032, all Ps ≤ 0.01 in mt‐S), may result from the shorter larval development than in other broadcast‐spawning corals, such as the genus Acropora. Nevertheless, intraspecific genetic diversity and connectivity have been maintained, and with both sexual and asexual reproduction, this species appears to have a potential for the recovery of populations after disturbance.     

Genetic structure and life history are key factors in species distribution models of spiny lobsters

AimWe incorporated genetic structure and life history phase in species distribution models (SDMs) constructed for a widespread spiny lobster, to reveal local adaptations specific to individual subspecies and predict future range shifts under the RCP 8.5 climate change scenario.LocationIndo‐West Pacific.MethodsMaxEnt was used to construct present‐day SDMs for the spiny lobster Panulirus homarus and individually for the three genetically distinct subspecies of which it comprises. SDMs incorporated both sea surface and benthic (seafloor) climate layers to recreate discrete influences of these habitats during the drifting larval and benthic juvenile and adult life history phases. Principle component analysis (PCA) was used to infer environmental variables to which individual subspecies were adapted. SDM projections of present‐day habitat suitability were compared with predictions for the year 2,100, under the RCP 8.5 climate change scenario.ResultsIn the PCA, salinity best explained P. h. megasculptus habitat suitability, compared with current velocity in P. h. rubellus and sea surface temperature in P. h. homarus. Drifting and benthic life history phases were adapted to different combinations of sea surface and benthic environmental variables considered. Highly suitable habitats for benthic phases were spatially enveloped within more extensive sea surface habitats suitable for drifting larvae. SDMs predicted that present‐day highly suitable habitats for P. homarus will decrease by the year 2,100.Main conclusionsIncorporating genetic structure in SDMs showed that individual spiny lobster subspecies had unique adaptations, which could not be resolved in species‐level models. The use of sea surface and benthic climate layers revealed the relative importance of environmental variables during drifting and benthic life history phases. SDMs that included genetic structure and life history were more informative in predictive models of climate change effects.     

Historic range dynamics in Kaiser's mountain newt (Neurergus kaiseri): Insights from phylogeographic analyses and species distribution modeling

Vulnerable Kaiser''s mountain newt, Neurergus kaiseri, is endemic to highland streams, springs, and pools of the southwestern Zagros mountain, Iran. The present study aimed to use an integration of phylogeographical and species distribution modeling (SDM) approaches to provide new insights into the evolutionary history of the species throughout Quaternary climate oscillations. The phylogeographical analysis was followed by analyzing two mitochondrial DNA (mt‐DNA) markers including 127 control region (D‐loop) and 72 NADH dehydrogenase 2 (ND2) sequences from 15 populations in the entire species range that were obtained from GenBank. Potential recent and past distribution (the Last Glacial Maximum, LGM, 21 Kya and the Mid‐Holocene, 6 Kya) reconstructed by ensemble SDM using nine algorithms with CCSM4, MIROC‐ESM, and MPI‐ESM‐P models. N. kaiseri displayed two distinct lineages in the northern and southern regions that diverged in the Early‐Pleistocene. The demographics analysis showed signs of a slight increase in effective population size for both northern and southern populations in the Mid‐Pleistocene. Biogeography analysis showed that both vicariance and dispersal events played an important role in the formation of recent species distribution of N. kaiseri. Based on SDM projection onto paleoclimatic data, N. kaiseri displayed a scenario of past range expansion that followed by postglacial contraction. The models showed that the distribution range of the species may have shifted to a lower altitude during LGM while with amelioration of climatic during Mid‐Holocene to recent conditions caused the species to shift to the higher altitude. The findings of the current study support the hypothesis that the Zagros mountains​ may be acting as climatic refugia and play an important role in the protection of isolated populations during climate oscillations.     

Epigenetic variation associated with responses to different habitats in the context of genetic divergence in Phragmites australis

The mechanisms underlying heritable phenotypic divergence associated with adaptation in response to environmental stresses may involve both genetic and epigenetic variations. Several prior studies have revealed even higher levels of epigenetic variation than genetic variation. However, few population‐level studies have explored the effects of epigenetic variation on species with high levels of genetic diversity distributed across different habitats. Using AFLP and methylation‐sensitive AFLP markers, we tested the hypothesis that epigenetic variation may contribute to differences in plants occupying different habitats when genetic variation alone cannot fully explain adaptation. As a cosmopolitan invasive species, Phragmites australis (common reed) together with high genetic diversity and remarkable adaptability has been suggested as a model for responses to global change and indicators of environmental fluctuations. We found high levels of genetic and epigenetic diversity and significant genetic/epigenetic structure within each of 12 studied populations sampled from four natural habitats of P. australis. Possible adaptive epigenetic variation was suggested by significant correlations between DNA methylation‐based epigenetic differentiation and adaptive genetic divergence in populations across the habitats. Meanwhile, various AMOVAs indicated that some epigenetic differences may respond to various local habitats. A partial Mantel test was used to tease out the correlations between genetic/epigenetic variation and habitat after controlling for the correlation between genetic and epigenetic variations. We found that epigenetic diversity was affected mostly by soil nutrient availability, suggesting that at least some epigenetic differentiation occurred independently of genetic variation. We also found stronger correlations between epigenetic variation and phenotypic traits than between genetic variation and such traits. Overall, our findings indicate that genetically based differentiation correlates with heterogeneous habitats, while epigenetic variation plays an important role in ecological differentiation in natural populations of P. australis. In addition, our results suggest that when assessing global change responses of plant species, intraspecific variation needs to be considered.     

A resurrection study reveals limited evolution of phenology in response to recent climate change across the geographic range of the scarlet monkeyflower

Premise of the studyAs global climate change alters drought regimes, rapid evolution of traits that facilitate adaptation to drought can rescue populations in decline. The evolution of phenological advancement can allow plant populations to escape drought, but evolutionary responses in phenology can vary across a species'' range due to differences in drought intensity and standing genetic variation.Methods Mimulus cardinalis, a perennial herb spanning a broad climatic gradient, recently experienced a period of record drought. Here, we used a resurrection study comparing flowering time and stem height at first flower of pre‐drought ancestors and post‐drought descendants from northern‐edge, central, and southern‐edge populations in a common environment to examine the evolution of drought escape across the latitudinal range.Key resultsContrary to the hypothesis of the evolution of advanced phenology in response to recent drought, flowering time did not advance between ancestors and descendants in any population, though storage condition and maternal effects could have impacted these results. Stem height was positively correlated with flowering time, such that plants that flowered earlier were shorter at first flower. This correlation could constrain the evolution of earlier flowering time if selection favors flowering early at a large size.ConclusionsThese findings suggest that rapid evolution of phenology will not rescue these populations from recent climate change. Future work is needed to examine the potential for the evolution of alternative drought strategies and phenotypic plasticity to buffer M. cardinalis populations from changing climate.     

Low genome‐wide divergence between two lizard populations with high adaptive phenotypic differentiation

Usually, adaptive phenotypic differentiation is paralleled by genetic divergence between locally adapted populations. However, adaptation can also happen in a scenario of nonsignificant genetic divergence due to intense gene flow and/or recent differentiation. While this phenomenon is rarely published, findings on incipient ecologically driven divergence or isolation by adaptation are relatively common, which could confound our understanding about the frequency at which they actually occur in nature. Here, we explore genome‐wide traces of divergence between two populations of the lacertid lizard Psammodromus algirus separated by a 600 m elevational gradient. These populations seem to be differentially adapted to their environments despite showing low levels of genetic differentiation (according to previously studies of mtDNA and microsatellite data). We performed a search for outliers (i.e., loci subject to selection) trying to identify specific loci with F_ST statistics significantly higher than those expected on the basis of overall, genome‐wide estimates of genetic divergence. We find that local phenotypic adaptation (in terms of a wide diversity of characters) was not accompanied by genome‐wide differentiation, even when we maximized the chances of unveiling such differentiation at particular loci with F_ST‐based outlier detection tests. Instead, our analyses confirmed the lack of genome‐wide differentiation on the basis of more than 70,000 SNPs, which is concordant with a scenario of local adaptation without isolation by environment. Our results add evidence to previous studies in which local adaptation does not lead to any kind of isolation (or early stages of ecological speciation), but maintains phenotypic divergence despite the lack of a differentiated genomic background.     

Uncovering the genetic diversity of yams (Dioscorea spp.) in China by combining phenotypic trait and molecular marker analyses

Yam is an important edible tuber and root plant worldwide; China as one of the native places of yams has many diverse local resources. The goal of this study was to clarify the genetic diversity of the commonly cultivated yam landraces and the genetic relationship between the main yam species in China. In this study, 26 phenotypic traits of 112 yam accessions from 21 provinces in China were evaluated, and 24 simple sequence repeat (SSR) and 29 sequence‐related amplified polymorphism (SRAP) markers were used for the genetic diversity analysis. Phenotypic traits revealed that Dioscorea opposita had the highest genetic diversity, followed by D. alata, D. persimilis, D. fordii, and D. esculenta. Among the 26 phenotypic traits, the Shannon diversity indexes of leaf shape, petiole color, and stem color were high, and the range in the variation of tuber‐related traits in the underground part was higher than that in the aboveground part. All accessions were divided into six groups by phenotypic trait clustering, which was also supported by principal component analysis (PCA). Molecular marker analysis showed that SSR and SRAP markers had good amplification effects and could effectively and accurately evaluate the genetic variation of yam. The unweighted pair‐group method with arithmetic means analysis based on SSR‐SRAP marker data showed that the 112 accessions were also divided into six groups, similar to the phenotypic trait results. The results of PCA and population structure analysis based on SSR‐SRAP data also produced similar results. In addition, the analysis of the origin and genetic relationship of yam indicated that the species D. opposita may have originated from China. These results demonstrate the genetic diversity and distinctness among the widely cultivated species of Chinese yam and provide a theoretical reference for the classification, breeding, germplasm innovation, utilization, and variety protection of Chinese yam resources.     

High genetic and epigenetic variation of transposable elements: Potential drivers to rapid adaptive evolution for the noxious invasive weed Mikania micrantha

Why invasive species can rapidly adapt to novel environments is a puzzling question known as the genetic paradox of invasive species. This paradox is explainable in terms of transposable elements (TEs) activity, which are theorized to be powerful mutational forces to create genetic variation. Mikania micrantha, a noxious invasive weed, in this sense provides an excellent opportunity to test the explanation. The genetic and epigenetic variation of 21 invasive populations of M. micrantha in southern China have been examined by using transposon display (TD) and transposon methylation display (TMD) techniques to survey 12 TE superfamilies. Our results showed that M. micrantha populations maintained an almost equally high level of TE‐based genetic and epigenetic variation and they have been differentiated into subpopulations genetically and epigenetically. A similar positive spatial genetic and epigenetic structure pattern was observed within 300 m. Six and seven TE superfamilies presented significant genetic and epigenetic isolation by distance (IBD) pattern. In total, 59 genetic and 86 epigenetic adaptive TE loci were identified. Of them, 51 genetic and 44 epigenetic loci were found to correlate with 25 environmental variables (including precipitation, temperature, vegetation coverage, and soil metals). Twenty‐five transposon‐inserted genes were sequenced and homology‐based annotated, which are found to be involved in a variety of molecular and cellular functions. Our research consolidates the importance of TE‐associated genetic and epigenetic variation in the rapid adaptation and invasion of M. micrantha.