Selection and validation of reference genes for quantitative RT-PCR expression studies of the non-model crop Musa

Reference genes for RT-qPCR analysis in Musa acuminata genotypes contrasting in resistance to Fusarium oxysporum f. sp. cubense subtropical race 4

Banana (Musa spp.) is the most widely consumed fruit globally. Fusarium wilt, caused by Fusarium oxysporum f. sp. cubense (Foc), is a highly threatening disease to banana production. Resistance genes to Foc exist in wild Musa genotypes such as Musa acuminata subsp. burmannicoides var. Calcutta 4. Whilst real-time PCR (RT-qPCR) is appropriate for accurate analysis of gene expression in pathways involved in host defence responses, reference genes with stable expression under specific biotic stress conditions and host tissue types are necessary for normalization of sample variation. In this context, the stability in potential host reference genes ACT1, APT, EF1α, GAPDH, αTUB, RAN, UBIQ1, UBIQ2, βTUB1, βTUB3, L2 and ACTA1 was evaluated in total RNA samples from root tissues in Calcutta 4 (resistant) and Musa sp. cultivar Prata-anã (susceptible) extracted during interaction with Foc subtropical race 4 (STR4). Expression stability was calculated using the algorithms geNorm, NormFinder and BestKeeper. βTUB3 and L2 were identified as the most stable in Calcutta 4, with ACTA1 and GAPDH the most stable in Prata-anã. These reference genes for analysis of gene expression modulation in the Musa-Foc STR4 pathosystem are fundamental for advancing understanding of host defence responses to this important pathogen.

In Silico Genome-Wide Analysis of B3 Transcription Factors in Cannabis sativa L

Introduction: The B3 transcription factor has been identified in Arabidopsis thaliana, Oryza sativa, and Solanum lycopersicum, among other species. This family of transcription factors regulates seed growth, development, and stress. Cannabis is a valuable crop with numerous applications; however, no B3 transcription factors have been identified in this plant. Materials and Methods: The cannabis B3 gene family was identified and analyzed using bioinformatics analysis tools, such as the NCBI database, plantTFDB website, TBtools, and MEGA software. Quantitative real-time polymerase chain reaction (qRT-PCR) experiments were used to confirm its function. Results: The cannabis B3 family contains 65 members spread across 10 chromosomes. The isoelectric point ranged from 10.03 to 4.65, and the molecular weight ranged from 99,542.88 to 14,310.9 Da. Most of the members were found in the nucleus. The upstream promoter region of the gene contains a variety of cis-acting elements related to the stress response. RNA-seq data and qRT-PCR results showed that CsB3 genes were expressed differently in five organs of female Diku plants and in glandular hairs of nine distinct types of female cannabis inflorescences. Collinearity analysis revealed that there were more homologous genes between cannabis and dicotyledons than monocotyledonous plants, which was consistent with the evolutionary relationship. Conclusions: Hormones and external environmental factors might influence CsB3 expression. Furthermore, some genes such as CsB3-02, CsB3-07, CsB3-50, CsB3-62, and CsB3-65 may participate in cannabis growth and development and play a role in secondary metabolite synthesis. This study provides a solid foundation for further research into the gene function of the cannabis B3 family.

Plastid DNA Barcoding and RtActin cDNA Fragment Isolation of Reutealis Trisperma: A Promising Bioresource for Biodiesel Production

Reutealis trisperma belonging to the family Euphorbiaceae is currently used for biodiesel production, and rapid development in plant-based biofuel production has led to its increasing demand. However, massive utilization of bio-industrial plants has led to conservation issues. Moreover, genetic information on R trisperma is still limited, which is crucial for developmental, physiological, and molecular studies. Studying gene expression is essential to explain plant physiological processes. Nonetheless, this technique requires sensitive and precise measurement of messenger RNA (mRNA). In addition, the presence of internal control genes is important to avoid bias. Therefore, collecting and preserving genetic data for R trisperma is indispensable. In this study, we aimed to evaluate the application of plastid loci, rbcL, and matK, to the DNA barcode of R trisperma for use in conservation programs. In addition, we isolated and cloned the RtActin (RtACT) gene fragment for use in gene expression studies. Sequence information was analyzed in silico by comparison with other Euphorbiaceae plants. For actin fragment isolation, reverse-transcription polymerase chain reaction was used. Molecular cloning of RtActin was performed using the pTA2 plasmid before sequencing. We successfully isolated and cloned 592 and 840 bp of RtrbcL and RtmatK fragment genes, respectively. The RtrbcL barcoding marker, rather than the RtmatK plastidial marker, provided discriminative molecular phylogenetic data for R Trisperma. We also isolated 986 bp of RtACT gene fragments. Our phylogenetic analysis demonstrated that R trisperma is closely related to the Vernicia fordii Actin gene (97% identity). Our results suggest that RtrbcL could be further developed and used as a barcoding marker for R trisperma. Moreover, the RtACT gene could be further investigated for use in gene expression studies of plant.

Genome-Wide Identification and Expression Profiles of bZIP Genes in Cannabis sativa L

Background: The bZIP gene family plays roles in biotic and abiotic stress, secondary metabolism, and other aspects in plants. They have been reported in Arabidopsis thaliana, Oryza sativa, Artemisia annua, and other plants, but their roles in Cannabis sativa have not been determined. Materials and Methods: In this study, we analyzed the genome-wide identification and expression profile of the bZIP gene family in C. sativa. Results: A total of 51 members of the bZIP gene family were identified based on the C. sativa genome and numbered in order from CsbZIP1 to CsbZIP51. Their phylogenetic relationships, cis-elements in promoter region, gene structures and motif compositions, physicochemical properties, chromosome locations, and expression profiles, were analyzed. The results showed that the 51 CsbZIPs were unevenly distributed on 10 chromosomes and could be clustered into 11 subfamilies. Furthermore, CsbZIPs located in the same subfamilies presented similar intron/exon organization and motif composition. The expression levels of CsbZIPs in various tissues (flowers, bracts, vegetative leaves, stems, and seeds) were determined using reverse transcription quantitative polymerase chain reaction. The expression levels of CsbZIPs were higher in flowers and bracts. The 51 CsbZIPs were explored, and their structure, evolution, and expression pattern in different tissues of C. sativa were characterized synthetically. The findings indicated that CsbZIPs are essential for the growth and development of C. sativa. Conclusions: These results provide a theoretical basis for subsequent research on hemp bZIP transcription factors and the cultivation of high-cannabidiol and low-tetrahydrocannabinol high-quality cannabis varieties.

Molecular, Histological and Histochemical Responses of Banana Cultivars Challenged with Fusarium oxysporum f. sp. cubense with Different Levels of Virulence

Fusarium wilt caused by Fusarium oxysporum f. sp. cubense (Foc) is the most limiting factor in the banana agribusiness worldwide. Therefore, studies regarding pathogen attack mechanisms, and especially host defense responses, in this pathosystem are of utmost importance for genetic breeding programs in the development of Foc-resistant banana cultivars. In this study, analysis at the molecular, histological and histochemical levels of the Musa spp. x Foc interaction was performed. Three Foc isolates representative of race 1 (R1), subtropical race 4 (ST4) and isolate 229A, which is a putative ST4, were inoculated in two Prata-type cultivars (Prata-Anã and BRS Platina) and one cultivar of the Cavendish type (Grand Naine). Of seven genes related to plant-pathogen interactions, five were overexpressed in 'BRS Platina' 12 h after inoculation (HAI) with Foc R1 and ST4 but had reduced or negative expression after inoculation with Foc 229A, according to RT-qPCR analyses. While hyphae, mycelia and spores of the Foc 229A isolate grow towards the central cylinder of the Grand Naine and Prata-Anã cultivars, culminating in the occlusion of the xylem vessels, the BRS Platina cultivar responds with increased presence of cellulose, phenolic compounds and calcium oxalate crystals, reducing colonization within 30 days after inoculation (DAI). In general, these data indicate that the cultivar BRS Platina has potential for use in banana-breeding programs focused on resistance to Foc tropical race 4 (TR4) and in aggregating information on the virulence relationships of the Foc pathogen and the defense responses of banana plants after infection.

Gene Expression, Histology and Histochemistry in the Interaction between Musa sp. and Pseudocercospora fijiensis

