Identification and validation of reference genes for normalization of transcripts from virus-infected Arabidopsis thaliana

Arabidopsis Calcium Dependent Protein Kinase 3, and Its Orthologues OsCPK1, OsCPK15, and AcCPK16, Are Involved in Biotic and Abiotic Stresses

Calcium-dependent protein kinases (CPKs) are plant proteins that directly bind calcium ions before phosphorylating substrates involved in biotic and abiotic stress responses, as well as development. Arabidopsis thaliana CPK3 (AtCPK3) is involved with plant signaling pathways such as stomatal movement regulation, salt stress response, apoptosis, seed germination and pathogen defense. In this study, AtCPK3 and its orthologues in relatively distant plant species such as rice (Oryza sativa, monocot) and kiwifruit (Actinidia chinensis, asterid eudicot) were analyzed in response to drought, bacteria, fungi, and virus infections. Two orthologues were studied in O. sativa, namely OsCPK1 and OsCPK15, while one orthologue-AcCPK16-was identified in A. chinensis. Reverse-transcriptase quantitative PCR (RT-qPCR) analysis revealed that OsCPK1 and AcCPK16 exhibit similar responses to stressors to AtCPK3. OsCPK15 responded differently, particularly in bacterial and fungal infections. An increase in expression was consistently observed among AtCPK3 and its orthologues in response to virus infection. Overexpression mutants in both Arabidopsis and kiwifruit showed slight tolerance to drought, while knockout mutants were slightly more susceptible or had little difference with wild-type plants. Overexpression mutants in Arabidopsis showed slight tolerance to virus infection. These findings highlight the importance of AtCPK3 and its orthologues in drought and pathogen responses and suggest such function must be conserved in its orthologues in a wide range of plants.

Functional characterisation of WRKY transcription factor CrWRKY48 involved in regulating seed abortion of Ponkan (Citrus reticulata)

Citrus fruits are one of the most important fruits in the world, and their seedless character is favored by consumers. WRKY is a plant-specific transcription factor family involved in all aspects of plant growth and development. However, the molecular mechanism of seedless fruit formation in citrus and the role of the WRKY gene family in seed abortion are still poorly understood. In this study, we identified 47 WRKY family genes in the citrus fruit Citrus reticulata and comprehensively characterized the WRKY gene family through gene structure and evolutionary relationships. The expression patterns and protein interaction networks of the WRKY gene family were analyzed based on citrus seed abortion transcriptome data, and several WRKY genes that may be involved in the seed abortion regulation were excavated. Furthermore, CrWRKY48 was verified to regulate seed abortion positively in Arabidopsis thaliana, and the rate of seed abortion caused by overexpression of CrWRKY48 reached 45.48%. Using terminal deoxynucleotidyl transferase-mediated dUTP nick-end labelling assays, DNA affinity purification sequencing and yeast-one-hybrid assays, we found that CrWRKY48 activated excessive programmed cell death by regulating the expression of programmed cell death-related genes such as SOBIR1. Our results show the potential regulation of the WRKY gene family for citrus seed abortion and provide novel insights into the role of CrWRKY48 in mediating citrus seed abortion by activating programmed cell death.

Affecting of Glyphosate Tolerance and Metabolite Content in Transgenic Arabidopsis thaliana Overexpressing EPSPS Gene from Eleusine indica

Long-term use of the global non-selective herbicide glyphosate for weed control has caused resistance in weeds. Overproducing of the target of glyphosate 5-enolpyruvylshikimate-3-phosphate synthase (EPSPS) is one of the resistance mechanisms in weeds. However, few studies have measured the effects on tolerance levels and metabolite content in model plant species overexpressing EPSPS from weeds. We assessed the resistance levels of transgenic Arabidopsis thaliana overexpressing EPSPS from Eleusine indica, and its effects on metabolite content using the method of both quasi-targeted and targeted metabolomics. The results showed that the average resistance index of the transgenic lines was 4.7 and the exogenous E. indica EPSPS expression levels were 265.3- to 532.0-fold higher than those in the wild-type (WT) line. The EPSPS protein ranged from 148.5 to 286.2 μg g-1, which was substantially higher than that in the WT line (9.1 μg g-1). 103 metabolites associated with flavone and flavonol biosynthesis, the metabolism of aromatic amino acids, energy metabolism, and auxin synthesis were significantly higher in the transgenic glyphosate-resistant individuals (R) than in the WT individuals. The results of quantitative analysis show that pyruvate, sedoheptulose 7-phosphate, and gluconic acid amounts in R plants were 1.1-, 1.6- and 1.3-fold higher than those in WT plants, respectively. However, both citric and glyceric acid levels were 0.9-fold lower than those in WT plants. The abundance of other metabolites in the glycolytic and pentose phosphate pathways of central carbon metabolism was similar in the WT and transgenic plants. Glutamic acid was significantly more abundant in the transgenic line than in the WT plants. In contrast, asparagine, glutamine, and lysine were less abundant. However, the concentration of other amino acids did not change significantly. Overexpression of E. indica EPSPS in A. thaliana conferred a moderate level of tolerance to glyphosate. Metabolites associated with flavone and flavonol biosynthesis, the metabolism of aromatic amino acids, and energy metabolism were significantly increased. The results of this study will be useful for evaluating the characterisation and risk assessment of transgenic plants, including identification of unintended effects of the respective transgenic modifications.

Two plant membrane-shaping reticulon-like proteins play contrasting complex roles in turnip mosaic virus infection

Positive-sense RNA viruses remodel cellular cytoplasmic membranes as the membranous sources for the formation of viral replication organelles (VROs) for viral genome replication. In plants, they traffic through plasmodesmata (PD), plasma membrane-lined pores enabling cytoplasmic connections between cells for intercellular movement and systemic infection. In this study, we employed turnip mosaic virus (TuMV), a plant RNA virus to investigate the involvement of RTNLB3 and RTNLB6, two ER (endoplasmic reticulum) membrane-bending, PD-located reticulon-like (RTNL) non-metazoan group B proteins (RTNLBs) in viral infection. We show that RTNLB3 interacts with TuMV 6K2 integral membrane protein and RTNLB6 binds to TuMV coat protein (CP). Knockdown of RTNLB3 promoted viral infection, whereas downregulation of RTNLB6 restricted viral infection, suggesting that these two RTNLs play contrasting roles in TuMV infection. We further demonstrate that RTNLB3 targets the α-helix motif 42LRKSM46 of 6K2 to interrupt 6K2 self-interactions and compromise 6K2-induced VRO formation. Moreover, overexpression of AtRTNLB3 apparently promoted the selective degradation of the ER and ER-associated protein calnexin, but not 6K2. Intriguingly, mutation of the α-helix motif of 6K2 that is required for induction of VROs severely affected 6K2 stability and abolished TuMV infection. Thus, RTNLB3 attenuates TuMV replication, probably through the suppression of 6K2 function. We also show that RTNLB6 promotes viral intercellular movement but does not affect viral replication. Therefore, the proviral role of RTNLB6 is probably by enhancing viral cell-to-cell trafficking. Taken together, our data demonstrate that RTNL family proteins may play diverse complex, even opposite, roles in viral infection in plants.

Calcium-Dependent Protein Kinase GhCDPK16 Exerts a Positive Regulatory Role in Enhancing Drought Tolerance in Cotton

Cotton is essential for the textile industry as a primary source of natural fibers. However, environmental factors like drought present significant challenges to its cultivation, adversely affecting both production levels and fiber quality. Enhancing cotton's drought resilience has the potential to reduce yield losses and support the growth of cotton farming. In this study, the cotton calcium-dependent protein kinase GhCDPK16 was characterized, and the transcription level of GhCDPK16 was significantly upregulated under drought and various stress-related hormone treatments. Physiological analyses revealed that the overexpression of GhCDPK16 improved drought stress resistance in Arabidopsis by enhancing osmotic adjustment capacity and boosting antioxidant enzyme activities. In contrast, silencing GhCDPK16 in cotton resulted in increased dehydration compared with the control. Furthermore, reduced antioxidant enzyme activities and downregulation of ABA-related genes were observed in GhCDPK16-silenced plants. These findings not only enhanced our understanding of the biological functions of GhCDPK16 and the mechanisms underlying drought stress resistance but also underscored the considerable potential of GhCDPK16 in improving drought resilience in cotton.

Recent findings in molecular reactions of E. coli as exposed to alkylated, nano- and ordinary chitosans

