Asymmetric gene expression and cell-type-specific regulatory networks in the root of bread wheat revealed by single-cell multiomics analysis

Single-cell atlases reveal leaf cell-type-specific regulation of metal transporters in the hyperaccumulator Sedum alfredii under cadmium stress

Hyperaccumulation in plants is a complex and dynamic biological process. Sedum alfredii, the most studied Cd hyperaccumulator, can accumulate up to 9000 mg kg-1 Cd in its leaves without suffering toxicity. Although several studies have reported the molecular mechanisms of Cd hyperaccumulation, our understanding of the cell-type-specific transcriptional regulation induced by Cd remains limited. In this study, the first full-length transcriptome of S. alfredii was generated using the PacBio Iso-Seq technology. A total of 18,718,513 subreads (39.90 Gb) were obtained, with an average length of 2133 bp. The single-cell RNA sequencing was employed on leaves of S. alfredii grown under Cd stress. A total of 12,616 high-quality single cells were derived from the control and Cd-treatment samples of S. alfredii leaves. Based on cell heterogeneity and the expression profiles of previously reported marker genes, seven cell types with 12 transcriptionally distinct cell clusters were identified, thereby constructing the first single-cell atlas for S. alfredii leaves. Metal transporters such as CAX5, COPT5, ZIP5, YSL7, and MTP1 were up-regulated in different cell types of S. alfredii leaves under Cd stress. The distinctive gene expression patterns of metal transporters indicate special gene regulatory networks underlying Cd tolerance and hyperaccumulation in S. alfredii. Collectively, our findings are the first observation of the cellular and molecular responses of S. alfredii leaves under Cd stress and lay the cornerstone for future hyperaccumulator scRNA-seq investigations.

Chromatin loops gather targets of upstream regulators together for efficient gene transcription regulation during vernalization in wheat

BACKGROUND

Plants respond to environmental stimuli by altering gene transcription that is highly related with chromatin status, including histone modification, chromatin accessibility, and three-dimensional chromatin interaction. Vernalization is essential for the transition to reproductive growth for winter wheat. How wheat reshapes its chromatin features, especially chromatin interaction during vernalization, remains unknown.

RESULTS

Combinatory analysis of gene transcription and histone modifications in winter wheat under different vernalization conditions identifies 17,669 differential expressed genes and thousands of differentially enriched peaks of H3K4me3, H3K27me3, and H3K9ac. We find dynamic gene expression across the vernalization process is highly associated with H3K4me3. More importantly, the dynamic H3K4me3- and H3K9ac-associated chromatin-chromatin interactions demonstrate that vernalization leads to increased chromatin interactions and gene activation. Remarkably, spatially distant targets of master regulators like VRN1 and VRT2 are gathered together by chromatin loops to achieve efficient transcription regulation, which is designated as a "shepherd" model. Furthermore, by integrating gene regulatory network for vernalization and natural variation of flowering time, TaZNF10 is identified as a negative regulator for vernalization-related flowering time in wheat.

CONCLUSIONS

We reveal dynamic gene transcription network during vernalization and find that the spatially distant genes can be recruited together via chromatin loops associated with active histone mark thus to be more efficiently found and bound by upstream regulator. It provides new insights into understanding vernalization and response to environmental stimuli in wheat and other plants.

Molecular Mechanisms of Heterosis and Its Applications in Tree Breeding: Progress and Perspectives

