Regional differences in the abiotic environment contribute to genomic divergence within a wild tomato species

Deciphering salt stress responses in Solanum pimpinellifolium through high-throughput phenotyping

Soil salinity is a major environmental stressor affecting agricultural productivity worldwide. Understanding plant responses to salt stress is crucial for developing resilient crop varieties. Wild relatives of cultivated crops, such as wild tomato, Solanum pimpinellifolium, can serve as a useful resource to further expand the resilience potential of the cultivated germplasm, S. lycopersicum. In this study, we employed high-throughput phenotyping in the greenhouse and field conditions to explore salt stress responses of a S. pimpinellifolium diversity panel. Our study revealed extensive phenotypic variations in response to salt stress, with traits such as transpiration rate, shoot mass, and ion accumulation showing significant correlations with plant performance. We found that while transpiration was a key determinant of plant performance in the greenhouse, shoot mass strongly correlated with yield under field conditions. Conversely, ion accumulation was the least influential factor under greenhouse conditions. Through a Genome Wide Association Study, we identified candidate genes not previously associated with salt stress, highlighting the power of high-throughput phenotyping in uncovering novel aspects of plant stress responses. This study contributes to our understanding of salt stress tolerance in S. pimpinellifolium and lays the groundwork for further investigations into the genetic basis of these traits, ultimately informing breeding efforts for salinity tolerance in tomato and other crops.

Landscape genomics reveals genetic signals of environmental adaptation of African wild eggplants

Crop wild relatives (CWR) provide a valuable resource for improving crops. They possess desirable traits that confer resilience to various environmental stresses. To fully utilize crop wild relatives in breeding and conservation programs, it is important to understand the genetic basis of their adaptation. Landscape genomics associates environments with genomic variation and allows for examining the genetic basis of adaptation. Our study examined the differences in allele frequency of 15,416 single nucleotide polymorphisms (SNPs) generated through genotyping by sequencing approach among 153 accessions of 15 wild eggplant relatives and two cultivated species from Africa, the principal hotspot of these wild relatives. We also explored the correlation between these variations and the bioclimatic and soil conditions at their collection sites, providing a comprehensive understanding of the genetic signals of environmental adaptation in African wild eggplant. Redundancy analysis (RDA) results showed that the environmental variation explained 6% while the geographical distances among the collection sites explained 15% of the genomic variation in the eggplant wild relative populations when controlling for population structure. Our findings indicate that even though environmental factors are not the main driver of selection in eggplant wild relatives, it is influential in shaping the genomic variation over time. The selected environmental variables and candidate SNPs effectively revealed grouping patterns according to the environmental characteristics of sampling sites. Using four genotype-environment association methods, we detected 396 candidate SNPs (2.5% of the initial SNPs) associated with eight environmental factors. Some of these SNPs signal genes involved in pathways that help adapt to environmental stresses such as drought, heat, cold, salinity, pests, and diseases. These candidate SNPs will be useful for marker-assisted improvement and characterizing the germplasm of this crop for developing climate-resilient eggplant varieties. The study provides a model for applying landscape genomics to other crops' wild relatives.

Genetic and behavioral differences between above and below ground Culex pipiens bioforms

Efficiency of mosquito-borne disease transmission is dependent upon both the preference and fidelity of mosquitoes as they seek the blood of vertebrate hosts. While mosquitoes select their blood hosts through multi-modal integration of sensory cues, host-seeking is primarily an odor-guided behavior. Differences in mosquito responses to hosts and their odors have been demonstrated to have a genetic component, but the underlying genomic architecture of these responses has yet to be fully resolved. Here, we provide the first characterization of the genomic architecture of host preference in the polymorphic mosquito species, Culex pipiens. The species exists as two morphologically identical bioforms, each with distinct avian and mammalian host preferences. Cx. pipiens females with empirically measured host responses were prepared into reduced representation DNA libraries and sequenced to identify genomic regions associated with host preference. Multiple genomic regions associated with host preference were identified on all 3 Culex chromosomes, and these genomic regions contained clusters of chemosensory genes, as expected based on work in Anopheles gambiae complex mosquitoes and in Aedes aegypti. One odorant receptor and one odorant binding protein gene showed one-to-one orthologous relationships to differentially expressed genes in A. gambiae complex members with divergent host preferences. Overall, our work identifies a distinct set of odorant receptors and odorant binding proteins that may enable Cx. pipiens females to distinguish between their vertebrate blood host species, and opens avenues for future functional studies that could measure the unique contributions of each gene to host preference phenotypes.

Genomic diversity and structure of a Neotropical microendemic fig tree

