Application of qRT-PCR and RNA-Seq analysis for the identification of housekeeping genes useful for normalization of gene expression values during Striga hermonthica development

Effect of Trichomonacide 6-Nitro-1H-benzimidazole Derivative Compounds on Expression Level of Metabolic Genes in Trichomonas vaginalis

The parasite Trichomonas vaginalis is the etiologic agent of trichomoniasis, the most common non-viral sexually transmitted disease worldwide. This infection often remains asymptomatic and is related to several health complications. The traditional treatment for trichomoniasis is the use of drugs of the 5-nitroimidazole family, such as metronidazole; however, scientific reports indicate an increasing number of drug-resistant strains. Benzimidazole derivatives could offer an alternative in the search for new anti-trichomonas drugs. In this sense, two attractive candidates are the compounds O2N-BZM7 and O2N-BZM9 (1H-benzimidazole derivatives), since, through in vitro tests, they have shown a higher trichomonacide activity. In this study, we determined the effect on the expression level of metabolic genes in T. vaginalis. The results show that genes involved in redox balance (NADHOX, G6PD::6PGL) are overexpressed, as well as the gene that participates in the first reaction of glycolysis (CK); on the other hand, structural genes such as ACT and TUB are decreased in expression in trophozoites treated with the compound O2N-BZM9, which would probably affect its morphology, motility and virulence. These results align with the trichomonacidal activity of the compounds, with benzimidazole O2N-BZM9 being the most potent, with an IC50 value of 4.8 μM. These results are promising for potential future therapeutic applications.

Gibberellins Promote Seed Conditioning by Up-Regulating Strigolactone Receptors in the Parasitic Plant Striga hermonthica

Dormant seeds of the root parasitic plant Striga hermonthica sense strigolactones from host plants as environmental cues for germination. This process is mediated by a diversified member of the strigolactone receptors encoded by HYPOSENSITIVE TO LIGHT/KARRIKIN INSENSITIVE2 genes. It is known that warm and moist treatment during seed conditioning gradually makes dormant Striga seeds competent to respond to strigolactones, although the mechanism behind it is poorly understood. In this report, we show that plant hormone gibberellins increase strigolactone competence by up-regulating mRNA expression of the major strigolactone receptors during the conditioning period. This idea was supported by a poor germination phenotype in which gibberellin biosynthesis was depleted by paclobutrazol during conditioning. Moreover, live imaging with a fluorogenic strigolactone mimic, yoshimulactone green W, revealed that paclobutrazol treatment during conditioning caused aberrant dynamics of strigolactone perception after germination. These observations revealed an indirect role of gibberellins in seed germination in Striga, which contrasts with their roles as dominant germination-stimulating hormones in non-parasitic plants. We propose a model of how the role of gibberellins became indirect during the evolution of parasitism in plants. Our work also highlights the potential role for gibberellins in field applications, for instance, in elevating the sensitivity of seeds toward strigolactones in the current suicidal germination approach to alleviate the agricultural threats caused by this parasite in Africa.

Genome-enabled discovery of candidate virulence loci in Striga hermonthica, a devastating parasite of African cereal crops

Parasites have evolved proteins, virulence factors (VFs), that facilitate plant colonisation, however VFs mediating parasitic plant-host interactions are poorly understood. Striga hermonthica is an obligate, root-parasitic plant of cereal hosts in sub-Saharan Africa, causing devastating yield losses. Understanding the molecular nature and allelic variation of VFs in S. hermonthica is essential for breeding resistance and delaying the evolution of parasite virulence. We assembled the S. hermonthica genome and identified secreted proteins using in silico prediction. Pooled sequencing of parasites growing on a susceptible and a strongly resistant rice host allowed us to scan for loci where selection imposed by the resistant host had elevated the frequency of alleles contributing to successful colonisation. Thirty-eight putatively secreted VFs had very different allele frequencies with functions including host cell wall modification, protease or protease inhibitor and kinase activities. These candidate loci had significantly higher Tajima's D than the genomic background, consistent with balancing selection. Our results reveal diverse strategies used by S. hermonthica to overcome different layers of host resistance. Understanding the maintenance of variation at virulence loci by balancing selection will be critical to managing the evolution of virulence as part of a sustainable control strategy.

