Genetic and Chemical Diversity in Perovskia abrotanoides Kar. (Lamiaceae) Populations Based on ISSRs Markers and Essential Oils Profile

Genetic and the essential oil composition variability among twelve Perovskia abrotanoides populations (PAbPs) growing wild in Iran were assessed by ISSR markers, GC‐FID and GC/MS, respectively. Nine selected ISSR primers produced 119 discernible bands, of them 96 (80.7%) being polymorphic. Genetic similarity values among populations ranged between 0.07 and 0.79 which indicated a high level of genetic variation. Polymorphic information content, resolving power and marker index generated by ISSR primers were, 0.31, 6.14, and 3.32, respectively. UPGMA grouped PAbPs into four main clusters. Altogether, 38 chemical compounds were identified in the oils, and a relatively high variation in their contents was found. Camphor (11.9 – 27.5%), 1,8‐cineole (11.3 – 21.3%), α‐bisabolol (0.0 – 13.1%), α‐pinene (5.9 – 10.8%), and δ‐3‐carene (0.1 – 10.5%) were the major compounds. Oxygenated monoterpenes (32.1 – 35.8%) and monoterpene hydrocarbons (25.7 – 30.4%) were the main groups of compounds in the oils studied. Cluster analysis and principal‐component analysis were used to characterize the samples according to oil components. Four main chemotypes were found to be Chemotype I (camphor/1,8‐cineol), Chemotype II (1,8‐cineole/camphor), Chemotype III (camphor/1,8‐cineol/α‐bisabolol), and Chemotype IV (camphor/δ‐3‐carene/α‐bisabolol). The information, provided here on P. abrotanoides populations, will be useful to introduce this plant into agricultural systems.     

Influence of gene position on the expression divergence of oxidative response genes in intraspecific yeast

Phenotypic variation can arise from differences in the protein coding sequence and in the regulatory elements. However, little is known about the contribution of regulatory difference to the expression divergence, especially the cis and trans regulatory variation to the expression divergence in intraspecific populations. In this study, we used two different yeast strains, BY4743 and RM11‐1a/α, to study the regulatory variation to the expression divergence between BY and RM under oxidative stress condition. Our results indicated that the expression divergence of BY and RM is mainly due to trans regulatory variations under both normal and oxidative stress conditions. However, cis regulatory variation seems to play a very important role in oxidative stress response in yeast because 36% of genes showed an increase in cis regulatory variation effect compared with 13% of genes that showed an increase in trans regulatory variation effect after oxidative stress. Our data also indicated that genes located on the longer arm of the chromosomes are more susceptible to cis variation effect under oxidative stress than genes on the shorter arm of the chromosomes.     

Outlier analyses to test for local adaptation to breeding grounds in a migratory arctic seabird

Investigating the extent (or the existence) of local adaptation is crucial to understanding how populations adapt. When experiments or fitness measurements are difficult or impossible to perform in natural populations, genomic techniques allow us to investigate local adaptation through the comparison of allele frequencies and outlier loci along environmental clines. The thick‐billed murre (Uria lomvia) is a highly philopatric colonial arctic seabird that occupies a significant environmental gradient, shows marked phenotypic differences among colonies, and has large effective population sizes. To test whether thick‐billed murres from five colonies along the eastern Canadian Arctic coast show genomic signatures of local adaptation to their breeding grounds, we analyzed geographic variation in genome‐wide markers mapped to a newly assembled thick‐billed murre reference genome. We used outlier analyses to detect loci putatively under selection, and clustering analyses to investigate patterns of differentiation based on 2220 genomewide single nucleotide polymorphisms (SNPs) and 137 outlier SNPs. We found no evidence of population structure among colonies using all loci but found population structure based on outliers only, where birds from the two northernmost colonies (Minarets and Prince Leopold) grouped with birds from the southernmost colony (Gannet), and birds from Coats and Akpatok were distinct from all other colonies. Although results from our analyses did not support local adaptation along the latitudinal cline of breeding colonies, outlier loci grouped birds from different colonies according to their non‐breeding distributions, suggesting that outliers may be informative about adaptation and/or demographic connectivity associated with their migration patterns or nonbreeding grounds.     

Quantitative genetic variance in experimental fly populations evolving with or without environmental heterogeneity

