A single-step method for RNA isolation from tropical crops in the field

Effect of fieldwork-friendly coffee blender-based extraction methods and leaf tissue storage on the transcriptome of non-model plants

The adaptation of plants to environmental conditions involves a transcriptional response. "Field transcriptomics" is an emerging concept for studying plants in their natural habitat. However, this term includes studies in which cold storage was possible until further processing in a laboratory. Previous studies proposing onsite RNA extraction methods are limited to descriptions of RNA purity, quantity, and quality, and lack a thorough evaluation of transcriptome quality, and transcriptomic evaluations of RNA storage solutions in plants are, to our knowledge, only available for periods of less than a day. This issue is critical for studying plants in geographically difficult-to-access regions, where keeping the cold chain is unrealistic. In this study, the transcriptome of the non-model plant Helonias orientalis (order: Liliales) was evaluated before and after storage of the leaf tissue for one and fourteen days at 25 °C in RNAlater and TRIzol, respectively. Additionally, field-friendly protocols were similarly evaluated for onsite plant RNA extraction at ambient temperature with lightweight equipment that can run on a portable generator, including a guanidine isothiocyanate-free protocol that is compatible with the polyphenol-rich wild strawberry Fragaria vesca. The quality of the transcriptome assembly after 1-day storage and our optimized onsite methods had similar results to that of the state-of-the-art. However, in terms of differential expression analysis, onsite extraction methods performed better overall than the stored tissue samples. We expect that our onsite RNA extraction methods will provide valuable insights into the transcriptional regulation of plants in areas where research equipment is difficult to access.

Bioinformatics Identification and Expression Analysis of Acetyl-CoA Carboxylase Reveal Its Role in Isoflavone Accumulation during Soybean Seed Development

Isoflavones belong to the class of flavonoid compounds, which are important secondary metabolites that play a crucial role in plant development and defense. Acetyl-CoA carboxylase (ACCase) is a biotin-dependent enzyme that catalyzes the conversion of Acetyl-CoA into Malonyl-CoA in plants. It is a key enzyme in fatty acid synthesis and also catalyzes the production of various secondary metabolites. However, information on the ACC gene family in the soybean (Glycine max L. Merr.) genome and the specific members involved in isoflavone biosynthesis is still lacking. In this study, we identified 20 ACC family genes (GmACCs) from the soybean genome and further characterized their evolutionary relationships and expression patterns. Phylogenetic analysis showed that the GmACCs could be divided into five groups, and the gene structures within the same groups were highly conserved, indicating that they had similar functions. The GmACCs were randomly distributed across 12 chromosomes, and collinearity analysis suggested that many GmACCs originated from tandem and segmental duplications, with these genes being under purifying selection. In addition, gene expression pattern analysis indicated that there was functional divergence among GmACCs in different tissues. The GmACCs reached their peak expression levels during the early or middle stages of seed development. Based on the transcriptome and isoflavone content data, a weighted gene co-expression network was constructed, and three candidate genes (Glyma.06G105900, Glyma.13G363500, and Glyma.13G057400) that may positively regulate isoflavone content were identified. These results provide valuable information for the further functional characterization and application of GmACCs in isoflavone biosynthesis in soybean.

Dynamic changes in the plastid and mitochondrial genomes of the angiosperm Corydalis pauciovulata (Papaveraceae)

BACKGROUND

Corydalis DC., the largest genus in the family Papaveraceae, comprises > 465 species. Complete plastid genomes (plastomes) of Corydalis show evolutionary changes, including syntenic arrangements, gene losses and duplications, and IR boundary shifts. However, little is known about the evolution of the mitochondrial genome (mitogenome) in Corydalis. Both the organelle genomes and transcriptomes are needed to better understand the relationships between the patterns of evolution in mitochondrial and plastid genomes.

RESULTS

