Series-Spatial Transcriptome Profiling of Leafy Head Reveals the Key Transition Leaves for Head Formation in Chinese Cabbage

Advances in Single-Cell Transcriptome Sequencing and Spatial Transcriptome Sequencing in Plants

"Omics" typically involves exploration of the structure and function of the entire composition of a biological system at a specific level using high-throughput analytical methods to probe and analyze large amounts of data, including genomics, transcriptomics, proteomics, and metabolomics, among other types. Genomics characterizes and quantifies all genes of an organism collectively, studying their interrelationships and their impacts on the organism. However, conventional transcriptomic sequencing techniques target population cells, and their results only reflect the average expression levels of genes in population cells, as they are unable to reveal the gene expression heterogeneity and spatial heterogeneity among individual cells, thus masking the expression specificity between different cells. Single-cell transcriptomic sequencing and spatial transcriptomic sequencing techniques analyze the transcriptome of individual cells in plant or animal tissues, enabling the understanding of each cell's metabolites and expressed genes. Consequently, statistical analysis of the corresponding tissues can be performed, with the purpose of achieving cell classification, evolutionary growth, and physiological and pathological analyses. This article provides an overview of the research progress in plant single-cell and spatial transcriptomics, as well as their applications and challenges in plants. Furthermore, prospects for the development of single-cell and spatial transcriptomics are proposed.

Identification of YUC genes associated with leaf wrinkling trait in Tacai variety of Chinese cabbage

Chinese cabbage (Brassica campestris L. ssp. chinensis (L.) Makino) stands as a widely cultivated leafy vegetable in China, with its leaf morphology significantly influencing both quality and yield. Despite its agricultural importance, the precise mechanisms governing leaf wrinkling development remain elusive. This investigation focuses on 'Wutacai', a representative cultivar of the Tacai variety (Brassica campestris L. ssp. chinensis var. rosularis Tsen et Lee), renowned for its distinct leaf wrinkling characteristics. Within the genome of 'Wutacai', we identified a total of 18 YUCs, designated as BraWTC_YUCs, revealing their conservation within the Brassica genus, and their close homology to YUCs in Arabidopsis. Expression profiling unveiled that BraWTC_YUCs in Chinese Cabbage exhibited organ-specific and leaf position-dependent variation. Additionally, transcriptome sequencing data from the flat leaf cultivar 'Suzhouqing' and the wrinkled leaf cultivar 'Wutacai' revealed differentially expressed genes (DEGs) related to auxin during the early phases of leaf development, particularly the YUC gene. In summary, this study successfully identified the YUC gene family in 'Wutacai' and elucidated its potential function in leaf wrinkling trait, to provide valuable insights into the prospective molecular mechanisms that regulate leaf wrinkling in Chinese cabbage.

Cabbage (Brassica oleracea var. capitata) Development in Time: How Differential Parenchyma Tissue Growth Affects Leafy Head Formation

This study aims to categorize the morphological changes during cabbage (B. oleracea ssp. capitata) development, seedling, rosette, folding, and heading, and to elucidate the cellular mechanisms of the leaf curvature, essential for the formation of the leafy head. We followed the growth of two cabbage cultivars with distinct head shapes (round and pointed) and one non-heading collard cultivar; we phenotyped the size and volume of the whole plant as well as the size, shape, and curvature of the leaves during growth. By integrating these phenotypic data, we determined the four vegetative stages for both cabbages. The histological phenotypes of microtome sections from five distinct leaf positions of the rosette, folding, and heading leaves at two timepoints during leaf growth were quantified and revealed variations in cellular parameters among leaf types, between leaf positions, and between the adaxial and abaxial sides. We identified two synergistic cellular mechanisms contributing to the curvature of heading leaves: differential growth across the leaf blade, with increased growth at the leaf's center relative to the margins; and the increased expansion of the spongy parenchyma layer compared to the palisade parenchyma layer, resulting in the direction of the curvature, which is inwards. These two processes together contribute to the typical leafy heads of cabbages.

Opportunities and challenges in the application of single-cell and spatial transcriptomics in plants

