Skip to main content
  • Original article
  • Open access
  • Published:

Selection of reference genes from Shiraia bambusicola for RT-qPCR analysis under different culturing conditions

Abstract

Stable reference genes are necessary to analyse quantitative real-time reverse transcription PCR (qRT-PCR) data and determine the reliability of the final results. For further studies of the valuable fungus Shiraia bambusicola, the identification of suitable reference genes has become increasingly urgent. In this study, three conventional reference genes and nine novel candidates were evaluated under different light conditions (all-dark, all-light and 12-h light/dark) and in different media (rice medium, PD medium, and Czapek–Dox medium). Three popular software programs (geNorm, NormFinder and BestKeeper) were used to analyse these genes, and the final ranking was determined using RefFinder. SbLAlv9, SbJsn1, SbSAS1 and SbVAC55 displayed the best stability among the genes, while SbFYVE and SbPKI showed the worst. These emerging genes exhibited significantly better properties than the three existing genes under almost all conditions. Furthermore, the most reliable reference genes were identified separately under different nutrient and light conditions, which would help accessible to make the most of the existing data. In summary, a group of novel housekeeping genes from S. bambusicola with more stable properties than before was explored, and these results could also provide a practical approach for other filamentous fungi.

Introduction

Shiraia bambusicola is an important and valuable macrofungus in the medical and food industries. It is noteworthy for its hypocrellins, the main secondary metabolites of S. bambusicola, whose use has been proposed for disease treatments that involve anti-clinical strains and anti-inflammatory and anti-viral activity (Su et al. 2011; Jiang et al. 2011; Zhou et al. 2009). Currently, to break the bottleneck of product yield and improve the understanding of biosynthetic pathways, the analysis of functional genes is attracting increasing attention (Deng et al. 2016).

Quantitative real-time reverse transcription PCR (qRT-PCR) is generally regarded as a convenient and efficient tool to analyse gene expression, but the RNA quantification, the reverse transcription reaction efficiency and other uncontrolled factors may limit the accuracy and stability of the final results (Huang et al. 2014; Hao et al. 2014). Thus, it is necessary to apply reliable reference genes to normalize the data.

A large number of research papers have been published on reference genes under different stresses or from different organs in plants (Warzybok and Migocka 2013; Lin et al. 2014), and similar works in filamentous fungi are gradually being conducted (Zampieri et al. 2014; Zhou et al. 2011). It is striking that most of the applied reference genes for fungi were directly copied from the existing results in plants or animals, such as actin, tubulin and 18S rRNA (Fang and Bidochka 2006; Yan and Liou 2006). However, the divergent regulatory mechanisms and environmental stresses suggested that these reliable reference genes might not be suitable for the analysis of fungal gene expression. For example, the classical reference genes elongation factor (EF-1) and beta-tubulin (β-tubulin), which are widely used in plants, were not appropriate in Hemileia vastatrix (Vieira et al. 2011); another reliable reference gene, glyceraldehyde-3-phosphate dehydrogenase (GAPDH), used in several qPCR studies for the normalization of data, was not stable in Heterobasidion annosum (Raffaello and Asiegbu 2013). Thus, it is necessary to systemically rescreen the candidates in novel species from the fungal kingdom.

According to the current literature, exploring putative reference genes from the transcriptomic data is a feasible and easy assay in fungal species (Llanos et al. 2015; Nadai et al. 2015; Steiger et al. 2010). Likewise, there is a series of popular Excel-based software programs, such as geNorm, NormFinder and BestKeeper, that have been successfully used in previous studies (Vandesompele et al. 2009; Andersen et al. 2004; Pfaffl et al. 2004).

Thus far, research on S. bambusicola has focused on strain mutagenesis and product fermentation (Song et al. 2015). However, since the genome and transcriptome were published (Yang et al. 2014; Zhao et al. 2016), gene function and molecular regulation have become research hotspots. To make qRT-PCR results more reliable and the application of reference genes more widespread, we comprehensively rescreened the candidate reference genes for S. bambusicola under different light conditions and in different media.

Materials and Methods

Fungal strain and culture conditions

