MAC3 protein analogs in upland cotton genome
Using amino acid sequences of AtMAC3A and AtMAC3B as queries to search against cotton genome (https://www.cottongen.org/), six analogous proteins (Gh_A01G0193, Gh_D01G2267, Gh_A02G0417, Gh_D02G0470, Gh_A09G2001 and Gh_D09G2202) were found in upland cotton genome and named as GhMAC3a~GhMAC3f according to the chromosome on which each gene resides. Phylogenetic analysis revealed that cotton MAC3 homologous proteins belong to the same clade and are closely related to those of soybean and Arabidopsis, but relatively far away from MACs in rice, maize and wheat (Fig. 1A). Multiple sequence alignment analysis showed that GhMAC3e contain a conserved U-box domain, Prp19 domain, and WD40 domain, respectively, and the MAC3s of cotton, soybean and Arabidopsis exhibit a high degree of similarity, suggesting that the homologous proteins of MAC3 are evolutionarily conserved and possibly perform similar functions in dicotyledonous plants. GhMAC3e (Gh_A09G2001) consists of a 1, 575 bp open reading frame, encoding 524 amino acids.
This protein was originally identified in a fiberless cotton mutant by proteomics analysis, its differential accumulation in the wild type and fiberless mutant ovules suggested that this gene is relevant to fiber development(Liu, et al. 2012), and therefore was molecularly cloned for further function analysis.
Expression pattern of GhMAC3e
To evaluate the possible biological functions of MAC3 genes in cotton, we analyzed the spatiotemporal expression characteristics of GhMAC3s and changes in response to abiotic and biotic stresses, respectively, by using public RNA-seq database of G. hirsutum (Bio Project Accession: PRJNA490626). The results showed that GhMAC3a is predominantly expressed in stamen and the ovules at -3 days post-anthesis (DPA); GhMAC3d is highly expressed in stem; GhMAC3b, GhMAC3c, GhMAC3e, and GhMAC3f exhibited similar expression patterns between different tissues, with relatively higher levels in stem and the ovules before anthesis (Fig. 2A). Six GhMAC3s can be induced to express by salt (NaCl), PEG (Polyethylene Glycol), and VD strain Vd6, and roughly in a similar pattern (Fig. 2B). These results indicate that GhMAC3s may be involved in cotton growth, ovule development, and responses to different stresses.
As phytohormones play crucial roles in regulating plant growth and stress tolerance, we further analyzed if expression of GhMAC3e could be induced by exogenous GA3 (representing growth hormone), and stress hormones by using RT-qPCR (Fig. 3). The results showed GhMAC3e could be significantly induced by all above phytohormone treatments: The expression of GhMAC3s was up-regulated by GA3 after 3 hours post treatment (hpt) with a peak at 6 hpt (Fig. 3B); under ABA treatment, the abundance of GhMAC3e mRNAs at 1, 3, and 12 hpt dropped to a level that even lower than at 0 hpt, but peaked at 24 hpt to be 2.8 fold higher than that at 0 hpt (Fig. 3C); expression of GhMAC3e can be significantly induced by SA, MeJA, and ETH as early as 1 hpt, and up-regulated up to 27-, 23-, and 14- fold higher levels compared to that at 0 htp, respectively (Fig. 3D-F). These results indicate that GhMAC3e expression is swiftly response to both exogenous GA3 and stress hormones, suggesting that GhMAC3e may play regulatory roles both in growth and stress responses.
Silencing GhMAC3e compromised primary stem growth and biomass accumulation in cotton
To elucidate the biological function of GhMAC3e, we employed virus-induced gene silencing (VIGS) technology to knock down GhMAC3e in cotton. In this system, silencing GhCLA1 to produce an albino phenotype was set as a visible marker to monitor the efficiency of VIGS, cotton leaves began to display the expected albino phenotype at around 14 d post Agro-infiltration with the TRV1 and TRV2:GhCLA1 (Fig. S1A). So we harvested TRV2:00 (empty vector control, WT) and TRV2:GhMAC3e plants at 14 days post-infiltration to determine the expression level of GhMAC3e by semi-quantitative RT-PCR, respectively. The results showed that GhMAC3e mRNA level was significantly reduced in TRV2:GhMAC3 plants (Fig. S1B), indicating that GhMAC3e expression was effectively silenced. The plant heights of GhMAC3e-silenced plants were significantly reduced by 25.7% compared to the control plants (Fig.4A, B), notably, this significant reduction attributed to the main stem length loss (Fig.4C), because no significant difference was observed in hypocotyl length between GhMAC3e-silenced plants and the control plants (Fig. 4D). Moreover, both the seedling fresh weight and dry weight were also significantly reduced by silencing GhMAC3e (Fig. 4 E, F).
On the longitudinal section along the stem center of the epicotyl of the seedlings, it can be observed that after silencing GhMAC3e, the number of the cortical cells per unit length of stem was increased coupling with significantly reduced longitudinal length of each of them (Fig. 4 G, H). These observations suggest GhMAC3e is required for plant growth, especially for cell elongation.
To explore the physiological and molecular mechanisms underlying the arrested stem apical growth in GhMAC3e-silenced cotton seedlings, firstly, the auxin concentrations in stem tips were extracted and quantified, and the auxin distribution in the petiole of the top leaves was directly visualized using DR5-GFP reporter system. The results showed that silencing GhMAC3e significantly reduced auxin concentration in the stem tips of cotton seedlings, meanwhile, the green fluorescence in the petiole phloem of was dramatically impaired (Fig. 5A-C). These results suggest that both auxin biosynthesis in stem tips and auxin transportation between stem and young leaves through petiole phloem could be attenuated by the suppressed expression of GhMAC3e, thereupon then block cotton growth and biomass accumulation.
To explore the relationship between GhMAC3e, auxin, and cotton growth, we investigated the influences of silencing GhMAC3e on the expression of GhXTH16 (xyloglucan endo-glucan transglucosylase/hydrolase 16), GhGH3.6 (Gretchen Hagen3.6), and GhSAUR12 (small auxin up RNA 12), which plays a crucial role in plant cell wall extension, auxin inactivation, and auxin responsive cell elongation, respectively. Substantial upregulation of GhXTH16 and Gh3.6 expression and downregulation of GhSAUR12 were observed in GhMAC3e-silenced cotton (Fig. 5D). Our findings indicate that GhMAC3e is involved in homeostasis of active and inactive auxin, auxin signaling, and cell wall extension, which play important roles in plant growth.
Overexpression of GhMAC3e enhances disease resistance of Verticillium wilt in Arabidopsis
To further elucidate the function of GhMAC3e, the wild-type Col-0 Arabidopsis plants were genetically modified using Agrobacterium tumefaciens to create transgenic lines overexpressing GhMAC3e driven by the CaMV35S promoter. Among 17 independent transgenic lines, two pure lines, named 35S::GhMAC3e #23 and 35S::GhMAC3e #40, were selected for further analysis, because each line has one single transgene insertion and significantly elevated expression level of GhMAC3e (Fig. S2A, B). The 10-day-old transgenic and wild-type Arabidopsis plants (WT) were inoculated with VD strain Vd991 to evaluate VW disease resistance. At 10 days post-inoculation (dpi), the Arabidopsis leaves exhibited characteristic chlorosis, wilting and necrosis symptoms (Fig 6A). In comparison to the wild-type Arabidopsis plants, these symptoms were much weaker in overexpression transgenic lines coupling with significantly reduced fungal biomass (Fig. 6B), even though the typical disease symptoms became more pronounced in both the wild type and the transgenic lines as inoculation duration increased. So we examined the survival rates to evaluate durable resistance, the overexpression line 40# with the highest expression of GhMAC3e was selected for survival rate assessment. By 19 dpi, most leaves in the wild-type plants exhibited chlorosis and wilting, and even some individual plants died completely, in contrast, much fewer yellow leaves showed on the transgenic plants (Fig. 6C). The survival rate of the WT plants started to decline at 18 dpi and completely decreased to 0 at 27 dpi, while that of the transgenic plants began to decline at 22 dpi and decreased to 0 at 36 dpi. Overexpression of GhMAC3e can delay Arabidopsis death caused by infection with VD991 (Fig. 6D), thereby gaining more persistent tolerance to the fungus. Taken together, the ectopic overexpression of GhMAC3e conferred Arabidopsis enhanced resistance to VD. To explore the possible molecular mechanisms underlying this enhancement in VW resistance, the expression of AtNPR1, AtPR1, AtPR1, AtMYC2, AtVSP2, AtTGA2, AtERF1, and AtAAO3 of WT and the transgenic line 40# in response to VD infection (2 dpi) was analysed using RT-qPCR. As shown in Fig. 7, overexpression of GhMAC3e did not alter AtERF1 expression, implying that the advanced VW resistance may not employ ethylene-signaling regulation pathway. AtMYC2 and AtVSP2, the two important members in the same JA signaling branch, were down-regulated in the transgenic plants after VD inoculation; AtAAO3, a key gene in ABA biosynthesis pathway, was also down-regulated in transgenic plants both under regular and VD infection condition. The transgenic GhMAC3e can slightly decrease the expression level of AtNPR1 and AtTGA2 without VD infection, but did not influence their expressions after VD infection. Increased expression of the AtPR1 and AtPR2 genes marks the activation of SA signaling pathway, VD inoculation increased AtPR1 expression by more than 146 fold in WT, ectopic expression of GhMAC3e dramatically upregulated AtPR1 expression by around 15 fold and up to more than 262 fold under basal and VD inoculation condition, respectively. Whereas, AtPR2 expression was moderately induced by VD, but was compromised by overexpression of GhMAC3e. Together, ectopic expression of GhMAC3e seems to regulate SA rather than JA, ABA, ET, or reactive oxygen species (ROS) signaling in response to VD infection. Furthermore, GhMAC3e regulating SA signaling is independent upon GhNPR1 transcription level.