Plant science's next top models

BACKGROUND

Model organisms are at the core of life science research. Notable examples include the mouse as a model for humans, baker's yeast for eukaryotic unicellular life and simple genetics, or the enterobacteria phage λ in virology. Plant research was an exception to this rule, with researchers relying on a variety of non-model plants until the eventual adoption of Arabidopsis thaliana as primary plant model in the 1980s. This proved to be an unprecedented success, and several secondary plant models have since been established. Currently, we are experiencing another wave of expansion in the set of plant models.

SCOPE

Since the 2000s, new model plants have been established to study numerous aspects of plant biology, such as the evolution of land plants, grasses, invasive and parasitic plant life, adaptation to environmental challenges, and the development of morphological diversity. Concurrent with the establishment of new plant models, the advent of the 'omics' era in biology has led to a resurgence of the more complex non-model plants. With this review, we introduce some of the new and fascinating plant models, outline why they are interesting subjects to study, the questions they will help to answer, and the molecular tools that have been established and are available to researchers.

CONCLUSIONS

Understanding the molecular mechanisms underlying all aspects of plant biology can only be achieved with the adoption of a comprehensive set of models, each of which allows the assessment of at least one aspect of plant life. The model plants described here represent a step forward towards our goal to explore and comprehend the diversity of plant form and function. Still, several questions remain unanswered, but the constant development of novel technologies in molecular biology and bioinformatics is already paving the way for the next generation of plant models.

Genome-wide identification of internal reference genes for normalization of gene expression values during endosperm development in wheat

Internal reference genes that are stably expressed are essential for normalization in comparative expression analyses. However, gene expression varies significantly among species, organisms, tissues, developmental stages, stresses, and treatments. Therefore, identification of stably expressed reference genes in developmental endosperm of bread wheat is important for expression analysis of endosperm genes. As the first study to systematically screen for reference genes across different developmental stages of wheat endosperm, nine genes were selected from among 76 relatively stable genes based on high-throughput RNA sequencing data. The expression stability of these candidate genes and five traditional reference genes was assessed by real-time quantitative PCR combined with three independent algorithms: geNorm, NormFinder, and BestKeeper. The results showed that ATG8d was the most stable gene during wheat endosperm development, followed by Ta54227, while the housekeeping gene GAPDH, commonly used as an internal reference, was the least stable. ATG8d and Ta54227 together formed the optimal combination of reference genes. Comparative expression analysis of glutenin genes indicated that credible quantification could be achieved by normalization against ATG8d in developmental endosperm. The stably expressed gene characterized here can act as a proper internal reference for expression analysis of wheat endosperm genes, especially nutrient- and nutrient synthesis-related genes.

Identification of Suitable Reference Genes for RT-qPCR Assays in Liriodendron chinense (Hemsl.) Sarg

The precision and reliability of reverse transcription quantitative polymerase chain reaction (RT-qPCR) depend mainly on suitable reference genes; however, reference genes have not yet been identified for Liriodendron chinense (Hemsl.) Sarg. In this study, the expression stability of 15 candidate reference genes, ACT7, ACT97, UBQ1, eIF2, eIF3, HIS, BIG, AGD11, EFG, GAPDH, CYP, RPL25, UBC, RPB1, and TUB, was tested across multiple organs of L. chinense using four algorithms, geNorm, NormFinder, BestKeeper, and RefFinder. To understand the difference between the selected reference genes and the unsuitable candidate reference genes, the expression level of a target gene, LcPAT7, was normalized across various plant samples. ACT97 and eIF3 represented the best combination across all samples tested, while AGD11 and UBQ1 were unsuitable for normalization in this case. In the vegetative organ subset, ACT97, ACT7, and GAPDH showed the highest expression stability. For floral organs, UBC and eIF3 were the most stable reference genes. Unsuitable reference genes underestimated the expression levels of a target gene, LcPAT7. This study identified two reference genes (ACT97 and eIF3) for the precise and reliable normalization of L. chinense RT-qPCR data across various organs. Our work provides an effective framework for quantifying gene expression in L. chinense.