Heterogeneous environments are typically expected to maintain more genetic variation in fitness within populations than homogeneous environments. However, the accuracy of this claim depends on the form of heterogeneity as well as the genetic basis of fitness traits and how similar the assay environment is to the environment of past selection. Here, we measure quantitative genetic (QG) variance for three traits important for fitness using replicated experimental populations of Drosophila melanogaster evolving under four selective regimes: constant salt‐enriched medium (Salt), constant cadmium‐enriched medium (Cad), and two heterogeneous regimes that vary either temporally (Temp) or spatially (Spatial). As theory predicts, we found that Spatial populations tend to harbor more genetic variation than Temp populations or those maintained in a constant environment that is the same as the assay environment. Contrary to expectation, Salt populations tend to have more genetic variation than Cad populations in both assay environments. We discuss the patterns for QG variances across regimes in relation to previously reported data on genome‐wide sequence diversity. For some traits, the QG patterns are similar to the diversity patterns of ecological selected SNPs, whereas the QG patterns for some other traits resembled that of neutral SNPs.     

Polymorphisms in a desaturase 2 ortholog associate with cuticular hydrocarbon and male mating success variation in a natural population of Drosophila serrata

Elucidating the nature of genetic variation underlying both sexually selected traits and the fitness components of sexual selection is essential to understanding the broader consequences of sexual selection as an evolutionary process. To date, there have been relatively few attempts to connect the genetic variance in sexually selected traits with segregating DNA sequence polymorphisms. We set out to address this in a well‐characterized sexual selection system – the cuticular hydrocarbons (CHCs) of Drosophila serrata – using an indirect association study design that allowed simultaneous estimation of the genetic variance in CHCs, sexual fitness and single nucleotide polymorphism (SNP) effects in an outbred population. We cloned and sequenced an ortholog of the D. melanogaster desaturase 2 gene, previously shown to affect CHC biosynthesis in D. melanogaster, and associated 36 SNPs with minor allele frequencies > 0.02 with variance in CHCs and sexual fitness. Three SNPs had significant multivariate associations with CHC phenotype (q‐value < 0.05). At these loci, minor alleles had multivariate effects on CHCs that were weakly associated with the multivariate direction of sexual selection operating on these traits. Two of these SNPs had pleiotropic associations with male mating success, suggesting these variants may underlie responses to sexual selection due to this locus. There were 15 significant male mating success associations (q‐value < 0.1), and interestingly, we detected a nonrandom pattern in the relationship between allele frequency and direction of effect on male mating success. The minor‐frequency allele usually reduced male mating success, suggesting a positive association between male mating success and total fitness at this locus.     

Resistance gene enrichment sequencing (RenSeq) enables reannotation of the NB‐LRR gene family from sequenced plant genomes and rapid mapping of resistance loci in segregating populations

RenSeq is a NB‐LRR (nucleotide binding‐site leucine‐rich repeat) gene‐targeted, Resistance gene enrichment and sequencing method that enables discovery and annotation of pathogen resistance gene family members in plant genome sequences. We successfully applied RenSeq to the sequenced potato Solanum tuberosum clone DM, and increased the number of identified NB‐LRRs from 438 to 755. The majority of these identified R gene loci reside in poorly or previously unannotated regions of the genome. Sequence and positional details on the 12 chromosomes have been established for 704 NB‐LRRs and can be accessed through a genome browser that we provide. We compared these NB‐LRR genes and the corresponding oligonucleotide baits with the highest sequence similarity and demonstrated that ~80% sequence identity is sufficient for enrichment. Analysis of the sequenced tomato S. lycopersicum ‘Heinz 1706’ extended the NB‐LRR complement to 394 loci. We further describe a methodology that applies RenSeq to rapidly identify molecular markers that co‐segregate with a pathogen resistance trait of interest. In two independent segregating populations involving the wild Solanum species S. berthaultii (Rpi‐ber2) and S. ruiz‐ceballosii (Rpi‐rzc1), we were able to apply RenSeq successfully to identify markers that co‐segregate with resistance towards the late blight pathogen Phytophthora infestans. These SNP identification workflows were designed as easy‐to‐adapt Galaxy pipelines.     

Changes in concentrations of circulating heat‐shock proteins in House Finches in response to different environmental stressors

