Genome-wide identification, phylogeny, and expression analysis of pectin methylesterases reveal their major role in cotton fiber development

Analysis of the genetic basis of fiber-related traits and flowering time in upland cotton using machine learning

Cotton is an important crop for fiber production, but the genetic basis underlying key agronomic traits, such as fiber quality and flowering days, remains complex. While machine learning (ML) has shown great potential in uncovering the genetic architecture of complex traits in other crops, its application in cotton has been limited. Here, we applied five machine learning models-AdaBoost, Gradient Boosting Regressor, LightGBM, Random Forest, and XGBoost-to identify loci associated with fiber quality and flowering days in cotton. We compared two SNP dataset down-sampling methods for model training and found that selecting SNPs with an Fscale value greater than 0 outperformed randomly selected SNPs in terms of model accuracy. We further performed machine learning quantitative trait loci (mlQTLs) analysis for 13 traits related to fiber quality and flowering days. These mlQTLs were then compared to those identified through genome-wide association studies (GWAS), revealing that the machine learning approach not only confirmed known loci but also identified novel QTLs. Additionally, we evaluated the effect of population size on model accuracy and found that larger population sizes resulted in better predictive performance. Finally, we proposed candidate genes for the identified mlQTLs, including two argonaute 5 proteins, Gh_A09G104100 and Gh_A09G104400, for the FL3/FS2 locus, as well as GhFLA17 and Syntaxin-121 (Gh_D09G143700) for the FSD09_2/FED09_2 locus. Our findings demonstrate the efficacy of machine learning in enhancing the identification of genetic loci in cotton, providing valuable insights for improving cotton breeding strategies.

GhPME36 aggravates susceptibility to Liriomyza sativae by affecting cell wall biosynthesis in cotton leaves

BACKGROUND

Cotton is an important economic crop and a host of Liriomyza sativae. Pectin methylesterase (PME)-mediated pectin metabolism plays an indispensable role in multiple biological processes in planta. However, the pleiotropic functions of PME often lead to unpredictable effects on crop resistance to pests. Additionally, whether and how PME affects susceptibility to Liriomyza sativae remain unclear.

RESULTS

Here, we isolated GhPME36, which is located in the cell wall, from upland cotton (Gossypium hirsutum L.). Interestingly, the overexpression of GhPME36 in cotton caused severe susceptibility to Liriomyza sativae but increased leaf biomass in Arabidopsis. Cytological observations revealed that the cell wall was thinner with more demethylesterified pectins in GhPME36-OE cotton leaves than in WT leaves, whereas the soluble sugar content of GhPME36-OE cotton leaf cell walls was accordingly higher; both factors attracted Liriomyza sativae to feed on GhPME36-OE cotton leaves. Metabolomic analysis demonstrated that glucose was significantly differentially accumulated. Transcriptomic analysis further revealed DEGs enriched in glucose metabolic pathways when GhPME36 was overexpressed, suggesting that GhPME36 aggravates susceptibility to Liriomyza sativae by affecting both the structure and components of cell wall biosynthesis. Moreover, GhPME36 interacts with another pectin-modifying enzyme, GhC/VIF1, to maintain the dynamic stability of pectin methyl esterification.

CONCLUSIONS

Taken together, our results reveal the cytological and molecular mechanisms by which GhPME36 aggravates susceptibility to Liriomyza sativae. This study broadens the knowledge of PME function and provides new insights into plant resistance to pests and the safety of genetically modified plants.

Genome-wide identification of walnut (Juglans regia) PME gene family members and expression analysis during infection with Cryptosphaeria pullmanensis pathogens

Walnuts exhibit a higher resistance to diseases, though they are not completely immune. This study focuses on the Pectin methylesterase (PME) gene family to investigate whether it is involved in disease resistance in walnuts. These 21 genes are distributed across 12 chromosomes, with four pairs demonstrating homology. Variations in conserved motifs and gene structures suggest diverse functions within the gene family. Phylogenetic and collinear gene pairs of the PME family indicate that the gene family has evolved in a relatively stable way. The cis-acting elements and gene ontology enrichment of these genes, underscores their potential role in bolstering walnuts' defense mechanisms. Transcriptomic analyses were conducted under conditions of Cryptosphaeria pullmanensis infestation and verified by RT-qPCR. The results showed that certain JrPME family genes were activated in response, leading to the hypothesis that some members may confer resistance to the disease.

Nuclear factor Y-A3b binds to the SINGLE FLOWER TRUSS promoter and regulates flowering time in tomato

The control of flowering time is essential for reproductive success and has a major effect on seed and fruit yield and other important agricultural traits in crops. Nuclear factors Y (NF-Ys) are transcription factors that form heterotrimeric protein complexes to regulate gene expression required for diverse biological processes, including flowering time control in plants. However, to our knowledge, there has been no report on mutants of individual NF-YA subunits that promote early flowering phenotype in plants. In this study, we identified SlNF-YA3b, encoding a member of the NF-Y transcription factor family, as a key gene regulating flowering time in tomato. Knockout of NF-YA3b resulted in an early flowering phenotype in tomato, whereas overexpression of NF-YA3b delayed flowering in transgenic tomato plants. NF-YA3b was demonstrated to form heterotrimeric protein complexes with multiple NF-YB/NF-YC heterodimers in yeast three-hybrid assays. Biochemical evidence indicated that NF-YA3b directly binds to the CCAAT cis-elements of the SINGLE FLOWER TRUSS (SFT) promoter to suppress its gene expression. These findings uncovered a critical role of NF-YA3b in regulating flowering time in tomato and could be applied to the management of flowering time in crops.

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.

Comparative transcriptional and co-expression network analysis of two upland cotton accessions with extreme phenotypic differences reveals molecular mechanisms of fiber development

Introduction

Upland cotton (Gossypium hirsutum) is the main source of natural fiber in the global textile industry, and thus its fiber quality and yield are important parameters. In this study, comparative transcriptomics was used to analyze differentially expressed genes (DEGs) due to its ability to effectively screen candidate genes during the developmental stages of cotton fiber. However, research using this method is limited, particularly on fiber development. The aim of this study was to uncover the molecular mechanisms underlying the whole period of fiber development and the differences in transcriptional levels.

Methods

Comparative transcriptomes are used to analyze transcriptome data and to screen for differentially expressed genes. STEM and WGCNA were used to screen for key genes involved in fiber development. qRT-PCR was performed to verify gene expression of selected DEGs and hub genes.

Results