An explanatory model of temperature influence on flowering through whole-plant accumulation of FLOWERING LOCUS T in Arabidopsis thaliana

We assessed mechanistic temperature influence on flowering by incorporating temperature-responsive flowering mechanisms across developmental age into an existing model. Temperature influences the leaf production rate as well as expression of FLOWERING LOCUS T (FT), a photoperiodic flowering regulator that is expressed in leaves. The Arabidopsis Framework Model incorporated temperature influence on leaf growth but ignored the consequences of leaf growth on and direct temperature influence of FT expression. We measured FT production in differently aged leaves and modified the model, adding mechanistic temperature influence on FT transcription, and causing whole-plant FT to accumulate with leaf growth. Our simulations suggest that in long days, the developmental stage (leaf number) at which the reproductive transition occurs is influenced by day length and temperature through FT, while temperature influences the rate of leaf production and the time (in days) the transition occurs. Further, we demonstrate that FT is mainly produced in the first 10 leaves in the Columbia (Col-0) accession, and that FT accumulation alone cannot explain flowering in conditions in which flowering is delayed. Our simulations supported our hypotheses that: (i) temperature regulation of FT, accumulated with leaf growth, is a component of thermal time, and (ii) incorporating mechanistic temperature regulation of FT can improve model predictions when temperatures change over time.

Phloem-Conducting Cells in Haustoria of the Root-Parasitic Plant Phelipanche aegyptiaca Retain Nuclei and Are Not Mature Sieve Elements

Phelipanche aegyptiaca parasitizes a wide range of plants, including important crops, and causes serious damage to their production. P. aegyptiaca develops a specialized intrusive organ called a haustorium that establishes connections to the host’s xylem and phloem. In parallel with the development of xylem vessels, the differentiation of phloem-conducting cells has been demonstrated by the translocation of symplasmic tracers from the host to the parasite. However, it is unclear yet whether haustorial phloem-conducting cells are sieve elements. In this study, we identified phloem-conducting cells in haustoria by the host-to-parasite translocation of green fluorescent protein (GFP) from AtSUC2pro::GFP tomato sieve tubes. Haustorial GFP-conducting cells contained nuclei but not callose-rich sieve plates, indicating that phloem-conducting cells in haustoria differ from conventional sieve elements. To ascertain why the nuclei were not degenerated, expression of the P. aegyptiaca homologs NAC-domain containing transcription factor (NAC45), NAC45/86-dependent exonuclease-domain protein 1 (NEN1), and NEN4 was examined. However, these genes were more highly expressed in the haustorium than in tubercle protrusion, implying that nuclear degradation in haustoria may not be exclusively controlled by the NAC45/86-NEN regulatory pathway. Our results also suggest that the formation of plasmodesmata with large size exclusion limits is independent of nuclear degradation and callose deposition.

Selection and validation of appropriate reference genes for quantitative real-time PCR analysis in Salvia hispanica

