Comprehensive Analysis of the COBRA-Like (COBL) Gene Family in Gossypium Identifies Two COBLs Potentially Associated with Fiber Quality

A COBRA family protein, PtrCOB3, contributes to gelatinous layer formation of tension wood fibers in poplar

Angiosperm trees usually develop tension wood (TW) in response to gravitational stimulation. TW comprises abundant gelatinous (G-) fibers with thick G-layers primarily composed of crystalline cellulose. Understanding the pivotal factors governing G-layer formation in TW fiber remains elusive. This study elucidates the role of a Populus trichocarpa COBRA family protein, PtrCOB3, in the G-layer formation of TW fibers. PtrCOB3 expression was upregulated, and its promoter activity was enhanced during TW formation. Comparative analysis with wild-type trees revealed that ptrcob3 mutants, mediated by Cas9/gRNA gene editing, were incapable of producing G-layers within TW fibers and showed severely impaired stem lift. Fluorescence immunolabeling data revealed a dearth of crystalline cellulose in the tertiary cell wall (TCW) of ptrcob3 TW fibers. The role of PtrCOB3 in G-layer formation is contingent upon its native promoter, as evidenced by the comparative phenotypic assessments of pCOB11::PtrCOB3, pCOB3::PtrCOB3, and pCOB3::PtrCOB11 transgenic lines in the ptrcob3 background. Overexpression of PtrCOB3 under the control of its native promoter expedited G-layer formation within TW fibers. We further identified 3 transcription factors that bind to the PtrCOB3 promoter and positively regulate its transcriptional levels. Alongside the primary TCW synthesis genes, these findings enable the construction of a 2-layer transcriptional regulatory network for the G-layer formation of TW fibers. Overall, this study uncovers mechanistic insight into TW formation, whereby a specific COB protein executes the deposition of cellulose, and consequently, G-layer formation within TW fibers.

Co-Expression Network Analysis and Introgressive Gene Identification for Fiber Length and Strength Reveal Transcriptional Differences in 15 Cotton Chromosome Substitution Segment Lines and Their Upland and Sea Island Parents

Fiber length (FL) and strength (FS) are the core indicators for evaluating cotton fiber quality. The corresponding stages of fiber elongation and secondary wall thickening are of great significance in determining FL and FS formation, respectively. QTL mapping and high-throughput sequencing technology have been applied to dissect the molecular mechanism of fiber development. In this study, 15 cotton chromosome segment substitution lines (CSSLs) with significant differences in FL and FS, together with their recurrent parental Gossypium hirsutum line CCRI45 and donor parent G. barbadense line Hai1, were chosen to conduct RNA-seq on developing fiber samples at 10 days post anthesis (DPA) and 20 DPA. Differentially expressed genes (DEGs) were obtained via pairwise comparisons among all 24 samples (each one with three biological repeats). A total of 969 DEGs related to FL-high, 1285 DEGs to FS-high, and 997 DEGs to FQ-high were identified. The functional enrichment analyses of them indicated that the GO terms of cell wall structure and ROS, carbohydrate, and phenylpropanoid metabolism were significantly enriched, while the GO terms of glucose and polysaccharide biosynthesis, and brassinosteroid and glycosylphosphatidylinositol metabolism could make great contributions to FL and FS formation, respectively. Weighted gene co-expressed network analyses (WGCNA) were separately conducted for analyzing FL and FS traits, and their corresponding hub DEGs were screened in significantly correlated expression modules, such as EXPA8, XTH, and HMA in the fiber elongation and WRKY, TDT, and RAC-like 2 during secondary wall thickening. An integrated analysis of these hub DEGs with previous QTL identification results successfully identified a total of 33 candidate introgressive DEGs with non-synonymous mutations between the Gh and Gb species. A common DEG encoding receptor-like protein kinase 1 was reported to likely participate in fiber secondary cell thickening regulation by brassionsteroid signaling. Such valuable information was conducive to enlightening the developing mechanism of cotton fiber and also provided an abundant gene pool for further molecular breeding.

Genome-wide dissection of the genetic factors underlying food quality in boiled and pounded white Guinea yam

BACKGROUND

Food quality traits related to the genetics of yam influence the acceptability for its consumption. This study aimed at identifying genetic factors underlying sensory and textural quality attributes of boiled and pounded yam, the two dominant food products from white Guinea yam.

RESULTS

A genome-wide association study (GWAS) of a panel of 184 genotypes derived from five multi-parent crosses population was conducted. The panel was phenotyped for the qualities of boiled and pounded yam using sensory quality and instrument-based textural profile assays. The genotypes displayed significant variation for most of the attributes. Population differentiation and structure analysis using principal component analysis (PCA) and population structure-based Bayesian information criteria revealed the presence of four well-defined clusters. The GWAS results from a multi-random mixed linear model with kinship and PCA used as covariate identified 13 single-nucleotide polymorphic (SNP) markers significantly associated with the boiled and pounded yam food qualities. The associated SNP markers explained 7.51-13.04% of the total phenotypic variance with a limit of detection exceeding 4.

CONCLUSION

Regions on chromosomes 7 and 15 were found to be associated with boiled and pounded yam quality attributes from sensory and instrument-based assays. Gene annotation analysis for the regions of associated SNPs revealed co-localization of several known putative genes involved in glucose export, hydrolysis and glycerol metabolism. Our study is one of the first reports of genetic factors underlying the boiled and pounded yam food quality to pave the way for marker-assisted selection in white Guinea yam. © 2023 Society of Chemical Industry.

Integrated multi-omic approach reveals the effect of a Graminaceae-derived biostimulant and its lighter fraction on salt-stressed lettuce plants

