Genome scans for divergent selection in natural populations of the widespread hardwood species Eucalyptus grandis (Myrtaceae) using microsatellites

Testing for fitness epistasis in a transplant experiment identifies a candidate adaptive locus in Timema stick insects

Identifying the genetic basis of adaptation is a central goal of evolutionary biology. However, identifying genes and mutations affecting fitness remains challenging because a large number of traits and variants can influence fitness. Selected phenotypes can also be difficult to know a priori, complicating top-down genetic approaches for trait mapping that involve crosses or genome-wide association studies. In such cases, experimental genetic approaches, where one maps fitness directly and attempts to infer the traits involved afterwards, can be valuable. Here, we re-analyse data from a transplant experiment involving Timema stick insects, where five physically clustered single-nucleotide polymorphisms associated with cryptic body coloration were shown to interact to affect survival. Our analysis covers a larger genomic region than past work and revealed a locus previously not identified as associated with survival. This locus resides near a gene, Punch (Pu), involved in pteridine pigments production, implying that it could be associated with an unmeasured coloration trait. However, by combining previous and newly obtained phenotypic data, we show that this trait is not eye or body coloration. We discuss the implications of our results for the discovery of traits, genes and mutations associated with fitness in other systems, as well as for supergene evolution. This article is part of the theme issue 'Genetic basis of adaptation and speciation: from loci to causative mutations'.

Climate Adaptation, Drought Susceptibility, and Genomic-Informed Predictions of Future Climate Refugia for the Australian Forest Tree Eucalyptus globulus

Understanding the capacity of forest tree species to adapt to climate change is of increasing importance for managing forest genetic resources. Through a genomics approach, we modelled spatial variation in climate adaptation within the Australian temperate forest tree Eucalyptus globulus, identified putative climate drivers of this genomic variation, and predicted locations of future climate refugia and populations at-risk of future maladaptation. Using 812,158 SNPs across 130 individuals from 30 populations (i.e., localities) spanning the species’ natural range, a gradientForest algorithm found 1177 SNPs associated with locality variation in home-site climate (climate-SNPs), putatively linking them to climate adaptation. Very few climate-SNPs were associated with population-level variation in drought susceptibility, signalling the multi-faceted nature and complexity of climate adaptation. Redundancy analysis (RDA) showed 24% of the climate-SNP variation could be explained by annual precipitation, isothermality, and maximum temperature of the warmest month. Spatial predictions of the RDA climate vectors associated with climate-SNPs allowed mapping of genomically informed climate selective surfaces across the species’ range under contemporary and projected future climates. These surfaces suggest over 50% of the current distribution of E. globulus will be outside the modelled adaptive range by 2070 and at risk of climate maladaptation. Such surfaces present a new integrated approach for natural resource managers to capture adaptive genetic variation and plan translocations in the face of climate change.

The characteristics of the complete chloroplast genome for Eucalyptus robusta (Myrtaceae)

Eucalyptus robusta is a great tree of Myrtaceae. It is mainly distributed in Southeast provinces of China and Southeast Asian countries. There is no study on the genome of E. robusta far. Here we report the complete chloroplast genome of E. robusta, assembled from whole-genome high-throughput sequencing data, as a resource for future studies on the taxonomy and evolution of E. robusta. The complete chloroplast genome of E. robusta is 160,201 bp in length with a typical quadripartite structure, consisting of a large single-copy region (LSC, 88,905 bp), a single-copy region (SSC, 18,498 bp) and a pair of inverted repeats (IRs, 26,399 bp). It was predicted to contain a total of 128 genes, with an overall GC content of 36.86%. Phylogenetic analysis placed E. robusta closest to Eucalyptus globulus.

Population Genetic Diversity and Structure of an Endangered Salicaceae Species in Northeast China: Chosenia arbutifolia (Pall.) A. Skv

Chosenia arbutifolia (Pall.) A. Skv. is a unique and endangered species belonging to the Salicaceae family. It has great potential for ornamental and industrial use. However, human interference has led to a decrease in and fragmentation of its natural populations in the past two decades. To effectively evaluate, utilize, and conserve available resources, the genetic diversity and population structure of C. arbutifolia were analyzed in this study. A total of 142 individuals from ten provenances were sampled and sequenced. Moderate diversity was detected among these, with a mean expected heterozygosity and Shannon’s Wiener index of 0.3505 and 0.5258, respectively. The inbreeding coefficient was negative, indicating a significant excess of heterozygotes. The fixation index varied from 0.0068 to 0.3063, showing a varied genetic differentiation between populations. Analysis of molecular variance demonstrated that differentiation accounted for 82.23% of the total variation among individuals, while the remaining 17.77% variation was between populations. Furthermore, the results of population structure analysis indicated that the 142 individuals originated from three primitive groups. To provide genetic information and help design conservation and management strategies, landscape genomics analysis was performed by investigating loci associated with environmental variables. Eighteen SNP markers were associated with altitude and annual average temperature, of which five were ascribed with specific functions. In conclusion, the current study furthers the understanding of C. arbutifolia genetic architecture and provides insights for germplasm protection.