Quantitative real-time polymerase chain reaction (qRT-PCR) has become the most popular choice for gene expression studies. For accurate expression analysis, it is pertinent to select a stable reference gene to normalize the data. It is now known that the expression of internal reference genes varies considerably during developmental stages and under different experimental conditions. For Salvia hispanica, an economically important oilseed crop, there are no reports of stable reference genes till date. In this study, we chose 13 candidate reference genes viz. Actin11 (ACT), Elongation factor 1-alpha (EF1-α), Eukaryotic translation initiation factor 3E (ETIF3E), alpha tubulin (α-TUB), beta tubulin (β-TUB), Glyceraldehyde 3-phosphate dehydrogenase (GAPDH), Cyclophilin (CYP), Clathrin adaptor complex (CAC), Serine/threonine-protein phosphatase 2A (PP2A), FtsH protease (FtsH), 18S ribosomal RNA (18S rRNA), S-adenosyl methionine decarboxylase (SAMDC) and Rubisco activase (RCA) and the expression levels of these genes were assessed in a diverse set of tissue samples representing vegetative stages, reproductive stages and various abiotic stress treatments. Two of the widely used softwares, geNorm and Normfinder were used to evaluate the expression stabilities of these 13 candidate reference genes under different conditions. Results showed that GAPDH and CYP expression remain stable throughout in the different abiotic stress treatments, CAC and PP2A expression were relatively stable under reproductive stages and α-TUB, PP2A and ETIF3E were found to be stably expressed in vegetative stages. Further, the expression levels of Diacylglycerol acyltransferase (DGAT1), a key enzyme in triacylglycerol synthesis was analyzed to confirm the validity of reference genes identified in the study. This is the first systematic study of selection of reference genes in S. hispanica, and will benefit future expression studies in this crop.

Genome-Wide Constitutively Expressed Gene Analysis and New Reference Gene Selection Based on Transcriptome Data: A Case Study from Poplar/Canker Disease Interaction

A number of transcriptome datasets for differential expression (DE) genes have been widely used for understanding organismal biology, but these datasets also contain untapped information that can be used to develop more precise analytical tools. With the use of transcriptome data generated from poplar/canker disease interaction system, we describe a methodology to identify candidate reference genes from high-throughput sequencing data. This methodology will improve the accuracy of RT-qPCR and will lead to better standards for the normalization of expression data. Expression stability analysis from xylem and phloem of Populus bejingensis inoculated with the fungal canker pathogen Botryosphaeria dothidea revealed that 729 poplar transcripts (1.11%) were stably expressed, at a threshold level of coefficient of variance (CV) of FPKM < 20% and maximum fold change (MFC) of FPKM < 2.0. Expression stability and bioinformatics analysis suggested that commonly used house-keeping (HK) genes were not the most appropriate internal controls: 70 of the 72 commonly used HK genes were not stably expressed, 45 of the 72 produced multiple isoform transcripts, and some of their reported primers produced unspecific amplicons in PCR amplification. RT-qPCR analysis to compare and evaluate the expression stability of 10 commonly used poplar HK genes and 20 of the 729 newly-identified stably expressed transcripts showed that some of the newly-identified genes (such as SSU_S8e, LSU_L5e, and 20S_PSU) had higher stability ranking than most of commonly used HK genes. Based on these results, we recommend a pipeline for deriving reference genes from transcriptome data. An appropriate candidate gene should have a unique transcript, constitutive expression, CV value of expression < 20% (or possibly 30%) and MFC value of expression <2, and an expression level of 50-1,000 units. Lastly, when four of the newly identified HK genes were used in the normalization of expression data for 20 differential expressed genes, expression analysis gave similar values to Cufflinks output. The methods described here provide an alternative pathway for the normalization of transcriptome data, a process that is essential for integrating analyses of transcriptome data across environments, laboratories, sequencing platforms, and species.

Quinone oxidoreductase 2 is involved in haustorium development of the parasitic plant Phtheirospermum japonicum

The family Orobanchaceae includes many parasitic plant species. Parasitic plants invade host vascular tissues and form organs called haustoria, which are used to obtain water and nutrients. Haustorium formation is initiated by host-derived chemicals including quinones and flavonoids. Two types of quinone oxidoreductase (QR) are involved in signal transduction leading to haustorium formation; QR1 mediates single-electron transfers and QR2 mediates 2-electron transfers. In the facultative parasite Triphysaria versicolor, QR1 is involved in haustorium induction signaling, while this role is played by QR2 in the model plant Phtheirospermum japonicum. Our results suggest that there is functional diversification in haustorium signaling molecules among different species of the Orobanchaceae.

