Evolution and subfunctionalization of CIPK6 homologous genes in regulating cotton drought resistance.

The occurrence of whole-genome duplication or polyploidy may promote plant adaptability to harsh environments. Here, we clarify the evolutionary relationship of eight GhCIPK6 homologous genes in upland cotton (Gossypium hirsutum). Gene expression and interaction analyses indicate that GhCIPK6 homologous genes show significant functional changes after polyploidy. Among these, GhCIPK6D1 and GhCIPK6D3 are significantly up-regulated by drought stress. Functional studies reveal that high GhCIPK6D1 expression promotes cotton drought sensitivity, while GhCIPK6D3 expression promotes drought tolerance, indicating clear functional differentiation. Genetic and biochemical analyses confirm the synergistic negative and positive regulation of cotton drought resistance through GhCBL1A1-GhCIPK6D1 and GhCBL2A1-GhCIPK6D3, respectively, to regulate stomatal movement by controlling the directional flow of K+ in guard cells. These results reveal differentiated roles of GhCIPK6 homologous genes in response to drought stress in upland cotton following polyploidy. The work provides a different perspective for exploring the functionalization and subfunctionalization of duplicated genes in response to polyploidization.

The complex hexaploid oil-Camellia genome traces back its phylogenomic history and multi-omics analysis of Camellia oil biosynthesis.

Oil-Camellia (Camellia oleifera), belonging to the Theaceae family Camellia, is an important woody edible oil tree species. The Camellia oil in its mature seed kernels, mainly consists of more than 90% unsaturated fatty acids, tea polyphenols, flavonoids, squalene and other active substances, which is one of the best quality edible vegetable oils in the world. However, genetic research and molecular breeding on oil-Camellia are challenging due to its complex genetic background. Here, we successfully report a chromosome-scale genome assembly for a hexaploid oil-Camellia cultivar Changlin40. This assembly contains 8.80 Gb genomic sequences with scaffold N50 of 180.0 Mb and 45 pseudochromosomes comprising 15 homologous groups with three members each, which contain 135 868 genes with an average length of 3936 bp. Referring to the diploid genome, intragenomic and intergenomic comparisons of synteny indicate homologous chromosomal similarity and changes. Moreover, comparative and evolutionary analyses reveal three rounds of whole-genome duplication (WGD) events, as well as the possible diversification of hexaploid Changlin40 with diploid occurred approximately 9.06 million years ago (MYA). Furthermore, through the combination of genomics, transcriptomics and metabolomics approaches, a complex regulatory network was constructed and allows to identify potential key structural genes (SAD, FAD2 and FAD3) and transcription factors (AP2 and C2H2) that regulate the metabolism of Camellia oil, especially for unsaturated fatty acids biosynthesis. Overall, the genomic resource generated from this study has great potential to accelerate the research for the molecular biology and genetic improvement of hexaploid oil-Camellia, as well as to understand polyploid genome evolution.

Salinity stress alters plant-mediated interactions between above- and below-ground herbivores.

Below-ground herbivory impacts plant development and often induces systemic responses in plants that affect the performance and feeding behavior of above-ground herbivores. Meanwhile, pest-damaged root tissue can enhance a plant's susceptibility to abiotic stress such as salinity. Yet, the extent to which herbivore-induced plant defenses are modulated by such abiotic stress has rarely been studied. In this study, we examine whether root feeding by larvae of the turnip moth, Agrotis segetum (Lepidoptera: Noctuidae) affects the performance of the above-ground, sap-feeding aphid Aphis gossypii (Hemiptera: Aphididae) on cotton, and assess whether those interactions are modulated by salinity stress. In the absence of salinity stress, A. segetum root feeding does not affect A. gossypii development. On the other hand, under intense salinity stress (i.e., 600 mM NaCl), A. segetum root feeding decreases aphid development time by 16.1 % and enhances fecundity by 72.0 %. Transcriptome, metabolome and bioassay trials showed that root feeding and salinity stress jointly trigger the biosynthesis of amino acids in cotton leaves. Specifically, increased titers of valine in leaf tissue relate to an enhanced performance of A. gossypii. Taken together, salinity stress alters the interaction between above- and below-ground feeders by changing amino acid accumulation. Our findings advance our understanding of how plants cope with concurrent biotic and abiotic stressors, and may help tailor plant protection strategies to varying production contexts.