Plant biostimulants are widely applied in agriculture for their ability to improve plant fitness. In the present work, the impact of Graminaceae-derived protein hydrolysate (P) and its lighter molecular fraction F3 (< 1 kDa) on lettuce plants, subjected to either no salt or high salt conditions, was investigated through the combination of metabolomics and transcriptomics. The results showed that both treatments significantly modulated the transcriptome and metabolome of plants under salinity stress, highlighting an induction of the hormonal response. Nevertheless, P and F3 also displayed several peculiarities. F3 specifically modulated the response to ethylene and MAPK signaling pathway, whereas P treatment induced a down-accumulation of secondary metabolites, albeit genes controlling the biosynthesis of osmoprotectants and antioxidants were up-regulated. Moreover, according with the auxin response modulation, P promoted cell wall biogenesis and plasticity in salt-stressed plants. Notably, our data also outlined an epigenetic control of gene expression induced by P treatment. Contrarily, experimental data are just partially in agreement when not stressed plants, treated with P or F3, were considered. Indeed, the reduced accumulation of secondary metabolites and the analyses of hormone pathways modulation would suggest a preferential allocation of resources towards growth, that is not coherent with the down-regulation of the photosynthetic machinery, the CO2 assimilation rate and leaves biomass. In conclusion, our data demonstrate that, although they might activate different mechanisms, both the P and F3 can result in similar benefits, as far as the accumulation of protective osmolytes and the enhanced tolerance to oxidative stress are concerned. Notably, the F3 fraction exhibits slightly greater growth promotion effects under high salt conditions. Most importantly, this research further corroborates that biostimulants' mode of action is dependent on plants' physiological status and their composition, underscoring the importance of investigating the bioactivity of the different molecular components to design tailored applications for the agricultural practice.

A comprehensive overview of cotton genomics, biotechnology and molecular biological studies

Cotton is an irreplaceable economic crop currently domesticated in the human world for its extremely elongated fiber cells specialized in seed epidermis, which makes it of high research and application value. To date, numerous research on cotton has navigated various aspects, from multi-genome assembly, genome editing, mechanism of fiber development, metabolite biosynthesis, and analysis to genetic breeding. Genomic and 3D genomic studies reveal the origin of cotton species and the spatiotemporal asymmetric chromatin structure in fibers. Mature multiple genome editing systems, such as CRISPR/Cas9, Cas12 (Cpf1) and cytidine base editing (CBE), have been widely used in the study of candidate genes affecting fiber development. Based on this, the cotton fiber cell development network has been preliminarily drawn. Among them, the MYB-bHLH-WDR (MBW) transcription factor complex and IAA and BR signaling pathway regulate the initiation; various plant hormones, including ethylene, mediated regulatory network and membrane protein overlap fine-regulate elongation. Multistage transcription factors targeting CesA 4, 7, and 8 specifically dominate the whole process of secondary cell wall thickening. And fluorescently labeled cytoskeletal proteins can observe real-time dynamic changes in fiber development. Furthermore, research on the synthesis of cotton secondary metabolite gossypol, resistance to diseases and insect pests, plant architecture regulation, and seed oil utilization are all conducive to finding more high-quality breeding-related genes and subsequently facilitating the cultivation of better cotton varieties. This review summarizes the paramount research achievements in cotton molecular biology over the last few decades from the above aspects, thereby enabling us to conduct a status review on the current studies of cotton and provide strong theoretical support for the future direction.

Recent genome resequencing paraded COBRA-Like gene family roles in abiotic stress and wood formation in Poplar

A cell wall determines the mechanical properties of a cell, serves as a barrier against plant stresses, and allows cell division and growth processes. The COBRA-Like (COBL) gene family encodes a putative glycosylphosphatidylinositol (GPI)-anchored protein that controls cellulose deposition and cell progression in plants by contributing to the microfibril orientation of a cell wall. Despite being studied in different plant species, there is a dearth of the comprehensive global analysis of COBL genes in poplar. Poplar is employed as a model woody plant to study abiotic stresses and biomass production in tree research. Improved genome resequencing has enabled the comprehensive exploration of the evolution and functional capacities of PtrCOBLs (Poplar COBRA-Like genes) in poplar. Phylogeny analysis has discerned and classified PtrCOBLs into two groups resembling the Arabidopsis COBL family, and group I genes possess longer proteins but have fewer exons than group II. Analysis of gene structure and motifs revealed PtrCOBLs maintained a rather stable motif and exon-intron pattern across members of the same group. Synteny and collinearity analyses exhibited that the evolution of the COBL gene family was heavily influenced by gene duplication events. PtrCOBL genes have undergone both segmental duplication and tandem duplication, followed by purifying selection. Promotor analysis flaunted various phytohormone-, growth- and stress-related cis-elements (e.g., MYB, ABA, MeJA, SA, AuxR, and ATBP1). Likewise, 29 Ptr-miRNAs of 20 families were found targeting 11 PtrCOBL genes. PtrCOBLs were found localized at the plasma membrane and extracellular matrix, while gene ontology analysis showed their involvement in plant development, plant growth, stress response, cellulose biosynthesis, and cell wall biogenesis. RNA-seq datasets depicted the bulk of PtrCOBL genes expression being found in plant stem tissues and leaves, rendering mechanical strength and rejoinders to environmental cues. PtrCOBL2, 3, 10, and 11 manifested the highest expression in vasculature and abiotic stress, and resemblant expression trends were upheld by qRT-PCR. Co-expression network analysis identified PtrCOBL2 and PtrCOBL3 as hub genes across all abiotic stresses and wood developing tissues. The current study reports regulating roles of PtrCOBLs in xylem differentiating tissues, tension wood formation, and abiotic stress latency that lay the groundwork for future functional studies of the PtrCOBL genes in poplar breeding.

Multi-omic characterization of the maize GPI synthesis mutant gwt1 with defects in kernel development

BACKGROUND

Glycosylphosphatidylinositol (GPI) and GPI-anchored proteins (GAPs) are important for cell wall formation and reproductive development in Arabidopsis. However, monocot counterparts that function in kernel endosperm development have yet to be discovered. Here, we performed a multi-omic analysis to explore the function of GPI related genes on kernel development in maize.

RESULTS

In maize, 48 counterparts of human GPI synthesis and lipid remodeling genes were identified, in which null mutation of the glucosaminyl-phosphatidylinositol O-acyltransferase1 gene, ZmGWT1, caused a kernel mutant (named gwt1) with defects in the basal endosperm transport layer (BETL). We performed plasma membrane (PM) proteomics to characterize the potential GAPs involved in kernel development. In total, 4,981 proteins were successfully identified in 10-DAP gwt1 kernels of mutant and wild-type (WT), including 1,638 membrane-anchored proteins with different posttranslational modifications. Forty-seven of the 256 predicted GAPs were differentially accumulated between gwt1 and WT. Two predicted BETL-specific GAPs (Zm00001d018837 and Zm00001d049834), which kept similar abundance at general proteome but with significantly decreased abundance at membrane proteome in gwt1 were highlighted.

CONCLUSIONS

