Genomic analyses provide insights into sex differentiation of tetraploid strawberry (Fragaria moupinensis).

The strawberry genus, Fragaria, exhibits a wide range of sexual systems and natural ploidy variation. Nearly, all polyploid strawberry species exhibit separate sexes (dioecy). Research has identified the sex-determining sequences as roughly conserved but with repeatedly changed genomic locations across octoploid strawberries. However, it remains unclear whether tetraploid wild strawberries evolved dioecy independently or shared a common origin with octoploid strawberries. In this study, we investigated the sex determinants of F. moupinensis, a dioecious plant with heterogametic females (ZW). Utilizing a combination of haplotype-resolved genome sequencing of the female F. moupinensis, k-mer-based and coverage-based genome-wide association studies (GWAS), and transcriptomic analysis, we discovered a non-recombining, approximately 33.6 kb W-specific region on chromosome 2a. Within this region, only one candidate sex-determining gene (FmoAFT) was identified. Furthermore, an extensive resequencing of the entire Fragaria genus indicated that the W-specific region displays conservative female specificity across all tetraploid species. This observation suggests that dioecy evolved independently in tetraploid and octoploid strawberries. Moreover, employing virus-induced gene silencing (VIGS), we knocked down the expression of the FmoAFT homologue transcript in cultivated strawberries, revealing its potential role in promoting female functions during early carpel development. We also applied DNA affinity purification sequencing (DAP-seq) and yeast one-hybrid assays to identify potential direct targets of FmoAFT. These insights shed new light on the genetic basis and evolutionary history of sex determination in strawberries, thereby facilitating the formulation of strategies to manipulate sex determination in breeding programs.

Evolutionary history and pan-genome dynamics of strawberry (Fragaria spp.).

Strawberry (Fragaria spp.) has emerged as a model system for various fundamental and applied research in recent years. In total, the genomes of five different species have been sequenced over the past 10 y. Here, we report chromosome-scale reference genomes for five strawberry species, including three newly sequenced species' genomes, and genome resequencing data for 128 additional accessions to estimate the genetic diversity, structure, and demographic history of key Fragaria species. Our analyses obtained fully resolved and strongly supported phylogenies and divergence times for most diploid strawberry species. These analyses also uncovered a new diploid species (Fragaria emeiensis Jia J. Lei). Finally, we constructed a pan-genome for Fragaria and examined the evolutionary dynamics of gene families. Notably, we identified multiple independent single base mutations of the MYB10 gene associated with white pigmented fruit shared by different strawberry species. These reference genomes and datasets, combined with our phylogenetic estimates, should serve as a powerful comparative genomic platform and resource for future studies in strawberry.

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.

Genome of Crucihimalaya himalaica, a close relative of Arabidopsis, shows ecological adaptation to high altitude.

Crucihimalaya himalaica, a close relative of Arabidopsis and Capsella, grows on the Qinghai-Tibet Plateau (QTP) about 4,000 m above sea level and represents an attractive model system for studying speciation and ecological adaptation in extreme environments. We assembled a draft genome sequence of 234.72 Mb encoding 27,019 genes and investigated its origin and adaptive evolutionary mechanisms. Phylogenomic analyses based on 4,586 single-copy genes revealed that C. himalaica is most closely related to Capsella (estimated divergence 8.8 to 12.2 Mya), whereas both species form a sister clade to Arabidopsis thaliana and Arabidopsis lyrata, from which they diverged between 12.7 and 17.2 Mya. LTR retrotransposons in C. himalaica proliferated shortly after the dramatic uplift and climatic change of the Himalayas from the Late Pliocene to Pleistocene. Compared with closely related species, C. himalaica showed significant contraction and pseudogenization in gene families associated with disease resistance and also significant expansion in gene families associated with ubiquitin-mediated proteolysis and DNA repair. We identified hundreds of genes involved in DNA repair, ubiquitin-mediated proteolysis, and reproductive processes with signs of positive selection. Gene families showing dramatic changes in size and genes showing signs of positive selection are likely candidates for C. himalaica's adaptation to intense radiation, low temperature, and pathogen-depauperate environments in the QTP. Loss of function at the S-locus, the reason for the transition to self-fertilization of C. himalaica, might have enabled its QTP occupation. Overall, the genome sequence of C. himalaica provides insights into the mechanisms of plant adaptation to extreme environments.

