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

Phylogeny and evolutionary dynamics of the Rubia genus based on the chloroplast genome of Rubia tibetica

Rubia tibetica, a well-known medicinal plant, holds significant medicinal value. In this study, we sequenced and assembled the complete chloroplast genome of R. tibetica and conducted comparative analyses with four other Rubia species to provide insights into genome structure, selective pressure, and phylogenetic relationships. The chloroplast genome of R. tibetica was 154,901 bp in length and exhibited a typical quadripartite circular structure. Comparative analysis revealed a highly conserved genome structure among Rubia species, with genome sizes ranging from 153,555 bp to 155,108 bp and a GC content of approximately 37%. Codon usage analysis indicated a preference for A/U-ending codons. Selective pressure analysis identified relaxed selection in genes such as rps15 and petA, suggesting species-specific adaptive strategies. Two hypervariable regions, rps16-trnQ and psaJ-rps18, were experimentally validated as potential molecular markers through PCR amplification and sequencing. However, sequencing results of the rps16-trnQ region differed from the reference chloroplast genome, requiring further investigation. Phylogenetic analysis was performed using maximum likelihood (ML) based on the complete chloroplast genome sequences of five Rubia species and confirmed the monophyly of the Rubia genus, with R. tibetica occupying a basal position. This suggesting the retention of ancestral traits and adaptation to the unique environment of the Qinghai-Tibet Plateau. Molecular dating analysis further revealed divergence times consistent with geological and climatic events in the region. These findings provide critical insights into the evolutionary dynamics of Rubia species and offer valuable resources for their classification, conservation, and medicinal application.

Transcriptomic and Physiological Analyses for the Role of Hormones and Sugar in Axillary Bud Development of Wild Strawberry Stolon

Strawberries are mainly propagated by stolons, which can be divided into monopodial and sympodial types. Monopodial stolons consistently produce ramets at each node following the initial single dormant bud, whereas sympodial stolons develop a dormant bud before each ramet. Sympodial stolon encompasses both dormant buds and ramet buds, making it suitable for studying the formation mechanism of different stolon types. In this study, we utilized sympodial stolons from Fragaria nilgerrensis as materials and explored the mechanisms underlying sympodial stolon development through transcriptomic and phytohormonal analyses. The transcriptome results unveiled that auxin, cytokinin, and sugars likely act as main regulators. Endogenous hormone analysis revealed that the inactivation of auxin could influence bud dormancy. Exogenous cytokinin application primarily induced dormant buds to develop into secondary stolons, with the proportion of ramet formation being very low, less than 10%. Furthermore, weighted gene co-expression network analysis identified key genes involved in ramet formation, including auxin transport and response genes, the cytokinin activation gene LOG1, and glucose transport genes SWEET1 and SFP2. Consistently, in vitro cultivation experiments confirmed that glucose enhances the transition of dormant buds into ramets within two days. Collectively, cytokinin and glucose act as dormant breakers, with cytokinin mainly driving secondary stolon formation and glucose promoting ramet generation. This study improved our understanding of stolon patterning and bud development in the sympodial stolon of strawberries.

Comparative analysis of the microbiomes of strawberry wild species Fragaria nilgerrensis and cultivated variety Akihime using amplicon-based next-generation sequencing

Fragaria nilgerrensis is a wild strawberry species widely distributed in southwest China and has strong ecological adaptability. Akihime (F. × ananassa Duch. cv. Akihime) is one of the main cultivated strawberry varieties in China and is prone to infection with a variety of diseases. In this study, high-throughput sequencing was used to analyze and compare the soil and root microbiomes of F. nilgerrensis and Akihime. Results indicate that the wild species F. nilgerrensis showed higher microbial diversity in nonrhizosphere soil and rhizosphere soil and possessed a more complex microbial network structure compared with the cultivated variety Akihime. Genera such as Bradyrhizobium and Anaeromyxobacter, which are associated with nitrogen fixation and ammonification, and Conexibacter, which is associated with ecological toxicity resistance, exhibited higher relative abundances in the rhizosphere and nonrhizosphere soil samples of F. nilgerrensis compared with those of Akihime. Meanwhile, the ammonia-oxidizing archaea Candidatus Nitrososphaera and Candidatus Nitrocosmicus showed the opposite tendencies. We also found that the relative abundances of potential pathogenic genera and biocontrol bacteria in the Akihime samples were higher than those in the F. nilgerrensis samples. The relative abundances of Blastococcus, Nocardioides, Solirubrobacter, and Gemmatimonas, which are related to pesticide degradation, and genus Variovorax, which is associated with root growth regulation, were also significantly higher in the Akihime samples than in the F. nilgerrensis samples. Moreover, the root endophytic microbiomes of both strawberry species, especially the wild F. nilgerrensis, were mainly composed of potential biocontrol and beneficial bacteria, making them important sources for the isolation of these bacteria. This study is the first to compare the differences in nonrhizosphere and rhizosphere soils and root endogenous microorganisms between wild and cultivated strawberries. The findings have great value for the research of microbiomes, disease control, and germplasm innovation of strawberry.