Heterosis, or hybrid vigor, refers to the phenomenon where hybrid progenies outperform their parents in traits such as yield and resistance. This phenomenon has been widely applied in plant breeding. Recent advances in high-throughput genomics have significantly advanced our understanding of heterosis. This review systematically summarizes the genetic, molecular, and epigenetic mechanisms underlying heterosis. Furthermore, we discuss recent advances in predictive methods for heterosis and their applications in improving growth rate, resistance to abiotic stresses, and wood yield in tree species. We also explore the role of tree genomics in unraveling the mechanisms underlying heterosis, emphasizing the potential of integrating high-resolution genomics, single-cell sequencing, and spatial transcriptomics to achieve a comprehensive understanding of heterosis from the molecular to spatial levels. Building on this, CRISPR-based gene-editing technologies can be employed to precisely edit heterotic loci, enabling the study of allele function. Additionally, molecular marker-assisted selection (MAS) can be utilized to identify heterotic loci in parental lines, facilitating the selection of optimal hybrid combinations and significantly reducing the labor and time costs of hybrid breeding. Finally, we review the utilization of heterosis in tree breeding and provide a forward-looking perspective on future research directions, highlighting the potential of integrating multi-omics approaches and emerging gene-editing tools to revolutionize tree hybrid breeding.

Single same-cell multiome for dissecting key plant traits

Understanding molecular dynamics at the single cell level is crucial to understand plant traits. Recently, Liu et al. and Cui et al. reported multiome analysis in the same cell/nucleus to dissect the key plant traits (osmotic stress response and pod development). Their results provide novel insights into pathways and regulatory networks at a single cell resolution.

A novel root hair mutant, srh1, affects root hair elongation and reactive oxygen species levels in wheat

Background

Root hairs are single-celled projections on root surfaces, critical for water and nutrient uptake. Here, we describe the first short root hair mutant in wheat (Triticum aestivum L.), identified in a mutagenized population and termed here short root hair 1 (srh1).

Results

While the srh1 mutant can initiate root hair bulges, lack of subsequent extension results in very short root hairs. Due to its semi-dominant nature, heterozygous lines displayed intermediate root hair lengths compared to wild-type. Bulked segregant analysis in a BC1F3 segregating population genotyped via exome capture sequencing localized the genetic control of this mutant to a region on the long arm of chromosome 3A. Via RNA sequencing and bioinformatic analysis, we identified two promising candidate genes. The first was a respiratory burst oxidase homolog (RBOH) encoding gene TaNOX3-A, orthologous to RBOH genes controlling root hair elongation in rice (OsNOX3) and maize (ZmRTH5), that carries a missense mutation in a conserved region of the predicted protein. RBOHs are membrane bound proteins that produce reactive oxygen species (ROS) which trigger cell wall extensibility, allowing subsequent root hair elongation. Notably, reduced ROS levels were observed in srh1 root hair bulges compared to wild-type. The second candidate was the calreticulin-3 encoding gene TaCRT3-A, located within the wider srh1 interval and whose expression was significantly downregulated in srh1 root tissues.

Conclusions

The identification of a major effect gene controlling wheat root hair morphology provides an entry point for future optimization of root hair architecture best suited to future agricultural environments.

Evolution of plant cell-type-specific cis-regulatory elements

Cis-regulatory elements (CREs) are critical in regulating gene expression, and yet understanding of CRE evolution remains challenging. Here, we constructed a comprehensive single-cell atlas of chromatin accessibility in Oryza sativa, integrating data from 103,911 nuclei representing 126 discrete cell states across nine distinct organs. We used comparative genomics to compare cell-type resolved chromatin accessibility between O. sativa and 57,552 nuclei from four additional grass species (Zea mays, Sorghum bicolor, Panicum miliaceum, and Urochloa fusca). Accessible chromatin regions (ACRs) had different levels of conservation depending on the degree of cell-type specificity. We found a complex relationship between ACRs with conserved noncoding sequences, cell-type specificity, conservation, and tissue-specific switching. Additionally, we found that epidermal ACRs were less conserved compared to other cell types, potentially indicating that more rapid regulatory evolution has occurred in the L1-derived epidermal layer of these species. Finally, we identified and characterized a conserved subset of ACRs that overlapped the repressive histone modification H3K27me3, implicating them as potentially silencer-like CREs maintained by evolution. Collectively, this comparative genomics approach highlights the dynamics of plant cell-type-specific CRE evolution.

Root hairs: an underexplored target for sustainable cereal crop production

