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.

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.

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.

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.

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.

Lateral transfers lead to the birth of momilactone biosynthetic gene clusters in grass.

Momilactone A, an important plant labdane-related diterpenoid, functions as a phytoalexin against pathogens and an allelochemical against neighboring plants. The genes involved in the biosynthesis of momilactone A are found in clusters, i.e., momilactone A biosynthetic gene clusters (MABGCs), in the rice and barnyardgrass genomes. In addition, we know little about the origin and evolution of MABGCs. Here, we integrated results from comprehensive phylogeny and comparative genomic analyses of the core genes of MABGC-like clusters and MABGCs in 40 monocot plant genomes, providing convincing evidence for the birth and evolution of MABGCs in grass species. The MABGCs found in the PACMAD clade of the core grass lineage (including Panicoideae and Chloridoideae) originated from a MABGC-like cluster in Triticeae (BOP clade) via lateral gene transfer (LGT) and followed by recruitment of MAS1/2 and CYP76L1 genes. The MABGCs in Oryzoideae originated from PACMAD through another LGT event and lost CYP76L1 afterwards. The Oryza MABGC and another Oryza diterpenoid cluster c2BGC are two distinct clusters, with the latter originating from gene duplication and relocation within Oryzoideae. Further comparison of the expression patterns of the MABGC genes between rice and barnyardgrass in response to pathogen infection and allelopathy provides novel insights into the functional innovation of MABGCs in plants. Our results demonstrate LGT-mediated origination of MABGCs in grass and shed lights into the evolutionary innovation and optimization of plant biosynthetic pathways.

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.

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.

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.

Characteristics of plant circular RNAs.