Genetic diversity is a key component of evolution, and unraveling factors that promote genetic differentiation in space and time is a central question in evolutionary biology. One of the most diverse and ecologically important tree genera in tropical forests worldwide is Ficus (Moraceae). It has been suggested that, given the great dispersal capacity of pollinating fig wasps (Chalcidoidea; Agaonidae), the spatial genetic structure, particularly in monoecious fig species, should be weak. However, no studies have addressed the factors that determine the genetic structure of Ficus species in regions of high geological, geographic, and climatic complexity, such as the Mexican Transition Zone. Using nuclear single nucleotide polymorphisms (5311 SNPs) derived from low-coverage whole genomes and 17 populations, we analyzed the population genomics of Ficus pringlei to characterize neutral and adaptive genetic variation and structure and its association with geographic barriers such as the Trans-Mexican Volcanic Belt, environmental heterogeneity, and wind connectivity. From genomic data of 71 individuals, high genetic diversity, and the identification of three genomic lineages were recorded (North, South, and Churumuco). The results suggest that genetic variation is primarily determined by climatic heterogeneity. Ficus pringlei populations from the north and south of the Trans-Mexican Volcanic Belt also exhibited minimal genetic differentiation (F ST = 0.021), indicating that this mountain range may not act as an insurmountable barrier to gene flow. Wind connectivity is also highlighted in structuring putative adaptive genetic variation, underscoring the intricate complexity of the various factors influencing genetic variation in the species. This study provides information on the possible mechanisms underlying the genetic variation of endemic species of the tropical dry forest of Western Mexico, such as F. pringlei.

Epigenetic architecture of Pseudotaxus chienii: Revealing the synergistic effects of climate and soil variables

Whether conifers can withstand environmental changes especially temperature fluctuations has been controversial. Epigenetic analysis may provide new perspectives for solving the issue. Pseudotaxus chienii is an endangered gymnosperm species endemic to China. In this study, we have examined the genetic and epigenetic variations in its natural populations aiming to disentangle the synergistic effects of climate and soil on its population (epi)genetic differentiation by using amplified fragment length polymorphism (AFLP) and methylation-sensitive AFLP (MSAP) techniques. We identified 23 AFLP and 26, 7, and 5 MSAP outliers in P. chienii. Twenty-one of the putative adaptive AFLP loci were found associated with climate and/or soil variables including precipitation, temperature, K, Fe, Zn, and Cu, whereas 21, 7, and 4 MSAP outliers were significantly related to precipitation of wettest month (Bio13), precipitation driest of month (Bio14), percent tree cover (PTC), and soil Fe, Mn, and Cu compositions. Total precipitation and precipitation in the driest seasons were the most influential factors for genetic and epigenetic variation, respectively. In addition, a high full-methylation level and a strong correlation between genetic and epigenetic variation were detected in P. chienii. Climate is found of greater importance than soil in shaping adaptive (epi)genetic differentiation, and the synergistic effects of climate and climate-soil variables were also observed. The identified climate and soil variables should be considered when applying ex situ conservation.

Climate has contributed to population diversification of Daphnia galeata across Eurasia

Climate is a fundamental abiotic factor that plays a key role in driving the evolution, distribution and population diversification of species. However, there have been few investigations of genomic signatures of adaptation to local climatic conditions in cladocerans. Here, we have provided the first high-quality chromosome-level genome assembly (~143 Mb, scaffold N50 12.6 Mb) of the waterflea, Daphnia galeata, and investigated genomic variation in 22 populations from Central Europe and Eastern China. Our ecological-niche models suggested that the historic distribution of D. galeata in Eurasia was significantly affected by Quaternary climate fluctuations. We detected pronounced genomic and morphometric divergences between European and Chinese D. galeata populations. Such divergences could be partly explained by genomic signatures of thermal adaptation to distinct climate regimes: a set of candidate single-nucleotide polymorphisms (SNPs) potentially associated with climate were detected. These SNPs were in genes significantly enriched in the Gene ontology terms "determination of adult lifespan" and "translation repressor activity", and especially, mthl5 and SOD1 involved in the IIS pathway, and EIF4EBP2 involved in the target of the rapamycin signalling pathway. Our study indicates that certain alleles might be associated with particular temperature regimes, playing a functional role in shaping the population structure of D. galeata at a large geographical scale. These results highlight the potential role of molecular variation in the response to climate variation, in the context of global climate change.

Population genomic analysis reveals geographic structure and climatic diversification for Macrophomina phaseolina isolated from soybean and dry bean across the United States, Puerto Rico, and Colombia