Our results show the importance of GPI and GAPs for endosperm development and provide candidate genes for further investigation of the regulatory network in which ZmGWT1 participates.

Comparative transcriptome analysis reveals key pathways and genes involved in trichome development in tea plant (Camellia sinensis)

Trichomes, which develop from epidermal cells, are considered one of the important characteristics of the tea plant [Camellia sinensis (L.) O. Kuntze]. Many nutritional and metabolomic studies have indicated the important contributions of trichomes to tea products quality. However, understanding the regulation of trichome formation at the molecular level remains elusive in tea plants. Herein, we present a genome-wide comparative transcriptome analysis between the hairless Chuyeqi (CYQ) with fewer trichomes and the hairy Budiaomao (BDM) with more trichomes tea plant genotypes, toward the identification of biological processes and functional gene activities that occur during trichome development. In the present study, trichomes in both cultivars CYQ and BDM were unicellular, unbranched, straight, and soft-structured. The density of trichomes was the highest in the bud and tender leaf periods. Further, using the high-throughput sequencing method, we identified 48,856 unigenes, of which 31,574 were differentially expressed. In an analysis of 208 differentially expressed genes (DEGs) encoding transcription factors (TFs), five may involve in trichome development. In addition, on the basis of the Gene Ontology (GO) annotation and the weighted gene co-expression network analysis (WGCNA) results, we screened several DEGs that may contribute to trichome growth, including 66 DEGs related to plant resistance genes (PRGs), 172 DEGs related to cell wall biosynthesis pathway, 29 DEGs related to cell cycle pathway, and 45 DEGs related to cytoskeleton biosynthesis. Collectively, this study provided high-quality RNA-seq information to improve our understanding of the molecular regulatory mechanism of trichome development and lay a foundation for additional trichome studies in tea plants.

Comparative genomic analysis of the COBRA genes in six Rosaceae species and expression analysis in Chinese white pear (Pyrus bretschneideri)

COBRA-Like (COBL) genes encode a glycosylphosphatidylinositol (GPI) anchoring protein unique to plants. In current study, 87 COBRA genes were identified in 6 Rosaceae species, including Pyrus bretschneideri (16 genes), Malus domestica (22 genes), Fragaria vesca (13 genes), Prunus mume (11 genes), Rubus occidentalis (13 genes) and Prunus avium (12 genes). We revealed the evolution of the COBRA gene in six Rosaceae species by phylogeny, gene structure, conservative sequence, hydrophobicity analysis, gene replication events and sliding window analysis. In addition, based on the analysis of expression patterns in pear fruit combined with bioinformatics, we identified PbCOBL12 and PbCOBL13 as potential genes regulating secondary cell wall (SCW) formation during pear stone cell development. This study aimed to understand the evolutionary relationship of the COBRA gene in Rosaceae species, clarify the potential function of COBRA in pear fruit development, and provide essential theoretical basis and gene resources for improving pear fruit quality through genetical modification mechanism.

The SlSHN2 transcription factor contributes to cuticle formation and epidermal patterning in tomato fruit

Tomato (Solanum lycopersicum) is an established model for studying plant cuticle because of its thick cuticle covering and embedding the epidermal cells of the fruit. In this study, we screened an EMS mutant collection of the miniature tomato cultivar Micro-Tom for fruit cracking mutants and found a mutant displaying a glossy fruit phenotype. By using an established mapping-by-sequencing strategy, we identified the causal mutation in the SlSHN2 transcription factor that is specifically expressed in outer epidermis of growing fruit. The point mutation in the shn2 mutant introduces a K to N amino acid change in the highly conserved 'mm' domain of SHN proteins. The cuticle from shn2 fruit showed a ~ fivefold reduction in cutin while abundance and composition of waxes were barely affected. In addition to alterations in cuticle thickness and properties, epidermal patterning and polysaccharide composition of the cuticle were changed. RNAseq analysis further highlighted the altered expression of hundreds of genes in the fruit exocarp of shn2, including genes associated with cuticle and cell wall formation, hormone signaling and response, and transcriptional regulation. In conclusion, we showed that a point mutation in the transcriptional regulator SlSHN2 causes major changes in fruit cuticle formation and its coordination with epidermal patterning.

Genome-wide identification and expression profiling of the COBRA-like genes reveal likely roles in stem strength in rapeseed (Brassica napus L.)

The COBRA-like (COBL) genes play key roles in cell anisotropic expansion and the orientation of microfibrils. Mutations in these genes cause the brittle stem and induce pathogen responsive phenotypes in Arabidopsis and several crop plants. In this study, an in silico genome-wide analysis was performed to identify the COBL family members in Brassica. We identified 44, 20 and 23 COBL genes in B. napus and its diploid progenitor species B. rapa and B. oleracea, respectively. All the predicted COBL genes were phylogenetically clustered into two groups: the AtCOB group and the AtCOBL7 group. The conserved chromosome locations of COBLs in Arabidopsis and Brassica, together with clustering, indicated that the expansion of the COBL gene family in B. napus was primarily attributable to whole-genome triplication. Among the BnaCOBLs, 22 contained all the conserved motifs and derived from 9 of 12 subgroups. RNA-seq analysis was used to determine the tissue preferential expression patterns of various subgroups. BnaCOBL9, BnaCOBL35 and BnaCOBL41 were highly expressed in stem with high-breaking resistance, which implies these AtCOB subgroup members may be involved in stem development and stem breaking resistance of rapeseed. Our results of this study may help to elucidate the molecular properties of the COBRA gene family and provide informative clues for high stem-breaking resistance studies.

Characterization of new COBRA like (COBL) genes in wheat (Triticum aestivum) and their expression analysis under drought stress

BACKGROUND

The COBL genes encode a plant-specific glycosylphosphatidylinositol (GPI)-anchored protein. Recently identified COBRA genes are supposed as a key regulator of the orientation of cell expansion in the root indicating that COBRA gene family members are likely to be important players at the plasma membrane-cell wall interface.

METHODS AND RESULTS