Strain zzz816 of Shiraia bambusicola was deposited in the China General Microbiological Culture Collection Centre (CGMCC, No: 3135). Fungal isolates adhered to agars were cultured in potato dextrose agar (PDA) medium (potato 200 g/L, dextrose 20 g/L, agar 15 g/L) at 26 °C for 7 days, and then the fresh mycelia were inoculated onto different media (rice medium: rice 800 g/L, PD medium: potato 200 g/L, dextrose 20 g/L, and Czapek–Dox medium: NaNO3 3 g/L, K2HPO4 1 g/L, MgSO47H2O 0.5 g/L, KCl 0.5 g/L, FeSO4 0.01 g/L, sucrose 30 g/L) under three different light conditions, namely, 24 h of continuous darkness (all-dark), 24 h of continuous lighting with a light intensity of 1500 lx (all-light), and 12: a 12 h light photoperiod with 1500 lx light intensity (12-h light/dark).

Finally, each sample was collected at different time points: 2, 4 and 6 days.

RNA isolation and cDNA synthesis

The total RNA of each sample was extracted from frozen mycelia using an improved method in our laboratory (Song et al. 2015), and the total RNA was dissolved in RNase-free water.

The first strand of cDNA was synthesized by reverse transcribing 1 μg of RNA with TransScript ® All-in-One First-Strand cDNA Synthesis SuperMix for qPCR (One-Step gDNA Removal; TransGen, China). The quantity and quality of the total RNA extracted was determined using an Eppendorf BioPhotometer plus (Eppendorf, Ger), and the cDNA samples were stored at −20 °C.

Primer design and quantitative real-time PCR