Selection of housekeeping genes as internal controls for quantitative RT-PCR analysis of the veined rapa whelk (Rapana venosa)

Background

The veined rapa whelk Rapana venosa is an important commercial shellfish in China and quantitative real-time PCR (qRT-PCR) has become the standard method to study gene expression in R. venosa. For accurate and reliable gene expression results, qRT-PCR assays require housekeeping genes as internal controls, which display highly uniform expression in different tissues or stages of development. However, to date no studies have validated housekeeping genes in R. venosa for use as internal controls for qRT-PCR.

Methods

In this study, we selected the following 13 candidate genes for suitability as internal controls: elongation factor-1α (EF-1α), α-actin (ACT), cytochrome c oxidase subunit 1 (COX1), nicotinamide adenine dinucleotide dehydrogenase (ubiquinone) 1α subcomplex subunit 7 (NDUFA7), 60S ribosomal protein L5 (RL5), 60S ribosomal protein L28 (RL28), glyceraldehyde 3-phosphate dehydrogenase (GAPDH), β-tubulin (TUBB), 40S ribosomal protein S25 (RS25), 40S ribosomal protein S8 (RS8), ubiquitin-conjugating enzyme E2 (UBE2), histone H3 (HH3), and peptidyl-prolyl cis-trans isomerase A (PPIA). We measured the expression levels of these 13 candidate internal controls in eight different tissues and twelve larvae developmental stages by qRT-PCR. Further analysis of the expression stability of the tested genes was performed using GeNorm and RefFinder algorithms.

Results

Of the 13 candidate genes tested, we found that EF-1α was the most stable internal control gene in almost all adult tissue samples investigated with RL5 and RL28 as secondary choices. For the normalization of a single specific tissue, we suggested that EF-1α and NDUFA7 are the best combination in gonad, as well as COX1 and RL28 for intestine, EF-1α and RL5 for kidney, EF-1α and COX1 for gill, EF-1α and RL28 for Leiblein and mantle, EF-1α, RL5, and NDUFA7 for liver, GAPDH, PPIA, and RL28 for hemocyte. From a developmental perspective, we found that RL28 was the most stable gene in all developmental stages measured, and COX1 and RL5 were appropriate secondary choices. For the specific developmental stage, we recommended the following combination for normalization, PPIA, RS25, and RL28 for stage 1, RL5 and RL28 for stage 2 and 5, RL28 and NDUFA7 for stage 3, and PPIA and TUBB for stage 4.

Discussion

Our results are instrumental for the selection of appropriately validated housekeeping genes for use as internal controls for gene expression studies in adult tissues or larval development of R. venosa in the future.

Comparative physiological and transcriptomic analyses provide integrated insight into osmotic, cold, and salt stress tolerance mechanisms in banana

The growth, development, and production of banana plants are constrained by multiple abiotic stressors. However, it remains elusive for the tolerance mechanisms of banana responding to multiple abiotic stresses. In this study, we found that Fen Jiao (FJ) was more tolerant to osmotic, cold, and salt stresses than BaXi Jiao (BX) by phenotypic and physiological analyses. Comparative transcriptomic analyses highlighted stress tolerance genes that either specifically regulated in FJ or changed more than twofold in FJ relative to BX after treatments. In total, 933, 1644, and 133 stress tolerance genes were identified after osmotic, cold, and salt treatments, respectively. Further integrated analyses found that 30 tolerance genes, including transcription factor, heat shock protein, and E3 ubiquitin protein ligase, could be commonly regulated by osmotic, cold, and salt stresses. Finally, ABA and ROS signaling networks were found to be more active in FJ than in BX under osmotic, cold, and salt treatments, which may contribute to the strong stress tolerances of FJ. Together, this study provides new insights into the tolerance mechanism of banana responding to multiple stresses, thus leading to potential applications in the genetic improvement of multiple abiotic stress tolerances in banana.

Stem-Loop RT-qPCR as an Efficient Tool for the Detection and Quantification of Small RNAs in Giardia lamblia