Whole-genome resequencing reveals genetic diversity and adaptive evolution in Chinese honeybee (Apis cerana cerana) in Guizhou, China

Introduction: Guizhou Province, characterized by complex and diverse geographic and climatic environments, has rich genetic resources for the Chinese honeybee (Apis cerana cerana) and is one of the main bee-producing areas in China. However, research on the genetic diversity of Chinese honeybee in the Guizhou region is very limited, despite implications for conservation of biodiversity. Methods: In this study, we analyzed the genetic diversity, differentiation, and selection signals based on 116 Chinese honeybees from 12 regions in Guizhou Province using whole-genome sequencing. Results: We identified 1,400,430 high-quality SNPs across all samples. A population structure analysis revealed two independent genetic subgroups of Chinese honeybees in Guizhou, a Yunnan-Guizhou Plateau population in western Guizhou and a hilly-mountainous population in eastern Guizhou. The average nucleotide diversity (Pi) ranged from 0.00138 to 0.00161 and average expected heterozygosity (He) ranged from 0.2592 to 0.2604. The average genetic differentiation index (F ST) for Chinese honeybees in pairwise comparisons of 12 regions ranged from 0.0094 to 0.0293. There was clear genetic differentiation between the western plateau and the eastern hilly mountainous areas of Guizhou; however, F ST values between the eastern and western populations ranged from 0.0170 to 0.0293, indicating a low degree of differentiation. A genome-wide scan revealed a number of genes under selection in the Yunnan-Guizhou Plateau environment. These genes were related to growth and development, reproduction, and cold resistance, and several candidate genes involved in environmental adaptation were identified, including CTR, MAPK, MAST, HSF, and MKKK. Discussion: The results of the present study provide important theoretical bases for the conservation, evaluation, development, and utilization of genetic resources for Chinese honeybees in the Guizhou region and for further investigations of environmental adaptation and underlying mechanisms in the species.

Genomic signatures of strawberry domestication and diversification

Cultivated strawberry (Fragaria × ananassa) has a brief history of less than 300 yr, beginning with the hybridization of octoploids Fragaria chiloensis and Fragaria virginiana. Here we explored the genomic signatures of early domestication and subsequent diversification for different climates using whole-genome sequences of 289 wild, heirloom, and modern varieties from two major breeding programs in the United States. Four nonadmixed wild octoploid populations were identified, with recurrent introgression among the sympatric populations. The proportion of F. virginiana ancestry increased by 20% in modern varieties over initial hybrids, and the proportion of F. chiloensis subsp. pacifica rose from 0% to 3.4%. Effective population size rapidly declined during early breeding. Meanwhile, divergent selection for distinct environments reshaped wild allelic origins in 21 out of 28 chromosomes. Overlapping divergent selective sweeps in natural and domesticated populations revealed 16 convergent genomic signatures that may be important for climatic adaptation. Despite 20 breeding cycles since initial hybridization, more than half of loci underlying yield and fruit size are still not under artificial selection. These insights add clarity to the domestication and breeding history of what is now the most widely cultivated fruit in the world.

Selective sweep and GWAS provide insights into adaptive variation of Populus cathayana leaves

Leaf morphology plays a crucial role in predicting the productivity and environmental adaptability of forest trees, making it essential to understand the genetic mechanisms behind leaf variation. In natural populations of Populus cathayana, leaf morphology exhibits rich intraspecific variation due to long-term selection. However, there have been no studies that systematically reveal the genetic mechanisms of leaf variation in P. cathayana. To fill this gap and enhance our understanding of leaf variation in P. cathayana, we collected nine leaf traits from the P. cathayana natural population, consisting of 416 accessions, and conducted the preliminary classification of leaf types with four categories. Subsequently, we conducted an analysis of selective sweep and genome-wide association studies (GWAS) to uncover the genetic basis of leaf traits variation. Most of the leaf traits displayed significant correlations, with broad-sense trait heritability ranging from 0.38 to 0.74. In total, three selective sweep methods ultimately identified 278 positively selected candidate regions and 493 genes associated with leaf size. Single-trait and multi-trait GWAS methods detected 13 and 59 genes, respectively. By integrating the results of selective sweep and GWAS, we further identified a total of nine overlapping genes. These genes may play a role in the leaf development process and are closely associated with leaf size. In particular, the gene CBSCBSPB3 (Pca07G009100) located on chromosome 7, was associated with the response to light stimulation. This study will deepen our understanding of the genetic mechanism of leaf adaptive variation in P. cathayana and provide valuable gene resources.