Heat‐shock proteins (HSP) are molecular chaperones that play key roles in the maintenance of cellular homeostasis under variable environmental conditions. Although HSP are frequently used in studies of wild vertebrates as indicators of stress, no one to date has assessed responses of HSP60, HSP70, and HSP90 in the same species to different environmental stressors. We studied changes in the circulating concentrations of HSP60, HSP70, and HSP90 in wild‐caught House Finches (Haemorhous mexicanus) in response to multiple and sequential stressful environments, including high temperatures, transportation, and pathogen exposure. House Finches sampled during a period of low‐environmental stress with moderate ambient temperatures had low levels of HSP60 and modest levels of HSP70 and HSP90 compared to birds sampled during a presumably more stressful period with high temperatures. After exposure to high‐ambient temperatures, transportation in a vehicle, and exposure to Mycoplasma gallisepticum, captive finches were found to have increasingly higher levels of HSP60. HSP70 tended to rise in response to each stressor, but to drop in the weeks between stress challenges. HSP90 levels increased significantly only in response to pathogen challenge. Our observations suggest that HSP60 and HSP70 are indices of a range of stressors in House Finches, whereas HSP90 primarily reflects health state.     

Finding needles in a genomic haystack: targeted capture identifies clear signatures of selection in a nonmodel plant species

Teasing apart neutral and adaptive genomic processes and identifying loci that are targets of selection can be difficult, particularly for nonmodel species that lack a reference genome. However, identifying such loci and the factors driving selection have the potential to greatly assist conservation and restoration practices, especially for the management of species in the face of contemporary and future climate change. Here, we focus on assessing adaptive genomic variation within a nonmodel plant species, the narrow‐leaf hopbush (Dodonaea viscosa ssp. angustissima), commonly used for restoration in Australia. We used a hybrid‐capture target enrichment approach to selectively sequence 970 genes across 17 populations along a latitudinal gradient from 30°S to 36°S. We analysed 8462 single‐nucleotide polymorphisms (SNPs) for F_ST outliers as well as associations with environmental variables. Using three different methods, we found 55 SNPs with significant correlations to temperature and water availability, and 38 SNPs to elevation. Genes containing SNPs identified as under environmental selection were diverse, including aquaporin and abscisic acid genes, as well as genes with ontologies relating to responses to environmental stressors such as water deprivation and salt stress. Redundancy analysis demonstrated that only a small proportion of the total genetic variance was explained by environmental variables. We demonstrate that selection has led to clines in allele frequencies in a number of functional genes, including those linked to leaf shape and stomatal variation, which have been previously observed to vary along the sampled environmental cline. Using our approach, gene regions subject to environmental selection can be readily identified for nonmodel organisms.     

Heterogeneous genetic invasions of three insecticide resistance mutations in Indo‐Pacific populations of Aedes aegypti (L.)

Nations throughout the Indo‐Pacific region use pyrethroid insecticides to control Aedes aegypti, the mosquito vector of dengue, often without knowledge of pyrethroid resistance status of the pest or origin of resistance. Two mutations (V1016G + F1534C) in the sodium channel gene (Vssc) of Ae. aegypti modify ion channel function and cause target‐site resistance to pyrethroid insecticides, with a third mutation (S989P) having a potential additive effect. Of 27 possible genotypes involving these mutations, some allelic combinations are never seen whereas others predominate. Here, five allelic combinations common in Ae. aegypti from the Indo‐Pacific region are described and their geographical distributions investigated using genome‐wide SNP markers. We tested the hypothesis that resistance allele combinations evolved de novo in populations versus the alternative that dispersal of Ae. aegypti between populations facilitated genetic invasions of allele combinations. We used latent factor mixed‐models to detect SNPs throughout the genome that showed structuring in line with resistance allele combinations and compared variation at SNPs within the Vssc gene with genome‐wide variation. Mixed‐models detected an array of SNPs linked to resistance allele combinations, all located within or in close proximity to the Vssc gene. Variation at SNPs within the Vssc gene was structured by resistance profile, whereas genome‐wide SNPs were structured by population. These results demonstrate that alleles near to resistance mutations have been transferred between populations via linked selection. This indicates that genetic invasions have contributed to the widespread occurrence of Vssc allele combinations in Ae. aegypti in the Indo‐Pacific region, pointing to undocumented mosquito invasions between countries.     