Stem-loop quantitative reverse transcription PCR (RT-qPCR) is a molecular technique used for identification and quantification of individual small RNAs in cells. In this work, we used a Universal ProbeLibrary (UPL)-based design to detect-in a rapid, sensitive, specific, and reproducible way-the small nucleolar RNA (snoRNA) GlsR17 and its derived miRNA (miR2) of Giardia lamblia using a stem-loop RT-qPCR approach. Both small RNAs could be isolated from both total RNA and small RNA samples. Identification of the two small RNAs was carried out by sequencing the PCR-amplified small RNA products upon ligation into the pJET1.2/blunt vector. GlsR17 is constitutively expressed during the 72 h cultures of trophozoites, while the mature miR2 is present in 2-fold higher abundance during the first 48 h than at 72 h. Because it has been suggested that miRNAs in G. lamblia have an important role in the regulation of gene expression, the use of the stem-loop RT-qPCR method could be valuable for the study of miRNAs of G. lamblia. This methodology will be a powerful tool for studying gene regulation in G. lamblia, and will help to better understand the features and functions of these regulatory molecules and how they work within the RNA interference (RNAi) pathway in G. lamblia.

LbCML38 and LbRH52, two reference genes derived from RNA-Seq data suitable for assessing gene expression in Lycium barbarum L

For quantitative real-time PCR (qRT-PCR) analysis, the key prerequisite that determines result accuracy is the selection of appropriate reference gene(s). Goji (Lycium barbarum L.) is a multi-branched shrub belonging to the Solanaceae family. To date, no systematic screening or evaluation of reference gene(s) in Goji has been performed. In this work, we identified 18 candidate reference genes from the transcriptomic sequencing data of 14 samples of Goji at different developmental stages and under drought stress condition. The expression stability of these candidate genes was rigorously analyzed using qRT-PCR and four different statistical algorithms: geNorm, BestKeeper, NormFinder and RefFinder. Two novel reference genes LbCML38 and LbRH52 showed the most stable expression, whereas the traditionally used reference genes such as LbGAPDH, LbHSP90 and LbTUB showed unstable expression in the tested samples. Expression of a target gene LbMYB1 was also tested and compared using optimal reference genes LbCML38 and LbRH52, mediocre reference gene LbActin7, and poor reference gene LbHSP90 as normalization standards, respectively. As expected, calculation of the target gene expression by normalization against LbCML38, LbActin7 or LbHSP90 showed significant differences. Our findings suggest that LbCML38 and LbRH52 can be used as reference genes for gene expression analysis in Goji.

A single amino acid mutation of OsSBEIIb contributes to resistant starch accumulation in rice

Foods rich in resistant starch can help prevent various diseases, including diabetes, colon cancers, diarrhea, and chronic renal and hepatic diseases. Variations in starch biosynthesis enzymes could contribute to the high content of resistant starch in some cultivars of rice (Oryza sativa L.). Our previously published work indicated that the sbe3-rs gene in the rice mutant line, 'Jiangtangdao1' was a putative allele of the rice starch branching enzyme gene SBEIIb (previously known as SBE3); sbe3-rs might control the biosynthesis of the high resistant starch content in the rice line. Biomolecular analysis showed that the activity of SBEs was significantly lower in soluble extracts of immature seeds harvested from 'Jiangtangdao1' 15 days after flowering than in the extracts of the wild-type rice line 'Huaqingdao'. We performed gene complementation assays by introducing the wild-type OsSBEIIb into the sbe3-rs mutant 'Jiangtangdao1'. The genetically complemented lines demonstrated restored seed-related traits. The structures of endosperm amylopectin and the morphological and physicochemical properties of the starch granules in the transformants recovered to wild-type levels. This study provides evidence that sbe3-rs is a novel allele of OsSBEIIb, responsible for biosynthesis of high resistant starch in 'Jiangtangdao1'.