GhHAM regulates GoPGF-dependent gland development and contributes to broad-spectrum pest resistance in cotton.

Cotton is a globally cultivated crop, producing 87% of the natural fiber used in the global textile industry. The pigment glands, unique to cotton and its relatives, serve as a defense structure against pests and pathogens. However, the molecular mechanism underlying gland formation and the specific role of pigment glands in cotton's pest defense are still not well understood. In this study, we cloned a gland-related transcription factor GhHAM and generated the GhHAM knockout mutant using CRISPR/Cas9. Phenotypic observations, transcriptome analysis, and promoter-binding experiments revealed that GhHAM binds to the promoter of GoPGF, regulating pigment gland formation in cotton's multiple organs via the GoPGF-GhJUB1 module. The knockout of GhHAM significantly reduced gossypol production and increased cotton's susceptibility to pests in the field. Feeding assays demonstrated that more than 80% of the cotton bollworm larvae preferred ghham over the wild type. Furthermore, the ghham mutants displayed shorter cell length and decreased gibberellins (GA) production in the stem. Exogenous application of GA3 restored stem cell elongation but not gland formation, thereby indicating that GhHAM controls gland morphogenesis independently of GA. Our study sheds light on the functional differentiation of HAM proteins among plant species, highlights the significant role of pigment glands in influencing pest feeding preference, and provides a theoretical basis for breeding pest-resistant cotton varieties to address the challenges posed by frequent outbreaks of pests.

Robust CRISPR/Mb2Cas12a genome editing tools in cotton plants.

The efficiency and accuracy of the CRISPR/Mb2Cas12a system were demonstrated in cotton, achieving an efficiency of over 90% at target sites. Notably, Mb2Cas12a exhibited significant tolerance under different temperatures ranging from 22°C to 32°C. Additionally, the Mb2Cas12a system revealed effective editing at more relaxed VTTV PAM sites in the cotton genome, which expanded the genome editing range by approximately 2.6-fold than the wide-type LbCas12a. Finally, a multiplex genome editing system was also developed based on Mb2Cas12a, enabling simultaneous editing of eight target sites using a single crRNA cassette.

A dominant negative mutation of GhMYB25-like alters cotton fiber initiation, reducing lint and fuzz.

Cotton (Gossypium hirsutum) fibers, vital natural textile materials, are single-cell trichomes that differentiate from the ovule epidermis. These fibers are categorized as lint (longer fibers useful for spinning) or fuzz (shorter, less useful fibers). Currently, developing cotton varieties with high lint yield but without fuzz remains challenging due to our limited knowledge of the molecular mechanisms underlying fiber initiation. This study presents the identification and characterization of a naturally occurring dominant negative mutation GhMYB25-like_AthapT, which results in a reduced lint and fuzzless phenotype. The GhMYB25-like_AthapT protein exerts its dominant negative effect by suppressing the activity of GhMYB25-like during lint and fuzz initiation. Intriguingly, the negative effect of GhMYB25-like_AthapT could be alleviated by high expression levels of GhMYB25-like. We also uncovered the role of GhMYB25-like in regulating the expression of key genes such as GhPDF2 (PROTODERMAL FACTOR 2), CYCD3; 1 (CYCLIN D3; 1) and PLD (Phospholipase D), establishing its significance as a pivotal transcription factor in fiber initiation. We identified other genes within this regulatory network, expanding our understanding of the determinants of fiber cell fate. These findings offer valuable insights for cotton breeding and contribute to our fundamental understanding of fiber development.

Comprehensive analysis of MAPK gene family in upland cotton (Gossypium hirsutum) and functional characterization of GhMPK31 in regulating defense response to insect infestation.