The antibacterial effects of chitosan have been widely studied, but the underlying molecular mechanisms are not fully understood. We investigated the molecular responses of Escherichia coli MG1655 cell, a model gram-negative bacterium, upon exposure to chitosan (Cs), alkylated Cs (AlkCs), and chitosan nanoparticles (CsNPs). Nine target genes involved in relevant signaling pathways (ompF, ompC, ompA, mrcA, mrcB, mgtA, glnA, kdpA, lptA) were selected for analysis. A significant reduction in the expression of mrcA, mgtA, glnA, and lptA genes was observed in the cells treated with Cs. Those treated with Cs, AlkCs, and CsNPs revealed an increase in ompF gene expression, but the expression level was lower in the cells treated with AlkCs and CsNPs compared to Cs. This increase in porin expression suggests compromised membrane integrity and disrupted nutrient transport. In addition, the changes in the expression of mgtA, kdpA, and glnA are related to different effects on membrane permeability. The higher expression in the genes mrcA and mrcB is associated with morphological changes of cells treated with AlkCs and CsNPs. These findings contribute to our understanding of the molecular mechanisms underlying chitosan-induced stress responses and provide insights for the development of safer antimicrobial compounds in the future.

Mechanisms of metabolic adaptation in the duckweed Lemna gibba: an integrated metabolic, transcriptomic and flux analysis

BACKGROUND

Duckweeds are small, rapidly growing aquatic flowering plants. Due to their ability for biomass production at high rates they represent promising candidates for biofuel feedstocks. Duckweeds are also excellent model organisms because they can be maintained in well-defined liquid media, usually reproduce asexually, and because genomic resources are becoming increasingly available. To demonstrate the utility of duckweed for integrated metabolic studies, we examined the metabolic adaptation of growing Lemna gibba cultures to different nutritional conditions.

RESULTS

To establish a framework for quantitative metabolic research in duckweeds we derived a central carbon metabolism network model of Lemna gibba based on its draft genome. Lemna gibba fronds were grown with nitrate or glutamine as nitrogen source. The two conditions were compared by quantification of growth kinetics, metabolite levels, transcript abundance, as well as by 13C-metabolic flux analysis. While growing with glutamine, the fronds grew 1.4 times faster and accumulated more protein and less cell wall components compared to plants grown on nitrate. Characterization of photomixotrophic growth by 13C-metabolic flux analysis showed that, under both metabolic growth conditions, the Calvin-Benson-Bassham cycle and the oxidative pentose-phosphate pathway are highly active, creating a futile cycle with net ATP consumption. Depending on the nitrogen source, substantial reorganization of fluxes around the tricarboxylic acid cycle took place, leading to differential formation of the biosynthetic precursors of the Asp and Gln families of proteinogenic amino acids. Despite the substantial reorganization of fluxes around the tricarboxylic acid cycle, flux changes could largely not be associated with changes in transcripts.

CONCLUSIONS

Through integrated analysis of growth rate, biomass composition, metabolite levels, and metabolic flux, we show that Lemna gibba is an excellent system for quantitative metabolic studies in plants. Our study showed that Lemna gibba adjusts to different nitrogen sources by reorganizing central metabolism. The observed disconnect between gene expression regulation and metabolism underscores the importance of metabolic flux analysis as a tool in such studies.

Plant YTHDF proteins are direct effectors of antiviral immunity against an N6-methyladenosine-containing RNA virus

In virus-host interactions, nucleic acid-directed first lines of defense that allow viral clearance without compromising growth are of paramount importance. Plants use the RNA interference pathway as a basal antiviral immune system, but additional RNA-based mechanisms of defense also exist. The infectivity of a plant positive-strand RNA virus, alfalfa mosaic virus (AMV), relies on the demethylation of viral RNA by the recruitment of the cellular N6-methyladenosine (m6 A) demethylase ALKBH9B, but how demethylation of viral RNA promotes AMV infection remains unknown. Here, we show that inactivation of the Arabidopsis cytoplasmic YT521-B homology domain (YTH)-containing m6 A-binding proteins ECT2, ECT3, and ECT5 is sufficient to restore AMV infectivity in partially resistant alkbh9b mutants. We further show that the antiviral function of ECT2 is distinct from its previously demonstrated function in the promotion of primordial cell proliferation: an ect2 mutant carrying a small deletion in its intrinsically disordered region is partially compromised for antiviral defense but not for developmental functions. These results indicate that the m6 A-YTHDF axis constitutes a novel branch of basal antiviral immunity in plants.

Laser capture microdissection-based spatiotemporal transcriptomes uncover regulatory networks during seed abortion in seedless Ponkan (Citrus reticulata)

Seed abortion is an important process in the formation of seedless characteristics in citrus fruits. However, the molecular regulatory mechanism underlying citrus seed abortion is poorly understood. Laser capture microdissection-based RNA-seq combined with Pacbio-seq was used to profile seed development in the Ponkan cultivars 'Huagan No. 4' (seedless Ponkan) (Citrus reticulata) and 'E'gan No. 1' (seeded Ponkan) (C. reticulata) in two types of seed tissue across three developmental stages. Through comparative transcriptome and dynamic phytohormone analyses, plant hormone signal, cell division and nutrient metabolism-related processes were revealed to play critical roles in the seed abortion of 'Huagan No. 4'. Moreover, several genes may play indispensable roles in seed abortion of 'Huagan No. 4', such as CrWRKY74, CrWRKY48 and CrMYB3R4. Overexpression of CrWRKY74 in Arabidopsis resulted in severe seed abortion. By analyzing the downstream regulatory network, we further determined that CrWRKY74 participated in seed abortion regulation by inducing abnormal programmed cell death. Of particular importance is that a preliminary model was proposed to depict the regulatory networks underlying seed abortion in citrus. The results of this study provide novel insights into the molecular mechanism across citrus seed development, and reveal the master role of CrWRKY74 in seed abortion of 'Huagan No. 4'.

Differential Participation of Plant Ribosomal Proteins from the Small Ribosomal Subunit in Protein Translation under Stress

Upon exposure to biotic and abiotic stress, plants have developed strategies to adapt to the challenges imposed by these unfavorable conditions. The energetically demanding translation process is one of the main elements regulated to reduce energy consumption and to selectively synthesize proteins involved in the establishment of an adequate response. Emerging data have shown that ribosomes remodel to adapt to stresses. In Arabidopsis thaliana, ribosomes consist of approximately eighty-one distinct ribosomal proteins (RPs), each of which is encoded by two to seven genes. Recent research has revealed that a mutation in a given single RP in plants can not only affect the functions of the RP itself but can also influence the properties of the ribosome, which could bring about changes in the translation to varying degrees. However, a pending question is whether some RPs enable ribosomes to preferentially translate specific mRNAs. To reveal the role of ribosomal proteins from the small subunit (RPS) in a specific translation, we developed a novel approach to visualize the effect of RPS silencing on the translation of a reporter mRNA (GFP) combined to the 5'UTR of different housekeeping and defense genes. The silencing of genes encoding for NbRPSaA, NbRPS5A, and NbRPS24A in Nicotiana benthamiana decreased the translation of defense genes. The NbRACK1A-silenced plant showed compromised translations of specific antioxidant enzymes. However, the translations of all tested genes were affected in NbRPS27D-silenced plants. These findings suggest that some RPS may be potentially involved in the control of protein translation.

Selection and Validation of Reference Genes in Virus-Infected Sweet Potato Plants

Quantitative real-time PCR (qRT-PCR) in sweet potatoes requires accurate data normalization; however, there are insufficient studies on appropriate reference genes for gene expression analysis. We examined variations in the expression of eight candidate reference genes in the leaf and root tissues of sweet potatoes (eight nonvirus-infected or eight virus-infected samples). Parallel analyses with geNorm, NormFinder, and Best-Keeper show that different viral infections and origin tissues affect the expression levels of these genes. Based on the results of the evaluation of the three software, the adenosine diphosphate-ribosylation factor is suitable for nonvirus or virus-infected sweet potato leaves. Cyclophilin and ubiquitin extension proteins are suitable for nonvirus-infected sweet potato leaves. Phospholipase D1 alpha is suitable for virus-infected sweet potato leaves. Actin is suitable for roots of nonvirus-infected sweet potatoes. Glyceraldehyde-3-phosphate dehydrogenase is suitable for virus-infected sweet potato roots. The research provides appropriate reference genes for further analysis in leaf and root samples of viruses in sweet potatoes.

An integrase toolbox to record gene-expression during plant development

There are many open questions about the mechanisms that coordinate the dynamic, multicellular behaviors required for organogenesis. Synthetic circuits that can record in vivo signaling networks have been critical in elucidating animal development. Here, we report on the transfer of this technology to plants using orthogonal serine integrases to mediate site-specific and irreversible DNA recombination visualized by switching between fluorescent reporters. When combined with promoters expressed during lateral root initiation, integrases amplify reporter signal and permanently mark all descendants. In addition, we present a suite of methods to tune the threshold for integrase switching, including: RNA/protein degradation tags, a nuclear localization signal, and a split-intein system. These tools improve the robustness of integrase-mediated switching with different promoters and the stability of switching behavior over multiple generations. Although each promoter requires tuning for optimal performance, this integrase toolbox can be used to build history-dependent circuits to decode the order of expression during organogenesis in many contexts.