Two accessions of upland cotton with extreme phenotypic differences, namely EZ60 and ZR014121, were used to carry out RNA sequencing (RNA-seq) on fiber samples from different fiber development stages. The results identified 704, 376, 141, 269, 761, and 586 genes that were upregulated, and 1,052, 476, 355, 259, 702, and 847 genes that were downregulated at 0, 5, 10, 15, 20, and 25 days post anthesis, respectively. Similar expression patterns of DEGs were monitored using short time-series expression miner (STEM) analysis, and associated pathways of DEGs within profiles were investigated. In addition, weighted gene co-expression network analysis (WGCNA) identified five key modules in fiber development and screened 20 hub genes involved in the development of fibers.

Discussion

Through the annotation of the genes, it was found that the excessive expression of resistance-related genes in the early fiber development stages affects the fiber yield, whereas the sustained expression of cell elongation-related genes is critical for long fibers. This study provides new information that can be used to improve fibers in newly developed upland cotton genotypes.

Integrated Characterization of Cassava (Manihot esculenta) Pectin Methylesterase (MePME) Genes to Filter Candidate Gene Responses to Multiple Abiotic Stresses

Plant pectin methylesterases (PMEs) play crucial roles in regulating cell wall modification and response to various stresses. Members of the PME family have been found in several crops, but there is a lack of research into their presence in cassava (Manihot esculent), which is an important crop for world food security. In this research, 89 MePME genes were identified in cassava that were separated into two types (type-Ⅰ and type-Ⅱ) according to the existence or absence of a pro-region (PMEI domain). The MePME gene members were unevenly located on 17 chromosomes, with 19 gene pairs being identified that most likely arose via duplication events. The MePMEs could be divided into ten sub-groups in type-Ⅰ and five sub-groups in type-Ⅱ. The motif analysis revealed 11 conserved motifs in type-Ⅰ and 8 in type-Ⅱ MePMEs. The number of introns in the CDS region of type-Ⅰ MePMEs ranged between one and two, and the number of introns in type-Ⅱ MePMEs ranged between one and nine. There were 21 type-Ⅰ and 31 type-Ⅱ MePMEs that contained signal peptides. Most of the type-Ⅰ MePMEs had two conserved "RK/RLL" and one "FPSWVS" domain between the pro-region and the PME domain. Multiple stress-, hormone- and tissue-specific-related cis-acting regulatory elements were identified in the promoter regions of MePME genes. A total of five co-expressed genes (MePME1, MePME2, MePME27, MePME65 and MePME82) were filtered from different abiotic stresses via the use of UpSet Venn diagrams. The gene expression pattern analysis revealed that the expression of MePME1 was positively correlated with the degree of cassava postharvest physiological deterioration (PPD). The expression of this gene was also significantly upregulated by 7% PEG and 14 °C low-temperature stress, but slightly downregulated by ABA treatment. The tissue-specific expression analysis revealed that MePME1 and MePME65 generally displayed higher expression levels in most tissues than the other co-expressed genes. In this study, we obtain an in-depth understanding of the cassava PME gene family, suggesting that MePME1 could be a candidate gene associated with multiple abiotic tolerance.

Experimental animal models of chronic inflammation

Inflammation is a general term for a wide variety of both physiological and pathophysiological processes in the body which primarily prevents the body from diseases and helps to remove dead tissues. It has a crucial part in the body immune system. Tissue damage can recruit inflammatory cells and cytokines and induce inflammation. Inflammation can be classified as acute, sub-acute, and chronic. If it remained unresolved and lasted for prolonged periods, it would be considered as chronic inflammation (CI), which consequently exacerbates tissue damage in different organs. CI is the main pathophysiological cause of many disorders such as obesity, diabetes, arthritis, myocardial infarction, and cancer. Thus, it is critical to investigate different mechanisms involved in CI to understand its processes and to find proper anti-inflammatory therapeutic approaches for it. Animal models are one of the most useful tools for study about different diseases and mechanisms in the body, and are important in pharmacological studies to find proper treatments. In this study, we discussed the various experimental animal models that have been used to recreate CI which can help us to enhance the understanding of CI mechanisms in human and contribute to the development of potent new therapies.

The cotton pectin methyl esterase gene GhPME21 functions in microspore development and fertility in Gossypium hirsutum L

Pectin widely exists in higher plants' cell walls and intercellular space of higher plants and plays an indispensable role in plant growth and development. We identified 55 differentially expressed genes related to pectin degradation by transcriptomic analysis in the male sterile mutant, ms1. A gene encoding pectin methylesterase (GhPME21) was found to be predominantly expressed in the developing stamens of cotton but was significantly down-regulated in ms1 stamens. The tapetal layer of GhPME21 interfered lines (GhPME21i) was significantly thickened compared to that of WT at the early stage; anther compartment morphology of GhPME21i lines was abnormal, and the microspore wall was broken at the middle stage; Alexander staining showed that the pollen grains of GhPME21i lines differed greatly in volume at the late stage. The mature pollen surfaces of GhPME21i lines were deposited with discontinuous and broken sheets and prickles viewed under SEM. Fewer pollen tubes were observed to germinate in vitro in GhPME21i lines, while tiny of those in vivo were found to elongate to the ovary. The seeds harvested from GhPME21i lines as pollination donors were dry and hollow. The changes of phenotypes in GhPME21i lines at various stages illustrated that the GhPME21 gene played a vital role in the development of cotton stamens and controlled plant fertility by affecting stamen development, pollen germination, and pollen tube elongation. The findings of this study laid the groundwork for further research into the molecular mechanisms of PMEs involved in microspore formation and the creation of cotton male sterility materials.

Transcriptomic Analysis on the Peel of UV-B-Exposed Peach Fruit Reveals an Upregulation of Phenolic- and UVR8-Related Pathways

UV-B treatment deeply influences plant physiology and biochemistry, especially by activating the expression of responsive genes involved in UV-B acclimation through a UV-B-specific perception mechanism. Although the UV-B-related molecular responses have been widely studied in Arabidopsis, relatively few research reports deepen the knowledge on the influence of post-harvest UV-B treatment on fruit. In this work, a transcriptomic approach is adopted to investigate the transcriptional modifications occurring in the peel of UV-B-treated peach (Prunus persica L., cv Fairtime) fruit after harvest. Our analysis reveals a higher gene regulation after 1 h from the irradiation (88% of the differentially expressed genes-DEGs), compared to 3 h recovery. The overexpression of genes encoding phenylalanine ammonia-lyase (PAL), chalcone syntase (CHS), chalcone isomerase (CHI), and flavonol synthase (FLS) revealed a strong activation of the phenylpropanoid pathway, resulting in the later increase in the concentration of specific flavonoid classes, e.g., anthocyanins, flavones, dihydroflavonols, and flavanones, 36 h after the treatment. Upregulation of UVR8-related genes (HY5, COP1, and RUP) suggests that UV-B-triggered activation of the UVR8 pathway occurs also in post-harvest peach fruit. In addition, a regulation of genes involved in the cell-wall dismantling process (PME) is observed. In conclusion, post-harvest UV-B exposure deeply affects the transcriptome of the peach peel, promoting the activation of genes implicated in the biosynthesis of phenolics, likely via UVR8. Thus, our results might pave the way to a possible use of post-harvest UV-B treatments to enhance the content of health-promoting compounds in peach fruits and extending the knowledge of the UVR8 gene network.