Genomic signals of local adaptation in Picea crassifolia

BACKGROUND

Global climate change poses a grave threat to biodiversity and underscores the importance of identifying the genes and corresponding environmental factors involved in the adaptation of tree species for the purposes of conservation and forestry. This holds particularly true for spruce species, given their pivotal role as key constituents of the montane, boreal, and sub-alpine forests in the Northern Hemisphere.

RESULTS

Here, we used transcriptomes, species occurrence records, and environmental data to investigate the spatial genetic distribution of and the climate-associated genetic variation in Picea crassifolia. Our comprehensive analysis employing ADMIXTURE, principal component analysis (PCA) and phylogenetic methodologies showed that the species has a complex population structure with obvious differentiation among populations in different regions. Concurrently, our investigations into isolation by distance (IBD), isolation by environment (IBE), and niche differentiation among populations collectively suggests that local adaptations are driven by environmental heterogeneity. By integrating population genomics and environmental data using redundancy analysis (RDA), we identified a set of climate-associated single-nucleotide polymorphisms (SNPs) and showed that environmental isolation had a more significant impact than geographic isolation in promoting genetic differentiation. We also found that the candidate genes associated with altitude, temperature seasonality (Bio4) and precipitation in the wettest month (Bio13) may be useful for forest tree breeding.

CONCLUSIONS

Our findings deepen our understanding of how species respond to climate change and highlight the importance of integrating genomic and environmental data in untangling local adaptations.

Selective sweeps identification in distinct groups of cultivated rye (Secale cereale L.) germplasm provides potential candidate genes for crop improvement

BACKGROUND

During domestication and subsequent improvement plants were subjected to intensive positive selection for desirable traits. Identification of selection targets is important with respect to the future targeted broadening of diversity in breeding programmes. Rye (Secale cereale L.) is a cereal that is closely related to wheat, and it is an important crop in Central, Eastern and Northern Europe. The aim of the study was (i) to identify diverse groups of rye accessions based on high-density, genome-wide analysis of genetic diversity within a set of 478 rye accessions, covering a full spectrum of diversity within the genus, from wild accessions to inbred lines used in hybrid breeding, and (ii) to identify selective sweeps in the established groups of cultivated rye germplasm and putative candidate genes targeted by selection.

RESULTS

Population structure and genetic diversity analyses based on high-quality SNP (DArTseq) markers revealed the presence of three complexes in the Secale genus: S. sylvestre, S. strictum and S. cereale/vavilovii, a relatively narrow diversity of S. sylvestre, very high diversity of S. strictum, and signatures of strong positive selection in S. vavilovii. Within cultivated ryes we detected the presence of genetic clusters and the influence of improvement status on the clustering. Rye landraces represent a reservoir of variation for breeding, and especially a distinct group of landraces from Turkey should be of special interest as a source of untapped variation. Selective sweep detection in cultivated accessions identified 133 outlier positions within 13 sweep regions and 170 putative candidate genes related, among others, to response to various environmental stimuli (such as pathogens, drought, cold), plant fertility and reproduction (pollen sperm cell differentiation, pollen maturation, pollen tube growth), and plant growth and biomass production.

CONCLUSIONS

Our study provides valuable information for efficient management of rye germplasm collections, which can help to ensure proper safeguarding of their genetic potential and provides numerous novel candidate genes targeted by selection in cultivated rye for further functional characterisation and allelic diversity studies.

Genetic structure and first genome-wide insights into the adaptation of a wild relative of grapevine, Vitis berlandieri