Macrophomina phaseolina causes charcoal rot, which can significantly reduce yield and seed quality of soybean and dry bean resulting from primarily environmental stressors. Although charcoal rot has been recognized as a warm climate-driven disease of increasing concern under global climate change, knowledge regarding population genetics and climatic variables contributing to the genetic diversity of M. phaseolina is limited. This study conducted genome sequencing for 95 M. phaseolina isolates from soybean and dry bean across the continental United States, Puerto Rico, and Colombia. Inference on the population structure using 76,981 single nucleotide polymorphisms (SNPs) revealed that the isolates exhibited a discrete genetic clustering at the continental level and a continuous genetic differentiation regionally. A majority of isolates from the United States (96%) grouped in a clade with a predominantly clonal genetic structure, while 88% of Puerto Rican and Colombian isolates from dry bean were assigned to a separate clade with higher genetic diversity. A redundancy analysis (RDA) was used to estimate the contributions of climate and spatial structure to genomic variation (11,421 unlinked SNPs). Climate significantly contributed to genomic variation at a continental level with temperature seasonality explaining the most variation while precipitation of warmest quarter explaining the most when spatial structure was accounted for. The loci significantly associated with multivariate climate were found closely to the genes related to fungal stress responses, including transmembrane transport, glycoside hydrolase activity and a heat-shock protein, which may mediate climatic adaptation for M. phaseolina. On the contrary, limited genome-wide differentiation among populations by hosts was observed. These findings highlight the importance of population genetics and identify candidate genes of M. phaseolina that can be used to elucidate the molecular mechanisms that underly climatic adaptation to the changing climate.

Population structure analysis to explore genetic diversity and geographical distribution characteristics of wild tea plant in Guizhou Plateau

BACKGROUND

Tea, the second largest consumer beverage in the world after water, is widely cultivated in tropical and subtropical areas. However, the effect of environmental factors on the distribution of wild tea plants is unclear.

RESULTS

A total of 159 wild tea plants were collected from different altitudes and geological types of the Guizhou Plateau. Using the genotyping-by-sequencing method, 98,241 high-quality single nucleotide polymorphisms were identified. Genetic diversity, population structure analysis, principal component analysis, phylogenetic analysis, and linkage disequilibrium were performed. The genetic diversity of the wild tea plant population from the Silicate Rock Classes of Camellia gymnogyna was higher than that from the Carbonate Rock Classes of Camellia tachangensis. In addition, the genetic diversity of wild tea plants from the second altitude gradient was significantly higher than that of wild tea plants from the third and first altitude gradients. Two inferred pure groups (GP01 and GP02) and one inferred admixture group (GP03) were identified by population structure analysis and were verified by principal component and phylogenetic analyses. The highest differentiation coefficients were determined for GP01 vs. GP02, while the lowest differentiation coefficients were determined for GP01 vs. GP03.

CONCLUSIONS

This study revealed the genetic diversity and geographical distribution characteristics of wild tea plants in the Guizhou Plateau. There are significant differences in genetic diversity and evolutionary direction between Camellia tachangensis with Carbonate Rock Classes at the first altitude gradient and Camellia gymnogyna with Silicate Rock Classes at the third altitude gradient. Geological environment, soil mineral element content, soil pH, and altitude markedly contributed to the genetic differentiation between Camellia tachangensis and Camellia gymnogyna.

Genotype-environment associations to reveal the molecular basis of environmental adaptation

A fundamental goal in plant biology is to identify and understand the variation underlying plants' adaptation to their environment. Climate change has given new urgency to this goal, as society aims to accelerate adaptation of ecologically important plant species, endangered plant species, and crops to hotter, less predictable climates. In the pre-genomic era, identifying adaptive alleles was painstaking work, leveraging genetics, molecular biology, physiology, and ecology. Now, the rise of genomics and new computational approaches may facilitate this research. Genotype-environment associations (GEAs) use statistical associations between allele frequency and environment of origin to test the hypothesis that allelic variation at a given gene is adapted to local environments. Researchers may scan the genome for GEAs to generate hypotheses on adaptive genetic variants (environmental genome-wide association studies). Despite the rapid adoption of these methods, many important questions remain about the interpretation of GEA findings, which arise from fundamental unanswered questions on the genetic architecture of adaptation and limitations inherent to association-based analyses. We outline strategies to ground GEAs in the underlying hypotheses of genetic architecture and better test GEA-generated hypotheses using genetics and ecophysiology. We provide recommendations for new users who seek to learn about the molecular basis of adaptation. When combined with a rigorous hypothesis testing framework, GEAs may facilitate our understanding of the molecular basis of climate adaptation for plant improvement.