Genome-wide identification and expression analysis of the ftsH protein family and its response to abiotic stress in Nicotiana tabacum L

BACKGROUND

The filamentous temperature-sensitive H protease (ftsH) gene family plays an important role in plant growth and development. FtsH proteins belong to the AAA protease family. Studies have shown that it is a key gene for plant chloroplast development and photosynthesis regulation. In addition, the ftsH gene is also involved in plant response to stress. At present, the research and analysis of the ftsH gene family are conducted in microorganisms such as Escherichia coli and Oenococcus and various plants such as Arabidopsis, pear, rice, and corn. However, analysis reports on ftsH genes from tobacco (Nicotiana tabacum L.), an important model plant, are still lacking. Since ftsH genes regulate plant growth and development, it has become necessary to systematically study this gene in an economically important plant like tobacco.

RESULTS

This is the first study to analyze the ftsH gene from Nicotiana tabacum L. K326 (NtftsH). We identified 20 ftsH genes from the whole genome sequence, renamed them according to their chromosomal locations, and divided them into eight subfamilies. These 20 NtftsH genes were unevenly distributed across the 24 chromosomes. We found four pairs of fragment duplications. We further investigated the collinearity between these genes and related genes in five other species. Quantitative real-time polymerase chain reaction (qRT-PCR) analysis identified differential expression patterns of NtftsH in different tissues and under various abiotic stress conditions.

CONCLUSIONS

This study provides a comprehensive analysis of the NtftsH gene family. The exon-intron structure and motif composition are highly similar in NtftsH genes that belong to the same evolutionary tree branch. Homology analysis and phylogenetic comparison of ftsH genes from several different plants provide valuable clues for studying the evolutionary characteristics of NtftsH genes. The NtftsH genes play important roles in plant growth and development, revealed by their expression levels in different tissues as well as under different stress conditions. Gene expression and phylogenetic analyses will provide the basis for the functional analysis of NtftsH genes. These results provide a valuable resource for a better understanding of the biological role of the ftsH genes in the tobacco plant.

Genome-wide identification of PME gene family and expression of candidate genes associated with aluminum tolerance in tea plant (Camellia sinensis)

BACKGROUND

The major aluminum (Al) detoxication mechanism of tea plant (Camellia sinensis), as an Al hyperaccumulator plant, is the fixation of almost 70% of Al in the cell walls. Pectin is the primary constituent of cell walls, a degree of methylation of pectin polysaccharides regulated by the pectin methylesterase (PME) genes can greatly affect the Al binding capacity. The knowledge on PME gene family in tea plant is still poor.

RESULTS

We identified 66 (CsPME1-CsPME66) PME genes from C. sinensis genome. We studied their protein characterization, conserved motifs, gene structure, systematic evolution and gene expression under Al treatments, to establish a basis for in-depth research on the function of PMEs in tea plant. Gene structures analysis revealed that the majority of PME genes had 2-4 exons. Phylogenetic results pointed out that the PME genes from the same species displayed comparatively high sequence consistency and genetic similarity. Selective pressure investigation suggested that the Ka/Ks value for homologous genes of PME family was less than one. The expression of CsPMEs under three Al concentration treatments was tissue specific, eight PME genes in leaves and 15 in roots displayed a trend similar to of the Al contents and PME activities under Al concentration treatments, indicating that the degree of pectin de-esterification regulated by PME was crucial for Al tolerance of tea plant.

CONCLUSIONS

Sixty-six CsPME genes were identified for the first time in tea plant. The genome-wide identification, classification, evolutionary and transcription analyses of the PME gene family provided a new direction for further research on the function of PME gene in Al tolerance of tea plant.

Revealing Genetic Differences in Fiber Elongation between the Offspring of Sea Island Cotton and Upland Cotton Backcross Populations Based on Transcriptome and Weighted Gene Coexpression Networks

Fiber length is an important indicator of cotton fiber quality, and the time and rate of cotton fiber cell elongation are key factors in determining the fiber length of mature cotton. To gain insight into the differences in fiber elongation mechanisms in the offspring of backcross populations of Sea Island cotton Xinhai 16 and land cotton Line 9, we selected two groups with significant differences in fiber length (long-fiber group L and short-fiber group S) at different fiber development stages 0, 5, 10 and 15 days post-anthesis (DPA) for transcriptome comparison. A total of 171.74 Gb of clean data was obtained by RNA-seq, and eight genes were randomly selected for qPCR validation. Data analysis identified 6055 differentially expressed genes (DEGs) between two groups of fibers, L and S, in four developmental periods, and gene ontology (GO) term analysis revealed that these DEGs were associated mainly with microtubule driving, reactive oxygen species, plant cell wall biosynthesis, and glycosyl compound hydrolase activity. Kyoto encyclopedia of genes and genomes (KEGG) pathway analysis indicated that plant hormone signaling, mitogen-activated protein kinase (MAPK) signaling, and starch and sucrose metabolism pathways were associated with fiber elongation. Subsequently, a sustained upregulation expression pattern, profile 19, was identified and analyzed using short time-series expression miner (STEM). An analysis of the weighted gene coexpression network module uncovered 21 genes closely related to fiber development, mainly involved in functions such as cell wall relaxation, microtubule formation, and cytoskeletal structure of the cell wall. This study helps to enhance the understanding of the Sea Island–Upland backcross population and identifies key genes for cotton fiber development, and these findings will provide a basis for future research on the molecular mechanisms of fiber length formation in cotton populations.

Transcriptome Time-Course Analysis in the Whole Period of Cotton Fiber Development

Gossypium hirsutum and Gossypium barbadense are the widely cultivated tetraploid cottons around the world, which evolved great differences in the fiber yield and quality due to the independent domestication process. To reveal the genetic basis of the difference, we integrated 90 samples from ten time points during the fiber developmental period for investigating the dynamics of gene expression changes associated with fiber in G. hirsutum acc. TM-1 and G. barbadense cv. Hai7124 and acc. 3-79. Globally, 44,484 genes expressed in all three cultivars account for 61.14% of the total genes. About 61.39% (N = 3,412) of the cotton transcription factors were involved in fiber development, which consisted of 58 cotton TF families. The differential analysis of intra- and interspecies showed that 3 DPA had more expression changes. To discover the genes with temporally changed expression profiles during the whole fiber development, 1,850 genes predominantly expressed in G. hirsutum and 1,050 in G. barbadense were identified, respectively. Based on the weighted gene co-expression network and time-course analysis, several candidate genes, mainly involved in the secondary cell wall synthesis and phytohormones, were identified in this study, underlying possibly the transcriptional regulation and molecular mechanisms of the fiber quality differences between G. barbadense and G. hirsutum. The quantitative real-time PCR validation of the candidate genes was consistent with the RNA-seq data. Our study provides a strong rationale for the analysis of gene function and breeding of high-quality cotton.

