RNA Interference for Functional Genomics and Improvement of Cotton (Gossypium sp.)

Omics-driven exploration and mining of key functional genes for the improvement of food and fiber crops

The deployment of omics technologies has obtained an incredible boost over the past few decades with the advances in next-generation sequencing (NGS) technologies, innovative bioinformatics tools, and the deluge of available biological information. The major omics technologies in the limelight are genomics, transcriptomics, proteomics, metabolomics, and phenomics. These biotechnological advances have modernized crop breeding and opened new horizons for developing crop varieties with improved traits. The genomes of several crop species are sequenced, and a huge number of genes associated with crucial economic traits have been identified. These identified genes not only provide insights into the understanding of regulatory mechanisms of crop traits but also decipher practical grounds to assist in the molecular breeding of crops. This review discusses the potential of omics technologies for the acquisition of biological information and mining of the genes associated with important agronomic traits in important food and fiber crops, such as wheat, rice, maize, potato, tomato, cassava, and cotton. Different functional genomics approaches for the validation of these important genes are also highlighted. Furthermore, a list of genes discovered by employing omics approaches is being represented as potential targets for genetic modifications by the latest genome engineering methods for the development of climate-resilient crops that would in turn provide great impetus to secure global food security.

Genome Editing and Improvement of Abiotic Stress Tolerance in Crop Plants

Genome editing aims to revolutionise plant breeding and could assist in safeguarding the global food supply. The inclusion of a 12-40 bp recognition site makes mega nucleases the first tools utilized for genome editing and first generation gene-editing tools. Zinc finger nucleases (ZFNs) are the second gene-editing technique, and because they create double-stranded breaks, they are more dependable and effective. ZFNs were the original designed nuclease-based approach of genome editing. The Cys2-His2 zinc finger domain's discovery made this technique possible. Clustered regularly interspaced short palindromic repeats (CRISPR) are utilized to improve genetics, boost biomass production, increase nutrient usage efficiency, and develop disease resistance. Plant genomes can be effectively modified using genome-editing technologies to enhance characteristics without introducing foreign DNA into the genome. Next-generation plant breeding will soon be defined by these exact breeding methods. There is abroad promise that genome-edited crops will be essential in the years to come for improving the sustainability and climate-change resilience of food systems. This method also has great potential for enhancing crops' resistance to various abiotic stressors. In this review paper, we summarize the most recent findings about the mechanism of abiotic stress response in crop plants and the use of the CRISPR/Cas mediated gene-editing systems to improve tolerance to stresses including drought, salinity, cold, heat, and heavy metals.

NM, Jebakani K. S, Selvaraj S, Pramitha J. L, Selvaraj R, Petchiammal K. I, Kather Sheriff S, Thinakaran J, Rathinamoorthy S, Kumar P. R. Genetic manipulation of anti-nutritional factors in major crops for a sustainable diet in future

The consumption of healthy food, in order to strengthen the immune system, is now a major focus of people worldwide and is essential to tackle the emerging pandemic concerns. Moreover, research in this area paves the way for diversification of human diets by incorporating underutilized crops which are highly nutritious and climate-resilient in nature. However, although the consumption of healthy foods increases nutritional uptake, the bioavailability of nutrients and their absorption from foods also play an essential role in curbing malnutrition in developing countries. This has led to a focus on anti-nutrients that interfere with the digestion and absorption of nutrients and proteins from foods. Anti-nutritional factors in crops, such as phytic acid, gossypol, goitrogens, glucosinolates, lectins, oxalic acid, saponins, raffinose, tannins, enzyme inhibitors, alkaloids, β-N-oxalyl amino alanine (BOAA), and hydrogen cyanide (HCN), are synthesized in crop metabolic pathways and are interconnected with other essential growth regulation factors. Hence, breeding with the aim of completely eliminating anti-nutrition factors tends to compromise desirable features such as yield and seed size. However, advanced techniques, such as integrated multi-omics, RNAi, gene editing, and genomics-assisted breeding, aim to breed crops in which negative traits are minimized and to provide new strategies to handle these traits in crop improvement programs. There is also a need to emphasize individual crop-based approaches in upcoming research programs to achieve smart foods with minimum constraints in future. This review focuses on progress in molecular breeding and prospects for additional approaches to improve nutrient bioavailability in major crops.

