Metabolic GWAS-based dissection of genetic basis underlying nutrient quality variation and domestication of cassava storage root.

BACKGROUND: Metabolites play critical roles in regulating nutritional qualities of plants, thereby influencing their consumption and human health. However, the genetic basis underlying the metabolite-based nutrient quality and domestication of root and tuber crops remain largely unknown. RESULTS: We report a comprehensive study combining metabolic and phenotypic genome-wide association studies to dissect the genetic basis of metabolites in the storage root (SR) of cassava. We quantify 2,980 metabolic features in 299 cultivated cassava accessions. We detect 18,218 significant marker-metabolite associations via metabolic genome-wide association mapping and identify 12 candidate genes responsible for the levels of metabolites that are of potential nutritional importance. Me3GT, MeMYB4, and UGT85K4/UGT85K5, which are involved in flavone, anthocyanin, and cyanogenic glucoside metabolism, respectively, are functionally validated through in vitro enzyme assays and in vivo gene silencing analyses. We identify a cluster of cyanogenic glucoside biosynthesis genes, among which CYP79D1, CYP71E7b, and UGT85K5 are highly co-expressed and their allelic combination contributes to low linamarin content. We find MeMYB4 is responsible for variations in cyanidin 3-O-glucoside and delphinidin 3-O-rutinoside contents, thus controlling SR endothelium color. We find human selection affects quercetin 3-O-glucoside content and SR weight per plant. The candidate gene MeFLS1 is subject to selection during cassava domestication, leading to decreased quercetin 3-O-glucoside content and thus increased SR weight per plant. CONCLUSIONS: These findings reveal the genetic basis of cassava SR metabolome variation, establish a linkage between metabolites and agronomic traits, and offer useful resources for genetically improving the nutrition of cassava and other root crops.

Genome-wide identification and characterization of 14-3-3 gene family related to negative regulation of starch accumulation in storage root of Manihot esculenta.

The 14-3-3 protein family is a highly conservative member of the acid protein family and plays an important role in regulating a series of important biological activities and various signal transduction pathways. The role of 14-3-3 proteins in regulating starch accumulation still remains largely unknown. To investigate the properties of 14-3-3 proteins, the structures and functions involved in starch accumulation in storage roots were analyzed, and consequently, 16 Me14-3-3 genes were identified. Phylogenetic analysis revealed that Me14-3-3 family proteins are split into two groups (ε and non-ε). All Me14-3-3 proteins contain nine antiparallel α-helices. Me14-3-3s-GFP fusion protein was targeted exclusively to the nuclei and cytoplasm. In the early stage of starch accumulation in the storage root, Me14-3-3 genes were highly expressed in high-starch cultivars, while in the late stage of starch accumulation, Me14-3-3 genes were highly expressed in low-starch cultivars. Me14-3-3 I, II, V, and XVI had relatively high expression levels in the storage roots. The transgenic evidence from Me14-3-3II overexpression in Arabidopsis thaliana and the virus-induced gene silencing (VIGS) in cassava leaves and storage roots suggest that Me14-3-3II is involved in the negative regulation of starch accumulation. This study provides a new insight to understand the molecular mechanisms of starch accumulation linked with Me14-3-3 genes during cassava storage root development.

Cassava Foliage Effects on Antioxidant Capacity, Growth, Immunity, and Ruminal Microbial Metabolism in Hainan Black Goats.