Genome-wide identification of PME genes, evolution and expression analyses in soybean (Glycine max L.)

BACKGROUND

Pectin methylesterase (PME) is one of pectin-modifying enzyme that affects the pectin homeostasis in cell wall and regulates plant growth and diverse biological processes. The PME genes have been well explored and characterized in different plants. Nevertheless, systematic research on the soybean (Glycine max L.) PME genes remain lacking.

RESULTS

We identified 127 Glycine max PME genes (GmPME) from the soybean Wm82.a2.v1 genome, which unevenly distributed on 20 soybean chromosomes. Phylogenetic analysis classified the GmPME genes into four clades (Group I, Group II, Group III and Group IV). GmPME gene members in the same clades displayed similar gene structures and motif patterns. The gene family expansion analysis demonstrated that segmental duplication was the major driving force to acquire novel GmPME genes compared to the tandem duplication events. Further synteny and evolution analyses showed that the GmPME gene family experienced strong purifying selective pressures during evolution. The cis-element analyses together with the expression patterns of the GmPME genes in various tissues suggested that the GmPME genes broadly participate in distinct biological processes and regulate soybean developments. Importantly, based on the transcriptome data and quantitative RT-PCR validations, we examined the potential roles of the GmPME genes in regulating soybean flower bud development and seed germination.

CONCLUSION

In conclusion, we provided a comprehensive characterization of the PME genes in soybean, and our work laid a foundation for the functional study of GmPME genes in the future.

Co-expression network and comparative transcriptome analysis for fiber initiation and elongation reveal genetic differences in two lines from upland cotton CCRI70 RIL population

Upland cotton is the most widely planted for natural fiber around the world, and either lint percentage (LP) or fiber length (FL) is the crucial component tremendously affecting cotton yield and fiber quality, respectively. In this study, two lines MBZ70-053 and MBZ70-236 derived from G. hirsutum CCRI70 recombinant inbred line (RIL) population presenting different phenotypes in LP and FL traits were chosen to conduct RNA sequencing on ovule and fiber samples, aiming at exploring the differences of molecular and genetic mechanisms during cotton fiber initiation and elongation stages. As a result, 249/128, 369/206, 4296/1198 and 3547/2129 up-/down- regulated differentially expressed genes (DGEs) in L2 were obtained at -3, 0, 5 and 10 days post-anthesis (DPA), respectively. Seven gene expression profiles were discriminated using Short Time-series Expression Miner (STEM) analysis; seven modules and hub genes were identified using weighted gene co-expression network analysis. The DEGs were mainly enriched into energetic metabolism and accumulating as well as auxin signaling pathway in initiation and elongation stages, respectively. Meanwhile, 29 hub genes were identified as 14-3-3ω , TBL35, GhACS, PME3, GAMMA-TIP, PUM-7, etc., where the DEGs and hub genes revealed the genetic and molecular mechanisms and differences during cotton fiber development.

Genome-wide identification of the BASS gene family in four Gossypium species and functional characterization of GhBASSs against salt stress

Bile acid sodium symporter (BASS) family proteins encode a class of sodium/solute symporters. Even though the sodium transporting property of BASSs in mammals was well studied, their sodium transportability and functional roles in plant salt tolerance remained largely unknown. Here, BASS family members from 4 cotton species, as well as 30 other species were identified. Then, they were designated as members of BASS1 to BASS5 subfamilies according to their sequence similarity and phylogenetic relationships. There were 8, 11, 16 and 18 putative BASS genes in four cotton species. While whole-genome duplications (WGD) and segmental duplications rendered the expansion of the BASS gene family in cotton, BASS gene losses occurred in the tetraploid cotton during the evolution from diploids to allotetraploids. Concerning functional characterizations, the transcript profiling of GhBASSs revealed that they not only preferred tissue-specific expression but also were differently induced by various stressors and phytohormones. Gene silencing and overexpression experiments showed that GhBASS1 and GhBASS3 positively regulated, whereas GhBASS2, GhBASS4 and GhBASS5 negatively regulated plant salt tolerance. Taken together, BASS family genes have evolved before the divergence from the common ancestor of prokaryotes and eukaryotes, and GhBASSs are plastidial sodium-dependent metabolite co-transporters that can influence plant salt tolerance.

The cellulose synthase (CesA) gene family in four Gossypium species: phylogenetics, sequence variation and gene expression in relation to fiber quality in Upland cotton

Cellulose synthases (CesAs) are multi-subunit enzymes found on the plasma membrane of plant cells and play a pivotal role in cellulose production. The cotton fiber is mainly composed of cellulose, and the genetic relationships between CesA genes and cotton fiber yield and quality are not fully understood. Through a phylogenetic analysis, the CesA gene family in diploid Gossypium arboreum and Gossypium raimondii, as well as tetraploid Gossypium hirsutum ('TM-1') and Gossypium barbadense ('Hai-7124' and '3-79'), was divided into 6 groups and 15 sub-groups, with each group containing two to five homologous genes. Most CesA genes in the four species are highly collinear. Among the five cotton genomes, 440 and 1929 single nucleotide polymorphisms (SNPs) in the CesA gene family were identified in exons and introns, respectively, including 174 SNPs resulting in amino acid changes. In total, 484 homeologous SNPs between the A and D genomes were identified in diploids, while 142 SNPs were detected between the two tetraploids, with 32 and 82 SNPs existing within G. hirsutum and G. barbadense, respectively. Additionally, 74 quantitative trait loci near 18 GhCesA genes were associated with fiber quality. One to four GhCesA genes were differentially expressed (DE) in ovules at 0 and 3 days post anthesis (DPA) between two backcross inbred lines having different fiber lengths, but no DE genes were identified between these lines in developing fibers at 10 DPA. Twenty-seven SNPs in above DE CesA genes were detected among seven cotton lines, including one SNP in Ghi_A08G03061 that was detected in four G. hirsutum genotypes. This study provides the first comprehensive characterization of the cotton CesA gene family, which may play important roles in determining cotton fiber quality.

Silencing of GhKEA4 and GhKEA12 Revealed Their Potential Functions Under Salt and Potassium Stresses in Upland Cotton