To meet the demands of a rising human population, plant breeders will need to develop improved crop varieties that maximize yield in the face of increasing pressure on crop production. Historically, the optimization of crop root architecture has represented a challenging breeding target due to the inaccessibility of the root systems. Root hairs, single cell projections from the root epidermis, are perhaps the most overlooked component of root architecture traits. Root hairs play a central role in facilitating water, nutrient uptake, and soil cohesion. Current root hair architectures may be suboptimal under future agricultural production regimes, coupled with an increasingly variable climate. Here, we review the genetic control of root hair development in the world's three most important crops-rice, maize, and wheat-and highlight conservation of gene function between monocots and the model dicot species Arabidopsis. Advances in genomic techniques including gene editing combined with traditional plant breeding methods have the potential to overcome many inherent issues associated with the design of improved root hair architectures. Ultimately, this will enable detailed characterization of the effects of contrasting root hair morphology strategies on crop yield and resilience, and the development of new varieties better adapted to deliver future food security.

Single-Cell Transcriptomics Applied in Plants

Single-cell RNA sequencing (scRNA-seq) is a high-tech method for characterizing the expression patterns of heterogeneous cells in the same tissue and has changed our evaluation of biological systems by increasing the number of individual cells analyzed. However, the full potential of scRNA-seq, particularly in plant science, has not yet been elucidated. To explore the utilization of scRNA-seq technology in plants, we firstly conducted a comprehensive review of significant scRNA-seq findings in the past few years. Secondly, we introduced the research and applications of scRNA-seq technology to plant tissues in recent years, primarily focusing on model plants, crops, and wood. We then offered five databases that could facilitate the identification of distinct expression marker genes for various cell types. Finally, we analyzed the potential problems, challenges, and directions for applying scRNA-seq in plants, with the aim of providing a theoretical foundation for the better use of this technique in future plant research.

A Quantitative Computational Framework for Allopolyploid Single-Cell Data Integration and Core Gene Ranking in Development

Polyploidization drives regulatory and phenotypic innovation. How the merger of different genomes contributes to polyploid development is a fundamental issue in evolutionary developmental biology and breeding research. Clarifying this issue is challenging because of genome complexity and the difficulty in tracking stochastic subgenome divergence during development. Recent single-cell sequencing techniques enabled probing subgenome-divergent regulation in the context of cellular differentiation. However, analyzing single-cell data suffers from high error rates due to high dimensionality, noise, and sparsity, and the errors stack up in polyploid analysis due to the increased dimensionality of comparisons between subgenomes of each cell, hindering deeper mechanistic understandings. In this study, we develop a quantitative computational framework, called "pseudo-genome divergence quantification" (pgDQ), for quantifying and tracking subgenome divergence directly at the cellular level. Further comparing with cellular differentiation trajectories derived from single-cell RNA sequencing data allows for an examination of the relationship between subgenome divergence and the progression of development. pgDQ produces robust results and is insensitive to data dropout and noise, avoiding high error rates due to multiple comparisons of genes, cells, and subgenomes. A statistical diagnostic approach is proposed to identify genes that are central to subgenome divergence during development, which facilitates the integration of different data modalities, enabling the identification of factors and pathways that mediate subgenome-divergent activity during development. Case studies have demonstrated that applying pgDQ to single-cell and bulk tissue transcriptomic data promotes a systematic and deeper understanding of how dynamic subgenome divergence contributes to developmental trajectories in polyploid evolution.

Single-nucleus RNA and ATAC sequencing analyses provide molecular insights into early pod development of peanut fruit

