An optimized isolation protocol yields high-quality RNA from cassava tissues (Manihot esculenta Crantz)

Comparative transcriptomics analysis reveals defense mechanisms of Manihot esculenta Crantz against Sri Lanka Cassava MosaicVirus

BACKGROUND

Cassava mosaic disease (CMD), caused by Sri Lankan cassava mosaic virus (SLCMV) infection, has been identified as a major pernicious disease in Manihot esculenta Crantz (cassava) plantations. It is widespread in Southeast Asia, especially in Thailand, which is one of the main cassava supplier countries. With the aim of restricting the spread of SLCMV, we explored the gene expression of a tolerant cassava cultivar vs. a susceptible cassava cultivar from the perspective of transcriptional regulation and the mechanisms underlying plant immunity and adaptation.

RESULTS

Transcriptomic analysis of SLCMV-infected tolerant (Kasetsart 50 [KU 50]) and susceptible (Rayong 11 [R 11]) cultivars at three infection stages-that is, at 21 days post-inoculation (dpi) (early/asymptomatic), 32 dpi (middle/recovery), and 67 dpi (late infection/late recovery)-identified 55,699 expressed genes. Differentially expressed genes (DEGs) between SLCMV-infected KU 50 and R 11 cultivars at (i) 21 dpi to 32 dpi (the early to middle stage), and (ii) 32 dpi to 67 dpi (the middle stage to late stage) were then identified and validated by real-time quantitative PCR (RT-qPCR). DEGs among different infection stages represent genes that respond to and regulate the viral infection during specific stages. The transcriptomic comparison between the tolerant and susceptible cultivars highlighted the role of gene expression regulation in tolerant and susceptible phenotypes.

CONCLUSIONS

This study identified genes involved in epigenetic modification, transcription and transcription factor activities, plant defense and oxidative stress response, gene expression, hormone- and metabolite-related pathways, and translation and translational initiation activities, particularly in KU 50 which represented the tolerant cultivar in this study.

Integrative transcriptomics reveals association of abscisic acid and lignin pathways with cassava whitefly resistance

BACKGROUND

Whiteflies are a global threat to crop yields, including the African subsistence crop cassava (Manihot esculenta). Outbreaks of superabundant whitefly populations throughout Eastern and Central Africa in recent years have dramatically increased the pressures of whitefly feeding and virus transmission on cassava. Whitefly-transmitted viral diseases threaten the food security of hundreds of millions of African farmers, highlighting the need for developing and deploying whitefly-resistant cassava. However, plant resistance to whiteflies remains largely poorly characterized at the genetic and molecular levels. Knowledge of cassava-defense programs also remains incomplete, limiting characterization of whitefly-resistance mechanisms. To better understand the genetic basis of whitefly resistance in cassava, we define the defense hormone- and Aleurotrachelus socialis (whitefly)-responsive transcriptome of whitefly-susceptible (COL2246) and whitefly-resistant (ECU72) cassava using RNA-seq. For broader comparison, hormone-responsive transcriptomes of Arabidopsis thaliana were also generated.

RESULTS

Whitefly infestation, salicylic acid (SA), jasmonic acid (JA), ethylene (ET), and abscisic acid (ABA) transcriptome responses of ECU72 and COL2246 were defined and analyzed. Strikingly, SA responses were largely reciprocal between the two cassava genotypes and we suggest candidate regulators. While susceptibility was associated with SA in COL2246, resistance to whitefly in ECU72 was associated with ABA, with SA-ABA antagonism observed. This was evidenced by expression of genes within the SA and ABA pathways and hormone levels during A. socialis infestation. Gene-enrichment analyses of whitefly- and hormone-responsive genes suggest the importance of fast-acting cell wall defenses (e.g., elicitor recognition, lignin biosynthesis) during early infestation stages in whitefly-resistant ECU72. A surge of ineffective immune and SA responses characterized the whitefly-susceptible COL2246's response to late-stage nymphs. Lastly, in comparison with the model plant Arabidopsis, cassava's hormone-responsive genes showed striking divergence in expression.

CONCLUSIONS

