Heterologous expression of taxane genes confers resistance to fall armyworm in Nicotiana benthamiana.

KEY MESSAGE: Taxadiene synthase, taxadiene-5α-hydroxylase, and taxane 13α-hydroxylase genes were introduced into Nicotiana benthamiana, and the improved resistance to lepidoptera pest fall armyworm was reported. Fall armyworm (FAW) is a serious agricultural pest. Genetic engineering techniques have been used to create pest-resistant plant varieties for reducing pest damage. Paclitaxel is a diterpenoid natural metabolite with antineoplastic effects in medicine. However, the effects of taxanes on the growth and development of lepidoptera pests, such as the FAW, are unknown. Here, selected paclitaxel precursor biosynthesis pathway genes, taxadiene synthase, taxane 5α-hydroxylase, and taxane 13α-hydroxylase, were engineered in the heterologous host Nicotiana benthamiana plants. Bioassay experiments showed that the transgenic N. benthamiana plants displayed improved resistance to FAW infestation, with degeneration of gut tissues and induced expression of apoptosis-related genes. Cytotoxicity experiment showed that the paclitaxel precursor, 10-deacetylbaccatin III, is cytotoxic to Sf9 cells, causing cell cycle arrest at the G2/M phase and disorder of the cytoskeleton. Metabolome analysis showed that heterologous expression of taxane genes in N. benthamiana affected the digestive system, steroid hormone and purine metabolism pathways of FAW larvae. In summary, this study provides a candidate approach for FAW control.

Genome-wide characterization and expression analysis of Erf gene family in cotton.

BACKGROUND: AP2/ERF transcription factors are important in a variety of biological activities, including plant growth, development, and responses to biotic and abiotic stressors. However, little study has been done on cotton's AP2/ERF genes, although cotton is an essential fibre crop. We were able to examine the tissue and expression patterns of AP2/ERF genes in cotton on a genome-wide basis because of the recently published whole genome sequence of cotton. Genome-wide analysis of ERF gene family within two diploid species (G. arboreum & G. raimondii) and two tetraploid species (G. barbadense, G. hirsutum) was performed. RESULTS: A total of 118, 120, 213, 220 genes containing the sequence of single AP2 domain were identified in G. arboreum, G. raimondii, G. barbadense and G. hirsutum respectively. The identified genes were unevenly distributed across 13/26 chromosomes of A and D genomes of cotton. Synteny and collinearity analysis revealed that segmental duplications may have played crucial roles in the expansion of the cotton ERF gene family, as well as tandem duplications played a minor role. Cis-acting elements of the promoter sites of Ghi-ERFs genes predict the involvement in multiple hormone responses and abiotic stresses. Transcriptome and qRT-PCR analysis revealed that Ghi-ERF-2D.6, Ghi-ERF-12D.13, Ghi-ERF-6D.1, Ghi-ERF-7A.6 and Ghi-ERF-11D.5 are candidate genes against salinity tolerance in upland cotton. CONCLUSION: Overwhelmingly, the present study paves the way to better understand the evolution of cotton ERF genes and lays a foundation for future investigation of ERF genes in improving salinity stress tolerance in cotton.

HDAC3 Knockdown Dysregulates Juvenile Hormone and Apoptosis-Related Genes in Helicoverpa armigera.