Single-cell and spatial transcriptomics have diverted researchers' attention from the multicellular level to the single-cell level and spatial information. Single-cell transcriptomes provide insights into the transcriptome at the single-cell level, whereas spatial transcriptomes help preserve spatial information. Although these two omics technologies are helpful and mature, further research is needed to ensure their widespread applicability in plant studies. Reviewing recent research on plant single-cell or spatial transcriptomics, we compared the different experimental methods used in various plants. The limitations and challenges are clear for both single-cell and spatial transcriptomic analyses, such as the lack of applicability, spatial information, or high resolution. Subsequently, we put forth further applications, such as cross-species analysis of roots at the single-cell level and the idea that single-cell transcriptome analysis needs to be combined with other omics analyses to achieve superiority over individual omics analyses. Overall, the results of this review suggest that combining single-cell transcriptomics, spatial transcriptomics, and spatial element distribution can provide a promising research direction, particularly for plant research.

A Gene Correlation Measurement Method for Spatial Transcriptome Data Based on Partitioning and Distribution

Spatial transcriptome (ST) technology provides both the spatial location and transcriptional profile of spots, as well as tissue images. ST data can be utilized to construct gene regulatory networks, which can help identify gene modules that facilitate the understanding of biological processes such as cell communication. Correlation measurement is the core basis for constructing a gene regulatory network. However, due to the high noise and sparsity in ST data, common correlation measurement methods such as the Pearson correlation coefficient (PCC) and Spearman correlation coefficient (SPCC) are not suitable. In this work, a new gene correlation measurement method called STgcor is proposed. STgcor defines vertexes as spots in a two-dimensional coordinate plane consisting of axes X and Y from the gene pair (X and Y). The joint probability density of Gaussian distribution of the gene pair (X and Y) is calculated to identify and eliminate outliers. To overcome sparsity, the degree, trend, and location of the distribution of vertexes are used to measure the correlation between gene pairs (X, Y). To validate the performance of the STgcor method, it is compared with the PCC and SPCC in a weighted coexpression network analysis method using two ST datasets of breast cancer and prostate cancer. The gene modules identified by these methods are then compared and analyzed. The results show that the STgcor method detects some special gene modules and cancer-related pathways that cannot be detected by the other two methods.

Role of Brassica orphan gene BrLFM on leafy head formation in Chinese cabbage (Brassica rapa)

Brassica orphan gene BrFLM, identified by two allelic mutants, was involved in leafy head formation in Chinese cabbage. Leafy head formation is a unique agronomic trait of Chinese cabbage that determines its yield and quality. In our previous study, an EMS mutagenesis Chinese cabbage mutant library was constructed using the heading Chinese cabbage double haploid (DH) line FT as the wild-type. Here, we screened two extremely similar leafy head deficiency mutants lfm-1 and lfm-2 with geotropic growth leaves from the library to investigate the gene(s) related to leafy head formation. Reciprocal crossing results showed that these two mutants were allelic. We utilized lfm-1 to identify the mutant gene(s). Genetic analysis showed that the mutated trait was controlled by a single nuclear gene Brlfm. Mutmap analysis showed that Brlfm was located on chromosome A05, and BraA05g012440.3C or BraA05g021450.3C were the candidate gene. Kompetitive allele-specific PCR analysis eliminated BraA05g012440.3C from the candidates. Sanger sequencing identified an SNP from G to A at the 271st nucleotide on BraA05g021450.3C. The sequencing of lfm-2 detected another non-synonymous SNP (G to A) located at the 266st nucleotide on BraA05g021450.3C, which verified its function on leafy head formation. We blasted BraA05g021450.3C on database and found that it belongs to a Brassica orphan gene encoding an unknown 13.74 kDa protein, named BrLFM. Subcellular localization showed that BrLFM was located in the nucleus. These findings reveal that BrLFM is involved in leafy head formation in Chinese cabbage.

The Roles of Gibberellins in Regulating Leaf Development

Plant growth and development are correlated with many aspects, including phytohormones, which have specific functions. However, the mechanism underlying the process has not been well elucidated. Gibberellins (GAs) play fundamental roles in almost every aspect of plant growth and development, including cell elongation, leaf expansion, leaf senescence, seed germination, and leafy head formation. The central genes involved in GA biosynthesis include GA20 oxidase genes (GA20oxs), GA3oxs, and GA2oxs, which correlate with bioactive GAs. The GA content and GA biosynthesis genes are affected by light, carbon availability, stresses, phytohormone crosstalk, and transcription factors (TFs) as well. However, GA is the main hormone associated with BR, ABA, SA, JA, cytokinin, and auxin, regulating a wide range of growth and developmental processes. DELLA proteins act as plant growth suppressors by inhibiting the elongation and proliferation of cells. GAs induce DELLA repressor protein degradation during the GA biosynthesis process to control several critical developmental processes by interacting with F-box, PIFS, ROS, SCLl3, and other proteins. Bioactive GA levels are inversely related to DELLA proteins, and a lack of DELLA function consequently activates GA responses. In this review, we summarized the diverse roles of GAs in plant development stages, with a focus on GA biosynthesis and signal transduction, to develop new insight and an understanding of the mechanisms underlying plant development.

