Genomic variation reveals demographic history and biological adaptation of the ancient relictual, lotus (Nelumbo Adans)

Cell wall remodeling confers plant architecture with distinct wall structure in Nelumbo nucifera

Lotus (Nelumbo nucifera G.) is a perennial aquatic horticultural plant with diverse architectures. Distinct plant architecture (PA) has certain attractive and practical qualities, but its genetic morphogenesis in lotus remains elusive. In this study, we employ genome-wide association analysis (GWAS) for the seven traits of petiole length (PLL), leaf length (LL), leaf width (LW), peduncle length (PLF), flower diameter (FD), petal length (PeL), and petal width (PeW) in 301 lotus accessions. A total of 90 loci are identified to associate with these traits across 4 years of trials. Meanwhile, we perform RNA sequencing (RNA-seq) to analyze the differential expression of the gene (DEG) transcripts between large and small PA (LPA and SPA) of lotus stems (peduncles and petioles). As a result, eight key candidate genes are identified that are all primarily involved in plant cell wall remodeling significantly associated with PA traits by integrating the results of DEGs and GWAS. To verify this result, we compare the cell wall compositions and structures of LPA versus SPA in representative lotus germplasms. Intriguingly, compared with the SPA lotus, the LPA varieties have higher content of cellulose and hemicellulose, but less filling substrates of pectin and lignin. Additionally, we verified longer cellulose chains and higher cellulose crystallinity with less interference in LPA varieties. Taken together, our study illustrates how plant cell wall remodeling affects PA in lotus, shedding light on the genetic architecture of this significant ornamental trait and offering a priceless genetic resource for future genomic-enabled breeding.

Global characterization of somatic mutations and DNA methylation changes during vegetative propagation in strawberries

Somatic mutations arise and accumulate during tissue culture and vegetative propagation, potentially affecting various traits in horticultural crops, but their characteristics are still unclear. Here, somatic mutations in regenerated woodland strawberry derived from tissue culture of shoot tips under different conditions and 12 cultivated strawberry individuals are analyzed by whole genome sequencing. The mutation frequency of single nucleotide variants is significantly increased with increased hormone levels or prolonged culture time in the range of 3.3 × 10-8-3.0 × 10-6 mutations per site. CG methylation shows a stable reduction (0.71%-8.03%) in regenerated plants, and hypoCG-DMRs are more heritable after sexual reproduction. A high-quality haplotype-resolved genome is assembled for the strawberry cultivar "Beni hoppe." The 12 "Beni hoppe" individuals randomly selected from different locations show 4731-6005 mutations relative to the reference genome, and the mutation frequency varies among the subgenomes. Our study has systematically characterized the genetic and epigenetic variants in regenerated woodland strawberry plants and different individuals of the same strawberry cultivar, providing an accurate assessment of somatic mutations at the genomic scale and nucleotide resolution in plants.

Genome-wide analysis of the CCT gene family and functional characterization of SlCCT6 in response to drought stress in tomato

CCT transcription factors are important for photoperiod and abiotic stress regulation in Arabidopsis and rice. However, the CCT gene family has not been reported in tomato. Here, we systematically analyzed this. Thirty-one SlCCT genes were identified and divided into five groups (CMF, TIFY, PRR, S8, and COL), with members unevenly distributed across 12 chromosomes and the third chromosome exhibiting the most distribution. SlCCT was found to interact with an interacting protein (SlGI), transcription factor (MYB), and non-coding RNA (sly-miR156-5p) to jointly regulate the tomato stress response. cis-Acting element analysis of the SlCCT promoter region indicated large stress- and hormone-response elements in this family. Real-time PCR results indicated that SlPRR subfamily genes respond to various abiotic stresses and hormones. Tissue expression analysis revealed that several PRR subfamily genes are highly expressed in flowers, and subcellular localization analysis indicated an SlCCT6 nuclear location. Notably, SlCCT6 expression was significantly induced by drought, and its silencing reduced drought stress tolerance. Moreover, SlCCT6 overexpression enhanced tomato drought resistance by increasing antioxidant enzyme activity and activating stress-related genes, whereas SlCCT6 knockout decreased drought resistance. In conclusion, this provides valuable insights for future research on SlCCT functions.

