Population genetic structure of the bank vole Myodes glareolus within its glacial refugium in peninsular Italy

It is now well established that Southern European peninsulas have been major glacial refugia for temperate species during Pleistocene climatic oscillations. However, substantial environmental changes occurred also within these peninsulas throughout the Pleistocene, raising questions about the role and interplay of various microevolutionary processes in shaping patterns of intraspecific diversity within these areas. Here, we investigate the patterns of genetic variation in the bank vole Myodes glareolus within the Italian peninsula. By using a panel of 13 microsatellite loci, we found more intraspecific variation than expected based on previous assessments. Indeed, both Bayesian and ordination‐based clustering analyses of variation recovered five main geographic/genetic clusters along the peninsula, with three clusters geographically restricted to the southern portion of the peninsula. This clustering is supported by previous evidences of some morphological distinctiveness among these populations. This pattern can be explained by a refugia‐within‐refugia scenario, with the occurrence of multiple sub‐refugia for the bank vole within the Italian peninsula, likely promoted by the major palaeo‐environmental changes which affected forested habitats within this area during the Pleistocene. Moreover, our results support a scenario whereby the high levels of intraspecific diversity observed within major Pleistocene refugia are better explained by dynamic microevolutionary processes occurred within these areas, rather than by long‐term demographic stability of refugial population. Finally, the narrow and isolated distribution of some of the identified lineages suggests the need for future assessments of their conservation and taxonomic status.     

Borrelia burgdorferi sensu lato infection pressure shapes innate immune gene evolution in natural rodent populations across Europe

Although parasite-mediated selection is assumed to be the main driver of immune gene evolution, empirical evidence that parasites induce allele frequency changes at host immune genes in time and/or space remains scarce. Here, I show that the frequency of a protective gene variant of the innate immune receptor Toll-like receptor 2 in natural bank vole (Myodes glareolus) populations is positively associated with the strength of Borrelia burgdorferi sensu lato infection risk across the European continent. Thereby, this study provides rare evidence for the role of spatially variable infection pressures in moulding the vertebrate immune system.     

Toll‐like receptor 9 protects non‐immune cells from stress by modulating mitochondrial ATP synthesis through the inhibition of SERCA2

Toll‐like receptor 9 (TLR9) has a key role in the recognition of pathogen DNA in the context of infection and cellular DNA that is released from damaged cells. Pro‐inflammatory TLR9 signalling pathways in immune cells have been well investigated, but we have recently discovered an alternative pathway in which TLR9 temporarily reduces energy substrates to induce cellular protection from stress in cardiomyocytes and neurons. However, the mechanism by which TLR9 stimulation reduces energy substrates remained unknown. Here, we identify the calcium‐transporting ATPase, SERCA2 (also known as Atp2a2), as a key molecule for the alternative TLR9 signalling pathway. TLR9 stimulation reduces SERCA2 activity, modulating Ca²⁺ handling between the SR/ER and mitochondria, which leads to a decrease in mitochondrial ATP levels and the activation of cellular protective machinery. These findings reveal how distinct innate responses can be elicited in immune and non‐immune cells—including cardiomyocytes—using the same ligand‐receptor system.     

Oxidative damage and antioxidant defense are assay and tissue‐dependent both in captive and wild‐caught bank voles (Myodes glareolus) before and after reproduction

Reproduction is costly and life‐history theory predicts that current parental investment will result in lower survival or decreased future reproduction. The physiological mechanisms mediating the link between reproduction and survival are still under debate and elevated oxidative damage during reproduction has been proposed as a plausible candidate. Previous studies of oxidative stress during reproduction in animals under natural conditions have been restricted to analyses of blood. Herein, we measured the level of oxidative damage to lipids (tiobarbituric‐acid‐reactive substances) and proteins (carbonyls) in the liver, kidneys, heart and skeletal muscles in free‐living bank vole females from spring and autumn generations, before and after reproduction. Antioxidant defense in the liver and kidneys was also determined. We expected oxidative damage to tissues and hypothesized that the damage would be more uniform between tissues in wild animals compared to those breeding under laboratory conditions. Considering all combinations of markers/tissues/generations, oxidative damage in females did not differ before and after reproduction in 12 comparisons, was lower after reproduction in three comparisons, and was higher after breeding in one comparison. The total glutathione was significantly increased after reproduction only in the liver of the autumn generation and there was no change in catalase activity. Our results confirm—for the first time in the field—previous observations from laboratory studies that there is no simple link between oxidative stress and reproduction and that patterns depend on the tissue and marker being studied. Overall, however, our study does not support the hypothesis that the cost of reproduction in bank voles is mediated by oxidative stress in these tissues.     