Cassava (Manihot esculenta Crantz) foliage is a byproduct of cassava production characterized by high biomass and nutrient content. In this study, we investigated the effects of cassava foliage on antioxidant capacity, growth performance, and immunity status in goats, as well as rumen fermentation and microbial metabolism. Twenty-five Hainan black goats were randomly divided into five groups (n = 5 per group) and accepted five treatments: 0% (T1), 25% (T2), 50% (T3), 75% (T4), and 100% (T5) of the cassava foliage silage replaced king grass, respectively. The feeding experiment lasted for 70 d (including 10 d adaptation period and 60 d treatment period). Feeding a diet containing 50% cassava foliage resulted in beneficial effects for goat growth and health, as reflected by the higher average daily feed intake (ADFI), average daily gain (ADG) and better feed conversion rate (FCR), as well as by the reduced serum levels of alanine aminotransferase (ALT), aspartate aminotransferase (AST), creatinine (CRE), and triglycerides (TG). Meanwhile, cassava foliage improved antioxidant activity by increasing the level of glutathion peroxidase (GSH-Px), superoxide dismutase (SOD), and total antioxidant capacity (T-AOC) and lowering malondialdehyde (MDA). Moreover, feeding cassava foliage was also beneficial to immunity status by enhancing complement 3 (C3), complement 4 (C4), immunoglobulin A (IgA), immunoglobulin G (IgG), and immunoglobulin M (IgM). Furthermore, the addition of dietary cassava foliage also altered rumen fermentation, rumen bacterial community composition, and metabolism. The abundance of Butyrivibrio_2 and Prevotella_1 was elevated, as were the concentrations of beneficial metabolites such as butyric acid; there was a concomitant decline in metabolites that hindered nutrient metabolism and harmed host health. In summary, goats fed a diet containing 50% cassava foliage silage demonstrated a greater abundance of Butyrivibrio_2, which enhanced the production of butyric acid; these changes led to greater antioxidant capacity, growth performance, and immunity in the goats.

Spraying chitosan on cassava roots reduces postharvest deterioration by promoting wound healing and inducing disease resistance.

Postharvest damage makes cassava roots vulnerable to pathogen infections and decay, which significantly hinders the development of the cassava industry. The objective of this study was to assess the antibacterial properties of chitosan in vitro, as well as its effect on wound healing and resistance in cassava roots. The findings demonstrated that the bacteriostatic effect of chitosan became increasingly prominent as the concentration of chitosan enhanced. Chitosan at a concentration of 0.5 mg/mL was revealed to significantly inhibit the germination of P. palmivora spores and damage to their structure. Moreover, chitosan activated the transcription of crucial genes and enzyme activities associated with the phenylpropane metabolism pathway in cassava roots, thus promoting rapid lignin accumulation and resulting in the early formation of a fracture layer. Chitosan was also found to enhance cassava root resistance by promoting the expression of pathogenesis-related proteins and the accumulation of flavonoids and total phenols. After 48 h of inoculation, cassava roots treated with chitosan exhibited a 51.4 % and 53.4 % decrease in lesion area for SC9 and SC6 varieties, respectively. The findings of this study offer a novel approach for managing postharvest deterioration of cassava roots.

Integrative analysis of metabolome and transcriptome reveals the mechanism of color formation in cassava (Manihot esculenta Crantz) leaves.

Cassava (Manihot esculenta Crantz) leaves are often used as vegetables in Africa. Anthocyanins possess antioxidant, anti-inflammatory, anti-cancer, and other biological activities. They are poor in green leaves but rich in the purple leaves of cassava. The mechanism of anthocyanin's accumulation in cassava is poorly understood. In this study, two cassava varieties, SC9 with green leaves and Ziyehuangxin with purple leaves (PL), were selected to perform an integrative analysis using metabolomics and transcriptomics. The metabolomic analysis indicated that the most significantly differential metabolites (SDMs) belong to anthocyanins and are highly accumulated in PL. The transcriptomic analysis revealed that differentially expressed genes (DEGs) are enriched in secondary metabolites biosynthesis. The analysis of the combination of metabolomics and transcriptomics showed that metabolite changes are associated with the gene expressions in the anthocyanin biosynthesis pathway. In addition, some transcription factors (TFs) may be involved in anthocyanin biosynthesis. To further investigate the correlation between anthocyanin accumulation and color formation in cassava leaves, the virus-induced gene silencing (VIGS) system was used. VIGS-MeANR silenced plant showed the altered phenotypes of cassava leaves, partially from green to purple color, resulting in a significant increase of the total anthocyanin content and reduction in the expression of MeANR. These results provide a theoretical basis for breeding cassava varieties with anthocyanin-rich leaves.