Genomic analysis of field pennycress (Thlaspi arvense) provides insights into mechanisms of adaptation to high elevation.

BACKGROUND: Understanding how organisms evolve and adapt to extreme habitats is of crucial importance in evolutionary ecology. Altitude gradients are an important determinant of the distribution pattern and range of organisms due to distinct climate conditions at different altitudes. High-altitude regions often provide extreme environments including low temperature and oxygen concentration, poor soil, and strong levels of ultraviolet radiation, leading to very few plant species being able to populate elevation ranges greater than 4000 m. Field pennycress (Thlaspi arvense) is a valuable oilseed crop and emerging model plant distributed across an elevation range of nearly 4500 m. Here, we generate an improved genome assembly to understand how this species adapts to such different environments. RESULTS: We sequenced and assembled de novo the chromosome-level pennycress genome of 527.3 Mb encoding 31,596 genes. Phylogenomic analyses based on 2495 single-copy genes revealed that pennycress is closely related to Eutrema salsugineum (estimated divergence 14.32-18.58 Mya), and both species form a sister clade to Schrenkiella parvula and genus Brassica. Field pennycress contains the highest percentage (70.19%) of transposable elements in all reported genomes of Brassicaceae, with the retrotransposon proliferation in the Middle Pleistocene being likely responsible for the expansion of genome size. Moreover, our analysis of 40 field pennycress samples in two high- and two low-elevation populations detected 1,256,971 high-quality single nucleotide polymorphisms. Using three complementary selection tests, we detected 130 candidate naturally selected genes in the Qinghai-Tibet Plateau (QTP) populations, some of which are involved in DNA repair and the ubiquitin system and potential candidates involved in high-altitude adaptation. Notably, we detected a single base mutation causing loss-of-function of the FLOWERING LOCUS C protein, responsible for the transition to early flowering in high-elevation populations. CONCLUSIONS: Our results provide a genome-wide perspective of how plants adapt to distinct environmental conditions across extreme elevation differences and the potential for further follow-up research with extensive data from additional populations and species.

Cultivated hawthorn (Crataegus pinnatifida var. major) genome sheds light on the evolution of Maleae (apple tribe).

Cultivated hawthorn (Crataegus pinnatifida var. major) is an important medicinal and edible plant with a long history of use for health protection in China. Herein, we provide a de novo chromosome-level genome sequence of the hawthorn cultivar "Qiu Jinxing." We assembled an 823.41 Mb genome encoding 40 571 genes and further anchored the 779.24 Mb sequence into 17 pseudo-chromosomes, which account for 94.64% of the assembled genome. Phylogenomic analyses revealed that cultivated hawthorn diverged from other species within the Maleae (apple tribe) at approximately 35.4 Mya. Notably, genes involved in the flavonoid and triterpenoid biosynthetic pathways have been significantly amplified in the hawthorn genome. In addition, our results indicated that the Maleae share a unique ancient tetraploidization event; however, no recent independent whole-genome duplication event was specifically detected in hawthorn. The amplification of non-specific long terminal repeat retrotransposons contributed the most to the expansion of the hawthorn genome. Furthermore, we identified two paleo-sub-genomes in extant species of Maleae and found that these two sub-genomes showed different rearrangement mechanisms. We also reconstructed the ancestral chromosomes of Rosaceae and discussed two possible paleo-polyploid origin patterns (autopolyploidization or allopolyploidization) of Maleae. Overall, our study provides an improved context for understanding the evolution of Maleae species, and this new high-quality reference genome provides a useful resource for the horticultural improvement of hawthorn.

Proteasome-mediated degradation of FRIGIDA modulates flowering time in Arabidopsis during vernalization.

Winter-annual accessions of Arabidopsis thaliana require either exposure to cold stress or vernalization to initiate flowering via FRIGIDA (FRI). FRI acts as a scaffold protein to recruit several chromatin modifiers that epigenetically modify flowering genes. Here, we report that proteasome-mediated FRI degradation regulates flowering during vernalization in Arabidopsis. Our genetic and biochemical experiments demonstrate that FRI directly interacts with the BTB (Bric-a-Brac/Tramtrack/Broad Complex) proteins LIGHT-RESPONSE BTB1 (LRB1) and LRB2 as well as the CULLIN3A (CUL3A) ubiquitin-E3 ligase in vitro and in vivo, leading to proteasomal degradation of FRI during vernalization. The degradation of FRI is accompanied by an increase in the levels of the long noncoding RNA ColdAIR, which reduces the level of histone H3Lys4 trimethylation (H3K4me3) in FLOWERING LOCUS C chromatin to promote flowering. Furthermore, we found that the cold-induced WRKY34 transcription factor binds to the W-box in the promoter region of CUL3A to modulate CUL3A expression. Deficiency of WRKY34 suppressed CUL3A transcription to enhance FRI protein stability and led to late flowering after vernalization. Conversely, overexpression of WRK34 promoted FRI degradation and early flowering through inducing CUL3A accumulation. Together, these data suggest that WRKY34-induced and CUL3A-dependent proteolysis of FRI modulate flowering in response to vernalization.