Climate‐associated genetic variation in Fagus sylvatica and potential responses to climate change in the French Alps

Local adaptation patterns have been found in many plants and animals, highlighting the genetic heterogeneity of species along their range of distribution. In the next decades, global warming is predicted to induce a change in the selective pressures that drive this adaptive variation, forcing a reshuffling of the underlying adaptive allele distributions. For species with low dispersion capacity and long generation time such as trees, the rapidity of the change could impede the migration of beneficial alleles and lower their capacity to track the changing environment. Identifying the main selective pressures driving the adaptive genetic variation is thus necessary when investigating species capacity to respond to global warming. In this study, we investigate the adaptive landscape of Fagus sylvatica along a gradient of populations in the French Alps. Using a double‐digest restriction‐site‐associated DNA (ddRAD) sequencing approach, we identified 7,000 SNPs from 570 individuals across 36 different sites. A redundancy analysis (RDA)‐derived method allowed us to identify several SNPs that were strongly associated with climatic gradients; moreover, we defined the primary selective gradients along the natural populations of F. sylvatica in the Alps. Strong effects of elevation and humidity, which contrast north‐western and south‐eastern site, were found and were believed to be important drivers of genetic adaptation. Finally, simulations of future genetic landscapes that used these findings allowed identifying populations at risk for F. sylvatica in the Alps, which could be helpful for future management plans.     

Water availability drives signatures of local adaptation in whitebark pine (Pinus albicaulis Engelm.) across fine spatial scales of the Lake Tahoe Basin,USA

Patterns of local adaptation at fine spatial scales are central to understanding how evolution proceeds, and are essential to the effective management of economically and ecologically important forest tree species. Here, we employ single and multilocus analyses of genetic data (n = 116 231 SNPs) to describe signatures of fine‐scale adaptation within eight whitebark pine (Pinus albicaulis Engelm.) populations across the local extent of the environmentally heterogeneous Lake Tahoe Basin, USA. We show that despite highly shared genetic variation (F_ST = 0.0069), there is strong evidence for adaptation to the rain shadow experienced across the eastern Sierra Nevada. Specifically, we build upon evidence from a common garden study and find that allele frequencies of loci associated with four phenotypes (mean = 236 SNPs), 18 environmental variables (mean = 99 SNPs), and those detected through genetic differentiation (n = 110 SNPs) exhibit significantly higher signals of selection (covariance of allele frequencies) than could be expected to arise, given the data. We also provide evidence that this covariance tracks environmental measures related to soil water availability through subtle allele frequency shifts across populations. Our results replicate empirical support for theoretical expectations of local adaptation for populations exhibiting strong gene flow and high selective pressures and suggest that ongoing adaptation of many P. albicaulis populations within the Lake Tahoe Basin will not be constrained by the lack of genetic variation. Even so, some populations exhibit low levels of heritability for the traits presumed to be related to fitness. These instances could be used to prioritize management to maintain adaptive potential. Overall, we suggest that established practices regarding whitebark pine conservation be maintained, with the additional context of fine‐scale adaptation.     

In the presence of population structure: From genomics to candidate genes underlying local adaptation

Understanding the genomic signatures, genes, and traits underlying local adaptation of organisms to heterogeneous environments is of central importance to the field evolutionary biology. To identify loci underlying local adaptation, models that combine allelic and environmental variation while controlling for the effects of population structure have emerged as the method of choice. Despite being evaluated in simulation studies, there has not been a thorough investigation of empirical evidence supporting local adaptation across these alleles. To evaluate these methods, we use 875 Arabidopsis thaliana Eurasian accessions and two mixed models (GEMMA and LFMM) to identify candidate SNPs underlying local adaptation to climate. Subsequently, to assess evidence of local adaptation and function among significant SNPs, we examine allele frequency differentiation and recent selection across Eurasian populations, in addition to their distribution along quantitative trait loci (QTL) explaining fitness variation between Italy and Sweden populations and cis‐regulatory/nonsynonymous sites showing significant selective constraint. Our results indicate that significant LFMM/GEMMA SNPs show low allele frequency differentiation and linkage disequilibrium across locally adapted Italy and Sweden populations, in addition to a poor association with fitness QTL peaks (highest logarithm of odds score). Furthermore, when examining derived allele frequencies across the Eurasian range, we find that these SNPs are enriched in low‐frequency variants that show very large climatic differentiation but low levels of linkage disequilibrium. These results suggest that their enrichment along putative functional sites most likely represents deleterious variation that is independent of local adaptation. Among all the genomic signatures examined, only SNPs showing high absolute allele frequency differentiation (AFD) and linkage disequilibrium (LD) between Italy and Sweden populations showed a strong association with fitness QTL peaks and were enriched along selectively constrained cis‐regulatory/nonsynonymous sites. Using these SNPs, we find strong evidence linking flowering time, freezing tolerance, and the abscisic‐acid pathway to local adaptation.     