A novel tasi RNA-based micro RNA-induced gene silencing strategy to tackle multiple pests and pathogens in cotton (Gossypium hirsutum L.)

This study demonstrates the combinatorial management of multiple pests through a trans-acting siRNA (tasiRNA)-based micro RNA-induced gene silencing (MIGS) strategy. Transgenic cotton events demonstrated improved efficacy against cotton leaf curl disease, cotton leaf hopper and root-knot nematode. Cotton (Gossypium hirsutum L.), an important commercial crop grown worldwide is confronted by several pests and pathogens, thus reiterating interventions for their management. In this study, we report, the utility of a novel Arabidopsis miRNA173-directed trans-acting siRNA (tasiRNA)-based micro RNA-induced gene silencing (MIGS) strategy for the simultaneous management of cotton leaf curl disease (CLCuD), cotton leaf hopper (CLH; Amrasca biguttula biguttula) and root-knot nematode (RKN, Meloidogyne incognita). Cotton transgenics were developed with the MIGS construct targeting a total of 7 genes by an apical meristem-targeted in planta transformation strategy. Stable transgenics were selected using stringent selection pressure, molecular characterization and stress-specific bio-efficacy studies. We identified 8 superior events with 50-100% resistance against CLCuD, while reduction in the root-knot nematode multiplication factor in the range of 35-75% confirmed resistance to RKN. These transgenic cotton events were also detrimental to the growth and development of CLH, as only 43.3-62.5% of nymphs could survive. Based on the corroborating evidences obtained by all the bioefficacy analyses, 3 events viz., L-75-1, E-27-11, E-27-7 were found to be consistent in tackling the target pests. To the best of our knowledge, this report is the first of its kind demonstrating the possibility of combinatorial management of pests/diseases in cotton using MIGS approach. These identified events demonstrate immense utility of the strategy towards combinatorial stress management in cotton improvement programs.

Data mining of transcriptional biomarkers at different cotton fiber developmental stages

Advancement of the gene expression study provides comprehensive information on pivotal genes at different cotton fiber development stages. For the betterment of cotton fiber yield and their quality, genetic improvement is a major target point for the cotton community. Therefore, various studies were carried out to understand the transcriptional machinery of fiber leading to the detailed integrative as well as innovative study. Through data mining and statistical approaches, we identified and validated the transcriptional biomarkers for staged specific differentiation of fiber. With the unique mapping read matrix of ~ 200 cotton transcriptome data and sequential statistical analysis, we identified several important genes that have a deciding and specific role in fiber cell commitment, initiation and elongation, or secondary cell wall synthesis stage. Based on the importance score and validation analysis, IQ domain 26, Aquaporin, Gibberellin regulated protein, methionine gamma lyase, alpha/beta hydrolases, and HAD-like superfamily have shown the specific and determining role for fiber developmental stages. These genes are represented as transcriptional biomarkers that provide a base for molecular characterization for cotton fiber development which will ultimately determine the high yield.

Functional Markers for Precision Plant Breeding

Advances in molecular biology including genomics, high-throughput sequencing, and genome editing enable increasingly faster and more precise cultivar development. Identifying genes and functional markers (FMs) that are highly associated with plant phenotypic variation is a grand challenge. Functional genomics approaches such as transcriptomics, targeting induced local lesions in genomes (TILLING), homologous recombinant (HR), association mapping, and allele mining are all strategies to identify FMs for breeding goals, such as agronomic traits and biotic and abiotic stress resistance. The advantage of FMs over other markers used in plant breeding is the close genomic association of an FM with a phenotype. Thereby, FMs may facilitate the direct selection of genes associated with phenotypic traits, which serves to increase selection efficiencies to develop varieties. Herein, we review the latest methods in FM development and how FMs are being used in precision breeding for agronomic and quality traits as well as in breeding for biotic and abiotic stress resistance using marker assisted selection (MAS) methods. In summary, this article describes the use of FMs in breeding for development of elite crop cultivars to enhance global food security goals.