Five COBL gene namely, TaCOBL 1, TaCOBL 2, TaCOBL 3, TaCOBL 4 and TaCOBL 5 were identified using database search and domain predictions. Chromosomal location of each gene was mapped on karyotype. Structure of genes, promoter analysis and phylogenetic analysis were performed using different bioinformatics tools. Set of novel SNPs were also predicted. Gene ontologies were analyzed, and the processes and pathways were identified in which COBRA genes were involved. The molecular weight all the cobra proteins was in range of 50-75 KDa with 429-461 amino acid residues. The COBL genes were mapped on homeologous groups 2, 4, 5, 6 and 7. Gene ontology analysis revealed that TaCOBL genes were involved in cellulose microfibril organization, mucilage biosynthetic process involved in seed coat development, plant-type cell wall biogenesis plant-type cell wall cellulose biosynthetic process, seed coat development and seed development. Three drought responsive cis-elements (WRKY, ABRE and DRE) were found nearby COBL genes The qRT-PCR revealed TaCOBL genes are drought responsive and can be further explored to understand their role in drought tolerance in wheat.

CONCLUSION

The comprehensive annotation and expression profiling of COBL genes revealed that all five COBL genes are drought response. The promoter cis-regulatory element analysis revealed that COBL genes had stress related WRKY, ABRE and DRE cis-regulatory elements. This evidence suggest that TaCOBL genes are involved in drought stress tolerance.

OsBC1L1 and OsBC1L8 function in stomatal development in rice

Stomata that are bordered by pairs of guard cells are specialized for regulating gas exchange and transpiration in plants. The stomatal morphology of grass is unique, characterized by two dumbbell-shaped guard cells flanked by two lateral subsidiary cells. This morphology and developmental pattern enable grass stomata to respond to environmental signals efficiently. In this study, we demonstrated that knockout either OsBC1L1 or OsBC1L8, two close homologs of OsBC1L family causes no discernible defects in rice stomatal development, however, the double knockout mutant osbc1l1 osbc1l8 exhibits excess stomatal production and stomatal clustering. OsBC1L1 overexpression also causes abnormal stomatal patterning in rice. Moreover, osbc1l1 osbc1l8 has many defective stomata complexes with only one subsidiary cell. The expression of OsSPCH2 and OsFAMA, two genes key to stomatal development is both down-regulated in osbc1l1 osbc1l8. In contrast, overexpressing OsBC1L1 suppresses only the expression of OsSPCH2. Both OsBC1L1 and OsBC1L8 could be detected to be localized at the cell plate and plasma membrane during cell division of guard mother cells and subsidiary mother cells. Taken together, these results suggest that OsBC1L1 and OsBC1L8 play essential roles in the development of rice stomatal complex likely through their involvement in cell reproduction.

Maize Brittle Stalk2-Like3, encoding a COBRA protein, functions in cell wall formation and carbohydrate partitioning

Carbohydrate partitioning from leaves to sink tissues is essential for plant growth and development. The maize (Zea mays) recessive carbohydrate partitioning defective28 (cpd28) and cpd47 mutants exhibit leaf chlorosis and accumulation of starch and soluble sugars. Transport studies with 14C-sucrose (Suc) found drastically decreased export from mature leaves in cpd28 and cpd47 mutants relative to wild-type siblings. Consistent with decreased Suc export, cpd28 mutants exhibited decreased phloem pressure in mature leaves, and altered phloem cell wall ultrastructure in immature and mature leaves. We identified the causative mutations in the Brittle Stalk2-Like3 (Bk2L3) gene, a member of the COBRA family, which is involved in cell wall development across angiosperms. None of the previously characterized COBRA genes are reported to affect carbohydrate export. Consistent with other characterized COBRA members, the BK2L3 protein localized to the plasma membrane, and the mutants condition a dwarf phenotype in dark-grown shoots and primary roots, as well as the loss of anisotropic cell elongation in the root elongation zone. Likewise, both mutants exhibit a significant cellulose deficiency in mature leaves. Therefore, Bk2L3 functions in tissue growth and cell wall development, and this work elucidates a unique connection between cellulose deposition in the phloem and whole-plant carbohydrate partitioning.

Anther dehiscence is regulated by gibberellic acid in yellow lupine (Lupinus luteus L.)

BACKGROUND

Anther dehiscence resulting in the release of pollen grains is tightly regulated in a spatiotemporal manner by various factors. In yellow lupine (Lupinus luteus L.), a species that shows cleistogamy, the anthers split before the flowers open, but the course and regulation of this process are unknown. The specific control of anther development takes place via hormonal pathways, the wide action of which ensures reproductive success. In our previous research concerning flower and early pod development in yellow lupine, we showed that the lowest transcript level of LlDELLA1, a main repressor of gibberellin (GA) signalling, occurs approximately at the time of anther opening; therefore, the main purpose of this study was to precisely investigate the gibberellic acid (GA3)-dependent regulation of the anther dehiscence in this species.

RESULTS

In this paper, we showed the specific changes in the yellow lupine anther structure during dehiscence, including secondary thickening in the endothecium by lignocellulosic deposition, enzymatic cell wall breakdown at the septum/stomium and cell degeneration via programmed cell death (PCD), and identified several genes widely associated with this process. The expression profile of genes varied over time, with the most intense mRNA accumulation in the phases prior to or at the time of anther opening. The transcriptional activity also revealed that these genes are highly coexpressed and regulated in a GA-dependent manner. The cellular and tissue localization of GA3 showed that these molecules are present before anther opening, mainly in septum cells, near the vascular bundle and in the endothecium, and that they are subsequently undetectable. GA3 localization strongly correlates with the transcriptional activity of genes related to GA biosynthesis and deactivation. The results also suggest that GA3 controls LlGAMYB expression via an LlMIR159-dependent pathway.

CONCLUSIONS

The presented results show a clear contribution of GA3 in the control of the extensive anther dehiscence process in yellow lupine. Understanding the processes underlying pollen release at the hormonal and molecular levels is a significant aspect of controlling fertility in this economically important legume crop species and is of increasing interest to breeders.

Genome-wide identification and expression analysis of the plant specific LIM genes in Gossypium arboreum under phytohormone, salt and pathogen stress