We obtained complete plastid and mitochondrial genomes from Corydalis pauciovulata using a hybrid assembly of Illumina and Oxford Nanopore Technologies reads to assess the evolutionary parallels between the organelle genomes. The mitogenome and plastome of C. pauciovulata had sizes of 675,483 bp and 185,814 bp, respectively. Three ancestral gene clusters were missing from the mitogenome, and expanded IR (46,060 bp) and miniaturized SSC (202 bp) regions were identified in the plastome. The mitogenome and plastome of C. pauciovulata contained 41 and 67 protein-coding genes, respectively; the loss of genes was a plastid-specific event. We also generated a draft genome and transcriptome for C. pauciovulata. A combination of genomic and transcriptomic data supported the functional replacement of acetyl-CoA carboxylase subunit β (accD) by intracellular transfer to the nucleus in C. pauciovulata. In contrast, our analyses suggested a concurrent loss of the NADH-plastoquinone oxidoreductase (ndh) complex in both the nuclear and plastid genomes. Finally, we performed genomic and transcriptomic analyses to characterize DNA replication, recombination, and repair (DNA-RRR) genes in C. pauciovulata as well as the transcriptomes of Liriodendron tulipifera and Nelumbo nuicifera. We obtained 25 DNA-RRR genes and identified their structure in C. pauciovulata. Pairwise comparisons of nonsynonymous (dN) and synonymous (dS) substitution rates revealed that several DNA-RRR genes in C. pauciovulata have higher dN and dS values than those in N. nuicifera.

CONCLUSIONS

The C. pauciovulata genomic data generated here provide a valuable resource for understanding the evolution of Corydalis organelle genomes. The first mitogenome of Papaveraceae provides an example that can be explored by other researchers sequencing the mitogenomes of related plants. Our results also provide fundamental information about DNA-RRR genes in Corydalis and their related rate variation, which elucidates the relationships between DNA-RRR genes and organelle genome stability.

Intrageneric structural variation in organelle genomes from the genus Dystaenia (Apiaceae): genome rearrangement and mitochondrion-to-plastid DNA transfer

Introduction

During plant evolution, intracellular DNA transfer (IDT) occurs not only from organelles to the nucleus but also between organelles. To further comprehend these events, both organelle genomes and transcriptomes are needed.

Methods

In this study, we constructed organelle genomes and transcriptomes for two Dystaenia species and described their dynamic IDTs between their nuclear and mitochondrial genomes, or plastid and mitochondrial genomes (plastome and mitogenome).

Results and Discussion

We identified the putative functional transfers of the mitochondrial genes 5' rpl2, rps10, rps14, rps19, and sdh3 to the nucleus in both Dystaenia species and detected two transcripts for the rpl2 and sdh3 genes. Additional transcriptomes from the Apicaceae species also provided evidence for the transfers and duplications of these mitochondrial genes, showing lineage-specific patterns. Intrageneric variations of the IDT were found between the Dystaenia organelle genomes. Recurrent plastid-to-mitochondrion DNA transfer events were only identified in the D. takeshimana mitogenome, and a pair of mitochondrial DNAs of plastid origin (MIPTs) may generate minor alternative isoforms. We only found a mitochondrion-to-plastid DNA transfer event in the D. ibukiensis plastome. This event may be linked to inverted repeat boundary shifts in its plastome. We inferred that the insertion region involved an MIPT that had already acquired a plastid sequence in its mitogenome via IDT. We propose that the MIPT acts as a homologous region pairing between the donor and recipient sequences. Our results provide insight into the evolution of organelle genomes across the family Apiaceae.

Comparative Analysis and Identification of Terpene Synthase Genes in Convallaria keiskei Leaf, Flower and Root Using RNA-Sequencing Profiling

The 'Lilly of the Valley' species, Convallaria, is renowned for its fragrant white flowers and distinctive fresh and green floral scent, attributed to a rich composition of volatile organic compounds (VOCs). However, the molecular mechanisms underlying the biosynthesis of this floral scent remain poorly understood due to a lack of transcriptomic data. In this study, we conducted the first comparative transcriptome analysis of C. keiskei, encompassing the leaf, flower, and root tissues. Our aim was to investigate the terpene synthase (TPS) genes and differential gene expression (DEG) patterns associated with essential oil biosynthesis. Through de novo assembly, we generated a substantial number of unigenes, with the highest count in the root (146,550), followed by the flower (116,434) and the leaf (72,044). Among the identified unigenes, we focused on fifteen putative ckTPS genes, which are involved in the synthesis of mono- and sesquiterpenes, the key aromatic compounds responsible for the essential oil biosynthesis in C. keiskei. The expression of these genes was validated using quantitative PCR analysis. Both DEG and qPCR analyses revealed the presence of ckTPS genes in the flower transcriptome, responsible for the synthesis of various compounds such as geraniol, germacrene, kaurene, linalool, nerolidol, trans-ocimene and valencene. The leaf transcriptome exhibited genes related to the biosynthesis of kaurene and trans-ocimene. In the root, the identified unigenes were associated with synthesizing kaurene, trans-ocimene and valencene. Both analyses indicated that the genes involved in mono- and sesquiterpene biosynthesis are more highly expressed in the flower compared to the leaf and root. This comprehensive study provides valuable resources for future investigations aiming to unravel the essential oil-biosynthesis-related genes in the Convallaria genus.