Insect development requires genes to be expressed in strict spatiotemporal order. The dynamic regulation of genes involved in insect development is partly orchestrated by the histone acetylation-deacetylation via histone acetyltransferases (HATs) and histone deacetylases (HDACs). Although histone deacetylase 3 (HDAC3) is required for mice during early embryonic development, its functions in Helicoverpa armigera (H. armigera) and its potential to be used as a target of insecticides remain unclear. We treated H. armigera with HDAC3 siRNA and RGFP966, a specific inhibitor, examining how the HDAC3 loss-of-function affects growth and development. HDAC3 siRNA and RGFP966 treatment increased mortality at each growth stage and altered metamorphosis, hampering pupation and causing abnormal wing development, reduced egg production, and reduced hatching rate. We believe that the misregulation of key hormone-related genes leads to abnormal pupa development in HDAC3 knockout insects. RNA-seq analysis identified 2788 differentially expressed genes (≥two-fold change; p ≤ 0.05) between siHDAC3- and siNC-treated larvae. Krüppel homolog 1 (Kr-h1), was differentially expressed in HDAC3 knockdown larvae. Pathway-enrichment analysis revealed the significant enrichment of genes involved in the Hippo, MAPK, and Wnt signaling pathways following HDAC3 knockdown. Histone H3K9 acetylation was increased in H. armigera after siHDAC3 treatment. In conclusion, HDAC3 knockdown dysregulated juvenile hormone (JH)-related and apoptosis-related genes in H. armigera. The results showed that the HDAC3 gene is a potential target for fighting H. armigera.

Genome-Wide Identification and Characterization of the PERK Gene Family in Gossypium hirsutum Reveals Gene Duplication and Functional Divergence.

Proline-rich extensin-like receptor kinases (PERKs) are an important class of receptor kinases in plants. Receptor kinases comprise large gene families in many plant species, including the 15 PERK genes in Arabidopsis. At present, there is no comprehensive published study of PERK genes in G. hirsutum. Our study identified 33 PERK genes in G. hirsutum. Phylogenetic analysis of conserved PERK protein sequences from 15 plant species grouped them into four well defined clades. The GhPERK gene family is an evolutionarily advanced gene family that lost its introns over time. Several cis-elements were identified in the promoter regions of the GhPERK genes that are important in regulating growth, development, light responses and the response to several stresses. In addition, we found evidence for gene loss or addition through segmental or whole genome duplication in cotton. Gene duplication and synteny analysis identified 149 orthologous/paralogous gene pairs. Ka/Ks values show that most GhPERK genes experienced strong purifying selection during the rapid evolution of the gene family. GhPERK genes showed high expression levels in leaves and during ovule development. Furthermore, the expression of GhPERK genes can be regulated by abiotic stresses and phytohormone treatments. Additionally, PERK genes could be involved in several molecular, biological and physiological processes that might be the result of functional divergence.

Cotton polyamine oxidase is required for spermine and camalexin signalling in the defence response to Verticillium dahliae.

Verticillium dahliae is a destructive, soil-borne fungal pathogen that causes vascular wilt disease in many economically important crops worldwide. A polyamine oxidase (PAO) gene was identified and cloned by screening suppression subtractive hybridisation and cDNA libraries of cotton genotypes tolerant to Verticillium wilt and was induced early and strongly by inoculation with V. dahliae and application of plant hormone. Recombinant cotton polyamine oxidase (GhPAO) was found to catalyse the conversion of spermine (Spm) to spermidine (Spd) in vitro. Constitutive expression of GhPAO in Arabidopsis thaliana produced improved resistance to V. dahliae and maintained putrescine, Spd and Spm at high levels. Hydrogen peroxide (H2 O2 ), salicylic acid and camalexin (a phytoalexin) levels were distinctly increased in GhPAO-overexpressing Arabidopsis plants during V. dahliae infection when compared with wild-type plants, and Spm and camalexin efficiently inhibited growth of V. dahliae in vitro. Spermine promoted the accumulation of camalexin by inducing the expression of mitogen-activated protein kinases and cytochrome P450 proteins in Arabidopsis and cotton plants. The three polyamines all showed higher accumulation in tolerant cotton cultivars than in susceptible cotton cultivars after inoculation with V. dahliae. GhPAO silencing in cotton significantly reduced the Spd level and increased the Spm level, leading to enhanced susceptibility to infection by V. dahliae, and the levels of H2 O2 and camalexin were distinctly lower in GhPAO-silenced cotton plants after V. dahliae infection. Together, these results suggest that GhPAO contributes to resistance of the plant against V. dahliae through the mediation of Spm and camalexin signalling.