EST-SSR-based landscape genetics of Pseudotaxus chienii, a tertiary relict conifer endemic to China

Pseudotaxus chienii, belonging to the monotypic genus Pseudotaxus (Taxaceae), is a relict conifer endemic to China. Its populations are usually small and patchily distributed, having a low capacity of natural regeneration. To gain a clearer understanding of how landscape variables affect the local adaptation of P. chienii, we applied EST-SSR markers in conjunction with landscape genetics methods: (a) to examine the population genetic pattern and spatial genetic structure; (b) to perform genome scan and selection scan to identify outlier loci and the associated landscape variables; and (c) to model the ecological niche under climate change. As a result, P. chienii was found to have a moderate level of genetic variation and a high level of genetic differentiation. Its populations displayed a significant positive relationship between the genetic and geographical distance (i.e., "isolation by distance" pattern) and a strong fine-scale spatial genetic structure within 2 km. A putatively adaptive locus EMS6 (functionally annotated to cellulose synthase A catalytic subunit 7) was identified, which was found significantly associated with soil Cu, K, and Pb content and the combined effects of temperature and precipitation. Moreover, P. chienii was predicted to experience significant range contractions in future climate change scenarios. Our results highlight the potential of specific soil metal content and climate variables as the driving force of adaptive genetic differentiation in P. chienii. The data would also be useful to develop a conservation action plan for P. chienii.

Seascape Genomics of the Sugar Kelp Saccharina latissima along the North Eastern Atlantic Latitudinal Gradient

Temperature is one of the most important range-limiting factors for many seaweeds. Driven by the recent climatic changes, rapid northward shifts of species' distribution ranges can potentially modify the phylogeographic signature of Last Glacial Maximum. We explored this question in detail in the cold-tolerant kelp species Saccharina latissima, using microsatellites and double digest restriction site-associated DNA sequencing ( ddRAD-seq) derived single nucleotide polymorphisms (SNPs) to analyze the genetic diversity and structure in 11 sites spanning the entire European Atlantic latitudinal range of this species. In addition, we checked for statistical correlation between genetic marker allele frequencies and three environmental proxies (sea surface temperature, salinity, and water turbidity). Our findings revealed that genetic diversity was significantly higher for the northernmost locality (Spitsbergen) compared to the southern ones (Northern Iberia), which we discuss in light of the current state of knowledge on phylogeography of S. latissima and the potential influence of the recent climatic changes on the population structure of this species. Seven SNPs and 12 microsatellite alleles were found to be significantly associated with at least one of the three environmental variables. We speculate on the putative adaptive functions of the genes associated with the outlier markers and the importance of these markers for successful conservation and aquaculture strategies for S. latissima in this age of rapid global change.

Genomic evidence of introgression and adaptation in a model subtropical tree species, Eucalyptus grandis

The genetic consequences of adaptation to changing environments can be deciphered using population genomics, which may help predict species' responses to global climate change. Towards this, we used genome-wide SNP marker analysis to determine population structure and patterns of genetic differentiation in terms of neutral and adaptive genetic variation in the natural range of Eucalyptus grandis, a widely cultivated subtropical and temperate species, serving as genomic reference for the genus. We analysed introgression patterns at subchromosomal resolution using a modified ancestry mapping approach and identified provenances with extensive interspecific introgression in response to increased aridity. Furthermore, we describe potentially adaptive genetic variation as explained by environment-associated SNP markers, which also led to the discovery of what is likely a large structural variant. Finally, we show that genes linked to these markers are enriched for biotic and abiotic stress responses.

Modelling Current and Future Potential Habitats for Plantations of Eucalyptus grandis Hill ex Maiden and E. dunnii Maiden in Uruguay

Eucalyptus grandis and E. dunnii have high productive potential in the South of Brazil, Uruguay, and central Argentina. This is based on the similarity of the climate and soil of these areas, which form an eco-region called Campos. However, previous results show that these species have differences in their distribution caused by the prioritization of Uruguayan soils for forestry, explained by the particular conditions of each site. In this study, the site variables (climate, soil, and topography) that better explain the distribution of both species were identified, and prediction models of current and future distribution were adjusted for different climate change scenarios (years 2050 and 2070). The distribution of E. grandis was associated with soil parameters, whereas for E. dunnii a greater effect of the climatic variables was observed. The ensemble biomod2 model was the most precise with regard to predicting the habitat for both species with respect to the simple models evaluated. For E. dunnii, the average values of the AUC, Kappa, and TSS index were 0.98, 0.88, and 0.77, respectively. For E. grandis, their values were 0.97, 0.86, and 0.80, respectively. In the projections of climatic change, the distribution of E. grandis occurrence remains practically unchanged, even in the scenarios of temperature increase. However, current distribution of E. dunnii shows high susceptibility in a scenario of increased temperature, to the point that most of the area currently planted may be at risk. Our results might be useful to political government and foresters for decision making in terms of future planted areas.

