RNA-Sequencing, Physiological and RNAi Analyses Provide Insights into the Response Mechanism of the ABC-Mediated Resistance to Verticillium dahliae Infection in Cotton

SR45a plays a key role in enhancing cotton resistance to Verticillium dahliae by alternative splicing of immunity genes

Alternative splicing (AS) of pre-mRNAs increases the diversity of transcriptome and proteome and plays fundamental roles in plant development and stress responses. However, the prevalent changes in AS events and the regulating mechanisms of plants in response to pathogens remain largely unknown. Here, we show that AS changes are an important mechanism conferring cotton immunity to Verticillium dahliae (Vd). GauSR45a, encoding a serine/arginine-rich RNA binding protein, was upregulated expression and underwent AS in response to Vd infection in Gossypium australe, a wild diploid cotton species highly resistant to Vd. Silencing GauSR45a substantially reduced the splicing ratio of Vd-induced immune-associated genes, including GauBAK1 (BRI1-associated kinase 1) and GauCERK1 (chitin elicitor receptor kinase 1). GauSR45a binds to the GAAGA motif that is commonly found in the pre-mRNA of genes essential for PTI, ETI, and defense. The binding between GauSR45a and the GAAGA motif in the pre-mRNA of BAK1 was enhanced by two splicing factors of GauU2AF35B and GauU1-70 K, thereby facilitating exon splicing; silencing either AtU2AF35B or AtU1-70 K decreased the resistance to Vd in transgenic GauSR45a Arabidopsis. Overexpressing the short splicing variant of BAK1GauBAK1.1 resulted in enhanced Verticillium wilt resistance rather than the long one GauBAK1.2. Vd-induced far more AS events were in G. barbadense (resistant tetraploid cotton) than those in G. hirsutum (susceptible tetraploid cotton) during Vd infection, indicating resistance divergence in immune responses at a genome-wide scale. We provided evidence showing a fundamental mechanism by which GauSR45a enhances cotton resistance to Vd through global regulation of AS of immunity genes.

An insight into the gene expression evolution in Gossypium species based on the leaf transcriptomes

BACKGROUND

Gene expression pattern is associated with biological phenotype and is widely used in exploring gene functions. Its evolution is also crucial in understanding species speciation and divergence. The genus Gossypium is a bona fide model for studying plant evolution and polyploidization. However, the evolution of gene expression during cotton species divergence has yet to be extensively discussed.

RESULTS

Based on the seedling leaf transcriptomes, this work analyzed the transcriptomic content and expression patterns across eight cotton species, including six diploids and two natural tetraploids. Our findings indicate that, while the biological function of these cotton transcriptomes remains largely conserved, there has been significant variation in transcriptomic content during species divergence. Furthermore, we conducted a comprehensive analysis of expression distances across cotton species. This analysis lends further support to the use of G. arboreum as a substitute for the A-genome donor of natural cotton polyploids. Moreover, our research highlights the evolution of stress-responsive pathways, including hormone signaling, fatty acid degradation, and flavonoid biosynthesis. These processes appear to have evolved under lower selection pressures, presumably reflecting their critical role in the adaptations of the studied cotton species to diverse environments.

CONCLUSIONS

In summary, this study provided insights into the gene expression variation within the genus Gossypium and identified essential genes/pathways whose expression evolution was closely associated with the evolution of cotton species. Furthermore, the method of characterizing genes and pathways under unexpected high or slow selection pressure can also serve as a new strategy for gene function exploration.

Decoding the guardians of cotton resilience: A comprehensive exploration of the βCA genes and its role in Verticillium dahliae resistance