Phylogeography of Thlaspi arvense (Brassicaceae) in China Inferred from Chloroplast and Nuclear DNA Sequences and Ecological Niche Modeling.

Thlaspi arvense is a well-known annual farmland weed with worldwide distribution, which can be found from sea level to above 4000 m high on the Qinghai-Tibetan Plateau (QTP). In this paper, a phylogeographic history of T. arvense including 19 populations from China was inferred by using three chloroplast (cp) DNA segments (trnL-trnF, rpl32-trnL and rps16) and one nuclear (n) DNA segment (Fe-regulated transporter-like protein, ZIP). A total of 11 chloroplast haplotypes and six nuclear alleles were identified, and haplotypes unique to the QTP were recognized (C4, C5, C7 and N4). On the basis of molecular dating, haplotypes C4, C5 and C7 have separated from others around 1.58 Ma for cpDNA, which corresponds to the QTP uplift. In addition, this article suggests that the T. arvense populations in China are a mixture of diverged subpopulations as inferred by hT/vT test (hT ≤ vT, cpDNA) and positive Tajima's D values (1.87, 0.05 < p < 0.10 for cpDNA and 3.37, p < 0.01 for nDNA). Multimodality mismatch distribution curves and a relatively large shared area of suitable environmental conditions between the Last Glacial Maximum (LGM) as well as the present time recognized by MaxEnt software reject the sudden expansion population model.

Genome-wide and molecular evolution analysis of the subtilase gene family in Vitis vinifera.

BACKGROUND: Vitis vinifera (grape) is one of the most economically significant fruit crops in the world. The availability of the recently released grape genome sequence offers an opportunity to identify and analyze some important gene families in this species. Subtilases are a group of subtilisin-like serine proteases that are involved in many biological processes in plants. However, no comprehensive study incorporating phylogeny, chromosomal location and gene duplication, gene organization, functional divergence, selective pressure and expression profiling has been reported so far for the grape. RESULTS: In the present study, a comprehensive analysis of the subtilase gene family in V. vinifera was performed. Eighty subtilase genes were identified. Phylogenetic analyses indicated that these subtilase genes comprised eight groups. The gene organization is considerably conserved among the groups. Distribution of the subtilase genes is non-random across the chromosomes. A high proportion of these genes are preferentially clustered, indicating that tandem duplications may have contributed significantly to the expansion of the subtilase gene family. Analyses of divergence and adaptive evolution show that while purifying selection may have been the main force driving the evolution of grape subtilases, some of the critical sites responsible for the divergence may have been under positive selection. Further analyses of real-time PCR data suggested that many subtilase genes might be important in the stress response and functional development of plants. CONCLUSIONS: Tandem duplications as well as purifying and positive selections have contributed to the functional divergence of subtilase genes in V. vinifera. The data may contribute to a better understanding of the grape subtilase gene family.

Uncovering fruit flavor and genetic diversity across diploid wild Fragaria species via comparative metabolomics profiling.

Wild Fragaria resources exhibit extensive genetic diversity and desirable edible traits, such as high soluble solid content and flavor compounds. However, specific metabolites in different wild strawberry fruits remain unknown. In this study, we characterized 1008 metabolites covering 11 subclasses among 13 wild diploid resources representing eight species, including F. vesca, F. nilgerrensis, F. viridis, F. nubicola, F. pentaphylla, F. mandschurica, F. chinensis, and F. emeiensis. Fifteen potential metabolite biomarkers were identified to distinguish fruit flavors among the 13 diploid wild Fragaria accessions. A total of nine distinct modules were employed to explore key metabolites related to fruit quality through weighted gene co-expression module analysis, with significant enrichment in amino acid biosynthesis pathway. Notably, the identified significantly different key metabolites highlighted the close association of amino acids, sugars, and anthocyanins with flavor formation. These findings offer valuable resources for improving fruit quality through metabolome-assisted breeding.