Bananas are the main fruits responsible for feeding more than 500 million people in tropical and subtropical countries. Black Sigatoka, caused by the fungus Pseudocercospora fijiensis, is one of the most destructive disease for the crop. This fungus is mainly controlled with the use of fungicides; however, in addition to being harmful to human health, they are associated with a high cost. The development of resistant cultivars through crosses of susceptible commercial cultivars is one of the main focuses of banana breeding programs worldwide. Thus, the objective of the present study was to investigate the interaction between Musa sp. and P. fijiensis through the relative expression of candidate genes involved in the defence response to black Sigatoka in four contrasting genotypes (resistant: Calcutta 4 and Krasan Saichon; susceptible: Grand Naine and Akondro Mainty) using quantitative real-time PCR (RT-qPCR) in addition to histological and histochemical analyses to verify the defence mechanisms activated during the interaction. Differentially expressed genes (DEGs) related to the jasmonic acid and ethylene signalling pathway, GDSL-like lipases and pathogenesis-related proteins (PR-4), were identified. The number and distance between stomata were directly related to the resistance/susceptibility of each genotype. Histochemical tests showed the production of phenolic compounds and callosis as defence mechanisms activated by the resistant genotypes during the interaction process. Scanning electron microscopy (SEM) showed pathogenic structures on the leaf surface in addition to calcium oxalate crystals. The resistant genotype Krasan Saichon stood out in the analyses and has potential for use in breeding programs for resistance to black Sigatoka in banana and plantains.

Identification and Validation of Reliable Reference Genes for Gene Expression Studies in Koelreuteria paniculata

RT-qPCR is considered a rapid and reliable technique for analyzing gene expression. This technique is commonly used to analyze the expression of various genes at diverse transcriptional levels in different samples. However, few studies have characterized ornamental Koelreuteria species for reliable reference genes. In this study, eight reference genes were evaluated as controls in RT-qPCR with SYBR green to quantify gene expression in different Koelreuteria paniculata samples. All selected reference genes showed a broad range of Ct values in all samples, which was supportive of their variable expression. Our results showed significant variation in the stable expression of K. paniculata genes. Sample data, analyzed using geNorm, NormFinder, and BestKeeper, showed that phospholipase (PLA2) and β-actin (ACT) were the most suitable and statistically reliable reference genes, whereas ribosomal protein L13 (RPL13) and elongation factor 1-α (EF1α) were less stable and unsuitable for use as internal controls. To compare gene expression levels, two or more reference genes should be used for data normalization. Thus, the stability and expression of both PLA2 and ACT were believed to provide better normalization and quantification of the transcript levels for gene expression studies in K. paniculata.

Identification of Suitable Reference Genes for qRT-PCR Normalization in Kiwifruit

Reference genes are used for the correction of qRT-PCR data, and it is necessary to investigate the optimum reference gene under certain conditions. The expression levels of seven traditional reference genes ACT1, ACT2, GAPDH, 18S rRNA, UBQ, TUB and CYP were analyzed using qRT-PCR in different varieties, tissues, developmental stages and hormone (or pollen polysaccharide) treatments in kiwifruit. Gene expression stability was assessed with the help of three common software (geNorm, NormFinder, BestKeeper), and the minimum number of reference genes necessary for normalization was also determined. GAPDH, ACT1 and ACT2 were selected as reference genes for different genotypes of kiwifruit. GAPDH and UBQ were the best combinations of reference genes for root, stem, leaf, flower and fruit. GAPDH and ACT1 could be the preferred reference genes for normalization of qRT-PCR data during fruit development. The pairing of ACT1 and UBQ constituted the optimal combination of reference genes in kiwifruit treated with different hormones (or pollen polysaccharide). This study provides a new and reliable option for the use of reference genes in the analysis of gene expression patterns of interest in kiwifruit.

Evaluation of Reference Genes in the Polyploid Complex Dianthus broteri (Caryophyllaceae) Using qPCR

Dianthus broteri is an endemic complex which is considered the largest polyploid series within the Dianthus genus. This polyploid species involves four cytotypes (2×, 4×, 6× and 12×) with spatial and ecological segregation. The study of gene expression in polyploid species must be very rigorous because of the effects of duplications on gene regulation. In these cases, real-time polymerase chain reaction (qPCR) is the most appropriate technique for determining the gene expression profile because of its high sensitivity. The relative quantification strategy using qPCR requires genes with stable expression, known as reference genes, for normalization. In this work, we evaluated the stability of 13 candidate genes to be considered reference genes in leaf and petal tissues in Dianthus broteri. Several statistical analyses were used to determine the most stable candidate genes: Bayesian analysis, network analysis based on equivalence tests, geNorm and BestKeeper algorithms. In the leaf tissue, the most stable candidate genes were TIP41, TIF5A, PP2A and SAMDC. Similarly, the most adequate reference genes were H3.1, TIP41, TIF5A and ACT7 in the petal tissue. Therefore, we suggest that the best reference genes to compare different ploidy levels for both tissues in D. broteri are TIP41 and TIF5A.

Candidate reference genes for quantitative gene expression analysis in Lagerstroemia indica

Quantitative gene expression analysis by qPCR requires reference genes for normalization. Lagerstroemia indica (crape myrtle) is a popular ornamental plant in the world, but suitable endogenous reference genes are lacking. To find suitable reference genes, we evaluated the stabilities of nine candidate genes in six experimental data sets: six different tissues, three leaf colors, nine flower colors, and under three abiotic stresses (salt, drought, cold) using four statistical algorithms. A target gene LiMYB56 (homolog of Arabidopsis MYB56) was used to verify the authenticity and accuracy of the candidate reference genes. The results showed that the combination of two stably expressed reference genes, rather than a single reference gene, improved the accuracy of the qPCR. LiEF1α-2 + LiEF1α-3 was best for the tissue, salt treatment, and drought treatment sets; LiEF1α-2 + LiEF1α-1 was optimal for leaf color; LiEF1α-2 + LiACT7 was optimal for cold treatment; and LiUBC + LiEF1α-1 was best for the flower color set. Notably, LiEF1α-2 had high expression stability in all six experimental sets, implying it may be a good reference gene for expression studies in L. indica. Our results will facilitate future gene expression studies in L. indica.

Identification and functional analysis of the CorA/MGT/MRS2-type magnesium transporter in banana

Magnesium (Mg) plays an irreplaceable role in plant growth and development. Mg transporters, especially CorA/MGT/MRS2 family proteins, played a vital role in regulating Mg content in plant cells. Although extensive work has been conducted in model crops, such as Arabidopsis, rice, and maize, the relevant information is scarce in tropical crops. In this study, 10 MaMRS2 genes in banana (Musa acuminata) were isolated from its genome and classified into five distinct clades. The putative physiochemical properties, chromosome location, gene structure, cis-acting elements, and duplication relationships in between these members were analyzed. Complementary experiments revealed that three MaMRS2 gene members (MaMRS2-1, MaMRS2-4, MaMRS2-7), from three distinct phylogenetic branches, were capable of restoring the function of Mg transport in Salmonella typhimurium mutants. Semi-quantitative RT-PCR showed that MaMRS2 genes were differentially expressed in banana cultivar ‘Baxijiao’ (Musa spp. AAA Cavendish) seedlings. The result was confirmed by real-time PCR analysis, in addition to tissue specific expression, expression differences among MaMRS2 members were also observed under Mg deficiency conditions. These results showed that Mg transporters may play a versatile role in banana growth and development, and our work will shed light on the functional analysis of Mg transporters in banana.

Differential regulation of the banana stress NAC family by individual and combined stresses of drought and heat in susceptible and resistant genotypes

Banana, an important tropical fruit crop, often faces drought, heat and its combination during its growth, leading to decreased yields. The combined stresses caused 100% yield loss in Grand Nain (GN) as compared to only 46% in Hill Banana (HB). To understand the response of combined stresses, we studied the stress-responsive NAC gene sub-family under individual and combined drought/heat stresses under controlled and field conditions in the stress-sensitive GN (AAA genotype) and stress-tolerant HB (AAB genotype). Under drought, expression of most stress-NACs increased with progression of drought in either one or the other genotype with little overlap. Heat stress caused a continuous decline in expression of most genes in HB unlike in GN where many NACs were up-regulated although to a lesser scale than for drought. Combination of the two stresses elicited a very different response compared with individual stresses. GN responded strongly to the combined stress with up-regulation of most genes unlike that seen in drought. Surprisingly, NAC genes in HB did not respond much to the more severe combination of the stresses despite being up-regulated strongly by drought. The response of the NACs to combined field stress was similar to that under controlled conditions. Most of the stress-NACs were strongly up-regulated upon treatment with exogenous ABA within 30-60 min, the increase being more prominent in GN. The studies suggest that the B genome in the stress-tolerant HB may counter more drastic combined stresses without taking recourse to the expression of stress NACs.