Plant Carbonic anhydrases (Cas) have been shown to be stress-responsive enzymes that may play a role in adapting to adverse conditions. Cotton is a significant economic crop in China, with upland cotton (Gossypium hirsutum) being the most widely cultivated species. We conducted genome-wide identification of the βCA gene in six cotton species and preliminary analysis of the βCA gene in upland cotton. In total, 73 βCA genes from six cotton species were identified, with phylogenetic analysis dividing them into five subgroups. GHβCA proteins were predominantly localized in the chloroplast and cytoplasm. The genes exhibited conserved motifs, with motifs 1, 2, and 3 being prominent. GHβCA genes were unevenly distributed across chromosomes and were associated with stress-responsive cis-regulatory elements, including those responding to light, MeJA, salicylic acid, abscisic acid, cell cycle regulation, and defence/stress. Expression analysis indicated that GHβCA6, GHβCA7, GHβCA10, GHβCA15, and GHβCA16 were highly expressed under various abiotic stress conditions, whereas GHβCA3, GHβCA9, GHβCA10, and GHβCA18 had higher expression patterns under Verticillium dahliae infection at different time intervals. In Gossypium thurberi, GthβCA1, GthβCA2, and GthβCA4 showed elevated expression across stress conditions and tissues. Silencing GHβCA10 through VIGS increased Verticillium wilt severity and reduced lignin deposition compared to non-silenced plants. GHβCA10 is crucial for cotton's defense against Verticillium dahliae. Further research is needed to understand the underlying mechanisms and develop strategies to enhance resistance against Verticillium wilt.

Dynamic characteristics and functional analysis provide new insights into the role of GauERF105 for resistance against Verticillium dahliae in cotton

BACKGROUND

The cotton industry suffers significant yield losses annually due to Verticillium wilt, which is considered the most destructive disease affecting the crop. However, the precise mechanisms behind this disease in cotton remain largely unexplored.

METHODS

Our approach involved utilizing transcriptome data from G. australe which was exposed to Verticillium dahliae infection. From this data, we identified ethylene-responsive factors and further investigated their potential role in resistance through functional validations via Virus-induced gene silencing (VIGS) in cotton and overexpression in Arabidopsis.

RESULTS

A total of 23 ethylene response factors (ERFs) were identified and their expression was analyzed at different time intervals (24 h, 48 h, and 72 h post-inoculation). Among them, GauERF105 was selected based on qRT-PCR expression analysis for further investigation. To demonstrate the significance of GauERF105, VIGS was utilized, revealing that suppressing GauERF105 leads to more severe infections in cotton plants compared to the wild-type. Additionally, the silenced plants exhibited reduced lignin deposition in the stems compared to the WT plants, indicating that the silencing of GauERF105 also impacts lignin content. The overexpression of GauERF105 in Arabidopsis confirmed its pivotal role in conferring resistance against Verticillium dahliae infection. Our results suggest that WT possesses higher levels of the oxidative stress markers MDA and H2O2 as compared to the overexpressed lines. In contrast, the activities of the antioxidant enzymes SOD and POD were higher in the overexpressed lines compared to the WT. Furthermore, DAB and trypan staining of the overexpressed lines suggested a greater impact of the disease in the wild-type compared to the transgenic lines.

CONCLUSIONS

Our findings provide confirmation that GauERF105 is a crucial candidate in the defense mechanism of cotton against Verticillium dahliae invasion, and plays a pivotal role in this process. These results have the potential to facilitate the development of germplasm resistance in cotton.

Genome-Wide Analysis of Lipoxygenase (LOX) Genes in Angiosperms

Lipoxygenases (LOXs) are enzymes that catalyze the addition of an oxygen molecule to unsaturated fatty acids, thus forming hydroperoxides. In plants, these enzymes are encoded by a multigene family found in several organs with varying activity patterns, by which they are classified as LOX9 or LOX13. They are involved in several physiological functions, such as growth, fruit development, and plant defense. Despite several studies on genes of the LOX family in plants, most studies are restricted to a single species or a few closely related species. This study aimed to analyze the diversity, evolution, and expression of LOX genes in angiosperm species. We identified 247 LOX genes among 23 species of angiosperms and basal plants. Phylogenetic analyses identified clades supporting LOX13 and two main clades for LOX9: LOX9_A and LOX9_B. Eudicot species such as Tarenaya hassleriana, Capsella rubella, and Arabidopsis thaliana did not present LOX9_B genes; however, LOX9_B was present in all monocots used in this study. We identified that there were potential new subcellular localization patterns and conserved residues of oxidation for LOX9 and LOX13 yet unexplored. In summary, our study provides a basis for the further functional and evolutionary study of lipoxygenases in angiosperms.

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.

Genome-wide identification and functional analysis of class III peroxidases in Gossypium hirsutum