Peanut (Arachis hypogaea L.) is an important leguminous oil and economic crop that produces flowers aboveground and fruits belowground. Subterranean fruit-pod development, which significantly affects peanut production, involves complex molecular mechanisms that likely require the coordinated regulation of multiple genes in different tissues. To investigate the molecular mechanisms that underlie peanut fruit-pod development, we characterized the anatomical features of early fruit-pod development and integrated single-nucleus RNA-sequencing (snRNA-seq) and single-nucleus assay for transposase-accessible chromatin with sequencing (snATAC-seq) data at the single-cell level. We identified distinct cell types, such as meristem, embryo, vascular tissue, cuticular layer, and stele cells within the shell wall. These specific cell types were used to examine potential molecular changes unique to each cell type during pivotal stages of fruit-pod development. snRNA-seq analyses of differentially expressed genes revealed cell-type-specific insights that were not previously obtainable from transcriptome analyses of bulk RNA. For instance, we identified MADS-box genes that contributes to the formation of parenchyma cells and gravity-related genes that are present in the vascular cells, indicating an essential role for the vascular cells in peg gravitropism. Overall, our single-nucleus analysis provides comprehensive and novel information on specific cell types, gene expression, and chromatin accessibility during the early stages of fruit-pod development. This information will enhance our understanding of the mechanisms that underlie fruit-pod development in peanut and contribute to efforts aimed at improving peanut production.

The root of plant-plant interactions: Belowground special cocktails

Plants interact with each other via a multitude of processes among which belowground communication facilitated by specialized metabolites plays an important but overlooked role. Until now, the exact targets, modes of action, and resulting phenotypes that these metabolites induce in neighboring plants have remained largely unknown. Moreover, positive interactions driven by the release of root exudates are prevalent in both natural field conditions and controlled laboratory environments. In particular, intraspecific positive interactions suggest a genotypic recognition mechanism in addition to non-self perception in plant roots. This review concentrates on recent discoveries regarding how plants interact with one another through belowground signals in intra- and interspecific mixtures. Furthermore, we elaborate on how an enhanced understanding of these interactions can propel the field of agroecology forward.

Progressive meristem and single-cell transcriptomes reveal the regulatory mechanisms underlying maize inflorescence development and sex differentiation

Maize develops separate ear and tassel inflorescences with initially similar morphology but ultimately different architecture and sexuality. The detailed regulatory mechanisms underlying these changes still remain largely unclear. In this study, through analyzing the time-course meristem transcriptomes and floret single-cell transcriptomes of ear and tassel, we revealed the regulatory dynamics and pathways underlying inflorescence development and sex differentiation. We identified 16 diverse gene clusters with differential spatiotemporal expression patterns and revealed biased regulation of redox, programmed cell death, and hormone signals during meristem differentiation between ear and tassel. Notably, based on their dynamic expression patterns, we revealed the roles of two RNA-binding proteins in regulating inflorescence meristem activity and axillary meristem formation. Moreover, using the transcriptional profiles of 53 910 single cells, we uncovered the cellular heterogeneity between ear and tassel florets. We found that multiple signals associated with either enhanced cell death or reduced growth are responsible for tassel pistil suppression, while part of the gibberellic acid signal may act non-cell-autonomously to regulate ear stamen arrest during sex differentiation. We further showed that the pistil-protection gene SILKLESS 1 (SK1) functions antagonistically to the known pistil-suppression genes through regulating common molecular pathways, and constructed a regulatory network for pistil-fate determination. Collectively, our study provides a deep understanding of the regulatory mechanisms underlying inflorescence development and sex differentiation in maize, laying the foundation for identifying new regulators and pathways for maize hybrid breeding and improvement.

Recent advances in exploring transcriptional regulatory landscape of crops

Crop breeding entails developing and selecting plant varieties with improved agronomic traits. Modern molecular techniques, such as genome editing, enable more efficient manipulation of plant phenotype by altering the expression of particular regulatory or functional genes. Hence, it is essential to thoroughly comprehend the transcriptional regulatory mechanisms that underpin these traits. In the multi-omics era, a large amount of omics data has been generated for diverse crop species, including genomics, epigenomics, transcriptomics, proteomics, and single-cell omics. The abundant data resources and the emergence of advanced computational tools offer unprecedented opportunities for obtaining a holistic view and profound understanding of the regulatory processes linked to desirable traits. This review focuses on integrated network approaches that utilize multi-omics data to investigate gene expression regulation. Various types of regulatory networks and their inference methods are discussed, focusing on recent advancements in crop plants. The integration of multi-omics data has been proven to be crucial for the construction of high-confidence regulatory networks. With the refinement of these methodologies, they will significantly enhance crop breeding efforts and contribute to global food security.

