Prerpints.org logo
Preprint
Article

Functional Characterization of Candidate Genes, Gohir.D05G103700 and Gohir.D12G153600, Identified Through the Expression QTL Analysis using Virus Induced Gene Silencing in Upland Cotton (Gossypium hirsutum L.)

Altmetrics

Downloads

140

Views

91

Comments

0

A peer-reviewed article of this preprint also exists.

supplementary.docx (438.84KB )

This version is not peer-reviewed

Submitted:

07 April 2023

Posted:

07 April 2023

You are already at the latest version

Alerts
Abstract
Cotton (Gossypium spp.) is a major source of natural fiber and an important cash crop. The cotton growth habit and architecture determine its productivity and influence management strategies for commercial production. The GATA transcription factors (TFs) control various developmental processes in plants, such as flower, bract and embryo development, and petal differentiation. As stable transformation is still a bottleneck in many plant species, TRV-VIGS was used to manipulate gene expression in different plants, including Gossypium hirsutum L. In this study, we undertook the TRV-based VIGS to functionally characterize two candidate genes, Gohir.D05G103700 and Gohir.D12G153600, identified through the expression QTL analysis for five floral induction and meristem identity genes using the upland cotton mini-core collection. Virus-induced silencing of the Gohir.D05G103700 gene resulted in up to a 1.4-fold reduction in the transcript level in two inoculated plants, G3 and G4, and Gohir.D12G153600 gene resulted in up to a 2.3-fold reduction in transcript level in a single inoculated plant P05 relative to the mock-treated plant. The TRV2-Gohir.D05G103700 inoculated plants precisely G3 and G4 also exhibited loss of the supernumerary (fourth) floral bract in the squares, whereas the TRV2-Gohir.D12G153600 inoculated plants did not show any observable phenotypic change relative to the mock-treated plants.  Altogether, this study suggested that TRV-VIGS can be used to characterize genes in cotton relatively rapidly and the cotton Gohir.D05G103700 gene is a positive regulator of the indeterminate growth habit in cotton, which could be manipulated to obtain a cotton plant with architecture best suited for the cultivation area.
Keywords: 
Subject: Biology and Life Sciences  -   Agricultural Science and Agronomy

1. Introduction

The genus Gossypium emerged almost 12.5 million years ago (mya), and about 5-10 mya diverged into diploid A and D sub-genome species [1]. Subsequently, around 1-2 mya, the two-diploid species underwent trans-oceanic hybridization and chromosome doubling, resulting in an allopolyploid AD clade that radiated from the Americas and diversified into multiple species [1]. The allopolyploid cotton, Gossypium hirsutum L. (Upland cotton), and Gossypium barbadense L. (Egyptian cotton) have originated from two independent domestication events [2]. Today, G. hirsutum dominates the world’s cotton fiber production (90%). Whereas G. barbadense, G. arboretum, and G. herbaceum together contribute only 10% to it [3]. Cotton production generates an annual economic impact of around 600 billion US dollars globally, making it a major contributor to the world’s economy [4].
The plant architecture plays a crucial role in determining crop production and its management. It is determined by the meristematic cells and could be indeterminate, leading to continuous vegetative growth, to determinate, where meristematic cells are consumed in producing a terminal structure. Primary meristems are dynamically added in cotton-producing monopodial or sympodial branches, inflorescences, and terminal shoots, i.e., flowers. Similarly, the secondary meristems (cambium) add continuously and contribute to wood production. The position of organ-specific meristems called primordia contributes to the final plant architecture [5]. Flowering in cotton depends on the differentiation of flower primordia that depends on axillary bud primordia, which in turn depends on the shoot apical meristem (SAM) differentiation. Axillary bud primordia initiate the transition from vegetative to the reproductive phase in cotton by differentiating from leaf primordia to flower primordia at a specific time point during plant development.
Transcription factors (TF) regulate tissue-specific gene expression by binding to the promoter region of the target genes. Subsequently, the DNA-TF complex promotes the recruitment of the RNA polymerase onto the target sequence and promotes gene expression. The GATA TFs (one of the important family of transcription factors) in plants express at different developmental stages, especially the light-induced processes, such as embryo development, flowering, petal differentiation and expansion, and maturation [6,7,8]. The NTL1 was the first GATA TF gene identified in Nicotiana tabacum [9]. So far, over 30 GATA TFs have been identified and characterized in Arabidopsis thaliana [10]. Additionally, 64 GATA TFs have been identified in Glycine max, 30 GATA TFs in Gossypium spp., and 28 GATA TFs in Oryza sativa [10,11,12]. These studies highlight the importance of GATA TFs in plant growth and development.
Another diverse family of transcription factors, SQUMOSA PROMOTER BINDING PROTEIN-LIKE (SPL), plays a crucial role in plant growth and development [13,14]. The SPL transcription factors play various roles in the growth and development of plants, such as controlling vegetative to reproductive phase change in A. thaliana [15,16], ovary and fruit development in tomato [17], yield-related traits in wheat [18], and floral organ size and ovule production in cotton [19].
Virus-induced gene silencing (VIGS) is an important reverse genetics tool for in vivo gene function study in plants [20]. VIGS depends on posttranscriptional gene silencing (PTGS) machinery in a sequence-specific manner. Briefly, in this machinery, double-stranded RNAs (dsRNAs), such as the replication intermediates of plant RNA viruses and highly structured genomic RNA, are processed by Dicer-like proteins (DCLs) to form small interfering RNAs (siRNAs). The guide strand of siRNA gets loaded onto the RNA-induced silencing complex (RISC) and cleaves complementary messenger RNA (mRNA) [21]. For VIGS, fragments of the gene of interest are cloned into a viral vector and applied to plants where the endogenous RNA silencing machinery in the host has degraded the viral RNA and, in this process, produced complementary siRNAs that reduce the target gene expression. VIGS has been demonstrated using numerous plant-virus combinations [22,23]. In the past decade, several viral genomes have been modified as a powerful reverse genetic tool for the functional characterization of plant genes. Some examples of VIGS constructs are the Tobacco rattle virus (TRV) [24], Apple latent spherical virus (ALSV) [25,26], African cassava mosaic virus (ASMV) [27], Cucumber mosaic virus (CMV) [28], and Barley streak mosaic virus (BSMV) [29].
TRV consists of a bipartite genome, TRV1, and TRV2 [30], where TRV1 determines viral movement [31] and TRV2 codes for the coat protein and other nonstructural proteins. Also, it shows a great deal of variability among different viral isolates [30]. These nonstructural proteins are involved in the nematode transmission of this virus [32], but these proteins are not essential for plant infection. Therefore, TRV2’s nonstructural protein-coding genes were replaced by multiple cloning site where fragments of the target gene could be introduced [33]. TRV1 and TRV2 were later cloned into a binary vector for Agrobacterium-mediated plant delivery [24,33].
TRV-VIGS has been effectively used to manipulate gene expression in different plant species such as A. thaliana, Solanum lycopersicum, Nicotiana spp., G. arboreum, and Petunia hybrida [34]. TRV-based VIGS is one of the most effective methods to study gene function in plants due to several advantages, such as mild symptoms, spread to large cell patches, ability to migrate to the growing meristems, and systemic spread to new tissues and all plant parts [35,36,37].
Given the desirable attributes of TRV and its successful implementation in VIGS in different plant species, in this study, we undertook the TRV-based VIGS of two candidate genes, Gohir.D05G103700 and Gohir.D12G153600. Using upland cotton mini-core collection, we identified these genes earlier by expression QTL (eQTL) analysis for five floral induction and meristem identity genes (FT, SOC1, AP1, FUL, and LFY) using the upland cotton mini-core collection (Naveed et al. unpublished data). In this analysis, Gohir.D05G103700 showed an association with the APETALA 1 (AP1) gene expression trait and annotated as GATA transcription factor 11-like, and Gohir.D12G153600 showed an association with FRUITFUL (FUL) and SUPPRESSOR OF OVEREXPRESSION OF CONSTANS 1 (SOC1) gene expression traits and annotated as SQUAMOSA promoter binding-like transcription factor. This work is focused on characterizing the molecular function of these eQTLs and determining their role in cotton architecture development. Unveiling the molecular functions of these genes might have implications in breeding cotton for reduced regrowth after defoliation or annual growth habit with improved lint yield.