The K⁺ efflux antiporter (KEA) mediates intracellular K⁺ and H⁺ homeostasis to improve salt tolerance in plants. However, the knowledge of KEA gene family in cotton is largely absent. In the present study, 8, 8, 15, and 16 putative KEA genes were identified in Gossypium arboreum, G. raimondii, G. hirsutum, and G. barbadense, respectively. These KEA genes were classified into three subfamilies, and members from the same subfamilies showed similar motif compositions and gene structure characteristics. Some hormone response elements and stress response elements were identified in the upstream 2000 bp sequence of GhKEAs. Transcriptome data showed that most of the GhKEAs were highly expressed in roots and stems. The quantificational real-time polymerase chain reaction (qRT-PCR) results showed that most of the GhKEAs responded to low potassium, salt and drought stresses. Virus-induced gene silencing (VIGS) experiments demonstrated that under salt stress, after silencing genes GhKEA4 and GhKEA12, the chlorophyll content, proline content, soluble sugar content, peroxidase (POD) activity and catalase (CAT) activity were significantly decreased, and the Na⁺/K⁺ ratio was extremely significantly increased in leaves, leading to greater salt sensitivity. Under high potassium stress, cotton plants silenced for the GhKEA4 could still maintain a more stable Na⁺ and K⁺ balance, and the activity of transporting potassium ions from roots into leaves was reduced silenced for GhKEA12. Under low potassium stress, silencing the GhKEA4 increased the activity of transporting potassium ions to shoots, and silencing the GhKEA12 increased the ability of absorbing potassium ions, but accumulated more Na⁺ in leaves. These results provided a basis for further studies on the biological roles of KEA genes in cotton development and adaptation to stress conditions.

Genome-wide characterization of the cellulose synthase gene superfamily in Pyrus bretschneideri and reveal its potential role in stone cell formation

Members of the cellulose synthase (CesA) and cellulose synthase-like (Csl) families from the cellulose synthase gene superfamily participate in cellulose and hemicellulose synthesis in the plasma membrane. The members of this superfamily are vital for cell wall construction during plant growth and development. However, little is known about their function in pear fruit, a model for Rosaceae species and for fleshy fruit development. In our research, a total of 36 CesA/Csl family members were identified from the pear and were grouped into six subfamilies (CesA, CslB, CslC, CslD, CslE, and CslG) according to phylogenetic relationships. We performed a protein sequence physicochemical analysis, phylogenetic tree construction, a gene structure, a conserved domain, and chromosomal localization analysis. The results indicated that most of the CesA/Csl genes from pear are closely related to genes in Arabidopsis, but these families have unique characteristics in terms of their gene structure, chromosomal localization, phylogeny, and deduced protein sequences, suggesting that they have evolved through different processes. Tissue expression analysis results showed that most of the CesA/Csl genes were constitutively expressed at different levels in different organs. Furthermore, the expression levels of four genes (Pbr032894.2, Pbr016107.1, Pbr00518.1, and Pbr034218.1) tended to first increase and then decrease during fruit development, implying that these four genes may be involved in the development of stone cells of pear fruit. Our results may help elucidate the evolutionary history and functional differences of the CesA/Csl genes in pear and lay a foundation for further investigation of the CesA/Csl genes in pear and other Rosaceae species.

Characterization of the pectin methyl-esterase gene family and its function in controlling pollen tube growth in pear (Pyrus bretschneideri)

Pectin methyl-esterases (PMEs) play crucial roles in plant growth. In this study, we identified 81 PbrPMEs in pear. Whole-genome duplication and purifying selection drove the evolution of PbrPME gene family. The expression of 47 PbrPMEs was detected in pear pollen tube, which were assigned to 13 clusters by an expression tendency analysis. One of the 13 clusters presented opposite expression trends towards the changes of methyl-esterified pectins at the apical cell wall. PbrPMEs were localized in the cytoplasm and plasma membrane. Repression of PbrPME11, PbrPME44, and PbrPME59 resulted in the inhibition of pear pollen tube growth and abnormal deposition of methyl-esterified pectins at pollen tube tip. Pharmacological analysis confirmed that reduced PbrPME activities repressed the pollen tube growth. Overall, we have explored the evolutionary characteristics of PbrPME gene family and found the key PbrPME genes that control the growth of pollen tube, which deepened the understanding of pear fertility regulation.

Genome wide identification and functional characterization of strawberry pectin methylesterases related to fruit softening

BACKGROUND

Pectin methylesterase (PME) is a hydrolytic enzyme that catalyzes the demethylesterification of homogalacturonans and controls pectin reconstruction, being essential in regulation of cell wall modification. During fruit ripening stage, PME-mediated cell wall remodeling is an important process to determine fruit firmness and softening. Strawberry fruit is a soft fruit with a short postharvest life, due to a rapid loss of firm texture. Hence, preharvest improvement of strawberry fruit rigidity is a prerequisite for extension of fruit refreshing time. Although PME has been well characterized in model plants, knowledge regarding the functionality and evolutionary property of PME gene family in strawberry remain limited.

RESULTS

A total of 54 PME genes (FvPMEs) were identified in woodland strawberry (Fragaria vesca 'Hawaii 4'). Phylogeny and gene structure analysis divided these FvPME genes into four groups (Group 1-4). Duplicate events analysis suggested that tandem and dispersed duplications effectively contributed to the expansion of the PME family in strawberry. Through transcriptome analysis, we identified FvPME38 and FvPME39 as the most abundant-expressed PMEs at fruit ripening stages, and they were positively regulated by abscisic acid. Genetic manipulation of FvPME38 and FvPME39 by overexpression and RNAi-silencing significantly influences the fruit firmness, pectin content and cell wall structure, indicating a requirement of PME for strawberry fruit softening.

CONCLUSION

Our study globally analyzed strawberry pectin methylesterases by the approaches of phylogenetics, evolutionary prediction and genetic analysis. We verified the essential role of FvPME38 and FvPME39 in regulation of strawberry fruit softening process, which provided a guide for improving strawberry fruit firmness by modifying PME level.

Genome-wide quantitative trait loci reveal the genetic basis of cotton fibre quality and yield-related traits in a Gossypium hirsutum recombinant inbred line population