This study provides the first characterization of cassava's global transcriptome responses to whitefly infestation and defense hormone treatment. Our analyses of ECU72 and COL2246 uncovered possible whitefly resistance/susceptibility mechanisms in cassava. Comparative analysis of cassava and Arabidopsis demonstrated that defense programs in Arabidopsis may not always mirror those in crop species. More broadly, our hormone-responsive transcriptomes will also provide a baseline for the cassava community to better understand global responses to other yield-limiting pests/pathogens.

An optimized method for high-quality RNA extraction from distinctive intrinsic laryngeal muscles in the rat model

Challenges related to high-quality RNA extraction from post-mortem tissue have limited RNA-sequencing (RNA-seq) application in certain skeletal muscle groups, including the intrinsic laryngeal muscles (ILMs). The present study identified critical factors contributing to substandard RNA extraction from the ILMs and established a suitable method that permitted high-throughput analysis. Here, standard techniques for tissue processing were adapted, and an effective means to control confounding effects during specimen preparation was determined. The experimental procedure consistently provided sufficient intact total RNA (N = 68) and RIN ranging between 7.0 and 8.6, which was unprecedented using standard RNA purification protocols. This study confirmed the reproducibility of the workflow through repeated trials at different postnatal time points and across the distinctive ILMs. High-throughput diagnostics from 90 RNA samples indicated no sequencing alignment scores below 70%, validating the extraction strategy. Significant differences between the standard and experimental conditions suggest circumvented challenges and broad applicability to other skeletal muscles. This investigation remains ongoing given the prospect of therapeutic insights to voice, swallowing, and airway disorders. The present methodology supports pioneering global transcriptome investigations in the larynx previously unfounded in literature.

Characterization of cassava ORANGE proteins and their capability to increase provitamin A carotenoids accumulation

Cassava (Manihot esculenta Crantz) biofortification with provitamin A carotenoids is an ongoing process that aims to alleviate vitamin A deficiency. The moderate content of provitamin A carotenoids achieved so far limits the contribution to providing adequate dietary vitamin A levels. Strategies to increase carotenoid content focused on genes from the carotenoids biosynthesis pathway. In recent years, special emphasis was given to ORANGE protein (OR), which promotes the accumulation of carotenoids and their stability in several plants. The aim of this work was to identify, characterize and investigate the role of OR in the biosynthesis and stabilization of carotenoids in cassava and its relationship with phytoene synthase (PSY), the rate-limiting enzyme of the carotenoids biosynthesis pathway. Gene and protein characterization of OR, expression levels, protein amounts and carotenoids levels were evaluated in roots of one white (60444) and two yellow cassava cultivars (GM5309-57 and GM3736-37). Four OR variants were found in yellow cassava roots. Although comparable expression was found for three variants, significantly higher OR protein amounts were observed in the yellow varieties. In contrast, cassava PSY1 expression was significantly higher in the yellow cultivars, but PSY protein amount did not vary. Furthermore, we evaluated whether expression of one of the variants, MeOR_X1, affected carotenoid accumulation in cassava Friable Embryogenic Callus (FEC). Overexpression of maize PSY1 alone resulted in carotenoids accumulation and induced crystal formation. Co-expression with MeOR_X1 led to greatly increase of carotenoids although PSY1 expression was high in the co-expressed FEC. Our data suggest that posttranslational mechanisms controlling OR and PSY protein stability contribute to higher carotenoid levels in yellow cassava. Moreover, we showed that cassava FEC can be used to study the efficiency of single and combinatorial gene expression in increasing the carotenoid content prior to its application for the generation of biofortified cassava with enhanced carotenoids levels.

An optimized nucleic acid isolation protocol for virus diagnostics in cassava (Manihot esculenta Crantz.)

Our group works on the detection and characterization of cassava viruses, supporting projects that involve large scale pathogen surveillance activities and resistance screening assays in multiple and remote locations. In order to comply with these applications, nucleic acid isolation protocols need to be cost effective, adjusted for samples that will stand long distance transport and harsh storage conditions, while maximizing the yield and quality of the nucleic acid extracts obtained. The method we describe here has been widely used and validated using different downstream tests (including, but not limited to, Rolling Circle Amplification and Illumina and Nanopore sequencing), but is currently unpublished. The protocol begins with milligram amounts of dry leaf samples stored in silica gel, does not require liquid Nitrogen nor phenol extraction and produces an average of 2.11 µg of nucleic acids per mg of dry tissue.•DNA purity estimations reveal OD260/280 ratios above 2.0 and OD260/230 ratios above 1.7, even for samples stored in silica gel for several months.•The high quality of the extracts is suitable for detection of DNA and RNA viruses, with high efficiency.•We suggest this method could be used as part of a gold standard kit for virus detection in cassava.

