- Research note
- Open Access
Reference genes for gene expression studies targeting sugarcane infected with Sugarcane mosaic virus (SCMV)
BMC Research Notes volume 12, Article number: 149 (2019)
The selection of reference genes in sugarcane under Sugarcane mosaic virus (SCMV) infection has not been reported and is indispensable to get reliable reverse transcription quantitative PCR (RT-qPCR) results for validation of transcriptome analysis. In this regard, seven potential reference genes were tested by RT-qPCR and ranked according to their stability using BestKeeper, NormFinder and GeNorm algorithms, and RefFinder WEB-based software in an experiment performed with samples from two sugarcane cultivars contrasting for SCMV resistance, when mechanically inoculated with a severe SCMV strain and using mock inoculated plant controls.
The genes Uridylate kinase (UK) and Ubiquitin-conjugating enzyme 18 (UBC18) were the most stable according to GeNorm algorithm and the Pearson correlation coefficients with the BestKeeper index. On the other hand, ribosomal protein L35-4 (RPL1), Actin (ACT) and Ubiquitin1 (UBQ1) were the least stable genes for all algorithms tested.
Sugarcane mosaic disease (SMD) is widely distributed among sugarcane-growing countries and may be caused by different virus species of the genera Potyvirus and Poacevirus, family Potyviridae . In Brazil, Sugarcane mosaic virus (SCMV), Potyvirus, is one the main viruses affecting sugarcane and the only causal agent of SMD, to date [2, 3]. The disease is controlled by the use of resistant cultivars making the comprehension of molecular bases of resistance to these viruses of great concern for sugarcane breeding programs worldwide [1, 2]. Transcriptome analysis has been applied in sugarcane to identify differentially expressed genes associated with biological traits [4,5,6] yet, few studies have investigated changes in the sugarcane transcriptome under infection by mosaic-causing viruses [7, 8]. The validation of transcriptome results via reverse transcription quantitative PCR (RT-qPCR) requires a normalization step for reducing its uncertainties [9, 10], commonly attained by the use of endogenous reference genes [11, 12]. The choice of appropriate reference genes is an essential step, since improper selection of references genes may result in unreliable RT-qPCR results . Several algorithms are available for identification of reliable candidate reference genes , which is a necessary procedure whenever different experimental conditions and genotypes are involved [13, 14], however such studies for SMD are yet to be reported. In this regard, the present study aimed to evaluate seven candidate reference genes based on previous reports in sugarcane under drought stress  and in closely related monocot species under viral infection .
Plant material and experimental design
The biological samples used in this study proceed from a previous experiment performed by Medeiros et al. . It was used a completely randomized factorial design with three factors, under greenhouse conditions: (a) two sugarcane cultivars, IACSP95-5000, resistant to SCMV, and IAC91-1099, susceptible to SCMV, both from the “Sugarcane Breeding Program, Instituto Agronômico de Campinas, IAC”, Brazil (b) two treatments (SCMV inoculated, s.i; and mock inoculated, m.i), and (c) a time course experiment with three sampling time points of the + 1 leaf, 24, 48 and 72 h post inoculation (hpi). Briefly, 36 sugarcane plantlets of each cultivar were obtained by meristem tip culture and indexed as virus-free by reverse transcription PCR (RT-PCR) using specific primers for the SCMV capsid protein . At 1-month-old, 18 plantlets of each cultivar were submitted to the s.i treatment using a severe strain of SCMV (SCMV Rib-1)  and the remaining 18 were submitted to the m.i treatment, according to Bain method . Therefore, six biological replicates were used for each combination of experimental factors. Among the s.i and m.i samples indexed by RT-PCR respectively as virus infected and virus free, three biological replicates from 24 and 72 hpi for each genotype × treatment were selected for the stability assessment of the seven candidate reference genes. This choice was based on the higher number of differentially transcribed fragments (DTFs) observed in cDNA-AFLP analysis at these sampling time points .
Total RNA isolation and cDNA synthesis
Total RNA was extracted from the sugarcane + 1 leaf of each biological replicate with Trizol reagent (Invitrogen, Carlsbad, USA) following the manufacturer’s instructions, and stored at − 80 °C. RNA concentration was estimated in a spectrophotometer NanoDrop2000 (Thermo Fischer Scientific, Wilmington DE, USA), and RNA integrity was checked in 1.5% agarose’s gel. Firstly, 1 μg of total RNA was treated with DNase I, following manufacturer’s instructions (Promega, Fitchburg WI, USA), to remove genomic DNA. Reverse transcription of DNase treated RNA was then performed using the GoScript Reverse Transcription System (Promega) kit, according to manufacturer’s instructions.
Candidate reference genes and primer design
The sequence of reference genes reported in sorghum (Sorghum bicolor) infected with Brome mosaic virus (BMV, Bromovirus) and in maize (Zea mays) with Barley stripe mosaic virus (BSMV, Hordeivirus), Rice black-streaked dwarf virus (RBSDV, Fijivirus) and SCMV, namely Uridylate kinase (UK), SAND protein family (SAND), and Ubiquitin-conjugating enzyme 18 gene (UBC18) , were obtained in the DFCI gene index database . These sequences were used as queries to search within the SUCEST-FUN (Sugarcane Expressed Sequence Tag Functional Analysis) database  by using BlastN tool and adopting an E-value of 1e−5 as inferior threshold. The primer design was performed using PrimerQuest tool  and analyzed using Netprimer software . The other four candidate reference genes were selected based on sugarcane gene expression studies under drought stress described by Andrade et al. : glyceraldehyde-3-phosphate dehydrogenase (GAPDH), 60S ribosomal protein L35-4 (RPL1), Actin (ACT) and Ubiquitin1 (UBQ1).
Quantitative PCR conditions
The RT-qPCR reactions were performed on an Applied Biosystems StepOnePlus System (Foster City CA, USA). The reaction mixture consisted in 5 μL of SYBR Green Power Master Mix (Applied Biosystems), 3 μL of 1:10 diluted cDNA and 0.2 μM of each forward and reverse primers in a total volume of 10 μL. The reaction thermal profile consisted in an initial denaturation step at 95 °C for 20 s, followed by 40 cycles of denaturation 95 °C for 3 s; 60 °C for 30 s. At the end of RT-qPCR reaction, dissociation curve profiles (melting curves) were carried out for amplicon specificity analysis.
Stability evaluation and selection of reference genes
The cDNA from the three aforementioned biological replicates were pooled together, resulting in eight cDNA samples, which were used in three technical replicates for the gene stability assessment (n = 24 for each gene): IAC91-1099 24 hpi (m.i), IAC91-1099 24 hpi (s.i), IAC91-1099 72 hpi (m.i), IAC91-1099 72 hpi (s.i), IACSP95-5000 24 hpi (m.i), IACSP95-5000 24 hpi (s.i), IACSP95-5000 72 hpi (m.i) and IACSP95-5000 72 hpi (s.i). PCR product threshold cycle (Ct) and PCR reaction efficiency data provided by LinReg PCR analysis  were used to seek for the best reference gene or best gene pair with NormFinder , BestKeeper  free Excel based software packages and GeNorm from NormqPCR R package , whilst RefFinder WEB-based software  identified the best reference gene based on Ct data.
Homology of sugarcane ESTs to maize and sorghum candidate reference genes and efficiency of RT-qPCR
The candidate reference genes UBC18, SAND and UK showed identity ranging from 93 to 96% and a highly significant alignment (E-value = 0) with the sugarcane ESTs (Additional file 1: Table S1). The ranking of the mean Ct-value was GAPDH > UK > SAND > UBC18 > RPL1 > ACT > UBQ1, whilst the ranking of coefficient of variation (CV%) values was SAND > UK > GAPDH > UBC18 > RPL1 > UBQ1 > ACT (Additional file 2: Figure S1; Additional file 3: Tables S2 and S3).
The LinReg PCR program, which detects the exponential phase of amplification by fluorescence data plotting in a logarithmical scale , showed amplification efficiency ranging from 90.2 to 98.2%. The matching degree of the plotted data to the standard curve in the PCR reaction, revealed by correlation coefficients (R2), ranged from 0.998 to 0.999 (Additional file 3: Tables S2 and S4). Each pair of primers showed a unique peak of fluorescence in the melting curves (Additional file 4: Figure S2), indicating single fragment amplification during RT-qPCR. The newly designed pairs of primers for genes SAND, UK and UBC18 showed the predicted amplicon size in 1% agarose gel (Additional file 5: Figure S3).
Gene expression stability
According to NormFinder the seven candidate reference genes showed the following ranking, from the most to the least stable gene: SAND > UK > GAPDH > UBC18 > RPL1 > UBQ1 > ACT. Moreover, with stability value of 0.181, SAND/GAPDH is the best combination of two genes, which represents the minimal combined intra- and intergroup variation in gene expression. The BestKeeper analysis involved two approaches, with the first, the BestKeeper standard deviation (SD) statistics, presenting the same stability ranking from above. The second approach was performed stepwise, with successive exclusion of candidate reference genes based on the SD threshold of 1.0 established by Pfaffl et al. , and on low Pearson correlation coefficient (r) values with significance cutoff at 5% level (P < 0.05). The selected genes were further ranked based on the Pearson correlation coefficients with the BestKeeper index, leading to the statement of UBC18 and UK as the most stable genes, all significant at 1% of probability. Results generated by GeNorm algorithm analysis showed the following gene ranking: UBC18/UK > SAND > GAPDH > RPL1 > UBQ1 > ACT (Table 1).
The pairwise variation Vn/Vn+1 of two sequential normalization factors NFn and NFn+1 calculated by GeNorm showed a V2/3 value of 0.08, indicating that the inclusion of a third reference gene has no significant effect for normalization, considering a threshold value below 0.15  (Fig. 1).
Among the algorithms tested in RefFinder, SAND and GAPDH were the most stable genes according to NormFinder. UBC18 and UK were the best pair of genes according to GeNorm while SAND and UK were indicated as the most stable by DeltaCt and BestKeeper algoritms. The comprehensive ranking presented the same ranking of the DeltaCt algorithm: SAND > UK > UBC18 > GAPDH > RPL1 > UBQ1 > ACT (Table 2).
The most commonly employed algorithms for reference gene expression stability analysis are based on different mathematical approaches , and often result in dissimilar outcomes. In the present study, this was observed by different statements of the most stable genes, i.e. SAND and UK, SAND/GAPDH, UBC18/UK. The comprehensive rank provided by RefFinder, allows an overall assessment of gene stability based on these different mathematical approaches , but should be restricted as a complementary tool for reference gene stability assessment taking into account the strengths and weaknesses of each algorithm . Among the tested algorithms, BestKeeper is addressed as a “common sense” between the need of reference genes with low SD values and good correlation among them, assuming that the reference genes are not co-regulated . Considering this, the Bestkeeper algorithm indicates that UK and UBC18 genes have an acceptable low SD and high correlation between them, being in agreement with GeNorm output. On the other hand, the best combination of two genes SAND/GAPDH calculated by NormFinder had low Pearson correlation coefficient (c) values according to BestKeeper, which favors the choice of UK and UBC18 genes.
The statement of genes UK and UBC18 as the most stable in SCMV-infected sugarcane by GeNorm and BestKeeper algorithms resemble the reports of Zhang et al. , e.g. GeNorm and BestKeeper outputs for BMV and BSMV-infected barley (Hordeum vulgare), and NormFinder output for BMV-infected sorghum. The UBC18 gene stability is noteworthy since ubiquitin expression has been used for normalization in maize infected by different potyviruses . The SAND and GAPDH genes also were reliably stable when subjected to all algorithms used in the present study. Similarly, SAND was reported as the most stable in wheat (Triticum aestivum) infected by BSMV and RBSDV according to NormFinder, and in BMV-infected Sorghum according to GeNorm, while GAPDH was ranked by NormFinder and GeNorm as the most stable in BSMV-infected Brachypodium (Brachypodium distachyon) . Our results rank RPL1, UBQ1 and ACT genes amongst the least stable genes by all algorithms. The poor transcript stability of ACT is in agreement with previous studies [16, 34, 35], while the report of 60 s ribosomal protein in replication complexes during potyvirus infection  seem to corroborate our observations for RPL1. The sugarcane UBQ1 gene contrasted with UBC18, suggesting that potyviruses may interfere with pathways involving certain ubiquitin genes as reported by Cheng; Wang .
The results indicate that UBC18 and UK are the most stable sugarcane reference genes in leaves when the target is gene expression studies in search for resistance to SCMV by RT-qPCR approaches, and should also be considered as candidate reference genes for accurate normalization for other expression studies involving SMD.
It is necessary to reassess expression stability of candidate reference genes when different experimental conditions and genotypes are involved in SMD studies. In addition, an important step for the selection of reference genes is the validation by RT-qPCR analysis of a well-studied sugarcane gene responsive to SCMV infection, which information is lacking in literature. Therefore, studies with good candidate genes, e.g. recent reports in maize [38, 39] could provide useful data.
sugarcane mosaic disease
hours post inoculation
Zhai Y, Deng Y, Cheng G, Peng L, Zheng Y, Yang Y, Xu J. Sugarcane Elongin C is involved in infection by sugarcane mosaic disease pathogens. Biochem Biophys Res Commun. 2015;466(3):312–8. https://doi.org/10.1016/j.bbrc.2015.09.015.
Gonçalves MC, Pinto LR, Souza SC, Landell MGA. Virus diseases of sugarcane. A constant challenge to sugarcane breeding in Brazil. Funct Plant Sci Biotechnol. 2012;6:108–16.
Camelo-García VM, da Silva Andrade SC, Geering AD, Kitajima EW, Rezende JA. Genome organization and host range of a Brazilian isolate of Johnsongrass mosaic virus. Arch Virol. 2016;161(5):1335–41. https://doi.org/10.1007/s00705-016-2772-4.
Manners JM, Casu RE. Transcriptome analysis and functional genomics of sugarcane. Trop Plant Biol. 2011;4:9–21. https://doi.org/10.1007/s12042-011-9066-5.
Schaker PD, Palhares AC, Taniguti LM, Peters LP, Creste S, Aitken KS, Van Sluys MA, Kitajima JP, Vieira ML, Monteiro-Vitorello CB. RNAseq transcriptional profiling following whip development in sugarcane smut disease. PLoS ONE. 2016;11(9):e0162237. https://doi.org/10.1371/journal.pone.0162237.
Vicentini R, Bottcher A, Brito MOS, Dos Santos AB, Creste S, Landell MG, Cesarino I, Mazzafera P. Large-scale transcriptome analysis of two sugarcane genotypes contrasting for lignin content. PLoS ONE. 2015;10(8):e0134909. https://doi.org/10.1371/journal.pone.0134909.
Medeiros CNF, Goncalves MC, Harakava R, Creste S, Nobile PM, Pinto LR, Perecin D, Landell MGA. Sugarcane transcript profiling assessed by cDNA-AFLP analysis during the interaction with Sugarcane mosaic virus. Adv Microbiol. 2014;4:511–20. https://doi.org/10.4236/aim.2014.49057.
Dong M, Cheng G, Peng L, Xu Q, Yang Y, Xu J. Transcriptome analysis of sugarcane response to the infection by Sugarcane Streak Mosaic Virus (SCSMV). Trop Plant Biol. 2017;10:45. https://doi.org/10.1007/s12042-016-9183-2.
Nolan T, Hands RE, Bustin SA. Quantification of mRNA using real-time RT-PCR. Nat Protoc. 2006;1(3):1559–82. https://doi.org/10.1038/nprot.2006.236.
Sanders R, Mason DJ, Foy CA, Huggett JF. Considerations for accurate gene expression measurement by reverse transcription quantitative PCR when analysing clinical samples. Anal Bioanal Chem. 2014;406(26):6471–83. https://doi.org/10.1007/s00216-014-7857-x.
Gusberti M, Gessler C, Broggini GA. RNA-Seq analysis reveals candidate genes for ontogenic resistance in Malus-Venturia pathosystem. PLoS ONE. 2013;8(11):e78457. https://doi.org/10.1371/journal.pone.0078457.
Nakayama T, Okada N, Yoshikawa M, Asaka D, Kuboki A, Kojima H, Tanaka Y, Haruna SI. Assessment of suitable reference genes for RT-qPCR studies in chronic rhinosinusitis. Sci Rep. 2018;8(1):1568. https://doi.org/10.1038/s41598-018-19834-9.
Guénin S, Mauriat M, Pelloux J, Van Wuytswinkel O, Bellini C, Gutierrez L. Normalization of qRT-PCR data: the necessity of adopting a systematic, experimental conditions-specific, validation of references. J Exp Bot. 2009;60(2):487–93. https://doi.org/10.1093/jxb/ern305.
Lacerda AL, Fonseca LN, Blawid R, Boiteux LS, Ribeiro SG, Brasileiro AC. Reference gene selection for qPCR analysis in Tomato-Bipartite Begomovirus interaction and validation in additional Tomato-Virus Pathosystems. PLoS ONE. 2015;10(8):e0136820. https://doi.org/10.1371/journal.pone.0136820.
Andrade LM, Dos Santos Brito M, Fávero Peixoto Junior R, Marchiori PER, Nóbile PM, Martins APB, Ribeiro RV, Creste S. Reference genes for normalization of qPCR assays in sugarcane plants under water deficit. Plant Methods. 2017;13:28. https://doi.org/10.1186/s13007-017-0178-2.
Zhang K, Niu S, Di D, Shi L, Liu D, Cao X, Miao H, Wang X, Han C, Yu J, et al. Selection of reference genes for gene expression studies in virus-infected monocots using quantitative real-time PCR. J Biotechnol. 2013;168(1):7–14. https://doi.org/10.1016/j.jbiotec.2013.08.008.
Yang ZN, Mirkov TE. Sequence and relationships of sugarcane mosaic and sorghum mosaic virus strains and development of RT-PCR-based RFLPs for strain discrimination. Phytopathology. 1997;87(9):932–9. https://doi.org/10.1094/PHYTO.19184.108.40.2062.
Gonçalves MC, Santos AS, Maia LG, Chagas CM, Harakava R. Caracterização de um isolado do Sugarcane mosaic virus que quebra a resistência de variedades comerciais de cana-de-açúcar. Fitopatologia Brasileira. 2007;32:32–9. https://doi.org/10.1590/S0100-41582007000100004.
Bain DC. The use of abrasive for inoculating sugarcane seedlings with the mosaic virus. Phytopathology. 1944;34:844–5.
Quackenbush J, Antonescu C, Antonescu V, Cho J, Holt I, Karamycheva S, et al. The TGI databases; 2014. http://compbio.dfci.harvard.edu/tgi/. Accessed 27 Ago 2016.
Sugarcane Functional Genomics Database (SUCEST-FUN). http://sucest-fun.org/. Accessed 26 Jan 2017.
Integrated DNA Technologies (IDT): PrimerQuest Tool; 2012. https://www.idtdna.com/primerquest/Home/Index. Accessed 26 Jan 2017.
PREMIER Biosoft: Netprimer. http://www.premierbiosoft.com/netprimer/. Accessed 26 Jan 2017.
Ramakers C, Ruijter JM, Deprez RH, Moorman AF. Assumption-free analysis of quantitative real-time polymerase chain reaction (PCR) data. Neurosci Lett. 2003;339(1):62–6. https://doi.org/10.1016/S0304-3940(02)01423-4.
Andersen CL, Jensen JL, Ørntoft TF. Normalization of real-time 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. 2004;64(15):5245–50. https://doi.org/10.1158/0008-5472.CAN-04-0496.
Pfaffl MW, Tichopad A, Prgomet C, Neuvians TP. Determination of stable housekeeping genes, differentially regulated target genes and sample integrity: bestKeeper–Excel-based tool using pair-wise correlations. Biotechnol Lett. 2004;26(6):509–15. https://doi.org/10.1023/B:BILE.0000019559.84305.47.
Perkins JR, Dawes JM, McMahon SB, Bennett DL, Orengo C, Kohl M. ReadqPCR and NormqPCR: R packages for the reading, quality checking and normalisation of RT-qPCR quantification cycle (Cq) data. BMC Genomics. 2012;13:296. https://doi.org/10.1186/1471-2164-13-296.
Zhang B, Xie F: RefFinder; 2012. http://leonxie.esy.es/RefFinder/?type=reference. Accessed 18 Jul 2017.
Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, Speleman F. Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002;3(7):RESEARCH0034.
Robledo D, Hernández-Urcera J, Cal RM, Pardo BG, Sánchez L, Martínez P. Analysis of qPCR reference gene stability determination methods and a practical turbot (Scophthalmus maximus) gonad dataset. BMC Genomics. 2014;15:648. https://doi.org/10.1186/1471-2164-15-648.
Xie F, Xiao P, Chen D, Xu L, Zhang B. miRDeepFinder: a miRNA analysis tool for deep sequencing of plant small RNAs. Plant Mol Biol. 2012. https://doi.org/10.1007/s11103-012-9885-2.
De Spiegelaere W, Dern-Wieloch J, Weigel R, Schumacher V, Schorle H, Nettersheim D, Bergmann M, Brehm R, Kliesch S, Vandekerckhove L, et al. Reference gene validation for RT-qPCR, a note on different available software packages. PLoS ONE. 2015;10(3):e0122515. https://doi.org/10.1371/journal.pone.0122515.
Shi Y, Qin YH, Cao YY, Sun H, Zhou T, Hong YG, Fan ZF. Influence of an m-type thioredoxin in maize on potyviral infection. Eur J Plant Pathol. 2011;131:317–26. https://doi.org/10.1007/s10658-011-9810-6.
Lilly ST, Drummond RS, Pearson MN, MacDiarmid RM. Identification and validation of reference genes for normalization of transcripts from virus-infected Arabidopsis thaliana. Mol Plant Microbe Interact. 2011;24(3):294–304. https://doi.org/10.1094/MPMI-10-10-0236.
Liu D, Shi L, Han C, Yu J, Li D, Zhang Y. Validation of reference genes for gene expression studies in virus-infected Nicotiana benthamiana using quantitative real-time PCR. PLoS ONE. 2012;7(9):e46451. https://doi.org/10.1371/journal.pone.0046451.
Lõhmus A, Varjosalo M, Mäkinen K. Protein composition of 6K2-induced membrane structures formed during Potato virus A infection. Mol Plant Pathol. 2016;17(6):943–58. https://doi.org/10.1111/mpp.12341.
Cheng X, Wang A. The potyvirus silencing suppressor protein VPg mediates degradation of SGS3 via ubiquitination and autophagy pathways. J Virol. 2017. https://doi.org/10.1128/jvi.01478-16.
Liu Q, Liu H, Gong Y, Tao Y, Jiang L, Zuo W, Yang Q, Ye J, Lai J, Wu J, et al. An atypical thioredoxin imparts early resistance to sugarcane mosaic virus in Maize. Mol Plant. 2017;10(3):483–97. https://doi.org/10.1016/j.molp.2017.02.002.
Zhu M, Chen Y, Ding XS, Webb SL, Zhou T, Nelson RS, Fan Z. Maize Elongin C interacts with the viral genome-linked protein, VPg, of Sugarcane mosaic virus and facilitates virus infection. New Phytol. 2014;203(4):1291–304. https://doi.org/10.1111/nph.12890.
LRP, MCG, MSB, PMN and LMA were responsible for the design and conceptualization of the experiment. SC produced the sugarcane virus-free plantlets obtained by meristem tip culture. CNM conducted the greenhouse experiment.MFS designed the primers used for RT-qPCR, performed the RT-qPCR analysis and wrote the manuscript. MFS, LMA, MCG, MSB, LRP and PMN reviewed and edited the manuscript. LRP and MCG provided funding acquisition for the research. All authors read and approved the final manuscript.
The authors declare that they have no competing interest.
Availability of data and materials
The conclusions of this study were based on the datasets included within the article and in additional files.
Consent for publication
Ethics approval and consent to participate
This work was supported by “Conselho Nacional de Desenvolvimento Científico e Tecnológico” (CNPq—Edital Universal, Award Number 424050/2016-1) and the “Fundação de Amparo à Pesquisa do Estado de São Paulo” (FAPESP, Brazil-2012/15060-6). MFS was supported by a Doctor fellowship from “Coordenação de Aperfeiçoamento de Pessoal de Nível Superior” CAPES.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
da Silva, M.F., Gonçalves, M.C., dos Santos Brito, M. et al. Reference genes for gene expression studies targeting sugarcane infected with Sugarcane mosaic virus (SCMV). BMC Res Notes 12, 149 (2019). https://doi.org/10.1186/s13104-019-4168-5
- Saccharum spp.
- Sugarcane mosaic disease