Class III peroxidase (PRX) genes play essential roles in various processes, such as auxin catabolism, removal of H₂O₂, crosslinking cell wall components, and response to biotic and abiotic stresses. In this study, we identified 166, 78 and 89 PRX genes from G. hirsutum, G. arboretum and G. raimondii, respectively. These PRX genes were classified into seven subfamilies based on phylogenetic tree analysis and the classification of PRX genes in Arabidopsis. Segmental duplication and purifying selection were the major factors driving the evolution of GhPRXs. GO and KEGG enrichment analysis revealed that GhPRX genes were mainly associated with responding to oxidative stresses, peroxidase activities and phenylpropanoid biosynthesis pathways. Transcriptome data analysis showed that GhPRX genes expression were significantly different in microspore development between the sterility line-JinA and the maintainer line MB177. We confirmed the up-regulation of GhPRX107 and down-regulation of GhPRX128 in the sterile line compared to its maintainer line using qRT-PCR, suggesting their roles in pollen fertility. In addition, silencing GhPRX107 in cotton showed a significant decrease of the reactive oxygen species (ROS) levels of microsporocyte stage anthers compared to control. Overexpressing GhPRX107 in Arabidopsis significantly increased the ROS levels of anthers compared to wild type. In conclusion, we identified GhPRX107 as a determinant of ROS levels in anther. This work sets a foundation for PRX studies in pollen development.

Inventory of ATP-binding cassette proteins in Lithospermum erythrorhizon as a model plant producing divergent secondary metabolites

ATP-binding cassette (ABC) proteins are the largest membrane transporter family in plants. In addition to transporting organic substances, these proteins function as ion channels and molecular switches. The development of multiple genes encoding ABC proteins has been associated with their various biological roles. Plants utilize many secondary metabolites to adapt to environmental stresses and to communicate with other organisms, with many ABC proteins thought to be involved in metabolite transport. Lithospermum erythrorhizon is regarded as a model plant for studying secondary metabolism, as cells in culture yielded high concentrations of meroterpenes and phenylpropanoids. Analysis of the genome and transcriptomes of L. erythrorhizon showed expression of genes encoding 118 ABC proteins, similar to other plant species. The number of expressed proteins in the half-size ABCA and full-size ABCB subfamilies was ca. 50% lower in L. erythrorhizon than in Arabidopsis, whereas there was no significant difference in the numbers of other expressed ABC proteins. Because many ABCG proteins are involved in the export of organic substances, members of this subfamily may play important roles in the transport of secondary metabolites that are secreted into apoplasts.

ATP-binding cassette transporters in nonmodel plants

Knowledge about plant ATP-binding cassette (ABC) proteins is of great value for sustainable agriculture, economic yield, and the generation of high-quality products, especially under unfavorable growth conditions. We have learned much about ABC proteins in model organisms, notably Arabidopsis thaliana; however, the importance of research dedicated to these transporters extends far beyond Arabidopsis biology. Recent progress in genomic and transcriptomic approaches for nonmodel and noncanonical model plants allows us to look at ABC transporters from a wider perspective and consider chemodiversity and functionally driven adaptation as distinctive mechanisms during their evolution. Here, by considering several representatives from agriculturally important families and recent progress in functional characterization of nonArabidopsis ABC proteins, we aim to bring attention to understanding the evolutionary background, distribution among lineages and possible mechanisms underlying the adaptation of this versatile transport system for plant needs. Increasing the knowledge of ABC proteins in nonmodel plants will facilitate breeding and development of new varieties based on, for example, genetic variations of endogenous genes and/or genome editing, representing an alternative to transgenic approaches.

Differences in the Abundance of Auxin Homeostasis Proteins Suggest Their Central Roles for In Vitro Tissue Differentiation in Coffea arabica

Coffea arabica is one of the most important crops worldwide. In vitro culture is an alternative for achieving Coffea regeneration, propagation, conservation, genetic improvement, and genome editing. The aim of this work was to identify proteins involved in auxin homeostasis by isobaric tandem mass tag (TMT) and the synchronous precursor selection (SPS)-based MS3 technology on the Orbitrap Fusion™ Tribrid mass spectrometer™ in three types of biological materials corresponding to C. arabica: plantlet leaves, calli, and suspension cultures. Proteins included in the β-oxidation of indole butyric acid and in the signaling, transport, and conjugation of indole-3-acetic acid were identified, such as the indole butyric response (IBR), the auxin binding protein (ABP), the ATP-binding cassette transporters (ABC), the Gretchen-Hagen 3 proteins (GH3), and the indole-3-acetic-leucine-resistant proteins (ILR). A more significant accumulation of proteins involved in auxin homeostasis was found in the suspension cultures vs. the plantlet, followed by callus vs. plantlet and suspension culture vs. callus, suggesting important roles of these proteins in the cell differentiation process.