Cotton S-adenosylmethionine decarboxylase-mediated spermine biosynthesis is required for salicylic acid- and leucine-correlated signaling in the defense response to Verticillium dahliae.

MAIN CONCLUSION: Cotton S-adenosylmethionine decarboxylase-, rather than spermine synthase-, mediated spermine biosynthesis is required for salicylic acid- and leucine-correlated signaling in the defense response to Verticillium dahliae. Spermine (Spm) signaling is correlated with plant resistance to the fungal pathogen Verticillium dahliae. We identified genes for key rate-limiting enzymes in the biosynthesis of Spm, namely S-adenosylmethionine decarboxylase (GhSAMDC) and Spm synthase (GhSPMS). These were found by screening suppression subtractive hybridization and cDNA libraries of cotton (Gossypium) species tolerant to Verticillium wilt. Both were induced early and strongly by inoculation with V. dahliae and application of plant hormones. Silencing of GhSPMS or GhSAMDC in cotton leaves led to a significant accumulation of upstream substrates and, ultimately, enhanced plant susceptibility to Verticillium infection. Exogenous supplementation of Spm to the silenced cotton plants improved resistance. When compared with the wild type (WT), constitutive expression of GhSAMDC in Arabidopsis thaliana was associated with greater Verticillium wilt resistance and higher accumulations of Spm, salicylic acid, and leucine during the infection period. By contrast, transgenic Arabidopsis plants that over-expressed GhSPMS were unexpectedly more susceptible than the WT to V. dahliae and they also had impaired levels of putrescine (Put) and salicylic acid (SA). The susceptibility exhibited in GhSPMS-overexpressing Arabidopsis plants was partially reversed by the exogenous supply of Put or SA. In addition, the responsiveness of those two transgenic Arabidopsis lines to V. dahliae was associated with an alteration in transcripts of genes involved in plant resistance to epidermal penetrations and amino acid signaling. Together, these results suggest that GhSAMDC-, rather than GhSPMS-, mediated spermine biosynthesis contributes to plant resistance against V. dahliae through SA- and leucine-correlated signaling.

Cotton ACAULIS5 is involved in stem elongation and the plant defense response to Verticillium dahliae through thermospermine alteration.

KEY MESSAGE: Overexpression of GhACL5 , an ACAULIS5 from cotton, in Arabidopsis increased plant height and T-Spm level. Silencing of GhACL5 in cotton exhibited a dwarf phenotype and reduced resistance to Verticillium dahliae. The Arabidopsis thaliana gene ACAULIS5 (ACL5), for which inactivation causes a defect in stem elongation, encodes thermospermine (T-Spm) synthase. However, limited information is available about improvement in plant height by the overexpression of ACL5 gene, and the biological functions of ACL5 genes in response to biotic stress. Here, this study reports that constitutive expression of the cotton ACL5 gene (GhACL5) in Arabidopsis thaliana significantly increased plant height and elevated the level of T-Spm. Silencing of that gene in cotton reduced the amount of T-Spm and led to a severe dwarf phenotype. Expression of GhACL5 was induced upon treatment with the fungal pathogen Verticillium dahliae and plant hormones salicylic acid, jasmonic acid, and ethylene in resistant cotton plants, but gene silencing in cotton enhanced their susceptibility to V. dahliae infection. Furthermore, T-Spm exposure effectively inhibited V. dahliae growth in vitro. In summary, GhACL5 expression is related to in planta levels of T-Spm and is involved in stem elongation and defense responses against V. dahliae.

iJAZ-based approach to engineer lepidopteran pest resistance in multiple crop species.