Asiatic cotton (Gossypium arboreum) cultivated as ‘desi cotton’ in India, is renowned for its climate resilience and robustness against biotic and abiotic stresses. The genome of G. arboreum is therefore, considered as a valued reserve of information for discovering novel genes or gene functions for trait improvements in the present context of cotton cultivation world-wide. In the present study, we carried out genome-wide analysis of LIM gene family in desi cotton and identified twenty LIM domain proteins (GaLIMs) which include sixteen animals CRP-like GaLIMs and four plant specific GaLIMs with presence (GaDA1) or absence (GaDAR) of UIM (Ubiquitin Interacting Motifs). Among the sixteen CRP-like GaLIMs, eleven had two conventional LIM domains while, five had single LIM domain which was not reported in LIM gene family of the plant species studied, except in Brassica rapa. Phylogenetic analysis of these twenty GaLIM proteins in comparison with LIMs of Arabidopsis, chickpea and poplar categorized them into distinct αLIM1, βLIM1, γLIM2, δLIM2 groups in CRP-like LIMs, and GaDA1 and GaDAR in plant specific LIMs group. Domain analysis had revealed consensus [(C-X₂-C-X₁₇-H-X₂-C)-X₂-(C-X₂-C-X₁₇-C-X₂-H)] and [(C-X₂-C-X₁₇-H-X₂-C)-X₂-(C-X₄-C-X₁₅-C-X₂-H)] being conserved as first and/or second LIM domains of animal CRP-like GaLIMs, respectively. Interestingly, single LIM domain containing GaLIM15 was found to contain unique consensus with longer inter-zinc-motif spacer but shorter second zinc finger motif. All twenty GaLIMs showed variable spatio-temporal expression patterns and accordingly further categorized into distinct groups of αLIM1, βLIM1, γLIM2 δLIM2 and plant specific LIM (DA1/DAR). For the first time, response of GaDA1/DAR under the influence of biotic and abiotic stresses were studied in cotton, involving treatments with phytohormones (Jasmonic acid and Abscisic acid), salt (NaCl) and wilt causing pathogen (Fusarium oxysporum). Expressions patterns of GaDA1/DAR showed variable response and identified GaDA2 as a probable candidate gene for stress tolerance in G. arboreum.

Genome-Wide Analysis of the COBRA-Like Gene Family Supports Gene Expansion through Whole-Genome Duplication in Soybean (Glycine max)

The COBRA-like (COBL) gene family has been associated with the regulation of cell wall expansion and cellulose deposition. COBL mutants result in reduced levels and disorganized deposition of cellulose causing defects in the cell wall and inhibiting plant development. In this study, we report the identification of 24 COBL genes (GmCOBL) in the soybean genome. Phylogenetic analysis revealed that the COBL proteins are divided into two groups, which differ by about 170 amino acids in the N-terminal region. The GmCOBL genes were heterogeneously distributed in 14 of the 20 soybean chromosomes. This study showed that segmental duplication has contributed significantly to the expansion of the COBL family in soybean during all Glycine-specific whole-genome duplication events. The expression profile revealed that the expression of the paralogous genes is highly variable between organs and tissues of the plant. Only 20% of the paralogous gene pairs showed similar expression patterns. The high expression levels of some GmCOBLs suggest they are likely essential for regulating cell expansion during the whole soybean life cycle. Our comprehensive overview of the COBL gene family in soybean provides useful information for further understanding the evolution and diversification of COBL genes in soybean.

Auxin Directly Upregulates GhRAC13 Expression to Promote the Onset of Secondary Cell Wall Deposition in Cotton Fibers

Cotton fibers are single cells that show a relatively independent developmental process of cell differentiation, elongation, and secondary wall deposition. Auxin promotes fiber cell protrusion from the surface of the ovule. However, the role of auxin at other stages of cotton fiber development remains largely unknown. To gain a deeper insight into this aspect, we measured indoleacetic acid (IAA) content in developing fibers. Results showed an increase in IAA content at the transition stage from elongation to secondary cell wall deposition. Subsequently, we investigated the differences between two transgenic cottons that show upregulated and downregulated fiber auxin levels, respectively. In planta analysis revealed that, in addition to promoting cell elongation, auxin regulated the time of initiation of reactive oxygen species (ROS) production and secondary wall deposition in cotton fibers. This was closely correlated with the upregulated expression of GhRAC13, which regulates ROS-triggered cellulose synthesis. We found multiple putative auxin-responsive elements existed within the promoter region of GhRAC13, and IAA could induce proGhRAC13 activity. The dual-luciferase reporter assay further proved the activation of proGhRAC13 by GhARF5, an auxin-signaling activator. Altogether, our results suggest a role of auxin in promoting the onset of secondary growth by directly upregulating GhRAC13 expression in cotton fibers.

Multi-omics analysis of cellular pathways involved in different rapid growth stages of moso bamboo

Moso bamboo (Phyllostachys edulis (Carriere) J. Houzeau) is a rapidly growing grass of industrial and ecological importance. However, the molecular mechanisms of its remarkable growth are not well understood. In this study, we investigated the early-stage growth of moso bamboo shoots and defined three different growth stages based on histological and biochemical analyses, namely, starting of cell division (SD), rapid division (RD) and rapid elongation (RE). Further analyses on potentially relevant cellular pathways in these growth stages using multi-omics approaches such as transcriptomics and proteomics revealed the involvement of multiple cellular pathways, including DNA replication, repair and ribosome biogenesis. A total of 8045 differentially expressed genes (DEGs) and 1053 differentially expressed proteins (DEPs) were identified in our analyses. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses of detected DEGs identified several key biological pathways such as phytohormone metabolism, signal transduction, cell wall development and carbohydrate metabolism. The comparative analysis of proteins displayed that a total of 213 DEPs corresponded with DEGs and 3 significant expression profiles that could be promoting the fast growth of bamboo internodes. Moreover, protein-protein interaction network prediction analysis is suggestive of the involvement of five major proteins of signal transduction, DNA synthesis and RNA transcription, and may act as key elements responsible for the rapid shoot growth. Our work exploits multi-omics and bioinformatic approaches to unfurl the complexity of molecular networks involved in the rapid growth of moso bamboo and opens up questions related to the interactions between the functions played by individual molecular pathway.

Comprehensive transcriptomics, proteomics, and metabolomics analyses of the mechanisms regulating tiller production in low-tillering wheat

