Retrospect and prospect of Nicotiana tabacum genome sequencing.

Investigating plant genomes offers crucial foundational resources for exploring various aspects of plant biology and applications, such as functional genomics and breeding practices. With the development in sequencing and assembly technology, several Nicotiana tabacum genomes have been published. In this paper, we reviewed the progress on N. tabacum genome assembly and quality, from the initial draft genomes to the recent high-quality chromosome-level assemblies. The application of long-read sequencing, optical mapping, and Hi-C technologies has significantly improved the contiguity and completeness of N. tabacum genome assemblies, with the latest assemblies having a contig N50 size over 50 Mb. Despite these advancements, further improvements are still required and possible, particularly on the development of pan-genome and telomere-to-telomere (T2T) genomes. These new genomes will capture the genomic diversity and variations among different N. tabacum cultivars and species, and provide a comprehensive view of the N. tabacum genome structure and gene content, so to deepen our understanding of the N. tabacum genome and facilitate precise breeding and functional genomics.

Advances in Cotton Genomics, Genetics and Breeding.

The cotton is an industrial crop of global significance, providing its fibers for the predominant textile material and its seed accumulating abundant oil and protein for other utilizations [...].

Genome­wide analysis of cotton SCAMP genes and functional characterization of GhSCAMP2 and GhSCAMP4 in salt tolerance.

BACKGROUND: Secretory carrier membrane proteins (SCAMPs) form a family of integral membrane proteins and play a crucial role in mediating exocytosis in both animals and plants. While SCAMP genes have been studied in several plant species, their functions in cotton, particularly in response to abiotic stress, have not yet been reported. RESULTS: In this study, a total of 53 SCAMP genes were identified in G. arboreum, G. raimondii, G. hirsutum, and G. barbadense. These genes were classified into five groups based on a phylogenetic analysis with SCAMPs from Arabidopsis thaliana. The main factor driving the expansion of the SCAMP gene family in G. hirsutum is tandem and segmental duplication events. Using MEME, in addition to the conserved SCAMP domain, we identified 3-13 other domains in each GhSCAMP. The cis-element analysis suggested that GhSCAMPs were widely involved in cotton growth and development, and responses to abiotic stresses. RNA sequencing (RNA-Seq) and reverse transcription-quantitative polymerase chain reaction (RT-qPCR) results showed that most GhSCAMPs were expressed highly in many tissues and had differential expression responses to drought, cold, and heat stresses. Knock-down of GhSCAMP2 and GhSCAMP4 by virus-induced gene silencing (VIGS) lead to a salt-sensitive phenotype and had a lower content of CAT, POD, and SOD. CONCLUSIONS: This study identified SCAMP genes in four cotton species, enhancing our understanding of the potential biological functions of SCAMPs. Additionally, we demonstrated that GhSCAMP2 and GhSCAMP4 positively regulate cotton tolerance to salt stress.

Single-cell landscape of long and short glandular trichomes in Nicotiana tabacum leaves.

Glandular trichomes (GTs) play a crucial role in plant defenses and the synthesis of secondary metabolites. Understanding the developmental trajectory of GTs is essential for unraveling their functional significance and potential applications. Here we established a comprehensive single-cell atlas of Nicotiana tabacum leaves, a model plant for GT studies. The atlas included a total of 40,433 cells and successfully captured both long GTs (LGTs) and short GTs (SGTs) from Nicotiana leaves. The developmental trajectories of these trichomes were delineated, revealing potential disparities in epidermal development. Comparative analysis of Arabidopsis and Nicotiana trichome development indicated limited similarity between Arabidopsis epidermal non-glandular trichomes and Nicotiana LGTs and SGTs, implying the essentiality of studying the genes directly involved in the development of Nicotiana GTs for a proper and comprehensive understanding of GT biology. Overall, our results provide profound insights into the developmental intricacies of the specialized GTs.

Biosynthesis and Signaling of Strigolactones Act Synergistically With That of ABA and JA to Enhance Verticillium dahliae Resistance in Cotton (Gossypium hirsutum L.).