Flyway structure in the circumpolar greater white‐fronted goose

Dispersal and migratory behavior are influential factors in determining how genetic diversity is distributed across the landscape. In migratory species, genetic structure can be promoted via several mechanisms including fidelity to distinct migratory routes. Particularly within North America, waterfowl management units have been delineated according to distinct longitudinal migratory flyways supported by banding data and other direct evidence. The greater white‐fronted goose (Anser albifrons) is a migratory waterfowl species with a largely circumpolar distribution consisting of up to six subspecies roughly corresponding to phenotypic variation. We examined the rangewide population genetic structure of greater white‐fronted geese using mtDNA control region sequence data and microsatellite loci from 23 locales across North America and Eurasia. We found significant differentiation in mtDNA between sampling locales with flyway delineation explaining a significant portion of the observed genetic variation (~12%). This is concordant with band recovery data which shows little interflyway or intercontinental movements. However, microsatellite loci revealed little genetic structure suggesting a panmictic population across most of the Arctic. As with many high‐latitude species, Beringia appears to have played a role in the diversification of this species. A common Beringian origin of North America and Asian populations and a recent divergence could at least partly explain the general lack of structure at nuclear markers. Further, our results do not provide strong support for the various taxonomic proposals for this species except for supporting the distinctness of two isolated breeding populations within Cook Inlet, Alaska (A. a. elgasi) and Greenland (A. a. flavirostris), consistent with their subspecies status.     

Drift,not selection,shapes toll‐like receptor variation among oceanic island populations

Understanding the relative role of different evolutionary forces in shaping the level and distribution of functional genetic diversity among natural populations is a key issue in evolutionary and conservation biology. To do so accurately genetic data must be analysed in conjunction with an unambiguous understanding of the historical processes that have acted upon the populations. Here, we focused on diversity at toll‐like receptor (TLR) loci, which play a key role in the vertebrate innate immune system and, therefore, are expected to be under pathogen‐mediated selection. We assessed TLR variation within and among 13 island populations (grouped into three archipelagos) of Berthelot's pipit, Anthus berthelotii, for which detailed population history has previously been ascertained. We also compared the variation observed with that found in its widespread sister species, the tawny pipit, Anthus campestris. We found strong evidence for positive selection at specific codons in TLR1LA, TLR3 and TLR4. Despite this, we found that at the allele frequency level, demographic history has played the major role in shaping patterns of TLR variation in Berthelot's pipit. Levels of diversity and differentiation within and across archipelagos at all TLR loci corresponded very closely with neutral microsatellite variation and with the severity of the bottlenecks that occurred during colonization. Our study shows that despite the importance of TLRs in combating pathogens, demography can be the main driver of immune gene variation within and across populations, resulting in patterns of functional variation that can persist over evolutionary timescales.     

High genomic diversity in the bank vole at the northern apex of a range expansion: The role of multiple colonizations and end‐glacial refugia

