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.

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.

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.

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.

Poaceae Chloroplast Genome Sequencing: Great Leap Forward in Recent Ten Years.

The first complete chloroplast genome of rice (Oryza sativa) was published in 1989, ushering in a new era of studies of chloroplast genomics in Poaceae. Progresses in Next-Generation Sequencing (NGS) and Third-Generation Sequencing (TGS) technologiesand in the development of genome assembly software, have significantly advanced chloroplast genomics research. Poaceae is one of the most targeted families in chloroplast genome research because of its agricultural, ecological, and economic importance. Over the last 30 years, 2,050 complete chloroplast genome sequences from 40 tribes and 282 genera have been generated, most (97%) of them in the recent ten years. The wealth of data provides the groundwork for studies on species evolution, phylogeny, genetic transformation, and other aspects of Poaceae chloroplast genomes. As a result, we have gained a deeper understanding of the properties of Poaceae chloroplast genomes. Here, we summarize the achievements of the studies of the Poaceae chloroplast genomes and envision the challenges for moving the area ahead.

Acetylation of GhCaM7 enhances cotton resistance to Verticillium dahliae.

Protein lysine acetylation is an important post-translational modification mechanism involved in cellular regulation in eukaryotes. Calmodulin (CaM) is a ubiquitous Ca2+ sensor in eukaryotes and is crucial for plant immunity, but it is so far unclear whether acetylation is involved in CaM-mediated plant immunity. Here, we found that GhCaM7 is acetylated upon Verticillium dahliae (V. dahliae) infection and a positive regulator of V. dahliae resistance. Overexpressing GhCaM7 in cotton and Arabidopsis enhances V. dahliae resistance and knocking-down GhCaM7 makes cotton more susceptible to V. dahliae. Transgenic Arabidopsis plants overexpressing GhCaM7 with mutation at the acetylation site are more susceptible to V. dahliae than transgenics overexpressing the wild-type GhCaM7, implying the importance of the acetylated GhCaM7 in response to V. dahliae infection. Yeast two-hybrid, bimolecular fluorescent complementation, luciferase complementation imaging, and coimmunoprecipitation assays demonstrated interaction between GhCaM7 and an osmotin protein GhOSM34 that was shown to have a positive role in V. dahliae resistance. GhCaM7 and GhOSM34 are co-localized in the cell membrane. Upon V. dahliae infection, the Ca2+ content reduces almost instantly in plants with downregulated GhCaM7 or GhOSM34. Down regulating GhOSM34 enhances accumulation of Na+ and increases cell osmotic pressure. Comparative transcriptomic analyses between cotton plants with an increased or reduced expression level of GhCaM7 and wild-type plants indicate the involvement of jasmonic acid signaling pathways and reactive oxygen species in GhCaM7-enabled disease resistance. Together, these results demonstrate the involvement of CaM protein in the interaction between cotton and V. dahliae, and more importantly, the involvement of the acetylated CaM in the interaction.

Identification of sugar transporter genes and their roles in the pathogenicity of Verticillium dahliae on cotton.

