Genetic basis for glandular trichome formation in cotton

Isolation of Cottonseed Extracts That Affect Human Cancer Cell Growth

Cottonseeds are classified as glanded or glandless seeds depending on the presence or absence of gossypol glands. Glanded cottonseed has anticancer property and glandless cottonseed was reported to cause cancer in one animal study. It is important to investigate the effect of bioactive components from cottonseeds. Our objectives were to isolate ethanol extracts from cottonseeds and investigate their effects on human cancer cells. A protocol was developed for isolating bioactive extracts from seed coat and kernel of glanded and glandless cottonseeds. HPLC-MS analyzed the four ethanol extracts but only quercetin was identified in the glandless seed coat extract. Residual gossypol was detected in the glanded and glandless seed kernel extracts and but only in the glanded seed coat extract. Ethanol extracts were used to treat human cancer cells derived from breast and pancreas followed by MTT assay for cell viability. Ethanol extracts from glanded and glandless cottonseed kernels and gossypol significantly decreased breast cancer cell mitochondrial activity. Ethanol extract from glanded cottonseed kernel and gossypol also significantly decreased pancreas cancer cell mitochondrial activity. These results suggest that ethanol extracts from cottonseeds, like gossypol, contain anticancer activities.

Characterization of gossypol biosynthetic pathway

Gossypol and related sesquiterpene aldehydes in cotton function as defense compounds but are antinutritional in cottonseed products. By transcriptome comparison and coexpression analyses, we identified 146 candidates linked to gossypol biosynthesis. Analysis of metabolites accumulated in plants subjected to virus-induced gene silencing (VIGS) led to the identification of four enzymes and their supposed substrates. In vitro enzymatic assay and reconstitution in tobacco leaves elucidated a series of oxidative reactions of the gossypol biosynthesis pathway. The four functionally characterized enzymes, together with (+)-δ-cadinene synthase and the P450 involved in 7-hydroxy-(+)-δ-cadinene formation, convert farnesyl diphosphate (FPP) to hemigossypol, with two gaps left that each involves aromatization. Of six intermediates identified from the VIGS-treated leaves, 8-hydroxy-7-keto-δ-cadinene exerted a deleterious effect in dampening plant disease resistance if accumulated. Notably, CYP71BE79, the enzyme responsible for converting this phytotoxic intermediate, exhibited the highest catalytic activity among the five enzymes of the pathway assayed. In addition, despite their dispersed distribution in the cotton genome, all of the enzyme genes identified show a tight correlation of expression. Our data suggest that the enzymatic steps in the gossypol pathway are highly coordinated to ensure efficient substrate conversion.

Combinatorial Transcriptional Control of Plant Specialized Metabolism

Plants produce countless specialized compounds of diverse chemical nature and biological activities. Their biosynthesis often exclusively occurs either in response to environmental stresses or is limited to dedicated anatomical structures. In both scenarios, regulation of biosynthesis appears to be mainly controlled at the transcriptional level, which is generally dependent on a combined interplay of DNA-related mechanisms and the activity of transcription factors that may act in a combinatorial manner. How environmental and developmental cues are integrated into a coordinated cell type-specific stress response has only partially been unraveled so far. Building on the available examples from (metabolic) gene expression, here we propose theoretical models of how this integration of signals may occur at the level of transcriptional control.

Divergence and evolution of cotton bHLH proteins from diploid to allotetraploid

BACKGROUND

Polyploidy is considered a major driving force in genome expansion, yielding duplicated genes whose expression may be conserved or divergence as a consequence of polyploidization.

RESULTS

We compared the genome sequences of tetraploid cotton (Gossypium hirsutum) and its two diploid progenitors, G. arboreum and G. raimondii, and found that the bHLH genes were conserved over the polyploidization. Oppositely, the expression of the homeolgous gene pairs was diversified. The biased homeologous proportion for bHLH family is significantly higher (64.6%) than the genome wide homeologous expression bias (40%). Compared with cacao (T. cacao), orthologous genes only accounted for a small proportion (41.7%) of whole cotton bHLHs family. The further Ks analysis indicated that bHLH genes underwent at least two distinct episodes of whole genome duplication: a recent duplication (1.0-60.0 million years ago, MYA, 0.005 < Ks < 0.312) and an old duplication (> 60.0 MYA, 0.312 < Ks < 3.0). The old duplication event might have played a key role in the expansion of the bHLH family. Both recent and old duplicated pairs (68.8%) showed a divergent expression profile, indicating specialized functions. The expression diversification of the duplicated genes suggested it might be a universal feature of the long-term evolution of cotton.