The history of repeated northern glacial cycling and southern climatic stability has long dominated explanations for how genetic diversity is distributed within temperate species in Eurasia and North America. However, growing evidence indicates the importance of cryptic refugia for northern colonization dynamics. An important geographic region to assess this is Fennoscandia, where recolonization at the end of the last glaciation was restricted to specific routes and temporal windows. We used genomic data to analyse genetic diversity and colonization history of the bank vole (Myodes glareolus) throughout Europe (>800 samples) with Fennoscandia as the northern apex. We inferred that bank voles colonized Fennoscandia multiple times by two different routes; with three separate colonizations via a southern land‐bridge route deriving from a “Carpathian” glacial refugium and one via a north‐eastern route from an “Eastern” glacial refugium near the Ural Mountains. Clustering of genome‐wide SNPs revealed high diversity in Fennoscandia, with eight genomic clusters: three of Carpathian origin and five Eastern. Time estimates revealed that the first of the Carpathian colonizations occurred before the Younger Dryas (YD), meaning that the first colonists survived the YD in Fennoscandia. Results also indicated that introgression between bank and northern red‐backed voles (Myodes rutilus) took place in Fennoscandia just after end‐glacial colonization. Therefore, multiple colonizations from the same and different cryptic refugia, temporal and spatial separations and interspecific introgression have shaped bank vole genetic variability in Fennoscandia. Together, these processes drive high genetic diversity at the apex of the northern expansion in this emerging model species.     

Short‐term variations in gene flow related to cyclic density fluctuations in the common vole

In highly fluctuating populations with complex social systems, genetic patterns are likely to vary in space and time due to demographic and behavioural processes. Cyclic rodents are extreme examples of demographically instable populations that often exhibit strong social organization. In such populations, kin structure and spacing behaviour may vary with density fluctuations and impact both the composition and spatial structure of genetic diversity. In this study, we analysed the multiannual genetic structure of a cyclic rodent, Microtus arvalis, using a sample of 875 individuals trapped over three complete cycles (from 1999 to 2007) and genotyped at 10 microsatellite loci. We tested the predictions that genetic diversity and gene flow intensity vary with density fluctuations. We found evidences for both spatial scale‐dependant variations in genetic diversity and higher gene flow during high density. Moreover, investigation of sex‐specific relatedness patterns revealed that, although dispersal is biased toward males in this species, distances moved by both sexes were lengthened during high density. Altogether, these results suggest that an increase in migration with density allows to restore the local loss of genetic diversity occurring during low density. We then postulate that this change in migration results from local competition, which enhances female colonization of empty spaces and male dispersal among colonies.     

Contrasting patterns of polymorphism and selection in bacterial‐sensing toll‐like receptor 4 in two house mouse subspecies

Detailed investigation of variation in genes involved in pathogen recognition is crucial for understanding co‐evolutionary processes between parasites and their hosts. Triggering immediate innate response to invading microbes, Toll‐like receptors (TLRs) belong presently among the best‐studied receptors of vertebrate immunity. TLRs exhibit remarkable interspecific variation and also intraspecific polymorphism is well documented. In humans and laboratory mice, several studies have recently shown that single amino acid substitution may significantly alter receptor function. Unfortunately, data concerning polymorphism in free‐living species are still surprisingly scarce. In this study, we analyzed the polymorphism of Toll‐like receptor 4 (Tlr4) over the Palearctic range of house mouse (Mus musculus). Our results reveal contrasting evolutionary patterns between the two recently (0.5 million years ago) diverged house mouse subspecies: M. m. domesticus (Mmd) and M. m. musculus (Mmm). Comparison with cytochrome b indicates strong directional selection in Mmd Tlr4. Throughout the whole Mmd western Palaearctic region, a single variant of the ligand‐binding region is spread, encoded mainly by one dominant haplotype (71% of Mmd). In contrast, Tlr4 in Mmm is much more polymorphic with several haplotypes at intermediate frequencies. Moreover, we also found clear signals of recombination between two principal haplogroups in Mmm, and we identified eight sites under positive selection in our dataset. Our results suggest that observed differences in Tlr4 diversity may be attributed to contrasting parasite‐mediated selection acting in the two subspecies.     

Nucleotide variability and haplotype heterogeneity at the porcine fat mass and obesity‐associated (FTO) gene