Primers were designed using Primer3 software (http://primer3.ut.ee/), and the specificity of the product was verified using 1% agarose gel electrophoresis and melting curves. The efficiency of the validated primer pairs remained at approximately 100% (Table 1). The qRT-PCR reaction was performed in LightCycler®480 multiwell plates (Roche Applied Science, Indianapolis, IN, USA) with a LightCycler®480II/96 (Roche Applied Science, Indianapolis, IN, USA) real-time PCR system using LightCycler®480 SYBR Green I Master Mix (Roche Applied Science, Indianapolis, IN, USA). The reactions were performed according to the recommendations of the manufacturer: 95 °C for 5 min for initial denaturation, followed by 45 thermal cycles of 10 s at 95 °C, 10 s at 60 °C and 20 s at 72 °C. The melting curve was performed with slow heating from 65 to 97 °C with continuous measurement of fluorescence 5 acquisitions per 1 °C. All reactions were performed with three biological and two technical replicates with negative controls. The qRT-PCR data were directly analysed using the “second derivative maximum” function in the LightCycler®480 Software Version 1.5 (Roche Applied Science, Indianapolis, IN, USA).

Table 1 Description of candidate reference genes, and the details of primers and amplicons

Statistical analyses

To select suitable reference genes, three software packages were used to calculate the stability: NormFinder, geNorm, and BestKeeper. An additional web-based tool, RefFinder, (http://fulxie.0fees.us/) was applied to integrate and rank all candidate reference genes (Xie et al. 2012).

To calculate the PCR amplification efficiencies (E) and correlation coefficients (R2) of each primer pair, standard curves were prepared using a 10-fold serial diluting plasmid, into which the reference gene was cloned in PEASY-T3® (TransGen, China).

The efficiency (E) was calculated according to the equation E = (10(−1/slope)−1) × 100% (Radonić et al. 2004).

Results

Strategy for selecting reference gene candidates

In this study, twelve candidate reference genes appeared in the stabilization assay, and nine of them were involved in this test for the first time. The novel ones were sought out directly from the public transcriptome data sets (Zhao et al. 2016), and the selection was based on the ranking of the expression levels of each gene, expressed as reads per kilobase per million (RKPM). For the whole data set, we evaluated the coefficient of variation of RKPM, and the ones with relatively smaller coefficients were considered the genes of interest, including GTP-binding protein SAS1 (SbSAS1), 4A/4B type thioredoxin-like protein (SbTRX), t-SNARE (SbtS), RNA binding protein-like protein Jsn1 (SbJsn1), conserved hypothetical protein (SbCHP), LAlv9 family protein (SbLAlv9), Pkinase-domain-containing protein (SbPKI), FYVE-domain-containing protein (SbFYVE) and vacuolar protein sorting 55 (SbVAC55). Other three genes, ubiquitin-activating enzyme (UBI), vacuolar ATPase subunit 1 (VAC) and transcription factor TFIIIC (TFC), were generally used to normalize the qRT-PCR data and appeared to be especially reliable as reference genes in the preliminary study for Shiraia bambusicola (Song et al. 2015).

Before the analysis of expression stability, the gene-specific amplification of these genes was confirmed by the single-peak melting curves of the qRT-PCR products, and the primers provided a good reaction efficiency ranging between 96 and 102% (Table 1).

The crossing point (CP) values were collected from all tested samples under different conditions and are shown in the box-plot (Fig. 1). The value of gene expression studied indicated a compact distribution and a limited range between 21 and 33. Among the genes, TFC displayed a lower expression variation, which mainly depended on the different media (Song et al. 2015).

Fig. 1
figure 1

The range of CP values of 12 reference gene candidates for all samples. Each box indicates the 25th and 75th percentiles, and the caps represent the maximum and minimum values. The median is shown by the line across the box

Expression stability analysis

Three different applets, geNorm, NormFinder, and BestKeeper, were applied to measure and rank the stability of candidate reference genes. The program geNorm (Steiger et al. 2010) classifies genes according to the control gene stability measure (M value), which represents the average of the pair-wise variation of a gene with all other control genes. NormFinder (Vandesompele et al. 2009) examines the stability of each single candidate gene independently and not in relation to the other genes. The results from geNorm and NormFinder can be compared easily because they both use raw data (relative quantities) as input data. BestKeeper (Andersen et al. 2004) is another Excel-based tool that determines the optimal reference genes using a pair-wise correlation analysis (Pearson correlation coefficient) of all pairs of candidate genes. It uses CP values (instead of relative quantities) as input and employs a different measure of expression stability from geNorm and NormFinder. However, this software is not able to analyse more than ten reference genes together, and thus the first ten ranking genes should depend on geNorm and NormFinder. The rankings of these software programs relied on different algorithms and were expected to lead to distinct outputs. Therefore, another software program, RefFinder (Radonić et al. 2004), was used to integrate the currently available major computational programs (geNorm, NormFinder, BestKeeper, and the comparative ΔΔCt method) to compare and rank the tested candidate reference genes. Based on the rankings from each program, it could assign an appropriate weight to an individual gene and calculate the geometric mean of their weights for the overall final ranking.

Using geNorm, SbLAlv9, SbJsn1, SbCHP and SbSAS1 displayed the highest reliability overall (Table 2). In the rice medium with difference light conditions, SbLAlv9, SbSAS1, SbJsn1 and TFC had a better effect than other candidate reference genes. In the PD medium, SbLAlv9, SbJsn1, SbtS, and SbSAS1 ranked at the top positions, and in the Czapek-Dox medium, SbJsn1, VAC, SbLAlv9 and SbFYVE were found to be the ideal reference genes. In the different light conditions, the best reference candidates were divided into three groups, including SbLAlv9, SbJsn1, VAC and SbCHP under all-dark conditions, SbLAlv9, SbJsn1, SbCHP, SbtS under all-light conditions and SbLAlv9, SbJsn1, SbSAS1, SbtS under the 12-h light/dark conditions. SbFYVE and SbPKI have been classified as the least reliable reference genes in most of the conditions. Furthermore, geNorm provided an output allowing a set of reliable normalization for the pairwise variation (Vn/n+1) to help to determine the optimal number of reference genes (Fig. 2). Two reference genes were sufficient for most of the conditions, but the continuous darkness required a third reference gene. Four reference genes were suggested to ensure accurate all-round analysis of the nutrient and light conditions.

Table 2 Ranking of candidate reference genes calculated by geNorm according to different expression conditions
Fig. 2
figure 2

The geNorm-based results of the pairwise variation analysis. A sufficient number of genes n can be used for reliable normalization when Vn/n + 1<0.15

As shown in Table 3, SbLAlv9, SbJsn1 and SbSAS1 were identified as the best reference candidates by NormFinder, and SbCHP, SbtS and VAC displayed good properties for certain nutrient and light conditions. SbFYVE and SbPKI were ranked as the most unstable genes, and this result was also in agreement with the geNorm calculation.

Table 3 Ranking of candidate reference genes calculated by NormFinder according to different expression conditions

In contrast to geNorm and NormFinder, the BestKeeper algorithm is based on the coefficient of variance (CV) and the standard deviation (SD) calculated by the average CP value of each reaction. Moreover, we removed SbFYVE and SbPKI due to the higher variations based on geNorm and NormFinder. (Table 4). SbVAC55, SbSAS1, UBI and SbTRX were the four best reference genes among all the samples. Based on the different points, SbSAS1, TFC, UBI and SbTRX were the most stable on the different media, and SbVAC55, UBI, SbSAS1 and SbtS were identified as suitable reference genes for the different light conditions, including all-dark, all-light and 12-h light/dark.

Table 4 Ranking of candidate reference genes were calculated by BestKeeper according to different expression conditions

To integrate the results obtained from different algorithms, RefFinder was used to perform the final ranking. As illustrated in Table 5, SbLAlv9 was singled out as the top candidate in the global consideration of different nutritional and light conditions, followed by SbJsn1, SbSAS1 and SbVAC55. Nevertheless, under specific conditions, the gene order would be altered, and the reliability of proper reference genes might be reduced. Thus, we propose two different methods—one medium with different light conditions and one light condition with different media—to improve the utilization of reference genes in practical experiments.

Table 5 Ranking of the candidate reference genes calculated using RefFinder under different conditions

Application in detailed conditions

Although the appropriate reference genes had been presented as in Table 5, it was still inconvenient to search for a target effectively for one particular condition. To optimize the application of the ideal reference genes, we reprocessed the data and characterized the reference genes as a two-way table (Table 6). Therein, we could find the suitable reference genes more easily and accurately. For example, under continuous darkness (Table 5) and limited to rice medium, SbVAC55 or SbLAlv9 was sufficient for qRT-PCR analysis (Table 6).

Table 6 The most stable reference genes under different nutritional and light conditions

Discussion

Of all the candidate reference genes, three familiar genes (UBI, VAC and TFC) had been approved previously (Song et al. 2015), and another nine (SbSAS1, SbTRX, SbtS, SbJsn1, SbCHP, SbLAlv9, SbPKI, SbFYVE, SbVAC55) first emerged for the reliability of gene expression. Three new genes, SbLAlv9, SbJsn1 and SbSAS1, indicated appropriate properties for all the nutrient and light conditions. SbtS and SbVAC55 also separately showed high reliability for certain conditions. In contrast, only one old reference gene, VAC, displayed good activity for Czapek–Dox medium or all-dark conditions. As illustrated by these results, the exploration of new reference genes could indeed provide the more reliable genes and improve the assurance of the stability of gene expression in S. bambusicola. Furthermore, the discovery of novel reference genes might contribute to the use of reference genes in other species.

We chose typical nutrient and light conditions in this study that generally cover the culturing conditions of S. bambusicola or other species. Therefore, the reliable reference genes in this work were anticipated to be broadly used in the gene expression analysis of S. bambusicola and could also provide benefits for other fungal species. Furthermore, a table (Table 6) was constructed to facilitate the application of these results, with a two-way form to facilitate finding the stable reference genes for specific conditions. In this work, we tried our best to expand the scope of application and make the process more convenient.

In conclusion, nine novel candidate reference genes were introduced for the first time, and several of them (SbLAlv9, SbJsn1, SbSAS1 and SbVAC55) were validated as ideal ones. In addition, this work provided a set of reference genes under different nutrient and light conditions, and these housekeeping genes are expected to be available in a more extensive range.

Abbreviations

qRT-PCR:

quantitative real-time reverse transcription PCR

EF-1:

elongation factor

β-tubulin:

beta-tubulin

GAPDH:

glyceraldehyde-3-phosphate dehydrogenase

RKPM:

reads per kilobase per million

SbSAS1:

GTP-binding protein SAS1

SbTRX:

4A/4B type thioredoxin-like protein

SbtS:

t-SNARE

SbJsn1:

RNA binding protein-like protein Jsn1

SbCHP:

conserved hypothetical protein

SbLAlv9:

LAlv9 family protein

SbPKI:

Pkinase-domain-containing protein

SbFYVE:

FYVE-domain-containing protein

SbVAC55:

vacuolar protein sorting 55

UBI:

ubiquitin-activating enzyme

VAC:

vacuolar ATPase subunit 1

TFC:

transcription factor TFIIIC

CP:

crossing point

CV:

coefficient of variance

SD:

standard deviation

References

  • Andersen CL, Jensen JL, Ørntoft TF (2004) Normalization of realtime quantitative reverse transcription-PCR data: a model-based variance estimation approach to identify genes suited for normalization, applied to bladder and colon cancer data sets. Cancer Res 64(15):5245–5250

    Article  CAS  PubMed  Google Scholar 

  • Deng H, Gao R, Liao X, Cai Y (2016) Reference genes selection and relative expression analysis from Shiraia sp. SUPER-H168 productive of hypocrellin. Gene 580(1):67–72

    Article  CAS  PubMed  Google Scholar 

  • Fang W, Bidochka M (2006) Expression of genes involved in germination, conidiogenesis and pathogenesis in Metarhizium anisopliae using quantitative real-time RT-PCR. Mycol Res 110(10):1165–1171

    Article  CAS  PubMed  Google Scholar 

  • Hao X, Horvath D, Chao W, Yang Y, Wang X, Xiao B (2014) Identification and evaluation of reliable reference genes for quantitative real-time PCR analysis in Tea Plant (Camellia sinensis (L.) O. Kuntze). Int J Mol Sci 5(12):22155–22172

    Article  Google Scholar 

  • Huang L, Yan H, Jiang X, Zhang X, Zhang Y, Huang X, Zhang Y, Miao J, Xu B, Frazier T, Zhao B (2014) Evaluation of candidate reference genes for normalization of quantitative RT-PCR in Switchgrass under various abiotic stress conditions. BioEnergy Res 7(4):1201–1211

    Article  CAS  Google Scholar 

  • Jiang Y, Leung A, Xiang J, Xu C (2011) Blue light-activated hypocrellin B damages ovarian cancer cells. Laser Phys 22(1):300–305

    Article  Google Scholar 

  • Lin Y, Zhang C, Lan H, Gao S, Liu H, Liu J, Cao M, Pan G, Rong T, Zhang S (2014) Validation of potential reference genes for qPCR in Maize across abiotic stresses, hormone treatments, and tissue types. PLoS ONE 9(5):e95445

    Article  PubMed  PubMed Central  Google Scholar 

  • Llanos A, François J, Parrou J (2015) Tracking the best reference genes for RT-qPCR data normalization in filamentous fungi. BMC Genom 16(1):71

    Article  Google Scholar 

  • Nadai C, Campanaro S, Giacomini A, Corich V (2015) Selection and validation of reference genes for quantitative real-time PCR studies during Saccharomyces cerevisiae alcoholic fermentation in the presence of sulfite. Int J Food Microbiol 215:49–56

    Article  CAS  PubMed  Google Scholar 

  • Pfaffl M, Tichopad A, Prgomet C, Neuvians T (2004) Determination of stable housekeeping genes, differentially regulated target genes and sample integrity: BestKeeper—Excel-based tool using pair-wise correlations. Biotechnol Lett 26(6):509–515

    Article  CAS  PubMed  Google Scholar 

  • Radonić A, Thulke S, Mackay I, Landt O, Siegert W, Nitsche A (2004) Guideline to reference gene selection for quantitative real-time PCR. Biochem Biophys Res Commun 313(4):856–862

    Article  PubMed  Google Scholar 

  • Raffaello T, Asiegbu F (2013) Evaluation of potential reference genes for use in gene expression studies in the conifer pathogen (Heterobasidion annosum). Mol Biol Rep 40(7):4605–4611

    Article  CAS  PubMed  Google Scholar 

  • Song L, Li T, Fan L, Shen X, Hou C (2015) Identification and evaluation of reliable reference genes in the medicinal fungus Shiraia bambusicola. Curr Microbiol 72(4):444–449

    Article  PubMed  Google Scholar 

  • Steiger M, Mach R, Mach-Aigner A (2010) An accurate normalization strategy for RT-qPCR in Hypocrea jecorina (Trichoderma reesei). J Biotechnol 145(1):30–37

    Article  CAS  PubMed  Google Scholar 

  • Su Y, Sun J, Rao S, Cai Y, Yang Y (2011) Photodynamic antimicrobial activity of hypocrellin A. J Photoch Photobio B 103(1):29–34

    Article  CAS  Google Scholar 

  • Vandesompele J, Kubista M, Pfaffl MW (2009) Reference gene validation software for improved normalization. In: Logan J (ed) Real-time PCR: current technology and applications, 1st edn. Caister Academic, London, pp 47–64

    Google Scholar 

  • Vieira A, Talhinhas P, Loureiro A, Duplessis S, Fernandez D, Silva M, Paulo O, Azinheira H (2011) Validation of RT-qPCR reference genes for in planta expression studies in Hemileia vastatrix, the causal agent of coffee leaf rust. Fungal Biol 115(9):891–901

    Article  CAS  PubMed  Google Scholar 

  • Warzybok A, Migocka M (2013) Reliable reference genes for normalization of gene expression in cucumber grown under different nitrogen nutrition. PLoS ONE 8(9):e72887

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  • Xie F, Xiao P, Chen D, Xu L, Zhang B (2012) miRDeepFinder: a miRNA analysis tool for deep sequencing of plant small RNAs. Plant Mol Biol 80(1):75–84

    Article  CAS  Google Scholar 

  • Yan H, Liou R (2006) Selection of internal control genes for real-time quantitative RT-PCR assays in the oomycete plant pathogen Phytophthora parasitica. Fungal Genet Biol 43(6):430–438

    Article  CAS  PubMed  Google Scholar 

  • Yang H, Wang Y, Zhang Z, Yan R, Zhu D (2014) Whole-genome shotgun assembly and analysis of the genome of Shiraia sp. Strain Slf14, a novel endophytic fungus producing Huperzine A and Hypocrellin A. Genome Announc 2(1):11–14

    Article  Google Scholar 

  • Zampieri D, Nora L, Basso V, Camassola M, Dillon A (2014) Validation of reference genes in Penicillium echinulatum to enable gene expression study using real-time quantitative RT-PCR. Curr Genet 60(3):231–236

    Article  CAS  PubMed  Google Scholar 

  • Zhao N, Lin X, Qi S, Luo Z, Chen S, Yan S (2016) De Novo Transcriptome assembly in Shiraia bambusicola to investigate putative genes involved in the biosynthesis of Hypocrellin A. Int J Mol Sci 17(3):311

    Article  PubMed  PubMed Central  Google Scholar 

  • Zhou J, Wu X, Gu X, Zhou L, Song K, Wei S, Feng Y, Shen J (2009) Spectroscopic studies on the interaction of hypocrellin A and hemoglobin. Spectrochim Acta A 72(1):151–155

    Article  Google Scholar 

  • Zhou Y, Zhang Y, Luo Z, Fan Y, Tang G, Liu L, Pei Y (2011) Selection of optimal reference genes for expression analysis in the entomopathogenic fungus Beauveria bassiana during development, under changing nutrient conditions, and after exposure to abiotic stresses. Appl Microbiol Biotechnol 93(2):679–685

    Article  PubMed  Google Scholar 

Download references

Authors’ contributions

CZ, CLH and XYS designed research. CZ, TL and XYS preformed all experiments. CLH and XYS contributed data analysis. All authors were involved in data interpretation and the writing of the paper. All authors read and approved the final manuscript.

Competing interests

The authors declare that they have no competing interests.

Availability of data and materials

The datasets supporting the conclusions of this article are included within the article and shown in tables.

Ethics approval

This article does not contain any studies concerned with experiment on human or animals.

Funding

This study was funded by National Natural Science Foundation of China (Grant Numbers 31500015 and 31470145), Scientific Research Common Program of Beijing Municipal Commission of Education (Grant Number KM201510028009).

Author information

Authors and Affiliations

Authors

Corresponding authors

Correspondence to Cheng-Lin Hou or Xiao-Ye Shen.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhang, C., Li, T., Hou, CL. et al. Selection of reference genes from Shiraia bambusicola for RT-qPCR analysis under different culturing conditions. AMB Expr 7, 14 (2017). https://doi.org/10.1186/s13568-016-0314-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13568-016-0314-9

Keywords