Magnesium chelatase subunit D is not only required for chlorophyll biosynthesis and photosynthesis, but also affecting starch accumulation in Manihot esculenta Crantz.

BACKGROUND: Magnesium chelatase plays an important role in photosynthesis, but only a few subunits have been functionally characterized in cassava. RESULTS: Herein, MeChlD was successfully cloned and characterized. MeChlD encodes a magnesium chelatase subunit D, which has ATPase and vWA conservative domains. MeChlD was highly expressed in the leaves. Subcellular localization suggested that MeChlD:GFP was a chloroplast-localized protein. Furthermore, the yeast two-hybrid system and BiFC analysis indicated that MeChlD interacts with MeChlM and MePrxQ, respectively. VIGS-induce silencing of MeChlD resulted in significantly decreased chlorophyll content and reduction the expression of photosynthesis-related nuclear genes. Furthermore, the storage root numbers, fresh weight and the total starch content in cassava storage roots of VIGS-MeChlD plants was significantly reduced. CONCLUSION: Taken together, MeChlD located at the chloroplast is not only required for chlorophyll biosynthesis and photosynthesis, but also affecting the starch accumulation in cassava. This study expands our understanding of the biological functions of ChlD proteins.

A major QTL identification and candidate gene analysis of watermelon fruit cracking using QTL-seq and RNA-seq.

Fruit cracking decreases the total production and the commercial value of watermelon. The molecular mechanisms of fruit cracking are unknown. In this study, 164 recombinant inbred lines (RILs) of watermelon, derived from the crossing of the WQ1 (cracking-sensitive) and WQ2 (cracking-tolerant) lines, were sequenced using specific length amplified fragment sequencing (SLAF-seq). A high-density genetic linkage map was constructed with 3,335 markers spanning 1,322.74 cM, at an average 0.40 cM across whole-genome flanking markers. The cracking tolerance capacity (CTC), depth of fruit cracking (DFC), rind thickness (RT), and rind hardness (RH) were measured for quantitative trait locus (QTL) analysis. Of the four traits analyzed, one major QTL with high phenotypic variation (41.04%-61.37%) was detected at 76.613-76.919 cM on chromosome 2, which contained 104 annotated genes. Differential gene expression analysis with RNA sequencing (RNA-seq) data between the two parents identified 4,508 differentially expressed genes (DEGs). Comparison of the genes between the QTL region and the DEGs obtained eight coexisting genes. Quantitative real-time PCR (qRT-PCR) analysis revealed that these genes were significant differentially expressed between the two parents. These results provide new insights into the identification of QTLs or genes and marker-assisted breeding in watermelon.

Integrated Metabolomic and Transcriptomic Analyses Reveals Sugar Transport and Starch Accumulation in Two Specific Germplasms of Manihot esculenta Crantz.

As a starchy and edible tropical plant, cassava (Manihot esculenta Crantz) has been widely used as an industrial raw material and a dietary source. However, the metabolomic and genetic differences in specific germplasms of cassava storage root were unclear. In this study, two specific germplasms, M. esculenta Crantz cv. sugar cassava GPMS0991L and M. esculenta Crantz cv. pink cassava BRA117315, were used as research materials. Results showed that sugar cassava GPMS0991L was rich in glucose and fructose, whereas pink cassava BRA117315 was rich in starch and sucrose. Metabolomic and transcriptomic analysis indicated that sucrose and starch metabolism had significantly changing metabolites enrichment and the highest degree of differential expression genes, respectively. Sugar transport in storage roots may contribute to the activities of sugar, which will eventually be exported to transporters (SWEETs), such as (MeSWEET1a, MeSWEET2b, MeSWEET4, MeSWEET5, MeSWEET10b, and MeSWEET17c), which transport hexose to plant cells. The expression level of genes involved in starch biosynthesis and metabolism were altered, which may result in starch accumulation. These results provide a theoretical basis for sugar transport and starch accumulation and may be useful in improving the quality of tuberous crops and increasing yield.

Genome-wide analysis of the metallothionein gene family in cassava reveals its role in response to physiological stress through the regulation of reactive oxygen species.