Selection and Validation of Reference Genes for RT-qPCR Analysis of Gene Expression in Nicotiana benthamiana upon Single Infections by 11 Positive-Sense Single-Stranded RNA Viruses from Four Genera

Quantitative real-time PCR (RT-qPCR) is a widely used method for studying alterations in gene expression upon infections caused by diverse pathogens such as viruses. Positive-sense single-stranded (ss(+)) RNA viruses form a major part of all known plant viruses, and some of them are damaging pathogens of agriculturally important crops. Analysis of gene expression following infection by ss(+) RNA viruses is crucial for the identification of potential anti-viral factors. However, viral infections are known to globally affect gene expression and therefore selection and validation of reference genes for RT-qPCR is particularly important. In this study, the expression of commonly used reference genes for RT-qPCR was studied in Nicotiana benthamiana following single infection by 11 ss(+) RNA viruses, including five tobamoviruses, four potyviruses, one potexvirus and one polerovirus. Stability of gene expression was analyzed in parallel by four commonly used algorithms: geNorm, NormFinder, BestKeeper, and Delta CT, and RefFinder was finally used to summarize all the data. The most stably expressed reference genes differed significantly among the viruses, even when those viruses were from the same genus. Our study highlights the importance of the selection and validation of reference genes upon different viral infections.

Ecotype-specific blockage of tasiARF production by two different RNA viruses in Arabidopsis

Arabidopsis thaliana is one of the most studied model organisms of plant biology with hundreds of geographical variants called ecotypes. One might expect that this enormous genetic variety could result in differential response to pathogens. Indeed, we observed previously that the Bur ecotype develops much more severe symptoms (upward curling leaves and wavy leaf margins) upon infection with two positive-strand RNA viruses of different families (turnip vein-clearing virus, TVCV, and turnip mosaic virus, TuMV). To find the genes potentially responsible for the ecotype-specific response, we performed a differential expression analysis of the mRNA and sRNA pools of TVCV and TuMV-infected Bur and Col plants along with the corresponding mock controls. We focused on the genes and sRNAs that showed an induced or reduced expression selectively in the Bur virus samples in both virus series. We found that the two ecotypes respond to the viral infection differently, yet both viruses selectively block the production of the TAS3-derived small RNA specimen called tasiARF only in the virus-infected Bur plants. The tasiARF normally forms a gradient through the adaxial and abaxial parts of the leaf (being more abundant in the adaxial part) and post-transcriptionally regulates ARF4, a major leaf polarity determinant in plants. The lack of tasiARF-mediated silencing could lead to an ectopically expressed ARF4 in the adaxial part of the leaf where the misregulation of auxin-dependent signaling would result in an irregular growth of the leaf blade manifesting as upward curling leaf and wavy leaf margin. QTL mapping using Recombinant Inbred Lines (RILs) suggests that the observed symptoms are the result of a multigenic interaction that allows the symptoms to develop only in the Bur ecotype. The particular nature of genetic differences leading to the ecotype-specific symptoms remains obscure and needs further study.

AtGSTU19 and AtGSTU24 as Moderators of the Response of Arabidopsis thaliana to Turnip mosaic virus

Plants produce glutathione as a response to the intercellular redox state. Glutathione actively participates in the reactive oxygen species (ROS)-dependent signaling pathway, especially under biotic stress conditions. Most of the glutathione S-transferases (GSTs) are induced in cells during the defense response of plants not only through highly specific glutathione-binding abilities but also by participating in the signaling function. The tau class of GSTs has been reported to be induced as a response under stress conditions. Although several studies have focused on the role of the tau class of GSTs in plant–pathogen interactions, knowledge about their contribution to the response to virus inoculation is still inadequate. Therefore, in this study, the response of Atgstu19 and Atgstu24 knockout mutants to mechanical inoculation of Turnip mosaic virus (TuMV) was examined. The systemic infection of TuMV was more dynamically promoted in Atgstu19 mutants than in wild-type (Col-0) plants, suggesting the role of GSTU19 in TuMV resistance. However, Atgstu24 mutants displayed virus limitation and downregulation of the relative expression of TuMV capsid protein, accompanied rarely by TuMV particles only in vacuoles, and ultrastructural analyses of inoculated leaves revealed the lack of virus cytoplasmic inclusions. These findings indicated that Atgstu24 mutants displayed a resistance-like reaction to TuMV, suggesting that GSTU24 may suppress the plant resistance. In addition, these findings confirmed that GSTU1 and GSTU24 are induced and contribute to the susceptible reaction to TuMV in the Atgstu19–TuMV interaction. However, the upregulation of GSTU19 and GSTU13 highly correlated with virus limitation in the resistance-like reaction in the Atgstu24–TuMV interaction. Furthermore, the highly dynamic upregulation of GST and glutathione reductase (GR) activities resulted in significant induction (between 1 and 14 days post inoculation [dpi]) of the total glutathione pool (GSH + GSSG) in response to TuMV, which was accompanied by the distribution of active glutathione in plant cells. On the contrary, in Atgstu19, which is susceptible to TuMV interaction, upregulation of GST and GR activity only up to 7 dpi symptom development was reported, which resulted in the induction of the total glutathione pool between 1 and 3 dpi. These observations indicated that GSTU19 and GSTU24 are important factors in modulating the response to TuMV in Arabidopsis thaliana. Moreover, it was clear that glutathione is an important component of the regulatory network in resistance and susceptible response of A. thaliana to TuMV. These results help achieve a better understanding of the mechanisms regulating the Arabidopsis–TuMV pathosystem.

Selection and Validation of Reference Genes for Gene Expression Studies Using Quantitative Real-Time PCR in Prunus Necrotic Ringspot Virus-Infected Cucumis sativus

Several members of the genus Ilarvirus infect fruit trees and are distributed worldwide. Prunus necrotic ringspot virus (PNRSV) is one of the most prevalent viruses, causing significant losses. Cucumissativus can be infected by several ilarviruses, leading to obvious symptoms, including PNRSV, which suggests that cucumbers could be good hosts for the study of the pathogenesis of ilarviruses. Real-time quantitative PCR is an optimal choice for studying gene expression because of its simplicity and its fast and high sensitivity, while its accuracy is highly dependent on the stability of the reference genes. In this study, we assessed the stability of eleven reference genes with geNorm, NormFinder, ΔCt method, BestKeeper, and the ranking software, RefFinder. The results indicated that the combined use of EF1α and F-BOX was the most accurate normalization method. In addition, the host genes AGO1, AGO4, and RDR6 were selected to test the reliability of the reference genes. This study provides useful information for gene expression analysis during PNRSV infection and will facilitate gene expression studies associated with ilarvirus infection.

Identification of positive and negative regulators of antiviral RNA interference in Arabidopsis thaliana

Virus-host coevolution often drives virus immune escape. However, it remains unknown whether natural variations of plant virus resistance are enriched in genes of RNA interference (RNAi) pathway known to confer essential antiviral defense in plants. Here, we report two genome-wide association study screens to interrogate natural variation among wild-collected Arabidopsis thaliana accessions in quantitative resistance to the endemic cucumber mosaic virus (CMV). We demonstrate that the highest-ranked gene significantly associated with resistance from both screens acts to regulate antiviral RNAi in ecotype Columbia-0. One gene, corresponding to Reduced Dormancy 5 (RDO5), enhances resistance by promoting amplification of the virus-derived small interfering RNAs (vsiRNAs). Interestingly, the second gene, designated Antiviral RNAi Regulator 1 (VIR1), dampens antiviral RNAi so its genetic inactivation by CRISPR/Cas9 editing enhances both vsiRNA production and CMV resistance. Our findings identify positive and negative regulators of the antiviral RNAi defense that may play important roles in virus-host coevolution.

Selection and Validation of Reference Genes for RT-qPCR Analysis in Tibetan Medicinal Plant Saussurea Laniceps Callus under Abiotic Stresses and Hormone Treatments

Real-time quantitative PCR (RT-qPCR) is an important technique for studying gene expression analysis, but accurate and reliable results depend on the use of a stable reference gene. This study proposes to test the expression stability of candidate reference genes in the callus of Saussurea laniceps, a unique Tibetan medicinal plant. Based on the S. laniceps callus transcriptome, eleven candidate reference genes, including TUA2, TUB3, TUB8, TIF3B1, TIF3H1, ELF5A, PP2AA2, UEV1D, UBL5, UBC36, and SKIP1), were validated for RT-qPCR normalization in the callus under abiotic stress (salt, cold, and UV) and hormone treatments (abscisic acid, MeJA, and salicylic acid). The stability of the candidate genes was evaluated in all the samples of S. laniceps. Comprehensive analysis of all samples showed that the best reference genes were UBC36 and UBL5. ELF5A and TIF3B1 were ranked as the most stable genes in the sample sets under abiotic stress. For hormone stimulation, UBC36 and TIF3H1 genes had the best stability. This study provides useful guidelines and a starting point for reference gene selection for expression analysis using RT-qPCR techniques in S. laniceps.