Efficient method for isolation of high-quality RNA from Psidium guajava L. tissues

Acquiring high-quality RNA in sufficient amounts is crucial in plant molecular biology and genetic studies. Several methods for RNA extraction from plants are available in the literature, mainly due to the great biochemical diversity present in each species and tissue, which can complicate or prevent the extraction. Psidium guajava (Myrtaceae family) is a perennial fruit tree of medicinal and economic value; nevertheless, only a few molecular studies are available for the species. One reason is the difficulty in obtaining RNA due to the content of the samples, which are rich in polyphenols, polysaccharides, and secondary metabolites. Furthermore, there are few studies available for the isolation of RNA from guava or Psidium samples, which hampers advances in the study of the genus. Here, quality and yields of RNA isolates were compared using six extraction protocols: two protocols based on the application of cetyltrimethylammonium bromide (CTAB) lysis buffer, one protocol which uses the TRIzol reagent, one which applies guanidine thiocyanate lysis buffer followed by organic phase extraction, and two commercial kits (PureLink RNA Mini Kit and RNeasy Plant Mini Kit). The CTAB-based method provided the highest RNA yields and quality for five different tissues (flower bud, immature leaf, young leaf, mature leaf, and root), genotypes, and stress conditions. For the most efficient protocol, the average yield of RNA from guava leaves was 203.06 μg/g of tissue, and the A₂₆₀/A₂₈₀ and A₂₆₀/A₂₃₀ ratios were 2.1 and 2.2, respectively. RT-qPCR analysis demonstrated that the purity of the samples was sufficient for molecular biology experiments. CTAB-based methods for RNA isolation were found to be the most efficient, providing the highest RNA yields and quality for tissues from P. guajava. Additionally, they were compatible for downstream RNA-based applications, besides being simple and cost-effective.

Field-transcriptome analyses reveal developmental transitions during flowering in cassava (Manihot esculenta Crantz)

We characterized genes that function in the photoperiodic flowering pathway in cassava. Transcriptome analysis of field-grown plants revealed characteristic expression patterns of these genes, demonstrating that field-grown cassava experiences two distinct developmental transitions. Cassava is an important crop for both edible and industrial purposes. Cassava develops storage roots that accumulate starch, providing an important source of staple food in tropical regions. To facilitate cassava breeding, it is important to elucidate how flowering is controlled. Several important genes that control flowering time have been identified in model plants; however, comprehensive characterization of these genes in cassava is still lacking. In this study, we identified genes encoding central flowering time regulators and examined these sequences for the presence or absence of conserved motifs. We found that cassava shares conserved genes for the photoperiodic flowering pathway, including florigen, anti-florigen and its associated transcription factor (GIGANTEA, CONSTANS, FLOWERING LOCUS T, CENTRORADIALIS/TERMINAL FLOWER1 and FD) and florigen downstream genes (SUPRESSOR OF OVEREXPRESSION OF CONSTANS1 and APETALA1/FRUITFUL). We conducted RNA-seq analysis of field-grown cassava plants and characterized the expression of flowering control genes. Finally, from the transcriptome analysis we identified two distinct developmental transitions that occur in field-grown cassava.

Comprehensive genomic resources related to domestication and crop improvement traits in Lima bean

Lima bean (Phaseolus lunatus L.), one of the five domesticated Phaseolus bean crops, shows a wide range of ecological adaptations along its distribution range from Mexico to Argentina. These adaptations make it a promising crop for improving food security under predicted scenarios of climate change in Latin America and elsewhere. In this work, we combine long and short read sequencing technologies with a dense genetic map from a biparental population to obtain the chromosome-level genome assembly for Lima bean. Annotation of 28,326 gene models show high diversity among 1917 genes with conserved domains related to disease resistance. Structural comparison across 22,180 orthologs with common bean reveals high genome synteny and five large intrachromosomal rearrangements. Population genomic analyses show that wild Lima bean is organized into six clusters with mostly non-overlapping distributions and that Mesomerican landraces can be further subdivided into three subclusters. RNA-seq data reveal 4275 differentially expressed genes, which can be related to pod dehiscence and seed development. We expect the resources presented here to serve as a solid basis to achieve a comprehensive view of the degree of convergent evolution of Phaseolus species under domestication and provide tools and information for breeding for climate change resiliency.