Using landscape genomics to assess local adaptation and genomic vulnerability of a perennial herb Tetrastigma hemsleyanum (Vitaceae) in subtropical China

Understanding adaptive genetic variation of plant populations and their vulnerabilities to climate change are critical to preserve biodiversity and subsequent management interventions. To this end, landscape genomics may represent a cost-efficient approach for investigating molecular signatures underlying local adaptation. Tetrastigma hemsleyanum is, in its native habitat, a widespread perennial herb of warm-temperate evergreen forest in subtropical China. Its ecological and medicinal values constitute a significant revenue for local human populations and ecosystem. Using 30,252 single nucleotide polymorphisms (SNPs) derived from reduced-representation genome sequencing in 156 samples from 24 sites, we conducted a landscape genomics study of the T. hemsleyanum to elucidate its genomic variation across multiple climate gradients and genomic vulnerability to future climate change. Multivariate methods identified that climatic variation explained more genomic variation than that of geographical distance, which implied that local adaptation to heterogeneous environment might represent an important source of genomic variation. Among these climate variables, winter precipitation was the strongest predictor of the contemporary genetic structure. F _ST outlier tests and environment association analysis totally identified 275 candidate adaptive SNPs along the genetic and environmental gradients. SNP annotations of these putatively adaptive loci uncovered gene functions associated with modulating flowering time and regulating plant response to abiotic stresses, which have implications for breeding and other special agricultural aims on the basis of these selection signatures. Critically, modelling revealed that the high genomic vulnerability of our focal species via a mismatch between current and future genotype-environment relationships located in central-northern region of the T. hemsleyanum's range, where populations require proactive management efforts such as assistant adaptation to cope with ongoing climate change. Taken together, our results provide robust evidence of local climate adaption for T. hemsleyanum and further deepen our understanding of adaptation basis of herbs in subtropical China.

Effects of the neo-X chromosome on genomic signatures of hybridization in Rumex hastatulus

Natural hybrid zones provide opportunities for studies of the evolution of reproductive isolation in wild populations. Although recent investigations have found that the formation of neo-sex chromosomes is associated with reproductive isolation, the mechanisms remain unclear in most cases. Here, we assess the contemporary structure of gene flow in the contact zone between largely allopatric cytotypes of the dioecious plant Rumex hastatulus, a species with evidence of sex chromosome turn-over. Males to the west of the Mississippi river, USA, have an X and a single Y chromosome, whereas populations to the east of the river have undergone a chromosomal rearrangement giving rise to a larger X and two Y chromosomes. Using reduced-representation sequencing, we provide evidence that hybrids form readily and survive multiple backcross generations in the field, demonstrating the potential for ongoing gene flow between the cytotypes. Cline analysis of each chromosome separately captured no signals of difference in cline shape between chromosomes. However, principal component regression revealed a significant increase in the contribution of individual SNPs to inter-cytotype differentiation on the neo-X chromosome, but no correlation with recombination rate. Cline analysis revealed that the only SNPs with significantly steeper clines than the genome average were located on the neo-X. Our data are consistent with a role for neo-sex chromosomes in reproductive isolation between R. hastatulus cytotypes. Our investigation highlights the importance of studying plant hybrid zones for understanding the evolution of sex chromosomes.

Landscape Genomics Provides Evidence of Ecotypic Adaptation and a Barrier to Gene Flow at Treeline for the Arctic Foundation Species Eriophorum vaginatum

Global climate change has resulted in geographic range shifts of flora and fauna at a global scale. Extreme environments, like the Arctic, are seeing some of the most pronounced changes. This region covers 14% of the Earth's land area, and while many arctic species are widespread, understanding ecotypic variation at the genomic level will be important for elucidating how range shifts will affect ecological processes. Tussock cottongrass (Eriophorum vaginatum L.) is a foundation species of the moist acidic tundra, whose potential decline due to competition from shrubs may affect ecosystem stability in the Arctic. We used double-digest Restriction Site-Associated DNA sequencing to identify genomic variation in 273 individuals of E. vaginatum from 17 sites along a latitudinal gradient in north central Alaska. These sites have been part of 30 + years of ecological research and are inclusive of a region that was part of the Beringian refugium. The data analyses included genomic population structure, demographic models, and genotype by environment association. Genome-wide SNP investigation revealed environmentally associated variation and population structure across the sampled range of E. vaginatum, including a genetic break between populations north and south of treeline. This structure is likely the result of subrefugial isolation, contemporary isolation by resistance, and adaptation. Forty-five candidate loci were identified with genotype-environment association (GEA) analyses, with most identified genes related to abiotic stress. Our results support a hypothesis of limited gene flow based on spatial and environmental factors for E. vaginatum, which in combination with life history traits could limit range expansion of southern ecotypes northward as the tundra warms. This has implications for lower competitive attributes of northern plants of this foundation species likely resulting in changes in ecosystem productivity.