CONCLUSIONS

Overview of cotton bHLH proteins indicated a conserved and divergent evolution from diploids to allotetraploid. Our results provided an excellent example for studying the long-term evolution of polyploidy.

SSR-based association mapping of fiber quality in upland cotton using an eight-way MAGIC population

The quality of fiber is significant in the upland cotton industry. As complex quantitative traits, fiber quality traits are worth studying at a genetic level. To investigate the genetic architecture of fiber quality traits, we conducted an association analysis using a multi-parent advanced generation inter-cross (MAGIC) population developed from eight parents and comprised of 960 lines. The reliable phenotypic data for six major fiber traits of the MAGIC population were collected from five environments in three locations. Phenotypic analysis showed that the MAGIC lines have a wider variation amplitude and coefficient than the founders. A total of 284 polymorphic SSR markers among eight parents screened from a high-density genetic map were used to genotype the MAGIC population. The MAGIC population showed abundant genetic variation and fast linkage disequilibrium (LD) decay (0.76 cM, r² > 0.1), which revealed the advantages of high efficiency and power in QTL exploration. Association mapping via a mixed linear model identified 52 significant loci associated with six fiber quality traits; 14 of them were mapped in reported QTL regions with fiber-related or other agronomic traits. Nine markers demonstrated the pleiotropism that controls more than two fiber traits. Furthermore, two SSR markers, BNL1231 and BNL3452, were authenticated as hotspots that were mapped with multi-traits. In addition, we provided candidate regions and screened six candidate genes for identified loci according to the LD decay distance. Our results provide valuable QTL for further genetic mapping and will facilitate marker-based breeding for fiber quality in cotton.

Cotton functional genomics reveals global insight into genome evolution and fiber development

Due to the economic value of natural textile fiber, cotton has attracted much research attention, which has led to the publication of two diploid genomes and two tetraploid genomes. These big data facilitate functional genomic study in cotton, and allow researchers to investigate cotton genome structure, gene expression, and protein function on the global scale using high-throughput methods. In this review, we summarized recent studies of cotton genomes. Population genomic analyses revealed the domestication history of cultivated upland cotton and the roles of transposable elements in cotton genome evolution. Alternative splicing of cotton transcriptomes was evaluated genome-widely. Several important gene families like MYC, NAC, Sus and GhPLDα1 were systematically identified and classified based on genetic structure and biological function. High-throughput proteomics also unraveled the key functional proteins correlated with fiber development. Functional genomic studies have provided unprecedented insights into global-scale methods for cotton research.

Population structure and genetic basis of the agronomic traits of upland cotton in China revealed by a genome-wide association study using high-density SNPs

Gossypium hirsutum L. represents the largest source of textile fibre, and China is one of the largest cotton-producing and cotton-consuming countries in the world. To investigate the genetic architecture of the agronomic traits of upland cotton in China, a diverse and nationwide population containing 503 G. hirsutum accessions was collected for a genome-wide association study (GWAS) on 16 agronomic traits. The accessions were planted in four places from 2012 to 2013 for phenotyping. The CottonSNP63K array and a published high-density map based on this array were used for genotyping. The 503 G. hirsutum accessions were divided into three subpopulations based on 11 975 quantified polymorphic single-nucleotide polymorphisms (SNPs). By comparing the genetic structure and phenotypic variation among three genetic subpopulations, seven geographic distributions and four breeding periods, we found that geographic distribution and breeding period were not the determinants of genetic structure. In addition, no obvious phenotypic differentiations were found among the three subpopulations, even though they had different genetic backgrounds. A total of 324 SNPs and 160 candidate quantitative trait loci (QTL) regions were identified as significantly associated with the 16 agronomic traits. A network was established for multieffects in QTLs and interassociations among traits. Thirty-eight associated regions had pleiotropic effects controlling more than one trait. One candidate gene, Gh_D08G2376, was speculated to control the lint percentage (LP). This GWAS is the first report using high-resolution SNPs in upland cotton in China to comprehensively investigate agronomic traits, and it provides a fundamental resource for cotton genetic research and breeding.