Verticillium wilt (VW) caused by the soil-borne fungal pathogen Verticillium dahliae reduces cotton productivity and quality. Numerous studies have explored the genetic and molecular mechanisms regulating VW resistance in cotton, but the role and mechanism of strigolactone (SL) is still elusive. We investigated the function of SL in cotton's immune response to V. dahliae infection by exogenously applying SL analog, blocking or enhancing biosynthesis of endogenous SLs in combination with comparative transcriptome analysis and by exploring cross-talk between SL and other phytohormones. Silencing GhDWARF27 and applying the SL analog GR24 or overexpressing GhDWARF27 decreased and enhanced V. dahliae resistance, respectively. Transcriptome analysis revealed SL-mediated activation of abscisic acid (ABA) and jasmonic acid (JA) biosynthesis and signaling pathways. Enhanced ABA biosynthesis and signaling led to increased activity of antioxidant enzymes and reduced buildup of excess reactive oxygen species. Enhanced JA biosynthesis and signaling facilitated transcription of JA-dependent disease resistance genes. One of the components of the SL signal transduction pathway, GhD53, was found to interact with GhNCED5 and GhLOX2, the key enzymes of ABA and JA biosynthesis, respectively. We revealed the molecular mechanism underlying SL-enabled V. dahliae resistance and provided potential solutions for improving VW resistance in cotton.

HSP70 and APX1 play important roles in cotton male fertility by mediating ROS homeostasis.

Male sterility is used in the production of hybrid seeds and can improve the breeding efficiency of cotton hybrids. Reactive oxygen species is closely associated with the tapetum and pollen development, but their relationship in cotton male fertility remains unclear. In this study, we comprehensively compared the cytology and proteome of the anthers from an Upland cotton (Gossypium hirsutum) material, Shida 98 (WT), and its nearly-isogenic male sterile line Shida 98A (MS). Cytology indicated delayed PCD in the tapetum and defects in microspores in MS anthers. And further studies revealed disruption of ROS homeostasis. Proteomic analysis identified proteins with differential abundance mainly being related to redox homeostasis, protein folding, and apoptotic signaling pathways. GhAPX1 interacted with GhHSP70 and played a crucial role in the development of cotton anthers. Exogenous application of HSP70 inhibitor increased H2O2 content and decreased the activity of APX1 and pollen viability. The GhAPX1 mutants generated by CRISPR/Cas9-mediated gene editing exhibited premature degradation of the tapetum, significant decrease in pollen viability, and significant increase in H2O2 content. Altogether, our results imply HSP70 and APX1 being the key players jointly regulating male fertility by mediating ROS homeostasis. These results provide insights into the proteins associated with male fertility.

Transcription factor OsbZIP10 modulates rice grain quality by regulating OsGIF1.

Understanding and optimizing the process of grain filling helps the quest to maximize rice (Oryza sativa L.) seed yield and quality, yet the intricate mechanisms at play remain fragmented. Transcription factors (TFs) are major players in the gene networks underlying the grain filling process. Here, we employed grain incomplete filling (OsGIF1)/cell wall invertase 2, a key gene involved in grain filling, to explore its upstream TFs and identified a bZIP family TF, OsbZIP10, to be a transcriptional activator of OsGIF1. Rice grains of the knockouts of OsbZIP10 showed increased white-core rates but lower amylose content (AC), leading to better eating and cooking qualities in all genetic backgrounds investigated, though the impact of mutations in OsbZIP10 on grain weight depended on genetic background. Multi-omics analyses suggested that, in addition to OsGIF1, multiple genes involved in different biological processes contributing to grain filling were targeted by OsbZIP10, including OsAGPS1, a gene encoding the ADP-Glc pyrophosphorylase (AGPase) small subunit, and genes contributing to homeostasis of reactive oxygen species. Distinct genetic make-up was observed in OsbZIP10 between japonica and indica rice varieties, with the majority varieties of each subspecies belonging to two different haplotypes that were closely associated with AC. Overexpressing the haplotype linked to high-AC in the low-AC genetic background increased AC. Overall, this study sheds crucial light on the significance of the OsbZIP10-OsGIF1 module in the determination of rice grain quality, offering a potential avenue for genetic engineering of rice to produce seeds with tailored attributes.