Introduction: Verticillium wilt (VW) caused by Verticillium dahliae is a soil-borne vascular fungal disease that severely affects cotton yield and fiber quality. Sugar metabolism plays an important role in the growth and pathogenicity of V. dahliae. However, limited information is known about the sugar transporter genes and their roles in the growth and pathogenicity of V. dahliae. Method: In this study, genome-wide identification of sugar transporter genes in V. dahliae was conducted and the expression profiles of these genes in response to root exudates from cotton varieties susceptible or resistant to V. dahliae were investigated based on RNA-seq data. Tobacco Rattle Virus-based host-induced gene silencing (TRV-based HIGS) and artificial small interfering RNAs (asiRNAs) were applied to investigate the function of candidate genes involved in the growth and pathogenic process of V. dahliae. Results: A total of 65 putative sugar transporter genes were identified and clustered into 8 Clades. Of the 65 sugar transporter genes, 9 were found to be induced only by root exudates from the susceptible variety, including VdST3 and VdST12 that were selected for further functional study. Silencing of VdST3 or VdST12 in host plants by TRV-based HIGS reduced fungal biomass and enhanced cotton resistance against V. dahliae. Additionally, silencing of VdST12 and VdST3 by feeding asiRNAs targeting VdST12 (asiR815 or asiR1436) and VdST3 (asiR201 or asiR1238) inhibited fungal growth, exhibiting significant reduction in hyphae and colony diameter, with a more significant effect observed for the asiRNAs targeting VdST12. The inhibitory effect of asiRNAs on the growth of V. dahliae was enhanced with the increasing concentration of asiRNAs. Silencing of VdST12 by feeding asiR815+asiR1436 significantly decreased the pathogenicity of V. dahliae. Discussion: The results suggest that VdST3 and VdST12 are sugar transporter genes required for growth and pathogenicity of V. dahliae and that asiRNA is a valuable tool for functional characterization of V. dahliae genes.

Single-cell transcriptomic analysis reveals the developmental trajectory and transcriptional regulatory networks of pigment glands in Gossypium bickii.

Comprehensive utilization of cottonseeds is limited by the presence of pigment glands and its inclusion gossypol. The ideal cotton has glandless seeds but a glanded plant, a trait found in only a few Australian wild cotton species, including Gossypium bickii. Introgression of this trait into cultivated species has proved to be difficult. Understanding the biological processes toward pigment gland morphogenesis and the associated underlying molecular mechanisms will facilitate breeding of cultivated cotton varieties with the trait of glandless seeds and glanded plant. In this study, single-cell RNA sequencing (scRNA-seq) was performed on 12 222 protoplasts isolated from cotyledons of germinating G. bickii seeds 48 h after imbibition. Clustered into 14 distinct clusters unsupervisedly, these cells could be grouped into eight cell populations with the assistance of known cell marker genes. The pigment gland cells were well separated from others and could be separated into pigment gland parenchyma cells, secretory cells, and apoptotic cells. By integrating the pigment gland cell developmental trajectory, transcription factor regulatory networks, and core transcription factor functional validation, we established a model for pigment gland formation. In this model, light and gibberellin were verified to promote the formation of pigment glands. In addition, three novel genes, GbiERF114 (ETHYLENE RESPONSE FACTOR 114), GbiZAT11 (ZINC FINGER OF ARABIDOPSIS THALIANA 11), and GbiNTL9 (NAC TRANSCRIPTION FACTOR-LIKE 9), were found to affect pigment gland formation. Collectively, these findings provide new insights into pigment gland morphogenesis and lay the cornerstone for future cotton scRNA-seq investigations.

VdPT1 Encoding a Neutral Trehalase of Verticillium dahliae Is Required for Growth and Virulence of the Pathogen.

Verticillum dahliae is a soil-borne phytopathogenic fungus causing destructive Verticillium wilt disease. We previously found a trehalase-encoding gene (VdPT1) in V. dahliae being significantly up-regulated after sensing root exudates from a susceptible cotton variety. In this study, we characterized the function of VdPT1 in the growth and virulence of V. dahliae using its deletion-mutant strains. The VdPT1 deletion mutants (ΔVdPT1) displayed slow colony expansion and mycelial growth, reduced conidial production and germination rate, and decreased mycelial penetration ability and virulence on cotton, but exhibited enhanced stress resistance, suggesting that VdPT1 is involved in the growth, pathogenesis, and stress resistance of V. dahliae. Host-induced silencing of VdPT1 in cotton reduced fungal biomass and enhanced cotton resistance against V. dahliae. Comparative transcriptome analysis between wild-type and mutant identified 1480 up-regulated and 1650 down-regulated genes in the ΔVdPT1 strain. Several down-regulated genes encode plant cell wall-degrading enzymes required for full virulence of V. dahliae to cotton, and down-regulated genes related to carbon metabolism, DNA replication, and amino acid biosynthesis seemed to be responsible for the decreased growth of the ΔVdPT1 strain. In contrast, up-regulation of several genes related to glycerophospholipid metabolism in the ΔVdPT1 strain enhanced the stress resistance of the mutated strain.