Tiller development in low-tillering wheat is related to several differentially expressed genes, proteins, and metabolites, as determined by an integrated omics approach combining transcriptome analysis, iTRAQ, and HPLC-MS on multiple NILs. Tillering is an important aspect of plant morphology that affects spike number, thereby contributing to the final crop yield. However, the mechanisms inhibiting tiller production in low-tillering wheat are poorly characterized. To investigate this aspect of wheat biology, two pairs of near-isogenic lines were developed, and an integrated omics approach combining transcriptome analysis, isobaric tags for relative and absolute quantification, and high-performance liquid chromatography-mass spectrometry were used to compare the free-tillering and low-tillering caused by an allele at Qltn.sicau-2D in wheat samples. Overall, 474 genes, 166 proteins, and 28 metabolites were identified as tillering-associated differentially expressed genes, proteins, and metabolites (DEGs, DEPs, and DEMs, respectively). Functional analysis indicated that the abundance of DEGs/DEPs/DEMs was related to lignin and cellulose metabolism, cell division, cell cycle processes, and glycerophospholipid metabolism; three transcription factor families, GRAS, GRF, and REV, might be related to the decrease in tillering in low-tillering wheat. These findings contribute to improve our understanding of the mechanisms responsible for the inhibition of tiller development in low-tillering wheat cultivars.

Transcriptome analysis identified aberrant gene expression in pollen developmental pathways leading to CGMS in cotton (Gossypium hirsutum L.)

Male sterility (induced or natural) is a potential tool for commercial hybrid seed production in different crops. Despite numerous endeavors to understand the physiological, hereditary, and molecular cascade of events governing CMS in cotton, the exact biological process controlling sterility and fertility reconstruction remains obscure. During current study, RNA-Seq using Ion Torrent S5 platform is carried out to identify 'molecular portraits' in floral buds among the Cytoplasmic Genic Male Sterility (CGMS) line, its near-isogenic maintainer, and restorer lines. A total of 300, 438 and 455 genes were differentially expressed in CGMS, Maintainer, and Restorer lines respectively. The functional analysis using AgriGo revealed suppression in the pathways involved in biogenesis and metabolism of secondary metabolites which play an important role in pollen and anther maturation. Enrichment analysis showed dearth related to pollen and anther's development in sterile line, including anomalous expression of genes and transcription factors that have a role in the development of the reproductive organ, abnormal cytoskeleton formation, defects in cell wall formation. The current study found aberrant expression of DYT1, AMS and cytochrome P450 genes involved in tapetum formation, pollen development, pollen exine and anther cuticle formation associated to male sterility as well as fertility restoration of CGMS. In the current study, more numbers of DEGs were found on Chromosome D05 and A05 as compared to other chromosomes. Expression pattern analysis of fourteen randomly selected genes using qRT-PCR showed high concurrence with gene expression profile of RNA-Seq analysis accompanied by a strong correlation of 0.82. The present study provides an important support for future studies in identifying interaction between cyto-nuclear molecular portraits, to accelerate functional genomics and molecular breeding related to cytoplasmic male sterility studies in cotton.

Plant Cell Wall Proteomics: A Focus on Monocot Species, Brachypodium distachyon, Saccharum spp. and Oryza sativa

Plant cell walls mostly comprise polysaccharides and proteins. The composition of monocots' primary cell walls differs from that of dicots walls with respect to the type of hemicelluloses, the reduction of pectin abundance and the presence of aromatic molecules. Cell wall proteins (CWPs) differ among plant species, and their distribution within functional classes varies according to cell types, organs, developmental stages and/or environmental conditions. In this review, we go deeper into the findings of cell wall proteomics in monocot species and make a comparative analysis of the CWPs identified, considering their predicted functions, the organs analyzed, the plant developmental stage and their possible use as targets for biofuel production. Arabidopsis thaliana CWPs were considered as a reference to allow comparisons among different monocots, i.e., Brachypodium distachyon, Saccharum spp. and Oryza sativa. Altogether, 1159 CWPs have been acknowledged, and specificities and similarities are discussed. In particular, a search for A. thaliana homologs of CWPs identified so far in monocots allows the definition of monocot CWPs characteristics. Finally, the analysis of monocot CWPs appears to be a powerful tool for identifying candidate proteins of interest for tailoring cell walls to increase biomass yield of transformation for second-generation biofuels production.

Defining the Genetic Basis of Plant⁻Endophytic Bacteria Interactions

Endophytic bacteria, which interact closely with their host, are an essential part of the plant microbiome. These interactions enhance plant tolerance to environmental changes as well as promote plant growth, thus they have become attractive targets for increasing crop production. Numerous studies have aimed to characterise how endophytic bacteria infect and colonise their hosts as well as conferring important traits to the plant. In this review, we summarise the current knowledge regarding endophytic colonisation and focus on the insights that have been obtained from the mutants of bacteria and plants as well as ‘omic analyses. These show how endophytic bacteria produce various molecules and have a range of activities related to chemotaxis, motility, adhesion, bacterial cell wall properties, secretion, regulating transcription and utilising a substrate in order to establish a successful interaction. Colonisation is mediated by plant receptors and is regulated by the signalling that is connected with phytohormones such as auxin and jasmonic (JA) and salicylic acids (SA). We also highlight changes in the expression of small RNAs and modifications of the cell wall properties. Moreover, in order to exploit the beneficial plant-endophytic bacteria interactions in agriculture successfully, we show that the key aspects that govern successful interactions remain to be defined.

Transcriptomic analysis of cultivated cotton Gossypium hirsutum provides insights into host responses upon whitefly-mediated transmission of cotton leaf curl disease

Cotton is a commercial and economically important crop that generates billions of dollars in annual revenue worldwide. However, cotton yield is affected by a sap-sucking insect Bemisia tabaci (whitefly), and whitefly-borne cotton leaf curl disease (CLCuD). The causative agent of devastating CLCuD is led by the viruses belonging to the genus Begomovirus (family Geminiviridae), collectively called cotton leaf curl viruses. Unfortunately, the extensively cultivated cotton (Gossypium hirsutum) species are highly susceptible and vulnerable to CLCuD. Yet, the concomitant influence of whitefly and CLCuD on the susceptible G. hirsutum transcriptome has not been interpreted. In the present study we have employed an RNA Sequencing (RNA-Seq) transcriptomics approach to explore the differential gene expression in susceptible G. hirsutum variety upon infection with viruliferous whiteflies. Comparative RNA-Seq of control and CLCuD infected plants was done using Illumina HiSeq 2500. This study yielded 468 differentially expressed genes (DEGs). Among them, we identified 220 up and 248 downregulated DEGs involved in disease responses and pathogen defense. We selected ten genes for downstream RT-qPCR analyses on two cultivars, Karishma and MNH 786 that are susceptible to CLCuD. We observed a similar expression pattern of these genes in both susceptible cultivars that was also consistent with our transcriptome data further implying a wider application of our global transcription study on host susceptibility to CLCuD. We next performed weighted gene co-expression network analysis that revealed six modules. This analysis also identified highly co-expressed genes as well as 55 hub genes that co-express with ≥ 50 genes. Intriguingly, most of these hub genes are shown to be downregulated and enriched in cellular processes. Under-expression of such highly co-expressed genes suggests their roles in favoring the virus and enhancing plant susceptibility to CLCuD. We also discuss the potential mechanisms governing the establishment of disease susceptibility. Overall, our study provides a comprehensive differential gene expression analysis of G. hirsutum under whitefly-mediated CLCuD infection. This vital study will advance the understanding of simultaneous effect of whitefly and virus on their host and aid in identifying important G. hirsutum genes which intricate in its susceptibility to CLCuD.