Cotton is widely cultivated globally because it provides natural fibre for the textile industry and human use. To identify quantitative trait loci (QTLs)/genes associated with fibre quality and yield, a recombinant inbred line (RIL) population was developed in upland cotton. A consensus map covering the whole genome was constructed with three types of markers (8295 markers, 5197.17 centimorgans (cM)). Six fibre yield and quality traits were evaluated in 17 environments, and 983 QTLs were identified, 198 of which were stable and mainly distributed on chromosomes 4, 6, 7, 13, 21 and 25. Thirty-seven QTL clusters were identified, in which 92.8% of paired traits with significant medium or high positive correlations had the same QTL additive effect directions, and all of the paired traits with significant medium or high negative correlations had opposite additive effect directions. In total, 1297 genes were discovered in the QTL clusters, 414 of which were expressed in two RNA-Seq data sets. Many genes were discovered, 23 of which were promising candidates. Six important QTL clusters that included both fibre quality and yield traits were identified with opposite additive effect directions, and those on chromosome 13 (qClu-chr13-2) could increase fibre quality but reduce yield; this result was validated in a natural population using three markers. These data could provide information about the genetic basis of cotton fibre quality and yield and help cotton breeders to improve fibre quality and yield simultaneously.

Genome-wide identification, phylogeny and expression analysis of the PME and PMEI gene families in maize

Pectins, the major components of cell walls in plants, are synthesized and secreted to cell walls as highly methyl-esterified polymers and then demethyl-esterified by pectin methylesterases (PMEs). The PMEs are spatially regulated by pectin methylesterase inhibitors (PMEIs). In this study, 43 and 49 putative PME and PMEI genes were identified in maize, respectively. Gene structure and motif analysis revealed that members in the same paralogous pairs or in the same subgroup generally had common motif compositions and gene structure patterns, which indicates functional similarity between the closely related ZmPME/PMEI genes. Gene ontology annotation analysis showed that most of the ZmPME/PMEI genes are involved in cell wall modification and pectin catabolic process with molecular functions of pectinesterase or pectinesterase inhibitor activities. There are 35 ZmPME/PMEI genes expressed higher in anthers than in other tissues from the NimbleGen maize microarray data, and the semiq-RT-PCR assay revealed most of these ZmPME/PMEIs specially expressed in anthers and pollens, indicating they possibly had role in anther and pollen development. In addition, these ZmPME/PMEI genes were highly expressed in the fertile anthers, while lowly or no expressed in sterile anthers. This further indicated these genes might be involved in the development of anther and pollen.

Long non-coding RNAs and their potential functions in Ligon-lintless-1 mutant cotton during fiber development

BACKGROUND

Long non-coding RNAs (LncRNAs) are part of genes, which are not translated into proteins and play a vital role in plant growth and development. Nevertheless, the presence of LncRNAs and how they functions in Ligon-lintless-1 mutant during the early cessation of cotton fiber development are still not well understood. In order to investigate the function of LncRNAs in cotton fiber development, it is necessary and important to identify LncRNAs and their potential roles in fiber cell development.

RESULTS

In this work, we identified 18,333 LncRNAs, with the proportion of long intergenic noncoding RNAs (LincRNAs) (91.5%) and anti-sense LncRNAs (8.5%), all transcribed from Ligon-lintless-1 (Li1) and wild-type (WT). Expression differences were detected between Ligon-lintless-1 and wild-type at 0 and 8 DPA (day post anthesis). Pathway analysis and Gene Ontology based on differentially expressed LncRNAs on target genes, indicated fatty acid biosynthesis and fatty acid elongation being integral to lack of fiber in mutant cotton. The result of RNA-seq and RT-qPCR clearly singles out two potential LncRNAs, LNC_001237 and LNC_017085, to be highly down-regulated in the mutant cotton. The two LncRNAs were found to be destabilized or repressed by ghr-miR2950. Both RNA-seq analysis and RT-qPCR results in Ligon-lintless-1 mutant and wild-type may provide strong evidence of LNC_001237, LNC_017085 and ghr-miR2950 being integral molecular elements participating in various pathways of cotton fiber development.

CONCLUSION

The results of this study provide fundamental evidence for the better understanding of LncRNAs regulatory role in the molecular pathways governing cotton fiber development. Further research on designing and transforming LncRNAs will help not only in the understanding of their functions but will also in the improvement of fiber quality.

Genome-Wide Identification and Expression Profiling of the Polygalacturonase (PG) and Pectin Methylesterase (PME) Genes in Grapevine (Vitis vinifera L.)

In pectin regulation, polygalacturonases (PGs) and pectin methylesterases (PMEs) are critical components in the transformation, disassembly network, and remodeling of plant primary cell walls. In the current study, we identified 36 PG and 47 PME genes using the available genomic resources of grapevine. Herein, we provide a comprehensive overview of PGs and PMEs, including phylogenetic and collinearity relationships, motif and gene structure compositions, gene duplications, principal component analysis, and expression profiling during developmental stages. Phylogenetic analysis of PGs and PMEs revealed similar domain composition patterns with Arabidopsis. The collinearity analysis showed high conservation and gene duplications with purifying selection. The type of duplications also varied in terms of gene numbers in PGs (10 dispersed, 1 proximal, 12 tandem, and 13 segmental, respectively) and PMEs (23 dispersed, 1 proximal, 16 tandem, and 7 segmental, respectively). The tissue-specific response of PG and PME genes based on the reported transcriptomic data exhibited diverged expression patterns in various organs during different developmental stages. Among PGs, VvPG8, VvPG10, VvPG13, VvPG17, VvPG18, VvPG19, VvPG20, VvPG22, and VvPG23 showed tissue- or organ-specific expression in majority of the tissues during development. Similarly, in PMEs, VvPME3, VvPME4, VvPME5, VvPME6, VvPME19, VvPME21, VvPME23, VvPME29, VvPME31, and VvPME32 suggested high tissue-specific response. The gene ontology (GO), Kyoto Encyclopedia of Genes and Genomics (KEGG) enrichment, and cis-elements prediction analysis also suggested the putative functions of PGs and PMEs in plant development, such as pectin and carbohydrate metabolism, and stress activities. Moreover, qRT-PCR validation of 32 PG and PME genes revealed their role in various organs of grapevines (i.e., root, stem, tendril, inflorescence, flesh, skins, and leaves). Therefore, these findings will lead to novel insights and encourage cutting-edge research on functional characterization of PGs and PMEs in fruit crop species.

A high density SLAF-SNP genetic map and QTL detection for fibre quality traits in Gossypium hirsutum

Background

Upland Cotton (Gossypium hirsutum) is a very important cash crop known for its high quality natural fiber. Recent advances in sequencing technologies provide powerful tools with which to explore the cotton genome for single nucleotide polymorphism marker identification and high density genetic map construction toward more reliable quantitative trait locus mapping.

Results

In the present study, a RIL population was developed by crossing a Chinese high fiber quality cultivar (Yumian 1) and an American high fiber quality line (CA3084), with distinct genetic backgrounds. Specific locus amplified fragment sequencing (SLAF-seq) technology was used to discover SNPs, and a genetic map containing 6254 SNPs was constructed, covering 3141.72 cM with an average distance of 0.5 cM between markers. A total of 95 QTL were detected for fiber quality traits in three environments, explaining 5.5-24.6% of the phenotypic variance. Fifty-five QTL found in multiple environments were considered stable QTL. Nine of the stable QTL were found in all three environments. We identified 14 QTL clusters on 13 chromosomes, each containing one or more stable QTL.