Population genetic differentiation and genomic signatures of adaptation to climate in an abundant lizard

Species distributed across climatic gradients will typically experience spatial variation in selection, but gene flow can prevent such selection from causing population genetic differentiation and local adaptation. Here, we studied genomic variation of 415 individuals across 34 populations of the common wall lizard (Podarcis muralis) in central Italy. This species is highly abundant throughout this region and populations belong to a single genetic lineage, yet there is extensive phenotypic variation across climatic regimes. We used redundancy analysis to, first, quantify the effect of climate and geography on population genomic variation in this region and, second, to test if climate consistently sorts specific alleles across the landscape. Climate explained 5% of the population genomic variation across the landscape, about half of which was collinear with geography. Linear models and redundancy analyses identified loci that were significantly differentiated across climatic regimes. These loci were distributed across the genome and physically associated with genes putatively involved in thermal tolerance, regulation of temperature-dependent metabolism and reproductive activity, and body colouration. Together, these findings suggest that climate can exercise sufficient selection in lizards to promote genetic differentiation across the landscape in spite of high gene flow.

Weaving place‐based knowledge for culturally significant species in the age of genomics: Looking to the past to navigate the future

Relationships with place provide critical context for characterizing biocultural diversity. Yet, genetic and genomic studies are rarely informed by Indigenous or local knowledge, processes, and practices, including the movement of culturally significant species. Here, we show how place‐based knowledge can better reveal the biocultural complexities of genetic or genomic data derived from culturally significant species. As a case study, we focus on culturally significant southern freshwater kōura (crayfish) in Aotearoa me Te Waipounamu (New Zealand, herein Aotearoa NZ). Our results, based on genotyping‐by‐sequencing markers, reveal strong population genetic structure along with signatures of population admixture in 19 genetically depauperate populations across the east coast of Te Waipounamu. Environment association and differentiation analyses for local adaptation also indicate a role for hydroclimatic variables—including temperature, precipitation, and water flow regimes—in shaping local adaptation in kōura. Through trusted partnerships between community and researchers, weaving genomic markers with place‐based knowledge has both provided invaluable context for the interpretation of data and created opportunities to reconnect people and place. We envisage such trusted partnerships guiding future genomic research for culturally significant species in Aotearoa NZ and beyond.

Haplotype analyses reveal novel insights into tomato history and domestication driven by long-distance migrations and latitudinal adaptations

A novel haplotype-based approach that uses Procrustes analysis and automatic classification was used to provide further insights into tomato history and domestication. Agrarian societies domesticated species of interest by introducing complex genetic modifications. For tomatoes, two species, one of which had two botanical varieties, are thought to be involved in its domestication: the fully wild Solanum pimpinellifolium (SP), the wild and semi-domesticated Solanum lycopersicum var. cerasiforme (SLC) and the cultivated S. l. var. lycopersicum (SLL). The Procrustes approach showed that SP evolved into SLC during a gradual migration from the Peruvian deserts to the Mexican rainforests and that Peruvian and Ecuadorian SLC populations were the result of more recent hybridizations. Our model was supported by independent evidence, including ecological data from the accession collection site and morphological data. Furthermore, we showed that photosynthesis-, and flowering time-related genes were selected during the latitudinal migrations.

Genomic population structure and local adaptation of the wild strawberry Fragaria nilgerrensis

The crop wild relative, Fragaria nilgerrensis, is adapted to a variety of diverse habitats across its native range in China. Thus, discoveries made in this species could serve useful to guide the development of new superior strawberry cultivars that are resilient to new or variable environments. However, the genetic diversity and genetic architecture of traits in this species underlying important adaptive traits remain poorly understood. Here, we used whole-genome resequencing data from 193 F. nilgerrensis individuals spanning the distribution range in China to investigate the genetic diversity, population structure and the genomic basis of local adaptation. We identified four genetic groups, with the western group located in Hengduan Mountains exhibited the highest genetic diversity. Redundancy analysis suggests that both environment and geographic variables shaped a significant proportion of genomic variation. Our analyses revealed that the environmental difference explains more of the observed genetic variation than geographic distance. This suggests that adaptation to distinct habitats, unique combination of abiotic factors, likely drove genetic differentiation. Lastly, by implementing selective sweeps scans and genome-environment association analysis throughout the genome, we identified the genetic variation associated with local adaptation and investigated the functions of putative candidate genes in F. nilgerrensis.