Glycosylphosphatidylinositol-Anchored Proteins in Arabidopsis and One of Their Common Roles in Signaling Transduction

Diverse proteins are found modified with glycosylphosphatidylinositol (GPI) at their carboxyl terminus in eukaryotes, which allows them to associate with membrane lipid bilayers and anchor on the external surface of the plasma membrane. GPI-anchored proteins (GPI-APs) play crucial roles in various processes, and more and more GPI-APs have been identified and studied. In this review, previous genomic and proteomic predictions of GPI-APs in Arabidopsis have been updated, which reveal their high abundance and complexity. From studies of individual GPI-APs in Arabidopsis, certain GPI-APs have been found associated with partner receptor-like kinases (RLKs), targeting RLKs to their subcellular localization and helping to recognize extracellular signaling polypeptide ligands. Interestingly, the association might also be involved in ligand selection. The analyses suggest that GPI-APs are essential and widely involved in signal transduction through association with RLKs.

Plant glycosylphosphatidylinositol anchored proteins at the plasma membrane-cell wall nexus

Approximately 1% of plant proteins are predicted to be post-translationally modified with a glycosylphosphatidylinositol (GPI) anchor that tethers the polypeptide to the outer leaflet of the plasma membrane. Whereas the synthesis and structure of GPI anchors is largely conserved across eukaryotes, the repertoire of functional domains present in the GPI-anchored proteome has diverged substantially. In plants, this includes a large fraction of the GPI-anchored proteome being further modified with plant-specific arabinogalactan (AG) O-glycans. The importance of the GPI-anchored proteome to plant development is underscored by the fact that GPI biosynthetic null mutants exhibit embryo lethality. Mutations in genes encoding specific GPI-anchored proteins (GAPs) further supports their contribution to diverse biological processes, occurring at the interface of the plasma membrane and cell wall, including signaling, cell wall metabolism, cell wall polymer cross-linking, and plasmodesmatal transport. Here, we review the literature concerning plant GPI-anchored proteins, in the context of their potential to act as molecular hubs that mediate interactions between the plasma membrane and the cell wall, and their potential to transduce the signal into the protoplast and, thereby, activate signal transduction pathways.

Ectopic expression of GhCOBL9A, a cotton glycosyl-phosphatidyl inositol-anchored protein encoding gene, promotes cell elongation, thickening and increased plant biomass in transgenic Arabidopsis

Cellulose is a major component of plant cell walls and is necessary for plant morphogenesis and biomass. COBL (COBRA-Like) proteins have been shown to be key regulators in the orientation of cell expansion and cellulose crystallinity status. To clarify the role of a cotton COBL gene, GhCOBL9A, we conducted the ectopic expression and functional analysis in Arabidopsis. Previous study showed that GhCOBL9A was preferentially expressed during secondary cell wall biosynthesis in cotton fibers, and showed a significant co-expression pattern with cellulose synthase genes. Here, we detected that overexpression of GhCOBL9A induced the up-regulation of genes related to cellulose synthesis and enhanced the cellulose deposition. As a result, GhCOBL9A transgenic plants displayed increased hypocotyl and root lengths in early development, and cell wall thickening at the SCW stage. Notably, overexpression of GhCOBL9A led to an erect, robust-stature phenotype and brought higher biomass in mature plants. In addition, overexpression of GhCOBL9A in Arabidopsis AtCOBL4 mutants, a paralogous gene of GhCOBL9A, also led to a stronger growth potential, but the Atcobl4 mutant phenotype could not be rescued, implying the functional divergence of GhCOBL9A and AtCOBL4 paralogs. Taken together, these results suggest that overexpression of GhCOBL9A contributes to plant cell elongation and thickening, and increased biomass, which provides references for further utilizing GhCOBL9A to improve yield and quality traits in cotton and other species.

The phosphatidylinositol synthase gene (GhPIS) contributes to longer, stronger, and finer fibers in cotton

Cotton fibers are the most important natural raw material used in textile industries world-wide. Fiber length, strength, and fineness are the three major traits which determine the quality and economic value of cotton. It is known that exogenous application of phosphatidylinositols (PtdIns), important structural phospholipids, can promote cotton fiber elongation. Here, we sought to increase the in planta production of PtdIns to improve fiber traits. Transgenic cotton plants were generated in which the expression of a cotton phosphatidylinositol synthase gene (i.e., GhPIS) was controlled by the fiber-specific SCFP promoter element, resulting in the specific up-regulation of GhPIS during cotton fiber development. We demonstrate that PtdIns content was significantly enhanced in transgenic cotton fibers and the elevated level of PtdIns stimulated the expression of genes involved in PtdIns phosphorylation as well as promoting lignin/lignin-like phenolic biosynthesis. Fiber length, strength and fineness were also improved in the transgenic plants as compared to the wild-type cotton, with no loss in overall fiber yield. Our data indicate that fiber-specific up-regulation of PtdIns synthesis is a promising strategy for cotton fiber quality improvement.

Comparative transcriptome analysis of cotton fiber development of Upland cotton (Gossypium hirsutum) and Chromosome Segment Substitution Lines from G. hirsutum × G. barbadense

Background

How to develop new cotton varieties possessing high yield traits of Upland cotton and superior fiber quality traits of Sea Island cotton remains a key task for cotton breeders and researchers. While multiple attempts bring in little significant progresses, the development of Chromosome Segment Substitution Lines (CSSLs) from Gossypium barbadense in G. hirsutum background provided ideal materials for aforementioned breeding purposes in upland cotton improvement. Based on the excellent fiber performance and relatively clear chromosome substitution segments information identified by Simple Sequence Repeat (SSR) markers, two CSSLs, MBI9915 and MBI9749, together with the recurrent parent CCRI36 were chosen to conduct transcriptome sequencing during the development stages of fiber elongation and Secondary Cell Wall (SCW) synthesis (from 10DPA and 28DPA), aiming at revealing the mechanism of fiber development and the potential contribution of chromosome substitution segments from Sea Island cotton to fiber development of Upland cotton.

