To Be a Flower or Fruiting Branch: Insights Revealed by mRNA and Small RNA Transcriptomes from Different Cotton Developmental Stages

Quantitative Trait Locus Mapping for Plant Height and Branch Number in CCRI70 Recombinant Inbred Line Population of Upland Cotton (Gossypium hirsutum)

Upland cotton accounts for a high percentage (95%) of the world's cotton production. Plant height (PH) and branch number (BN) are two important agronomic traits that have an impact on improving the level of cotton mechanical harvesting and cotton yield. In this research, a recombinant inbred line (RIL) population with 250 lines developed from the variety CCRI70 was used for constructing a high-density genetic map and identification of quantitative trait locus (QTL). The results showed that the map harbored 8298 single nucleotide polymorphism (SNP) markers, spanning a total distance of 4876.70 centimorgans (cMs). A total of 69 QTLs for PH (9 stable) and 63 for BN (11 stable) were identified and only one for PH was reported in previous studies. The QTLs for PH and BN harbored 495 and 446 genes, respectively. Combining the annotation information, expression patterns and previous studies of these genes, six genes could be considered as potential candidate genes for PH and BN. The results could be helpful for cotton researchers to better understand the genetic mechanism of PH and BN development, as well as provide valuable genetic resources for cotton breeders to manipulate cotton plant architecture to meet future demands.

Regulation of endogenous hormone and miRNA in leaves of alfalfa (Medicago sativa L.) seedlings under drought stress by endogenous nitric oxide

BACKGROUND

Alfalfa (Medicago sativa. L) is one of the best leguminous herbage in China and even in the world, with high nutritional and ecological value. However, one of the drawbacks of alfalfa is its sensitivity to dry conditions, which is a global agricultural problem. The objective of this study was to investigate the regulatory effects of endogenous nitric oxide (NO) on endogenous hormones and related miRNAs in alfalfa seedling leaves under drought stress. The effects of endogenous NO on endogenous hormones such as ABA, GA3, SA, and IAA in alfalfa leaves under drought stress were studied. In addition, high-throughput sequencing technology was used to identify drought-related miRNAs and endogenous NO-responsive miRNAs in alfalfa seedling leaves under drought stress.

RESULT

By measuring the contents of four endogenous hormones in alfalfa leaves, it was found that endogenous NO could regulate plant growth and stress resistance by inducing the metabolism levels of IAA, ABA, GA3, and SA in alfalfa, especially ABA and SA in alfalfa. In addition, small RNA sequencing technology and bioinformatics methods were used to analyze endogenous NO-responsive miRNAs under drought stress. It was found that most miRNAs were enriched in biological pathways and molecular functions related to hormones (ABA, ETH, and JA), phenylpropane metabolism, and plant stress tolerance.

CONCLUSION

In this study, the analysis of endogenous hormone signals and miRNAs in alfalfa leaves under PEG and PEG + cPTIO conditions provided an important basis for endogenous NO to improve the drought resistance of alfalfa at the physiological and molecular levels. It has important scientific value and practical significance for endogenous NO to improve plant drought resistance.

Phylogenetic and functional analysis of tiller angle control homeologs in allotetraploid cotton

Introduction

Plants can adapt their growth to optimize light capture in competitive environments, with branch angle being a crucial factor influencing plant phenotype and physiology. Decreased branch angles in cereal crops have been shown to enhance productivity in high-density plantings. The Tiller Angle Control (TAC1) gene, known for regulating tiller inclination in rice and corn, has been found to control branch angle in eudicots. Manipulating TAC1 in field crops like cotton offers the potential for improving crop productivity.

Methods

Using a homolog-based methodology, we examined the distribution of TAC1-related genes in cotton compared to other angiosperms. Furthermore, tissue-specific qPCR analysis unveiled distinct expression patterns of TAC1 genes in various cotton tissues. To silence highly expressed specific TAC1 homeologs in the stem, we applied CRISPR-Cas9 gene editing and Agrobacterium-mediated transformation, followed by genotyping and subsequent phenotypic validation of the mutants.

Results