Genotyping-by-Sequencing of Gossypium hirsutum Races and Cultivars Uncovers Novel Patterns of Genetic Relationships and Domestication Footprints

Determining the genetic rearrangement and domestication footprints in Gossypium hirsutum cultivars and primitive race genotypes are essential for effective gene conservation efforts and the development of advanced breeding molecular markers for marker-assisted breeding. In this study, 94 accessions representing the 7 primitive races of G hirsutum, along with 9 G hirsutum and 12 Gossypium barbadense cultivated accessions were evaluated. The genotyping-by-sequencing (GBS) approach was employed and 146 558 single nucleotide polymorphisms (SNP) were generated. Distinct SNP signatures were identified through the combination of selection scans and association analyses. Phylogenetic analyses were also conducted, and we concluded that the Latifolium, Richmondi, and Marie-Galante race accessions were more genetically related to the G hirsutum cultivars and tend to cluster together. Fifty-four outlier SNP loci were identified by selection-scan analysis, and 3 SNPs were located in genes related to the processes of plant responding to stress conditions and confirmed through further genome-wide signals of marker-phenotype association analysis, which indicate a clear selection signature for such trait. These results identified useful candidate gene locus for cotton breeding programs.

Thirteen years under arid conditions: exploring marker-trait associations in Eucalyptus cladocalyx for complex traits related to flowering, stem form and growth

We present an association analysis for seven key traits related to flowering, stem form and growth in Eucalyptus cladocalyx, a tree species suitable for low rainfall sites, using a long-term progeny trial with 49 open-pollinated maternal families in the southern Atacama Desert, Chile. The progeny trial was carried out in an arid environment with a mean annual rainfall of 152 mm. Simple sequence repeats (SSR) from a full consensus map of Eucalyptus were used for genotyping 245 individual trees. Twenty-three significant marker-trait associations were identified, explaining between 5.9 and 23.7% of the phenotypic variance. The marker EMBRA101 located on LG10 at 56.5 cM was concomitantly associated with diameter at breast height and tree height. Nine SSR were significantly associated with stem forking and stem straightness, explaining between 5.9 and 14.8% of the phenotypic variation. To our knowledge, this is the first study reporting a SSR-based association mapping analysis for stem form traits in Eucalyptus. These results provide novel and valuable information for understanding the genetic base of key traits in E. cladocalyx for breeding purposes under arid conditions.

The first identification of genomic loci in plants associated with resistance to galling insects: a case study in Eucalyptus L'Hér. (Myrtaceae)

Genomic loci related with resistance to gall-inducing insects have not been identified in any plants. Here, association mapping was used to identify molecular markers for resistance to the gall wasp Leptocybe invasa in two Eucalyptus species. A total of 86 simple sequence repeats (SSR) markers were screened out from 839 SSRs and used for association mapping in E. grandis. By applying the mixed linear model, seven markers were identified to be associated significantly (P ≤ 0.05) with the gall wasp resistance in E. grandis, including two validated with a correction of permutation test (P ≤ 0.008). The proportion of the variance in resistance explained by a significant marker ranged from 3.3% to 37.8%. Four out of the seven significant associations in E. grandis were verified and also validated (P ≤ 0.073 in a permutation test) in E. tereticornis, with the variation explained ranging from 24.3% to 48.5%. Favourable alleles with positive effect were also mined from the significant markers in both species. These results provide insight into the genetic control of gall wasp resistance in plants and have great potential for marker-assisted selection for resistance to L. invasa in the important tree genus Eucalyptus.

The first identification of genomic loci in plants associated with resistance to galling insects: a case study in Eucalyptus L'Hér. (Myrtaceae)

Genomic loci related with resistance to gall-inducing insects have not been identified in any plants. Here, association mapping was used to identify molecular markers for resistance to the gall wasp Leptocybe invasa in two Eucalyptus species. A total of 86 simple sequence repeats (SSR) markers were screened out from 839 SSRs and used for association mapping in E. grandis. By applying the mixed linear model, seven markers were identified to be associated significantly (P ≤ 0.05) with the gall wasp resistance in E. grandis, including two validated with a correction of permutation test (P ≤ 0.008). The proportion of the variance in resistance explained by a significant marker ranged from 3.3% to 37.8%. Four out of the seven significant associations in E. grandis were verified and also validated (P ≤ 0.073 in a permutation test) in E. tereticornis, with the variation explained ranging from 24.3% to 48.5%. Favourable alleles with positive effect were also mined from the significant markers in both species. These results provide insight into the genetic control of gall wasp resistance in plants and have great potential for marker-assisted selection for resistance to L. invasa in the important tree genus Eucalyptus.