KEY MESSAGE: The transcriptomic, phenotypic and metabolomic analysis of transgenic plants overexpressing GhMPK31 in upland cotton revealed the regulation of H2O2 burst and the synthesis of defensive metabolites by GhMPK31. Mitogen-activated protein kinases (MAPKs) are a crucial class of protein kinases, which play an essential role in various biological processes in plants. Upland cotton (G. hirsutum) is the most widely cultivated cotton species with high economic value. To gain a better understanding of the role of the MAPK gene family, we conducted a comprehensive analysis of the MAPK gene family in cotton. In this study, a total of 55 GhMPK genes were identified from the whole genome of G. hirsutum. Through an investigation of the expression patterns under diverse stress conditions, we discovered that the majority of GhMPK family members demonstrated robust responses to abiotic stress, pathogen stress and pest stress. Furthermore, the overexpression of GhMPK31 in cotton leaves led to a hypersensitive response (HR)-like cell death phenotype and impaired the defense capability of cotton against herbivorous insects. Transcriptome and metabolomics data analysis showed that overexpression of GhMPK31 enhanced the expression of H2O2-related genes and reduced the accumulation of defensive related metabolites. The direct evidence of GhMPK31 interacting with GhRBOHB (H2O2-generating protein) were found by Y2H, BiFC, and LCI. Therefore, we propose that the increase of H2O2 content caused by overexpression of GhMPK31 resulted in HR-like cell death in cotton leaves while reducing the accumulation of defensive metabolites, ultimately leading to a decrease in the defense ability of cotton against herbivorous insects. This study provides valuable insights into the function of MAPK genes in plant resistance to herbivorous insects.

Precise fine-turning of GhTFL1 by base editing tools defines ideal cotton plant architecture.

BACKGROUND: CRISPR/Cas-derived base editor enables precise editing of target sites and has been widely used for basic research and crop genetic improvement. However, the editing efficiency of base editors at different targets varies greatly. RESULTS: Here, we develop a set of highly efficient base editors in cotton plants. GhABE8e, which is fused to conventional nCas9, exhibits 99.9% editing efficiency, compared to GhABE7.10 with 64.9%, and no off-target editing is detected. We further replace nCas9 with dCpf1, which recognizes TTTV PAM sequences, to broaden the range of the target site. To explore the functional divergence of TERMINAL FLOWER 1 (TFL1), we edit the non-coding and coding regions of GhTFL1 with 26 targets to generate a comprehensive allelic population including 300 independent lines in cotton. This allows hidden pleiotropic roles for GhTFL1 to be revealed and allows us to rapidly achieve directed domestication of cotton and create ideotype germplasm with moderate height, shortened fruiting branches, compact plant, and early-flowering. Further, by exploring the molecular mechanism of the GhTFL1L86P and GhTFL1K53G+S78G mutations, we find that the GhTFL1L86P mutation weakens the binding strength of the GhTFL1 to other proteins but does not lead to a complete loss of GhTFL1 function. CONCLUSIONS: This strategy provides an important technical platform and genetic information for the study and creation of ideal plant architecture.

Plant organellar genomes: much done, much more to do.

Plastids and mitochondria are the only organelles that possess genomes of endosymbiotic origin. In recent decades, advances in sequencing technologies have contributed to a meteoric rise in the number of published organellar genomes, and have revealed greatly divergent evolutionary trajectories. In this review, we quantify the abundance and distribution of sequenced plant organellar genomes across the plant tree of life. We compare numerous genomic features between the two organellar genomes, with an emphasis on evolutionary trajectories, transfers, the current state of organellar genome editing by transcriptional activator-like effector nucleases (TALENs), transcription activator-like effector (TALE)-mediated deaminase, and clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein (Cas), as well as genetic transformation. Finally, we propose future research to understand these different evolutionary trajectories, and genome-editing strategies to promote functional studies and eventually improve organellar genomes.

Construction of Host Plant Insect-Resistance Mutant Library by High-Throughput CRISPR/Cas9 System and Identification of A Broad-Spectrum Insect Resistance Gene.

Insects pose significant challenges in cotton-producing regions. Here, they describe a high-throughput CRISPR/Cas9-mediated large-scale mutagenesis library targeting endogenous insect-resistance-related genes in cotton. This library targeted 502 previously identified genes using 968 sgRNAs, generated ≈2000 T0 plants and achieved 97.29% genome editing with efficient heredity, reaching upto 84.78%. Several potential resistance-related mutants (10% of 200 lines) their identified that may contribute to cotton-insect molecular interaction. Among these, they selected 139 and 144 lines showing decreased resistance to pest infestation and targeting major latex-like protein 423 (GhMLP423) for in-depth study. Overexpression of GhMLP423 enhanced insect resistance by activating the plant systemic acquired resistance (SAR) of salicylic acid (SA) and pathogenesis-related (PR) genes. This activation is induced by an elevation of cytosolic calcium [Ca2+ ]cyt flux eliciting reactive oxygen species (ROS), which their demoted in GhMLP423 knockout (CR) plants. Protein-protein interaction assays revealed that GhMLP423 interacted with a human epidermal growth factor receptor substrate15 (EPS15) protein at the cell membrane. Together, they regulated the systemically propagating waves of Ca2+ and ROS, which in turn induced SAR. Collectively, this large-scale mutagenesis library provides an efficient strategy for functional genomics research of polyploid plant species and serves as a solid platform for genetic engineering of insect resistance.