Gene duplication events of TAC1 specific to the Gossypium lineage were identified, with 3 copies in diploid progenitors and 6 copies in allotetraploid cottons. Sequence analysis of the TAC1 homeologs in Gossypium hirsutum revealed divergence from other angiosperms with 1-2 copies, suggesting possible neo- or sub-functionalization for the duplicated copies. These TAC1 homeologs exhibited distinct gene expression patterns in various tissues over developmental time, with elevated expression of A11G109300 and D11G112200, specifically in flowers and stems, respectively. CRISPR-mediated loss of these TAC1 homeologous genes resulted in a reduction in branch angle and altered petiole angles, and a 5 to 10-fold reduction in TAC1 expression in the mutants, confirming their role in controlling branch and petiole angles. This research provides a promising strategy for genetically engineering branch and petiole angles in commercial cotton varieties, potentially leading to increased productivity.

Comprehensive transcriptome and WGCNA analysis reveals the potential function of anthocyanins in low-temperature resistance of a red flower mutant tobacco

The anthocyanin is a protective substance in various plants, and plays important roles in resisting to low-temperature. Here, we explored transcriptome analysis of pink flower (as CK) and the natural mutant red flower (as research objects) under low-temperature conditions, and aimed to reveal the potential functions of anthocyanins and anthocyanin-related regulatory factors in resistance to low-temperature. Our results showed that most of the differentially expressed genes (DEGs) encoding key enzymes in the late stage of anthocyanin metabolism in the mutant were significantly up-regulated. Meanwhile, several genes significantly differentially expressed in CK or mutant were obtained by classification and analysis of transcription factors (TFs), phytohormones and osmoregulators. Additionally, WGCNA was carried out to mine hub genes resistanted to low-temperature stress in flavonoid pathway. Finally, one UFGT family gene, three MYB and one bHLH were obtained as the future hub genes of this study. Combined with the above information, we concluded that the ability of the red flower mutant to grow and develop normally at low-temperatures was the result of a combination of flavonoids and cold resistance genes.

COTTONOMICS: a comprehensive cotton multi-omics database