The Potential of Molecular Indicators of Plant Virus Infection: Are Plants Able to Tell Us They Are Infected?

To our knowledge, there are no reports that demonstrate the use of host molecular markers for the purpose of detecting generic plant virus infection. Two approaches involving molecular indicators of virus infection in the model plant Arabidopsis thaliana were examined: the accumulation of small RNAs (sRNAs) using a microfluidics-based method (Bioanalyzer); and the transcript accumulation of virus-response related host plant genes, suppressor of gene silencing 3 (AtSGS3) and calcium-dependent protein kinase 3 (AtCPK3) by reverse transcriptase-quantitative PCR (RT-qPCR). The microfluidics approach using sRNA chips has previously demonstrated good linearity and good reproducibility, both within and between chips. Good limits of detection have been demonstrated from two-fold 10-point serial dilution regression to 0.1 ng of RNA. The ratio of small RNA (sRNA) to ribosomal RNA (rRNA), as a proportion of averaged mock-inoculation, correlated with known virus infection to a high degree of certainty. AtSGS3 transcript decreased between 14- and 28-days post inoculation (dpi) for all viruses investigated, while AtCPK3 transcript increased between 14 and 28 dpi for all viruses. A combination of these two molecular approaches may be useful for assessment of virus-infection of samples without the need for diagnosis of specific virus infection.

Studying the Function of Phytoplasma Effector Proteins Using a Chemical-Inducible Expression System in Transgenic Plants

Phytoplasmas are bacterial pathogens that live mainly in the phloem of their plant hosts. They dramatically manipulate plant development by secreting effector proteins that target developmental proteins of their hosts. Traditionally, the effects of individual effector proteins have been studied by ectopic overexpression using strong, ubiquitously active promoters in transgenic model plants. However, the impact of phytoplasma infection on the host plants depends on the intensity and timing of infection with respect to the developmental stage of the host. To facilitate investigations addressing the timing of effector protein activity, we have established chemical-inducible expression systems for the three most well-characterized phytoplasma effector proteins, SECRETED ASTER YELLOWS WITCHES’ BROOM PROTEIN 11 (SAP11), SAP54 and TENGU in transgenic Arabidopsis thaliana. We induced gene expression either continuously, or at germination stage, seedling stage, or flowering stage. mRNA expression was determined by quantitative reverse transcription PCR, protein accumulation by confocal laser scanning microscopy of GFP fusion proteins. Our data reveal tight regulation of effector gene expression and strong upregulation after induction. Phenotypic analyses showed differences in disease phenotypes depending on the timing of induction. Comparative phenotype analysis revealed so far unreported similarities in disease phenotypes, with all three effector proteins interfering with flower development and shoot branching, indicating a surprising functional redundancy of SAP54, SAP11 and TENGU. However, subtle but mechanistically important differences were also observed, especially affecting the branching pattern of the plants.

Arabidopsis pavement cell morphogenesis requires FERONIA binding to pectin for activation of ROP GTPase signaling

Sensing and signaling of cell wall status and dynamics regulate many processes in plants, such as cell growth and morphogenesis, but the underpinning mechanisms remain largely unknown. Here, we demonstrate that the CrRLK1L receptor kinase FERONIA (FER) binds the cell wall pectin, directly leading to the activation of the ROP6 guanosine triphosphatase (GTPase) signaling pathway that regulates the formation of the puzzle piece shape of pavement cells in Arabidopsis. The extracellular malectin domain of FER binds demethylesterified pectin in vivo and in vitro. Both loss-of-FER mutations and defects in pectin demethylesterification caused similar changes in pavement cell shape and ROP6 GTPase signaling. FER is required for the activation of ROP6 by demethylesterified pectin and physically and genetically interacts with the ROP6 activator, RopGEF14. Thus, our findings elucidate a signaling pathway that directly connects the cell wall pectin to cellular morphogenesis via the cell surface receptor FER.

Severe Stunting Symptoms upon Nepovirus Infection Are Reminiscent of a Chronic Hypersensitive-like Response in a Perennial Woody Fruit Crop

Virus infection of plants can result in various degrees of detrimental impacts and disparate symptom types and severities. Although great strides have been made in our understanding of the virus-host interactions in herbaceous model plants, the mechanisms underlying symptom development are poorly understood in perennial fruit crops. Grapevine fanleaf virus (GFLV) causes variable symptoms in most vineyards worldwide. To better understand GFLV-grapevine interactions in relation to symptom development, field and greenhouse trials were conducted with a grapevine genotype that exhibits distinct symptoms in response to a severe and a mild strain of GFLV. After validation of the infection status of the experimental vines by high-throughput sequencing, the transcriptomic and metabolomic profiles in plants infected with the two viral strains were tested and compared by RNA-Seq and LC-MS, respectively, in the differentiating grapevine genotype. In vines infected with the severe GFLV strain, 1023 genes, among which some are implicated in the regulation of the hypersensitive-type response, were specifically deregulated, and a higher accumulation of resveratrol and phytohormones was observed. Interestingly, some experimental vines restricted the virus to the rootstock and remained symptomless. Our results suggest that GFLV induces a strain- and cultivar-specific defense reaction similar to a hypersensitive reaction. This type of defense leads to a severe stunting phenotype in some grapevines, whereas others are resistant. This work is the first evidence of a hypersensitive-like reaction in grapevine during virus infection.

GTP binding by Arabidopsis extra-large G protein 2 is not essential for its functions

The extra-large guanosine-5'-triphosphate (GTP)-binding protein 2, XLG2, is an unconventional Gα subunit of the Arabidopsis (Arabidopsis thaliana) heterotrimeric GTP-binding protein complex with a major role in plant defense. In vitro biochemical analyses and molecular dynamic simulations show that affinity of XLG2 for GTP is two orders of magnitude lower than that of the conventional Gα, AtGPA1. Here we tested the physiological relevance of GTP binding by XLG2. We generated an XLG2(T476N) variant with abolished GTP binding, as confirmed by in vitro GTPγS binding assay. Yeast three-hybrid, bimolecular fluorescence complementation, and split firefly-luciferase complementation assays revealed that the nucleotide-depleted XLG2(T476N) retained wild-type XLG2-like interactions with the Gβγ dimer and defense-related receptor-like kinases. Both wild-type and nucleotide-depleted XLG2(T476N) restored the defense responses against Fusarium oxysporum and Pseudomonas syringae compromised in the xlg2 xlg3 double mutant. Additionally, XLG2(T476N) was fully functional restoring stomatal density, root growth, and sensitivity to NaCl, but failed to complement impaired germination and vernalization-induced flowering. We conclude that XLG2 is able to function in a GTP-independent manner and discuss its possible mechanisms of action.

Arabidopsis Ca2+-ATPases 1, 2, and 7 in the endoplasmic reticulum contribute to growth and pollen fitness

Generating cellular Ca2+ signals requires coordinated transport activities from both Ca2+ influx and efflux pathways. In Arabidopsis (Arabidopsis thaliana), multiple efflux pathways exist, some of which involve Ca2+-pumps belonging to the Autoinhibited Ca2+-ATPase (ACA) family. Here, we show that ACA1, 2, and 7 localize to the endoplasmic reticulum (ER) and are important for plant growth and pollen fertility. While phenotypes for plants harboring single-gene knockouts (KOs) were weak or undetected, a triple KO of aca1/2/7 displayed a 2.6-fold decrease in pollen transmission efficiency, whereas inheritance through female gametes was normal. The triple KO also resulted in smaller rosettes showing a high frequency of lesions. Both vegetative and reproductive phenotypes were rescued by transgenes encoding either ACA1, 2, or 7, suggesting that all three isoforms are biochemically redundant. Lesions were suppressed by expression of a transgene encoding NahG, an enzyme that degrades salicylic acid (SA). Triple KO mutants showed elevated mRNA expression for two SA-inducible marker genes, Pathogenesis-related1 (PR1) and PR2. The aca1/2/7 lesion phenotype was similar but less severe than SA-dependent lesions associated with a double KO of vacuolar pumps aca4 and 11. Imaging of Ca2+ dynamics triggered by blue light or the pathogen elicitor flg22 revealed that aca1/2/7 mutants display Ca2+ transients with increased magnitudes and durations. Together, these results indicate that ER-localized ACAs play important roles in regulating Ca2+ signals, and that the loss of these pumps results in male fertility and vegetative growth deficiencies.

Reference gene selection for miRNA and mRNA normalization in potato in response to potato virus Y