The chromosome-scale reference genome of mirid bugs (Adelphocoris suturalis) genome provides insights into omnivory, insecticide resistance, and survival adaptation.

BACKGROUND: Adelphocoris suturalis (Hemiptera: Miridae) is a notorious agricultural pest, which causes serious economic losses to a diverse range of agricultural crops around the world. The poor understanding of its genomic characteristics has seriously hindered the establishment of sustainable and environment-friendly agricultural pest management through biotechnology and biological insecticides. RESULTS: Here, we report a chromosome-level assembled genome of A. suturalis by integrating Illumina short reads, PacBio, 10x Chromium, and Hi-C mapping technologies. The resulting 1.29 Gb assembly contains twelve chromosomal pseudomolecules with an N50 of 1.4 and 120.6 Mb for the contigs and scaffolds, respectively, and carries 20,010 protein-coding genes. The considerable size of the A. suturalis genome is predominantly attributed to a high amount of retrotransposons, especially long interspersed nuclear elements (LINEs). Transcriptomic and phylogenetic analyses suggest that A. suturalis-specific candidate effectors, and expansion and expression of gene families associated with omnivory, insecticide resistance and reproductive characteristics, such as digestion, detoxification, chemosensory receptors and long-distance migration likely contribute to its strong environmental adaptability and ability to damage crops. Additionally, 19 highly credible effector candidates were identified and transiently overexpressed in Nicotiana benthamiana for functional assays and potential targeting for insect resistance genetic engineering. CONCLUSIONS: The high-quality genome of A. suturalis provides an important genomic landscape for further investigations into the mechanisms of omnivory, insecticide resistance and survival adaptation, and for the development of integrated management strategies.

Origin, evolution, and diversification of the wall-associated kinase gene family in plants.

KEY MESSAGE: The study of the origin, evolution, and diversification of the wall-associated kinase gene family in plants facilitates their functional investigations in the future. Wall-associated kinases (WAKs) make up one subfamily of receptor-like kinases (RLKs), and function directly in plant cell elongation and responses to biotic and abiotic stresses. The biological functions of WAKs have been extensively characterized in angiosperms; however, the origin and evolutionary history of the WAK family in green plants remain unclear. Here, we performed a comprehensive analysis of the WAK family to reveal its origin, evolution, and diversification in green plants. In total, 1061 WAK genes were identified in 37 species from unicellular algae to multicellular plants, and the results showed that WAK genes probably originated before bryophyte differentiation and were widely distributed in land plants, especially angiosperms. The phylogeny indicated that the land plant WAKs gave rise to five clades and underwent lineage-specific expansion after species differentiation. Cis-acting elements and expression patterns analyses of WAK genes in Arabidopsis and rice demonstrated the functional diversity of WAK genes in these two species. Many gene gains and losses have occurred in angiosperms, leading to an increase in the number of gene copies. The evolutionary trajectory of the WAK family during polyploidization was uncovered using Gossypium species. Our results provide insights into the evolution of WAK genes in green plants, facilitating their functional investigations in the future.

Single-cell resolution analysis reveals the preparation for reprogramming the fate of stem cell niche in cotton lateral meristem.

BACKGROUND: Somatic embryogenesis is a major process for plant regeneration. However, cell communication and the gene regulatory network responsible for cell reprogramming during somatic embryogenesis are still largely unclear. Recent advances in single-cell technologies enable us to explore the mechanism of plant regeneration at single-cell resolution. RESULTS: We generate a high-resolution single-cell transcriptomic landscape of hypocotyl tissue from the highly regenerable cotton genotype Jin668 and the recalcitrant TM-1. We identify nine putative cell clusters and 23 cluster-specific marker genes for both cultivars. We find that the primary vascular cell is the major cell type that undergoes cell fate transition in response to external stimulation. Further developmental trajectory and gene regulatory network analysis of these cell clusters reveals that a total of 41 hormone response-related genes, including LAX2, LAX1, and LOX3, exhibit different expression patterns in the primary xylem and cambium region of Jin668 and TM-1. We also identify novel genes, including CSEF, PIS1, AFB2, ATHB2, PLC2, and PLT3, that are involved in regeneration. We demonstrate that LAX2, LAX1 and LOX3 play important roles in callus proliferation and plant regeneration by CRISPR/Cas9 editing and overexpression assay. CONCLUSIONS: This study provides novel insights on the role of the regulatory network in cell fate transition and reprogramming during plant regeneration driven by somatic embryogenesis.