Heteroplasmy Is Rare in Plant Mitochondria Compared with Plastids despite Similar Mutation Rates

Plant cells harbor two membrane-bound organelles containing their own genetic material-plastids and mitochondria. Although the two organelles coexist and coevolve within the same plant cells, they differ in genome copy number, intracellular organization, and mode of segregation. How these attributes affect the time to fixation or, conversely, loss of neutral alleles is currently unresolved. Here, we show that mitochondria and plastids share the same mutation rate, yet plastid alleles remain in a heteroplasmic state significantly longer compared with mitochondrial alleles. By analyzing genetic variants across populations of the marine flowering plant Zostera marina and simulating organelle allele dynamics, we examine the determinants of allele segregation and allele fixation. Our results suggest that the bottlenecks on the cell population, e.g. during branching or seeding, and stratification of the meristematic tissue are important determinants of mitochondrial allele dynamics. Furthermore, we suggest that the prolonged plastid allele dynamics are due to a yet unknown active plastid partition mechanism. The dissimilarity between plastid and mitochondrial novel allele fixation at different levels of organization may manifest in differences in adaptation processes. Our study uncovers fundamental principles of organelle population genetics that are essential for further investigations of long-term evolution and molecular dating of divergence events.

The genome of Stephania japonica provides insights into the biosynthesis of cepharanthine

Stephania japonica is an early-diverging eudicotyledon plant with high levels of cepharanthine, proven to be effective in curing coronavirus infections. Here, we report a high-quality S. japonica genome. The genome size is 688.52 Mb, and 97.37% sequences anchor to 11 chromosomes. The genome comprises 67.46% repetitive sequences and 21,036 genes. It is closely related to two Ranunculaceae species, which diverged from their common ancestor 55.90-71.02 million years ago (Mya) with a whole-genome duplication 85.59-96.75 Mya. We further reconstruct ancestral karyotype of Ranunculales. Several cepharanthine biosynthesis genes are identified and verified by western blot. Two genes (Sja03G0243 and Sja03G0241) exhibit catalytic activity as shown by liquid chromatography-mass spectrometry. Then, cepharanthine biosynthesis genes, transcription factors, and CYP450 family genes are used to construct a comprehensive network. Finally, we construct an early-diverging eudicotyledonous genome resources (EEGR) database. As the first genome of the Menispermaceae family to be released, this study provides rich resources for genomic studies.

Exploring domestication pattern in lotus: insights from dispensable genome assembly

Lotus (Nelumbo nucifera Gaertn.), an important aquatic plant in horticulture and ecosystems, has been cultivated for more than 7000 years and domesticated into three different subgroups: flower lotus, rhizome lotus, and seed lotus. To explore the domesticated regions of each subgroup, re-sequencing data of 371 lotus accessions collected from the public database were aligned to the genome of 'China-Antique (CA)'. Unmapped reads were used to build the dispensable genome of each subgroup using a metagenome-like assembly strategy. More than 27 Mb of the dispensable genome in these three subgroups and the wild group was assembled, of which 11,761 genes were annotated. Some of the contigs in the dispensable genome were similar to the genomic segments of other lotus accessions other than 'CA'. The annotated genes in each subgroup played essential roles in specific developmental processes. Dissection of selective signals in three cultivated subgroups also demonstrated that subgroup-specific metabolic pathways, such as the brassinosteroids metabolism enrichment in FL, associated with these selected genes in each subgroup and the contigs in dispensable genome nearly located in the domesticated regions of each subgroup, respectively. Our data presented a valuable resource for facilitating lotus genomic studies, complemented the helpful information to the reference genome, and shed light on the selective signals of domesticated subgroups.

Environment as a limiting factor of the historical global spread of mungbean