Conclusion

A high-density genetic map of Gossypium hirsutum developed by using specific locus amplified fragment sequencing technology provides detailed mapping of fiber quality QTL, and identification of ‘stable QTL’ found in multiple environments. A marker-rich genetic map provides a foundation for fine mapping, candidate gene identification and marker-assisted selection of favorable alleles at stable QTL in breeding programs.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5294-5) contains supplementary material, which is available to authorized users.

Pectin Methylesterases: Cell Wall Remodeling Proteins Are Required for Plant Response to Heat Stress

Heat stress (HS) is expected to be of increasing worldwide concern in the near future, especially with regard to crop yield and quality as a consequence of rising or varying temperatures as a result of global climate change. HS response (HSR) is a highly conserved mechanism among different organisms but shows remarkable complexity and unique features in plants. The transcriptional regulation of HSR is controlled by HS transcription factors (HSFs) which allow the activation of HS-responsive genes, among which HS proteins (HSPs) are best characterized. Cell wall remodeling constitutes an important component of plant responses to HS to maintain overall function and growth; however, little is known about the connection between cell wall remodeling and HSR. Pectin controls cell wall porosity and has been shown to exhibit structural variation during plant growth and in response to HS. Pectin methylesterases (PMEs) are present in multigene families and encode isoforms with different action patterns by removal of methyl esters to influencing the properties of cell wall. We aimed to elucidate how plant cell walls respond to certain environmental cues through cell wall-modifying proteins in connection with modifications in cell wall machinery. An overview of recent findings shed light on PMEs contribute to a change in cell-wall composition/structure. The fine-scale modulation of apoplastic calcium ions (Ca2+) content could be mediated by PMEs in response to abiotic stress for both the assembly and disassembly of the pectic network. In particular, this modulation is prevalent in guard cell walls for regulating cell wall plasticity as well as stromal aperture size, which comprise critical determinants of plant adaptation to HS. These insights provide a foundation for further research to reveal details of the cell wall machinery and stress-responsive factors to provide targets and strategies to facilitate plant adaptation.

Methanol induces cytosolic calcium variations, membrane depolarization and ethylene production in arabidopsis and tobacco

Background and Aims

Methanol is a volatile organic compound released from plants through the action of pectin methylesterases (PMEs), which demethylesterify cell wall pectins. Plant PMEs play a role in developmental processes but also in responses to herbivory and infection by fungal or bacterial pathogens. However, molecular mechanisms that explain how methanol could affect plant defences remain poorly understood.

Methods

Using cultured cells and seedlings from Arabidopsis thaliana and tobacco BY2 expressing the apoaequorin gene, allowing quantification of cytosolic Ca2+, a reactive oxygen species (ROS) probe (CLA, Cypridina luciferin analogue) and electrophysiological techniques, we followed early plant cell responses to exogenously supplied methanol applied as a liquid or as volatile.

Key Results

Methanol induces cytosolic Ca2+ variations that involve Ca2+ influx through the plasma membrane and Ca2+ release from internal stores. Our data further suggest that these Ca2+ variations could interact with different ROS and support a signalling pathway leading to well known plant responses to pathogens such as plasma membrane depolarization through anion channel regulation and ethylene synthesis.

Conclusions

Methanol is not only a by-product of PME activities, and our data suggest that [Ca2+]cyt variations could participate in signalling processes induced by methanol upstream of plant defence responses.

Genome-Wide Identification and Characterization of Pectin Methylesterase Inhibitor Genes in Brassica oleracea

The activities of pectin methylesterases (PMEs) are regulated by pectin methylesterase inhibitors (PMEIs), which consequently control the pectin methylesterification status. However, the role of PMEI genes in Brassica oleracea, an economically important vegetable crop, is poorly understood. In this study, 95 B. oleracea PMEI (BoPMEI) genes were identified. A total of 77 syntenic ortholog pairs and 10 tandemly duplicated clusters were detected, suggesting that the expansion of BoPMEI genes was mainly attributed to whole-genome triplication (WGT) and tandem duplication (TD). During diploidization after WGT, BoPMEI genes were preferentially retained in accordance with the gene balance hypothesis. Most homologous gene pairs experienced purifying selection with ω (Ka/Ks) ratios lower than 1 in evolution. Five stamen-specific BoPMEI genes were identified by expression pattern analysis. By combining the analyses of expression and evolution, we speculated that nonfunctionalization, subfunctionalization, neofunctionalization, and functional conservation can occur in the long evolutionary process. This work provides insights into the characterization of PMEI genes in B. oleracea and contributes to the further functional studies of BoPMEI genes.

Genome-wide identification and expression analyses of the pectate lyase (PEL) gene family in cotton (Gossypium hirsutum L.)

Background

Pectin is a major component and structural polysaccharide of the primary cell walls and middle lamella of higher plants. Pectate lyase (PEL, EC 4.2.2.2), a cell wall modification enzyme, degrades de-esterified pectin for cell wall loosening, remodeling and rearrangement. Nevertheless, there have been few studies on PEL genes and no comprehensive analysis of the PEL gene family in cotton.

Results

We identified 53, 42 and 83 putative PEL genes in Gossypium raimondii (D5), Gossypium arboreum (A2), and Gossypium hirsutum (AD1), respectively. These PEL genes were classified into five subfamilies (I-V). Members from the same subfamilies showed relatively conserved gene structures, motifs and protein domains. An analysis of gene chromosomal locations and gene duplication revealed that segmental duplication likely contributed to the expansion of the GhPELs. The 2000 bp upstream sequences of all the GhPELs contained auxin response elements. A transcriptomic data analysis showed that 62 GhPELs were expressed in various tissues. Notably, most (29/32) GhPELs of subfamily IV were preferentially expressed in the stamen, and five GhPELs of subfamily V were prominently expressed at the fiber elongation stage. In addition, qRT-PCR analysis revealed the expression characteristics of 24 GhPELs in four pollen developmental stages and significantly different expression of some GhPELs between long- and short-fiber cultivars. Moreover, some members were responsive to IAA treatment. The results indicate that GhPELs play significant and functionally diverse roles in the development of different tissues.

Conclusions

In this study, we comprehensively analyzed PELs in G. hirsutum, providing a foundation to better understand the functions of GhPELs in different tissues and pathways, especially in pollen, fiber and the auxin signaling pathway.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5047-5) contains supplementary material, which is available to authorized users.

QTL delineation for five fiber quality traits based on an intra-specific Gossypium hirsutum L. recombinant inbred line population