Are miRNAs applicable for balancing crop growth and defense trade-off?

Securing agricultural supplies for the increasing population without negative impacts on environment demands new crop varieties with higher yields, better quality, and stronger stress resilience. But breeding such super crop varieties is restrained by growth-defense (G-D) trade-off. MicroRNAs (miRNAs) are versatile regulators of plant growth and immune responses, with several being demonstrated to simultaneously regulate crop growth and defense against biotic stresses and to balance G-D trade-off. Increasing evidence also links miRNAs to the metabolism and signaling of phytohormones, another type of master regulator of plant growth and defense. Here, we synthesize the reported functions of miRNAs in crop growth, development, and responses to bio-stressors, summarize the regulatory scenarios of miRNAs based on their relationship with target(s), and discuss how miRNAs, particularly those involved in crosstalk with phytohormones, can be applied in balancing G-D trade-off in crops. We also propose several open questions to be addressed for adopting miRNAs in balancing crop G-D trade-off.

Spatiotemporal transcriptomic landscape of rice embryonic cells during seed germination.

Characterizing cellular features during seed germination is crucial for understanding the complex biological functions of different embryonic cells in regulating seed vigor and seedling establishment. We performed spatially enhanced resolution omics sequencing (Stereo-seq) and single-cell RNA sequencing (scRNA-seq) to capture spatially resolved single-cell transcriptomes of germinating rice embryos. An automated cell-segmentation model, employing deep learning, was developed to accommodate the analysis requirements. The spatial transcriptomes of 6, 24, 36, and 48 h after imbibition unveiled both known and previously unreported embryo cell types, including two unreported scutellum cell types, corroborated by in situ hybridization and functional exploration of marker genes. Temporal transcriptomic profiling delineated gene expression dynamics in distinct embryonic cell types during seed germination, highlighting key genes involved in nutrient metabolism, biosynthesis, and signaling of phytohormones, reprogrammed in a cell-type-specific manner. Our study provides a detailed spatiotemporal transcriptome of rice embryo and presents a previously undescribed methodology for exploring the roles of different embryonic cells in seed germination.

Quantitative genomics-enabled selection for simultaneous improvement of lint yield and seed traits in cotton (Gossypium hirsutum L.).

KEY MESSAGE: A Bayesian linkage disequilibrium-based multiple-locus mixed model identified QTLs for fibre, seed and oil traits and predicted breeding worthiness of test lines, enabling their simultaneous improvement in cotton. Improving cotton seed and oil yields has become increasingly important while continuing to breed for higher lint yield. In this study, a novel Bayesian linkage disequilibrium-based multiple-locus mixed model was developed for QTL identification and genomic prediction (GP). A multi-parent population consisting of 256 recombinant inbred lines, derived from four elite cultivars with distinct combinations of traits, was used in the analysis of QTLs for lint percentage, seed index, lint index and seed oil content and their interrelations. All four traits were moderately heritable and correlated but with no large influence of genotype × environment interactions across multiple seasons. Seven to ten major QTLs were identified for each trait with many being adjacent or overlapping for different trait pairs. A fivefold cross-validation of the model indicated prediction accuracies of 0.46-0.62. GP results based on any two-season phenotypes were strongly correlated with phenotypic means of a pooled analysis of three-season experiments (r = 0.83-0.92). When used for selection of improvement in lint, seed and oil yields, GP captured 40-100% of individuals with comparable lint yields of those selected based on the three-season phenotypic results. Thus, this quantitative genomics-enabled approach can not only decipher the genomic variation underlying lint, seed and seed oil traits and their interrelations, but can provide predictions for their simultaneous improvement. We discuss future breeding strategies in cotton that will enhance the entire value of the crop, not just its fibre.