Transcriptome Analysis of a Cotton Cultivar Provides Insights into the Differentially Expressed Genes Underlying Heightened Resistance to the Devastating Verticillium Wilt

Cotton is an important economic crop worldwide. Verticillium wilt (VW) caused by Verticillium dahliae (V. dahliae) is a serious disease in cotton, resulting in massive yield losses and decline of fiber quality. Breeding resistant cotton cultivars is an efficient but elaborate method to improve the resistance of cotton against V. dahliae infection. However, the functional mechanism of several excellent VW resistant cotton cultivars is poorly understood at present. In our current study, we carried out RNA-seq to discover the differentially expressed genes (DEGs) in the roots of susceptible cotton Gossypium hirsutum cultivar Junmian 1 (J1) and resistant cotton G.hirsutum cultivar Liaomian 38 (L38) upon Vd991 inoculation at two time points compared with the mock inoculated control plants. The potential function of DEGs uniquely expressed in J1 and L38 was also analyzed by GO enrichment and KEGG pathway associations. Most DEGs were assigned to resistance-related functions. In addition, resistance gene analogues (RGAs) were identified and analyzed for their role in the heightened resistance of the L38 cultivar against the devastating Vd991. In summary, we analyzed the regulatory network of genes in the resistant cotton cultivar L38 during V. dahliae infection, providing a novel and comprehensive insight into VW resistance in cotton.

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.

Dynamic characteristics and functional analysis provide new insights into long non-coding RNA responsive to Verticillium dahliae infection in Gossypium hirsutum

BACKGROUND

Verticillium wilt is a widespread and destructive disease, which causes serious loss of cotton yield and quality. Long non-coding RNA (lncRNA) is involved in many biological processes, such as plant disease resistance response, through a variety of regulatory mechanisms, but their possible roles in cotton against Verticillium dahliae infection remain largely unclear.

RESULTS

Here, we measured the transcriptome of resistant G. hirsutum following infection by V. dahliae and 4277 differentially expressed lncRNAs (delncRNAs) were identified. Localization and abundance analysis revealed that delncRNAs were biased distribution on chromosomes. We explored the dynamic characteristics of disease resistance related lncRNAs in chromosome distribution, induced expression profiles, biological function, and these lncRNAs were divided into three categories according to their induced expression profiles. For the delncRNAs, 687 cis-acting pairs and 14,600 trans-acting pairs of lncRNA-mRNA were identified, which indicated that trans-acting was the main way of Verticillium wilt resistance-associated lncRNAs regulating target mRNAs in cotton. Analyzing the regulation pattern of delncRNAs revealed that cis-acting and trans-acting lncRNAs had different ways to influence target genes. Gene Ontology (GO) enrichment analysis revealed that the regulatory function of delncRNAs participated significantly in stimulus response process, kinase activity and plasma membrane components. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis indicated that delncRNAs participated in some important disease resistance pathways, such as plant-pathogen interaction, alpha-linolenic acid metabolism and plant hormone signal transduction. Additionally, 21 delncRNAs and 10 target genes were identified as being involved in alpha-linolenic acid metabolism associated with the biosynthesis of jasmonic acid (JA). Subsequently, we found that GhlncLOX3 might regulate resistance to V. dahliae through modulating the expression of GhLOX3 implicated in JA biosynthesis. Further functional analysis showed that GhlncLOX3-silenced seedlings displayed a reduced resistance to V. dahliae, with down-regulated expression of GhLOX3 and decreased content of JA.

CONCLUSION

This study shows the dynamic characteristics of delncRNAs in multiaspect, and suggests that GhlncLOX3-GhLOX3-JA network participates in response to V. dahliae invasion. Our results provide novel insights for genetic improvement of Verticillium wilt resistance in cotton using lncRNAs.

Transcriptome Analysis and RNA Interference Reveal GhGDH2 Regulating Cotton Resistance to Verticillium Wilt by JA and SA Signaling Pathways