Physical geography, isolation by distance and environmental variables shape genomic variation of wild barley (Hordeum vulgare L. ssp. spontaneum) in the Southern Levant

Determining the extent of genetic variation that reflects local adaptation in crop-wild relatives is of interest for the purpose of identifying useful genetic diversity for plant breeding. We investigated the association of genomic variation with geographical and environmental factors in wild barley (Hordeum vulgare L. ssp. spontaneum) populations of the Southern Levant using genotyping by sequencing (GBS) of 244 accessions in the Barley 1K+ collection. The inference of population structure resulted in four genetic clusters that corresponded to eco-geographical habitats and a significant association between lower gene flow rates and geographical barriers, e.g. the Judaean Mountains and the Sea of Galilee. Redundancy analysis (RDA) revealed that spatial autocorrelation explained 45% and environmental variables explained 15% of total genomic variation. Only 4.5% of genomic variation was solely attributed to environmental variation if the component confounded with spatial autocorrelation was excluded. A synthetic environmental variable combining latitude, solar radiation, and accumulated precipitation explained the highest proportion of genomic variation (3.9%). When conditioned on population structure, soil water capacity was the most important environmental variable explaining 1.18% of genomic variation. Genome scans with outlier analysis and genome-environment association studies were conducted to identify adaptation signatures. RDA and outlier methods jointly detected selection signatures in the pericentromeric regions, which have reduced recombination, of the chromosomes 3H, 4H, and 5H. However, selection signatures mostly disappeared after correction for population structure. In conclusion, adaptation to the highly diverse environments of the Southern Levant over short geographical ranges had a limited effect on the genomic diversity of wild barley. This highlighted the importance of nonselective forces in genetic differentiation.

Limited dispersal and local adaptation promote allopatric speciation in a biodiversity hotspot

Recently diverged or diverging populations can offer unobstructed insights into early barriers to gene flow during the initial stages of speciation. The current study utilised a novel insect system (order Mantophasmatodea) to shed light on the early drivers of speciation. The members of this group have limited dispersal abilities, small allopatric distributions and strong habitat associations in the Cape Floristic Region biodiversity hotspot in South Africa. Sister taxa from the diverse family Austrophasmatidae were chosen as focal species (Karoophasma biedouwense, K. botterkloofense). Population genetics and Generalized Dissimilarity Modelling (GDM) were used to characterise spatial patterns of genetic variation and evaluate the contribution of environmental factors to population divergence and speciation. Extensive sampling confirmed the suspected allopatry of these taxa. However, hybrids were identified in a narrow region occurring between the species' distributions. Strong population structure was found over short geographic distances; particularly in K. biedouwense in which geographic distance accounted for 32% of genetic variation over a scale of 50 km (r = .56, p < .001). GDM explained 42%-78% of the deviance in observed genetic dissimilarities. Geographic distance was consistently indicated to be important for between species and within population differentiation, suggesting that limited dispersal ability may be an important neutral driver of divergence. Temperature, altitude, precipitation and vegetation were also indicated as important factors, suggesting the possible role of adaptation to local environmental conditions for species divergence. The discovery of the hybrid-zone, and the multiple allopatric species pairs in Austrophasmatidae support the idea that this could be a promising group to further our understanding of speciation modes.

Genome-Wide Identification and Evolutionary Analysis of the SRO Gene Family in Tomato

SRO (SIMILAR TO RCD ONE) is a family of plant-specific small molecule proteins that play an important role in plant growth and development and environmental responses. However, SROs still lack systematic characterization in tomato. Based on bioinformatics methods, SRO family genes were identified and characterized from cultivated tomatoes and several wild tomatoes. qRT-PCR was used to study the expression of SRO gene in cultivated tomatoes. Phylogenetic and evolutionary analyses showed that SRO genes in angiosperms share a common ancestor and that the number of SRO family members changed as plants diverged and evolved. Cultivated tomato had six SRO members, five of which still shared some degree of identity with the ancestral SRO genes. Genetic structure and physicochemical properties showed that tomato SRO genes were highly conserved with chromosomal distribution. They could be divided into three groups based on exon-intron structure, and cultivated tomato contained only two of these subclades. A number of hormonal, light and abiotic stress-responsive cis-regulatory elements were identified from the promoter of the tomato SRO gene, and they also interacted with a variety of stress-responsive proteins and microRNAs. RNA-seq analysis showed that SRO genes were widely expressed in different tissues and developmental stages of tomato, with significant tissue-specific features. Expression analysis also showed that SRO genes respond significantly to high temperature and salt stress and mediate the tomato hormone regulatory network. These results provide a theoretical basis for further investigation of the functional expression of tomato SRO genes and provide potential genetic resources for tomato resistance breeding.