Stable reference genes for RT-qPCR analysis of gene expression in the Musa acuminata-Pseudocercospora musae interaction

Leaf pathogens are limiting factors in banana (Musa spp.) production, with Pseudocercospora spp. responsible for the important Sigatoka disease complex. In order to investigate cellular processes and genes involved in host defence responses, quantitative real-time PCR (RT-qPCR) is an analytical technique for gene expression quantification. Reliable RT-qPCR data, however, requires that reference genes for normalization of mRNA levels in samples are validated under the conditions employed for expression analysis of target genes. We evaluated the stability of potential reference genes ACT1, α-TUB, UBQ1, UBQ2, GAPDH, EF1α, APT and RAN. Total RNA was extracted from leaf tissues of Musa acuminata genotypes Calcutta 4 (resistant) and Cavendish Grande Naine (susceptible), both subjected to P. musae infection. Expression stability was determined with NormFinder, BestKeeper, geNorm and RefFinder algorithms. UBQ2 and RAN were the most stable across all M. acuminata samples, whereas when considering inoculated and non-inoculated leaf samples, APT and UBQ2 were appropriate for normalization in Calcutta 4, with RAN and α-TUB most stable in Cavendish Grande Naine. This first study of reference genes for relative quantification of target gene expression in the M. acuminata-P. musae interaction will enable reliable analysis of gene expression in this pathosystem, benefiting elucidation of disease resistance mechanisms.

Effect of paleopolyploidy and allopolyploidy on gene expression in banana

Background

Bananas (Musa spp.) are an important crop worldwide. Most modern cultivars resulted from a complex polyploidization history that comprised three whole genome duplications (WGDs) shaping the haploid Musa genome, followed by inter- and intra-specific crosses between Musa acuminata and M. balbisiana (A and B genome, respectively). Unresolved hybridizations finally led to banana diversification into several autotriploid (AAA) and allotriploid cultivars (AAB and ABB). Using transcriptomic data, we investigated the impact of the genome structure on gene expression patterns in roots of 12 different triploid genotypes covering AAA, AAB and ABB subgenome constitutions.

Results

We demonstrate that (i) there are different genome structures, (ii) expression patterns go beyond the predicted genomic groups, and (iii) the proportion of the B genome influences the gene expression.

The presence of the B genome is associated with a higher expression of genes involved in flavonoid biosynthesis, fatty acid metabolism, amino sugar and nucleotide sugar metabolism and oxidative phosphorylation. There are cultivar-specific chromosome regions with biased B:A gene expression ratios that demonstrate homoeologous exchanges (HE) between A and B sub-genomes. In two cultivars, aneuploidy was detected. We identified 3674 genes with a different expression level between allotriploid and autotriploid with ~ 57% having recently duplicated copies (paralogous). We propose a Paralog Inclusive Expression (PIE) analysis that appears to be suitable for genomes still in a downsizing and fractionation process following whole genome duplications. Our approach allows highlighting the genes with a maximum likelihood to affect the plant phenotype.

Conclusions

This study on banana is a good case to investigate the effects of alloploidy in crops. We conclude that allopolyploidy triggered changes in the genome structure of a crop and it clearly influences the gene.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5618-0) contains supplementary material, which is available to authorized users.

Comprehensive evaluation of candidate reference genes for real-time quantitative PCR (RT-qPCR) data normalization in nutri-cereal finger millet [Eleusine Coracana (L.)]

Finger millet (Eleusine coracana L.) is an annual herbaceous self-pollinating C4 cereal crop of the arid and semi-arid regions of the world. Finger millet is a food security crop proven to have resilience to changing climate and scores very high in nutrition. In the current study, we have assessed sixteen candidate reference genes for their appropriateness for the normalization studies in finger millet subjected to experimental regimes and treatments. Ten candidate reference genes (GAPDH, β-TUB, CYP, EIF4α, TIP41, UBC, G6PD, S24, MACP and MDH) were cloned and six (ACT, ELF1α, PP2A, PT, S21 and TFIID) were mined from the NCBI database as well as from the literature. Expression stability ranking of the finger millet reference genes was validated using four different statistical tools i.e., geNorm, NormFinder, BestKeeper, ΔCt and RefFinder. From the study, we endorse MACP, CYP, EIF4α to be most stable candidate reference genes in all 'tissues', whereas PT, TFIID, MACP ranked high across genotypes, β-TUB, CYP, ELF1α were found to be best under abiotic stresses and 'all samples set'. The study recommends using minimum of two reference genes for RT-qPCR data normalizations in finger millet. All in all, CYP, β-TUB, and EF1α, in combination were found to be best for robust normalizations under most experimental conditions. The best and the least stable genes were validated for confirmation by assessing their appropriateness for normalization studies using EcNAC1 gene. The report provides the first comprehensive list of suitable stable candidate reference genes for nutritional rich cereal finger millet that will be advantageous to gene expression studies in this crop.

Selection and Evaluation of Appropriate Reference Genes for RT-qPCR Normalization of Volvariella volvacea Gene Expression under Different Conditions

Volvariella volvacea (V. volvacea), commonly referred to as Chinese (paddy straw) mushroom, is a basidiomycete with a protein-rich volva and pileus. Selecting appropriate reference genes is a crucial step in the normalization of quantitative real-time PCR data. Therefore, 12 candidate reference genes were selected from the V. volvacea transcriptome based on previous studies and then BestKeeper, geNorm, and NormFinder were used to identify reference genes stably expressed during different developmental stages and conditions. Of the 12 candidate reference genes, SPRY domain protein (SPRYp), alpha-tubulin (TUBα), cyclophilin (CYP), L-asparaginase (L-asp), and MSF1-domain-containing protein (MSF1) were the most stably expressed under different experimental conditions, while 18S ribosomal RNA (18S), 28S ribosomal RNA (28S), and beta-actin (ACTB) were the least stably expressed. This investigation not only revealed potential factors influencing the suitability of reference genes, but also identified optimal reference genes from a pool of candidate genes under a wide range of conditions.

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.

A banana NAC transcription factor (MusaSNAC1) impart drought tolerance by modulating stomatal closure and H2O2 content

MusaSNAC1 function in H₂O₂ mediated stomatal closure and promote drought tolerance by directly binding to CGT[A/G] motif in regulatory region of multiple stress-related genes. Drought is a abiotic stress-condition, causing reduced plant growth and diminished crop yield. Guard cells of the stomata control photosynthesis and transpiration by regulating CO₂ exchange and water loss, thus affecting growth and crop yield. Roles of NAC (NAM, ATAF1/2 and CUC2) protein in regulation of stress-conditions has been well documented however, their control over stomatal aperture is largely unknown. In this study we report a banana NAC protein, MusaSNAC1 which induced stomatal closure by elevating H₂O₂ content in guard cells during drought stress. Overexpression of MusaSNAC1 in banana resulted in higher number of stomata closure causing reduced water loss and thus elevated drought-tolerance. During drought, expression of GUS (β-glucuronidase) under P _MusaSNAC1 was remarkably elevated in guard cells of stomata which correlated with its function as a transcription factor regulating stomatal aperture closing. MusaSNAC1 is a transcriptional activator belonging to SNAC subgroup and its 5'-upstream region contain multiple Dof1 elements as well as stress-associated cis-elements. Moreover, MusaSNAC1 also regulate multiple stress-related genes by binding to core site of NAC-proteins CGT[A/G] in their 5'-upstream region. Results indicated an interesting mechanism of drought tolerance through stomatal closure by H₂O₂ generation in guard cells, regulated by a NAC-protein in banana.

Selection and validation of reference genes for qRT-PCR analysis of gene expression in Microsporum canis growing under different adhesion-inducing conditions