ZmDWF1 regulates leaf angle in maize

Leaf angle (LA) is a critical agronomic trait enhancing grain yield under high-density planting in maize. A number of researches have been conducted in recent years to investigate the quantitative trait loci/genes responsible for LA variation, while only a few genes were identified through map-based cloning. Here we cloned the ZmDWF1 gene, which was previously reported to encode Δ24-sterol reductase in the brassinosteroids (BRs) biosynthesis pathway. Overexpression of ZmDWF1 resulted in enlarged LA, indicating that ZmDWF1 is a positive regulator of LA in maize. To reveal the regulatory framework of ZmDWF1, we conducted RNA-Sequencing and yeast-two hybrid (Y2H) screening analysis. RNA-Sequencing analyzing results indicate ZmDWF1 mainly affected expression level of genes involved in cell wall associated metabolism and hormone metabolism including BR, gibberellin, and auxin. Y2H screening with Bi-FC assay confirmed three proteins (ZmPP2C-1, ZmROF1, and ZmTWD1) interacting with ZmDWF1. We revealed a new regulatory network of ZmDWF1 gene in controlling plant architecture in maize.

An optimized method for RNA extraction from the polyurethane oligomer degrading strain Pseudomonas capeferrum TDA1 growing on aromatic substrates such as phenol and 2,4-diaminotoluene

Bacterial degradation of xenobiotic compounds is an intense field of research already for decades. Lately, this research is complemented by downstream applications including Next Generation Sequencing (NGS), RT-PCR, qPCR, and RNA-seq. For most of these molecular applications, high-quality RNA is a fundamental necessity. However, during the degradation of aromatic substrates, phenolic or polyphenolic compounds such as polycatechols are formed and interact irreversibly with nucleic acids, making RNA extraction from these sources a major challenge. Therefore, we established a method for total RNA extraction from the aromatic degrading Pseudomonas capeferrum TDA1 based on RNAzol® RT, glycogen and a final cleaning step. It yields a high-quality RNA from cells grown on TDA1 and on phenol compared to standard assays conducted in the study. To our knowledge, this is the first report tackling the problem of polyphenolic compound interference with total RNA isolation in bacteria. It might be considered as a guideline to improve total RNA extraction from other bacterial species.

Transcriptome-Wide Identification and Quantification of Caffeoylquinic Acid Biosynthesis Pathway and Prediction of Its Putative BAHDs Gene Complex in A. spathulifolius

The phenylpropanoid pathway is a major secondary metabolite pathway that helps plants overcome biotic and abiotic stress and produces various byproducts that promote human health. Its byproduct caffeoylquinic acid is a soluble phenolic compound present in many angiosperms. Hydroxycinnamate-CoA shikimate/quinate transferase is a significant enzyme that plays a role in accumulating CQA biosynthesis. This study analyzed transcriptome-wide identification of the phenylpropanoid to caffeoylquinic acid biosynthesis candidate genes in A. spathulifolius flowers and leaves. Transcriptomic analyses of the flowers and leaves showed a differential expression of the PPP and CQA biosynthesis regulated unigenes. An analysis of PPP-captive unigenes revealed a major duplication in the following genes: PAL, 120 unigenes in leaves and 76 in flowers; C3′H, 169 unigenes in leaves and 140 in flowers; 4CL, 41 unigenes in leaves and 27 in flowers; and C4H, 12 unigenes in leaves and 4 in flowers. The phylogenetic analysis revealed 82 BAHDs superfamily members in leaves and 72 in flowers, among which five unigenes encode for HQT and three for HCT. The three HQT are common to both leaves and flowers, whereas the two HQT were specialized for leaves. The pattern of HQT synthesis was upregulated in flowers, whereas HCT was expressed strongly in the leaves of A. spathulifolius. Overall, 4CL, C4H, and HQT are expressed strongly in flowers and CAA and HCT show more expression in leaves. As a result, the quantification of HQT and HCT indicates that CQA biosynthesis is more abundant in the flowers and synthesis of caffeic acid in the leaves of A. spathulifolius.