In grafted plants, such as grapevine, increasing the diversity of rootstocks available to growers is an ideal strategy for helping plants to adapt to climate change. The rootstocks used for grapevine are hybrids of various American Vitis, including V. berlandieri. The rootstocks currently use in vineyards are derived from breeding programs involving very small numbers of parental individuals. We investigated the structure of a natural population of V. berlandieri and the association of genetic diversity with environmental variables. In this study, we collected seeds from 78 wild V. berlandieri plants in Texas after open fertilization. We genotyped 286 individuals to describe the structure of the population, and environmental information collected at the sampling site made it possible to perform genome-environment association analysis (GEA). De novo long-read whole-genome sequencing was performed on V. berlandieri and a STRUCTURE analysis was performed. We identified and filtered 104,378 SNPs. We found that there were two subpopulations associated with differences in elevation, temperature, and rainfall between sampling sites. GEA identified three QTL for elevation and 15 QTL for PCA coordinates based on environmental parameter variability. This original study is the first GEA study to be performed on a population of grapevines sampled in natural conditions. Our results shed new light on rootstock genetics and could open up possibilities for introducing greater diversity into genetic improvement programs for grapevine rootstocks.

Integrated transcriptome and methylome analyses reveal the molecular regulation of drought stress in wild strawberry (Fragaria nilgerrensis)

BACKGROUND

Fragaria nilgerrensis, which is a diploid wild strawberry with excellent drought-resistance, would provide useful candidate genes for improving drought resistance of cultivated strawberry. So far, its molecular regulatory networks involved in drought stress are unclear. We therefore investigated the drought response regulatory networks of F. nilgerrensis based on the integrated analysis of DNA methylation, transcriptome and physiological traits during four time points under drought stress.  RESULTS: The most differentially expressed genes and the physiological changes were found at 8 days (T8) compared with 0 day (T0, control). Methylome analysis revealed slight dynamic changes in genome-wide mC levels under drought conditions, while the most hypomethylated and hypermethylated regions were identified at T4 and T8. Association analysis of the methylome and transcriptome revealed that unexpressed genes exhibited expected hypermethylation levels in mCHG and mCHH contexts, and highly expressed genes exhibited corresponding hypomethylation levels in the gene body, but mCG contexts showed the opposite trend. Then, 835 differentially methylated and expressed genes were identified and grouped into four clustering patterns to characterize their functions. The genes with either negative or positive correlation between methylation and gene expression were mainly associated with kinases, Reactive Oxygen Species (ROS) synthesis, scavenging, and the abscisic acid (ABA) signal pathway. Consistently, weighted gene co-expression network analysis (WGCNA) revealed Hub genes including NCED, CYP707A2, PP2Cs and others that play important roles in the ABA signaling pathway.

CONCLUSION

F. nilgerrensis drought is dominated by ABA-dependent pathways, possibly accompanied by ABA-independent crosstalk. DNA methylation may affect gene expression, but their correlation was more subtle and multiple types of association exist. Maintaining the balance between ROS regeneration and scavenging is an important factor in drought resistance in F. nilgerrensis. These results deepen our understanding of drought resistance and its application in breeding in strawberry plants.

Transcriptomic and Metabolomic Analyses Provide Insights into the Formation of the Peach-like Aroma of Fragaria nilgerrensis Schlecht. Fruits

Fragaria nilgerrensis Schlecht. is a wild diploid strawberry species. The intense peach-like aroma of its fruits makes F. nilgerrensis an excellent resource for strawberry breeding programs aimed at enhancing flavors. However, the formation of the peach-like aroma of strawberry fruits has not been comprehensively characterized. In this study, fruit metabolome and transcriptome datasets for F. nilgerrensis (HA; peach-like aroma) and its interspecific hybrids PA (peach-like aroma) and NA (no peach-like aroma; control) were compared. In total, 150 differentially accumulated metabolites were detected. The K-means analysis revealed that esters/lactones, including acetic acid, octyl ester, δ-octalactone, and δ-decalactone, were more abundant in HA and PA than in NA. These metabolites may be important for the formation of the peach-like aroma of F. nilgerrensis fruits. The significantly enriched gene ontology terms assigned to the differentially expressed genes (DEGs) were fatty acid metabolic process and fatty acid biosynthetic process. Twenty-seven DEGs were predicted to be associated with ester and lactone biosynthesis, including AAT, LOX, AOS, FAD, AIM1, EH, FAH, ADH, and cytochrome P450 subfamily genes. Thirty-five transcription factor genes were predicted to be associated with aroma formation, including bHLH, MYB, bZIP, NAC, AP2, GATA, and TCPfamily members. Moreover, we identified differentially expressed FAD, AOS, and cytochrome P450 family genes and NAC, MYB, and AP2 transcription factor genes that were correlated with δ-octalactone and δ-decalactone. These findings provide key insights into the formation of the peach-like aroma of F. nilgerrensis fruits, with implications for the increased use of wild strawberry resources.