Evaluation and selection of internal reference genes from two- and six-row U.S. malting barley varieties throughout micromalting for use in RT-qPCR

Selection of Reference Genes for Expression Normalization by RT-qPCR in Dracocephalum moldavica L

Dracocephalum moldavica is widely used as an ornamental, medicine, and perfume in industry. Real-time fluorescence quantitative polymerase chain reaction (RT-qPCR) is widely and accurately utilized for gene expression evaluations. Selecting optimal reference genes is essential for normalizing RT-qPCR results. However, the identification of suitable reference genes in D. moldavica has not been documented. A total of 12 reference genes in D. moldavica were identified by PEG6000 (15%) treatment under hypertonia conditions in different tissues (roots, stem, leaves, flower, seeds and sepal) and during three stages of flower development, then used to validate the expression stability. There were four algorithms (delta Ct, geNorm, NormFinder, and BestKeeper) used to analyze the stability. Finally, the RefFinder program was employed to evaluate the candidate reference genes' stability. The results showed that ACTIN, glyceraldehyde-3-phosphate dehydrogenase (GAPDH), and EF1α (elongation factor-1α) were stable reference genes under the PEG6000 treatment. Heat shock protein 70 (HSP70) was the most stable gene across different flower development stages. ADP-ribosylation factor (ARF) was the most stable gene in different tissues and total samples. This study provides reliable gene expression studies for future research in D. moldavica.

Exogenous spike-in mouse RNAs for accurate differential gene expression analysis in barley using RT-qPCR

Reverse transcriptase quantitative polymerase chain reaction (RT-qPCR) followed by the 2-ΔΔCt method is the most common way to measure transcript levels for relative gene expression assays. The quality of an RT-qPCR assay is dependent upon the identification and validation of reference genes to normalize gene expression data. The so-called housekeeping genes are commonly used as internal reference genes because they are assumed to be ubiquitously expressed at stable levels. Commonly, researchers do not validate their reference genes but rely on historical reference genes or previously validated genes from an unrelated experiment. Using previously validated reference genes to assess gene expression changes occurring during malting resulted in extensive variability. Therefore, a new method was tested and validated to circumvent the use of internal reference genes. Total mouse RNA was chosen as the external reference RNA and a suite of primer sets to putatively stable mouse genes was created to identify stably expressed genes for use as an external reference gene. cDNA was created by co-amplifying total mouse RNA, as an RNA spike-in, and barley RNA. When using the external reference genes to normalize malting gene expression data, standard deviations were significantly reduced and significant differences in transcript abundance were observed, whereas when using the internal reference genes, standard deviations were larger with no significant differences seen. Furthermore, external reference genes were more accurate at assessing expression levels in malting and developing grains, whereas the internal reference genes overestimated abundance in developing grains and underestimated abundance in malting grains.

Selection of suitable internal control gene for assaying gene expression in rice through qRT-PCR during sheath blight infection

qRT-PCR is a globally accepted technique for assaying gene expression in relative terms which compares the difference between critical threshold (C_T) values of a gene calculated form two independently isolated RNA samples. Independent RNA isolations, however, include error due to batch effect which must be normalized for error-free calculation of relative gene expression. Hence, C_T values of internal control (IC) genes are used for normalization during the calculation of expression fold-change in gene expression analysis. The expression of ICs genes expected to be stable in all the experimental conditions. However, it is almost impossible to find such a gene which do not depict expression fluctuation in response to the changes in experimental conditions. Hence, it is necessary to identify suitable IC gene(s) for any given experimental condition before conducting any particular gene expression study. Here, we examined the suitability of eight candidate IC genes, namely glyceraldehyde 3-phosphate dehydrogenase (GAPDH), eukaryotic elongation factor-1 (eEF-1α), 25 S rRNA (25 S), 18 S rRNA (18 S), ubiquitin C E2 ligase (UBC), Actin (Act), ubiquitin 5 (UBQ5) and ubiquitin 10 (UBQ10), for assaying gene expression in rice during sheath blight infection. Our analysis suggest that GAPDH might be the IC of choice when expression studies include contrasting genotypes differing in their tolerance to sheath blight pathogen as well as progressive infection time. While if expression analysis have to be performed only in one genotype but under progressive sheath blight infection, UBQ5 might be chosen as IC because of its high expression stability under the proposed experimental setup.

