Global gene expression profile and functional analysis reveal the conservation of reproduction-associated gene networks in Gossypium hirsutum.

KEY MESSAGE: Lastly, the bZIP gene family encompasses genes that have been reported to play a role in flower development, such as bZIP14 (FD). Notably, bZIP14 is essential for Flowering Locus T (FT) initiation of floral development in Arabidopsis (Abe et al. 2005). Cotton (Gossypium hirsutum L.) is the world's most extensively cultivated fiber crop. However, its reproductive development is poorly characterized at the molecular level. Thus, this study presents a detailed transcriptomic analysis of G. hirsutum at three different reproductive stages. We provide evidence that more than 64,000 genes are active in G. hirsutum during flower development, among which 94.33% have been assigned to functional terms and specific pathways. Gene set enrichment analysis (GSEA) revealed that the biological process categories of floral organ development, pollen exine formation, and stamen development were enriched among the genes expressed during the floral development of G. hirsutum. Furthermore, we identified putative Arabidopsis homologs involved in the G. hirsutum gene regulatory network (GRN) of pollen and flower development, including transcription factors such as WUSCHEL (WUS), INNER NO OUTER (INO), AGAMOUS-LIKE 66 (AGL66), SPOROCYTELESS/NOZZLE (SPL/NZZ), DYSFUNCTIONAL TAPETUM 1 (DYT1), ABORTED MICROSPORES (AMS), and ASH1-RELATED 3 (ASHR3), which are known crucial genes for plant reproductive success. The cotton MADS-box protein-protein interaction pattern resembles the previously described patterns for AGAMOUS (AG), SEEDSTICK (STK), SHATTERPROOF (SHP), and SEPALLATA3 (SEP3) homolog proteins from Arabidopsis. In addition to serving as a resource for comparative flower development studies, this work highlights the changes in gene expression profiles and molecular networks underlying stages that are valuable for cotton breeding improvement.

In vivo and in silico comparison analyses of Cry toxin activities toward the sugarcane giant borer.

The sugarcane giant borer, Telchin licus licus, is an insect pest that causes significant losses in sugarcane crops and in the sugar-alcohol sector. Chemical and manual control methods are not effective. As an alternative, in the current study, we have screened Bacillus thuringiensis (Bt) Cry toxins with high toxicity against this insect. Bioassays were conducted to determine the activity of four Cry toxins (Cry1A (a, b, and c) and Cry2Aa) against neonate T. licus licus larvae. Notably, the Cry1A family toxins had the lowest LC50 values, in which Cry1Ac presented 2.1-fold higher activity than Cry1Aa, 1.7-fold larger than Cry1Ab, and 9.7-fold larger than Cry2Aa toxins. In silico analyses were performed as a perspective to understand putative interactions between T. licus licus receptors and Cry1A toxins. The molecular dynamics and docking analyses for three putative aminopeptidase N (APN) receptors (TlAPN1, TlAPN3, and TlAPN4) revealed evidence for the amino acids that may be involved in the toxin-receptor interactions. Notably, the properties of Cry1Ac point to an interaction site that increases the toxin's affinity for the receptor and likely potentiate toxicity. The interacting amino acid residues predicted for Cry1Ac in this work are probably those shared by the other Cry1A toxins for the same region of APNs. Thus, the presented data extend the existing knowledge of the effects of Cry toxins on T. licus licus and should be considered in further development of transgenic sugarcane plants resistant to this major occurring insect pest in sugarcane fields.

Functional characterization of the pUceS8.3 promoter and its potential use for ectopic gene overexpression.