Dermatophytes are the group of filamentous fungi infecting keratinized structures such as skin, hair, and nails. Knowledge about genes and molecular mechanisms responsible for pathogenicity, as well as other biological properties of Microsporum canis is still relatively poor. The qRT-PCR is a reliable technique for quantifying gene expression across various biological processes, and choosing a set of suitable reference genes to normalize the expression data is a crucial step of this technique. We investigated the suitability of nine candidate reference genes: β-act, β-tub, adp-rf, ef1-α, sdha, rpl2, mbp1, psm1, and rGTPa for gene expression analysis in the dermatophyte M. canis in response to different carbon sources, phosphate levels, and pH shifts - factors that are extremely important and necessary for growth of dermatophyte in the host tissue. The transcription stability of these genes was evaluated using NormFinder, geNorm, BestKeeper, and RefFinder software. Regarding expression stability, mbp1, β-act, and sdha were the most stable housekeeping genes which we recommend for future qRT-PCR studies on M. canis strains. To the best of our knowledge this is the first study on selection and validation of reference genes for qRT-PCR data normalization in M. canis growth in culture media which promote adhesion-inducing conditions.

Reference genes identification for normalization of qPCR under multiple stresses in Hordeum brevisubulatum

BACKGROUND

Real-time quantitative PCR has been widely used as the most reliable method to measure gene expression, due to its high accuracy and specificity. Wild barley (Hordeum brevisubulatum (Trin.) Link) is a wild relative species in Triticeae that has strong tolerance to abiotic stresses and extremely wide adaptation. However, suitable references gene have not been documented for standardization of gene expression in wild barley under abiotic stress.

RESULTS

Here we report the first systematic and comprehensive analysis of reference genes for quantitative real-time PCR standardization in wild barley. We selected 11 genes, including ACT (Actin), ADP (ADP-ribosylation factor 1), CYP2 (Cyclophilin 2), EF-1α (Elongation factor 1-alpha), GAPDH (Glyceraldehyde 3-phosphate dehydrogenase), HSP90 (Heat shock protein 90), TUBα (Alpha-tubulin), TUBβ6 (Beta-tubulin 6), UBI (Ubiquitin), 18SrRNA-1 (guanine1575-N7-methyltransferase) and 18SrRNA-3 (adenine1779-N6-dimethyltransferase) from a wild barley transcriptome database and analyzed their expression stabilities in shoots and roots of wild barley seedling under various stress conditions using comparative ΔCt, BestKeeper, Normfinder and geNorm software. The results demonstrated that ADP was the most suitable reference gene in salt stress while UBI showed peak stability under mannitol and ABA stress; EF-1α was the most appropriate reference gene for PEG, GA₃, ethylene and heat stress; 18SrRNA-3 was the best choice for cold stress; and TUBα was the first stable gene across different tissues.

CONCLUSIONS

Our main contribution was to identify reference genes with suitable and stable expression in wild barley under various stress conditions and in different tissues to provide a useful resource for future studies. The results demonstrate the importance of transcriptome data as a useful resource for the screening of candidate reference genes and highlight the need for specific reference genes for specific conditions. Furthermore, these findings will provide valuable information for wild barley and relative species for future research.

Validation of Reference Genes in mRNA Expression Analysis Applied to the Study of Asthma

The quantitative Polymerase Chain Reaction is the most used technique for the study of gene expression. To correct putative experimental errors of this technique is necessary normalizing the expression results of the gene of interest with the obtained for reference genes. Here, we describe an example of the process to select reference genes. In this particular case, we select reference genes for expression studies in the peripheral blood mononuclear cells of asthmatic patients.

Evaluation of reference genes for reverse transcription quantitative real-time PCR (RT-qPCR) studies in Silene vulgaris considering the method of cDNA preparation

Accurate gene expression measurements are essential in studies of both crop and wild plants. Reverse transcription quantitative real-time PCR (RT-qPCR) has become a preferred tool for gene expression estimation. A selection of suitable reference genes for the normalization of transcript levels is an essential prerequisite of accurate RT-qPCR results. We evaluated the expression stability of eight candidate reference genes across roots, leaves, flower buds and pollen of Silene vulgaris (bladder campion), a model plant for the study of gynodioecy. As random priming of cDNA is recommended for the study of organellar transcripts and poly(A) selection is indicated for nuclear transcripts, we estimated gene expression with both random-primed and oligo(dT)-primed cDNA. Accordingly, we determined reference genes that perform well with oligo(dT)- and random-primed cDNA, making it possible to estimate levels of nucleus-derived transcripts in the same cDNA samples as used for organellar transcripts, a key benefit in studies of cyto-nuclear interactions. Gene expression variance was estimated by RefFinder, which integrates four different analytical tools. The SvACT and SvGAPDH genes were the most stable candidates across various organs of S. vulgaris, regardless of whether pollen was included or not.

Banana NAC transcription factor MusaNAC042 is positively associated with drought and salinity tolerance

Banana is an important fruit crop and its yield is hampered by multiple abiotic stress conditions encountered during its growth. The NAC (NAM, ATAF, and CUC) transcription factors are involved in plant response to biotic and abiotic stresses. In the present study, we studied the induction of banana NAC042 transcription factor in drought and high salinity conditions and its overexpression in transgenic banana to improve drought and salinity tolerance. MusaNAC042 expression was positively associated with stress conditions like salinity and drought and it encoded a nuclear localized protein. Transgenic lines of banana cultivar Rasthali overexpressing MusaNAC042 were generated by Agrobacterium-mediated transformation of banana embryogenic cells and T-DNA insertion was confirmed by PCR and Southern blot analysis. Our results using leaf disc assay indicated that transgenic banana lines were able to tolerate drought and high salinity stress better than the control plants and retained higher level of total chlorophyll and lower level of MDA content (malondialdehyde). Transgenic lines analyzed for salinity (250 mM NaCl) and drought (Soil gravimetric water content 0.15) tolerance showed higher proline content, better Fv/Fm ratio, and lower levels of MDA content than control suggesting that MusaNAC042 may be involved in responses to higher salinity and drought stresses in banana. Expression of several abiotic stress-related genes like those coding for CBF/DREB, LEA, and WRKY factors was altered in transgenic lines indicating that MusaNAC042 is an efficient modulator of abiotic stress response in banana.

Histone Acetylation Modifications Affect Tissue-Dependent Expression of Poplar Homologs of C4 Photosynthetic Enzyme Genes

Histone modifications play important roles in regulating the expression of C₄ photosynthetic genes. Given that all enzymes required for the C₄ photosynthesis pathway are present in C₃ plants, it has been hypothesized that this expression regulatory mechanism has been conserved. However, the relationship between histone modification and the expression of homologs of C₄ photosynthetic enzyme genes has not been well determined in C₃ plants. In the present study, we cloned nine hybrid poplar (Populus simonii × Populus nigra) homologs of maize (Zea mays) C₄ photosynthetic enzyme genes, carbonic anhydrase (CA), pyruvate orthophosphate dikinase (PPDK), phosphoenolpyruvate carboxykinase (PCK), and phosphoenolpyruvate carboxylase (PEPC), and investigated the correlation between the expression levels of these genes and the levels of promoter histone acetylation modifications in four vegetative tissues. We found that poplar homologs of C₄ homologous genes had tissue-dependent expression patterns that were mostly well-correlated with the level of histone acetylation modification (H3K9ac and H4K5ac) determined by chromatin immunoprecipitation assays. Treatment with the histone deacetylase inhibitor trichostatin A further confirmed the role of histone acetylation in the regulation of the nine target genes. Collectively, these results suggest that both H3K9ac and H4K5ac positively regulate the tissue-dependent expression pattern of the PsnCAs, PsnPPDKs, PsnPCKs, and PsnPEPCs genes and that this regulatory mechanism seems to be conserved among the C₃ and C₄ species. Our findings provide new insight that will aid efforts to modify the expression pattern of these homologs of C₄ genes to engineer C₄ plants from C₃ plants.

Identification of Differentially-Expressed Genes in Response to Mycosphaerella fijiensis in the Resistant Musa Accession 'Calcutta-4' Using Suppression Subtractive Hybridization