Introgression shapes fruit color convergence in invasive Galápagos tomato

Invasive species represent one of the foremost risks to global biodiversity. Here, we use population genomics to evaluate the history and consequences of an invasion of wild tomato-Solanum pimpinellifolium-onto the Galápagos Islands from continental South America. Using >300 archipelago and mainland collections, we infer this invasion was recent and largely the result of a single event from central Ecuador. Patterns of ancestry within the genomes of invasive plants also reveal post-colonization hybridization and introgression between S. pimpinellifolium and the closely related Galápagos endemic Solanum cheesmaniae. Of admixed invasive individuals, those that carry endemic alleles at one of two different carotenoid biosynthesis loci also have orange fruits-characteristic of the endemic species-instead of typical red S. pimpinellifolium fruits. We infer that introgression of two independent fruit color loci explains this observed trait convergence, suggesting that selection has favored repeated transitions of red to orange fruits on the Galápagos.

Population transcriptomic sequencing reveals allopatric divergence and local adaptation in Pseudotaxus chienii (Taxaceae)

BACKGROUND

Elucidating the effects of geography and selection on genetic variation is critical for understanding the relative importance of adaptation in driving differentiation and identifying the environmental factors underlying its occurrence. Adaptive genetic variation is common in tree species, especially widely distributed long-lived species. Pseudotaxus chienii can occupy diverse habitats with environmental heterogeneity and thus provides an ideal material for investigating the process of population adaptive evolution. Here, we characterize genetic and expression variation patterns and investigate adaptive genetic variation in P. chienii populations.

RESULTS

We generated population transcriptome data and identified 13,545 single nucleotide polymorphisms (SNPs) in 5037 unigenes across 108 individuals from 10 populations. We observed lower nucleotide diversity (π = 0.000701) among the 10 populations than observed in other gymnosperms. Significant negative correlations between expression diversity and nucleotide diversity in eight populations suggest that when the species adapts to the surrounding environment, gene expression and nucleotide diversity have a reciprocal relationship. Genetic structure analyses indicated that each distribution region contains a distinct genetic group, with high genetic differentiation among them due to geographical isolation and local adaptation. We used FST outlier, redundancy analysis, and latent factor mixed model methods to detect molecular signatures of local adaptation. We identified 244 associations between 164 outlier SNPs and 17 environmental variables. The mean temperature of the coldest quarter, soil Fe and Cu contents, precipitation of the driest month, and altitude were identified as the most important determinants of adaptive genetic variation. Most candidate unigenes with outlier signatures were related to abiotic and biotic stress responses, and the monoterpenoid biosynthesis and ubiquitin-mediated proteolysis KEGG pathways were significantly enriched in certain populations and deserve further attention in other long-lived trees.

CONCLUSIONS

Despite the strong population structure in P. chienii, genomic data revealed signatures of divergent selection associated with environmental variables. Our research provides SNPs, candidate unigenes, and biological pathways related to environmental variables to facilitate elucidation of the genetic variation in P. chienii in relation to environmental adaptation. Our study provides a promising tool for population genomic analyses and insights into the molecular basis of local adaptation.

Allelic Diversity at Abiotic Stress Responsive Genes in Relationship to Ecological Drought Indices for Cultivated Tepary Bean, Phaseolus acutifolius A. Gray, and Its Wild Relatives

Some of the major impacts of climate change are expected in regions where drought stress is already an issue. Grain legumes are generally drought susceptible. However, tepary bean and its wild relatives within Phaseolus acutifolius or P. parvifolius are from arid areas between Mexico and the United States. Therefore, we hypothesize that these bean accessions have diversity signals indicative of adaptation to drought at key candidate genes such as: Asr2, Dreb2B, and ERECTA. By sequencing alleles of these genes and comparing to estimates of drought tolerance indices from climate data for the collection site of geo-referenced, tepary bean accessions, we determined the genotype x environmental association (GEA) of each gene. Diversity analysis found that cultivated and wild P. acutifolius were intermingled with var. tenuifolius and P. parvifolius, signifying that allele diversity was ample in the wild and cultivated clade over a broad sense (sensu lato) evaluation. Genes Dreb2B and ERECTA harbored signatures of directional selection, represented by six SNPs correlated with the environmental drought indices. This suggests that wild tepary bean is a reservoir of novel alleles at genes for drought tolerance, as expected for a species that originated in arid environments. Our study corroborated that candidate gene approach was effective for marker validation across a broad genetic base of wild tepary accessions.