Comparative transcriptome analysis between inbred and hybrids reveals molecular insights into yield heterosis of upland cotton

BACKGROUND

Utilization of heterosis has greatly improved the productivity of many crops worldwide. Understanding the potential molecular mechanism about how hybridization produces superior yield in upland cotton is critical for efficient breeding programs.

RESULTS

In this study, high, medium, and low hybrids varying in the level of yield heterosis were screened based on field experimentation of different years and locations. Phenotypically, high hybrid produced a mean of 14% more seed cotton yield than its better parent. Whole-genome RNA sequencing of these hybrids and their four inbred parents was performed using different tissues of the squaring stage. Comparative transcriptomic differences in each hybrid parent triad revealed a higher percentage of differentially expressed genes (DEGs) in each tissue. Expression level dominance analysis identified majority of hybrids DEGs were biased towards parent like expressions. An array of DEGs involved in ATP and protein binding, membrane, cell wall, mitochondrion, and protein phosphorylation had more functional annotations in hybrids. Sugar metabolic and plant hormone signal transduction pathways were most enriched in each hybrid. Further, these two pathways had most mapped DEGs on known seed cotton yield QTLs. Integration of transcriptome, QTLs, and gene co-expression network analysis discovered genes Gh_A03G1024, Gh_D08G1440, Gh_A08G2210, Gh_A12G2183, Gh_D07G1312, Gh_D08G1467, Gh_A03G0889, Gh_A08G2199, and Gh_D05G0202 displayed a complex regulatory network of many interconnected genes. qRT-PCR of these DEGs was performed to ensure the accuracy of RNA-Seq data.

CONCLUSIONS

Through genome-wide comparative transcriptome analysis, the current study identified nine key genes and pathways associated with biological process of yield heterosis in upland cotton. Our results and data resources provide novel insights and will be useful for dissecting the molecular mechanism of yield heterosis in cotton.

Role of MicroRNAs and small RNAs in regulation of developmental processes and agronomic traits in Gossypium species

Small RNAs (sRNAs) are short, non-coding, 17-24 nucleotides long RNA molecules that play vital roles in regulating gene expression in every known organism investigated to date including cotton (Gossypium ssp.). These tiny RNA molecules target diverse categories of genes from different bioliogical and metabolic processes and have been reported in the three domains of life. Small RNAs, including miRNAs, are involved in ovule and fiber development, biotic and abiotic stresses, fertility, and other biochemical processes in cotton species. Also, sRNAs are the critical components in RNA interference pathway. In this article, we have reviewed the research efforts related to the isolation and characterization of miRNAs using molecular and genomic approaches. The progress made in understanding the functional roles of miRNAs in regulation, alteration, and inactivation of fundamental plant processes and traits of importance in cotton are presented here.

ABP9, a maize bZIP transcription factor, enhances tolerance to salt and drought in transgenic cotton