While the domestication process has been investigated in many crops, the detailed route of cultivation range expansion and factors governing this process received relatively little attention. Here, using mungbean (Vigna radiata var. radiata) as a test case, we investigated the genomes of more than 1000 accessions to illustrate climatic adaptation's role in dictating the unique routes of cultivation range expansion. Despite the geographical proximity between South and Central Asia, genetic evidence suggests mungbean cultivation first spread from South Asia to Southeast, East and finally reached Central Asia. Combining evidence from demographic inference, climatic niche modeling, plant morphology, and records from ancient Chinese sources, we showed that the specific route was shaped by the unique combinations of climatic constraints and farmer practices across Asia, which imposed divergent selection favoring higher yield in the south but short-season and more drought-tolerant accessions in the north. Our results suggest that mungbean did not radiate from the domestication center as expected purely under human activity, but instead, the spread of mungbean cultivation is highly constrained by climatic adaptation, echoing the idea that human commensals are more difficult to spread through the south-north axis of continents.

Comprehensive genome-wide identification and functional characterization of MAPK cascade gene families in Nelumbo

Mitogen-activated protein kinase (MAPK) cascade signaling pathway plays pivotal roles in various plant biological processes. However, systematic study of MAPK cascade gene families is yet to be conducted in lotus. Herein, 198 putative MAPK genes, including 152 MAP3Ks, 15 MKKs, and 31 MPKs genes were identified in Nelumbo. Segmental duplication was identified as the predominant factor driving MAPK cascade gene family expansion in lotus. MAPK cascade genes in N. nucifera and N. lutea shared high degree of sequence homologies, with 84, 9, and 19 homologous MAP3K, MKK, and MPK gene pairs being detected between the two species, respectively, with most genes predominantly undergoing purifying selection. Gene expression profiling indicated that NnMAPK cascade genes were extensively involved in plant development and submergence stress response. Co-expression analysis revealed potential interaction between transcription factors (TFs) and NnMAPK cascade genes in various biological processes. NnMKK showed predicted interactions with multiple NnMAP3K or NnMPK proteins, which suggested that functional diversity of MAPK cascade genes could be as a result of their complex protein interaction mechanisms. This first systematic analysis of MAPK cascade families in lotus provides deeper insights into their evolutionary dynamics and functional properties, which potentially could be crucial for lotus genetic improvement.

Identification of the NAC Transcription Factors and Their Function in ABA and Salinity Response in Nelumbo nucifera

Nelumbo nucifera Gaertn. is an important perennial aquatic herb that has high ornamental, edible, medicinal, and economic value, being widely distributed and used in China. The NAC superfamily (NAM, ATAF1/2, CUC2) plays critical roles in plant growth, development, and response to abiotic and biotic stresses. Though there have been a few reports about NAC genes in lotus, systematic analysis is still relatively lacking. The present study aimed to characterize all the NAC genes in the lotus and obtain better insights on the NnNACs in response to salt stress by depending on ABA signaling. Here, 97 NAC genes were identified by searching the whole lotus genome based on the raw HMM models of the conserved NAM domain and NAC domain. They were characterized by bioinformatics analysis and divided into 18 subgroups based on the phylogenetic tree. Cis-element analysis demonstrated that NAC genes are responsive to biotic and abiotic stresses, light, low temperature, and plant hormones. Meanwhile, NAC genes had tissue expression specificity. qRT-PCR analysis indicated that NAC genes could be upregulated or downregulated by NaCl treatment, ABA, and fluoridone. In addition, NAC016, NAC025, and NAC070, whose encoding genes were significantly induced by NaCl and ABA, were located in the nucleus. Further analysis showed the three NAC proteins had transcriptional activation capabilities. The co-expression network analysis reflected that NAC proteins may form complexes with other proteins to play a role together. Our study provides a theoretical basis for further research to be conducted on the regulatory mechanisms of salinity resistance in the lotus.

Studies on Lotus Genomics and the Contribution to Its Breeding

Lotus (Nelumbo nucifera), under the Nelumbonaceae family, is one of the relict plants possessing important scientific research and economic values. Because of this, much attention has been paid to this species on both its biology and breeding among the scientific community. In the last decade, the genome of lotus has been sequenced, and several high-quality genome assemblies are available, which have significantly facilitated functional genomics studies in lotus. Meanwhile, re-sequencing of the natural and genetic populations along with different levels of omics studies have not only helped to classify the germplasm resources but also to identify the domestication of selected regions and genes controlling different horticultural traits. This review summarizes the latest progress of all these studies on lotus and discusses their potential application in lotus breeding.