A large number of studies have confirmed that variants within the fat mass and obesity‐associated (FTO) gene are associated with higher obesity risk in humans. We and others have shown that FTO polymorphisms are associated with fat deposition and related traits in several pig populations, thus confirming the role of this gene in fatness across species. However, some differences observed in different pig populations may be derived, at least in part, from genetic heterogeneity at this locus. Here, we characterise the nucleotide variability and haplotype diversity of the porcine FTO gene in breeds having different predispositions to fat deposition traits. We resequenced 4749 bp of coding and non‐coding regions of the porcine FTO gene in 44 pigs of eight different breeds and identified 27 single nucleotide polymorphisms (SNPs) and four insertions/deletions. A positive Tajima's D‐value (P < 0.10) obtained in Italian Duroc pigs may be compatible with putative balancing selection. From the sequenced pig panel, 20 haplotypes were inferred, some of which clustered according to the breed of origin (Meishan and Italian Duroc). Genetic heterogeneity at this locus could complicate the dissection of the effects of this gene on fat deposition and production traits in pigs. This situation resembles, to some extent, what has been reported in humans, thus making the study of the porcine FTO gene variability especially interesting, as it could be used as a model to understand the complex and elusive role of this gene in human obesity.     

Contrasting genetic responses to population fragmentation in a coevolving fig and fig wasp across a mainland–island archipelago

Interacting species of pollinator–host systems, especially the obligate ones, are sensitive to habitat fragmentation, due to the nature of mutual dependence. Comparative studies of genetic structure can provide insights into how habitat fragmentation contributes to patterns of genetic divergence among populations of the interacting species. In this study, we used microsatellites to analyse genetic variation in Chinese populations of a typical mutualistic system – Ficus pumila and its obligate pollinator Wiebesia sp. 1 – in a naturally fragmented landscape. The plants and wasps showed discordant patterns of genetic variation and geographical divergence. There was no significant positive relationship in genetic diversity between the two species. Significant isolation‐by‐distance (IBD) patterns occurred across the populations of F. pumila and Wiebesia sp. 1 as whole, and IBD also occurred among island populations of the wasps, but not the plants. However, there was no significant positive relationship in genetic differentiation between them. The pollinator populations had significantly lower genetic variation in small habitat patches than in larger patches, and three island pollinator populations showed evidence of a recent bottleneck event. No effects of patch size or genetic bottlenecks were evident in the plant populations. Collectively, the results indicate that, in more fragmented habitats, the pollinators, but not the plants, have experienced reduced genetic variation. The contrasting patterns have multiple potential causes, including differences in longevity and hence number of generations experiencing fragmentation; different dispersal patterns, with the host's genes dispersed as seeds as well as a result of pollen dispersal via the pollinator; asymmetrical responses to fluctuations in partner populations; and co‐existence of a rare second pollinating wasp on some islands. These results indicate that strongly interdependent species may respond in markedly different ways to habitat fragmentation.     

Overexpression of the leucine‐rich receptor‐like kinase gene LRK2 increases drought tolerance and tiller number in rice

Drought represents a key limiting factor of global crop distribution. Receptor‐like kinases play major roles in plant development and defence responses against stresses such as drought. In this study, LRK2, which encodes a leucine‐rich receptor‐like kinase, was cloned and characterized and found to be localized on the plasma membrane in rice. Promoter–GUS analysis revealed strong expression in tiller buds, roots, nodes and anthers. Transgenic plants overexpressing LRK2 exhibited enhanced tolerance to drought stress due to an increased number of lateral roots compared with the wild type at the vegetative stage. Moreover, ectopic expression of LRK2 seedlings resulted in increased tiller development. Yeast two‐hybrid screening and bimolecular fluorescence complementation (BiFC) indicated a possible interaction between LRK2 and elongation factor 1 alpha (OsEF1A) in vitro. These results suggest that LRK2 functions as a positive regulator of the drought stress response and tiller development via increased branch development in rice. These findings will aid our understanding of branch regulation in other grasses and support improvements in rice genetics.     

Local adaptation,genetic divergence,and experimental selection in a foundation grass across the US Great Plains’ climate gradient