Aster spathulifolius Maxim. a leaf transcriptome provides an overall functional characterization, discovery of SSR marker and phylogeny analysis

Aster spathulifolius Maxim. is belongs to the Asteraceae family, which is distributed only in Korea and Japan. The species is traditionally a medicinal plant and is economically valuable in the ornamental field. On the other hand, the Aster genus, among the Asteraceae family, lacks genomic resources and its molecular functions. Therefore, in our study the high-throughput RNA-sequencing transcriptome data of A. spathulifolius were obtained to identify the molecular functions and its characterization. The de novo assembly produced 98660 uniqueness with an N50 value of 1126bp. Total unigenes were procure to analyze the functional annotation against databases like non-redundant protein, Pfam, Uniprot, KEGG and Gene ontology. The overall percentage of functional annotation to the nr database (43.71%), uniprotein database (49.97%), Pfam (39.94%), KEGG (42.3%) and to GO (30.34%) were observed. Besides, 377 unigenes were found to be involved in the terpenoids pathway and 666 unigenes were actively engaged in other secondary metabolites synthesis, given that 261 unigenes were within phenylpropanoid pathway and 81 unigenes to flavonoid pathway. A further prediction of stress resistance (9,513) unigenes and transcriptional factor (3,027) unigenes in 53 types were vastly regulated in abiotic stress respectively in salt, heat, MAPK and hormone signal transduction pathway. This study discovered 29,692 SSR markers that assist the genotyping approaches and the genetic diversity perspectives. In addition, eight Asteraceae species as in-group together with one out-group were used to construct the phylogenetic relationship by employing their plastid genome and single-copy orthologs genes. Among 50 plastid protein-coding regions, A. spathulifolius is been closely related to A. annua and by 118 single copy orthologs genes, O. taihangensis is more neighboring species to A. spathulifolius. Apart from this, A. spathulifolius and O. taihangensis, genera have recently diverged from other species. Overall, this research gains new insights into transcriptome data by revealing and exposing the secondary metabolite compounds for drug development, the stress-related genes for producing resilient crops and an ortholog gene of A. spathulifolius for the robustness of phylogeny reconstruction among Asteraceae genera.

Identification of a small set of genes commonly regulated in rice roots in response to beneficial rhizobacteria

Rhizosphere bacteria, whether phytopathogenic or phytobeneficial, are thought to be perceived by the plant as a threat. Plant Growth-Promoting Rhizobacteria (PGPR), such as many strains of the Azospirillum genus known as the main phytostimulator of cereals, cooperate with host plants and favorably affect their growth and health. An earlier study of rice root transcriptome, undertaken with two rice cultivars and two Azospirillum strains, revealed a strain-dependent response during the rice-Azospirillum association and showed that only a few genes, including some implicated in plant defense, were commonly regulated in all tested conditions. Here, a set of genes was selected from previous studies and their expression was monitored by qRT-PCR in rice roots inoculated with ten PGPR strains isolated from various plants and belonging to various genera (Azospirillum, Herbaspirillum, Paraburkholderia). A common expression pattern was highlighted for four genes that are proposed to be markers of the rice-PGPR interaction: two genes involved in diterpenoid phytoalexin biosynthesis (OsDXS3 and OsDTC2) and one coding for an uncharacterized protein (Os02g0582900) were significantly induced by PGPR whereas one defense-related gene encoding a pathogenesis-related protein (PR1b, Os01g0382000) was significantly repressed. Interestingly, exposure to a rice bacterial pathogen also triggered the expression of OsDXS3 while the expression of Os02g0582900 and PR1b was down-regulated, suggesting that these genes might play a key role in rice-bacteria interactions. Integration of these results with previous data led us to propose that the jasmonic acid signaling pathway might be triggered in rice roots upon inoculation with PGPR.

A common metabolomic signature is observed upon inoculation of rice roots with various rhizobacteria