A universal method for high-quality RNA extraction from plant tissues rich in starch, proteins and fiber

Using existing protocols, RNA extracted from seeds rich in starch often results in poor quality RNA, making it inappropriate for downstream applications. Though some methods are proposed for extracting RNA from plant tissue rich in starch and other polysaccharides, they invariably yield less and poor quality RNA. In order to obtain high yield and quality RNA from seeds and other plant tissues including roots a modified SDS-LiCl method was compared with existing methods, including TRIZOL kit (Invitrogen), Plant RNeasy mini kit (Qiagen), Furtado (2014) method, and CTAB-LiCl method. Modifications in the extraction buffer and solutions used for RNA precipitation resulted in a robust method for extracting RNA in seeds and roots, where extracting quality RNA is challenging. The modified SDS-LiCl method revealed intense RNA bands through gel electrophoresis and a nanodrop spectrophotometer detected ratios of ≥ 2 and 1.8 for A260/A230 and A260/A280, respectively. The absence of starch co-precipitation during RNA extraction resulted in enhanced yield and quality of RNA with RIN values of 7-9, quantified using a bioanalyzer. The high-quality RNA obtained was demonstrated to be suitable for downstream applications, such as cDNA synthesis, gene amplification, and RT-qPCR. The method was also effective in extracting RNA from seeds of other cereals including field-grown sorghum and corn. The modified SDS-LiCl method is a robust and highly reproducible RNA extraction method for plant tissues rich in starch and other secondary metabolites. The modified SDS-LiCl method successfully extracted high yield and quality RNA from mature, developing, and germinated seeds, leaves, and roots exposed to different abiotic stresses.

Genome-wide analyses of cassava Pathogenesis-related (PR) gene families reveal core transcriptome responses to whitefly infestation, salicylic acid and jasmonic acid

Background

Whiteflies are a threat to cassava (Manihot esculenta), an important staple food in many tropical/subtropical regions. Understanding the molecular mechanisms regulating cassava’s responses against this pest is crucial for developing control strategies. Pathogenesis-related (PR) protein families are an integral part of plant immunity. With the availability of whole genome sequences, the annotation and expression programs of the full complement of PR genes in an organism can now be achieved. An understanding of the responses of the entire complement of PR genes during biotic stress and to the defense hormones, salicylic acid (SA) and jasmonic acid (JA), is lacking. Here, we analyze the responses of cassava PR genes to whiteflies, SA, JA, and other biotic aggressors.

Results

The cassava genome possesses 14 of the 17 plant PR families, with a total of 447 PR genes. A cassava PR gene nomenclature is proposed. Phylogenetic relatedness of cassava PR proteins to each other and to homologs in poplar, rice and Arabidopsis identified cassava-specific PR gene family expansions. The temporal programs of PR gene expression in response to the whitefly (Aleurotrachelus socialis) in four whitefly-susceptible cassava genotypes showed that 167 of the 447 PR genes were regulated after whitefly infestation. While the timing of PR gene expression varied, over 37% of whitefly-regulated PR genes were downregulated in all four genotypes. Notably, whitefly-responsive PR genes were largely coordinately regulated by SA and JA. The analysis of cassava PR gene expression in response to five other biotic stresses revealed a strong positive correlation between whitefly and Xanthomonas axonopodis and Cassava Brown Streak Virus responses and negative correlations between whitefly and Cassava Mosaic Virus responses. Finally, certain associations between PR genes in cassava expansions and response to biotic stresses were observed among PR families.

Conclusions

This study represents the first genome-wide characterization of PR genes in cassava. PR gene responses to six biotic stresses and to SA and JA are demonstrably different to other angiosperms. We propose that our approach could be applied in other species to fully understand PR gene regulation by pathogens, pests and the canonical defense hormones SA and JA.