Plant Glandular Trichomes: Natural Cell Factories of High Biotechnological Interest

Multicellular glandular trichomes are epidermal outgrowths characterized by the presence of a head made of cells that have the ability to secrete or store large quantities of specialized metabolites. Our understanding of the transcriptional control of glandular trichome initiation and development is still in its infancy. This review points to some central questions that need to be addressed to better understand how such specialized cell structures arise from the plant protodermis. A key and unique feature of glandular trichomes is their ability to synthesize and secrete large amounts, relative to their size, of a limited number of metabolites. As such, they qualify as true cell factories, making them interesting targets for metabolic engineering. In this review, recent advances regarding terpene metabolic engineering are highlighted, with a special focus on tobacco (Nicotiana tabacum). In particular, the choice of transcriptional promoters to drive transgene expression and the best ways to sink existing pools of terpene precursors are discussed. The bioavailability of existing pools of natural precursor molecules is a key parameter and is controlled by so-called cross talk between different biosynthetic pathways. As highlighted in this review, the exact nature and extent of such cross talk are only partially understood at present. In the future, awareness of, and detailed knowledge on, the biology of plant glandular trichome development and metabolism will generate new leads to tap the largely unexploited potential of glandular trichomes in plant resistance to pests and lead to the improved production of specialized metabolites with high industrial or pharmacological value.

Rapid mapping and cloning of the virescent-1 gene in cotton by bulked segregant analysis-next generation sequencing and virus-induced gene silencing strategies

Map-based gene cloning is a vital strategy for the identification of the quantitative trait loci or genes underlying important agronomic traits. The conventional map-based cloning method is powerful but generally time-consuming and labor-intensive. In this context, we introduce an improved bulked segregant analysis method in combination with a virus-induced gene silencing (VIGS) strategy for rapid and reliable gene mapping, identification and functional verification. This method was applied to a multiple recessive marker line of upland cotton, Texas 582 (T582), and identified unique genomic positions harboring mutant loci, showing the reliability and efficacy of this method. The v1 locus was further fine-mapped. Only one gene, GhCHLI, which encodes one of the subunits of Mg chelatase, was differentially down-regulated in T582 compared with TM-1. A point mutation occurred in the AAA+ conserved region of GhCHLI and led to an amino acid substitution. Suppression of its expression by VIGS in TM-1 resulted in a yellow blade phenotype that was similar to T582. This integrated approach provides a paradigm for the rapid mapping and identification of the candidate genes underlying the genetic traits in plants with large and complex genomes in the future.

Genome-wide comparative transcriptome analysis of CMS-D2 and its maintainer and restorer lines in upland cotton

BACKGROUND

Cytoplasmic male sterility (CMS) conferred by the cytoplasm from Gossypium harknessii (D2) is an important system for hybrid seed production in Upland cotton (G. hirsutum). The male sterility of CMS-D2 (i.e., A line) can be restored to fertility by a restorer (i.e., R line) carrying the restorer gene Rf1 transferred from the D2 nuclear genome. However, the molecular mechanisms of CMS-D2 and its restoration are poorly understood.

RESULTS

In this study, a genome-wide comparative transcriptome analysis was performed to identify differentially expressed genes (DEGs) in flower buds among the isogenic fertile R line and sterile A line derived from a backcross population (BC8F1) and the recurrent parent, i.e., the maintainer (B line). A total of 1464 DEGs were identified among the three isogenic lines, and the Rf1-carrying Chr_D05 and its homeologous Chr_A05 had more DEGs than other chromosomes. The results of GO and KEGG enrichment analysis showed differences in circadian rhythm between the fertile and sterile lines. Eleven DEGs were selected for validation using qRT-PCR, confirming the accuracy of the RNA-seq results.

CONCLUSIONS

Through genome-wide comparative transcriptome analysis, the differential expression profiles of CMS-D2 and its maintainer and restorer lines in Upland cotton were identified. Our results provide an important foundation for further studies into the molecular mechanisms of the interactions between the restorer gene Rf1 and the CMS-D2 cytoplasm.

Salicylic acid-related cotton (Gossypium arboreum) ribosomal protein GaRPL18 contributes to resistance to Verticillium dahliae

BACKGROUND

Verticillium dahliae is a phytopathogenic fungal pathogen that causes vascular wilt diseases responsible for considerable decreases in cotton yields. The complex mechanism underlying cotton resistance to Verticillium wilt remains uncharacterized. Identifying an endogenous resistance gene may be useful for controlling this disease.