Structural basis of unique ligand specificity of KAI2-like protein from parasitic weed Striga hermonthica

The perception of two plant germination inducers, karrikins and strigolactones, are mediated by the proteins KAI2 and D14. Recently, KAI2-type proteins from parasitic weeds, which are possibly related to seed germination induced by strigolactone, have been classified into three clades characterized by different responses to karrikin/strigolactone. Here we characterized a karrikin-binding protein in Striga (ShKAI2iB) that belongs to intermediate-evolving KAI2 and provided the structural bases for its karrikin-binding specificity. Binding assays showed that ShKAI2iB bound karrikins but not strigolactone, differing from other KAI2 and D14. The crystal structures of ShKAI2iB and ShKAI2iB-karrikin complex revealed obvious structural differences in a helix located at the entry of its ligand-binding cavity. This results in a smaller closed pocket, which is also the major cause of ShKAI2iB’s specificity of binding karrikin. Our structural study also revealed that a few non-conserved amino acids led to the distinct ligand-binding profile of ShKAI2iB, suggesting that the evolution of KAI2 resulted in its diverse functions.

Validation of housekeeping genes as an internal control for gene expression studies in Giardia lamblia using quantitative real-time PCR

The analysis of transcript levels of specific genes is important for understanding transcriptional regulation and for the characterization of gene function. Real-time quantitative reverse transcriptase PCR (RT-qPCR) has become a powerful tool to quantify gene expression. The objective of this study was to identify reliable housekeeping genes in Giardia lamblia. Twelve genes were selected for this purpose, and their expression was analyzed in the wild type WB strain and in two strains with resistance to nitazoxanide (NTZ) and metronidazole (MTZ), respectively. RefFinder software analysis showed that the expression of the genes is different in the three strains. The integrated data from the four analyses showed that the NADH oxidase (NADH) and aldolase (ALD) genes were the most steadily expressed genes, whereas the glyceraldehyde-3-phosphate dehydrogenase gene was the most unstable. Additionally, the relative expression of seven genes were quantified in the NTZ- and MTZ-resistant strains by RT-qPCR, using the aldolase gene as the internal control, and the results showed a consistent differential pattern of expression in both strains. The housekeeping genes found in this work will facilitate the analysis of mRNA expression levels of other genes of interest in G. lamblia.

Evaluation of Sorghum [Sorghum bicolor (L.)] Reference Genes in Various Tissues and under Abiotic Stress Conditions for Quantitative Real-Time PCR Data Normalization

Accurate and reliable gene expression data from qPCR depends on stable reference gene expression for potential gene functional analyses. In this study, 15 reference genes were selected and analyzed in various sample sets including abiotic stress treatments (salt, cold, water stress, heat, and abscisic acid) and tissues (leaves, roots, seedlings, panicle, and mature seeds). Statistical tools, including geNorm, NormFinder and RefFinder, were utilized to assess the suitability of reference genes based on their stability rankings for various sample groups. For abiotic stress, PP2A and CYP were identified as the most stable genes. In contrast, EIF4α was the most stable in the tissue sample set, followed by PP2A; PP2A was the most stable in all the sample set, followed by EIF4α. GAPDH, and UBC1 were the least stably expressed in the tissue and all the sample sets. These results also indicated that the use of two candidate reference genes would be sufficient for the optimization of normalization studies. To further verify the suitability of these genes for use as reference genes, SbHSF5 and SbHSF13 gene expression levels were normalized using the most and least stable sorghum reference genes in root and water stressed-leaf tissues of five sorghum varieties. This is the first systematic study of the selection of the most stable reference genes for qPCR-related assays in Sorghum bicolor that will potentially benefit future gene expression studies in sorghum and other closely related species.

Comparative transcriptome analyses reveal core parasitism genes and suggest gene duplication and repurposing as sources of structural novelty