A chromosome-level genome assembly of the Knoxia roxburghii (Rubiaceae).

Knoxia roxburghii is a well-known medicinal plant that is widely distributed in southern China and Southeast Asia. Its dried roots, known as hongdaji in traditional Chinese medicine, are used to treat a range of diseases, including cancers, carbuncles, and ascites. In this study, we report a de novo chromosome-level genome sequence for this diploid plant, which has a length of approximately 446.30 Mb with a contig N50 size of 42.26 Mb and scaffold N50 size of 44.38 Mb. Approximately 99.78% of the assembled sequences were anchored to 10 pseudochromosomes and 3 gapless assembled chromosomes were included in this assembly. A total of 24,507 genes were annotated, along with 68.92% of repetitive elements. Overall, our results will facilitate further active component biosynthesis for K. roxburghii and provide insights for future functional genomic studies and DNA-informed breeding.

Integrative analysis of the metabolome and transcriptome reveals the potential mechanism of fruit flavor formation in wild hawthorn (Crataegus chungtienensis).

Hawthorns are important medicinal and edible plants with a long history of health protection in China. Besides cultivated hawthorn, other wild hawthorns may also have excellent medicinal and edible value, such as Crataeguschungtienensis, an endemic species distributed in the Southwest of China. In this study, by integrating the flavor-related metabolome and transcriptome data of the ripening fruit of C. chungtienensis, we have developed an understanding of the formation of hawthorn fruit quality. The results show that a total of 849 metabolites were detected in the young and mature fruit of C. chungtienensis, of which flavonoids were the most detected metabolites. Among the differentially accumulated metabolites, stachyose, maltotetraose and cis-aconitic acid were significantly increased during fruit ripening, and these may be important metabolites affecting fruit flavor change. Moreover, several flavonoids and terpenoids were reduced after fruit ripening compared with young fruit. Therefore, using the unripe fruit of C. chungtienensis may allow us to obtain more medicinal active ingredients such as flavonoids and terpenoids. Furthermore, we screened out some differentially expressed genes (DEGs) related to fruit quality formation, which had important relationships with differentially accumulated sugars, acids, flavonoids and terpenoids. Our study provides new insights into flavor formation in wild hawthorn during fruit development and ripening, and at the same time this study lays the foundation for the improvement of hawthorn fruit flavor.

Comparative Genomic and Phylogenetic Analysis of Chloroplast Genomes of Hawthorn (Crataegus spp.) in Southwest China.

The hawthorns (Crataegus spp.) are widely distributed and famous for their edible and medicinal values. There are ∼18 species and seven varieties of hawthorn in China distributed throughout the country. We now report the chloroplast genome sequences from C. scabrifolia, C. chungtienensis and C. oresbia, from the southwest of China and compare them with the previously released six species in Crataegus and four species in Rosaceae. The chloroplast genome structure of Crataegus is typical and can be divided into four parts. The genome sizes are between 159,654 and 159,898bp. The three newly sequenced chloroplast genomes encode 132 genes, including 85 protein-coding genes, 37 tRNA genes, and eight rRNA genes. Comparative analysis of the chloroplast genomes revealed six divergent hotspot regions, including ndhA, rps16-trnQ-UUG, ndhF-rpl32, rps16-psbK, trnR-UCU-atpA and rpl32-trnL-UAG. According to the correlation and co-occurrence analysis of repeats with indels and SNPs, the relationship between them cannot be ignored. The phylogenetic tree constructed based on the complete chloroplast genome and intergenic region sequences indicated that C. scabrifolia has a different origin from C. chungtienensis and C. oresbia. We support the placement of C. hupehensis, C. cuneata, C. scabrifolia in C. subg. Crataegus and C. kansuensis, C. oresbia, C. kansuensis in C. subg. Sanguineae. In addition, based on the morphology, geographic distribution and phylogenetic relationships of C. chungtienensis and C. oresbia, we speculate that these two species may be the same species. In conclusion, this study has enriched the chloroplast genome resources of Crataegus and provided valuable information for the phylogeny and species identification of this genus.

How to survive in the world's third poplar: Insights from the genome of the highest altitude woody plant, Hippophae tibetana (Elaeagnaceae).