Gossypium hirsutum L. is the most important fiber crop worldwide and contributes to more than 95% of global cotton production. Marker-assisted selection (MAS) is an effective approach for improving fiber quality, and quantitative trait loci (QTL) mapping of fiber quality traits is important for cotton breeding. In this study, a permanent intra-specific recombinant inbred line (RIL) population containing 137 families was used for fiber quality testing. Based on a previously reported high-density genetic map with an average marker distance of 0.63 cM, 186 additive QTLs were obtained for five fiber quality traits over five consecutive years, including 39 for fiber length (FL), 36 for fiber strength (FS), 50 for fiber uniformity (FU), 33 for micronaire (MC) and 28 for fiber elongation (FE). Three stable QTLs, qMC-A4-1, qMC-D2-3 and qFS-D9-1, were detected in four datasets, and another eight stable QTLs, qMC-A4-2, qMC-D11-2, qFU-A9-1, qFU-A10-4, qFS-D11-1, qFL-D9-2, qFL-D11-1 and qFE-A3-2, were detected in three datasets. The annotated genes in these 11 stable QTLs were collected, and these genes included many transcription factors with functions during fiber development. 33 QTL coincidence regions were found, and these involved nearly half of the total QTLs. Four chromosome regions containing at least 6 QTLs were promising for fine mapping. In addition, 41 pairs of epistatic QTLs (e-QTLs) were screened, including 6 for FL, 30 for FS, 2 for FU and 3 for MC. The identification of stable QTLs adds valuable information for further QTL fine mapping and gene positional cloning for fiber quality genetic detection and provides useful markers for further molecular breeding in enhancing fiber quality.

Genome-wide identification and analysis of the evolution and expression patterns of the cellulose synthase gene superfamily in Gossypium species

The cellulose synthase gene superfamily, which includes the cellulose synthase (Ces) and cellulose synthase-like (Csl) families, is involved in the synthesis of cellulose and hemicellulose. This superfamily is critical for cotton fiber development in Gossypium species. Applying a series of bioinformatic methods, we identified 228 Ces/Csl genes from four Gossypium species (G. hirsutum, G. barbadense, G. arboreum, and G. raimondii). These genes were then grouped into 11 subfamilies based on phylogenetic relationships. A subsequent analysis of gene evolution revealed sites in CSLG and CSLJ genes that were under long-term positive selection pressure, with a posterior probability >0.95. Moreover, the d_N:d_S value for the CSLJ clade was 1.305, suggesting this subfamily was under positive selection pressure. Our data indicated that the d_N:d_S value ranged from 0.0084 to 0.9693 among the homologous Ces/Csl genes, implying they were under purifying selection pressure. Our transcriptome and qRT-PCR analyses revealed that CesA genes were more highly expressed in tetraploids than in diploids. However, the Csl expression levels exhibited the opposite trend. Furthermore, changes to promoter sequences may have influenced the expression of homologous Ces/Csl genes. Our findings may provide novel insights into the evolutionary relationships and expression patterns of the Ces/Csl genes in Gossypium species.

Genome-wide identification of CBL family and expression analysis of CBLs in response to potassium deficiency in cotton

Calcineurin B-like (CBL) proteins, as calcium sensors, play pivotal roles in plant responses to diverse abiotic stresses and in growth and development through interaction with CBL-interacting protein kinases (CIPKs). However, knowledge about functions and evolution of CBLs in Gossypium plants is scarce. Here, we conducted a genome-wide survey and identified 13, 13 and 22 CBL genes in the progenitor diploid Gossypium arboreum and Gossypium raimondii, and the cultivated allotetraploid Gossypium hirsutum, respectively. Analysis of physical properties, chromosomal locations, conserved domains and phylogeny indicated rather conserved nature of CBLs among the three Gossypium species. Moreover, these CBLs have closer genetic evolutionary relationship with the CBLs from cocoa than with those from other plants. Most CBL genes underwent evolution under purifying selection in the three Gossypium plants. Additionally, nearly all G. hirsutum CBL (GhCBL) genes were expressed in the root, stem, leaf, flower and fiber. Many GhCBLs were preferentially expressed in the flower while several GhCBLs were mainly expressed in roots. Expression patterns of GhCBL genes in response to potassium deficiency were also studied. The expression of most GhCBLs were moderately induced in roots after treatments with low-potassium stress. Yeast two-hybrid experiments indicated that GhCBL1-2, GhCBL1-3, GhCBL4-4, GhCBL8, GhCBL9 and GhCBL10-3 interacted with GhCIPK23, respectively. Our results provided a comprehensive view of the CBLs and valuable information for researchers to further investigate the roles and functional mechanisms of the CBLs in Gossypium.

The Intergenic Interplay between Aldose 1-Epimerase-Like Protein and Pectin Methylesterase in Abiotic and Biotic Stress Control

The mechanical damage that often precedes the penetration of a leaf by a pathogen promotes the activation of pectin methylesterase (PME); the activation of PME leads to the emission of methanol, resulting in a "priming" effect on intact leaves, which is accompanied by an increased sensitivity to Tobacco mosaic virus (TMV) and resistance to bacteria. In this study, we revealed that mRNA levels of the methanol-inducible gene encoding Nicotiana benthamiana aldose 1-epimerase-like protein (NbAELP) in the leaves of intact plants are very low compared with roots. However, stress and pathogen attack increased the accumulation of the NbAELP mRNA in the leaves. Using transiently transformed plants, we obtained data to support the mechanism underlying AELP/PME-related negative feedback The insertion of the NbAELP promoter sequence (proNbAELP) into the N. benthamiana genome resulted in the co-suppression of the natural NbAELP gene expression, accompanied by a reduction in the NbAELP mRNA content and increased PME synthesis. Knockdown of NbAELP resulted in high activity of PME in the cell wall and a decrease in the leaf glucose level, creating unfavorable conditions for Agrobacterium tumefaciens reproduction in injected leaves. Our results showed that NbAELP is capable of binding the TMV movement protein (MP_TMV) in vitro and is likely to affect the cellular nucleocytoplasmic transport, which may explain the sensitivity of NbAELP knockdown plants to TMV. Although NbAELP was primarily detected in the cell wall, the influence of this protein on cellular PME mRNA levels might be associated with reduced transcriptional activity of the PME gene in the nucleus. To confirm this hypothesis, we isolated the N. tabacum PME gene promoter (proNtPME) and showed the inhibition of proNtPME-directed GFP and GUS expression in leaves when co-agroinjected with the NbAELP-encoding plasmid. We hypothesized that plant wounding and/or pathogen attack lead to PME activation and increased methanol emission, followed by increased NbAELP expression, which results in reversion of PME mRNA level and methanol emission to levels found in the intact plant.