The evolutionary landscape and expression pattern of plant lincRNAs.

Long intergenic non-coding RNAs (lincRNAs) are important regulators of cellular processes, including development and stress response. Many lincRNAs have been bioinformatically identified in plants, but their evolutionary dynamics and expression characteristics are still elusive. Here, we systematically identified thousands of lincRNAs in 26 plant species, including 6 non-flowering plants, investigated the conservation of the identified lincRNAs in different levels of plant lineages based on sequence and/or synteny homology and explored characteristics of the conserved lincRNAs during plant evolution and their co-expression relationship with protein-coding genes (PCGs). In addition to confirmation of the features well documented in literature for lincRNAs, such as species-specific, fewer exons, tissue-specific expression patterns and less abundantly expressed, we revealed that histone modification signals and/or binding sites of transcription factors were enriched in the conserved lincRNAs, implying their biological functionalities, as demonstrated by identifying conserved lincRNAs related to flower development in both the Brassicaceae and grass families and ancient lincRNAs potentially functioning in meristem development of non-flowering plants. Compared to PCGs, lincRNAs are more likely to be associated with transposable elements (TEs), but with different characteristics in different evolutionary lineages, for instance, the types of TEs and the variable level of association in lincRNAs with different conservativeness. Together, these results provide a comprehensive view on the evolutionary landscape of plant lincRNAs and shed new insights on the conservation and functionality of plant lincRNAs.

Efficient breeding of low glutelin content rice germplasm by simultaneous editing multiple glutelin genes via CRISPR/Cas9.

Chronic kidney disease (CKD) and phenylketonuria (PKU) patients need to eat rice with low glutelin content. Therefore, breeding low glutelin content rice varieties with high yield and delicious taste is one of the major goals of rice breeders due to the high demand for the product. In this study, we designed three sgRNAs targeting nine glutelin genes and generated nine T-DNA-free homozygous editing lines with reduced glutelin content compared with the wild-type due to simultaneous mutation(s) in 5-7 glutelin genes. The glutelin content of two lines is even significantly lower than that of the low glutelin content cultivar, LGC-1. Compared to the wild-type, these low glutelin lines showed similar agronomic traits, including yield components and viscosity properties, and can be used as new varieties or parental materials for further breeding.

Comparative Metabolome and Transcriptome Analysis of Anthocyanin Biosynthesis in White and Pink Petals of Cotton (Gossypium hirsutum L.).

Upland cotton (Gossypium hirsutum L.) is one of the important fiber crops. Cotton flowers usually appear white (or cream-colored) without colored spots at the petal base, and turn pink on the next day after flowering. In this study, using a mutant showing pink petals with crimson spots at their base, we conducted comparative metabolome and transcriptome analyses to investigate the molecular mechanism of coloration in cotton flowers. Metabolic profiling showed that cyanidin-3-O-glucoside and glycosidic derivatives of pelargonidins and peonidins are the main pigments responsible for the coloration of the pink petals of the mutant. A total of 2443 genes differentially expressed (DEGs) between the white and pink petals were identified by RNA-sequencing. Many DEGs are structural genes and regulatory genes of the anthocyanin biosynthesis pathway. Among them, MYB21, UGT88F3, GSTF12, and VPS32.3 showed significant association with the accumulation of cyanidin-3-O-glucoside in the pink petals. Taken together, our study preliminarily revealed the metabolites responsible for the pink petals and the key genes regulating the biosynthesis and accumulation of anthocyanins in the pink petals. The results provide new insights into the biochemical and molecular mechanism underlying anthocyanin biosynthesis in upland cotton.