The origin of novel traits is recognized as an important process underlying many major evolutionary radiations. We studied the genetic basis for the evolution of haustoria, the novel feeding organs of parasitic flowering plants, using comparative transcriptome sequencing in three species of Orobanchaceae. Around 180 genes are upregulated during haustorial development following host attachment in at least two species, and these are enriched in proteases, cell wall modifying enzymes, and extracellular secretion proteins. Additionally, about 100 shared genes are upregulated in response to haustorium inducing factors prior to host attachment. Collectively, we refer to these newly identified genes as putative “parasitism genes.” Most of these parasitism genes are derived from gene duplications in a common ancestor of Orobanchaceae and Mimulus guttatus, a related nonparasitic plant. Additionally, the signature of relaxed purifying selection and/or adaptive evolution at specific sites was detected in many haustorial genes, and may play an important role in parasite evolution. Comparative analysis of gene expression patterns in parasitic and nonparasitic angiosperms suggests that parasitism genes are derived primarily from root and floral tissues, but with some genes co-opted from other tissues. Gene duplication, often taking place in a nonparasitic ancestor of Orobanchaceae, followed by regulatory neofunctionalization, was an important process in the origin of parasitic haustoria.

Selection of reference genes as internal controls for gene expression in tissues of red abalone Haliotis rufescens (Mollusca, Vetigastropoda; Swainson, 1822)

The red abalone Haliotis rufescens is one of the most important species for aquaculture in Baja California, México, and despite this, few gene expression studies have been done in tissues such as gill, head and gonad. For this purpose, reverse transcription and quantitative real time PCR (RT-qPCR) is a powerful tool for gene expression evaluation. For a reliable analysis, however, it is necessary to select and validate housekeeping genes that allow proper transcription quantification. Stability of nine housekeeping genes (ACTB, BGLU, TUBB, CY, GAPDH, HPRTI, RPL5, SDHA and UBC) was evaluated in different tissues of red abalone (gill, head and gonad/digestive gland). Four-fold serial dilutions of cDNA (from 25 ngμL(-1) to 0.39 ngμL(-1)) were used to prepare the standard curve, and it showed gene efficiencies between 0.95 and 0.99, with R(2)=0.99. geNorm and NormFinder analysis showed that RPL5 and CY were the most stable genes considering all tissues, whereas in gill HPRTI and BGLU were most stable. In gonad/digestive gland, RPL5 and TUBB were the most stable genes with geNorm, while SDHA and HPRTI were the best using NormFinder. Similarly, in head the best genes were RPL5 and UBC with geNorm, and GAPDH and CY with NormFinder. The technical variability analysis with RPL5 and abalone gonad/digestive gland tissue indicated a high repeatability with a variation coefficient within groups ≤ 0.56% and between groups ≤ 1.89%. These results will help us for further research in reproduction, thermoregulation and endocrinology in red abalone.

Identification of internal reference genes for gene expression normalization between the two sexes in dioecious white Campion

Quantitative real time (qRT)-PCR is a precise and efficient method for studying gene expression changes between two states of interest, and is frequently used for validating interesting gene expression patterns in candidate genes initially identified in genome-wide expression analyses, such as RNA-seq experiments. For an adequate normalisation of qRT-PCR data, it is essential to have reference genes available whose expression intensities are constant among the different states of interest. In this study we present and validate a catalogue of traditional and newly identified reference genes that were selected from RNA-seq data from multiple individuals from the dioecious plant Silene latifolia with the aim of studying gene expression differences between the two sexes in both reproductive and vegetative tissues. The catalogue contains more than 15 reference genes with both stable expression intensities and a range of expression intensities in flower buds and leaf tissues. These reference genes were used to normalize expression differences between reproductive and vegetative tissues in eight candidate genes with sex-biased expression. Our results suggest a trend towards a reduced sex-bias in sex-linked gene expression in vegetative tissues. In this study, we report on the systematic identification and validation of internal reference genes for adequate normalization of qRT-PCR-based analyses of gene expression differences between the two sexes in S. latifolia. We also show how RNA-seq data can be used efficiently to identify suitable reference genes in a wide diversity of species.