Verticillium wilt, caused by Verticillium dahliae, is one of the most damaging and widespread soil-borne cotton diseases. The molecular mechanisms underlying the cotton defense against V. dahliae remain largely elusive. Here, we compared the transcriptional differences between Upland cotton cultivars: one highly resistant (HR; Shidalukang 1) and one highly susceptible (HS; Junmian 1). This was done at multiple time points after V. dahliae inoculation, which identified 2010 and 1275 differentially expressed genes (DEGs) in HR and HS, respectively. Plant hormone signal transduction-related genes were enriched in HR, whereas genes related to lignin biosynthesis were enriched in both HR and HS. Weighted gene co-expression network analysis (WGCNA) using the 2868 non-redundant genes differentially expressed between the V. dahliae infected and uninfected samples in HR or HS identified 10 different gene network modules and 22 hub genes with a potential role in regulating cotton defense against V. dahliae infection. GhGDH2, encoding glutamate dehydrogenase (GDH), was selected for functional characterization. Suppressing the expression level of GhGDH2 by virus-induced gene silencing (VIGS) in HS led to inhibition of the salicylic acid (SA) biosynthesis/signaling pathways and activation of the jasmonic acid (JA) biosynthesis/signaling pathways, which resulted in an increase of 42.1% JA content and a reduction of 78.9% SA content in cotton roots, and consequently enhanced V. dahliae resistance. Our finding provides new insights on the molecular mechanisms of cotton resistance to V. dahliae infection and candidate genes for breeding V. dahliae resistance cotton cultivars by genetic modification.

Genome-Wide Identification and Expression Profile Analysis of the PHT1 Gene Family in Gossypium hirsutum and Its Two Close Relatives of Subgenome Donor Species

Phosphate transporter (PHT) is responsible for plant phosphorus (P) absorption and transport. PHT1 is a component of the high-affinity phosphate transporter system and plays pivotal roles in P absorption under P starvation conditions. However, in cotton, the number and identity of PHT1 genes that are crucial for P absorption from soil remain unclear. Here, genome-wide identification detected twelve PHT1 genes in Gossypium hirsutum and seven and eight PHT1 genes in two close relatives of the G. hirsutum genome-G. arboreum and G. raimondii, respectively. In addition, under low-phosphate treatment, the expressions of GaPHT1;3, GaPHT1;4, and GaPHT1;5 in roots were upregulated after 3 h of induction, and GhPHT1;3-At, GhPHT1;4-At, GhPHT1;5-At, GhPHT1;3-Dt, GhPHT1;4-Dt, and GhPHT1;5-Dt in the roots began to respond after 1 h of induction. Homologous pairs-GaPHT1;4 and GhPHT1;4-At; GaPHT1;5 and GhPHT1;5-At; GrPHT1;4 and GhPHT1;4-Dt, with GhPHT1;5-Dt and GhPHT1;5-At being syntenic-were all highly expressed in the roots under normal conditions. Among the genes highly expressed in the roots, GhPHT1;4-At, GhPHT1;5-At, GhPHT1;4-Dt and GhPHT1;5-Dt were continuously upregulated by P starvation. Therefore, it is concluded that these four genes might be key genes for P uptake in cotton roots. The results of this study provide insights into the mechanisms of P absorption and transport in cotton.

Dynamic infection of Verticillium dahliae in upland cotton

Verticillium wilt, an infection caused by the soilborne fungus Verticillium dahliae, is one of the most serious diseases in cotton. No effective control method against V. dahliae has been established, and the infection mechanism of V. dahliae in upland cotton remains unknown. GFP-tagged V. dahliae isolates with different pathogenic abilities were used to analyse the colonisation and infection of V. dahliae in the roots and leaves of different upland cotton cultivars, the relationships among infection processes, the immune responses and the resistance ability of different cultivars against V. dahliae. Here, we report a new infection model for V. dahliae in upland cotton plants. V. dahliae can colonise and infect any organ of upland cotton plants and then spread to the entire plant from the infected organ through the surface and interior of the organ. Vascular tissue was found to not be the sole transmission route of V. dahliae in cotton plants. In addition, the rate of infection of a V. dahliae isolate with strong pathogenicity was notably faster than that of an isolate with weak pathogenicity. The resistance of upland cotton to Verticillium wilt was related to the degree of the immune response induced in plants infected with V. dahliae. These results provide a theoretical basis for studying the mechanism underlying the interaction between V. dahliae and upland cotton. These results provide a theoretical basis for studying the mechanism underlying the interaction between V. dahliae and upland cotton.