BACKGROUND: Cassava (Manihot esculenta Crantz) is widely planted in tropical and several subtropical regions in which drought, high temperatures, and other abiotic stresses occur. Metallothionein (MT) is a group of conjugated proteins with small molecular weight and rich in cysteine. These proteins play a substantial role in response to physiological stress through the regulation of reactive oxygen species (ROS). However, the biological functions of MT genes in cassava are unknown. RESULTS: A total of 10 MeMT genes were identified in the cassava genome. The MeMTs were divided into 3 groups (Types 2-4) based on the contents and distribution of Cys residues. The MeMTs exhibited tissue-specific expression and located on 7 chromosomes. The MeMT promoters contain some hormones regulatory and stresses responsiveness elements. MeMTs were upregulated under hydrogen peroxide (H2O2) treatment and in respond to post-harvest physiological deterioration (PPD). The results were consistent with defense-responsive cis-acting elements in the MeMT promoters. Further, four of MeMTs were selected and silenced by using the virus-induced gene silencing (VIGS) method to evaluate their functional characterization. The results of gene-silenced cassava suggest that MeMTs are involved in oxidative stress resistance, as ROS scavengers. CONCLUSION: We identified the 10 MeMT genes, and explore their evolutionary relationship, conserved motif, and tissue-specific expression. The expression profiles of MeMTs under three kinds of abiotic stresses (wounding, low-temperature, and H2O2) and during PPD were analyzed. The tissue-specific expression and the response to abiotic stresses revealed the role of MT in plant growth and development. Furthermore, silenced expression of MeMTs in cassava leaves decreased its tolerance to ROS, consistent with its predicted role as ROS scavengers. In summary, our results suggest an important role of MeMTs in response to physiological stress as well as species adaptation via the regulation of ROS homeostasis.

Overexpression of leucoanthocyanidin reductase or anthocyanidin reductase elevates tannins content and confers cassava resistance to two-spotted spider mite.

The two-spotted spider mite (TSSM) is a destructive cassava pest. Intensive demonstration of resistance mechanism greatly facilitates the creation of TSSM-resistant cassava germplasm. Gene to metabolite network plays a crucial role in modulating plant resistance, but little is known about the genes and related metabolites which are responsible for cassava resistance to TSSM. Here, a highly resistant (HR) and a highly susceptible (HS) cassava cultivar were used, integrative and comparative transcriptomic and metabolomic analyses between these two cultivars after TSSM infestation revealed that several genes and metabolites were closely related and significantly different in abundance. In particular, the expression of leucoanthocyanidin reductase (LAR) and anthocyanidin reductase (ANR) genes showed a high positive correlation with most of the metabolites in the tannin biosynthesis pathway. Furthermore, transgenic cassava lines overexpressing either of the genes elevated tannin concentrations and conferred cassava resistance to TSSM. Additionally, different forms of tannins possessed distinct bioactivity on TSSM, of which total condensed tannins (LC50 = 375.68 mg/l) showed maximum lethal effects followed by procyanidin B1 (LC50 = 3537.10 mg/l). This study accurately targets LAR, ANR and specific tannin compounds as critical genes and metabolites in shaping cassava resistance to TSSM, which could be considered as biomarkers for evaluation and creation of pest-resistant cassava germplasm.

Systematic Analysis of bHLH Transcription Factors in Cassava Uncovers Their Roles in Postharvest Physiological Deterioration and Cyanogenic Glycosides Biosynthesis.