The TCP transcription factor HvTB2 heterodimerizes with VRS5 and controls spike architecture in barley

Understanding the molecular network, including protein-protein interactions, of VRS5 provide new routes towards the identification of other key regulators of plant architecture in barley. The TCP transcriptional regulator TEOSINTE BRANCHED 1 (TB1) is a key regulator of plant architecture. In barley, an important cereal crop, HvTB1 (also referred to as VULGARE SIX-ROWED spike (VRS) 5), inhibits the outgrowth of side shoots, or tillers, and grains. Despite its key role in barley development, there is limited knowledge on the molecular network that is utilized by VRS5. In this work, we performed protein-protein interaction studies of VRS5. Our analysis shows that VRS5 potentially interacts with a diverse set of proteins, including other class II TCP's, NF-Y TF, but also chromatin remodelers. Zooming in on the interaction capacity of VRS5 with other TCP TFs shows that VRS5 preferably interacts with other class II TCP TFs in the TB1 clade. Induced mutagenesis through CRISPR-Cas of one of the putative VRS5 interactors, HvTB2 (also referred to as COMPOSITUM 1 and BRANCHED AND INDETERMINATE SPIKELET 1), resulted in plants that have lost their characteristic unbranched spike architecture. More specifically, hvtb2 mutants exhibited branches arising at the main spike, suggesting that HvTB2 acts as inhibitor of branching. Our protein-protein interaction studies of VRS5 resulted in the identification of HvTB2 as putative interactor of VRS5, another key regulator of spike architecture in barley. The study presented here provides a first step to underpin the protein-protein interactome of VRS5 and to identify other, yet unknown, key regulators of barley plant architecture.

Abscisic acid activation of oleosin gene HvOle3 expression prevents the coalescence of protein storage vacuoles in barley aleurone cells

Protein storage vacuoles (PSVs) in aleurone cells coalesce during germination, and this process is highly coupled with mobilization of PSV reserves, allowing de novo synthesis of various hydrolases in aleurone cells for endosperm degradation. Here we show that in barley (Hordeum vulgare L.) oleosins, the major integral proteins of oleosomes are encoded by four genes (HvOle1 to 4), and the expression of HvOle1 and HvOle3 is strongly up-regulated by abscisic acid (ABA), which shows antagonism to gibberellic acid. In aleurone cells, all HvOLEs were subcellularly targeted to the tonoplast of PSVs. Gain-of-function analyses revealed that HvOLE3 effectively delayed PSV coalescence, whereas HvOLE1 only had a moderate effect, with no notable effect of HvOLE2 and 4. With regard to longevity, HvOLE3 chiefly outperformed other HvOLEs, followed by HvOLE1. Experiments swapping the N- and C-terminal domain between HvOLE3 and other HvOLEs showed that the N-terminal region of HvOLE3 is mainly responsible, with some positive effect by the C-terminal region, for mediating the specific preventive effect of HvOLE3 on PSV coalescence. Three ACGT-core elements and the RY-motif were responsible for ABA induction of HvOle3 promoter activity. Transient expression assays using aleurone protoplasts demonstrated that transcriptional activation of the HvOle3 promoter was mediated by transcription factors HvABI3 and HvABI5, which acted downstream of protein kinase HvPKABA1.

Description and functional analysis of the transcriptome from malting barley