The Conservation of Long Intergenic Non-Coding RNAs and Their Response to Verticillium dahliae Infection in Cotton.

Long intergenic non-coding RNAs (lincRNAs) have been demonstrated to be vital regulators of diverse biological processes in both animals and plants. While many lincRNAs have been identified in cotton, we still know little about the repositories and conservativeness of lincRNAs in different cotton species or about their role in responding to biotic stresses. Here, by using publicly available RNA-seq datasets from diverse sources, including experiments of Verticillium dahliae (Vd) infection, we identified 24,425 and 17,713 lincRNAs, respectively, in Gossypium hirsutum (Ghr) and G. barbadense (Gba), the two cultivated allotetraploid cotton species, and 6933 and 5911 lincRNAs, respectively, in G. arboreum (Gar) and G. raimondii (Gra), the two extant diploid progenitors of the allotetraploid cotton. While closely related subgenomes, such as Ghr_At and Gba_At, tend to have more conserved lincRNAs, most lincRNAs are species-specific. The majority of the synthetic and transcribed lincRNAs (78.2%) have a one-to-one orthologous relationship between different (sub)genomes, although a few of them (0.7%) are retained in all (sub)genomes of the four species. The Vd responsiveness of lincRNAs seems to be positively associated with their conservation level. The major functionalities of the Vd-responsive lincRNAs seem to be largely conserved amongst Gra, Ghr, and Gba. Many Vd-responsive Ghr-lincRNAs overlap with Vd-responsive QTL, and several lincRNAs were predicted to be endogenous target mimicries of miR482/2118, with a pair being highly conserved between Ghr and Gba. On top of the confirmation of the feature characteristics of the lincRNAs previously reported in cotton and other species, our study provided new insights into the conservativeness and divergence of lincRNAs during cotton evolution and into the relationship between the conservativeness and Vd responsiveness of lincRNAs. The study also identified candidate lincRNAs with a potential role in disease response for functional characterization.

Impact of water deficiency on leaf cuticle lipids and gene expression networks in cotton (Gossypium hirsutum L.).

BACKGROUND: Water deficit (WD) has serious effect on the productivity of crops. Formation of cuticular layer with increased content of wax and cutin on leaf surfaces is closely related to drought tolerance. Identification of drought tolerance associated wax components and cutin monomers and the genes responsible for their biosynthesis is essential for understanding the physiological and genetic mechanisms underlying drought tolerance and improving crop drought resistance. RESULT: In this study, we conducted comparative phenotypic and transcriptomic analyses of two Gossypium hirsutum varieties that are tolerant (XL22) or sensitive (XL17) to drought stress. XL17 consumed more water than XL22, particularly under the WD conditions. WD significantly induced accumulation of most major wax components (C29 and C31 alkanes) and cutin monomers (palmitic acid and stearic acid) in leaves of both XL22 and XL17, although accumulation of the major cutin monomers, i.e., polyunsaturated linolenic acid (C18:3n-3) and linoleic acid (C18:2n-6), were significantly repressed by WD in both XL22 and XL17. According to the results of transcriptome analysis, although many genes and their related pathways were commonly induced or repressed by WD in both XL22 and XL17, WD-induced differentially expressed genes specific to XL22 or XL17 were also evident. Among the genes that were commonly induced by WD were the GhCER1 genes involved in biosynthesis of alkanes, consistent with the observation of enhanced accumulation of alkanes in cotton leaves under the WD conditions. Interestingly, under the WD conditions, several GhCYP86 genes, which encode enzymes catalyzing the omega-hydroxylation of fatty acids and were identified to be the hub genes of one of the co-expression gene modules, showed a different expression pattern between XL22 and XL17 that was in agreement with the WD-induced changes of the content of hydroxyacids or fatty alcohols in these two varieties. CONCLUSION: The results contribute to our comprehending the physiological and genetic mechanisms underlying drought tolerance and provide possible solutions for the difference of drought resistance of different cotton varieties.