Technology-enabled great leap in deciphering plant genomes.

Plant genomes provide essential and vital basic resources for studying many aspects of plant biology and applications (for example, breeding). From 2000 to 2020, 1,144 genomes of 782 plant species were sequenced. In the past three years (2021-2023), 2,373 genomes of 1,031 plant species, including 793 newly sequenced species, have been assembled, representing a great leap. The 2,373 newly assembled genomes, of which 63 are telomere-to-telomere assemblies and 921 have been generated in pan-genome projects, cover the major phylogenetic clades. Substantial advances in read length, throughput, accuracy and cost-effectiveness have notably simplified the achievement of high-quality assemblies. Moreover, the development of multiple software tools using different algorithms offers the opportunity to generate more complete and complex assemblies. A database named N3: plants, genomes, technologies has been developed to accommodate the metadata associated with the 3,517 genomes that have been sequenced from 1,575 plant species since 2000. We also provide an outlook for emerging opportunities in plant genome sequencing.

A reference genome of Commelinales provides insights into the commelinids evolution and global spread of water hyacinth (Pontederia crassipes).

Commelinales belongs to the commelinids clade, which also comprises Poales that includes the most important monocot species, such as rice, wheat, and maize. No reference genome of Commelinales is currently available. Water hyacinth (Pontederia crassipes or Eichhornia crassipes), a member of Commelinales, is one of the devastating aquatic weeds, although it is also grown as an ornamental and medical plant. Here, we present a chromosome-scale reference genome of the tetraploid water hyacinth with a total length of 1.22 Gb (over 95% of the estimated size) across 8 pseudochromosome pairs. With the representative genomes, we reconstructed a phylogeny of the commelinids, which supported Zingiberales and Commelinales being sister lineages of Arecales and shed lights on the controversial relationship of the orders. We also reconstructed ancestral karyotypes of the commelinids clade and confirmed the ancient commelinids genome having 8 chromosomes but not 5 as previously reported. Gene family analysis revealed contraction of disease-resistance genes during polyploidization of water hyacinth, likely a result of fitness requirement for its role as a weed. Genetic diversity analysis using 9 water hyacinth lines from 3 continents (South America, Asia, and Europe) revealed very closely related nuclear genomes and almost identical chloroplast genomes of the materials, as well as provided clues about the global dispersal of water hyacinth. The genomic resources of P. crassipes reported here contribute a crucial missing link of the commelinids species and offer novel insights into their phylogeny.

Gh4CL20/20A involved in flavonoid biosynthesis is essential for male fertility in cotton (Gossypium hirsutum L.).

Flavonoids have been shown to play an essential role in plant growth and fertility. 4-Coumarate CoA ligase (4CL) is one of the indispensable enzymes involved in the biosynthesis of flavonoids. However, the role of 4CL and flavonoids in impact on cotton fertility is still unknown. In this study, on the basis of identification of an additional Gh4CL gene, Gh4CL20A, by using an updated G. hirsutum genome, we found that Gh4CL20A and its homologous Gh4CL20 were preferentially expressed in petals and stamens. The petals of the loss-of-function Gh4CL20/Gh4CL20A mutant generated by CRISPR/Cas9 gene editing remained white until wilting. Notably, the mutant showed indehiscent anthers, reduced number of pollen grains and pollen viability, leading to male sterility. Histological analysis revealed that abnormal degradation of anther tapetum at the tetrad stage and abnormal pollen grain development at the mature stage caused male sterility of the gene editing mutant. Analysis of the anther transcriptome identified a total of 10574 and 11962 genes up- and down-regulated in the mutant, respectively, compared to the wild-type. GO, KEGG, and WGCNA analyses linked the abnormality of the mutant anthers to the defective flavonoid biosynthetic pathway, leading to decreased activity of 4CL and chalcone isomerase (CHI) and reduced accumulation of flavonoids in the mutant. These results imply a role of Gh4CL20/Gh4CL20A in assuring proper development of cotton anthers by regulating flavonoid metabolism. This study elucidates a molecular mechanism underlying cotton anther development and provides candidate genes for creating cotton male sterile germplasm that has the potential to be used in production of hybrid seeds.