Many prior studies have uncovered evidence for local adaptation using reciprocal transplant experiments. However, these studies are rarely conducted for a long enough time to observe succession and competitive dynamics in a community context, limiting inferences for long‐lived species. Furthermore, the genetic basis of local adaptation and genetic associations with climate has rarely been identified. Here, we report on a long‐term (6‐year) experiment conducted under natural conditions focused on Andropogon gerardii, the dominant grass of the North American Great Plains tallgrass ecosystem. We focus on this foundation grass that comprises 80% of tallgrass prairie biomass and is widely used in 20,000 km² of restoration. Specifically, we asked the following questions: (a) Whether ecotypes are locally adapted to regional climate in realistic ecological communities. (b) Does adaptive genetic variation underpin divergent phenotypes across the climate gradient? (c) Is there evidence of local adaptation if the plants are exposed to competition among ecotypes in mixed ecotype plots? Finally, (d) are local adaptation and genetic divergence related to climate? Reciprocal gardens were planted with 3 regional ecotypes (originating from dry, mesic, wet climate sources) of Andropogon gerardii across a precipitation gradient (500–1,200 mm/year) in the US Great Plains. We demonstrate local adaptation and differentiation of ecotypes in wet and dry environments. Surprisingly, the apparent generalist mesic ecotype performed comparably under all rainfall conditions. Ecotype performance was underpinned by differences in neutral diversity and candidate genes corroborating strong differences among ecotypes. Ecotype differentiation was related to climate, primarily rainfall. Without long‐term studies, wrong conclusions would have been reached based on the first two years. Further, restoring prairies with climate‐matched ecotypes is critical to future ecology, conservation, and sustainability under climate change.     

Microsatellite‐based analysis of genetic structure and gene flow of Mythimna separata (Walker) (Lepidoptera: Noctuidae) in China

The oriental armyworm, Mythimna separata, is a serious agricultural pest in China. Seasonal and roundtrip migration has recently led to sudden, localized outbreaks and crop losses. To evaluate genetic differentiation between populations in eastern and western China and elucidate gene flow, the genetic structure of 20 natural populations from nine provinces was examined using seven microsatellite markers. The results indicated high genetic diversity. However, little to moderate (0 < F_ST < 0.15) genetic differentiation was detected, and there was no correlation between genetic distance and geographical distance. Bayesian clustering analysis identified three groups whereas discriminant analysis of principal components identified ten clusters that were considered as two clear‐cut clusters and one admixed group. Gene flow occurred frequently in most population pairs, and an asymmetrical migration rate was detected in several pairwise population comparisons. The bottleneck test showed that few populations had experienced recent bottlenecks. Correspondingly, large‐scale and long‐distance migration of M. separata has caused low genetic differentiation and frequent gene exchange. Our findings are important for studying genetic evolution and help to improve predictions of M. separata outbreaks in China.     

Range‐wide analysis of genetic structure in a widespread,highly mobile species (Odocoileus hemionus) reveals the importance of historical biogeography

Highly mobile species that thrive in a wide range of habitats are expected to show little genetic differentiation across their range. A limited but growing number of studies have revealed that patterns of broad‐scale genetic differentiation can and do emerge in vagile, continuously distributed species. However, these patterns are complex and often shaped by both historical and ecological factors. Comprehensive surveys of genetic variation at a broad scale and at high resolution are useful for detecting cryptic spatial genetic structure and for investigating the relative roles of historical and ecological processes in structuring widespread, highly mobile species. In this study, we analysed 10 microsatellite loci from over 1900 samples collected across the full range of mule deer (Odocoileus hemionus), one of the most widely distributed and abundant of all large mammal species in North America. Through both individual‐ and population‐based analyses, we found evidence for three main genetic lineages, one corresponding to the ‘mule deer’ morphological type and two to the ‘black‐tailed deer’ type. Historical biogeographic events likely are the primary drivers of genetic divergence in this species; boundaries of the three lineages correspond well with predictions based on Pleistocene glacial cycles, and substructure within each lineage demonstrates island vicariance. However, across large geographic areas, including the entire mule deer lineage, we found that genetic variation fit an isolation‐by‐distance pattern rather than discrete clusters. A lack of genetic structure across wide geographic areas of the continental west indicates that ecological processes have not resulted in restrictions to gene flow sufficient for spatial genetic structure to emerge. Our results have important implications for our understanding of evolutionary mechanisms of divergence, as well as for taxonomy, conservation and management.     

Single residues in the LRR domain of the wheat PM3A immune receptor can control the strength and the spectrum of the immune response