Telomere-to-telomere reference genome for Panax ginseng highlights the evolution of saponin biosynthesis

Ginseng (Panax ginseng) is a representative of Chinese traditional medicine, also used worldwide, while the triterpene saponin ginsenoside is the most important effective compound within it. Ginseng is an allotetraploid, with complex genetic background, making the study of its metabolic evolution challenging. In this study, we assembled a telomere-to-telomere ginseng reference genome, constructed of 3.45 Gb with 24 chromosomes and 77 266 protein-coding genes. Additionally, the reference genome was divided into two subgenomes, designated as subgenome A and B. Subgenome A contains a larger number of genes, whereas subgenome B has a general expression advantage, suggesting that ginseng subgenomes experienced asymmetric gene loss with biased gene expression. The two subgenomes separated approximately 6.07 million years ago, and subgenome B shows the closest relation to Panax vietnamensis var. fuscidiscus. Comparative genomics revealed an expansion of gene families associated with ginsenoside biosynthesis in both ginseng subgenomes. Furthermore, both tandem duplications and proximal duplications play crucial roles in ginsenoside biosynthesis. We also screened functional genes identified in previous research and found that some of these genes located in colinear regions between subgenomes have divergence functions, revealing an unbalanced evolution in both subgenomes and the saponin biosynthesis pathway in ginseng. Our work provides important resources for future genetic studies and breeding programs of ginseng, as well as the biosynthesis of ginsenosides.

Multiome in the Same Cell Reveals the Impact of Osmotic Stress on Arabidopsis Root Tip Development at Single-Cell Level

Cell-specific transcriptional regulatory networks (TRNs) play vital roles in plant development and response to environmental stresses. However, traditional single-cell mono-omics techniques are unable to directly capture the relationships and dynamics between different layers of molecular information within the same cells. While advanced algorithm facilitates merging scRNA-seq and scATAC-seq datasets, accurate data integration remains a challenge, particularly when investigating cell-type-specific TRNs. By examining gene expression and chromatin accessibility simultaneously in 16,670 Arabidopsis root tip nuclei, the TRNs are reconstructed that govern root tip development under osmotic stress. In contrast to commonly used computational integration at cell-type level, 12,968 peak-to-gene linkage is captured at the bona fide single-cell level and construct TRNs at an unprecedented resolution. Furthermore, the unprecedented datasets allow to more accurately reconstruct the coordinated changes of gene expression and chromatin states during cellular state transition. During root tip development, chromatin accessibility of initial cells precedes gene expression, suggesting that changes in chromatin accessibility may prime cells for subsequent differentiation steps. Pseudo-time trajectory analysis reveal that osmotic stress can shift the functional differentiation of trichoblast. Candidate stress-related gene-linked cis-regulatory elements (gl-cCREs) as well as potential target genes are also identified, and uncovered large cellular heterogeneity under osmotic stress.

Dissecting the molecular basis of spike traits by integrating gene regulatory networks and genetic variation in wheat

Spike architecture influences both grain weight and grain number per spike, which are the two major components of grain yield in bread wheat (Triticum aestivum L.). However, the complex wheat genome and the influence of various environmental factors pose challenges in mapping the causal genes that affect spike traits. Here, we systematically identified genes involved in spike trait formation by integrating information on genomic variation and gene regulatory networks controlling young spike development in wheat. We identified 170 loci that are responsible for variations in spike length, spikelet number per spike, and grain number per spike through genome-wide association study and meta-QTL analyses. We constructed gene regulatory networks for young inflorescences at the double ridge stage and the floret primordium stage, in which the spikelet meristem and the floret meristem are predominant, respectively, by integrating transcriptome, histone modification, chromatin accessibility, eQTL, and protein-protein interactome data. From these networks, we identified 169 hub genes located in 76 of the 170 QTL regions whose polymorphisms are significantly associated with variation in spike traits. The functions of TaZF-B1, VRT-B2, and TaSPL15-A/D in establishment of wheat spike architecture were verified. This study provides valuable molecular resources for understanding spike traits and demonstrates that combining genetic analysis and developmental regulatory networks is a robust approach for dissection of complex traits.