This was the first report on evaluating candidate reference genes for quantifying the expression profiles of both coding (e.g., mRNA) and non-coding (e.g., miRNA) genes in potato response to potato virus Y (PVY) inoculation. The reverse transcription-quantitative real-time polymerase chain reaction (RT-qPCR) method was employed to quantify the expression profiles of eight selected candidate reference genes; their expression stability was analyzed by four statistical algorithms, i.e., geNorm, BestKeeper, NormFinder and RefFinder. The most stable reference genes were sEF1a, sTUBb and seIF5 with a high stability. The least stable ones were sPP2A, sSUI1 and sGAPDH. The same reference gene allows for normalization of both miRNA and mRNA levels from a single RNA sample using cDNAs synthesized in a single RT reaction, in which a stem-loop primer was used for miRNAs and the oligo (dT) for mRNAs.

Identification, characterization, and expression analysis of cowpea (Vigna unguiculata [L.] Walp.) miRNAs in response to cowpea severe mosaic virus (CPSMV) challenge

Cowpea miRNAs and Argonaute genes showed differential expression patterns in response to CPSMV challenge Several biotic stresses affect cowpea production and yield. CPSMV stands out for causing severe negative impacts on cowpea. Plants have two main induced immune systems. In the basal system (PTI, PAMP-triggered immunity), plants recognize and respond to conserved molecular patterns associated with pathogens (PAMPs). The second type (ETI, Effector-triggered immunity) is induced after plant recognition of specific factors from pathogens. RNA silencing is another important defense mechanism in plants. Our research group has been using biochemical and proteomic approaches to learn which proteins and pathways are involved and could explain why some cowpea genotypes are resistant whereas others are susceptible to CPSMV. This current study was conducted to determine the role of cowpea miRNA in the interaction between a resistant cowpea genotype (BRS-Marataoã) and CPSMV. Previously identified and deposited plant microRNA sequences were used to find out all possible microRNAs in the cowpea genome. This search detected 617 mature microRNAs, which were distributed in 89 microRNA families. Next, 4 out of these 617 miRNAs and their possible target genes that encode the proteins Kat-p80, DEAD-Box, GST, and SPB9, all involved in the defense response of cowpea to CPSMV, had their expression compared between cowpea leaves uninoculated and inoculated with CPSMV. Additionally, the differential expression of genes that encode the Argonaute (AGO) proteins 1, 2, 4, 6, and 10 is reported. In summary, the studied miRNAs and AGO 2 and AGO4 associated genes showed differential expression patterns in response to CPSMV challenge, which indicate their role in cowpea defense.

Identification and validation of reference genes for real-time RT-PCR in Aphelenchoides besseyi

Fragments of four candidate reference genes of Aphelenchoides besseyi, including actin, glyceraldehyde-3-phosphate dehydrogenase (GAPDH), ubiquitin conjugating-3 enzyme (UBC) and alpha-tubulin (α-tubulin) were cloned from the transcriptome database of A. besseyi. The expression level of these four candidate reference genes and a commonly used reference gene of A. besseyi (18S rRNA) in three experimental conditions, including the four life stages (female, male, juvenile and egg) of two populations and the mixed-stage nematodes of four populations with different origins and hosts were analyzed by RT-qPCR. The expression stability of the five candidate reference genes under the three experimental conditions was analyzed by ΔCt, geNorm, NormFinder and RefFinder respectively. The analysis results of ΔCt, geNorm, NormFinder and RefFinder all indicated that UBC was the gene with the highest average ranking of stability. In conclusion, the expression stability of UBC was optimal under the three experimental conditions, indicating that UBC could be used as a suitable reference gene instead of 18S rRNA in the RT-qPCR analysis for A. besseyi.

Selection and Validation of Reference Genes for Normalisation of Gene Expression in Glehnia littoralis

Glehnia littoralis is an important medicinal halophyte—the dried root of which is used as Chinese herbal medicine. However, the use, selection and stability of reference genes are rarely verified in studies of G. littoralis, which hampers investigation of its salt tolerance and metabolism. In this study, we selected 13 candidate reference genes from the transcriptome data of G. littoralis—serine/threonine-protein phosphatase PP2A (PP2A), polyubiquitin 10 (UBQ10), actin (ACT), elongation factor 1-α (EF1-α), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), α-tubulin (α-TUB), β-tubulin (β-TUB), polypyrimidine tract-binding protein 1 (PTBP1), expressed protein 1 (EXP1), expressed protein 2 (EXP2), TIP41-like (TIP41), SAND family (SAND), and cyclophilin 2 (CYP2), and used qRT-PCR to analyse their expression levels in roots of G. littoralis treated with NaCl, polyethylene glycol (PEG), abscisic acid (ABA), and methyl jasmonate (MeJA), as well as in various organs of G. littoralis. The ΔCt, geNorm, NormFinder, and BestKeeper algorithms were used to assess the expression stability of the candidate reference genes and the results were then used to generate a comprehensive rank list with the RankAggreg R package. The most stable reference genes for normalisation were EXP1 and PP2A in response to NaCl, EXP2 and PP2A in response to ABA, CYP2 and α-TUB in response to MeJA, and ACT and EXP1 in the PEG and the organ subsets. GAPDH, β-TUB, and UBQ10 exhibited low stability and so were unsuitable for normalisation. This study is the first systematic analysis of candidate reference genes in G. littoralis and will facilitate further investigation of normalisation of gene expression in G. littoralis.

A conserved motif in three viral movement proteins from different genera is required for host factor recruitment and cell-to-cell movement

Due to their minimal genomes, plant viruses are forced to hijack specific cellular pathways to ensure host colonization, a condition that most frequently involves physical interaction between viral and host proteins. Among putative viral interactors are the movement proteins, responsible for plasmodesma gating and genome binding during viral transport. Two of them, DGBp1 and DGBp2, are required for alpha-, beta- and gammacarmovirus cell-to-cell movement, but the number of DGBp-host interactors identified at present is limited. By using two different approaches, yeast two-hybrid and bimolecular fluorescence complementation assays, we found three Arabidopsis factors, eIF3g1, RPP3A and WRKY36, interacting with DGBp1s from each genus mentioned above. eIF3g1 and RPP3A are mainly involved in protein translation initiation and elongation phases, respectively, while WRKY36 belongs to WRKY transcription factor family, important regulators of many defence responses. These host proteins are not expected to be associated with viral movement, but knocking out WRKY36 or silencing either RPP3A or eIF3g1 negatively affected Arabidopsis infection by Turnip crinkle virus. A highly conserved FNF motif at DGBp1 C-terminus was required for protein-protein interaction and cell-to-cell movement, suggesting an important biological role.

Nucleotide exchange-dependent and nucleotide exchange-independent functions of plant heterotrimeric GTP-binding proteins

Heterotrimeric guanine nucleotide-binding proteins (G proteins), which are composed of α, β, and γ subunits, are versatile, guanine nucleotide-dependent, molecular on-off switches. In animals and fungi, the exchange of GDP for GTP on Gα controls G protein activation and is crucial for normal cellular responses to diverse extracellular signals. The model plant Arabidopsis thaliana has a single canonical Gα subunit, AtGPA1. We found that, in planta, the constitutively active, GTP-bound AtGPA1(Q222L) mutant and the nucleotide-free AtGPA1(S52C) mutant interacted with Gβγ1 and Gβγ2 dimers with similar affinities, suggesting that G protein heterotrimer formation occurred independently of nucleotide exchange. In contrast, AtGPA1(Q222L) had a greater affinity than that of AtGPA1(S52C) for Gβγ3, suggesting that the GTP-bound conformation of AtGPA1(Q222L) is distinct and tightly associated with Gβγ3. Functional analysis of transgenic lines expressing either AtGPA1(S52C) or AtGPA1(Q222L) in the gpa1-null mutant background revealed various mutant phenotypes that were complemented by either AtGPA1(S52C) or AtGPA1(Q222L). We conclude that, in addition to the canonical GDP-GTP exchange-dependent mechanism, plant G proteins can function independently of nucleotide exchange.

Interactions between soybean, Bradyrhizobium japonicum and Soybean mosaic virus: the effects depend on the interaction sequence

The symbiotic interaction between soybean and nitrogen-fixing rhizobia can lead to plant growth promotion and induced systemic responses. Symbiotic interactions may increase tolerance/resistance to abiotic/biotic stress conditions, but are also sensitive to environmental conditions. Soybean mosaic virus (SMV), which is transmitted by seed and aphids, severely affects crop yields in many areas of the world, consequently virus infection may precede rhizobium infection or vice versa in the field. With the hypothesis that sequence of interaction is a key determinant of the resulting responses; growth, primary metabolism and defence responses were evaluated in different interaction sequences. Results showed that vegetative growth was promoted by Bradyrhizobium japonicum (Bj) inoculation and drastically impaired by SMV infection. The negative effect of SMV single infection on soybean growth parameters was correlated with photosynthesis decrease, sugar accumulation, oxidative damage, and increases in salicylic acid levels. Bj inoculation partially reversed virus-induced symptoms, mainly at Bj-SMV sequence. However, this symptom attenuation did not correlate with less virus accumulation. Nodulation was negatively affected by SMV, particularly when virus infection was previous to Bj inoculation (SMV-Bj). Defence related hormones (salicylic acid (SA)/jasmonic acid (JA)) and the expression of defence-related genes were dependent on the sequence of tripartite interaction. The present study showed that the sequence of the tripartite interaction among soybean, Bj and SMV determinates the tolerance/susceptibility to SMV infection, through changes in the defence mechanism and metabolic alteration.