Tapping Diversity From the Wild: From Sampling to Implementation

The diversity observed among crop wild relatives (CWRs) and their ability to flourish in unfavorable and harsh environments have drawn the attention of plant scientists and breeders for many decades. However, it is also recognized that the benefit gained from using CWRs in breeding is a potential rose between thorns of detrimental genetic variation that is linked to the trait of interest. Despite the increased interest in CWRs, little attention was given so far to the statistical, analytical, and technical considerations that should guide the sampling design, the germplasm characterization, and later its implementation in breeding. Here, we review the entire process of sampling and identifying beneficial genetic variation in CWRs and the challenge of using it in breeding. The ability to detect beneficial genetic variation in CWRs is strongly affected by the sampling design which should be adjusted to the spatial and temporal variation of the target species, the trait of interest, and the analytical approach used. Moreover, linkage disequilibrium is a key factor that constrains the resolution of searching for beneficial alleles along the genome, and later, the ability to deplete linked deleterious genetic variation as a consequence of genetic drag. We also discuss how technological advances in genomics, phenomics, biotechnology, and data science can improve the ability to identify beneficial genetic variation in CWRs and to exploit it in strive for higher-yielding and sustainable crops.

Predicting Thermal Adaptation by Looking Into Populations' Genomic Past

Molecular evolution offers an insightful theory to interpret the genomic consequences of thermal adaptation to previous events of climate change beyond range shifts. However, disentangling often mixed footprints of selective and demographic processes from those due to lineage sorting, recombination rate variation, and genomic constrains is not trivial. Therefore, here we condense current and historical population genomic tools to study thermal adaptation and outline key developments (genomic prediction, machine learning) that might assist their utilization for improving forecasts of populations' responses to thermal variation. We start by summarizing how recent thermal-driven selective and demographic responses can be inferred by coalescent methods and in turn how quantitative genetic theory offers suitable multi-trait predictions over a few generations via the breeder's equation. We later assume that enough generations have passed as to display genomic signatures of divergent selection to thermal variation and describe how these footprints can be reconstructed using genome-wide association and selection scans or, alternatively, may be used for forward prediction over multiple generations under an infinitesimal genomic prediction model. Finally, we move deeper in time to comprehend the genomic consequences of thermal shifts at an evolutionary time scale by relying on phylogeographic approaches that allow for reticulate evolution and ecological parapatric speciation, and end by envisioning the potential of modern machine learning techniques to better inform long-term predictions. We conclude that foreseeing future thermal adaptive responses requires bridging the multiple spatial scales of historical and predictive environmental change research under modern cohesive approaches such as genomic prediction and machine learning frameworks.

Comparative Analyses of Phytochemical Variation Within and Between Congeneric Species of Willow Herb, Epilobium hirsutum and E. parviflorum: Contribution of Environmental Factors

The plants in the Epilobium genus are considered to have several important medicinal properties due to their unique chemical composition. Although metabolic profiles of medicinal plants are mainly controlled by genetic factors, their production is also to some degree influenced by environmental factors, thus, variations in the levels of phytochemicals may represent long-term ecological and evolutionary interactions. In order to depict the magnitude of natural variation in level of chemical compounds among conspecific populations of Epilobium hirsutum (n = 31) and E. parviflorum (n = 16), metabolite profiling of aerial parts of plants was performed with gas chromatography/mass spectrometry analysis. Putative identification and structure annotation revealed the presence of 74 compounds including 46 compounds considered secondary metabolites categorized into flavonoids (n = 8), phenolic acids (n = 26), steroids (n = 3), and terpenes (n = 5) across all populations. Although there was a considerable natural variation among conspecific populations, principal component analysis revealed a clear separation of populations of each species based on the second main principal component which was highly correlated with eight secondary metabolites. The level of secondary metabolites was significantly correlated between species (r = 0.91), suggesting shared metabolic pathways underlying the production of chemical compounds. In addition, redundancy and variance partitioning analyses by including bioclimatic variables and altitude revealed a significant contribution of elevation in explaining the total variation of secondary metabolites in E. hirsutum. Two-thirds of all secondary metabolites were significantly correlated with altitude in E. hirsutum. The large-scale geographic analyses of populations revealed additionally detected flavonoids and terpenes (E. hirsutum and E. parviflorum) and steroids (E. hirsutum) for the first time. This study provides significant information on additional chemical compounds found across the distribution range of the two ecologically important species of willow herb and emphasizes the importance of geographic-wide sampling as a valuable strategy to depict intraspecific and interspecific variability in chemical traits.