Hippophae tibetana (Tibetan sea-buckthorn) is one of the highest distributed woody plants in the world (3,000-5,200 meters a.s.l.). It is characterized by adaptation to extreme environment and important economic values. Here, we combined PacBio Hifi platform and Hi-C technology to assemble a 1,452.75 Mb genome encoding 33,367 genes with a Contig N50 of 74.31 Mb, and inferred its sexual chromosome. Two Hippophae-specific whole-genome duplication events (18.7-21.2 million years ago, Ma; 28.6-32.4 Ma) and long terminal repeats retroelements (LTR-RTs) amplifications were detected. Comparing with related species at lower altitude, Ziziphus jujuba (<1, 700 meters a.s.l.), H. tibetana had some significantly rapid evolving genes involved in adaptation to high altitude habitats. However, comparing with Hippophae rhamnoides (<3, 700 meters a.s.l.), no rapid evolving genes were found except microtubule and microtubule-based process genes, H. tibetana has a larger genome, with extra 2, 503 genes (7.5%) and extra 680.46 Mb transposable elements (TEs) (46.84%). These results suggest that the changes in the copy number and regulatory pattern of genes play a more important role for H. tibetana adapting to more extreme and variable environments at higher altitude by more TEs and more genes increasing genome variability and expression plasticity. This suggestion was supported by two findings: nitrogen-fixing genes of H. tibetana having more copies, and intact TEs being significantly closer genes than fragmentary TEs. This study provided new insights into the evolution of alpine plants.

Transcriptional memory of gene expression across generations participates in transgenerational plasticity of field pennycress in response to cadmium stress.

Transgenerational plasticity (TGP) occurs when maternal environments influence the expression of traits in offspring, and in some cases may increase fitness of offspring and have evolutionary significance. However, little is known about the extent of maternal environment influence on gene expression of offspring, and its relationship with trait variations across generations. In this study, we examined TGP in the traits and gene expression of field pennycress (Thlaspi arvense) in response to cadmium (Cd) stress. In the first generation, along with the increase of soil Cd concentration, the total biomass, individual height, and number of seeds significantly decreased, whereas time to flowering, superoxide dismutase (SOD) activity, and content of reduced glutathione significantly increased. Among these traits, only SOD activity showed a significant effect of TGP; the offspring of Cd-treated individuals maintained high SOD activity in the absence of Cd stress. According to the results of RNA sequencing and bioinformatic analysis, 10,028 transcripts were identified as Cd-responsive genes. Among them, only 401 were identified as transcriptional memory genes (TMGs) that maintained the same expression pattern under normal conditions in the second generation as in Cd-treated parents in the first generation. These genes mainly participated in Cd tolerance-related processes such as response to oxidative stress, cell wall biogenesis, and the abscisic acid signaling pathways. The results of weighted correlation network analysis showed that modules correlated with SOD activity recruited more TMGs than modules correlated with other traits. The SOD-coding gene CSD2 was found in one of the modules correlated with SOD activity. Furthermore, several TMGs co-expressed with CSD2 were hub genes that were highly connected to other nodes and critical to the network's topology; therefore, recruitment of TMGs in offspring was potentially related to TGP. These findings indicated that, across generations, transcriptional memory of gene expression played an important role in TGP. Moreover, these results provided new insights into the trait evolution processes mediated by phenotypic plasticity.

Pulvinatusia (Brassicaceae), a new cushion genus from China and its systematic position.

The new genus and species Pulvinatusiaxuegulaensis (Brassicaceae) are described and illustrated. The species is a cushion plant collected from Xuegu La, Xizang, China. Its vegetative parts are most similar to those of Arenariabryophylla (Caryophyllaceae) co-occurring in the same region, while its leaves and fruits closely resemble those of Xerodrabapatagonica (Brassicaceae) from Patagonian Argentina and Chile. Family-level phylogenetic analyses based on both nuclear ITS and plastome revealed that it is a member of the tribe Crucihimalayeae, but the infra-/intergeneric relationships within the tribe are yet to be resolved.

Chloroplast phylogeography of Terminalia franchetii (Combretaceae) from the eastern Sino-Himalayan region and its correlation with historical river capture events.