The fall armyworm (FAW) poses a significant threat to global crop production. Here we showed that overexpression of jasmonate ZIM-domain (JAZ) protein GhJAZ24 confers resistance to cotton bollworm and FAW, while also causing sterility in transgenic cotton by recruiting TOPLESS and histone deacetylase 6. We identified the NGR motif of GhJAZ24 that recognizes and binds the aminopeptidase N receptor, enabling GhJAZ24 to enter cells and disrupt histone deacetylase 3, leading to cell death. To overcome plant sterility associated with GhJAZ24 overexpression, we developed iJAZ (i, induced), an approach involving damage-induced expression and a switch from intracellular to extracellular localization of GhJAZ24. iJAZ transgenic cotton maintained fertility and showed insecticidal activity against cotton bollworm and FAW. In addition, iJAZ transgenic rice, maize and tobacco plants showed insecticidal activity against their lepidopteran pests, resulting in an iJAZ-based approach for generating alternative insecticidal proteins with distinctive mechanisms of action, thus holding immense potential for future crop engineering.

Arginine Decarboxylase Gene ADC2 Regulates Fiber Elongation in Cotton.

Cotton is an important agro-industrial crop providing raw material for the textile industry. Fiber length is the key factor that directly affects fiber quality. ADC, arginine decarboxylase, is the key rate-limiting enzyme in the polyamine synthesis pathway; whereas, there is no experimental evidence that ADC is involved in fiber development in cotton yet. Our transcriptome analysis of the fiber initiation material of Gossypium arboreum L. showed that the expression profile of GaADC2 was induced significantly. Here, GhADC2, the allele of GaADC2 in tetraploid upland cotton Gossypium hirsutum L., exhibited up-regulated expression pattern during fiber elongation in cotton. Levels of polyamine are correlated with fiber elongation; especially, the amount of putrescine regulated by ADC was increased. Scanning electron microscopy showed that the fiber length was increased with exogenous addition of an ADC substrate or product putrescine; whereas, the fiber density was decreased with exogenous addition of an ADC specific inhibitor. Next, genome-wide transcriptome profiling of fiber elongation with exogenous putrescine addition was performed to determine the molecular basis in Gossypium hirsutum. A total of 3163 differentially expressed genes were detected, which mainly participated in phenylpropanoid biosynthesis, fatty acid elongation, and sesquiterpenoid and triterpenoid biosynthesis pathways. Genes encoding transcription factors MYB109, WRKY1, and TCP14 were enriched. Therefore, these results suggested the ADC2 and putrescine involvement in the development and fiber elongation of G. hirsutum, and provides a basis for cotton fiber development research in future.

Heterologous expression of cry3Bb1 and cry3 genes for enhanced resistance against insect pests in cotton.

Transgenic technology played a crucial role in developing insect-resistant plants resulting in the reduced application of pesticides. This article reports the expression of two cry proteins (Cry3Bb1 and Cry3) in cotton for enhanced resistance against chewing insect pests. The aforementioned genes were synthetically developed and were cloned under appropriate regulatory sequences followed by transformation into Eagle-2 genotype (Gossypium hirsutum) of cotton through shoot apex-cut Agro-infiltration. The transgene integration was validated by polymerase chain reaction using primers flanking the aforementioned cry genes. Transgene expression was assessed by qRT-PCR using GADPH as a reference gene. The relative fold expression analyses revealed the highest expression of the transgene(s) in M1 plants, which is a 4.5-fold expression (Cry3 + Cry3Bb1) followed by M3 (fold expression, 3.0) (Cry3Bb1) and M2 (fold expression, 2.5) (Cry3) transformants of cotton. The confirmed transgenic plants were exposed to insect pests, pink bollworm (Pectinophora gossypiella), and army bollworm (Helicoverpa armigera). Bioassay results revealed that 60% mortality was observed against pink bollworm, and 75% mortality was observed against army bollworm in transgenic plants containing both Cry3Bb1 and Cry3 genes (M1 transgenic plants). In M2 transgenic plants containing only the Cry3Bb1 gene, the mortality was observed to be 40% in the pink bollworm population, whereas 45% mortality was observed in the army bollworm population. In the case of M3 transgenic plants containing single gene-Cry3, the mortality was 20% in the pink bollworm population, whereas 30% mortality was observed in the army bollworm population. Almost no mortality was observed in non-transgenic Eagle-2 control plants. Hence, the developed cotton transformants have improved resistance against chewing insect pests.