MAIN CONCLUSION: The pUceS8.3 is a constitutive gene promoter with potential for ectopic and strong genes overexpression or active biomolecules in plant tissues attacked by pests, including nematode-induced giant cells or galls. Soybean (Glycine max) is one of the most important agricultural commodities worldwide and a major protein and oil source. Herein, we identified the soybean ubiquitin-conjugating (E2) enzyme gene (GmUBC4; Glyma.18G216000), which is significantly upregulated in response to Anticarsia gemmatalis attack and Meloidogyne incognita-induced galls during plant parasitism by plant nematode. The GmUBC4 promoter sequence and its different modules were functionally characterized in silico and in planta using transgenic Arabidopsis thaliana and G. max lines. Its full-length transcriptional regulatory region (promoter and 5´-UTR sequences, named pUceS8.3 promoter) was able to drive higher levels of uidA (ß-glucuronidase) gene expression in different tissues of transgenic A. thaliana lines compared to its three shortened modules and the p35SdAMV promoter. Notably, higher ß-glucuronidase (GUS) enzymatic activity was shown in M. incognita-induced giant cells when the full pUceS8.3 promoter drove the expression of this reporter gene. Furthermore, nematode-specific dsRNA molecules were successfully overexpressed under the control of the pUceS8.3 promoter in transgenic soybean lines. The RNAi gene construct used here was designed to post-transcriptionally downregulate the previously characterized pre-mRNA splicing factor genes from Heterodera glycines and M. incognita. A total of six transgenic soybean lines containing RNAi gene construct were selected for molecular characterization after infection with M. incognita pre-parasitic second-stage (ppJ2) nematodes. A strong reduction in the egg number produced by M. incognita after parasitism was observed in those transgenic soybean lines, ranging from 71 to 92% compared to wild-type control plants. The present data demonstrated that pUceS8.3 is a gene promoter capable of effectively driving dsRNA overexpression in nematode-induced giant cells of transgenic soybean lines and can be successfully applied as an important biotechnological asset to generate transgenic crops with improved resistance to root-knot nematodes as well as other pests.

Discovery and functional characterization of novel cotton promoters with potential application to pest control.

KEY MESSAGE: pGhERF105 and pGhNc-HARBI1 promoters are highly responsive to CBW infestation and exhibit strong activity in vegetative and reproductive tissues, increasing their potential application in GM crop plants for pest control. The main challenge to cotton (Gossypium hirsutum) crop productivity is the constant attack of several pests, including the cotton boll weevil (CBW, Anthonomus grandis), which uses cotton floral buds for feeding and egg-laying. The endophytic nature of the early developmental stages of CBW makes conventional pesticide-based control poorly efficient. Most biotechnological assets used for pest control are based on Bacillus thurigiensis insecticidal Cry toxins or the silencing of insect-pest essential genes using RNA-interference technology. However, suitable plant promoter sequences are required to efficiently drive insecticidal molecules to the target plant tissue. This study selected the Ethylene Responsive Factor 105 (GhERF105) and Harbinger transposase-derived nuclease (GhNc-HARBI1) genes based on available transcriptome-wide data from cotton plants infested by CBW larvae. The GhERF105 and GhNc-HARBI1 genes showed induction kinetics from 2 to 96 h under CBW's infestation in cotton floral buds, uncovering the potential application of their promoters. Therefore, the promoter regions (1,500 base pairs) were assessed and characterized using Arabidopsis thaliana transgenic plants. The pGhERF105 and pGhNc-HARBI1 promoters showed strong activity in plant vegetative (leaves and roots) and reproductive (flowers and fruits) tissues, encompassing higher GUS transcriptional activity than the viral-constitutive Cauliflower Mosaic Virus 35S promoter (pCaMV35S). Notably, pGhERF105 and pGhNc-HARBI1 promoters demonstrated more efficiency in driving reporter genes in flowers than other previously characterized cotton flower-specific promoters. Overall, the present study provides a new set of cotton promoters suitable for biotechnological application in cotton plants for pest resistance.

Characterization of floral morphoanatomy and identification of marker genes preferentially expressed during specific stages of cotton flower development.