Circular RNA (circRNA) is a kind of covalently closed single-stranded RNA molecules that have been proved to play important roles in transcriptional regulation of genes in diverse species. With the rapid development of bioinformatics tools, a huge number (95143) of circRNAs have been identified from different plant species, providing an opportunity for uncovering the overall characteristics of plant circRNAs. Here, based on publicly available circRNAs, we comprehensively analyzed characteristics of plant circRNAs with the help of various bioinformatics tools as well as in-house scripts and workflows, including the percentage of coding genes generating circRNAs, the frequency of alternative splicing events of circRNAs, the non-canonical splicing signals of circRNAs and the networks involving circRNAs, miRNAs and mRNAs. All this information has been integrated into an upgraded online database, PlantcircBase 3.0 (http://ibi.zju.edu.cn/plantcircbase/). In this database, we provided browse, search and visualization tools as well as a web-based blast tool, BLASTcirc, for prediction of circRNAs from query sequences based on searching against plant genomes and transcriptomes.

Recent origination of circular RNAs in plants.

Circular RNA (circRNA) is a kind of new regulatory RNA with diverse biological functions. Numerous circRNAs have been identified in many plant species; however, evolution of plant circRNAs remains largely unknown. In this study, we assembled full-length sequences of 6519 rice (Oryza sativa) circRNAs and analyzed their conservation in another 46 plant species based on comparison of sequences and expression patterns. We found that, at the genomic level, 8.7% of the 6519 circRNAs were conserved in dicotyledonous plants and 49.1% in Oryza genus. Meanwhile, 57.8% of parental protein-coding genes of the rice circRNAs originated recently after divergence of monocotyledonous plants, implying recent origin of the majority of rice circRNAs, a conclusion further supported by the results based on analysis of 4663 full-length circRNAs in Arabidopsis thaliana. Accordingly, we proposed three models to address the origination of different types of circRNAs. Taken together, the results obtained in this study provide new insights for the evolutionary dynamics of plant circRNAs and candidate circRNAs for further functional exploration.

Genomic prediction of cotton fibre quality and yield traits using Bayesian regression methods.

Genomic selection or genomic prediction (GP) has increasingly become an important molecular breeding technology for crop improvement. GP aims to utilise genome-wide marker data to predict genomic breeding value for traits of economic importance. Though GP studies have been widely conducted in various crop species such as wheat and maize, its application in cotton, an essential renewable textile fibre crop, is still significantly underdeveloped. We aim to develop a new GP-based breeding system that can improve the efficiency of our cotton breeding program. This article presents a GP study on cotton fibre quality and yield traits using 1385 breeding lines from the Commonwealth Scientific and Industrial Research Organisation (CSIRO, Australia) cotton breeding program which were genotyped using a high-density SNP chip that generated 12,296 informative SNPs. The aim of this study was twofold: (1) to identify the models and data sources (i.e. genomic and pedigree) that produce the highest prediction accuracies; and (2) to assess the effectiveness of GP as a selection tool in the CSIRO cotton breeding program. The prediction analyses were conducted under various scenarios using different Bayesian predictive models. Results highlighted that the model combining genomic and pedigree information resulted in the best cross validated prediction accuracies: 0.76 for fibre length, 0.65 for fibre strength, and 0.64 for lint yield. Overall, this work represents the largest scale genomic selection studies based on cotton breeding trial data. Prediction accuracies reported in our study indicate the potential of GP as a breeding tool for cotton. The study highlighted the importance of incorporating pedigree and environmental factors in GP models to optimise the prediction performance.

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.

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.

Transcriptome Profiling Provides New Insights into the Molecular Mechanism Underlying the Sensitivity of Cotton Varieties to Mepiquat Chloride.

Mepiquat chloride (MC) is a plant growth regulator widely used in cotton production to control vegetative overgrowth of cotton plants to achieve ideal plant architecture required for high yielding. Cotton varieties respond differently to MC application, but there is little information about the molecular mechanisms underlying the varietal difference. In this study, comparative transcriptome analysis was conducted by using two Upland cotton varieties with different sensitivity (XLZ74, insensitive; SD1068, sensitive) to MC treatment, aiming to understand the molecular mechanisms responsible for varietal difference of MC sensitivity. RNA-seq data were generated from the two varieties treated with MC or water at three time points, 1, 3 and 6 days post-spray (dps). Genes differentially expressed between the MC and mock treatments of XLZ74 (6252) and SD1068 (6163) were subjected to Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses to compare the enriched GO terms and KEGG pathways between the two varieties. Signal transduction of phytohormones, biosynthesis of gibberellins (GAs) and brassinosteroids (BRs) and profiles of transcription factors (TFs) seemed to be differentially affected by MC in the two varieties. The transcriptomic results were further consolidated with the content changes of phytohormones in young stem. Several GA catabolic genes, GA2ox, were highly induced by MC in both varieties especially in SD1068, consistent with a more significant decrease in GA4 in SD1068. Several AUX/IAA and SAUR genes and CKX genes were induced by MC in both varieties, but with a more profound effect observed in SD1068 that showed a significant reduction in indole-3-acetic acid (IAA) and a significant increase in cytokinin (CTK) at 6 days post-spray (dps). BR biosynthesis-related genes were downregulated in SD1068, but not in XLZ74. Additionally, more downregulated TFs were observed in MC-treated SD1068 than in MC-treated XLZ74, and the two varieties had very different profiles of genes involved in starch and sucrose metabolism, with those of SD1068 and XLZ74 being downregulated and upregulated by MC treatment, respectively. Together, these results indicate that although the same or similar biological pathways are affected by MC treatment in cotton varieties showing different MC sensitivity, the extent of effect is variable, leading to their different phenotypic outcomes. How the quantitative effect of MC on the biological processes associated with growth retardation is regulated is still an open question.

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.

Expansion of MIR482/2118 by a class-II transposable element in cotton.

Some plant microRNA (miRNA) families contain multiple members generating identical or highly similar mature miRNA variants. Mechanisms underlying the expansion of miRNA families remain elusive, although tandem and/or segmental duplications have been proposed. In this study of two tetraploid cottons, Gossypium hirsutum and Gossypium barbadense, and their extant diploid progenitors, Gossypium arboreum and Gossypium raimondii, we investigated the gain and loss of members of the miR482/2118 superfamily, which modulates the expression of nucleotide-binding site leucine-rich repeat (NBS-LRR) disease resistance genes. We found significant expansion of MIR482/2118d in G. barbadense, G. hirsutum and G. raimondii, but not in G. arboreum. Several newly expanded MIR482/2118d loci have mutated to produce different miR482/2118 variants with altered target-gene specificity. Based on detailed analysis of sequences flanking these MIR482/2118 loci, we found that this expansion of MIR482/2118d and its derivatives resulted from an initial capture of an MIR482/2118d by a class-II DNA transposable element (TE) in G. raimondii prior to the tetraploidization event, followed by transposition to new genomic locations in G. barbadense, G. hirsutum and G. raimondii. The 'GosTE' involved in the capture and proliferation of MIR482/2118d and its derivatives belongs to the PIF/Harbinger superfamily, generating a 3-bp target site duplication upon insertion at new locations. All orthologous MIR482/2118 loci in the two diploids were retained in the two tetraploids, but mutation(s) in miR482/2118 were observed across all four species as well as in different cultivars of both G. barbadense and G. hirsutum, suggesting a dynamic co-evolution of miR482/2118 and its NBS-LRR targets. Our results provide fresh insights into the mechanisms contributing to MIRNA proliferation and enrich our knowledge on TEs.

The Cotton Lignin Biosynthetic Gene Gh4CL30 Regulates Lignification and Phenolic Content and Contributes to Verticillium Wilt Resistance.

Verticillium wilt is a vascular disease causing tremendous damage to cotton production worldwide. However, our knowledge of the mechanisms of cotton resistance or susceptibility to this disease is very limited. In this study, we compared the defense transcriptomes of cotton (Gossypium hirsutum) cultivars Shidalukang 1 (Verticillium dahliae resistant, HR) and Junmian 1 (V. dahliae susceptible, HS) before and after V. dahliae infection, identified hub genes of the network associated with responses to V. dahliae infection, and functionally characterized one of the hub genes involved in biosynthesis of lignin and phenolics. We identified 6,831 differentially expressed genes (DEGs) between the basal transcriptomes of HR and HS; 3,685 and 3,239 of these DEGs were induced in HR and HS, respectively, at different time points after V. dahliae infection. KEGG pathway analysis indicated that DEGs were enriched for genes involved in lignin biosynthesis. In all, 23 hub genes were identified based on a weighted gene coexpression network analysis of the 6,831 DEGs and their expression profiles at different time points after V. dahliae infection. Knockdown of Gh4CL30, one of the hub genes related to the lignin biosynthesis pathway, by virus-induced gene silencing, led to a decreased content of flavonoids, lignin, and S monomer but an increased content of G monomer, G/S lignin monomer, caffeic acid, and ferulic acid, and enhanced cotton resistance to V. dahliae. These results suggest that Gh4CL30 is a key gene modulating the outputs of different branches of the lignin biosynthesis pathway, and provide new insights into cotton resistance to V. dahliae.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Effects of Sample Size on Plant Single-Cell RNA Profiling.

Single-cell RNA (scRNA) profiling or scRNA-sequencing (scRNA-seq) makes it possible to parallelly investigate diverse molecular features of multiple types of cells in a given plant tissue and discover cell developmental processes. In this study, we evaluated the effects of sample size (i.e., cell number) on the outcome of single-cell transcriptome analysis by sampling different numbers of cells from a pool of ~57,000 Arabidopsis thaliana root cells integrated from five published studies. Our results indicated that the most significant principal components could be achieved when 20,000-30,000 cells were sampled, a relatively high reliability of cell clustering could be achieved by using ~20,000 cells with little further improvement by using more cells, 96% of the differentially expressed genes could be successfully identified with no more than 20,000 cells, and a relatively stable pseudotime could be estimated in the subsample with 5000 cells. Finally, our results provide a general guide for optimizing sample size to be used in plant scRNA-seq studies.

Explore the gene network regulating the composition of fatty acids in cottonseed.

BACKGROUND: Cottonseed is one of the major sources of vegetable oil. Analysis of the dynamic changes of fatty acid components and the genes regulating the composition of fatty acids of cottonseed oil is of great significance for understanding the biological processes underlying biosynthesis of fatty acids and for genetic improving the oil nutritional qualities. RESULTS: In this study, we investigated the dynamic relationship of 13 fatty acid components at 12 developmental time points of cottonseed (Gossypium hirsutum L.) and generated cottonseed transcriptome of the 12 time points. At 5-15 day post anthesis (DPA), the contents of polyunsaturated linolenic acid (C18:3n-3) and saturated stearic acid (C18:0) were higher, while linoleic acid (C18:2n-6) was mainly synthesized after 15 DPA. Using 5 DPA as a reference, 15,647 non-redundant differentially expressed genes were identified in 10-60 DPA cottonseed. Co-expression gene network analysis identified six modules containing 3275 genes significantly associated with middle-late seed developmental stages and enriched with genes related to the linoleic acid metabolic pathway and α-linolenic acid metabolism. Genes (Gh_D03G0588 and Gh_A02G1788) encoding stearoyl-ACP desaturase were identified as hub genes and significantly up-regulated at 25 DPA. They seemed to play a decisive role in determining the ratio of saturated fatty acids to unsaturated fatty acids. FAD2 genes (Gh_A13G1850 and Gh_D13G2238) were highly expressed at 25-50 DPA, eventually leading to the high content of C18:2n-6 in cottonseed. The content of C18:3n-3 was significantly decreased from 5 DPA (7.44%) to 25 DPA (0.11%) and correlated with the expression characteristics of Gh_A09G0848 and Gh_D09G0870. CONCLUSIONS: These results contribute to our understanding on the relationship between the accumulation pattern of fatty acid components and the expression characteristics of key genes involved in fatty acid biosynthesis during the entire period of cottonseed development.