The reorganization of major river drainages in the Sino-Himalayan region of Southwest China was caused primarily by river separation and capture events following the most recent uplift of the Tibetan Plateau (≤ 3.4 Ma). Terminalia franchetii is a shrub or small tree species endemic to the river valleys of this region. Based on a range-wide sampling comprising 28 populations and 258 individuals, we investigated the relationship between the modern phylogeographic structure of T. franchetii and geological changes in drainage patterns, using chloroplast DNA sequences (trnL-F, petL-psbE). T. franchetii was found to harbor high haplotype diversity (h(T)=0.784) but low average within-population diversity (h(S)=0.124). Mismatch distribution and neutrality tests provided no evidence of recent demographic population growth. Two (out of five) population groups identified exhibited a disjunctive distribution of dominant haplotypes between northern and southern valleys, corresponding to the geography of past drainage systems. We conclude that the modern disjunctive distribution of T. franchetii, and associated patterns of cpDNA haplotype variation, result from vicariance caused by several historical river separation and capture events. Overall, our inferred timings of these events (mostly mid-to-late Pleistocene) agree with previous time estimates of drainage re-arrangements in the Sino-Himalayan region.

Phylogeographic structure of Terminalia franchetii (combretaceae) in southwest China and its implications for drainage geological history.

Following the rapid uplift of the Qinghai-Tibetan Plateau, the reorganization of the major river drainages in southwest China was primarily caused by river capture events. However, the impact of these past changes in drainage patterns on the current distribution and genetic structure of the endemic flora of this region remains largely unknown. Here we report a survey of amplified fragment length polymorphism (AFLP) in Terminalia franchetii, an endemic shrub or small tree of the deep and dry-hot river valleys of this region. We surveyed AFLP variation within and among 21 populations (251 individuals) of T. franchetii, distributed disjunctively between northern and southern drainage systems. Using STRUCTURE, principal coordinates analysis, and genetic distance methods, we identified two main population genetic groups (I and II) and four subgroups within the species, as follows: (I) the Upper Jinshajiang Valley (subgroup I((north))) and the Honghe drainage area (subgroup I((south))); (II) the Middle and Lower Jinshajiang and Yalongjiang Valleys (subgroup II((north))) and the Nanpanjiang drainage area (subgroup II((south))). Genetic diversity was lower in group I than in group II. According to the genetic diversity and genetic structure results, we suggest that the modern disjunctive distribution and associated patterns of genetic structure of T. franchetii result from vicariance caused by several historical drainage capture events, involving the separation of the Upper Jinshajiang, Yalongjiang and Daduhe from the Honghe or Nanpanjiang in southwest China.

Spatial genetic and epigenetic structure of Thlaspi arvense (field pennycress) in China.

Thlaspi arvense (field pennycress) is widespread in temperate regions of the northern hemisphere. We estimated the genetic and epigenetic structure of eight T. arvense populations (131 individuals) in China using amplified fragment length polymorphism and methylation-sensitive amplified polymorphism molecular-marker techniques. We detected low diversity at both genetic (mean = 0.03; total = 0.07) and epigenetic (mean = 0.04; total = 0.07) levels, while significant genetic (FST = 0.42, P < 0.001) and epigenetic (FST = 0.32, P < 0.001) divergence was found across the distribution range. Using Mantel testing, we found spatial genetic and epigenetic differentiation, consistent with isolation-by-distance models. We also identified a strong correlation between genetic and epigenetic differentiation (r = 0.7438, P < 0.001), suggesting genetic control of the epigenetic variation. Our results indicate that mating system, natural selection and gene flow events jointly structure spatial patterns of genetic and epigenetic variation. Moreover, epigenetic variation may serve as a basis of natural selection and ecological evolution to enable species to adapt to heterogeneous habitats. Our study provides novel clues for the adaptation of T. arvense.

The complete chloroplast genome and phylogenetic analysis of Saussurea wettsteiniana (Compositae).

Saussurea wettsteiniana is a medicinally important herb endemic to Hengduan Mountains. Here, we report and characterize the complete chloroplast genome sequence of S. wettsteiniana to provide genomic resources useful for future study. The complete chloroplast genome is 152,631 bp in length, consisting of a large single copy and a small single copy of 83,552 bp and 18,637 bp, which were separated by a pair of inverted repeats of 25,221 bp. Totally 133 genes were annotated, including 87 protein-coding genes, 36 tRNA genes, and eight rRNA genes. We also detected two pseudo-genes (ycf1 and rps19). The overall GC content of the whole genome is 37.7%. The phylogenetic tree based on 23 complete plastomes indicated that S. wettsteiniana was closely related to S. involucrata of Compositae.