The rapid advancement of sequencing technology, including next-generation sequencing (NGS), has greatly improved sequencing efficiency and decreased cost. Consequently, huge amounts of genomic, transcriptomic and epigenetic data concerning cotton species have been generated and released. These large-scale data provide immense opportunities for the study of cotton genomic structure and evolution, population genetic diversity and genome-wide mining of excellent genes for important traits. However, the complexity of NGS data also causes distress, as it cannot be utilized easily. Here, we presented the cotton omics data platform COTTONOMICS (http://cotton.zju.edu.cn/), an easily accessible web database that integrates 32.5 TB of omics data including seven assembled genomes, resequencing data from 1180 allotetraploid cotton accessions and RNA-sequencing (RNA-seq), small RNA-sequencing (smRNA-seq), Chromatin Immunoprecipitation sequencing (ChIP-seq), DNase hypersensitive sites sequencing (DNase-seq) and Bisulfite sequencing (BS-seq). COTTONOMICS allows users to employ various search scenarios and retrieve information concerning the cotton genomes, genomic variation (Single nucleotide polymorphisms (SNPs) and Insertion and Deletion (InDels)), gene expression, smRNA expression, epigenetic regulation and quantitative trait locus (QTLs). The user-friendly web interface offers a variety of modules for storing, retrieving, analyzing and visualizing cotton multi-omics data to diverse ends, thereby enabling users to decipher cotton population genetics and identify potential novel genes that influence agronomically beneficial traits.

Database URL: http://cotton.zju.edu.cn

Insights Into MicroRNA-Mediated Regulation of Flowering Time in Cotton Through Small RNA Sequencing

The timing of flowering is a key determinant for plant reproductive. It has been demonstrated that microRNAs (miRNAs) play an important role in transition from the vegetative to reproductive stage in cotton; however, knowledge remains limited about the regulatory role of miRNAs involved in flowering time regulation in cotton. To elucidate the molecular basis of miRNAs in response to flowering time in cotton, we performed high-throughput small RNA sequencing at the fifth true leaf stage. We identified 56 and 43 miRNAs that were significantly up- and downregulated in two elite early flowering cultivars (EFC) compared with two late flowering cultivars (LFC), respectively. The miRNA targets by RNA sequencing analysis showed that GhSPL4 in SBP transcription factor family targeted by GhmiR156 was significantly upregulated in EFCs. Co-expression regulatory network analysis (WGCNA) revealed that GhSOC1, GhAP1, GhFD, GhCOL3, and GhAGL16 act as node genes in the auxin- and gibberellin-mediated flowering time regulatory networks in cotton. Therefore, elucidation of miRNA-mediated flowering time regulatory network will contribute to our understanding of molecular mechanisms underlying flowering time in cotton.

Genomic mapping and identification of candidate genes encoding nulliplex-branch trait in sea-island cotton (Gossypium barbadense L.) by multi-omics analysis

Nulliplex branch is a key architectural trait in sea-island cotton (Gossypium barbadense L.), but its genetic basis is not well understood. Here we investigated the genetic basis of the nulliplex-branch trait in cotton by combining newly created bulked segregant analysis (BSA)-seq data, published RNA-seq data, and published whole-genome resequencing (WGR) data. We delimited the nulliplex-branch locus (qD07-NB) to D07, region 14.8-17.1 Mb, using various BSA methods and markers. We integrated our BSA data with WGR data of sea-island cotton varieties and detected a missense single-nucleotide polymorphism in the candidate gene (Gbar_D07G011870) of qD07-NB. This gene was under positive selection during sea-island cotton breeding in the Xinjiang Uygur Autonomous Region, China. Notably, the nulliplex-branch varieties possessed a better fiber quality than the long-branch varieties, and a set of high-quality molecular markers was identified for molecular breeding of the nulliplex-branch trait in cotton. We combined BSA-seq and RNA-seq data to compare gene expression profiles between two elite sea-island cotton varieties during three developmental stages. We identified eleven relevant candidate genes, five downregulated and six upregulated, in the qD07-NB locus. This research will expand our understanding of the genetic basis of the nulliplex-branch trait and provide guidance for architecture-focused breeding in cotton.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-021-01229-w.

Biocidal radiuses of cycloxaprid, imidacloprid and lambda-cyhalothrin droplets controlling against cotton aphid (Aphis gossypii) using an unmanned aerial vehicle

BACKGROUND

Unmanned aerial vehicles (UAVs) are a recently advanced aerial spraying technology. However, the median lethal number of droplets (LN₅₀ ) and biocidal radiuses (r₅₀ ) of insecticides droplets sprayed by UAVs are still unknown. Therefore, we evaluated the LN₅₀ and r₅₀ of cycloxaprid, imidacloprid and lambda-cyhalothrin droplets associated with adjuvant controlling against Aphis gossypii.

RESULTS

A small UAV and Potter spray tower (PST) were used to generate different size of droplets (Dv0.5 = 185 ± 5 and 43 ± 2 μm). The mortality of A. gossypii showed a droplet density-dependent process. At the concentration of 5 g L^-1 , the LN₅₀ of lambda-cyhalothrin, imidacloprid and cycloxaprid droplets sprayed by the UAV were only 49.2, 34.6 and 19.7 droplets cm^-2 , respectively, and the r₅₀ were 0.57, 0.68 and 0.90 mm, respectively. The LN₅₀ values were negatively correlated with insecticide concentrations, but the r₅₀ values increased with the increasing concentrations. Although cycloxaprid is less toxic to A. gossypii than lambda-cyhalothrin, cycloxaprid had a larger r₅₀ than lambda-cyhalothrin due to its translocation ability. Furthermore, cycloxaprid had a relatively larger r₅₀ than imidacloprid because it is more toxic to A. gossypii. Moreover, adjuvant silwet DRS-60 can significantly increase the r₅₀ of droplets.

CONCLUSION

Our results revealed that the r₅₀ of a droplet was higher than its droplet size at tested concentrations. Smaller droplets generated by the PST had relatively higher insecticidal potential. The median lethal dose (LD₅₀ ) and translocation ability of insecticides and spray adjuvant contributed to their r₅₀ . Therefore, ultra-low-volume spray with UAVs is feasible to control cotton aphids. © 2020 Society of Chemical Industry.

Genome-wide sRNA and mRNA transcriptomic profiling insights into dynamic regulation of taproot thickening in radish (Raphanus sativus L.)

BACKGROUND

Taproot is the main edible organ and ultimately determines radish yield and quality. However, the precise molecular mechanism underlying taproot thickening awaits further investigation in radish. Here, RNA-seq was performed to identify critical genes involved in radish taproot thickening from three advanced inbred lines with different root size.

RESULTS

A total of 2606 differentially expressed genes (DEGs) were shared between 'NAU-DY' (large acicular) and 'NAU-YB' (medium obovate), which were significantly enriched in 'phenylpropanoid biosynthesis', 'glucosinolate biosynthesis', and 'starch and sucrose metabolism' pathway. Meanwhile, a total of 16 differentially expressed miRNAs (DEMs) were shared between 'NAU-DY' and 'NAU-YH' (small circular), whereas 12 miRNAs exhibited specific differential expression in 'NAU-DY'. Association analysis indicated that miR393a-bHLH77, miR167c-ARF8, and miR5658-APL might be key factors to biological phenomenon of taproot type variation, and a putative regulatory model of taproot thickening and development was proposed. Furthermore, several critical genes including SUS1, EXPB3, and CDC5 were characterized and profiled by RT-qPCR analysis.

CONCLUSION

This integrated study on the transcriptional and post-transcriptional profiles could provide new insights into comprehensive understanding of the molecular regulatory mechanism underlying taproot thickening in root vegetable crops.

Genome-wide analysis and comparison of the DNA-binding one zinc finger gene family in diploid and tetraploid cotton (Gossypium)

The Dof (DNA-binding one zinc finger) transcription factor family is a representative of plant-specific classes of transcription factors. In this study, we performed a genome-wide screening and characterization of the Dof gene family within two tetraploid species Gossypium barbadense, Gossypium hirsutum, and two diploid species Gossypium arboreum, Gossypium raimondii. 115, 116, 55 and 56 Dof genes were identified respectively and all of the genes contain a sequence encoding the Dof DNA-binding domain. Those genes were unevenly distributed across 13/26 chromosomes of the cotton. Genome comparison revealed that segmental duplication may have played crucial roles in the expansion of the cotton Dof gene family, and tandem duplication also played a minor role. Analysis of RNA-Seq data indicated that cotton Dof gene expression levels varied across different tissues and in response to different abiotic stress. Overall, our results could provide valuable information for better understanding the evolution of cotton Dof genes, and lays a foundation for future investigation in cotton.

Genome sequencing of the Australian wild diploid species Gossypium australe highlights disease resistance and delayed gland morphogenesis

The diploid wild cotton species Gossypium australe possesses excellent traits including resistance to disease and delayed gland morphogenesis, and has been successfully used for distant breeding programmes to incorporate disease resistance traits into domesticated cotton. Here, we sequenced the G. australe genome by integrating PacBio, Illumina short read, BioNano (DLS) and Hi-C technologies, and acquired a high-quality reference genome with a contig N50 of 1.83 Mb and a scaffold N50 of 143.60 Mb. We found that 73.5% of the G. australe genome is composed of various repeat sequences, differing from those of G. arboreum (85.39%), G. hirsutum (69.86%) and G. barbadense (69.83%). The G. australe genome showed closer collinear relationships with the genome of G. arboreum than G. raimondii and has undergone less extensive genome reorganization than the G. arboreum genome. Selection signature and transcriptomics analyses implicated multiple genes in disease resistance responses, including GauCCD7 and GauCBP1, and experiments revealed induction of both genes by Verticillium dahliae and by the plant hormones strigolactone (GR24), salicylic acid (SA) and methyl jasmonate (MeJA). Experiments using a Verticillium-resistant domesticated G. barbadense cultivar confirmed that knockdown of the homologues of these genes caused a significant reduction in resistance against Verticillium dahliae. Moreover, knockdown of a newly identified gland-associated gene GauGRAS1 caused a glandless phenotype in partial tissues using G. australe. The G. australe genome represents a valuable resource for cotton research and distant relative breeding as well as for understanding the evolutionary history of crop genomes.

Light- and temperature-entrainable circadian clock in soybean development

In plants, the spatiotemporal expression of circadian oscillators provides adaptive advantages in diverse species. However, the molecular basis of circadian clock in soybean is not known. In this study, we used soybean hairy roots expression system to monitor endogenous circadian rhythms and the sensitivity of circadian clock to environmental stimuli. We discovered in experiments with constant light and temperature conditions that the promoters of clock genes GmLCLb2 and GmPRR9b1 drive a self-sustained, robust oscillation of about 24-h in soybean hairy roots. Moreover, we demonstrate that circadian clock is entrainable by ambient light/dark or temperature cycles. Specifically, we show that light and cold temperature pulses can induce phase shifts of circadian rhythm, and we found that the magnitude and direction of phase responses depends on the specific time of these two zeitgeber stimuli. We obtained a quadruple mutant lacking the soybean gene GmLCLa1, LCLa2, LCLb1, and LCLb2 using CRISPR, and found that loss-of-function of these four GmLCL orthologs leads to an extreme short-period circadian rhythm and late-flowering phenotype in transgenic soybean. Our study establishes that the morning-phased GmLCLs genes act constitutively to maintain circadian rhythmicity and demonstrates that their absence delays the transition from vegetative growth to reproductive development.

Functional analysis of the GbDWARF14 gene associated with branching development in cotton

Plant architecture, including branching pattern, is an important agronomic trait of cotton crops. In recent years, strigolactones (SLs) have been considered important plant hormones that regulate branch development. In some species such as Arabidopsis, DWARF14 is an unconventional receptor that plays an important role in the SL signaling pathway. However, studies on SL receptors in cotton are still lacking. Here, we cloned and analysed the structure of the GbD14 gene in Gossypium barbadense and found that it contains the domains necessary for a SL receptor. The GbD14 gene was expressed primarily in the roots, leaves and vascular bundles, and the GbD14 protein was determined via GFP to localize to the cytoplasm and nucleus. Gene expression analysis revealed that the GbD14 gene not only responded to SL signals but also was differentially expressed between cotton plants whose types of branching differed. In particular, GbD14 was expressed mainly in the axillary buds of normal-branching cotton, while it was expressed the most in the leaves of nulliplex-branch cotton. In cotton, the GbD14 gene can be induced by SL and other plant hormones, such as indoleacetic acid, abscisic acid, and jasmonic acid. Compared with wild-type Arabidopsis, GbD14-overexpressing Arabidopsis responded more rapidly to SL signals. Moreover, we also found that GbD14 can rescue the multi-branched phenotype of Arabidopsis Atd14 mutants. Our results indicate that the function of GbD14 is similar to that of AtD14, and GbD14 may be a receptor for SL in cotton and involved in regulating branch development. This research provides a theoretical basis for a profound understanding of the molecular mechanism of branch development and ideal plant architecture for cotton breeding improvements.

MicroRNAome Profile of Euphorbia kansui in Response to Methyl Jasmonate

miRNAs play vital regulatory roles in different plant developmental stages and in plant response to biotic and abiotic stresses. However, information is limited on the miRNA regulatory mechanism to methyl jasmonate (MeJA). In this study, we used the microRNAome profile to illustrate the relevant regulatory mechanisms of Euphorbia kansui in response to methyl jasmonate (MeJA) through Illumina RNA-Seq. As a result, we identified 875 miRNAs corresponding to 11,277 target mRNAs, among them, 168 known miRNA families representing 6019 target mRNAs sequences were obtained. 452 miRNA-mRNA pairs presented an anti-correlationship (Cor < −0.50 and p-value of correlation ≤ 0.05). The miRNA with a fold change ≥ 2 and a p (p-Value) < 0.05 in pairwise comparison were identified as significant differentially expressed miRNAs (DEMs). The DEMs in MeJA treatment of 0, 24, 36 and 48 h were compared by using Short Time Expression Miner (STEM) cluster and 4 significant gene profiles (p-value ≤ 0.02) were identified. Through the kyoto encyclopedia of genes and genomes (KEGG) pathway and gene ontology (GO) enrichment analysis on all miRNA targets, we identified 33 mRNAs in terpenoid biosynthesis, which were regulated by miRNAs under MeJA treatment, so the miRNA maybe involved in the response of E. kansui plant to exogenous MeJA and the results would provide very useful information on illustrating the regulatory mechanism of E. kansui and also provide an overall view of the miRNAs response to MeJA stress of a non-model plant.

Nulliplex-branch, a TERMINAL FLOWER 1 ortholog, controls plant growth habit in cotton

Nulliplex-branch (nb) mutants in cotton display a specific architecture. The gene responsible for the nb phenotype was identified, and its modulation mode was further studied. Plant architecture is an important agronomic factor influencing various traits such as yield and variety adaptability in crop plants. Cotton (Gossypium) simultaneously displays monopodial and sympodial growth. Nulliplex-branch (nb) mutants showing determinate sympodial shoots have been reported in both G. hirsutum (Ghnb) and G. barbadense (Gbnb). In this study, the gene responsible for the nb phenotype was identified. GhNB and GbNB were found to be allelic loci and are TERMINAL FLOWER 1 orthologs on the Dt subgenome, though the At copies remain native. Sequencing and association analyses identified four (Gh-nb1-Gh-nb4) and one (Gb-nb1) type of point mutation in the coding sequences of Ghnb and Gbnb, respectively. The NB gene was mainly expressed in the root and shoot apex, and expression rhythms were also observed in these tissues, suggesting that the expression of the NB gene could be regulated by photoperiod. Constitutive overexpression of GhNB suppresses the differentiation of the reproductive shoots. Knockout of both copies of GhNB caused the main and lateral shoots to terminate in flowers, which is a more determinate architecture than that of the nb mutants and implies that its function might be dosage dependent. A protein lipid overlay assay indicated that the amino acid substitutions in Gh-nb1 and Gb-nb1 weaken the ligand-binding activity of the NB protein in vitro. These findings suggest that the NB gene plays crucial roles in regulating the determinacy of shoots, and the modulation of this gene should constitute an effective crop improvement approach through adjusting the growth habit of cotton.

Natural variation in a CENTRORADIALIS homolog contributed to cluster fruiting and early maturity in cotton

BACKGROUND

Plant architecture and the vegetative-reproductive transition have major impacts on the agronomic success of crop plants, but genetic mechanisms underlying these traits in cotton (Gossypium spp.) have not been identified.

RESULTS

We identify four natural mutations in GoCEN-D_t associated with cluster fruiting (cl) and early maturity. The situ hybridization shows that GhCEN is preferentially expressed in cotton shoot apical meristems (SAM) of the main stem and axillary buds. Constitutive GhCEN-Dt overexpression suppresses the transition of the cotton vegetative apex to a reproductive shoot. Silencing GoCEN leads to early flowering and determinate growth, and in tetraploids causes the main stem to terminate in a floral bud, a novel phenotype that exemplifies co-adaptation of polyploid subgenomes and suggests new research and/or crop improvement approaches. Natural cl variations are enriched in cottons adapted to high latitudes with short frost-free periods, indicating that mutants of GoCEN have been strongly selected for early maturity.

CONCLUSION

We show that the cotton gene GoCEN-Dt, a homolog of Antirrhinum CENTRORADIALIS, is responsible for determinate growth habit and cluster fruiting. Insight into the genetic control of branch and flower differentiation offers new approaches to develop early maturing cultivars of cotton and other crops with plant architecture appropriate for mechanical harvesting.

Changes in DNA methylation assessed by genomic bisulfite sequencing suggest a role for DNA methylation in cotton fruiting branch development

Cotton plant architecture, including fruit branch formation and flowering pattern, influences plant light exploitation, cotton yield and planting cost. DNA methylation has been widely observed at different developmental stages in both plants and animals and is associated with regulation of gene expression, chromatin remodelling, genome protection and other functions. Here, we investigated the global epigenetic reprogramming during the development of fruiting branches and floral buds at three developmental stages: the seedling stage, the pre-squaring stage and the squaring stage. We first identified 22 cotton genes which potentially encode DNA methyltransferases and demethylases. Among them, the homologous genes of CMT, DRM2 and MET1 were upregulated at pre-squaring and squaring stages, suggesting that DNA methylation is involved in the development of floral buds and fruit branches. Although the global methylation at all of three developmental stages was not changed, the CHG-type methylation of non-expressed genes was higher than those of expressed genes. In addition, we found that the expression of the homologous genes of the key circadian rhythm regulators, including CRY, LHY and CO, was associated with changes of DNA methylation at three developmental stages.