Genetic Mapping and Characterization of Verticillium Wilt Resistance in a Recombinant Inbred Population of Upland Cotton.

Verticillium wilt (VW) is an important and widespread disease of cotton and once established is long-lived and difficult to manage. In Australia, the non-defoliating pathotype of Verticillium dahliae is the most common, and extremely virulent. Breeding cotton varieties with increased VW resistance is the most economical and effective method of controlling this disease and is greatly aided by understanding the genetics of resistance. This study aimed to investigate VW resistance in 240 F7 recombinant inbred lines (RIL) derived from a cross between MCU-5, which has good resistance, and Siokra 1-4, which is susceptible. Using a controlled environment bioassay, we found that resistance based on plant survival or shoot biomass was complex but with major contributions from chromosomes D03 and D09, with genomic prediction analysis estimating a prediction accuracy of 0.73 based on survival scores compared to 0.36 for shoot biomass. Transcriptome analysis of MCU-5 and Siokra 1-4 roots uninfected or infected with V. dahliae revealed that the two cultivars displayed very different root transcriptomes and responded differently to V. dahliae infection. Ninety-nine differentially expressed genes were located in the two mapped resistance regions and so are potential candidates for further identifying the genes responsible for VW resistance.

Population genomic analysis unravels the evolutionary roadmap of pericarp color in rice.

Pigmented rice stands out for its nutritional value and is gaining more and more attention. Wild rice, domesticated red rice, and weedy rice all have a red pericarp and a comprehensive genetic background in terms of the red-pericarp phenotype. We performed population genetic analyses using 5104 worldwide rice accessions, including 2794 accessions with red or black pericarps, 85 of which were newly sequenced in this study. The results suggested an evolutionary trajectory of red landraces originating from wild rice, and the split times of cultivated red and white rice populations were estimated to be within the past 3500 years. Cultivated red rice was found to feralize to weedy rice, and weedy rice could be further re-domesticated to cultivated red rice. A genome-wide association study based on the 2794 accessions with pigmented pericarps revealed several new candidate genes associated with the red-pericarp trait for further functional characterization. Our results provide genomic evidence for the origin of pigmented rice and a valuable genomic resource for genetic investigation and breeding of pigmented rice.

Dissecting the major genetic components underlying cotton lint development.

Numerous genetic loci and several functionally characterized genes have been linked to determination of lint percentage (lint%), one of the most important cotton yield components, but we still know little about the major genetic components underlying lint%. Here, we first linked the genetic loci containing MYB25-like_At and HD1_At to the fiberless seed trait of 'SL1-7-1' and found that MYB25-like_At and HD1_At were very lowly expressed in 'SL1-7-1' ovules during fiber initiation. We then dissected the genetic components involved in determination of lint% using segregating populations derived from crosses of fuzzless mutants and intermediate segregants with different lint%, which not only confirmed the HD1_At locus but identified the HD1_Dt locus as being the major genetic components contributing to fiber initiation and lint%. The segregating populations also allowed us to evaluate the relative contributions of MYB25-like_At, MYB25-like_Dt, HD1_At, and HD1_Dt to lint%. Haplotype analysis of an Upland cotton (Gossypium hirsutum) population with 723 accessions (including 81 fuzzless seed accessions) showed that lint% of the accessions with the LP allele (higher lint%) at MYB25-like_At, MYB25-like_Dt, or HD1_At was significantly higher than that with the lp allele (lower lint%). The lint% of the Upland cotton accessions with 3 or 4 LP alleles at MYB25-like and HD1 was significantly higher than that with 2 LP alleles. The results prompted us to propose a strategy for breeding high-yielding cotton varieties, i.e. pyramiding the LP alleles of MYB25-like and HD1 with new lint% LP alleles without negative impact on seed size and fiber quality.