A near-complete genome assembly of Catharanthus roseus and insights into its vinblastine biosynthesis and high susceptibility to the Huanglongbing pathogen.

This study reports the assembly of a near-complete genome of Catharanthus roseus, consisting of 561.7 Mb scaffolded into 8 pseudochromosomes with a contig N50 of 24.7 Mb and a scaffold N50 of 71.1 Mb. The assembly enables the construction of a gene regulatory network of the vinblastine biosynthetic pathway and provides insights into the high susceptibility of C. roseus to the Huanglongbing pathogen.

Cotton 4-coumarate-CoA ligase 3 enhanced plant resistance to Verticillium dahliae by promoting jasmonic acid signaling-mediated vascular lignification and metabolic flux.

Lignins and their antimicrobial-related polymers cooperatively enhance plant resistance to pathogens. Several isoforms of 4-coumarate-coenzyme A ligases (4CLs) have been identified as indispensable enzymes involved in lignin and flavonoid biosynthetic pathways. However, their roles in plant-pathogen interaction are still poorly understood. This study uncovers the role of Gh4CL3 in cotton resistance to the vascular pathogen Verticillium dahliae. The cotton 4CL3-CRISPR/Cas9 mutant (CR4cl) exhibited high susceptibility to V. dahliae. This susceptibility was most probably due to the reduction in the total lignin content and the biosynthesis of several phenolic metabolites, e.g., rutin, catechin, scopoletin glucoside, and chlorogenic acid, along with jasmonic acid (JA) attenuation. These changes were coupled with a significant reduction in 4CL activity toward p-coumaric acid substrate, and it is likely that recombinant Gh4CL3 could specifically catalyze p-coumaric acid to form p-coumaroyl-coenzyme A. Thus, overexpression of Gh4CL3 (OE4CL) showed increasing 4CL activity that augmented phenolic precursors, cinnamic, p-coumaric, and sinapic acids, channeling into lignin and flavonoid biosyntheses and enhanced resistance to V. dahliae. Besides, Gh4CL3 overexpression activated JA signaling that instantly stimulated lignin deposition and metabolic flux in response to pathogen, which all established an efficient plant defense response system, and inhibited V. dahliae mycelium growth. Our results propose that Gh4CL3 acts as a positive regulator for cotton resistance against V. dahliae by promoting JA signaling-mediated enhanced cell wall rigidity and metabolic flux.

A comprehensive overview of cotton genomics, biotechnology and molecular biological studies.

Cotton is an irreplaceable economic crop currently domesticated in the human world for its extremely elongated fiber cells specialized in seed epidermis, which makes it of high research and application value. To date, numerous research on cotton has navigated various aspects, from multi-genome assembly, genome editing, mechanism of fiber development, metabolite biosynthesis, and analysis to genetic breeding. Genomic and 3D genomic studies reveal the origin of cotton species and the spatiotemporal asymmetric chromatin structure in fibers. Mature multiple genome editing systems, such as CRISPR/Cas9, Cas12 (Cpf1) and cytidine base editing (CBE), have been widely used in the study of candidate genes affecting fiber development. Based on this, the cotton fiber cell development network has been preliminarily drawn. Among them, the MYB-bHLH-WDR (MBW) transcription factor complex and IAA and BR signaling pathway regulate the initiation; various plant hormones, including ethylene, mediated regulatory network and membrane protein overlap fine-regulate elongation. Multistage transcription factors targeting CesA 4, 7, and 8 specifically dominate the whole process of secondary cell wall thickening. And fluorescently labeled cytoskeletal proteins can observe real-time dynamic changes in fiber development. Furthermore, research on the synthesis of cotton secondary metabolite gossypol, resistance to diseases and insect pests, plant architecture regulation, and seed oil utilization are all conducive to finding more high-quality breeding-related genes and subsequently facilitating the cultivation of better cotton varieties. This review summarizes the paramount research achievements in cotton molecular biology over the last few decades from the above aspects, thereby enabling us to conduct a status review on the current studies of cotton and provide strong theoretical support for the future direction.