ABP9 , encoding a bZIP transcription factor from maize, enhances tolerance to multiple stresses and may participate in the ABA signaling pathway in transgenic cotton by altering physiological and biochemical processes and stress-related gene expression. Abiotic stresses, such as soil salinity and drought, negatively affect growth, development, and yield in cotton. Gene ABP9, which encodes a bZIP transcription factor, binds to the abscisic acid (ABA)-responsive-element (ABRE2) motif of the maize catalase1 gene. Its expression significantly improves tolerance in Arabidopsis to multiple abiotic stresses, but little is known about its role in cotton. In the present study, the ABP9 gene was introduced into upland cotton (Gossypium hirsutum L.) cultivar R15 by Agrobacterium tumefaciens-mediated transformation, and 12 independent transgenic cotton lines were obtained. Cotton plants over-expressing ABP9 have enhanced tolerance to salt and osmotic stress. Under stress, they developed better root systems in a greenhouse and higher germination, reduced stomatal aperture, and stomatal density in a growth chamber. Under drought conditions, survival rate and relative water content (RWC) of transgenic cotton were higher than those of R15 plants. Under salt and osmotic stresses, chlorophyll, proline, and soluble sugar contents significantly increased in transgenic cotton leaves and the malondialdehyde (MDA) content was lower than in R15. Overexpression of ABP9 also enhanced oxidative stress tolerance, reduced cellular levels of reactive oxygen species (ROS) through increased activities of antioxidative enzymes, and alleviated oxidative damage to cell. Interestingly, ABP9 over-expressing cotton was more sensitive to exogenous ABA than R15 at seed germination, root growth, stomatal aperture, and stomatal density. Moreover, ABP9 overexpression upregulated significantly the transcription levels of stress-related genes such as GhDBP2, GhNCED2, GhZFP1, GhERF1, GhHB1, and GhSAP1 under salt treatment. Conjointly, these results showed that overexpression of ABP9 conferred enhanced tolerance to multiple abiotic stresses in cotton. The stress-tolerant transgenic lines provide valuable resources for cotton breeding.

RNAi of Grp78 may disturb the fusion of ICR mouse palate cultured in vitro

RNA interference (RNAi) is a powerful tool to silence or minimize gene expression, and palate culture in vitro is an important technique for study of the palate development. Our previous study demonstrated that the gene expression of glucose-regulated protein-78 (Grp78) was downregulation in the all-trans retinoic acid-induced mouse models of cleft palate (CP) during embryogenesis. To find the role of Grp78, the small interfering RNA (siRNA) of this gene carried by fluorescent vector was injected with a microinjector, through which about 30 pmol siRNA was injected into the Institute of Cancer Research (ICR) mouse palate explants. After 6, 12, 24, 48, and 72 h, these palate explants were removed from culture to observe their fluorescent and Alcian blue-staining phenotypes, and the expression of the unfolded protein response (UPR) key members (Grp78, Inositol-responsive enzyme 1, protein kinase RNA-like endoplasmic reticulum kinase, activating transcription factor-6 and X-box binding protein-1) was measured. After cultured for 72 h, the partially or completely fused bilateral palates were observed in the control siRNA group, while CPs were found in the Grp78 siRNA group. In the Grp78 siRNA group, the relatively mRNA abundance of the key genes belonged to UPR at each time point was lower than that of the control siRNA group, and their protein expression also displayed the same change. By the system of RNAi strategies with mouse palate culture, we found the siRNA of Grp78 disturbed the fusion of mouse palate cultured in vitro.

Genome-wide identification and characterization of microRNAs differentially expressed in fibers in a cotton phytochrome A1 RNAi line

Cotton fiber is an important commodity throughout the world. Fiber property determines fiber quality and commercial values. Previous studies showed that silencing phytochrome A1 gene (PHYA1) by RNA interference in Upland cotton (Gossypium hirsutum L. cv. Coker 312) had generated PHYA1 RNAi lines with simultaneous improvements in fiber quality (longer, stronger and finer fiber) and other key agronomic traits. Characterization of the altered molecular processes in these RNAi genotypes and its wild-type controls is a great interest to better understand the PHYA1 RNAi phenotypes. In this study, a total of 77 conserved miRNAs belonging to 61 families were examined in a PHYA1 RNAi line and its parental Coker 312 genotype by using multiplex sequencing. Of these miRNAs, seven (miR7503, miR7514, miR399c, miR399d, miR160, miR169b, and miR2950) were found to be differentially expressed in PHYA1 RNAi cotton. The target genes of these differentially expressed miRNAs were involved in the metabolism and signaling pathways of phytohormones, which included Gibberellin, Auxin and Abscisic Acid. The expression of several MYB transcription factors was also affected by miRNAs in RNAi cotton. In addition, 35 novel miRNAs (novel miR1-novel miR35) were identified in fibers for the first time in this study. Target genes of vast majority of these novel miRNAs were also predicted. Of these, nine novel miRNAs (novel-miR1, 2, 16, 19, 26, 27, 28, 31 and 32) were targeted to cytochrome P450-like TATA box binding protein (TBP). The qRT-PCR confirmed expression levels of several differentially regulated miRNAs. Expression patterns of four miRNAs-targets pairs were also examined via RNA deep sequencing. Together, the results imply that the regulation of miRNA expression might confer to the phenotype of the PHYA1 RNAi line(s) with improved fiber quality.