Studying plant vascular development using single-cell approaches

Vascular cells form a highly complex and heterogeneous tissue. Its composition, function, shape, and arrangement vary with the developmental stage and between organs and species. Understanding the transcriptional regulation underpinning this complexity thus requires a high-resolution technique that is capable of capturing rapid events during vascular cell formation. Single-cell and single-nucleus RNA sequencing (sc/snRNA-seq) approaches provide powerful tools to extract transcriptional information from these lowly abundant and dynamically changing cell types, which allows the reconstruction of developmental trajectories. Here, we summarize and reflect on recent studies using single-cell transcriptomics to study vascular cell types and discuss current and future implementations of sc/snRNA-seq approaches in the field of vascular development.

Advances in the Application of Single-Cell Transcriptomics in Plant Systems and Synthetic Biology

Plants are complex systems hierarchically organized and composed of various cell types. To understand the molecular underpinnings of complex plant systems, single-cell RNA sequencing (scRNA-seq) has emerged as a powerful tool for revealing high resolution of gene expression patterns at the cellular level and investigating the cell-type heterogeneity. Furthermore, scRNA-seq analysis of plant biosystems has great potential for generating new knowledge to inform plant biosystems design and synthetic biology, which aims to modify plants genetically/epigenetically through genome editing, engineering, or re-writing based on rational design for increasing crop yield and quality, promoting the bioeconomy and enhancing environmental sustainability. In particular, data from scRNA-seq studies can be utilized to facilitate the development of high-precision Build-Design-Test-Learn capabilities for maximizing the targeted performance of engineered plant biosystems while minimizing unintended side effects. To date, scRNA-seq has been demonstrated in a limited number of plant species, including model plants (e.g., Arabidopsis thaliana), agricultural crops (e.g., Oryza sativa), and bioenergy crops (e.g., Populus spp.). It is expected that future technical advancements will reduce the cost of scRNA-seq and consequently accelerate the application of this emerging technology in plants. In this review, we summarize current technical advancements in plant scRNA-seq, including sample preparation, sequencing, and data analysis, to provide guidance on how to choose the appropriate scRNA-seq methods for different types of plant samples. We then highlight various applications of scRNA-seq in both plant systems biology and plant synthetic biology research. Finally, we discuss the challenges and opportunities for the application of scRNA-seq in plants.

Genomic asymmetric epigenetic modification of transposable elements is involved in gene expression regulation of allopolyploid Brassica napus

Polyploids are common and have a wide geographical distribution and environmental adaptability. Allopolyploidy may lead to the activation of transposable elements (TE). However, the mechanism of epigenetic modification of TEs in the establishment and evolution of allopolyploids remains to be explored. We focused on the TEs of model allopolyploid Brassica napus (An An Cn Cn ), exploring the TE characteristics of the genome, epigenetic modifications of TEs during allopolyploidization, and regulation of gene expression by TE methylation. In B. napus, approximately 50% of the genome was composed of TEs. TEs increased with proximity to genes, especially DNA transposons. TE methylation levels were negatively correlated with gene expression, and changes in TE methylation levels were able to regulate the expression of neighboring genes related to responses to light intensity and stress, which promoted powerful adaptation of allopolyploids to new environments. TEs can be synergistically regulated by RNA-directed DNA methylation pathways and histone modifications. The epigenetic modification levels of TEs tended to be similar to those of the diploid parents during the genome evolution of B. napus. The TEs of the An subgenome were more likely to be modified, and the imbalance in TE number and epigenetic modification level in the An and Cn subgenomes may lead to the establishment of subgenome dominance. Our study analyzed the characteristics of TE location, DNA methylation, siRNA, and histone modification in B. napus and highlighted the importance of TE epigenetic modifications during the allopolyploidy process, providing support for revealing the mechanism of allopolyploid formation and evolution.