RESULTS

We cloned the ribosomal protein L18 (GaRPL18) gene, which mediates resistance to Verticillium wilt, from a wilt-resistant cotton species (Gossypium arboreum). We then characterized the function of this gene in cotton and Arabidopsis thaliana plants. GaRPL18 encodes a 60S ribosomal protein subunit important for intracellular protein biosynthesis. However, previous studies revealed that some ribosomal proteins are also inhibitory toward oncogenesis and congenital diseases in humans and play a role in plant disease defense. Here, we observed that V. dahliae infections induce GaRPL18 expression. Furthermore, we determined that the GaRPL18 expression pattern is consistent with the disease resistance level of different cotton varieties. GaRPL18 expression is upregulated by salicylic acid (SA) treatments, suggesting the involvement of GaRPL18 in the SA signal transduction pathway. Virus-induced gene silencing technology was used to determine whether the GaRPL18 expression level influences cotton disease resistance. Wilt-resistant cotton species in which GaRPL18 was silenced became more susceptible to V. dahliae than the control plants because of a significant decrease in the abundance of immune-related molecules. We also transformed A. thaliana ecotype Columbia (Col-0) plants with GaRPL18 according to the floral dip method. The plants overexpressing GaRPL18 were more resistant to V. dahliae infections than the wild-type Col-0 plants. The enhanced resistance of transgenic A. thaliana plants to V. dahliae is likely mediated by the SA pathway.

CONCLUSION

Our findings provide new insights into the role of GaRPL18, indicating that it plays a crucial role in resistance to cotton "cancer", also known as Verticillium wilt, mainly regulated by an SA-related signaling pathway mechanism.

HOMEODOMAIN PROTEIN 1 is required for jasmonate-mediated glandular trichome initiation in Artemisia annua

Glandular trichomes are generally considered biofactories that produce valuable chemicals. Increasing glandular trichome density is a very suitable way to improve the productivity of these valuable metabolites, but little is known about the regulation of glandular trichome formation. Phytohormone jasmonate (JA) promotes glandular trichome initiation in various plants, but its mechanism is also unknown. By searching transcription factors regulated by JA in Artemisia annua, we identified a novel homeodomain-leucine zipper transcription factor, HOMEODOMAIN PROTEIN 1 (AaHD1), which positively controls both glandular and nonglandular trichome initiations. Overexpression of AaHD1 in A. annua significantly increased glandular trichome density without harming plant growth. Consequently, the artemisinin content was improved. AaHD1 interacts with A. annua jasmonate ZIM-domain 8 (AaJAZ8), which is a repressor of JA, thereby resulting in decreased transcriptional activity. AaHD1 knockdown lines show decreased sensitivity to JA on glandular trichome initiation, which indicates that AaHD1 plays an important role in JA-mediated glandular trichome initiation. We identified a new transcription factor that promotes A. annua glandular trichome initiation and revealed a novel molecular mechanism by which a homeodomain protein transduces JA signal to promote glandular trichome initiation. Our results also suggested a connection between glandular and nonglandular trichome formations.

Targeting insect mitochondrial complex I for plant protection

Plant engineered to express double-stranded RNA (dsRNA) can target the herbivorous insect gene for silencing. Although mounting evidence has emerged to support feasibility of this new pest control technology, field application is slow largely due to lack of potent targets. Here, we show that suppression of the gene encoding NDUFV2, a subunit of mitochondrial complex I that catalyses NADH dehydrogenation in respiratory chain, was highly lethal to insects. Feeding cotton bollworm (Helicoverpa armigera) larvae with transgenic cotton tissues expressing NDUFV2 dsRNA led to mortality up to 80% within 5 days, and almost no larvae survived after 7 days of feeding, due to the altered mitochondrial structure and activity. Transcriptome comparisons showed a drastic repression of dopa decarboxylase genes. Reciprocal assays with Asian corn borer (Ostrinia furnacalis), another lepidopteran species, revealed the sequence-specific effect of NDUFV2 suppression. Furthermore, the hemipteran lugus Apolygus lucorum was also liable to NDUFV2 repression. These data demonstrate that the mitochondrial complex I is a promising target with both sequence specificity and wide applicability for the development of new-generation insect-proof crops.