The development of improved plant nucleotide‐binding, leucine‐rich repeat (LRR) immune receptors (NLRs) has mostly been based on random mutagenesis or on structural information available for specific receptors complexed with the recognized pathogen effector. Here, we use a targeted mutagenesis approach based on the natural diversity of the Pm3 powdery mildew resistance alleles present in different wheat (Triticum aestivum) genotypes. In order to understand the functional importance of the amino acid polymorphisms between the active immune receptor PM3A and the inactive ancestral variant PM3CS, we exchanged polymorphic regions and residues in the LRR domain of PM3A with the corresponding segments of PM3CS. These novel variants were functionally tested for recognition of the corresponding AVRPM3^A2/F2 avirulence protein in Nicotiana benthamiana. We identified polymorphic residues in four regions of PM3A that enhance the immune response, but also residues that reduce it or result in complete loss of function. We found that the identified critical residues in PM3A modify its activation threshold towards different protein variants of AVRPM3^A2/F2. PM3A variants with a lowered threshold gave a stronger overall response and gained an extended recognition spectrum. One of these variant proteins with a single amino acid change was stably transformed into wheat, where it conferred race‐specific resistance to mildew. This is a proof of concept that improved PM3A variants with an enlarged recognition spectrum can be engineered based on natural diversity by exchanging single or multiple residues that modulate resistance function.     

Copy number variations of CBF genes at the Fr‐A2 locus are essential components of winter hardiness in wheat

Winter hardiness is important for the adaptation of wheat to the harsh winter conditions in temperate regions and is thus also an important breeding goal. Here, we employed a panel of 407 European winter wheat cultivars to dissect the genetic architecture of winter hardiness. We show that copy number variation (CNV) of CBF (C‐repeat Binding Factor) genes at the Fr‐A2 locus is the essential component for winter survival, with CBF‐A14 CNV being the most likely causal polymorphism, accounting for 24.3% of the genotypic variance. Genome‐wide association mapping identified several markers in the Fr‐A2 chromosomal region, which even after accounting for the effects of CBF‐A14 copy number explained approximately 15% of the genotypic variance. This suggests that additional, as yet undiscovered, polymorphisms are present at the Fr‐A2 locus. Furthermore, CNV of Vrn‐A1 explained an additional 3.0% of the genotypic variance. The allele frequencies of all loci associated with winter hardiness were found to show geographic patterns consistent with their role in adaptation. Collectively, our results from the candidate gene analysis, association mapping and genome‐wide prediction show that winter hardiness in wheat is a quantitative trait, but with a major contribution of the Fr‐A2 locus.     

Climate‐driven range shifts explain the distribution of extant gene pools and predict future loss of unique lineages in a marine brown alga

The climate‐driven dynamics of species ranges is a critical research question in evolutionary ecology. We ask whether present intraspecific diversity is determined by the imprint of past climate. This is an ongoing debate requiring interdisciplinary examination of population genetic pools and persistence patterns across global ranges. Previously, contrasting inferences and predictions have resulted from distinct genomic coverage and/or geographical information. We aim to describe and explain the causes of geographical contrasts in genetic diversity and their consequences for the future baseline of the global genetic pool, by comparing present geographical distribution of genetic diversity and differentiation with predictive species distribution modelling (SDM) during past extremes, present time and future climate scenarios for a brown alga, Fucus vesiculosus. SDM showed that both atmospheric and oceanic variables shape the global distribution of intertidal species, revealing regions of persistence, extinction and expansion during glacial and postglacial periods. These explained the distribution and structure of present genetic diversity, consisting of differentiated genetic pools with maximal diversity in areas of long‐term persistence. Most of the present species range comprises postglacial expansion zones and, in contrast to highly dispersive marine organisms, expansions involved only local fronts, leaving distinct genetic pools at rear edges. Besides unravelling a complex phylogeographical history and showing congruence between genetic diversity and persistent distribution zones, supporting the hypothesis of niche conservatism, range shifts and loss of unique genetic diversity at the rear edge were predicted for future climate scenarios, impoverishing the global gene pool.