Development and preliminary evaluation of a genomewide single nucleotide polymorphisms resource generated by RAD‐seq for the small yellow croaker (Larimichthys polyactis)

Recent advances in high‐throughput sequencing technologies have offered the possibility to generate genomewide sequence data to delineate previously unidentified genetic structure, obtain more accurate estimates of demographic parameters and to evaluate potential adaptive divergence. Here, we identified 27 556 single nucleotide polymorphisms for the small yellow croaker (Larimichthys polyactis) using restriction‐site‐associated DNA (RAD) sequencing of 24 individuals from two populations. Significant sources of genetic variation were identified, with an average nucleotide diversity (π) of 0.00105 ± 0.000425 across individuals, and long‐term effective population size was thus estimated to range between 26 172 and 261 716. According to the results, no differentiation between the two populations was detected based on the SNP data set of top quality score per contig or neutral loci. However, the two analysed populations were highly differentiated based on SNP data set of both top F_ST value per contig and the outlier SNPs. Moreover, local adaptation was highlighted by an F_ST‐based outlier tests implemented in LOSITAN and a total of 538 potentially locally selected SNPs were identified. blast2go annotation of contigs containing the outlier SNPs yielded hits for 37 (66%) of 56 significant blastx matches. Candidate genes for local adaptation constituted a wide array of biological functions, including cellular response to oxidative stress, actin filament binding, ion transmembrane transport and synapse assembly. The generated SNP resources in this study provided a valuable tool for future population genetics and genomics studies of L. polyactis.     

Association between host's genetic diversity and parasite burden in damselflies

Recent research indicates that low genetic variation in individuals can increase susceptibility to parasite infection, yet evidence from natural invertebrate populations remains scarce. Here, we studied the relationship between genetic heterozygosity, measured as AFLP‐based inbreeding coefficient f_AFLP, and gregarine parasite burden from eleven damselfly, Calopteryx splendens, populations. We found that in the studied populations, 5–92% of males were parasitized by endoparasitic gregarines (Apicomplexa: Actinocephalidae). Number of parasites ranged from none to 47 parasites per male, and parasites were highly aggregated in a few hosts. Mean individual f_AFLP did not differ between populations. Moreover, we found a positive association between individual's inbreeding coefficient and parasite burden. In other words, the more homozygous the individual, the more parasites it harbours. Thus, parasites are likely to pose strong selection pressure against inbreeding and homozygosity. Our results support the heterozygosity‐fitness correlation hypothesis, which suggests the importance of heterozygosity for an individual's pathogen resistance.     

Whole‐genome resequencing uncovers molecular signatures of natural and sexual selection in wild bighorn sheep

The identification of genes influencing fitness is central to our understanding of the genetic basis of adaptation and how it shapes phenotypic variation in wild populations. Here, we used whole‐genome resequencing of wild Rocky Mountain bighorn sheep (Ovis canadensis) to >50‐fold coverage to identify 2.8 million single nucleotide polymorphisms (SNPs) and genomic regions bearing signatures of directional selection (i.e. selective sweeps). A comparison of SNP diversity between the X chromosome and the autosomes indicated that bighorn males had a dramatically reduced long‐term effective population size compared to females. This probably reflects a long history of intense sexual selection mediated by male–male competition for mates. Selective sweep scans based on heterozygosity and nucleotide diversity revealed evidence for a selective sweep shared across multiple populations at RXFP2, a gene that strongly affects horn size in domestic ungulates. The massive horns carried by bighorn rams appear to have evolved in part via strong positive selection at RXFP2. We identified evidence for selection within individual populations at genes affecting early body growth and cellular response to hypoxia; however, these must be interpreted more cautiously as genetic drift is strong within local populations and may have caused false positives. These results represent a rare example of strong genomic signatures of selection identified at genes with known function in wild populations of a nonmodel species. Our results also showcase the value of reference genome assemblies from agricultural or model species for studies of the genomic basis of adaptation in closely related wild taxa.     