Bananas and plantains are considered an important crop around the world. Banana production is affected by several constraints, of which Black Sigatoka Disease, caused by the fungus Mycosphaerella fijiensis, is considered one of the most important diseases in banana plantations. The banana accession ‘Calcutta-4’ has a natural resistance to Black Sigatoka; however, the fruit is not valuable for commercialization. Gene identification and expression studies in ‘Calcutta-4’ might reveal possible gene candidates for resistant to the disease and elucidate mechanisms for resistance. A subtracted cDNA library was generated from leaves after 6, 9 and 12 days inoculated with M. fijiensis conidia on greenhouse banana plants of the accession ‘Calcutta-4’. Bioinformatic analysis revealed 99 good quality sequences. Blast2go analysis revealed that 31% of the sequences could not be categorized and, according to the Biological Process Category, 32 and 28 ESTs are related to general metabolic and cellular processes, respectively; while 10 ESTs response to stimulus. Seven sequences were redundant and one was similar to genes that may be involved in pathogen resistance including the putative disease resistance protein RGA1. Genes encoding zinc finger domains were identified and may play an important role in pathogen resistance by inducing the expression of downstream genes. Expression analysis of four selected genes was performed using RT-qPCR during the early stage of the disease development at 6, 9, 12 and 15 days post inoculation showing a peak of up regulation at 9 or 12 days post inoculation. Three of the four genes showed an up-regulation of expression in ‘Calcutta-4’ when compared to ‘Williams’ after inoculation with M. fijiensis, suggesting a fine regulation of specific gene candidates that may lead to a resistance response. The genes identified in early responses in a plant-pathogen interaction may be relevant for the resistance response of ‘Calcutta-4’ to Black Sigatoka. Genes with different functions may play a role in plant response to the disease. The present study suggests a fine up regulation of these genes that might be needed to perform an incompatible interaction. Further gene functional studies need to be performed to validate their use as candidate resistance genes in susceptible banana cultivars.

Screening and Validation of Housekeeping Genes of the Root and Cotyledon of Cunninghamia lanceolata under Abiotic Stresses by Using Quantitative Real-Time PCR

Cunninghamia lanceolata (Chinese fir) is a fast-growing and commercially important conifer of the Cupressaceae family. Due to the unavailability of complete genome sequences and relatively poor genetic background information of the Chinese fir, it is necessary to identify and analyze the expression levels of suitable housekeeping genes (HKGs) as internal reference for precise analysis. Based on the results of database analysis and transcriptome sequencing, we have chosen five candidate HKGs (Actin, GAPDH, EF1a, 18S rRNA, and UBQ) with conservative sequences in the Chinese fir and related species for quantitative analysis. The expression levels of these HKGs in roots and cotyledons under five different abiotic stresses in different time intervals were measured by qRT-PCR. The data were statistically analyzed using the following algorithms: NormFinder, BestKeeper, and geNorm. Finally, RankAggreg was applied to merge the sequences generated from three programs and rank these according to consensus sequences. The expression levels of these HKGs showed variable stabilities under different abiotic stresses. Among these, Actin was the most stable internal control in root, and GAPDH was the most stable housekeeping gene in cotyledon. We have also described an experimental procedure for selecting HKGs based on the de novo sequencing database of other non-model plants.

Analysis of gene expression and histone modification between C4 and non-C4 homologous genes of PPDK and PCK in maize

More efficient photosynthesis has allowed C4 plants to adapt to more diverse ecosystems (such as hot and arid conditions) than C3 plants. To better understand C4 photosynthesis, we investigated the expression patterns of C4 genes (C4PPDK and PCK1) and their non-C4 homologous genes (CyPPDK1, CyPPDK2, and PCK2) in the different organs of maize (Zea mays). Both C4 genes and non-C4 genes showed organ-dependent expression patterns. The mRNA levels of C4 genes were more abundant in leaf organ than in seeds at 25 days after pollination (DAP), while non-C4 genes were mainly expressed in developing seeds. Further, acetylation of histone H3 lysine 9 (H3K9ac) positively correlates with mRNA levels of C4 genes (C4PPDK and PCK1) in roots, stems, leaves, and seeds at 25 DAP, acetylation of histone H4 lysine 5 (H4K5ac) in the promoter regions of both C4 (C4PPDK and PCK1) and non-C4 genes (CyPPDK1, CyPPDK2, and PCK2) correlated well with their transcripts abundance in stems. In photosynthetic organs (stems and leaves), dimethylation of histone H3 lysine 9 (H3K9me2) negatively correlated with mRNA levels of both C4 and non-C4 genes. Taken together, our data suggest that histone modification was involved in the transcription regulation of both C4 genes and non-C4 genes, which might provide a clue of the functional evolution of C4 genes.

Identification of Reference Genes for Quantitative Gene Expression Studies in a Non-Model Tree Pistachio (Pistacia vera L.)

The tree species, Pistacia vera (P. vera) is an important commercial product that is salt-tolerant and long-lived, with a possible lifespan of over one thousand years. Gene expression analysis is an efficient method to explore the possible regulatory mechanisms underlying these characteristics. Therefore, having the most suitable set of reference genes is required for transcript level normalization under different conditions in P. vera. In the present study, we selected eight widely used reference genes, ACT, EF1α, α-TUB, β-TUB, GAPDH, CYP2, UBQ10, and 18S rRNA. Using qRT-PCR their expression was assessed in 54 different samples of three cultivars of P. vera. The samples were collected from different organs under various abiotic treatments (cold, drought, and salt) across three time points. Several statistical programs (geNorm, NormFinder, and BestKeeper) were applied to estimate the expression stability of candidate reference genes. Results obtained from the statistical analysis were then exposed to Rank aggregation package to generate a consensus gene rank. Based on our results, EF1α was found to be the superior reference gene in all samples under all abiotic treatments. In addition to EF1α, ACT and β-TUB were the second best reference genes for gene expression analysis in leaf and root. We recommended β-TUB as the second most stable gene for samples under the cold and drought treatments, while ACT holds the same position in samples analyzed under salt treatment. This report will benefit future research on the expression profiling of P. vera and other members of the Anacardiaceae family.

Defense Gene Expression Associated with Biotrophic Phase of Mycosphaerella fijiensis M. Morelet Infection in Banana

Banana black leaf streak, caused by Mycosphaerella fijiensis M. Morelet, is a primary phytosanitary concern in export of this fruit around the world. To develop improved cultivars resistant to this disease, an understanding of host response to infection is necessary. In this study, we obtained expression data on 14,872 genes by microarray analysis in the resistant genotype Musa acuminata subsp. burmannicoides 'Calcutta 4' after inoculation with Mycosphaerella fijiensis. From these data, 16 genes were analyzed as potential reference genes and 12 genes were identified as potential early indicators of the onset of the host defense response. Subsequently, these genes were analyzed by quantitative reverse-transcription polymerase chain reaction in susceptible 'Williams' and resistant Calcutta 4. The 18S and 26S ribosomal subunit genes in both cultivars showed the best characteristics as reference genes. In all, 5 of the 12 defense genes expressed shortly after infection (peroxidase, pathogenesis-related [PR]-4, PR-10, phenylalanine ammonia-liase, and disease resistance response 1) showed overexpression in Calcutta 4 between 6 and 24 h after inoculation as opposed to Williams, which did not show overexpression after 144 h. Early induction of defense-related genes could be a key component of the resistance of the Calcutta 4 genotype against M. fijiensis. In addition, these five genes could be used as indicators of the activation of defense responses in the interaction between banana and M. fijiensis.

Selection of Reference Genes for Gene Expression Normalization in Peucedanum praeruptorum Dunn under Abiotic Stresses, Hormone Treatments and Different Tissues

Peucedanum praeruptorum Dunn is one of the main traditional Chinese medicines producing coumarins and plenty of literatures are focused on the biosynthesis of coumarins. Quantitative real-time reverse transcription PCR (qRT-PCR) is a widely used method in studying the biosynthesis pathway and the selection of reference genes plays a crucial role in accurate normalization. To facilitate biosynthesis study of coumarins, twelve candidate reference genes were selected from the transcriptome database of P. praeruptorum according to previous studies. Then, BestKeeper, geNoFrm and NormFinder were used for selecting stably expressed reference genes in different tissues and under various stress treatments. The results indicated that, among the twelve candidate reference genes, the SAND family protein (SAND), actin 2 (ACT2), ubiquitin-conjugating enzyme 9 (UBC9), protein phosphatase 2A gene (PP2A) and polypyrimidine tract-binding protein (PTBP1) were the most stable reference genes under different experimental treatments, while glyceraldehyde 3-phosphate dehydrogenase (GAPDH) and tubulin beta-6 (TUB6) were the least stable genes. In addition, the suitability of SAND, TIP41-like protein (TIP41), UBC9, ACT2, TUB6 and their combination as reference genes were confirmed by normalizing the expression of 1-aminocyclopropane-1-carboxylate oxidase (ACO) in different treatments. This work is the first survey of the stability of reference genes in P. praeruptorum and provides guidelines to obtain more accurate qRT-PCR results in P. praeruptorum and other plant species.