QTL mapping for flowering-time and photoperiod insensitivity of cotton Gossypium darwinii Watt

Most wild and semi-wild species of the genus Gossypium are exhibit photoperiod-sensitive flowering. The wild germplasm cotton is a valuable source of genes for genetic improvement of modern cotton cultivars. A bi-parental cotton population segregating for photoperiodic flowering was developed by crossing a photoperiod insensitive irradiation mutant line with its pre-mutagenesis photoperiodic wild-type G. darwinii Watt genotype. Individuals from the F2 and F3 generations were grown with their parental lines and F1 hybrid progeny in the long day and short night summer condition (natural day-length) of Uzbekistan to evaluate photoperiod sensitivity, i.e., flowering-time during the seasons 2008-2009. Through genotyping the individuals of this bi-parental population segregating for flowering-time, linkage maps were constructed using 212 simple-sequence repeat (SSR) and three cleaved amplified polymorphic sequence (CAPS) markers. Six QTLs directly associated with flowering-time and photoperiodic flowering were discovered in the F2 population, whereas eight QTLs were identified in the F3 population. Two QTLs controlling photoperiodic flowering and duration of flowering were common in both populations. In silico annotations of the flanking DNA sequences of mapped SSRs from sequenced cotton (G. hirsutum L.) genome database has identified several potential 'candidate' genes that are known to be associated with regulation of flowering characteristics of plants. The outcome of this research will expand our understanding of the genetic and molecular mechanisms of photoperiodic flowering. Identified markers should be useful for marker-assisted selection in cotton breeding to improve early flowering characteristics.

Development, genetic mapping and QTL association of cotton PHYA, PHYB, and HY5-specific CAPS and dCAPS markers

BACKGROUND

Among SNP markers that become increasingly valuable in molecular breeding of crop plants are the CAPS and dCAPS markers derived from the genes of interest. To date, the number of such gene-based markers is small in polyploid crop plants such as allotetraploid cotton that has A- and D-sub-genomes. The objective of this study was to develop and map new CAPS and dCAPS markers for cotton developmental-regulatory genes that are important in plant breeding programs.

RESULTS

Gossypium hirsutum and G. barbadense, are the two cultivated allotetraploid cotton species. These have distinct fiber quality and other agronomic traits. Using comparative sequence analysis of characterized GSTs of the PHYA1, PHYB, and HY5 genes of G. hirsutum and G. barbadense one PHYA1-specific Mbo I/Dpn II CAPS, one PHYB-specific Alu I dCAPS, and one HY5-specific Hinf I dCAPS cotton markers were developed. These markers have successfully differentiated the two allotetraploid genomes (AD₁ and AD₂) when tested in parental genotypes of 'Texas Marker-1' ('TM-1'), 'Pima 3-79' and their F₁ hybrids. The genetic mapping and chromosome substitution line-based deletion analyses revealed that PHYA1 gene is located in A-sub-genome chromosome 11, PHYB gene is in A-sub-genome chromosome 10, and HY5 gene is in D-sub-genome chromosome 24, on the reference 'TM-1' x 'Pima 3-79' RIL genetic map. Further, it was found that genetic linkage map regions containing phytochrome and HY5-specific markers were associated with major fiber quality and flowering time traits in previously published QTL mapping studies.

CONCLUSION

This study detailed the genome mapping of three cotton phytochrome genes with newly developed CAPS and dCAPS markers. The proximity of these loci to fiber quality and other cotton QTL was demonstrated in two A-subgenome and one D-subgenome chromosomes. These candidate gene markers will be valuable for marker-assisted selection (MAS) programs to rapidly introgress G. barbadense phytochromes and/or HY5 gene (s) into G. hirsutum cotton genotypes or vice versa.