Variations in Essential Oil Composition and Antioxidant Activity in Perovskia abrotanoides Kar. Collected from Different Regions in Iran

Essential oil (EO) composition, phenolic content, and antioxidant activity were investigated in 17 P. abrotanoides populations collected from different geographical regions in Iran. The highest (3.61%) and lowest (1.25%) essential oil yields were measured in populations from Semnan Province (PSESM₂) and PISKS from Isfahan Province, respectively. GC/MS analysis identified camphor (4.05 – 35.94%), 1,8‐cineole (7.15 – 24.34%), borneol (0 – 21.75%), and α‐pinene (2.05 – 10.33%) as the main constituents of Perovskia essential oil. Cluster analysis classified the studied populations into four different groups: (I) high camphene, (II) high camphor/1,8‐cineole, (III) high borneol/δ‐3‐carene, and (IV) high α‐cadinol/trans‐caryophyllene. The highest flavonoid and phenolic contents were detected in PISAK from Isfahan Province (4.09 ± 0.05 mgQE/gDW, 58.51 ± 1.63 mgGAE/gDW) and PKRGS from Khorasan Province (3.80 ± 0.002 mgQE/gDW, 66.86 ± 0.002 mgGAE/gDW). DPPH and reducing power activity model systems identified PMASA and PKRKL as the populations with the highest antioxidant activity. Finally, the data obtained represented valuable information for introducing elite populations with EO components favorable to pharmaceutical and industrial applications.     

Effects of Competition and Temporal Variation on the Evolutionary Potential of Two Native Bunchgrass Species

The capacity of restored plant populations to adapt to new environmental challenges depends on within‐population genetic variation. We examined how much genetic and environmentally based variation for fitness‐associated traits exists within populations of two native grasses commonly used for restoration in California. We were also interested in understanding how phenotypic expression of genetic variation for these traits varies with growth environment. Thirty maternal families of Elymus glaucus (Blue wild rye) and Nassella pulchra (Purple needlegrass) were sampled from both coastal and interior populations and reciprocally transplanted into three replicated common gardens with and without interspecific competition at each site. Reproductive output of families differed both among years and with competition treatments. Phenotypic expression of genetic variation in culm production differed among populations and was very low when families were grown with interspecific competition. Without interspecific competition, the degree of genetic determination peaked in year two in both species (8.4 and 15.1% in E. glaucus and N. pulchra, respectively). Significant genetic differences in reproduction and phenotypic plasticity were found among N. pulchra subpopulations sampled less than 3 km apart, further highlighting the importance of thoroughly sampling available genetic variation in populations used for restoration. The variable and generally low expression of genetic variation indicates that rates of adaptation in restored populations of these native grasses may vary temporally and may be especially slow within competitive environments.     

SNPs across time and space: population genomic signatures of founder events and epizootics in the House Finch (Haemorhous mexicanus)

Identifying genomic signatures of natural selection can be challenging against a background of demographic changes such as bottlenecks and population expansions. Here, we disentangle the effects of demography from selection in the House Finch (Haemorhous mexicanus) using samples collected before and after a pathogen‐induced selection event. Using ddRADseq, we genotyped over 18,000 SNPs across the genome in native pre‐epizootic western US birds, introduced birds from Hawaii and the eastern United States, post‐epizootic eastern birds, and western birds sampled across a similar time span. We found 14% and 7% reductions in nucleotide diversity, respectively, in Hawaiian and pre‐epizootic eastern birds relative to pre‐epizootic western birds, as well as elevated levels of linkage disequilibrium and other signatures of founder events. Despite finding numerous significant frequency shifts (outlier loci) between pre‐epizootic native and introduced populations, we found no signal of reduced genetic diversity, elevated linkage disequilibrium, or outlier loci as a result of the epizootic. Simulations demonstrate that the proportion of outliers associated with founder events could be explained by genetic drift. This rare view of genetic evolution across time in an invasive species provides direct evidence that demographic shifts like founder events have genetic consequences more widespread across the genome than natural selection.