Involvement of OpsLTP1 from Opuntia streptacantha in abiotic stress adaptation and lipid metabolism

Plant lipid transfer proteins (LTPs) exhibit the ability to transfer lipids between membranes in vitro, and have been implicated in diverse physiological processes associated to plant growth, reproduction, development, biotic and abiotic stress responses. However, their mode of action is not yet fully understood. To explore the functions of the OpsLTP1 gene encoding a LTP from cactus pear Opuntia streptacantha Lem., we generated transgenic Arabidopsis thaliana (L.) Heynh. plants to overexpress OpsLTP1 and contrasted our results with the loss-of-function mutant ltp3 from A. thaliana under abiotic stress conditions. The ltp3 mutant seeds showed impaired germination under salt and osmotic treatments, in contrast to OpsLTP1 overexpressing lines that displayed significant increases in germination rate. Moreover, stress recovery assays showed that ltp3 mutant seedlings were more sensitive to salt and osmotic treatments than wild-type plants suggesting that AtLTP3 is required for stress-induced responses, while the OpsLTP1 overexpressing line showed no significant differences. In addition, OpsLTP1 overexpressing and ltp3 mutant seeds stored lower amount of total lipids compared with wild-type seeds, showing changes primarily on 16C and 18C fatty acids. However, ltp3 mutant also lead changes in lipid profile and no over concrete lipids which may suggest a compensatory activation of other LTPs. Interestingly, linoleic acid (18:2ω6) was consistently increased in neutral, galactoglycerolipids and phosphoglycerolipids of OpsLTP1 overexpressing line indicating a role of OpsLTP1 in the modulation of lipid composition in A. thaliana.

Selection of suitable reference genes for quantitive real-time PCR normalization in Miscanthus lutarioriparia

Miscanthus lutarioriparia, which is found widespread in China, has attracted great attention as a most potential bioenergy plant for years. The quantitative real time PCR (RT-qPCR) has appeared as a sensitive and powerful technique to measure gene expression in living organisms during different development stages. In this study, we evaluated ten candidate genes, including 25S ribosomal RNA gene (25S rRNA), actin1 gene (ACT1), carotenoid-binding protein 20 gene (CBP20), glyceraldehyde-3-phosphate dehydrogenase gene (GAPDH), Ubiquitin gene (UBQ), eukaryotic elongation factor 1-αgene (eEF-1α), α-tubulin gene (α-TUB), β-tubulin gene (β-TUB), eukaryotic translation initiation factor 4α-1 gene (eIF-4α) and NAC domain protein gene(NAC) in a series of 30 M. lutarioriparia samples followed by statistical algorithms geNorm and Normfinder to analyze the gene expression stability. The results indicated that eIF-4αand UBQ were the most stable expressed genes while CBP20 showed as the least stable among all the samples. Based on above research, we recommend that at least two top-ranked reference genes should be employed for expression data normalization. The best genes selected in this study will provide a starting point to select reference genes in the future in other tissues and under other experimental conditions in this energy crop candidate.

Reference genes for gene expression studies targeting sugarcane infected with Sugarcane mosaic virus (SCMV)

Objective

The selection of reference genes in sugarcane under Sugarcane mosaic virus (SCMV) infection has not been reported and is indispensable to get reliable reverse transcription quantitative PCR (RT-qPCR) results for validation of transcriptome analysis. In this regard, seven potential reference genes were tested by RT-qPCR and ranked according to their stability using BestKeeper, NormFinder and GeNorm algorithms, and RefFinder WEB-based software in an experiment performed with samples from two sugarcane cultivars contrasting for SCMV resistance, when mechanically inoculated with a severe SCMV strain and using mock inoculated plant controls.

Results

The genes Uridylate kinase (UK) and Ubiquitin-conjugating enzyme 18 (UBC18) were the most stable according to GeNorm algorithm and the Pearson correlation coefficients with the BestKeeper index. On the other hand, ribosomal protein L35-4 (RPL1), Actin (ACT) and Ubiquitin1 (UBQ1) were the least stable genes for all algorithms tested.

Electronic supplementary material

The online version of this article (10.1186/s13104-019-4168-5) contains supplementary material, which is available to authorized users.

Screening of reference genes in real-time PCR for Radopholus similis

Six candidate reference genes were chosen from the transcriptome database of Radopholus similis using the bioinformatics method, including four conventional reference genes (actin, Eukaryotic translation initiation factor 5A (eIF5A), Tubulin alpha (a-tubulin), ubiquitin (UBI)) and two new candidate reference genes (Ribosomal protein S21 (Rps21) and Serine/threonine protein phosphatase PP1-β catalytic subunit (β-PP1)). In addition, a traditional reference gene 18S ribosomal RNA (18S rRNA) obtained from NCBI databases was also added to the analysis. Real-time PCR was used to detect the expression of seven candidate reference genes in six populations of R. similis and four developmental stages (female, male, larva and egg) of a population. The stability of the expression of candidate genes was evaluated by three software programs, BestKeeper, geNorm and NormFinder. The results showed that eIF5A is the most suitable reference gene for gene functional research of different populations, while both Rps21 and eIF5A are the most suitable reference genes for different developmental stages of a population. Therefore, eIF5A is the best reference gene for studying R. similis. However, one defect of this study is that only seven candidate reference genes were analyzed; ideally, more genes should be tested.

Screening internal controls for expression analyses involving numerous treatments by combining statistical methods with reference gene selection tools

Real-time PCR is always the method of choice for expression analyses involving comparison of a large number of treatments. It is also the favored method for final confirmation of transcript levels followed by high throughput methods such as RNA sequencing and microarray. Our analysis comprised 16 different permutation and combinations of treatments involving four different Agrobacterium strains and three time intervals in the model plant Arabidopsis thaliana. The routinely used reference genes for biotic stress analyses in plants showed variations in expression across some of our treatments. In this report, we describe how we narrowed down to the best reference gene out of 17 candidate genes. Though we initiated our reference gene selection process using common tools such as geNorm, Normfinder and BestKeeper, we faced situations where these software-selected candidate genes did not completely satisfy all the criteria of a stable reference gene. With our novel approach of combining simple statistical methods such as t test, ANOVA and post hoc analyses, along with the routine software-based analyses, we could perform precise evaluation and we identified two genes, UBQ10 and PPR as the best reference genes for normalizing mRNA levels in the context of 16 different conditions of Agrobacterium infection. Our study emphasizes the usefulness of applying statistical analyses along with the reference gene selection software for reference gene identification in experiments involving the comparison of a large number of treatments.

Effects of 5-Aminolevulinic Acid on Gene Expression, Immunity, and ATP Levels in Pacific White Shrimp, Litopenaeus vannamei

With the emergence of several infectious diseases in shrimp aquaculture, there is a growing interest in the use of feed additives to enhance shrimp immunity. Recently, the use of 5-aminolevulinic acid (5-ALA), a non-protein amino acid that plays a rate-limiting role in heme biosynthesis, has received attention for its positive effect on immunity in livestock animals. To evaluate the effect of 5-ALA in the Pacific white shrimp, Litopenaeus vannamei, we conducted microarray analysis, a Vibrio parahaemolyticus immersion challenge test, an ATP level assay, and gene expression analysis of some hemoproteins and genes associated with heme synthesis and degradation. Out of 15,745 L. vannamei putative genes on the microarray, 101 genes were differentially expressed by more than fourfold (p < 0.05) between 5-ALA-supplemented and control shrimp hepatopancreas. 5-ALA upregulated 99 of the 101 genes, 41 of which were immune- and defense-related genes based on sequence homology. Compared to the control, the 5-ALA-supplemented group had a higher survival rate in the challenge test, higher transcript levels of porphobilinogen synthase, ferrochelatase, catalase, nuclear receptor E75, and heme oxygenase-1 and higher levels of ATP. These findings suggest that dietary 5-ALA enhanced the immune response of L. vannamei to V. parahaemolyticus, upregulated immune- and defense-related genes, and enhanced aerobic energy metabolism, respectively. Further studies are needed to elucidate the extent of 5-ALA use in shrimp culture.

Making a Better Home: Modulation of Plant Defensive Response by Brevipalpus Mites