The basic helix-loop-helix (bHLH) proteins are a large superfamily of transcription factors, and play a central role in a wide range of metabolic, physiological, and developmental processes in higher organisms. However, systematic investigation of bHLH gene family in cassava (Manihot esculenta Crantz) has not been reported. In the present study, we performed a genome-wide survey and identified 148 MebHLHs genes were unevenly harbored in 18 chromosomes. Through phylogenetic analyses along with Arabidopsis counterparts, these MebHLHs genes were divided into 19 groups, and each gene contains a similar structure and conserved motifs. Moreover, many cis-acting regulatory elements related to various defense and stress responses showed in MebHLH genes. Interestingly, transcriptome data analyses unveiled 117 MebHLH genes during postharvest physiological deterioration (PPD) process of cassava tuberous roots, while 65 MebHLH genes showed significantly change. Meanwhile, the relative quantitative analysis of 15 MebHLH genes demonstrated that they were sensitive to PPD, suggesting they may involve in PPD process regulation. Cyanogenic glucosides (CGs) biosynthesis during PPD process was increased, silencing of MebHLH72 and MebHLH114 showed that linamarin content was significantly decreased in the leaves. To summarize, the genome-wide identification and expression profiling of MebHLH candidates pave a new avenue for uderstanding their function in PPD and CGs biosynthesis, which will accelerate the improvement of PPD tolerance and decrease CGs content in cassava tuberous roots.

The transcription factor MebHLH18 in cassava functions in decreasing low temperature-induced leaf abscission to promote low-temperature tolerance.

The reactive oxygen species (ROS) signal regulates stress-induced leaf abscission in cassava. The relationship between the function of the cassava transcription factor bHLH gene and low temperature-induced leaf abscission is still unclear. Here, we report that MebHLH18, a transcription factor, involved in regulating low temperature-induced leaf abscission in cassava. The expression of the MebHLH18 gene was significantly related to low temperature-induced leaf abscission and POD level. Under low temperatures, the levels of ROS scavengers in different cassava genotypes were significantly different in the low temperature-induced leaf abscission process. Cassava gene transformation showed that MebHLH18 overexpression significantly decreased the low temperature-induced leaf abscission rate. Simultaneously, interference expression increased the rate of leaf abscission under the same conditions. ROS analysis showed a connection between the decrease in the low temperature-induced leaf abscission rate caused by MebHLH18 expression and the increase in antioxidant activity. A Genome-wide association studies analysis showed a relationship between the natural variation of the promoter region of MebHLH18 and low temperature-induced leaf abscission. Furthermore, studies showed that the change in MebHLH18 expression was caused by a single nucleotide polymorphism variation in the promoter region upstream of the gene. The high expression of MebHLH18 led to a significant increase in POD activity. The increased POD activity decreased the accumulation of ROS at low temperatures and the rate of leaf abscission. It indicates that the natural variation in the promoter region of MebHLH18 increases antioxidant levels under low temperatures and slows down low temperature-induced leaf abscission.

Different Fertilizers Applied Alter Fungal Community Structure in Rhizospheric Soil of Cassava (Manihot esculenta Crantz) and Increase Crop Yield.

Soil microbes play an important role in the ecosystem and have a relationship with plant growth, development, and production. There are only a few reports on the effects of planting patterns of cassava on the microbial community structure in the rhizospheric soil. Here, we investigated the effects of different fertilization on the microbial community structure in the cassava rhizospheric soil. SC205 cultivar was used in this study as the experimental material. Compound fertilizer (CF) and reduced fertilizer (RF) were applied to the soil prior to planting. Soil samples were collected before harvest, and fungi were analyzed using IonS5TMXL sequencing platform. Results showed that CF and RF treatments significantly increased cassava yield. Amplicon sequencing result indicated that the fungi richness in rhizospheric soil of cassava was increased after CF was applied, and the diversity was decreased. However, the fungal diversity and richness were decreased in rhizospheric soil after RF was applied. The most dominant fungal phylum was Ascomycota, which increased after fertilization. In addition, the abundance of beneficial fungi such as Chaetomium increased after fertilization, while that of pathogenic fungi such as Fusarium solani was decreased. The composition of the fungal community in rhizospheric soil with CF and RF applied was similar, but the richness and diversity of fungi were different. Canonical correspondence analysis (CCA) indicates there was a positive correlation between soil nutrition and fungal community structure. Overall, our results indicate that fertilization alters the fungal community structure of cassava rhizospheric soil, such that the abundance of potentially beneficial fungi increased, while that of potentially pathogenic fungi decreased, thereby significantly promoting plant growth and yield of cassava. Thus, during actual production, attention should be paid to maintain the stability of cassava rhizospheric soil micro-ecology.