MAIN CONCLUSION: Characterization of anther and ovule developmental programs and expression analyses of stage-specific floral marker genes in Gossypium hirsutum allowed to build a comprehensive portrait of cotton flower development before fiber initiation. Gossypium hirsutum is the most important cotton species that is cultivated worldwide. Although cotton reproductive development is important for fiber production, since fiber is formed on the epidermis of mature ovules, cotton floral development remains poorly understood. Therefore, this work aims to characterize the cotton floral morphoanatomy by performing a detailed description of anther and ovule developmental programs and identifying stage-specific floral marker genes in G. hirsutum. Using light microscopy and scanning electron microscopy, we analyzed anther and ovule development during 11 stages of flower development. To better characterize the ovule development in cotton, we performed histochemical analyses to evaluate the accumulation of phenolic compounds, pectin, and sugar in ovule tissues. After identification of major hallmarks of floral development, three key stages were established in G. hirsutum floral development: in stage 1 (early-EF), sepal, petal, and stamen primordia were observed; in stage 2 (intermediate-IF), primordial ovules and anthers are present, and the differentiating archesporial cells were observed, marking the beginning of microsporogenesis; and in stage 6 (late-LF), flower buds presented initial anther tapetum degeneration and microspore were released from the tetrad, and nucellus and both inner and outer integuments are developing. We used transcriptome data of cotton EF, IF and LF stages to identify floral marker genes and evaluated their expression by real-time quantitative PCR (qPCR). Twelve marker genes were preferentially expressed in a stage-specific manner, including the putative homologs for AtLEAFY, AtAPETALA 3, AtAGAMOUS-LIKE 19 and AtMALE STERILITY 1, which are crucial for several aspects of reproductive development, such as flower organogenesis and anther and petal development. We also evaluated the expression profile of B-class MADS-box genes in G. hirsutum floral transcriptome (EF, IF, and LF). In addition, we performed a comparative analysis of developmental programs between Arabidopsis thaliana and G. hirsutum that considered major morphoanatomical and molecular processes of flower, anther, and ovule development. Our findings provide the first detailed analysis of cotton flower development.

Functional characterization of AGAMOUS-subfamily members from cotton during reproductive development and in response to plant hormones.

KEY MESSAGE: Expression analysis of the AG -subfamily members from G. hirsutum during flower and fruit development. Reproductive development in cotton, including the fruit and fiber formation, is a complex process; it involves the coordinated action of gene expression regulators, and it is highly influenced by plant hormones. Several studies have reported the identification and expression of the transcription factor family MADS-box members in cotton ovules and fibers; however, their roles are still elusive during the reproductive development in cotton. In this study, we evaluated the expression profiles of five MADS-box genes (GhMADS3, GhMADS4, GhMADS5, GhMADS6 and GhMADS7) belonging to the AGAMOUS-subfamily in Gossypium hirsutum. Phylogenetic and protein sequence analyses were performed using diploid (G. arboreum, G. raimondii) and tetraploid (G. barbadense, G. hirsutum) cotton genomes, as well as the AG-subfamily members from Arabidopsis thaliana, Petunia hybrida and Antirrhinum majus. qPCR analysis showed that the AG-subfamily genes had high expression during flower and fruit development in G. hirsutum. In situ hybridization analysis also substantiates the involvement of AG-subfamily members on reproductive tissues of G. hirsutum, including ovule and ovary. The effect of plant hormones on AG-subfamily genes expression was verified in cotton fruits treated with gibberellin, auxin and brassinosteroid. All the genes were significantly regulated in response to auxin, whereas only GhMADS3, GhMADS4 and GhMADS7 genes were also regulated by brassinosteroid treatment. In addition, we have investigated the GhMADS3 and GhMADS4 overexpression effects in Arabidopsis plants. Interestingly, the transgenic plants from both cotton AG-like genes in Arabidopsis significantly altered the fruit size compared to the control plants. This alteration suggests that cotton AG-like genes might act regulating fruit formation. Our results demonstrate that members of the AG-subfamily in G. hirsutum present a conserved expression profile during flower development, but also demonstrate their expression during fruit development and in response to phytohormones.