2. Materials and Methods

2.1. Plant husbandry

Cotton cultivar Coker 201 was selected for the VIGS experiment, as this cultivar was used in the past for genetic transformation experiments using Agrobacterium, which is also an intended method for delivering the TRV constructs here [38]. Seeds of cotton cultivar Coker 201 were germinated in 0.5-gallon plastic pots filled with an autoclaved 3:1 mixture of Fafard 3B and playground sand. The pots were placed in a walk-in growth chamber (GR96, Conviron, Pembina, ND) set at 30°C daytime and 20°C nighttime temperature with a 16h photoperiod. Three sets of six plants were placed in the growth chamber (one set each for two genes of interest and one as a mock control). All plants were treated alike regarding watering and fertilization for the entire length of the study, and the plant ruminants were discarded after autoclaving at the end of the study. The whole experiment was repeated once.

2.2. Construction of the gene silencing constructs

We selected two candidate genes, Gohir.D05G103700 and Gohir.D12G153600 (Table 1), which showed association with the cotton AP1 and FUL/SOC1 gene expression traits, respectively, for the functional characterization using the TRV-based VIGS system. Target sites for PCR primer designing from the cotton Gohir.D05G103700 and Gohir.D12G153600 genes were selected based on the conservation of sequences between the homoeologous genes and were determined using multiple sequence alignment (MSA) on the Clustal Omega software. Primers flanking the conserved regions in these genes were designed using the Primer 3 software (https://primer3.ut.ee/). The primers were used to amplify a 175 bp and a 166 bp fragment for Gohir.D05G103700 and Gohir.D12G153600 genes from cotton cultivar Coker 201 (for the list of primers, see Supplementary Table S1). For PCR cloning of these fragments in TRV2, primers were supplemented with the XbaI and KpnI restriction sites. After PCR, the product was digested with XbaI and KpnI and ligated with the TRV2 plasmid pYL156 linearized formerly using the same restriction enzymes. Subsequently, the construct was transformed into E. coli DH5α cells by heat shock method. The plasmid DNA was isolated and validated by PCR, restriction digestion, and DNA sequencing (McLAB, San Francisco, CA), and later whole plasmid sequence was obtained from Plasmidsaurus (Eugene, OR) (Supplemental Figure S1).

2.3. Agrobacterium tumefaciens transfection with TRV1 and TRV2 constructs and plant inoculation

The detailed structure of TRV1 and TRV2 constructs amenable for Agrobacterium-mediated gene silencing has been described previously [24,33,34]. The TRV2 plasmid containing the Gohir.D05G103700 gene sequence (pYL156-Gohir.D05G103700), hereafter referred to as TRV2-GATA-TF and pTRV1 were transfected into Agrobacterium tumefaciens strain GV3101 by the heat shock method and selected on plates supplemented with kanamycin (50μg/mL), gentamycin (50μg/mL), and rifampicin (10μg/mL). Likewise, a TRV2-SPL-TF and an empty viral control TRV2-MCS (mock control) were developed and transfected into A. tumefaciens strain GV3101. Subsequently, the leaves were agroinfiltrated with a 1:1 mixture of TRV1 and recombinant TRV2 constructs [34]. A schematic representation of these steps is presented in Figure 1. The leaf tissues were collected from the infiltrated plants 14 days after inoculation in liquid nitrogen and stored at -80°C until used.

2.4. RNA extraction and cDNA conversion

Total RNA was extracted from leaf tissue using the Spectrum™ Plant Total RNA Kit (Sigma-Aldrich, USA) following the manufacturer’s instructions with minor modifications. RNA concentrations were determined using NanoDrop™ (Thermo Scientific, USA). The total RNA (1μg) was converted to cDNA using the Revert Aid first strand cDNA synthesis kit (Thermo Scientific, USA) following the manufacturer’s instructions using the oligo d(T) primer.

2.5. Quantitative Real-Time PCR (qRT-PCR)

The quantitative real-time PCR (qRT-PCR) analysis was performed using the iTaq™ Universal SYBR® Green Supermix (BioRad, USA) following the manufacturer’s instructions on iCycler iQ™ (BioRad, USA). The following PCR profile was used for the gene expression analysis: initial denaturation at 95°C for 3 min followed by 50 cycles at 95°C for 10s and 56°C (Gohir.D05G103700 and ACTIN1) or 59.5°C (Gohir.D12G153600) for 30 s and 50 cycles of melt curve analysis at 56°C to 93.5°C for 10 s and final incubation at 4°C. The cotton ACTIN1 gene was used as an internal control for data normalization. The gene expression data were analyzed using the 2–∆∆Ct method and presented as log 2 of the fold change relative to the reference gene [39].

2.6. Phenotypic data analysis

The data on plant height, number of leaves per plant, number of squares per plant, and number of bracts per square were collected from the inoculated and mock control plants in the walk-in growth chamber (Supplementary Table S2). The statistical analysis of the phenotypic data was performed using the SAS package.

3. Results

3.1. Sequence analysis for the antisense construct development

The DNA sequence of the cotton Gohir.D05G103700 and Gohir.D12G153600 genes identified as eQTLs on chromosomes D05 (SNP markers i09222Gh and i00443Gh at 8730.32 Mb) and D12 (SNP marker i08185Gh at 48443.48 Mb) for the cotton AP1 and FUL/SOC1 gene expression traits respectively were retrieved from NCBI (National Center for Biotechnology Information) website. Both genes existed as homoeologous gene pairs. The former gene mapped to upland cotton chromosomes A05 (LOC107906238; location 8729.803-8733.656 Mb) and D05 (LOC107959021; location 9724.645-9728.563 Mb) and coded for a GATA transcription factor and the latter gene mapped to chromosomes A12 (LOC107934496; location 91532.20-91536.64 Mb) and D12 (LOC107946275; location 48443.14-48445.05 Mb) and coded for a SQUAMOSA PROMOTER BINDING LIKE protein (SPL). When the homoeologous gene sequences were aligned, these showed 95-96% similarity. A 175 bp antisense construct for GATA-TF and a 166 bp antisense construct for the SPL-TF genes were designed from sites showing 99-100% sequence similarity or conservation (Figure 2). The primer pairs were designed to amplify and PCR-clone products of expected size from the upland cotton cultivar ‘Coker-201’ in the TRV2 vector. The resulting plasmids were confirmed via DNA sequencing, and the maps were constructed. After sequence confirmation, the construct was transformed to the Agrobacterium tumefaciens strain GV3101. Likewise, the TRV1 construct was also transformed to the Agrobacterium tumefaciens GV3101 and bacteria with TRV1 and TRV2 were mixed at 1:1 proportion, and the plants were inoculated with the medley [34].

3.2. Agrobacterium mediated-VIGS of Gohir.D05G103700 and Gohir.D12G153600 genes in the upland Cotton Cultivar ‘Coker-201’

Three-week-old plants were inoculated with TRV2-GATA-TF plus TRV1 and in a separate experiment with TRV2-SPL-TF plus TRV1. Leaf samples were collected from inoculated and control plants two weeks after the treatment for RNA extraction and qRT-PCR analysis. The qRT-PCR was performed for the cotton Gohir.D05G103700 and Gohir.D12G153600 genes, and the cotton housekeeping ACTIN1 gene was used for the normalization of the gene expression data. The TRV2-GATA-TF showed an about 1.4-fold reduced abundance of the cotton Gohir.D05G103700 transcript in treated plants, G3 and G4 relative to the uninoculated control, hence, suggested downregulation of the gene (Figure 3A). Other inoculated plants, G2, G5, and G6 except G1 showed a negligible reduction in Gohir.D05G103700 transcript abundance, whereas, no silencing was observed in G1.
Likewise, qRT-PCR analysis of TRV2-SPL-TF inoculated plants showed up to a 2.3-fold reduction in Gohir.D12G153600 transcript abundance in P05 relative to uninoculated control (Figure 3B). Whereas, about 0.8-fold reduction in transcript abundance was observed in rest of the treated plants, P01, P04, and P06, except P02 that exhibited wild type expression level.

3.3. Expression profiling of the cotton AP1 genes in TRV2-GATA-TF inoculated plants and the FUL and SOC1 genes in TRV2-SPL-TF inoculated plants

Since the cotton Gohir.D05G103700 and Gohir.D12G153600 genes were identified as eQTLs for the cotton AP1 and FUL/SOC1 genes, respectively, we studied the effect of silencing these genes on the transcription levels of the cotton AP1 gene in TRV2-GATA-TF inoculated plants and the cotton FUL and SOC1 genes in TRV2-SPL-TF inoculated plants. Among the TRV2-GATA-TF inoculated plants, G3 exhibited a ~1.4-fold reduction in the transcript level of the cotton AP1 gene. In contrast, no reduction was observed in the case of G1, G4, G5, and G6, and a slight over-accumulation was observed in the case of G2 (Figure 3C).
Similarly, among the TRV2-SPL-TF inoculated plants, P01 and P03 showed up to a 1-fold reduction in the cotton FUL transcript level, P04 and P05 exhibited negligible reductions, and P02 showed no reduction. On the other hand, P06 showed over-accumulation (Figure 3D). The transcription level of the cotton SOC1, TRV2-SPL-TF inoculated plants showed up to a 1.2-fold reduction in P05 and ~ 1-fold reduction in P03. A negligible reduction in transcript abundance of cotton SOC1 was observed in P04. On the other hand, P01, P02, and P06 did not show any reduction in gene expression (Figure 3E).

3.4. Phenotypic observations

Data on various traits such as plant height, leaf and flower number per plant, and number and arrangement of floral organs (including the number of floral bracts per square) were recorded from the inoculated and mock control plants grown in the walk-in growth chamber (Supplementary Table S2). Interestingly, we observed loss of the supernumerary (fourth) floral bract in the squares of the TRV2-GATA-TF inoculated plants, G3 and G4 with a relatively higher level of Gohir.D05G103700 suppression (Figure 4). None of the other TRV2-GATA-TF inoculated plants and mock-treated plants exhibited suppression of the fourth bract in the floral buds of cotton cv. Coker-201. On the other hand, the TRV2-SPL-TF inoculated plants showed no observable phenotypic change relative to the mock-treated plants.

4. Discussion

Upland cotton an allopolyploid has a large and complex genome. Genetic transformation in this crop was successful and was the source of several herbicide tolerant and insect pest resistant genotypes available commercially [41]. However, genetic transformation is still not commonplace in cotton and is highly genotype dependent and time consuming [42]; hence cannot be applied to test gene functions. Virus-induced gene silencing offers a solution to this problem, as the experiments are relatively less time consuming. However, the gene-silencing effects are broadly transient, which has both pros and cons [43]. There are two viral systems available for gene silencing and gene delivery in upland cotton, Cotton leaf crumple virus (CLCrV) and Tobacco rattle virus (TRV) [43]. Both systems were shown to work well; however, TRV-based system is more commonly used and hence was adopted in the present study for functional characterization of Gohir.D05G103700 and Gohir.D12G153600.
Despite the fact that we used an upland cotton cultivar ‘Coker 201’ that earlier showed susceptibility to Agroinfiltration [38], hence consistent delivery and systemic spread of the recombinant virus, we observed plant-to-plant variation in the silencing efficiency in our experiments. For instance, five of the six plants showed some level of suppression; however, the level of suppression ranged from about 0.4-fold to 1.4-fold for Gohir.D05G103700 and from about 0.8-fold to a little over 2.2-fold for Gohir.D12G153600 in two and a single plant, respectively. A review of the literature suggested that variation in the level of gene suppression is not uncommon in VIGS experiments, as was documented in several earlier studies performed in different crop plants [44]. Several reasons were proposed to contribute to the observed variation in the level of gene suppression, from the effect of abiotic factors to genetic factors [45,46]. Further, in this study, we observed incomplete gene silencing for either of the genes. This incomplete silencing could be attributed to multiple factors, including the plant genotype, ploidy level, the plant development stages, temperature, and inoculation method [47,48,49]. Further, it was demonstrated that TRV-based gene silencing was more efficient in the diploid (G. arboreum and G. herbaceum) than in tetraploid cotton (G. hirsutum) [49]. Furthermore, the inoculation method is a critical factor for the successful delivery of the VIGS construct and subsequent steps involving systemic spread and gene silencing [29].
In addition, the two genes Gohir.D05G103700 and Gohir.D12G153600, differed in the level of suppression. There could be several possible reasons, including VIGS is a random process and may not lead to a similar level of suppression in two different attempts, as the long antisense RNA after complementation with the target RNAs process is somewhat random fashion. Also, the two homoeologous copies of a gene do not carry perfect similarity; hence, they do not exhibit similar levels of suppression. Unfortunately, it is quite difficult to develop homoeologue-specific primers to study the level of suppression of individual homoeologues. Hence, a cumulative level of gene suppression was studied in most of the studies [50,51]. In some cases, the antisense was able to silence both copies equally well, and some only one of the two copies. Additionally, often times a compensatory increase in the expression level of one of the homoeologues was observed upon silencing the other homoeologue [45].
Phenotypic analysis revealed that Coker 201 tends to produce a supernumerary bract (fourth bract) in each flower (visible in the mock control), which was absent from GATA-VIGS inoculated plants (Figure 4A). The earlier research suggested that cotton tends to produce a fourth bract in certain environmental conditions when the meristem in the axil of the true leaves is not consumed fully in flower production see Figure 4B [40]. It was proposed that the three bracts in the cotton bud are modifications of one true leaf and two stipules, and the fourth bract represents the second leaf (indicative of indeterminate growth). These results suggested that the GATA transcription factor identified as an eQTL for the cotton AP1 gene has a target site in the AP1 promoter, where it binds and occupies a site not letting the FT/FD assembly bind and induce/promote its expression, which triggers flowering and hence consumption of the meristem in floral development [12,52,53,54]. We hypothesized that VIGS of the GATA transcription factor released the site from suppression, allowing FT/FD-induced activation of the AP1 gene. The qRT-PCR analysis, however, did not exhibit AP1 over accumulation. Furthermore, the fourth bract in cotton is an undesirable trait [40], as its enlargement creates a gap allowing thrips and other insects to enter and feed on the developing fruit. Functional characterization of the Gohir.D05G103700 gene identified as the cause of the undesirable fourth floral bract and made available the target to select for this trait. However, we did not see any observable phenotypic changes in the TRV2-SPL-TF plants, which suggested the functional redundancy of this transcription factor in cotton or a need to completely silence the gene to have a notable phenotypic effect. Quantitative RT-PCR analysis of the AP1 and the SOC1/FUL genes established a connection between these genes and eQTLs, Gohir.D05G103700 and Gohir.D12G153600. In sum, this study reconfirmed the utility of TRV-based VIGS in cotton by functional characterization of an eQTL target for the AP1 expression trait. Also, the study reflected on the transcriptional regulation of the cotton AP1 gene and provided a target to breed for thrip resistance in cotton cultivars with a supernumerary bract (fourth bract).

Supplementary Materials

The following supporting information can be downloaded at the website of this paper posted on Preprints.org. Figure S1: Plasmid maps of TRV2-GATA-TF, A and TRV2-SPL-TF, B.; Table S1: List of primers used in this study; Table S2: Phenotypic data collected from Gohir.D05G103700, Gohir.D12G153600, and mock treated plants; Table S3: List of the Tobacco rattle virus (TRV) based Virus Induced Gene Silencing (VIGS) constructs.

Author Contributions

Conceptualization, S.R.; methodology, S.N.; experiment, S.N.; data analysis, S.N.; writing-original draft preparation, S.N.; visualization, S.N., and S.R.; Supervision, S.R.; writing-review and editing, S.N., and S.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Cotton Incorporated, grant number 2015969.

Data Availability Statement

The data presented in this paper are available on request from the corresponding author.

Acknowledgments

The authors would like to acknowledge the support for this work from the SC Cotton Board and Cotton Incorporated.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Paterson, A. H.; Wendel, J. F.; Gundlach, H.; Guo, H.; Jenkins, J.; Jin, D.; Llewellyn, D.; Showmaker, K. C.; Shu, S.; Udall, J.; Yoo, M.; Byers, R.; Chen, W.; Doron-Faigenboim, A.; Duke, M. V.; Gong, L.; Grimwood, J.; Grover, C.; Grupp, K.; Hu, G. Repeated Polyploidization of Gossypium Genomes and the Evolution of Spinnable Cotton Fibres. Nature 2012, 492, 423–427. [Google Scholar] [CrossRef]
  2. Renny-Byfield, S.; Page, J. T.; Udall, J. A.; Sanders, W. S.; Peterson, D. G.; Arick, M. A.; Grover, C. E.; Wendel, J. F. Independent Domestication of Two Old World Cotton Species. Genome Biology and Evolution 2016, 8, 1940–1947. [Google Scholar] [CrossRef]
  3. Zhang, H.-B.; Li, Y.; Wang, B.; Chee, P. W. Recent Advances in Cotton Genomics. International Journal of Plant Genomics 2008, 2008, 1–20. [Google Scholar] [CrossRef]
  4. Ahmad, S. (2014). Cotton Production and Uses Agronomy. Crop Protection and Postharvest Technologies. Springer Nature Singapore 2014. [Google Scholar] [CrossRef]
  5. Sussex, I. M.; Kerk, N. M. The Evolution of Plant Architecture. Current Opinion in Plant Biology 2001, 4, 33–37. [Google Scholar] [CrossRef]
  6. Borello, U.; Ceccarelli, E.; Giuliano, G. Constitutive, Light-Responsive and Circadian Clock-Responsive Factors Compete for the Different I Box Elements in Plant Light-Regulated Promoters. The Plant Journal 1993, 4, 611–619. [Google Scholar] [CrossRef]
  7. Terzaghi, W. B.; Cashmore, A. R. Light-Regulated Transcription. Annual Review of Plant Physiology and Plant Molecular Biology 1995, 46, 445–474. [Google Scholar] [CrossRef]
  8. Luo, X.-M.; Lin, W.-H.; Zhu, S.; Zhu, J.-Y.; Sun, Y.; Fan, X.-Y.; Cheng, M.; Hao, Y.; Oh, E.; Tian, M.; Liu, L.; Zhang, M.; Xie, Q.; Chong, K.; Wang, Z.-Y. Integration of Light- and Brassinosteroid-Signaling Pathways by a GATA Transcription Factor in Arabidopsis. Developmental Cell 2010, 19, 872–883. [Google Scholar] [CrossRef]
  9. Daniel-Vedele, F.; Caboche, M. A Tobacco CDNA Clone Encoding a GATA-1 Zinc Finger Protein Homologous to Regulators of Nitrogen Metabolism in Fungi. Molecular and General Genetics MGG 1993, 240, 365–373. [Google Scholar] [CrossRef]
  10. Reyes J., C.; Muro-Pastor, M. I.; Florencio, F. J. The GATA Family of Transcription Factors in Arabidopsis and Rice. Plant Physiology 2004, 134, 1718–1732. [Google Scholar] [CrossRef] [PubMed]
  11. Zhang, C.; Hou, Y.; Hao, Q.; Chen, H.; Chen, L.; Yuan, S.; Shan, Z.; Zhang, X.; Yang, Z.; Qiu, D.; Zhou, X.; Huang, W. Genome-Wide Survey of the Soybean GATA Transcription Factor Gene Family and Expression Analysis under Low Nitrogen Stress. PLOS ONE 2015, 10, e0125174. [Google Scholar] [CrossRef]
  12. Zhang, Z.; Zou, X.; Huang, Z.; Fan, S.; Qun, G.; Liu, A.; Gong, J.; Li, J.; Gong, W.; Shi, Y.; Fan, L.; Zhang, Z.; Liu, R.; Jiang, X.; Lei, K.; Shang, H.; Xu, A.; Yuan, Y. Genome-Wide Identification and Analysis of the Evolution and Expression Patterns of the GATA Transcription Factors in Three Species of Gossypium Genus. Gene 2019, 680, 72–83. [Google Scholar] [CrossRef]
  13. Klein, J.; Saedler, H.; Huijser, P. A New Family of DNA Binding Proteins Includes Putative Transcriptional Regulators of TheAntirrhinum Majus Floral Meristem Identity GeneSQUAMOSA. Molecular and General Genetics MGG 1996, 250, 7–16. [Google Scholar] [CrossRef]
  14. Yang, Z.; Wang, X.; Gu, S.; Hu, Z.; Xu, H.; Xu, C. Comparative Study of SBP-Box Gene Family in Arabidopsis and Rice. Gene 2007, 407, 1–11. [Google Scholar] [CrossRef] [PubMed]
  15. Wang, H.; Wang, H. The MiR156/SPL Module, a Regulatory Hub and Versatile Toolbox, Gears up Crops for Enhanced Agronomic Traits. Molecular Plant 2015, 8, 677–688. [Google Scholar] [CrossRef] [PubMed]
  16. Xu, M.; Hu, T.; Zhao, J.; Park, M.-Y.; Earley, K. W.; Wu, G.; Yang, L.; Poethig, R. S. Developmental Functions of MiR156-Regulated SQUAMOSA PROMOTER BINDING PROTEIN-like (SPL) Genes in Arabidopsis Thaliana. PLOS Genetics 2016, 12, e1006263. [Google Scholar] [CrossRef]
  17. Silva, G. F. F. e; Silva, E. M.; da Silva Azevedo, M.; Guivin, M. A. C.; Ramiro, D. A.; Figueiredo, C. R.; Carrer, H.; Peres, L. E. P.; Nogueira, F. T. S. MicroRNA156-Targeted SPL/SBP Box Transcription Factors Regulate Tomato Ovary and Fruit Development. The Plant Journal 2014, 78, 604–618. [Google Scholar] [CrossRef] [PubMed]
  18. Zhang, B.; Xu, W.; Liu, X.; Mao, X.; Li, A.; Wang, J.; Chang, X.; Zhang, X.; Jing, R. Functional Conservation and Divergence among Homoeologs of TaSPL20 and TaSPL21, Two SBP-Box Genes Governing Yield-Related Traits in Hexaploid Wheat1[OPEN]. Plant Physiology 2017, 174, 1177–1191. [Google Scholar] [CrossRef]
  19. Liu, N.; Tu, L.; Wang, L.; Hu, H.; Xu, J.; Zhang, X. MicroRNA 157-Targeted SPL Genes Regulate Floral Organ Size and Ovule Production in Cotton. BMC Plant Biology 2017, 17. [Google Scholar] [CrossRef]
  20. Hayward, A.; Padmanabhan, M.; Dinesh-Kumar, S. P. Virus-Induced Gene Silencing in Nicotiana Benthamiana and Other Plant Species. Methods in Molecular Biology (Clifton, N.J.) 2011, 678, 55–63. [Google Scholar] [CrossRef]
  21. Baulcombe, D. RNA Silencing in Plants. Nature 2004, 431, 356–363. [Google Scholar] [CrossRef] [PubMed]
  22. Burch-Smith, T. M.; Anderson, J. C.; Martin, G. B.; Dinesh-Kumar, S. P. Applications and Advantages of Virus-Induced Gene Silencing for Gene Function Studies in Plants. The Plant Journal 2004, 39, 734–746. [Google Scholar] [CrossRef] [PubMed]
  23. Kant, R.; Dasgupta, I. Gene Silencing Approaches through Virus-Based Vectors: Speeding up Functional Genomics in Monocots. Plant Molecular Biology 2019, 100, 3–18. [Google Scholar] [CrossRef] [PubMed]
  24. Liu, Y.; Schiff, M.; Marathe, R.; Dinesh-Kumar, S. P. Tobacco Rar1, EDS1 and NPR1/NIM1 like Genes Are Required for N-Mediated Resistance to Tobacco Mosaic Virus. The Plant Journal 2002, 30, 415–429. [Google Scholar] [CrossRef] [PubMed]
  25. Gedling, C. R.; Ali, E. M.; Gunadi, A.; Finer, J. J.; Xie, K.; Liu, Y.; Yoshikawa, N.; Qu, F.; Dorrance, A. E. Improved Apple Latent Spherical Virus-Induced Gene Silencing in Multiple Soybean Genotypes through Direct Inoculation of Agro-Infiltrated Nicotiana Benthamiana Extract. Plant Methods 2018, 14. [Google Scholar] [CrossRef]
  26. Igarashi, A.; Yamagata, K.; Sugai, T.; Takahashi, Y.; Sugawara, E.; Tamura, A.; Yaegashi, H.; Yamagishi, N.; Takahashi, T.; Isogai, M.; Takahashi, H.; Yoshikawa, N. Apple Latent Spherical Virus Vectors for Reliable and Effective Virus-Induced Gene Silencing among a Broad Range of Plants Including Tobacco, Tomato, Arabidopsis Thaliana, Cucurbits, and Legumes. Virology 2009, 386, 407–416. [Google Scholar] [CrossRef]
  27. Lentz, E. M.; Kuon, J.-E.; Alder, A.; Mangel, N.; Zainuddin, I. M.; McCallum, E. J.; Anjanappa, R. B.; Gruissem, W.; Vanderschuren, H. Cassava Geminivirus Agroclones for Virus-Induced Gene Silencing in Cassava Leaves and Roots. Plant Methods 2018, 14. [Google Scholar] [CrossRef]
  28. Tzean, Y.; Lee, M.-C.; Jan, H.-H.; Chiu, Y.-S.; Tu, T.-C.; Hou, B.-H.; Chen, H.-M.; Chou, C.-N.; Yeh, H.-H. Cucumber Mosaic Virus-Induced Gene Silencing in Banana. Scientific Reports 2019, 9. [Google Scholar] [CrossRef]
  29. Rustgi, S.; Naveed, S.; Windham, J.; Zhang, H.; Demirer, G. S. Plant Biomacromolecule Delivery Methods in the 21st Century. Frontiers in Genome Editing 2022, 4. [Google Scholar] [CrossRef]
  30. MacFarlane, S. A. Molecular Biology of the Tobraviruses. Journal of General Virology 1999, 80, 2799–2807. [Google Scholar] [CrossRef]
  31. Ziegler-Graff, V.; Guilford, P. J.; Baulcombe, D. C. Tobacco Rattle Virus RNA-1 29K Gene Product Potentiates Viral Movement and Also Affects Symptom Induction in Tobacco. Virology 1991, 182, 145–155. [Google Scholar] [CrossRef] [PubMed]
  32. Visser, P. B.; Bol, J. F. Nonstructural Proteins of Tobacco Rattle Virus Which Have a Role in Nematode-Transmission: Expression Pattern and Interaction with Viral Coat Protein. Journal of General Virology 1999, 80, 3273–3280. [Google Scholar] [CrossRef]
  33. Ratcliff, F.; Martin-Hernandez, A. M.; Baulcombe, D. C. Technical Advance: Tobacco Rattle Virus as a Vector for Analysis of Gene Function by Silencing. The Plant Journal 2001, 25, 237–245. [Google Scholar] [CrossRef] [PubMed]
  34. Senthil-Kumar, M.; Mysore, K. S. Tobacco Rattle Virus–Based Virus-Induced Gene Silencing in Nicotiana Benthamiana. Nature Protocols 2014, 9, 1549–1562. [Google Scholar] [CrossRef] [PubMed]
  35. Saedler, R. Virus-Induced Gene Silencing of Jasmonate-Induced Direct Defences, Nicotine and Trypsin Proteinase-Inhibitors in Nicotiana Attenuata. Journal of Experimental Botany 2003, 55, 151–157. [Google Scholar] [CrossRef]
  36. Galis, I.; Schuman, M. C.; Gase, K.; Hettenhausen, C.; Hartl, M.; Dinh, S. T.; Wu, J.; Bonaventure, G.; Baldwin, I. T. The Use of VIGS Technology to Study Plant-Herbivore Interactions. Methods in Molecular Biology (Clifton, N.J.) 2013, 975, 109–137. [Google Scholar] [CrossRef] [PubMed]
  37. Atamian, H.S. Virus induced gene silencing optimization in plants: things to be considered. Virus 2014, 2. [Google Scholar] [CrossRef]
  38. Jin, S.; Zhang, X.; Nie, Y.; Guo, X.; Liang, S.; Zhu, H. Identification of a Novel Elite Genotype for in Vitro Culture and Genetic Transformation of Cotton. Biologia plantarum 2006, 50, 519–524. [Google Scholar] [CrossRef]
  39. Ye, K.; Teng, T.; Yang, T.; Zhao, D.; Zhao, Y. Transcriptome Analysis Reveals the Effect of Grafting on Gossypol Biosynthesis and Gland Formation in Cotton. BMC Plant Biology 2023, 23. [Google Scholar] [CrossRef]
  40. Mauney, J. Chapter 2 VEGETATIVE GROWTH AND DEVELOPMENT OF FRUITING SITES. Available online: https://www.semanticscholar.org/paper/Chapter-2-VEGETATIVE-GROWTH-AND-DEVELOPMENT-OF-Mauney/7f0dcf43b3d71eaac0317c6690f8a9aedcb3989c.
  41. Ribeiro, T. P.; Lourenço-Tessutti, I. T.; de Melo, B. P.; Morgante, C. V.; Filho, A. S.; Lins, C. B. J.; Ferreira, G. F.; Mello, G. N.; Macedo, L. L. P.; Lucena, W. A.; Silva, M. C. M.; Oliveira-Neto, O. B.; Grossi-de-Sa, M. F. Improved Cotton Transformation Protocol Mediated by Agrobacterium and Biolistic Combined-Methods. Planta 2021, 254, 20. [Google Scholar] [CrossRef]
  42. Rajasekaran, K. Biolistic Transformation of Cotton Embryogenic Cell Suspension Cultures. Methods in Molecular Biology 2018, 55–66. [Google Scholar] [CrossRef]
  43. Zhang, J.; Wang, F.; Zhang, C.; Zhang, J.; Chen, Y.; Liu, G.; Zhao, Y.; Hao, F.; Zhang, J. A Novel VIGS Method by Agroinoculation of Cotton Seeds and Application for Elucidating Functions of GhBI-1 in Salt-Stress Response. Plant Cell Reports 2018, 37, 1091–1100. [Google Scholar] [CrossRef] [PubMed]
  44. Bennypaul, H. S.; Mutti, J. S.; Rustgi, S.; Kumar, N.; Okubara, P. A.; Gill, K. S. Virus-Induced Gene Silencing (VIGS) of Genes Expressed in Root, Leaf, and Meiotic Tissues of Wheat. Functional & Integrative Genomics 2012, 12, 143–156. [Google Scholar] [CrossRef]
  45. Chaudhary, B.; Flagel, L.; Stupar, R. M.; Udall, J. A.; Verma, N.; Springer, N. M.; Wendel, J. F. Reciprocal Silencing, Transcriptional Bias and Functional Divergence of Homeologs in Polyploid Cotton (Gossypium). Genetics 2009, 182, 503–517. [Google Scholar] [CrossRef] [PubMed]
  46. Coate, J. E.; Doyle, J. J. Quantifying Whole Transcriptome Size, a Prerequisite for Understanding Transcriptome Evolution across Species: An Example from a Plant Allopolyploid. Genome Biology and Evolution 2010, 2, 534–546. [Google Scholar] [CrossRef]
  47. Gao, X.; Britt Jr., R.C.; Shan, L.; He, P. Agrobacterium-Mediated Virus-Induced Gene Silencing Assay in Cotton. Journal of Visualized Experiments 2011, 54. [Google Scholar] [CrossRef]
  48. Tuttle, J. R.; Idris, A. M.; Brown, J. K.; Haigler, C. H.; Robertson, D. Geminivirus-Mediated Gene Silencing from Cotton Leaf Crumple Virus Is Enhanced by Low Temperature in Cotton. Plant Physiology 2008, 148, 41–50. [Google Scholar] [CrossRef]
  49. Mustafa, R.; Shafiq, M.; Mansoor, S.; Briddon, R. W.; Scheffler, B. E.; Scheffler, J.; Amin, I. Virus-Induced Gene Silencing in Cultivated Cotton (Gossypium Spp.) Using Tobacco Rattle Virus. Molecular Biotechnology 2015, 58, 65–72. [Google Scholar] [CrossRef]
  50. Lloyd, A.; Blary, A.; Charif, D.; Charpentier, C.; Tran, J.; Balzergue, S.; Delannoy, E.; Rigaill, G.; Jenczewski, E. Homoeologous Exchanges Cause Extensive Dosage-Dependent Gene Expression Changes in an Allopolyploid Crop. New Phytologist 2017, 217, 367–377. [Google Scholar] [CrossRef]
  51. Glover, N. M.; Redestig, H.; Dessimoz, C. Homoeologs: What Are They and How Do We Infer Them? Trends in Plant Science 2016, 21, 609–621. [Google Scholar] [CrossRef]
  52. Zhao, Y.; Medrano, L.; Ohashi, K.; Fletcher, J. C.; Yu, H.; Sakai, H.; Meyerowitz, E. M. HANABA TARANU Is a GATA Transcription Factor That Regulates Shoot Apical Meristem and Flower Development in Arabidopsis[W]. The Plant Cell 2004, 16, 2586–2600. [Google Scholar] [CrossRef] [PubMed]
  53. Wang, L.; Yin, H.; Qian, Q.; Yang, J.; Huang, C.; Hu, X.; Luo, D. NECK LEAF 1, a GATA Type Transcription Factor, Modulates Organogenesis by Regulating the Expression of Multiple Regulatory Genes during Reproductive Development in Rice. Cell Research 2009, 19, 598–611. [Google Scholar] [CrossRef] [PubMed]
  54. Whipple, C. J.; Hall, D. H.; DeBlasio, S.; Taguchi-Shiobara, F.; Schmidt, R. J.; Jackson, D. P. A Conserved Mechanism of Bract Suppression in the Grass Family. The Plant Cell 2010, 22, 565–578. [Google Scholar] [CrossRef] [PubMed]
Preprints 70654 g001
Figure 1. Schematic representation of the VIGS process, from gene cloning to plant infiltration (A). Diagrammatic representation of the TRV1 and TRV2 plasmids, specifically highlighting the restriction sites used to clone the gene fragments (B). Agarose gels showing the PCR products (175 bp Gohir.D05G103700 fragment and 166 bp Gohir.D12G153600 fragment) amplified using gene-specific primers (Table S1) for cloning in the TRV2 plasmid (C).
Figure 1. Schematic representation of the VIGS process, from gene cloning to plant infiltration (A). Diagrammatic representation of the TRV1 and TRV2 plasmids, specifically highlighting the restriction sites used to clone the gene fragments (B). Agarose gels showing the PCR products (175 bp Gohir.D05G103700 fragment and 166 bp Gohir.D12G153600 fragment) amplified using gene-specific primers (Table S1) for cloning in the TRV2 plasmid (C).
Preprints 70654 g001
Preprints 70654 g002
Figure 2. Multiple sequence alignment of the target sites in the cotton Gohir.D05G103700 (A), Gohir.D12G153600 (B) gene sequences showing conservation between homoeologues.
Figure 2. Multiple sequence alignment of the target sites in the cotton Gohir.D05G103700 (A), Gohir.D12G153600 (B) gene sequences showing conservation between homoeologues.
Preprints 70654 g002
Preprints 70654 g003
Figure 3. A quantitative real-time PCR (qRT-PCR) analysis of the Gohir.D05G103700 gene (A), and the GhAP1 gene (B) in TRV2-GATA-TF inoculated plants, and the Gohir.D12G153600 (C), GhFUL (D) and GhSOC1 (E) genes in TRV2-SPL-TF inoculated plants. The expression data were analyzed using the 2–∆∆Ct method and normalized to the cotton ACTIN1 gene. Note: Different letters above error bars indicate significant differences (P<0.05).
Figure 3. A quantitative real-time PCR (qRT-PCR) analysis of the Gohir.D05G103700 gene (A), and the GhAP1 gene (B) in TRV2-GATA-TF inoculated plants, and the Gohir.D12G153600 (C), GhFUL (D) and GhSOC1 (E) genes in TRV2-SPL-TF inoculated plants. The expression data were analyzed using the 2–∆∆Ct method and normalized to the cotton ACTIN1 gene. Note: Different letters above error bars indicate significant differences (P<0.05).
Preprints 70654 g003
Preprints 70654 g004
Figure 4. Virus-induced gene silencing of the cotton Gohir.D05G103700 gene, identified as an expression QTL for the AP1 expression trait, A. Bottom left - diagrammatic representation of the cotton meristems that lead to the flower bud production, B. M = meristem, B = bracts [40]. Bottom middle, a picture of three and four bract squares. The square on the right is a typical three bract square. The square on the left has four well form bracts. The fourth bract (arrow) on the square in the center is not well-formed (very similar to the condition in Coker 201). Bottom right - the hypothetical models proposed to explain the phenotype observed for the GATA-VIGS silenced plants. Note: the presence of the fourth bract is an undesirable trait, as the gap produced by the expansion of the fourth bract allows thrips to enter and feed on the developing ovary/fruit [40].
Figure 4. Virus-induced gene silencing of the cotton Gohir.D05G103700 gene, identified as an expression QTL for the AP1 expression trait, A. Bottom left - diagrammatic representation of the cotton meristems that lead to the flower bud production, B. M = meristem, B = bracts [40]. Bottom middle, a picture of three and four bract squares. The square on the right is a typical three bract square. The square on the left has four well form bracts. The fourth bract (arrow) on the square in the center is not well-formed (very similar to the condition in Coker 201). Bottom right - the hypothetical models proposed to explain the phenotype observed for the GATA-VIGS silenced plants. Note: the presence of the fourth bract is an undesirable trait, as the gap produced by the expansion of the fourth bract allows thrips to enter and feed on the developing ovary/fruit [40].
Preprints 70654 g004
Table 1. List of the candidate genes identified earlier by expression QTL analysis for FT, SOC1, AP1, FUL, and LFY genes using upland cotton mini-core collection (Naveed et al. unpublished data).
Table 1. List of the candidate genes identified earlier by expression QTL analysis for FT, SOC1, AP1, FUL, and LFY genes using upland cotton mini-core collection (Naveed et al. unpublished data).
Sr. No. Molecular markers Gene ID Sub-genome Gene Name Position on genome (Mb) Annotation
1 i02927Gh Gohir.A01G208700 A01 AP1, FT, LFY 117.19592 Trihelix transcription factor PTL
2 i43992Gh Gohir.A08G034500 A08 FT 4.35425 MYB3-like transcription factor
3 i13158Gh Gohir.A13G050400 A13 FT, LFY 7.06423 GATA transcription factor 28-like
4* i09222Gh; i00443Gh Gohir.D05G103700 D05 AP1 8.73032 GATA transcription factor 11-like
5* i08185Gh Gohir.D12G153600 D12 FUL, SOC1 48.44348 SQUAMOSA promoter binding-like transcription factor
6 i13848Gh; i13851Gh Gohir.D13G236200 D13 LFY 64.67230 Homeobox-leucine zipper protein REVOLUTA-like
Note: The candidate genes selected for virus induced gene silencing (VIGS) are marked with an asterisk (*).
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Functional Characterization of Candidate Genes, Gohir.D05G103700 and Gohir.D12G153600, Identified Through the Expression QTL Analysis using Virus Induced Gene Silencing in Upland Cotton (Gossypium hirsutum L.)

Salman Naveed

et al.

,

2023

Comparative Transcriptomics of Non-Embryogenic and Embryogenic Callus in Semi-Recalcitrant and Non-Recalcitrant Upland Cotton Lines

Sonika Kumar

et al.

,

2021

Genome Wide Identification of Cotton C-Repeat Binding Factor (CBF) and Overexpression of Gthu17439 (GthCBF4) Gene Confer Cold Stress Tolerance in Arabidopsis Thaliana

Jiangna Liu

et al.

,

2021

Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

facebook logotwitter logolinkedin logo
MDPI Initiatives
Important Links
Subscribe

Disclaimer

Terms of Use

Privacy Policy

© 2024 MDPI (Basel, Switzerland) unless otherwise stated