Results

In total, 15 RNA-seq libraries were constructed and sequenced separately, generating 705.433 million clean reads with mean GC content of 45.13% and average Q30 of 90.26%. Through multiple comparisons between libraries, 1801 differentially expressed genes (DEGs) were identified, of which the 902 up-regulated DEGs were mainly involved in cell wall organization and response to oxidative stress and auxin, while the 898 down-regulated ones participated in translation, regulation of transcription, DNA-templated and cytoplasmic translation based on GO annotation and KEGG enrichment analysis. Subsequently, STEM software was performed to explicate the temporal expression pattern of DEGs. Two peroxidases and four flavonoid pathway-related genes were identified in the “oxidation-reduction process”, which could play a role in fiber development and quality formation. Finally, the reliability of RNA-seq data was validated by quantitative real-time PCR of randomly selected 20 genes.

Conclusions

The present report focuses on the similarities and differences of transcriptome profiles between the two CSSLs and the recurrent parent CCRI36 and provides novel insights into the molecular mechanism of fiber development, and into further exploration of the feasible contribution of G. barbadense substitution segments to fiber quality formation, which will lay solid foundation for simultaneously improving fiber yield and quality of upland cotton through CSSLs.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4077-8) contains supplementary material, which is available to authorized users.

Molecular characterization of GhPLDα1 and its relationship with secondary cell wall thickening in cotton fibers

Phospholipase D (PLD) hydrolyzes phospholipids to generate a free polar head group (e.g., choline) and a second messenger phosphatidic acid and plays diverse roles in plant growth and development, including seed germination, leaf senescence, root hair growth, and hypocotyl elongation. However, the function of PLD in cotton remains largely unexplored. Here, the comprehensive molecular characterization of GhPLDα1 was explored with its role in upland cotton (Gossypium hirsutum) fiber development. The GhPLDα1 gene was cloned successfully, and a sequence alignment showed that GhPLDα1 contains one C2 domain and two HKD (HxKxxxxD) domains. Quantitative reverse transcriptase-polymerase chain reaction measured the expression of GhPLDα1 in various cotton tissues with the highest level in fibers at 20 days post anthesis (d.p.a.). Fluorescent microscopy and immunoblotting in tobacco epidermis showed the GhPLDα1 distribution in both cell membranes and the cytoplasm. An activity assay indicated changes in PLDα enzyme activity in developing fiber cells with a peak level at 20 d.p.a., coinciding with the onset of cellulose accumulation and the increased H₂O₂ content during fiber development. Furthermore, the inhibition of PLDα activity obviously decreased the cellulose and H₂O₂ contents of in vitro-cultured cotton fibers. These results provide important evidence explaining the relationship of GhPLDα1 with secondary cell wall thickening in cotton fibers in that GhPLDα1 may correlate with the increased H₂O₂ content at the onset of secondary cell wall thickening, ultimately promoting cellulose biosynthesis.

A Cotton Annexin Affects Fiber Elongation and Secondary Cell Wall Biosynthesis Associated with Ca2+ Influx, ROS Homeostasis, and Actin Filament Reorganization

Annexins play pivotal roles in a variety of cellular processes as well as in fiber development; however, the functional mechanisms of their activities are unclear. Here, an annexin gene that is preferentially expressed in fibers, GhFAnnxA, was found to be significantly associated with various cotton (Gossypium hirsutum) fiber traits. Transgenic analysis demonstrated that GhFAnnxA affected cotton fiber elongation and was involved in secondary cell wall (SCW) biosynthesis. Functional studies demonstrated that GhFAnnxA may act as a Ca(2+) conductance regulator and that reactive oxygen species (ROS) produced by Rbohs in a Ca(2+)-dependent manner may determine fiber elongation caused by elevated intracellular turgor and cell wall loosening. However, excessive hydrogen peroxide (H2O2) inhibited cotton fiber elongation in vitro. We speculate that a positive feedback loop involving ROS and Ca(2+) is regulated by GhCDPK1 and regulates fiber cell elongation. Furthermore, the convergence of actin filaments is altered by their interaction with GhFAnnxA, and this also may contribute to fiber elongation. Moreover, GhFAnnxA may affect SCW biosynthesis through changes in cell wall components caused by an increase in H2O2 levels. These results not only provide new insights into the signaling pathways of GhFAnnxA in fiber development but also clarify the role of ROS in fiber development.

The Glycerol-3-Phosphate Acyltransferase GPAT6 from Tomato Plays a Central Role in Fruit Cutin Biosynthesis

The thick cuticle covering and embedding the epidermal cells of tomato (Solanum lycopersicum) fruit acts not only as a protective barrier against pathogens and water loss but also influences quality traits such as brightness and postharvest shelf-life. In a recent study, we screened a mutant collection of the miniature tomato cultivar Micro-Tom and isolated several glossy fruit mutants in which the abundance of cutin, the polyester component of the cuticle, was strongly reduced. We employed a newly developed mapping-by-sequencing strategy to identify the causal mutation underlying the cutin deficiency in a mutant thereafter named gpat6-a (for glycerol-3-phosphate acyltransferase6). To this end, a backcross population (BC1F2) segregating for the glossy trait was phenotyped. Individuals displaying either a wild-type or a glossy fruit trait were then pooled into bulked populations and submitted to whole-genome sequencing prior to mutation frequency analysis. This revealed that the causal point mutation in the gpat6-a mutant introduces a charged amino acid adjacent to the active site of a GPAT6 enzyme. We further showed that this mutation completely abolished the GPAT activity of the recombinant protein. The gpat6-a mutant showed perturbed pollen formation but, unlike a gpat6 mutant of Arabidopsis (Arabidopsis thaliana), was not male sterile. The most striking phenotype was observed in the mutant fruit, where cuticle thickness, composition, and properties were altered. RNA sequencing analysis highlighted the main processes and pathways that were affected by the mutation at the transcriptional level, which included those associated with lipid, secondary metabolite, and cell wall biosynthesis.