Resequencing of 388 cassava accessions identifies valuable loci and selection for variation in heterozygosity.

BACKGROUND: Heterozygous genomes are widespread in outcrossing and clonally propagated crops. However, the variation in heterozygosity underlying key agronomic traits and crop domestication remains largely unknown. Cassava is a staple crop in Africa and other tropical regions and has a highly heterozygous genome. RESULTS: We describe a genomic variation map from 388 resequenced genomes of cassava cultivars and wild accessions. We identify 52 loci for 23 agronomic traits through a genome-wide association study. Eighteen allelic variations in heterozygosity for nine candidate genes are significantly associated with seven key agronomic traits. We detect 81 selective sweeps with decreasing heterozygosity and nucleotide diversity, harboring 548 genes, which are enriched in multiple biological processes including growth, development, hormone metabolisms and responses, and immune-related processes. Artificial selection for decreased heterozygosity has contributed to the domestication of the large starchy storage root of cassava. Selection for homozygous GG allele in MeTIR1 during domestication contributes to increased starch content. Selection of homozygous AA allele in MeAHL17 is associated with increased storage root weight and cassava bacterial blight (CBB) susceptibility. We have verified the positive roles of MeTIR1 in increasing starch content and MeAHL17 in resistance to CBB by transient overexpression and silencing analysis. The allelic combinations in MeTIR1 and MeAHL17 may result in high starch content and resistance to CBB. CONCLUSIONS: This study provides insights into allelic variation in heterozygosity associated with key agronomic traits and cassava domestication. It also offers valuable resources for the improvement of cassava and other highly heterozygous crops.

Cassava breeding and agronomy in Asia: 50 years of history and future directions.

In Asia, cassava (Manihot esculenta) is cultivated by more than 8 million farmers, driving the rural economy of many countries. The International Center for Tropical Agriculture (CIAT), in partnership with national agricultural research institutes (NARIs), instigated breeding and agronomic research in Asia, 1983. The breeding program has successfully released high-yielding cultivars resulting in an average yield increase from 13.0 t ha-1 in 1996 to 21.3 t ha-1 in 2016, with significant economic benefits. Following the success in increasing yields, cassava breeding has turned its focus to higher-value traits, such as waxy cassava, to reach new market niches. More recently, building resistance to invasive pests and diseases has become a top priority due to the emergent threat of cassava mosaic disease (CMD). The agronomic research involves driving profitability with advanced technologies focusing on better agronomic management practices thereby maintaining sustainable production systems. Remote sensing technologies are being tested for trait discovery and large-scale field evaluation of cassava. In summary, cassava breeding in Asia is driven by a combination of food and market demand with technological innovations to increase the productivity. Further, exploration in the potential of data-driven agriculture is needed to empower researchers and producers for sustainable advancement.

Endophytic bacterial and fungal microbiota in different cultivars of cassava (Manihot esculenta Crantz).

Endophytes colonize tissues of healthy host plants and play a crucial role in plant growth and development. However, little attention has been paid to the endophytes of tuber crops such as cassava, which is used as a staple food by approximately 800 million people worldwide. This study aimed to elucidate the diversity and composition of endophytic bacterial and fungal communities in different cassava cultivars using high-throughput sequencing. Although no significant differences in richness or diversity were observed among the different cassava cultivars, the community compositions were diverse. Two cultivars (SC124 and SC205) tolerant to root rot exhibited similar community compositions, while two other cultivars (SC10 and SC5), which are moderately and highly susceptible to root rot, respectively, harboured similar community compositions. Proteobacteria, Firmicutes, and Ascomycota dominated the endophyte assemblages, with Weissella, Serratia, Lasiodiplodia, Fusarium, and Diaporthe being the predominant genera. The differentially abundant taxonomic clades between the tolerant and susceptible cultivars were mainly rare taxa, such as Lachnoclostridium_5, Rhizobium, Lampropedia, and Stenotrophomonas. These seemed to be key genera that affected the susceptibility of cassava to root rot. Moreover, the comparison of KEGG functional profiles revealed that 'Environmental adaptation' category was significantly enriched in the tolerant cultivars, while 'Infectious diseases: Parasitic' category was significantly enriched in the susceptible cultivars. The present findings open opportunities for further studies on the roles of endophytes in the susceptibility of plants to diseases.