False-spider mites of the genus Brevipalpus are highly polyphagous pests that attack hundreds of plant species of distinct families worldwide. Besides causing direct damage, these mites may also act as vectors of many plant viruses that threaten high-value ornamental plants like orchids and economically important crops such as citrus and coffee. To better understand the molecular mechanisms behind plant-mite interaction we used an RNA-Seq approach to assess the global response of Arabidopsis thaliana (Arabidopsis) plants along the course of the infestation with Brevipalpus yothersi, the main vector species within the genus. Mite infestation triggered a drastic transcriptome reprogramming soon at the beginning of the interaction and throughout the time course, deregulating 1755, 3069 and 2680 genes at 6 hours after infestation (hai), 2 days after infestation (dai), and 6 dai, respectively. Gene set enrichment analysis revealed a clear modulation of processes related to the plant immune system. Co-expressed genes correlated with specific classes of transcription factors regulating defense pathways and developmental processes. Up-regulation of defensive responses correlated with the down-regulation of growth-related processes, suggesting the triggering of the growth-defense crosstalk to optimize plant fitness. Biological processes (BPs) enriched at all time points were markedly related to defense against herbivores and other biotic stresses involving the defense hormones salicylic acid (SA) and jasmonic acid (JA). Levels of both hormones were higher in plants challenged with mites than in the non-infested ones, supporting the simultaneous induction of genes from both pathways. To further clarify the functional relevance of the plant hormonal pathways on the interaction, we evaluated the mite performance on Arabidopsis mutants impaired in SA- or JA-mediated response. Mite oviposition was lower on mutants defective in SA biosynthesis (sid2) and signaling (npr1), showing a function for SA pathway in improving the mite reproduction, an unusual mechanism compared to closely-related spider mites. Here we provide the first report on the global and dynamic plant transcriptome triggered by Brevipalpus feeding, extending our knowledge on plant-mite interaction. Furthermore, our results suggest that Brevipalpus mites manipulate the plant defensive response to render the plant more susceptible to their colonization by inducing the SA-mediated pathway.

QTLs and eQTLs mapping related to citrandarins' resistance to citrus gummosis disease

Background

Phytophthora nicotianae Breda de Haan (Phytophthora parasitica Dastur) causes severe damage to citrus crops worldwide. A population of citrandarins was created from the cross between the susceptible parent Citrus sunki Hort. Ex Tan. and the resistant parent Poncirus trifoliata (L.) Raf. cv. Rubidoux, both parents and two reference rootstocks (Rangpur lime and Swingle citrumelo) were grafted in a greenhouse on Rangpur lime. Inoculations were performed at 10 cm and 15 cm above the grafting region and the resulting lesions were evaluated by measuring the lesion length 60 days after inoculation. As control, non-inoculated plants of each genotype were used. In addition, we evaluated the expression of 19 candidate genes involved in citrus defense response 48 h after pathogen infection by quantitative Real-Time PCR (qPCR). We mapped genomic regions of Quantitative Trait Loci (QTLs) and Expression Quantitative Trait Loci (eQTLs) associated with resistance to P. parasitica in the linkage groups (LGs) of the previously constructed maps of C. sunki and P. trifoliata.

Results

We found disease severity differences among the generated hybrids, with lesion lengths varying from 1.15 to 11.13 mm. The heritability of the character was 65%. These results indicate that there is a great possibility of success in the selection of resistant hybrids within this experiment. The analysis of gene expression profile demonstrated a great variation of responses regarding the activation of plant defense pathways, indicating that citrandarins have several defense strategies to control oomycete infection. The information of the phenotypic and gene expression data made possible to detect genomic regions associated with resistance. Three QTLs and 84 eQTLs were detected in the linkage map of P. trifoliata, while one QTL and 110 eQTLs were detected in C. sunki.

Conclusions

This is the first study to use eQTLs mapping in the Phytophthora-citrus interaction. Our results from the QTLs and eQTLs mapping allow us to conclude that the resistance of some citrandarins to the infection by P. parasitica is due to a favorable combination of QTLs and eQTLs transmitted by both parents.

Electronic supplementary material

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

Identification and evaluation of PCR reference genes for host and pathogen in sugarcane-Sporisorium scitamineum interaction system

BACKGROUND

Sugarcane (Saccharum L. plant) is an important crop for sugar and bio-energy production around the world. Among sugarcane diseases, smut caused by Sporisorium scitamineum is one of the major fungal diseases causing severe losses to the sugarcane industry. The use of PCR reference genes is essential to the normalization of data on gene expression involving the sugarcane-S. scitamineum interaction system; however, no report that addresses criteria in selecting these reference genes has been published to date.

RESULTS

In this study, 10 sugarcane genes and eight S. scitamineum genes were selected as candidate PCR reference genes in the sugarcane-S. scitamineum interaction system. The stability and reliability of these 18 candidate genes were analyzed in smut-resistant (NCo376) and -susceptible (YC71-374) genotypes using the statistical algorithms geNorm, NormFinder, BestKeeper, and deltaCt method. Subsequently, the relative expression levels of the sugarcane chitinase I-3 gene and S. scitamineum chorismate mutase gene were determined to validate the applicability of these sugarcane and S. scitamineum PCR reference genes, respectively. We finally found that the acyl-CoA dehydrogenase gene (ACAD), serine/arginine repetitive matrix protein 1 gene (SARMp1), or their combination (ACAD + SARMp1) could be utilized as the most suitable reference genes for normalization of sugarcane gene expression in sugarcane bud tissues after S. scitamineum infection. Similarly, the inosine 5'-monophosphate dehydrogenase gene (S10), the SEC65-signal recognition particle subunit gene (S11), or their combination (S10 + S11) were suitable for normalization of S. scitamineum gene expression in sugarcane bud tissues.

CONCLUSIONS

The PCR reference genes ACAD, SARMp1, S10, and S11 may be employed in gene transcriptional studies involving the sugarcane-S. scitamineum interaction system.

Selection and validation of reference genes for quantitative real-time PCR in Artemisia sphaerocephala based on transcriptome sequence data

Artemisia sphaerocephala, a dicotyledonous perennial semi-shrub belonging to the Artemisia genus of the Compositae family, is widely distributed in northwestern China. This shrub is one of the most important pioneer plants which is capable of protecting rangelands from wind erosion. It therefore plays a vital role in maintaining desert ecosystem stability. In addition, to its use as a forage grass, it has excellent prospective applications as a source of plant oil and as a plant-based fuel. The use of internal genes is the basis for accurately assessing Real time quantitative PCR. In this study, based on transcriptome data of A. sphaerocephala, we analyzed 21 candidate internal genes to determine the optimal internal genes in this shrub. The stabilities of candidate genes were evaluated in 16 samples of A. sphaerocephala. Finally, UBC9 and TIP41-like were determined as the optimal reference genes in A. sphaerocephala by Delta Ct and three various programs. There were GeNorm, NormFinder and BestKeeper.

Reference gene identification for reliable normalisation of quantitative RT-PCR data in Setaria viridis

BACKGROUND

Quantitative real-time polymerase chain reaction (RT-qPCR) is the key platform for the quantitative analysis of gene expression in a wide range of experimental systems and conditions. However, the accuracy and reproducibility of gene expression quantification via RT-qPCR is entirely dependent on the identification of reliable reference genes for data normalisation. Green foxtail (Setaria viridis) has recently been proposed as a potential experimental model for the study of C₄ photosynthesis and is closely related to many economically important crop species of the Panicoideae subfamily of grasses, including Zea mays (maize), Sorghum bicolor (sorghum) and Sacchurum officinarum (sugarcane). Setaria viridis (Accession 10) possesses a number of key traits as an experimental model, namely; (i) a small sized, sequenced and well annotated genome; (ii) short stature and generation time; (iii) prolific seed production, and; (iv) is amendable to Agrobacterium tumefaciens-mediated transformation. There is currently however, a lack of reference gene expression information for Setaria viridis (S. viridis). We therefore aimed to identify a cohort of suitable S. viridis reference genes for accurate and reliable normalisation of S. viridis RT-qPCR expression data.

RESULTS

Eleven putative candidate reference genes were identified and examined across thirteen different S. viridis tissues. Of these, the geNorm and NormFinder analysis software identified SERINE/THERONINE-PROTEIN PHOSPHATASE 2A (PP2A), 5'-ADENYLYLSULFATE REDUCTASE 6 (ASPR6) and DUAL SPECIFICITY PHOSPHATASE (DUSP) as the most suitable combination of reference genes for the accurate and reliable normalisation of S. viridis RT-qPCR expression data. To demonstrate the suitability of the three selected reference genes, PP2A, ASPR6 and DUSP, were used to normalise the expression of CINNAMYL ALCOHOL DEHYDROGENASE (CAD) genes across the same tissues.

CONCLUSIONS

This approach readily demonstrated the suitably of the three selected reference genes for the accurate and reliable normalisation of S. viridis RT-qPCR expression data. Further, the work reported here forms a highly useful platform for future gene expression quantification in S. viridis and can also be potentially directly translatable to other closely related and agronomically important C₄ crop species.