Unraveling Heat Tolerance in Upland Cotton (Gossypium hirsutum L.) Using Univariate and Multivariate Analysis.

The ever-changing global environment currently includes an increasing ambient temperature that can be a devastating stress for organisms. Plants, being sessile, are adversely affected by heat stress in their physiology, development, growth, and ultimately yield. Since little is known about the response of biochemical traits to high-temperature ambiance, we evaluated eight parental lines (five lines and three testers) and their 15 F1 hybrids under normal and high-temperature stress to assess the impact of these conditions over 2 consecutive years. The research was performed under a triplicate randomized complete block design including a split-plot arrangement. Data were recorded for agronomic, biochemical, and fiber quality traits. Mean values of agronomic traits were significantly reduced under heat stress conditions, while hydrogen peroxide, peroxidase, total soluble protein, superoxide dismutase, catalase (CAT), carotenoids, and fiber strength displayed higher mean values under heat stress conditions. Under both conditions, high genetic advance and high heritability were observed for seed cotton yield (SCY), CAT, micronaire value, plant height, and chlorophyll-a and b content, indicating that an additive type of gene action controls these traits under both the conditions. For more insights into variation, Pearson correlation analysis and principal component analysis (PCA) were performed. Significant positive associations were observed among agronomic, biochemical, and fiber quality-related traits. The multivariate analyses involving hierarchical clustering and PCA classified the 23 experimental genotypes into four groups under normal and high-temperature stress conditions. Under both conditions, the F1 hybrid genotype FB-SHAHEEN × JSQ WHITE GOLD followed by Ghuari-1, CCRI-24, Eagle-2 × FB-Falcon, Ghuari-1 × JSQ White Gold, and Eagle-2 exhibited better performance in response to high-temperature stress regarding the agronomic and fiber quality-related traits. The mentioned genotypes could be utilized in future cotton breeding programs to enhance heat tolerance and improve cotton yield and productivity through resistance to environmental stressors.

Target of Rapamycin Regulates Genome Methylation Reprogramming to Control Plant Growth in Arabidopsis.

DNA methylation is an indispensable epigenetic modification that dynamically regulates gene expression and genome stability during cell growth and development processes. The target of rapamycin (TOR) has emerged as a central regulator to regulate many fundamental cellular metabolic processes from protein synthesis to autophagy in all eukaryotic species. However, little is known about the functions of TOR in DNA methylation. In this study, the synergistic growth inhibition of Arabidopsis seedlings can be observed when DNA methylation inhibitor azacitidine was combined with TOR inhibitors. Global DNA methylation level was evaluated using whole-genome bisulfite sequencing (WGBS) under TOR inhibition. Hypomethylation level of whole genome DNA was observed in AZD-8055 (AZD), rapamycin (RAP) and AZD + RAP treated Arabidopsis seedlings. Based on functional annotation and KEGG pathway analysis of differentially methylated genes (DMGs), most of DMGs were enriched in carbon metabolism, biosynthesis of amino acids and other metabolic processes. Importantly, the suppression of TOR caused the change in DNA methylation of the genes associated with plant hormone signal transduction, indicating that TOR played an important role in modulating phytohormone signals in Arabidopsis. These observations are expected to shed light on the novel functions of TOR in DNA methylation and provide some new insights into how TOR regulates genome DNA methylation to control plant growth.

Identification of Histone H3 (HH3) Genes in Gossypium hirsutum Revealed Diverse Expression During Ovule Development and Stress Responses.