Differential root transcriptomics in a polyploid non-model crop: the importance of respiration during osmotic stress

To explore the transcriptomic global response to osmotic stress in roots, 18 mRNA-seq libraries were generated from three triploid banana genotypes grown under mild osmotic stress (5% PEG) and control conditions. Illumina sequencing produced 568 million high quality reads, of which 70-84% were mapped to the banana diploid reference genome. Using different uni- and multivariate statistics, 92 genes were commonly identified as differentially expressed in the three genotypes. Using our in house workflow to analyze GO enriched and underlying biochemical pathways, we present the general processes affected by mild osmotic stress in the root and focus subsequently on the most significantly overrepresented classes associated with: respiration, glycolysis and fermentation. We hypothesize that in fast growing and oxygen demanding tissues, mild osmotic stress leads to a lower energy level, which induces a metabolic shift towards (i) a higher oxidative respiration, (ii) alternative respiration and (iii) fermentation. To confirm the mRNA-seq results, a subset of twenty up-regulated transcripts were further analysed by RT-qPCR in an independent experiment at three different time points. The identification and annotation of this set of genes provides a valuable resource to understand the importance of energy sensing during mild osmotic stress.

Introduction to the Gene Expression Analysis

In 1941, Beadle and Tatum published experiments that would explain the basis of the central dogma of molecular biology, whereby the DNA through an intermediate molecule, called RNA, results proteins that perform the functions in cells. Currently, biomedical research attempts to explain the mechanisms by which develops a particular disease, for this reason, gene expression studies have proven to be a great resource. Strictly, the term "gene expression" comprises from the gene activation until the mature protein is located in its corresponding compartment to perform its function and contribute to the expression of the phenotype of cell.The expression studies are directed to detect and quantify messenger RNA (mRNA) levels of a specific gene. The development of the RNA-based gene expression studies began with the Northern Blot by Alwine et al. in 1977. In 1969, Gall and Pardue and John et al. independently developed the in situ hybridization, but this technique was not employed to detect mRNA until 1986 by Coghlan. Today, many of the techniques for quantification of RNA are deprecated because other new techniques provide more information. Currently the most widely used techniques are qPCR, expression microarrays, and RNAseq for the transcriptome analysis. In this chapter, these techniques will be reviewed.

Reference Gene Selection for qPCR Normalization of Kosteletzkya virginica under Salt Stress

Kosteletzkya virginica (L.) is a newly introduced perennial halophytic plant. Presently, reverse transcription quantitative real-time PCR (qPCR) is regarded as the best choice for analyzing gene expression and its accuracy mainly depends on the reference genes which are used for gene expression normalization. In this study, we employed qPCR to select the most stable reference gene in K. virginica which showed stable expression profiles under our experimental conditions. The candidate reference genes were 18S ribosomal RNA (18SrRNA), β-actin (ACT), α-tubulin (TUA), and elongation factor (EF). We tracked the gene expression profiles of the candidate genes and analyzed their stabilities through BestKeeper, geNorm, and NormFinder software programs. The results of the three programs were identical and 18SrRNA was assessed to be the most stable reference gene in this study. However, TUA was identified to be the most unstable. Our study proved again that the traditional reference genes indeed displayed a certain degree of variations under given experimental conditions. Importantly, our research also provides guidance for selecting most suitable reference genes and lays the foundation for further studies in K. virginica.

Constitutive and stress-inducible overexpression of a native aquaporin gene (MusaPIP2;6) in transgenic banana plants signals its pivotal role in salt tolerance

High soil salinity constitutes a major abiotic stress and an important limiting factor in cultivation of crop plants worldwide. Here, we report the identification and characterization of a aquaporin gene, MusaPIP2;6 which is involved in salt stress signaling in banana. MusaPIP2;6 was firstly identified based on comparative analysis of stressed and non-stressed banana tissue derived EST data sets and later overexpression in transgenic banana plants was performed to study its tangible functions in banana plants. The overexpression of MusaPIP2;6 in transgenic banana plants using constitutive or inducible promoter led to higher salt tolerance as compared to equivalent untransformed control plants. Cellular localization assay performed using transiently transformed onion peel cells indicated that MusaPIP2;6 protein tagged with green fluorescent protein was translocated to the plasma membrane. MusaPIP2;6-overexpressing banana plants displayed better photosynthetic efficiency and lower membrane damage under salt stress conditions. Our results suggest that MusaPIP2;6 is involved in salt stress signaling and tolerance in banana.

Reference genes selection and normalization of oxidative stress responsive genes upon different temperature stress conditions in Hypericum perforatum L