A mobile loop near the active site acts as a switch between the dual activities of a viral protease/deubiquitinase

The positive-strand RNA virus Turnip yellow mosaic virus (TYMV) encodes an ovarian tumor (OTU)-like protease/deubiquitinase (PRO/DUB) protein domain involved both in proteolytic processing of the viral polyprotein through its PRO activity, and in removal of ubiquitin chains from ubiquitylated substrates through its DUB activity. Here, the crystal structures of TYMV PRO/DUB mutants and molecular dynamics simulations reveal that an idiosyncratic mobile loop participates in reversibly constricting its unusual catalytic site by adopting "open", "intermediate" or "closed" conformations. The two cis-prolines of the loop form a rigid flap that in the most closed conformation zips up against the other side of the catalytic cleft. The intermediate and closed conformations also correlate with a reordering of the TYMV PRO/DUB catalytic dyad, that then assumes a classical, yet still unusually mobile, OTU DUB alignment. Further structure-based mutants designed to interfere with the loop's mobility were assessed for enzymatic activity in vitro and in vivo, and were shown to display reduced DUB activity while retaining PRO activity. This indicates that control of the switching between the dual PRO/DUB activities resides prominently within this loop next to the active site. Introduction of mutations into the viral genome revealed that the DUB activity contributes to the extent of viral RNA accumulation both in single cells and in whole plants. In addition, the conformation of the mobile flap was also found to influence symptoms severity in planta. Such mutants now provide powerful tools with which to study the specific roles of reversible ubiquitylation in viral infection.

Viruses have much smaller genomes than their hosts. Consequently, they often encode proteins which are multifunctional. For instance, some viral proteases have a dual function, being also deubiquitinases, i.e. enzymes capable of removing ubiquitin tags grafted onto proteins and that often target them for destruction. The protease and deubiquitinase activities share a single active site that is used alternately for one function or the other, but how this switch between activities may be regulated is presently unknown. To answer this question, we studied a simple plant virus that is a useful model system for these complex molecular biology phenomena, and that encodes a simplified protease/deubiquitinase. Here, thanks to a combination of structural and functional analyses, we managed to decouple the two activities, killing the deubiquitinase activity while preserving the protease one. This successful decoupling relies on our discovery that a loop inserted next to the active site is mobile, and can thus act as a switch between the two activities. This result allowed us to demonstrate the importance of the specific deubiquinase activity in viral multiplication. In addition, viral symptoms were also severely affected by mutations affecting the loop mobility. Our data provide powerful tools for further studies, that may also be relevant for more complex or medically relevant viruses.

A spiking strategy facilitates housekeeping selection for RT-qPCR analysis under different biotic stresses in eggplant

Endogenous housekeeping genes are traditionally employed to normalize the expression of target genes in RT-qPCR studies. Assuming that a perfect housekeeping suitable for every condition does not exist, expression stability of the chosen reference gene should be evaluated at every new experiment. The housekeeping selection process reveals furthermore complicated and time-consuming when different conditions have to be compared in the same experimental dataset. As an alternative strategy, we spiked an external reference transcript (ERT) into all RNA samples of our dataset (eggplant roots subjected to different biotic stresses), and used it to normalize the expression levels of native candidate housekeeping. ERT expression resulted highly stable across all samples and enabled to indicate glyceraldehyde 3-phosphate dehydrogenase (GAPDH) as the most stable endogenous housekeeping. This result was confirmed by the use of GeNorm, Normfinder, and BestKeeper algorithms. This method might be generally applied to expedite the selection process of the best reference gene.

Validation of reference genes for quantifying changes in gene expression in virus-infected tobacco

To facilitate quantification of gene expression changes in virus-infected tobacco plants, eight housekeeping genes were evaluated for their stability of expression during infection by one of three systemically-infecting viruses (cucumber mosaic virus, potato virus X, potato virus Y) or a hypersensitive-response-inducing virus (tobacco mosaic virus; TMV) limited to the inoculated leaf. Five reference-gene validation programs were used to establish the order of the most stable genes for the systemically-infecting viruses as ribosomal protein L25 > β-Tubulin > Actin, and the least stable genes Ubiquitin-conjugating enzyme (UCE) < PP2A < GAPDH. For local infection by TMV, the most stable genes were EF1α > Cysteine protease > Actin, and the least stable genes were GAPDH < PP2A < UCE. Using two of the most stable and the two least stable validated reference genes, three defense responsive genes were examined to compare their relative changes in gene expression caused by each virus.

Identification and comprehensive evaluation of reference genes for RT-qPCR analysis of host gene-expression in Brassica juncea-aphid interaction using microarray data

Brassica juncea is a chief oil yielding crop in many parts of the world including India. With advancement of molecular techniques, RT-qPCR based study of gene-expression has become an integral part of experimentations in crop breeding. In RT-qPCR, use of appropriate reference gene(s) is pivotal. The virtue of the reference genes, being constant in expression throughout the experimental treatments, needs to be validated case by case. Appropriate reference gene(s) for normalization of gene-expression data in B. juncea during the biotic stress of aphid infestation is not known. In the present investigation, 11 reference genes identified from microarray database of Arabidopsis-aphid interaction at a cut off FDR ≤0.1, along with two known reference genes of B. juncea, were analyzed for their expression stability upon aphid infestation. These included 6 frequently used and 5 newly identified reference genes. Ranking orders of the reference genes in terms of expression stability were calculated using advanced statistical approaches such as geNorm, NormFinder, delta Ct and BestKeeper. The analysis suggested CAC, TUA and DUF179 as the most suitable reference genes. Further, normalization of the gene-expression data of STP4 and PR1 by the most and the least stable reference gene, respectively has demonstrated importance and applicability of the recommended reference genes in aphid infested samples of B. juncea.

Lipid flippases promote antiviral silencing and the biogenesis of viral and host siRNAs in Arabidopsis

Dicer-mediated processing of virus-specific dsRNA into short interfering RNAs (siRNAs) in plants and animals initiates a specific antiviral defense by RNA interference (RNAi). In this study, we developed a forward genetic screen for the identification of host factors required for antiviral RNAi in Arabidopsis thaliana Using whole-genome sequencing and a computational pipeline, we identified aminophospholipid transporting ATPase 2 (ALA2) and the related ALA1 in the type IV subfamily of P-type ATPases as key components of antiviral RNAi. ALA1 and ALA2 are flippases, which are transmembrane lipid transporter proteins that transport phospholipids across cellular membranes. We found that the ala1/ala2 single- and double-mutant plants exhibited enhanced disease susceptibility to cucumber mosaic virus when the virus-encoded function to suppress RNAi was disrupted. Notably, the antiviral activity of both ALA1 and ALA2 was abolished by a single amino acid substitution known to inactivate the flippase activity. Genetic analysis revealed that ALA1 and ALA2 acted to enhance the amplification of the viral siRNAs by RNA-dependent RNA polymerase (RdRP) 1 (RDR1) and RDR6 and of the endogenous virus-activated siRNAs by RDR1. RNA virus replication by plant viral RdRPs occurs inside vesicle-like membrane invaginations induced by the recruitment of the viral RdRP and host factors to subcellular membrane microdomains enriched with specific phospholipids. Our results suggest that the phospholipid transporter activity of ALA1/ALA2 may be necessary for the formation of similar invaginations for the synthesis of dsRNA precursors of highly abundant viral and host siRNAs by the cellular RdRPs.

Identification and Molecular Characterization of the Chloroplast Targeting Domain of Turnip yellow mosaic virus Replication Proteins

Turnip yellow mosaic virus (TYMV) is a positive-strand RNA virus infecting plants. The TYMV 140K replication protein is a key organizer of viral replication complex (VRC) assembly, being responsible for recruitment of the viral polymerase and for targeting the VRCs to the chloroplast envelope where viral replication takes place. However, the structural requirements determining the subcellular localization and membrane association of this essential viral protein have not yet been defined. In this study, we investigated determinants for the in vivo chloroplast targeting of the TYMV 140K replication protein. Subcellular localization studies of deletion mutants identified a 41-residue internal sequence as the chloroplast targeting domain (CTD) of TYMV 140K; this sequence is sufficient to target GFP to the chloroplast envelope. The CTD appears to be located in the C-terminal extension of the methyltransferase domain-a region shared by 140K and its mature cleavage product 98K, which behaves as an integral membrane protein during infection. We predicted the CTD to fold into two amphipathic α-helices-a folding that was confirmed in vitro by circular dichroism spectroscopy analyses of a synthetic peptide. The importance for subcellular localization of the integrity of these amphipathic helices, and the function of 140K/98K, was demonstrated by performing amino acid substitutions that affected chloroplast targeting, membrane association and viral replication. These results establish a short internal α-helical peptide as an unusual signal for targeting proteins to the chloroplast envelope membrane, and provide new insights into membrane targeting of viral replication proteins-a universal feature of positive-strand RNA viruses.