Comparative transcriptome, physiological and biochemical analyses reveal response mechanism mediated by CBF4 and ICE2 in enhancing cold stress tolerance in Gossypium thurberi

Low temperature is one of the key environmental stresses that impair plant growth and significantly restricts the productivity and spatial distribution of crop plants. Gossypium thurberi, a wild diploid cotton species, has adapted to a wide range of temperatures and exhibits a better tolerance to chilling stress. Here, we compared phenotypes and physiochemical changes in G. thurberi under cold stress and found this species indeed showed better cold tolerance. Therefore, to understand the molecular mechanisms of the cold tolerance in G. thurberi, we compared transcription changes in leaves of G. thurberi under cold stress by high-throughput transcriptome sequencing. In total, 35 617 unigenes were identified in the whole-genome transcription profile, and 4226 differentially expressed genes (DEGs) were discovered in the leaves upon cold treatment. Gene Ontology (GO) classification analyses showed that the majority of DEGs belonged to categories of signal transduction, transcription factors (TFs) and carbohydrate transport and metabolism. The expression of several cold-responsive genes such as ICE1, CBF4, RAP2-7 and abscisic acid (ABA) biosynthesis genes involved in different signalling pathways were induced after G. thurberi seedlings were exposed to cold stress. Furthermore, cold sensitivity was increased in CBF4 and ICE2 virus-induced gene silencing (VIGS) plants, and high level of malondialdehyde (MDA) showed that the CBF4 and ICE2 silenced plants were under oxidative stress compared to their wild types, which relatively had higher levels of antioxidant enzyme activity, as evident by high levels of proline and superoxide dismutase (SOD) content. In conclusion, our findings reveal a new regulatory network of cold stress response in G. thurberi and broaden our understanding of the cold tolerance mechanism in cotton, which might accelerate functional genomics studies and genetic improvement for cold stress tolerance in cultivated cotton.

Transcriptome sequencing of Verticillium dahliae from a cotton farm reveals positive correlation between virulence and tolerance of sugar-induced hyperosmosis

Verticillium dahliae causes disease symptoms in its host plants; however, due to its rapid variability, V. dahliae is difficult to control. To analyze the reason for this pathogenic differentiation, 22 V. dahliae strains with different virulence were isolated from a cotton farm. The genetic diversity of cotton varieties make cotton cultivars have different Verticillium wilt resistance, so the Xinluzao 7 (susceptible to V. dahliae), Zhongmian 35 (tolerant), and Xinluzao 33 (resistant) were used to investigate the pathogenicity of the strains in a green house. Vegetative compatibility groups (VCGs) assays, Internal Transcribed Spacer (ITS) PCR, and pathogenicity analysis showed that SHZ-4, SHZ-5, and SHZ-9 had close kinship and significantly different pathogenicity. Transcriptome sequencing of the three strains identified 19 of 146 unigenes in SHZ-4_vs_ SHZ-5, SHZ-5_vs_ SHZ-9, and SHZ-4_vs_ SHZ-9. In these unigenes, three proteinase and four polysaccharide degrading hydrolases were found to be associated with the pathogenicity. However, due to a number of differentially expressed genes in the transport, these unigenes not only played a role in nutrition absorption but might also contribute to the resistance of sugar-induced hyperosmosis. Moreover, the tolerance ability was positively related to the pathogenicity of V. dahliae. This resistance to sugar-induced hyperosmosis might help V. dahliae to access the nutrition of the host. The pathogenicity of V. dahliae correlated with the resistance of sugar-induced-hyperosmosis, which provides clues for the cultivation of V. dahliae resistant varieties.

Polyamine Oxidases Play Various Roles in Plant Development and Abiotic Stress Tolerance

Polyamines not only play roles in plant growth and development, but also adapt to environmental stresses. Polyamines can be oxidized by copper-containing diamine oxidases (CuAOs) and flavin-containing polyamine oxidases (PAOs). Two types of PAOs exist in the plant kingdom; one type catalyzes the back conversion (BC-type) pathway and the other catalyzes the terminal catabolism (TC-type) pathway. The catabolic features and biological functions of plant PAOs have been investigated in various plants in the past years. In this review, we focus on the advance of PAO studies in rice, Arabidopsis, and tomato, and other plant species.