Identification of NPF Family Genes in Brassica rapa Reveal Their Potential Functions in Pollen Development and Response to Low Nitrate Stress

Nitrate Transporter 1/Peptide Transporter Family (NPF) genes encode membrane transporters involved in the transport of diverse substrates. However, little is known about the diversity and functions of NPFs in Brassica rapa. In this study, 85 NPFs were identified in B. rapa (BrNPFs) which comprised eight subfamilies. Gene structure and conserved motif analysis suggested that BrNFPs were conserved throughout the genus. Stress and hormone-responsive cis-acting elements and transcription factor binding sites were identified in BrNPF promoters. Syntenic analysis suggested that tandem duplication contributed to the expansion of BrNPFs in B. rapa. Transcriptomic profiling analysis indicated that BrNPF2.6, BrNPF2.15, BrNPF7.6, and BrNPF8.9 were expressed in fertile floral buds, suggesting important roles in pollen development. Thirty-nine BrNPFs were responsive to low nitrate availability in shoots or roots. BrNPF2.10, BrNPF2.19, BrNPF2.3, BrNPF5.12, BrNPF5.16, BrNPF5.8, and BrNPF6.3 were only up-regulated in roots under low nitrate conditions, indicating that they play positive roles in nitrate absorption. Furthermore, many genes were identified in contrasting genotypes that responded to vernalization and clubroot disease. Our results increase understanding of BrNPFs as candidate genes for genetic improvement studies of B. rapa to promote low nitrate availability tolerance and for generating sterile male lines based on gene editing methods.

BrAN contributes to leafy head formation by regulating leaf width in Chinese cabbage (Brassica rapa L. ssp. pekinensis)

Leafy head is an important agronomic trait that determines the yield and quality of Chinese cabbage. The molecular mechanism underlying heading in Chinese cabbage has been the focus of research, and wide leaves are a prerequisite for leafy head formation. In our study, two allelic leafy heading-deficient mutants (lhd1 and lhd2) with narrow leaf phenotypes were screened in an ethyl methanesulfonate mutagenized population from a heading Chinese cabbage double haploid line 'FT'. Genetic analysis revealed that the mutant trait was controlled by a recessive nuclear gene, which was found to be BraA10g000480.3C by MutMap and Kompetitive allele-specific PCR analyses. As BraA10g000480.3C was the ortholog of ANGUSTIFOLIA in Arabidopsis, which has been found to regulate leaf width by controlling cortical microtubule arrangement and pavement cell shape, we named it BrAN. BrAN in mutant lhd1 carried an SNP (G to A) on intron 2 that co-segregated with the mutant phenotype, and disrupted the exon-intron splice junction generating intron retention and a putative truncated protein. BrAN in mutant lhd2 carried an SNP (G to A) on exon 4 leading to a premature stop codon. The ectopic overexpression of BrAN restored normal leaf phenotype due to abnormal cortical microtubule arrangement and pavement cell shape in the Arabidopsis an-t1 mutant. However, transformation of Bran did not rescue the an-t1 phenotype. These results indicate that BrAN contributes to leafy head formation of Chinese cabbage.

Transposable element insertion: a hidden major source of domesticated phenotypic variation in Brassica rapa

Transposable element (TE) is prevalent in plant genomes. However, studies on their impact on phenotypic evolution in crop plants are relatively rare, because systematically identifying TE insertions within a species has been a challenge. Here, we present a novel approach for uncovering TE insertion polymorphisms (TIPs) using pan-genome analysis combined with population-scale resequencing, and we adopt this pipeline to retrieve TIPs in a Brassica rapa germplasm collection. We found that 23% of genes within the reference Chiifu-401-42 genome harbored TIPs. TIPs tended to have large transcriptional effects, including modifying gene expression levels and altering gene structure by introducing new introns. Among 524 diverse accessions, TIPs broadly influenced genes related to traits and acted a crucial role in the domestication of B. rapa morphotypes. As examples, four specific TIP-containing genes were found to be candidates that potentially involved in various climatic conditions, promoting the formation of diverse vegetable crops in B. rapa. Our work reveals the hitherto hidden TIPs implicated in agronomic traits and highlights their widespread utility in studies of crop domestication.