Histone acts as the core for nucleosomes and is a key protein component of chromatin. Among different histone variants, histone H3 (HH3) variants have been reported to play vital roles in plant development. However, biological information and evolutionary relationships of HH3 genes in cotton remain to be elucidated. The current study identified 34 HH3 genes in Gossypium hirsutum. Phylogenetic analysis classified HH3 genes of 19 plant species into eight distinct clades. Sequence logos analysis among Arabidopsis, rice, and G. hirsutum amino acid residues showed higher conservation in amino acids. Using collinearity analysis, we identified 81 orthologous/paralogous gene pairs among the four genomes (A, D, At, and Dt) of cotton. Further, orthologous/paralogous and the Ka/Ks ratio demonstrated that cotton HH3 genes experienced strong purifying selection pressure with restricted functional divergence resulting from segmental and whole genome duplication. Expression pattern analysis indicated that GhHH3 genes were preferentially expressed in cotton ovule tissues. Additionally, GhHH3 gene expression can be regulated by abiotic stresses (cold, heat, sodium chloride (NaCl), and polyethylene glycol (PEG)) and phytohormonal (brassinolide (BL), gibberellic acid (GA), indole-3-acetic acid (IAA), salicylic acid (SA), and methyl jasmonate (MeJA)) treatments, suggesting that GhHH3 genes might play roles in abiotic and hormone stress resistance. Taken together, this work provides important information to decipher complete molecular and physiological functions of HH3 genes in cotton.

Transcriptome profiling of Gossypium arboreum during fiber initiation and the genome-wide identification of trihelix transcription factors.

Cotton fiber initiation is the first step in fiber development, and it determines the yield. Here, genome-wide transcriptome profiling of Gossypium arboreum was performed to determine the molecular basis of cotton fiber initiation. A comparison of the transcriptomes of fiber-bearing ovules at -0.5, 0, 0.5, 1, 1.5, 2, 2.5 and 3 d post-anthesis detected 12,049 differentially expressed genes that mainly participated in ribosome, carbon metabolism and amino acid biosynthesis pathways. Genes encoding alcohol dehydrogenase 1 and hydroxycinnamoyl-CoA shikimate/quinate hydroxycinnamoyl transferase, involving in fatty acid degradation and flavonoid biosynthesis, were enriched. Furthermore, 1049 differentially expressed transcription factors were identified. Among these, 17 were trihelix family transcription factors, which play important roles in plant development and responses to biotic and abiotic stresses. In total, 52 full-length trihelix genes, named as GaGTs, were identified in G. arboreum and located in 12 of the 13 cotton chromosomes. Transcriptomic data and a quantitative real-time PCR analysis indicated that several GaGTs were significantly induced during fiber initiation in G. arboreum. Thus, the genome-wide comprehensive analysis of gene expression in G. arboreum fiber initiation will serve as a useful resource for unraveling the functions of specific genes. The phylogenetic relationships and expression analyses of the G. arboreum trihelix genes established a solid foundation for future comprehensive functional analyses of the GaGTs.

Extensive intraspecific gene order and gene structural variations in upland cotton cultivars.

Multiple cotton genomes (diploid and tetraploid) have been assembled. However, genomic variations between cultivars of allotetraploid upland cotton (Gossypium hirsutum L.), the most widely planted cotton species in the world, remain unexplored. Here, we use single-molecule long read and Hi-C sequencing technologies to assemble genomes of the two upland cotton cultivars TM-1 and zhongmiansuo24 (ZM24). Comparisons among TM-1 and ZM24 assemblies and the genomes of the diploid ancestors reveal a large amount of genetic variations. Among them, the top three longest structural variations are located on chromosome A08 of the tetraploid upland cotton, which account for ~30% total length of this chromosome. Haplotype analyses of the mapping population derived from these two cultivars and the germplasm panel show suppressed recombination rates in this region. This study provides additional genomic resources for the community, and the identified genetic variations, especially the reduced meiotic recombination on chromosome A08, will help future breeding.