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 LC₅₀ 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.

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

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.

A nonproteinaceous Fusarium cell wall extract triggers receptor-like protein-dependent immune responses in Arabidopsis and cotton

Fusarium wilt caused by the ascomycete fungus Fusarium oxysporum is a devastating disease of many economically important crops. The mechanisms underlying plant responses to F. oxysporum infections remain largely unknown. We demonstrate here that a water-soluble, heat-resistant and nonproteinaceous F. oxysporum cell wall extract (FoCWE) component from multiple F. oxysporum isolates functions as a race-nonspecific elicitor, also termed pathogen-associated molecular pattern (PAMP). FoCWE triggers several demonstrated immune responses, including mitogen-activated protein (MAP) kinase phosphorylation, reactive oxygen species (ROS) burst, ethylene production, and stomatal closure, in cotton and Arabidopsis. Pretreated FoCWE protects cotton seeds against infections by virulent F. oxysporum f. sp. vasinfectum (Fov), and Arabidopsis plants against the virulent bacterium, Pseudomonas syringae, suggesting the potential application of FoCWEs in crop protection. Host-mediated responses to FoCWE do not appear to require LYKs/CERK1, BAK1 or SOBIR1, which are commonly involved in PAMP perception and/or signalling. However, FoCWE responses and Fusarium resistance in cotton partially require two receptor-like proteins, GhRLP20 and GhRLP31. Transcriptome analysis suggests that FoCWE preferentially activates cell wall-mediated defence, and Fov has evolved virulence mechanisms to suppress FoCWE-induced defence. These findings suggest that FoCWE is a classical PAMP that is potentially recognised by a novel pattern-recognition receptor to regulate cotton resistance to Fusarium infections.

Identification of Players Controlling Meristem Arrest Downstream of the FRUITFULL-APETALA2 Pathway

The end of the reproductive phase in monocarpic plants is determined by a coordinated arrest of all active meristems, a process known as global proliferative arrest (GPA). GPA is linked to the correlative control exerted by developing seeds and, possibly, the establishment of strong source-sink relationships. It has been proposed that the meristems that undergo arrest at the end of the reproductive phase behave at the transcriptomic level as dormant meristems, with low mitotic activity and high expression of abscisic acid response genes. Meristem arrest is also controlled genetically. In Arabidopsis (Arabidopsis thaliana), the MADS-box transcription factor FRUITFULL induces GPA by directly repressing genes of the APETALA2 (AP2) clade. The AP2 genes maintain shoot apical meristem (SAM) activity in part by keeping WUSCHEL expression active, but the mechanisms downstream of this pathway remain elusive. To identify target genes, we performed a transcriptomic analysis, inducing AP2 activity in meristems close to arrest. Our results suggest that AP2 controls meristem arrest by repressing genes related to axillary bud dormancy in the SAM and negative regulators of cytokinin signaling. In addition, our analysis indicates that genes involved in the response to environmental signals also respond to AP2, suggesting that it could modulate the end of flowering by controlling responses to both endogenous and exogenous signals. Our results support the previous observation that at the end of the reproductive phase the arrested SAM behaves as a dormant meristem, and they strongly support AP2 as a master regulator of this process.

Physiological and molecular responses of Setaria viridis to osmotic stress

Drought-tolerant species, such as Setaria viridis, a C₄ model plant, make physiological and biochemical adjustments water limitation and recover from the stress upon its release. We investigated S. viridis (A10.1 accession) responses to continuing osmotic stress. The osmotic stress was imposed using polyethylene glycol (PEG) 8000 (7.5%) for 10 days. Morphological traits and stomatal conductance were measured daily for the 10 days. On days 6 and 10, the following traits were measured separately for root and shoot: relative water content (RWC), osmotic potential (OP), electrolytic leakage (EL), and proline content. qPCR analysis was used to evaluate the expression of five selected genes in roots (SvLEA, SvDREB1C, SvPIP2-1, SvHSP20, and SvP5CS2), and chlorophyll a fluorescence was measured on three key days. The morphological data demonstrated a drastic reduction in shoot biomass as an effect of water deficit caused by the osmotic stress. Shoot biomass reduction could be associated with putative ABA-dependent signaling involved in SvDREB1C expression. Stomatal conductance and photosynthesis were severely affected up until day 6, however, stomatal conductance and some photosynthetic parameters such as F_V/F_M, ABS/RC, and DI₀/RC showed total or slight recovery on day 10. Root EL decreased in treated plants suggesting an investment in membrane protection by osmoregulator expression such as dehydrin (SvLEA) and proline (SvP5CS2) genes. Our data suggest that S. viridis exhibited a partial recovery from an imposed and constant osmotic stress within 10 days.

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

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.

Genome-wide analysis of the MADS-box gene family in polyploid cotton (Gossypium hirsutum) and in its diploid parental species (Gossypium arboreum and Gossypium raimondii)

The MADS-box gene family encodes transcription factors that share a highly conserved domain known to bind to DNA. Members of this family control various processes of development in plants, from root formation to fruit ripening. In this work, a survey of diploid (Gossypium raimondii and Gossypium arboreum) and tetraploid (Gossypium hirsutum) cotton genomes found a total of 147, 133 and 207 MADS-box genes, respectively, distributed in the MIKC, Mα, Mβ, Mγ, and Mδ subclades. A comparative phylogenetic analysis among cotton species, Arabidopsis, poplar and grapevine MADS-box homologous genes allowed us to evaluate the evolution of each MADS-box lineage in cotton plants and identify sequences within well-established subfamilies. Chromosomal localization and phylogenetic analysis revealed that G. raimondii and G. arboreum showed a conserved evolution of the MIKC subclade and a distinct pattern of duplication events in the Mα, Mγ and Mδ subclades. Additionally, G. hirsutum showed a combination of its parental subgenomes followed by a distinct evolutionary history including gene gain and loss in each subclade. qPCR analysis revealed the expression patterns of putative homologs in the AP1, AP3, AGL6, SEP4, AGL15, AG, AGL17, TM8, SVP, SOC and TT16 subfamilies of G. hirsutum. The identification of putative cotton orthologs is discussed in the light of evolution and gene expression data from other plants. This analysis of the MADS-box genes in Gossypium species opens an avenue to understanding the origin and evolution of each gene subfamily within diploid and polyploid species and paves the way for functional studies in cotton species.

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

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.