Application and prospects of single-cell and spatial omics technologies in woody plants

Over the past decade, high-throughput sequencing and high-resolution single-cell transcriptome sequencing technologies have undergone rapid development, leading to significant breakthroughs. Traditional molecular biology methods are limited in their ability to unravel cellular-level heterogeneity within woody plant tissues. Consequently, techniques such as single-cell transcriptomics, single-cell epigenetics, and spatial transcriptomics are rapidly gaining popularity in the study of woody plants. In this review, we provide a comprehensive overview of the development of these technologies, with a focus on their applications and the challenges they present in single-cell transcriptome research in woody plants. In particular, we delve into the similarities and differences among the results of current studies and analyze the reasons behind these differences. Furthermore, we put forth potential solutions to overcome the challenges encountered in single-cell transcriptome applications in woody plants. Finally, we discuss the application directions of these techniques to address key challenges in woody plant research in the future.

Deciphering the evolution and complexity of wheat germplasm from a genomic perspective

Bread wheat provides an essential fraction of the daily calorific intake for humanity. Due to its huge and complex genome, progress in studying on the wheat genome is substantially trailed behind those of the other two major crops, rice and maize, for at least a decade. With rapid advances in genome assembling and reduced cost of high-throughput sequencing, emerging de novo genome assemblies of wheat and whole-genome sequencing data are leading to a paradigm shift in wheat research. Here, we review recent progress in dissecting the complex genome and germplasm evolution of wheat since the release of the first high-quality wheat genome. New insights have been gained in the evolution of wheat germplasm during domestication and modern breeding progress, genomic variations at multiple scales contributing to the diversity of wheat germplasm, and complex transcriptional and epigenetic regulations of functional genes in polyploid wheat. Genomics databases and bioinformatics tools meeting the urgent needs of wheat genomics research are also summarized. The ever-increasing omics data, along with advanced tools and well-structured databases, are expected to accelerate deciphering the germplasm and gene resources in wheat for future breeding advances.

Cell type-specific cytonuclear coevolution in three allopolyploid plant species

Cytonuclear disruption may accompany allopolyploid evolution as a consequence of the merger of different nuclear genomes in a cellular environment having only one set of progenitor organellar genomes. One path to reconcile potential cytonuclear mismatch is biased expression for maternal gene duplicates (homoeologs) encoding proteins that target to plastids and/or mitochondria. Assessment of this transcriptional form of cytonuclear coevolution at the level of individual cells or cell types remains unexplored. Using single-cell (sc-) and single-nucleus (sn-) RNAseq data from eight tissues in three allopolyploid species, we characterized cell type-specific variations of cytonuclear coevolutionary homoeologous expression and demonstrated the temporal dynamics of expression patterns across development stages during cotton fiber development. Our results provide unique insights into transcriptional cytonuclear coevolution in plant allopolyploids at the single-cell level.

scPlant: A versatile framework for single-cell transcriptomic data analysis in plants

Single-cell transcriptomics has been fully embraced in plant biological research and is revolutionizing our understanding of plant growth, development, and responses to external stimuli. However, single-cell transcriptomic data analysis in plants is not trivial, given that there is currently no end-to-end solution and that integration of various bioinformatics tools involves a large number of required dependencies. Here, we present scPlant, a versatile framework for exploring plant single-cell atlases with minimum input data provided by users. The scPlant pipeline is implemented with numerous functions for diverse analytical tasks, ranging from basic data processing to advanced demands such as cell-type annotation and deconvolution, trajectory inference, cross-species data integration, and cell-type-specific gene regulatory network construction. In addition, a variety of visualization tools are bundled in a built-in Shiny application, enabling exploration of single-cell transcriptomic data on the fly.