The present study aimed to establish an early model of the malting barley transcriptome, which describes the expression of genes and their ontologies, identify the period during malting with the largest dynamic shift in gene expression for future investigation, and to determine the expression patterns of all starch degrading enzyme genes relevant to the malting and brewing industry. Large dynamic increases in gene expression occurred early in malting with differential expressed genes enriched for cell wall and starch hydrolases amongst many malting related categories. Twenty-five of forty starch degrading enzyme genes were differentially expressed in the malting barley transcriptome including eleven α-amylase genes, six β-amylase genes, three α-glucosidase genes, and all five starch debranching enzyme genes. Four new or novel α-amylase genes, one β-amylase gene (Bmy3), three α-glucosidase genes, and two isoamylase genes had appreciable expression that requires further exploration into their potential relevance to the malting and brewing industry.

The Barley (Hordeum vulgare ssp. vulgare) Respiratory Burst Oxidase Homolog (HvRBOH) Gene Family and Their Plausible Role on Malting Quality

Controlled generation of reactive oxygen species (ROS) is pivotal for normal plant development and adaptation to changes in the external milieu. One of the major enzymatic sources of ROS in plants are the plasma-membrane localized NADPH oxidases, also called as Respiratory Burst Oxidase Homologs (RBOH). In addition to the six previously reported, seven new members of RBOH gene family were identified in barley using in silico analysis. Conservation of genomic structure and key residues important for catalytic activity and co-factor binding was observed in barley RBOH genes. Phylogenetic analysis of plant RBOHs revealed distinct clades for monocot and dicot RBOH proteins. Hence, we propose to use the rice nomenclature for naming barley RBOH genes. Temporal changes in ROS profiles were observed during barley malting and was accompanied by changes in protein carbonylation, lipid peroxidation, and antioxidant capacity. Among the nine differentially expressed HvRBOHs during various malting stages, HvRBOHA and HvRBOHC showed most significant sustained changes in expression. RNAi knockdown lines with reduced expression of HvRBOHA/C gene exhibited genetic compensation via inducible expression of other gene family members during malting. However, the physiological consequence of reduced expression of HvRBOHA/C manifested as a poor malting quality profile attributable to low alpha-amylase activity and high levels of beta-glucan. We propose that the HvRBOHs play a critical role in modulating the redox milieu during the early stages of malting, which in turn can significantly impact carbohydrate metabolism.

Identification of qPCR reference genes suitable for normalizing gene expression in the mdx mouse model of Duchenne muscular dystrophy

The mdx mouse is the most widely-used animal model of the human disease Duchenne muscular dystrophy, and quantitative PCR analysis of gene expression in the muscles of this animal plays a key role in the study of pathogenesis and disease progression and in evaluation of potential therapeutic interventions. Normalization to appropriate stably-expressed reference genes is essential for accurate quantitative measurement, but determination of such genes is challenging: healthy and dystrophic muscles present very different transcriptional environments, further altering with disease progression and muscle use, raising the possibility that no single gene or combination of genes may be stable under all experimental comparative scenarios. Despite the pedigree of this animal model, this problem remains unaddressed. The aim of this work was therefore to comprehensively assess reference gene suitability in the muscles of healthy and dystrophic mice, identifying reference genes appropriate for specific experimental comparisons, and determining whether an essentially universally-applicable set of genes exists. Using a large sample collection comprising multiple muscles (including the tibialis anterior, diaphragm and heart muscles) taken from healthy and mdx mice at three disease-relevant ages, and a panel of sixteen candidate reference genes (FBXO38, FBXW2, MON2, ZFP91, HTATSF1, GAPDH, ACTB, 18S, CDC40, SDHA, RPL13a, CSNK2A2, AP3D1, PAK1IP1, B2M and HPRT1), we used the geNorm, BestKeeper and Normfinder algorithms to identify genes that were stable under multiple possible comparative scenarios. We reveal that no single gene is stable under all conditions, but a normalization factor derived from multiple genes (RPL13a, CSNK2A2, AP3D1 and the widely-used ACTB) appears suitable for normalizing gene expression in both healthy and dystrophic mouse muscle regardless of muscle type or animal age. We further show that other popular reference genes, including GAPDH, are markedly disease- or muscle-type correlated. This study demonstrates the importance of empirical reference gene identification, and should serve as a valuable resource for investigators wishing to study gene expression in mdx mice.