Plant growth-promoting rhizobacteria (PGPR), whose growth is stimulated by root exudates, are able to improve plant growth and health. Among those, bacteria of the genus Azospirillum were shown to affect root secondary metabolite content in rice and maize, sometimes without visible effects on root architecture. Transcriptomic studies also revealed that expression of several genes involved in stress and plant defense was affected, albeit with fewer genes when a strain was inoculated onto its original host cultivar. Here, we investigated, via a metabolic profiling approach, whether rice roots responded differently and with gradual intensity to various PGPR, isolated from rice or not. A common metabolomic signature of nine compounds was highlighted, with the reduced accumulation of three alkylresorcinols and increased accumulation of two hydroxycinnamic acid amides (HCAA), identified as N-p-coumaroylputrescine and N-feruloylputrescine. This was accompanied by the increased transcription of two genes involved in the N-feruloylputrescine biosynthetic pathway. Interestingly, exposure to a rice bacterial pathogen triggered a reduced accumulation of these HCAA in roots, a result contrasting with previous reports of increased HCAA content in leaves upon pathogen infection. Accumulation of HCAA, that are potential antimicrobial compounds, might be considered as a primary reaction of plant to bacterial perception.

Full moonlight-induced circadian clock entrainment in Coffea arabica

BACKGROUND

It is now well documented that moonlight affects the life cycle of invertebrates, birds, reptiles, and mammals. The lunisolar tide is also well-known to alter plant growth and development. However, although plants are known to be very photosensitive, few studies have been undertaken to explore the effect of moonlight on plant physiology.

RESULTS

Here for the first time we report a massive transcriptional modification in Coffea arabica genes under full moonlight conditions, particularly at full moon zenith and 3 h later. Among the 3387 deregulated genes found in our study, the main core clock genes were affected.

CONCLUSIONS

Moonlight also negatively influenced many genes involved in photosynthesis, chlorophyll biosynthesis and chloroplast machinery at the end of the night, suggesting that the full moon has a negative effect on primary photosynthetic machinery at dawn. Moreover, full moonlight promotes the transcription of major rhythmic redox genes and many heat shock proteins, suggesting that moonlight is perceived as stress. We confirmed this huge impact of weak light (less than 6 lx) on the transcription of circadian clock genes in controlled conditions mimicking full moonlight.

An Altered Circadian Clock Coupled with a Higher Photosynthesis Efficiency Could Explain the Better Agronomic Performance of a New Coffee Clone When Compared with a Standard Variety

In a context where climate change is threatening coffee productivity, the management of coffee leaf rust is a challenging issue. Major resistant genes, which have been used for many years, are systematically being overcome by pathogens. Developing healthy plants, able to defend themselves and be productive even when attacked by the pathogen, should be part of a more sustainable alternative approach. We compared one hybrid (GPFA124), selected for its good health in various environments including a reduced rust incidence, and the cv. 'Caturra', considered as a standard in terms of productivity and quality but highly susceptible to rust, for phenotypic variables and for the expression of genes involved in the circadian clock and in primary photosynthetic metabolism. The GPFA124 hybrid showed increased photosynthetic electron transport efficiency, better carbon partitioning, and higher chlorophyll content. A strong relationship exists between chlorophyll a fluorescence and the expression of genes related to the photosynthetic electron transport chain. We also showed an alteration of the amplitude of circadian clock genes in the clone. Our work also indicated that increased photosynthetic electron transport efficiency is related to the clone's better performance. Chlorophyll a fluorescence measurement is a good indicator of the coffee tree's physiological status for the breeder. We suggest a connection between the circadian clock and carbon metabolism in coffee tree.

Differential fine-tuning of gene expression regulation in coffee leaves by CcDREB1D promoter haplotypes under water deficit

Despite the importance of the DREB1D gene (also known as CBF4) in plant responses to water deficit and cold stress, studies analysing its regulation by transgenic approaches are lacking. In the current work, a functional study of three CcDREB1D promoter haplotypes (named HP15, HP16 and HP17) isolated from drought-tolerant and drought-sensitive clones of Coffea canephora was carried out in plants of C. arabica stably transformed by Agrobacterium tumefaciens by analysing their ability to regulate the expression of the uidA reporter gene in response to water deficit mimicked by polyethylene glycol (-2.0 MPa) and low relative humidity treatments. A deletion analysis of their corresponding 5'-upstream regions revealed increased specificity of β-glucuronidase activity in the polyethylene glycol and low relative humidity treatments, with high expression in leaf mesophyll and guard cells in full-length constructs. RT-qPCR assays also revealed that the HP16 haplotype (specific to clone tolerant to water deficit) had stronger and earlier activity compared with the HP15 and HP17 haplotypes. As most of the cis-regulatory elements involved in ABA-dependent and -independent networks, tissue specificity and light regulation are common to these haplotypes, we propose that their organization, as well as the nucleic acid polymorphisms present outside these boxes, may play a role in modulating activities of DREB1D promoters in guard cells.