Relevance of Class I α-Mannosidases to Cassava Postharvest Physiological Deterioration.

Class I α-mannosidases (MNSs) play important roles in protein N-glycosylation. However, no data are currently available about MNSs in cassava (Manihot esculenta), of which the functions are therefore not known, particularly in relevance to postharvest physiological deterioration (PPD). A total of seven genes were identified from the cassava genome in the present study. Two (MeMNS2 and MeMNS6) of the seven genes may be pseudogenes, as indicated by sequence alignment and exon-intron organizations. Five MNSs could be classified into three subfamilies. Tissue-specific expression analysis revealed that MNS genes have distinct expression patterns in different tissues between sugar cassava and cultivated cassava varieties, indicating their functional diversity. A PPD response and defense model was proposed based on the transcription data of MNSs and genes involved in reactive oxygen species, signal transduction, and cell wall remodeling. The findings help in the understanding of PPD responses in cassava.

Protein Cross-Interactions for Efficient Photosynthesis in the Cassava Cultivar SC205 Relative to Its Wild Species.

The underlying mechanisms of the higher photosynthetic efficiency of cultivated cassava relative to its wild species are poorly understood. In the present study, proteins in leaves and chloroplasts were analyzed to compare the differences among the cultivar SC205, its wild ancestor W14, and the related species Glaziovii. The functions of differential proteins are associated with 10 ontology groups including photosynthesis, carbohydrate and energy metabolism, as well as potential signal pathway. The protein-protein networks among 41 differential proteins showed that PGK1 is a hub protein and protein cross-interactions affected the differentiation of photosynthetic rate. Anatomy patterns and PEPC detection suggested that SC205 has more C4 photosynthesis characteristics than Glaziovii and W14. Finally, a mechanism model of the efficient photosynthesis was proposed based on the remarkable variations in photosynthetic parameters and protein functions in the domestic cultivars.

Comparative physiology and transcriptome analysis allows for identification of lncRNAs imparting tolerance to drought stress in autotetraploid cassava.

BACKGROUND: Polyploidization, pervasive among higher plant species, enhances adaptation to water deficit, but the physiological and molecular advantages need to be investigated widely. Long non-coding RNAs (lncRNAs) are involved in drought tolerance in various crops. RESULTS: Herein, we demonstrate that tetraploidy potentiates tolerance to drought stress in cassava (Manihot esculenta Crantz). Autotetraploidy reduces transpiration by lesser extent increasing of stomatal density, smaller stomatal aperture size, or greater stomatal closure, and reducing accumulation of H2O2 under drought stress. Transcriptome analysis of autotetraploid samples revealed down-regulation of genes involved in photosynthesis under drought stress, and less down-regulation of subtilisin-like proteases involved in increasing stomatal density. UDP-glucosyltransferases were increased more or reduced less in dehydrated leaves of autotetraploids compared with controls. Strand-specific RNA-seq data (validated by quantitative real time PCR) identified 2372 lncRNAs, and 86 autotetraploid-specific lncRNAs were differentially expressed in stressed leaves. The co-expressed network analysis indicated that LNC_001148 and LNC_000160 in autotetraploid dehydrated leaves regulated six genes encoding subtilisin-like protease above mentioned, thereby result in increasing the stomatal density to a lesser extent in autotetraploid cassava. Trans-regulatory network analysis suggested that autotetraploid-specific differentially expressed lncRNAs were associated with galactose metabolism, pentose phosphate pathway and brassinosteroid biosynthesis, etc. CONCLUSION: Tetraploidy potentiates tolerance to drought stress in cassava, and LNC_001148 and LNC_000160 mediate drought tolerance by regulating stomatal density in autotetraploid cassava.