Reverse transcription-quantitative real-time PCR (RT-qPCR) is a widely used technique for gene expression analysis. The reliability of this method depends largely on the suitable selection of stable reference genes for accurate data normalization. Hypericum perforatum L. (St. John's wort) is a field growing plant that is frequently exposed to a variety of adverse environmental stresses that can negatively affect its productivity. This widely known medicinal plant with broad pharmacological properties (anti-depressant, anti-tumor, anti-inflammatory, antiviral, antioxidant, anti-cancer, and antibacterial) has been overlooked with respect to the identification of reference genes suitable for RT-qPCR data normalization. In this study, 11 candidate reference genes were analyzed in H. perforatum plants subjected to cold and heat stresses. The expression stability of these genes was assessed using GeNorm, NormFinder and BestKeeper algorithms. The results revealed that the ranking of stability among the three algorithms showed only minor differences within each treatment. The best-ranked reference genes differed between cold- and heat-treated samples; nevertheless, TUB was the most stable gene in both experimental conditions. GSA and GAPDH were found to be reliable reference genes in cold-treated samples, while GAPDH showed low expression stability in heat-treated samples. 26SrRNA and H2A had the highest stabilities in the heat assay, whereas H2A was less stable in the cold assay. Finally, AOX1, AOX2, CAT1 and CHS genes, associated with plant stress responses and oxidative stress, were used as target genes to validate the reliability of identified reference genes. These target genes showed differential expression profiles over time in treated samples. This study not only is the first systematic analysis for the selection of suitable reference genes for RT-qPCR studies in H. perforatum subjected to temperature stress conditions, but may also provide valuable information about the roles of genes associated with temperature stress responses.

Assessment of RNAi-induced silencing in banana (Musa spp.)

BACKGROUND

In plants, RNA- based gene silencing mediated by small RNAs functions at the transcriptional or post-transcriptional level to negatively regulate target genes, repetitive sequences, viral RNAs and/or transposon elements. Post-transcriptional gene silencing (PTGS) or the RNA interference (RNAi) approach has been achieved in a wide range of plant species for inhibiting the expression of target genes by generating double-stranded RNA (dsRNA). However, to our knowledge, successful RNAi-application to knock-down endogenous genes has not been reported in the important staple food crop banana.

RESULTS

Using embryogenic cell suspension (ECS) transformed with ß-glucuronidase (GUS) as a model system, we assessed silencing of gusAINT using three intron-spliced hairpin RNA (ihpRNA) constructs containing gusAINT sequences of 299-nt, 26-nt and 19-nt, respectively. Their silencing potential was analysed in 2 different experimental set-ups. In the first, Agrobacterium-mediated co-transformation of banana ECS with a gusAINT containing vector and an ihpRNA construct resulted in a significantly reduced GUS enzyme activity 6-8 days after co-cultivation with either the 299-nt and 19-nt ihpRNA vectors. In the second approach, these ihpRNA constructs were transferred to stable GUS-expressing ECS and their silencing potential was evaluated in the regenerated in vitro plants. In comparison to control plants, transgenic plants transformed with the 299-nt gusAINT targeting sequence showed a 4.5 fold down-regulated gusA mRNA expression level, while GUS enzyme activity was reduced by 9 fold. Histochemical staining of plant tissues confirmed these findings. Northern blotting used to detect the expression of siRNA in the 299-nt ihpRNA vector transgenic in vitro plants revealed a negative relationship between siRNA expression and GUS enzyme activity. In contrast, no reduction in GUS activity or GUS mRNA expression occurred in the regenerated lines transformed with either of the two gusAINT oligo target sequences (26-nt and 19-nt).

CONCLUSIONS

RNAi-induced silencing was achieved in banana, both at transient and stable level, resulting in significant reduction of gene expression and enzyme activity. The success of silencing was dependent on the targeted region of the target gene. The successful generation of transgenic ECS for second transformation with (an)other construct(s) can be of value for functional genomics research in banana.

Reference genes in real-time PCR

This paper aims to discuss various aspects of the use of reference genes in qPCR technique used in the thousands of present studies. Most frequently, these are housekeeping genes and they must meet several criteria so that they can lay claim to the name. Lots of papers report that in different conditions, for different organisms and even tissues the basic assumption—the constant level of the expression is not maintained for many genes that seem to be perfect candidates. Moreover, their transcription can not be affected by experimental factors. Sounds simple and clear but a great number of designed protocols and lack of consistency among them brings confusion on how to perform experiment properly. Since during selection of the most stable normalizing gene we can not use any reference gene, different ways and algorithms for their selection were developed. Such methods, including examples of best normalizing genes in some specific cases and possible mistakes are presented based on available sources. Numerous examples of reference genes applications, which are usually in too few numbers in relevant articles not allowing to make a solid fundament for a reader, will be shown along with instructive compilations to make an evidence for presented statements and an arrangement of future qPCR experiments. To include all the pitfalls and problems associated with the normalization methods there is no way not to begin from sample preparation and its storage going through candidate gene selection, primer design and statistical analysis. This is important because numerous short reviews available cover the topic only in lesser extent at the same time giving the reader false conviction of complete topic recognition.

Evaluation of reference genes for quantitative real-time PCR in oil palm elite planting materials propagated by tissue culture

BACKGROUND

The somatic embryogenesis tissue culture process has been utilized to propagate high yielding oil palm. Due to the low callogenesis and embryogenesis rates, molecular studies were initiated to identify genes regulating the process, and their expression levels are usually quantified using reverse transcription quantitative real-time PCR (RT-qPCR). With the recent release of oil palm genome sequences, it is crucial to establish a proper strategy for gene analysis using RT-qPCR. Selection of the most suitable reference genes should be performed for accurate quantification of gene expression levels.

RESULTS

In this study, eight candidate reference genes selected from cDNA microarray study and literature review were evaluated comprehensively across 26 tissue culture samples using RT-qPCR. These samples were collected from two tissue culture lines and media treatments, which consisted of leaf explants cultures, callus and embryoids from consecutive developmental stages. Three statistical algorithms (geNorm, NormFinder and BestKeeper) confirmed that the expression stability of novel reference genes (pOP-EA01332, PD00380 and PD00569) outperformed classical housekeeping genes (GAPDH, NAD5, TUBULIN, UBIQUITIN and ACTIN). PD00380 and PD00569 were identified as the most stably expressed genes in total samples, MA2 and MA8 tissue culture lines. Their applicability to validate the expression profiles of a putative ethylene-responsive transcription factor 3-like gene demonstrated the importance of using the geometric mean of two genes for normalization.

CONCLUSIONS

Systematic selection of the most stably expressed reference genes for RT-qPCR was established in oil palm tissue culture samples. PD00380 and PD00569 were selected for accurate and reliable normalization of gene expression data from RT-qPCR. These data will be valuable to the research associated with the tissue culture process. Also, the method described here will facilitate the selection of appropriate reference genes in other oil palm tissues and in the expression profiling of genes relating to yield, biotic and abiotic stresses.

Comprehensive selection of reference genes for gene expression normalization in sugarcane by real time quantitative rt-PCR

The increasingly used real time quantitative reverse transcription-PCR (qRT-PCR) method for gene expression analysis requires one or several reference gene(s) acting as normalization factor(s). In order to facilitate gene expression studies in sugarcane (Saccharum officinarum), a non-model plant with limited genome information, the stability of 13 candidate reference genes was evaluated. The geNorm, NormFinder and deltaCt methods were used for selecting stably expressed internal controls across different tissues and under various experimental treatments. These results revealed that, among these 13 candidate reference genes, GAPDH, eEF-1a and eIF-4α were the most stable and suitable for use as normalization factors across all various experimental samples. In addition, APRT could be a candidate for examining the relationship between gene copy number and transcript levels in sugarcane tissue samples. According to the results evaluated by geNorm, combining CUL and eEF-1α in hormone treatment experiments; CAC and CUL in abiotic stress tests; GAPDH, eEF-1a and CUL in all treatment samples plus CAC, CUL, APRT and TIPS-41 in cultivar tissues as groups for normalization would lead to more accurate and reliable expression quantification in sugarcane. This is the first systematic validation of reference genes for quantification of transcript expression profiles in sugarcane. This study should provide useful information for selecting reference genes for more accurate quantification of gene expression in sugarcane and other plant species.

Selection of suitable reference genes for qPCR normalization under abiotic stresses in Oenanthe javanica (BI.) DC

Accurate normalization of gene expression data is an absolute prerequisite to obtain reliable results in qPCR analysis. Oenanthe javanica, an aquatic perennial herb, belongs to the Oenanthe genus in Apiaceae family, with known medicinal properties. In the current study, O. javanica was subjected to hormone stimuli (gibberellin, salicylic acid, methyl jasmonate, and abscisic acid) and abiotic stresses (heat, cold, salt, and drought), and the expression of nine candidate reference genes (eIF-4α, ACT7, TIP41, GAPDH, SAND, EF-1α, PP2A, TBP, and TUB) was evaluated. Stability of the genes was assessed using geNorm, NormFinder and BestKeeper. All the genes presented distinct expression profiles under the experimental conditions analyzed. Under abiotic stress conditions, ACT7 and PP2A genes displayed the maximum stability; PP2A and SAND were the most stable genes under hormone stimuli. Even though PP2A gene was most stable across all the samples, individual analysis revealed changes in expression profile. To further validate the suitability of the reference genes identified in this study, the expression level of M6PR gene under salt treatment was studied. Based on our data, we propose that it is essential to normalize the target gene expression with specific reference genes under different experimental conditions for most accurate results. To our knowledge, this is the first systematic analysis for reference genes under abiotic stress and hormone stimuli conditions in O. javanica. This will be beneficial for future studies on O. javanica and other plants in Apiaceae family at molecular level.

Validation of reference genes aiming accurate normalization of qRT-PCR data in Dendrocalamus latiflorus Munro

BACKGROUND

Dendrocalamus latiflorus Munro distributes widely in subtropical areas and plays vital roles as valuable natural resources. The transcriptome sequencing for D. latiflorus Munro has been performed and numerous genes especially those predicted to be unique to D. latiflorus Munro were revealed. qRT-PCR has become a feasible approach to uncover gene expression profiling, and the accuracy and reliability of the results obtained depends upon the proper selection of stable reference genes for accurate normalization. Therefore, a set of suitable internal controls should be validated for D. latiflorus Munro.

RESULTS

In this report, twelve candidate reference genes were selected and the assessment of gene expression stability was performed in ten tissue samples and four leaf samples from seedlings and anther-regenerated plants of different ploidy. The PCR amplification efficiency was estimated, and the candidate genes were ranked according to their expression stability using three software packages: geNorm, NormFinder and Bestkeeper. GAPDH and EF1α were characterized to be the most stable genes among different tissues or in all the sample pools, while CYP showed low expression stability. RPL3 had the optimal performance among four leaf samples. The application of verified reference genes was illustrated by analyzing ferritin and laccase expression profiles among different experimental sets. The analysis revealed the biological variation in ferritin and laccase transcript expression among the tissues studied and the individual plants.

CONCLUSIONS

geNorm, NormFinder, and BestKeeper analyses recommended different suitable reference gene(s) for normalization according to the experimental sets. GAPDH and EF1α had the highest expression stability across different tissues and RPL3 for the other sample set. This study emphasizes the importance of validating superior reference genes for qRT-PCR analysis to accurately normalize gene expression of D. latiflorus Munro.

Reference genes selection for quantitative real-time PCR using RankAggreg method in different tissues of Capra hircus

Identification of reference genes with stable levels of gene expression is an important prerequisite for obtaining reliable results in analysis of gene expression data using quantitative real time PCR (RT-qPCR). Since the underlying assumption of reference genes is that expressed at the exact same level in all sample types, in this study, we evaluated the expression stability of nine most commonly used endogenous controls (GAPDH, ACTB, 18S rRNA, RPS18, HSP-90, ALAS, HMBS, ACAC, and B2M) in four different tissues of the domestic goat, Capra hircus, including liver, visceral, subcutaneous fat and longissimus muscles, across different experimental treatments (a standard diet prepared using the NRC computer software as control and the same diet plus one mg chromium/day). We used six different software programs for ranking of reference genes and found that individual rankings of the genes differed among them. Additionally, there was a significant difference in ranking patterns of the studied genes among different tissues. A rank aggregation method was applied to combine the ranking lists of the six programs to a consensus ranking. Our results revealed that HSP-90 was nearly always among the two most stable genes in all studied tissues. Therefore, it is recommended for accurate normalization of RT-qPCR data in goats, while GAPDH, ACTB, and RPS18 showed the most varied expressions and should be avoided as reference genes.

Transgenic banana plants overexpressing a native plasma membrane aquaporin MusaPIP1;2 display high tolerance levels to different abiotic stresses

Water transport across cellular membranes is regulated by a family of water channel proteins known as aquaporins (AQPs). As most abiotic stresses like suboptimal temperatures, drought or salinity result in cellular dehydration, it is imperative to study the cause-effect relationship between AQPs and the cellular consequences of abiotic stress stimuli. Although plant cells have a high isoform diversity of AQPs, the individual and integrated roles of individual AQPs in optimal and suboptimal physiological conditions remain unclear. Herein, we have identified a plasma membrane intrinsic protein gene (MusaPIP1;2) from banana and characterized it by overexpression in transgenic banana plants. Cellular localization assay performed using MusaPIP1;2::GFP fusion protein indicated that MusaPIP1;2 translocated to plasma membrane in transformed banana cells. Transgenic banana plants overexpressing MusaPIP1;2 constitutively displayed better abiotic stress survival characteristics. The transgenic lines had lower malondialdehyde levels, elevated proline and relative water content and higher photosynthetic efficiency as compared to equivalent controls under different abiotic stress conditions. Greenhouse-maintained hardened transgenic plants showed faster recovery towards normal growth and development after cessation of abiotic stress stimuli, thereby underlining the importance of these plants in actual environmental conditions wherein the stress stimuli is often transient but severe. Further, transgenic plants where the overexpression of MusaPIP1;2 was made conditional by tagging it with a stress-inducible native dehydrin promoter also showed similar stress tolerance characteristics in in vitro and in vivo assays. Plants developed in this study could potentially enable banana cultivation in areas where adverse environmental conditions hitherto preclude commercial banana cultivation.

Proteomic analysis of Fusarium oxysporum f. sp. cubense tropical race 4-inoculated response to Fusarium wilts in the banana root cells

BACKGROUND

Fusarium wilt of banana is one of the most destructive diseases in the world. This disease has caused heavy losses in major banana production areas. Except for molecular breeding methods based on plant defense mechanisms, effective methods to control the disease are still lacking. Dynamic changes in defense mechanisms between susceptible, moderately resistant, and highly resistant banana and Fusarium oxysporum f. sp. cubense tropical race 4 (Foc4) at the protein level remain unknown. This research reports the proteomic profile of three banana cultivars in response to Foc4 and transcriptional levels correlated with their sequences for the design of disease control strategies by molecular breeding.

RESULTS

Thirty-eight differentially expressed proteins were identified to function in cell metabolism. Most of these proteins were positively regulated after Foc4 inoculation. These differentially regulated proteins were found to have important functions in banana defense response. Functional categories implicated that these proteins were associated with pathogenesis-related (PR) response; isoflavonoid, flavonoid, and anthocyanin syntheses; cell wall strengthening; cell polarization; reactive oxygen species production and scavenging; jasmonic acid-, abscisic acid-, and auxin-mediated signaling conduction; molecular chaperones; energy; and primary metabolism. By comparing the protein profiles of resistant and susceptible banana cultivars, many proteins showed obvious distinction in their defense mechanism functions. PR proteins in susceptible 'Brazil' were mainly involved in defense. The proteins related to PR response, cell wall strengthening and antifungal compound synthesis in moderately resistant 'Nongke No.1' were mainly involved in defense. The proteins related to PR response, cell wall strengthening, and antifungal compound synthesis in highly resistant 'Yueyoukang I' were mainly involved in defense. 12 differentially regulated genes were selected to validate through quantitative real time PCR method. Quantitative RT-PCR analyses of these selected genes corroborate with their respective protein abundance after pathogen infection.

CONCLUSIONS

This report is the first to use proteomic profiling to study the molecular mechanism of banana roots infected with Foc4. The differentially regulated proteins involved in different defense pathways are likely associated with different resistant levels of the three banana cultivars.

Reliable reference genes for normalization of gene expression in cucumber grown under different nitrogen nutrition

In plants, nitrogen is the most important nutritional factor limiting the yield of cultivated crops. Since nitrogen is essential for synthesis of nucleotides, amino acids and proteins, studies on gene expression in plants cultivated under different nitrogen availability require particularly careful selection of suitable reference genes which are not affected by nitrogen limitation. Therefore, the objective of this study was to select the most reliable reference genes for qPCR analysis of target cucumber genes under varying nitrogen source and availability. Among twelve candidate cucumber genes used in this study, five are highly homologous to the commonly used internal controls, whereas seven novel candidates were previously identified through the query of the cucumber genome. The expression of putative reference genes and the target CsNRT1.1 gene was analyzed in roots, stems and leaves of cucumbers grown under nitrogen deprivation, varying nitrate availability or different sources of nitrogen (glutamate, glutamine or NH3). The stability of candidate genes expression significantly varied depending on the tissue type and nitrogen supply. However, in most of the outputs genes encoding CACS, TIP41, F-box protein and EFα proved to be the most suitable for normalization of CsNRT1.1 expression. In addition, our results suggest the inclusion of 3 or 4 references to obtain highly reliable results of target genes expression in all cucumber organs under nitrogen-related stress.

Identification and validation of reference genes for transcript normalization in strawberry (Fragaria × ananassa) defense responses

Strawberry (Fragaria spp) is an emerging model for the development of basic genomics and recombinant DNA studies among rosaceous crops. Functional genomic and molecular studies involve relative quantification of gene expression under experimental conditions of interest. Accuracy and reliability are dependent upon the choice of an optimal reference control transcript. There is no information available on validated endogenous reference genes for use in studies testing strawberry-pathogen interactions. Thirteen potential pre-selected strawberry reference genes were tested against different tissues, strawberry cultivars, biotic stresses, ripening and senescent conditions, and SA/JA treatments. Evaluation of reference candidate's suitability was analyzed by five different methodologies, and information was merged to identify best reference transcripts. A combination of all five methods was used for selective classification of reference genes. The resulting superior reference genes, FaRIB413, FaACTIN, FaEF1α and FaGAPDH2 are strongly recommended as control genes for relative quantification of gene expression in strawberry. This report constitutes the first systematic study to identify and validate optimal reference genes for accurate normalization of gene expression in strawberry plant defense response studies.

MusaSAP1, a A20/AN1 zinc finger gene from banana functions as a positive regulator in different stress responses

A20/AN1 zinc finger domain containing Stress Associated Proteins (SAP) are involved in diverse stress response pathways in plants. In the present study, a novel banana SAP gene, MusaSAP1, was identified from banana EST database and was subsequently characterized by overexpression in transgenic banana plants. Expression profiling in native banana plants showed that MusaSAP1 was up-regulated by drought, salt, cold, heat and oxidative stress as well as by treatment with abscisic acid. Cellular localization assay carried out by making a MusaSAP1::GFP fusion protein indicated that MusaSAP1 is incompletely translocated to nucleus. Copy number analysis performed using real time PCR and Southern blotting indicated that MusaSAP1 occurs in the banana genome in a single copy per 11 chromosome set. Transgenic banana plants constitutively overexpressing MusaSAP1 displayed better stress endurance characteristics as compared to controls in both in vitro and ex vivo assays. Lesser membrane damage as indicated by reduced malondialdehyde levels in transgenic leaves subjected to drought, salt or oxidative stress pointed towards significant role for MusaSAP1 in stress amelioration pathways of banana. Strong up-regulation of a polyphenol oxidase (PPO) coding transcript in MusaSAP1 overexpressing plants together with induction of MusaSAP1 by wounding and methyl jasmonate treatment indicated possible involvement of MusaSAP1 in biotic stress responses where PPOs perform major functions in multiple defense pathways.