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.

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.

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.

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.

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.

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.