Abstract
In this study, we quantify and document the distribution and organization of cpSSRs in the chloroplast genome of 31 Eucalyptus species. Our sample included all previously sequenced plastomes of Eucalyptus species available from the NCBI online database. We processed the complete cpDNA sequences and identified mono-, di-, tri-, tetra-, penta-, and hexanucleotide cpSSRs, with the majority of cpSSRs classified as mononucleotide. After genome microsatellite selection, we evaluated the microsatellites for coding and non-coding regions and cpSSRs were predominantly identified in non-coding regions of cpDNA for all nucleotide types. Penta- and hexanucleotide cpSSRs were the least frequent types of microsatellites. We also developed and virtually amplified 60 primers pairs that can be used in studies of Eucalyptus species. Thus, these cpSSR regions can be used in studies assessing the ecology, breeding, and conservation of the genus.
Key words:
Chloroplast DNA; DNA fingerprint; genetic diversity; species differentiation.
INTRODUCTION
Eucalyptus L'Hér. is the most dominant genus in Australian flora and with approximately 700 species, it is the largest genus in the Myrtaceae family (Smith et al. 2003). The genus originates from the region in and around Australia, including Timor, Indonesia, Papua New Guinea, the Molucas Islands, Irian Jaya, and the Southern Philippines, between 9° N and 44° S latitude (Eldridge et al. 1993). In 2015, Eucalyptus forest plantations represented about 10% of the approximately 291 million hectares of planted forests worldwide (FAO 2016, Sands 2005).
Because of its extensive use in forestry, plant ecologists, breeders and geneticists have significant interest in establishing methods for fast and robust Eucalyptus identification or species distinction. Currently, the use of marker-assisted approaches in Eucalyptus breeding programs is common (Mahajan and Gupta 2012, Fuchs et al. 2015). These approaches are powerful tools for phylogenetic and population research (Doorduin et al. 2011) and can be used to assist ecologists and breeders in their analyses (Steane et al. 1998, Freeman et al. 2001, McKinnon et al. 2001, McKinnon et al. 2010, Melotto-Passarin et al. 2011).
The development of molecular markers for Eucalyptus species improved dramatically after the complete sequencing of its nuclear genome (nDNA) (Myburg et al. 2014, Fuchs et al. 2015). Currently, microsatellite markers (SSRs) are widely used in. Although many nSSRs have been developed for the Eucalyptus genus (Kirst et al. 2005), the use of cpSSRs is still limited. Nevertheless, as with nDNA, cpDNA (i.e., cpSSR) has demonstrated great potential for use in phylogeography, analyses of DNA fingerprints, hybridization and progeny analyses in tree species (Steane et al. 2005, Delgado et al. 2007, McKinnon et al. 2010). Because the majority of cpDNA evolves slowly, its overall utility for evolutionary and population genetic studies can be limited, particularly at lower taxonomic levels (Steane et al. 1991, Steane et al. 1998, Freeman et al. 2001, Shaw et al. 2005). As a result, many researchers have sought to use hypervariable regions of cpDNA to characterize the genetic variation of species (Shaw et al. 2007). In this context, cpSSRs have become a significant resource for genetic analyses (Weising and Gardner 1999) focused on ecology, systematics, conservation, phylogenetic studies, and in breeding programs. As such, chloroplast genomes (cpDNA) can provide important data that can be used in ecological (Latouche-Hallé et al. 2003) and breeding studies of Eucalyptus species, such as estimating the phylogenetic relationships within and between species (Whitock et al. 2003).
In this study, we aim to quantify and document the distribution and organization of cpSSRs in all 31 species of Eucalyptus that have an available plastome sequence. We also developed 60 potential primer pairs that can be used by researchers in a variety of different fields of study.
MATERIAL AND METHODS
In silico identification of microsatellite
This study is based on complete plastid genome sequences available from the National Center for Biotechnology Information (NCBI; https://www.ncbi.nlm.nih.gov/genomes/GenomesGroup.cgi?opt=plastid&taxid=2759). From this database, all 31 chloroplast genomes of Eucalyptus spp. were used (Table 1). The complete cpDNA sequences were collected from NCBI and processed using the FastPCR 6.5.40.0 software (Kalendar et al. 2016) to identify the microsatellite regions (cpSSR) as mono-, di-, tri-, tetra-, penta-, or hexanucleotide. We considered only those repeats in which the motifs repeated as follows: mononucleotide repeats with a repeat length ≥ 8; dinucleotide with a repeat length ≥ 6; and tri-, tetra-, penta-, and hexanucleotide with a repeat length ≥ 3. After identification, we assessed microsatellites for coding (gene) and non-coding (intergenic) regions according the information available in the NCBI database for each species.
Primer design and in silico PCR amplification
The complete chloroplast genome of Eucalyptus grandis (GenBank NC_014570) was screened for cpSSRs and 60 cpSSR flanking hypervariable regions were used to develop primer pairs. Primer pairs were designed using Primer3 (v. 0.4.0) (Untergasser et al. 2012) based on the following criteria: annealing temperature between 57 and 63 °C; final amplification product of 100 - 350bp; GC content of 20 - 80%; and primer size of 18 - 27bp.
Virtual amplification was performed using the Sequence Manipulation Suite (http://www.bioinformatics.org/sms2/pcr_products.html; Stothard 2000) for all 31 Eucalyptus species. The Suite determines the expected product size for each primer pair.
RESULTS AND DISCUSSION
All microsatellite regions in the cpDNA for the 31 studied Eucalyptus species were evaluated (Table 2). The total number of microsatellites observed for each species varied from 71 for E. melliodora to 135 for E. aromaphloia.
While we identified mono-, tri-, and tetranucleotide microsatellites for all species (Table 2), the majority of Eucalyptus cpSSRs were classified as mononucleotide. Currently there are no other studies that assess the distribution and organization of cpSSRs in Eucalyptus species; therefore, we provide a general comparison between our results and published studies on other plant species. For example, Vieira et al. (2015) analyzed the occurrence, type, and distribution of cpSSRs for 20 Bambuseae species and their results showed an average total of 141.8 cpSSRs, with a predominance of mono- and dinucleotide repeats. Similarly, Vieira et al. (2016a), studying Retrophyllum piresii Silba C.N. Page, reported that 94.5% of total cpSSRs were mono- and dinucleotides. Furthermore, George et al. (2015), studying the abundance and distribution of SSRs in 164 sequenced plastomes from a wide range of plants, demonstrate that plant species predominantly present mononucleotide repeats, as observed herein for Eucalyptus species.
Chloroplasts are well understood to originate from endosymbiosis between cyanobacteria and higher plants (McFadden 2001, Raven and Allen 2003). The theory of endosymbiosis argues that cyanobacteria, also called blue-green algae, are the ancestors of chloroplasts, despite the fact that cyanobacteria are prokaryotes (Clegg et al. 1994). For the studied Eucalyptus species, most mononucleotide cpSSRs were found in non-coding regions of cpDNA, with the exception of E. nitens and E. marginata which present a more balanced proportion of mononucleotide microsatellites in coding and non-coding regions (Table 2). Gandhi et al. (2010) evaluated cpSSR dynamics in 12 species of the Brassicaceae family using data collected from GenBank. Their results revealed that, although 51% of cpDNA occurs in coding regions, the majority of the total number of SSRs was in non-coding regions that contain twice the number of SSRs than those in coding cpDNA regions.
Herein, dinucleotide microsatellites were found in cpDNA of 15 species, all in non-coding regions (Table 2). Tambarussi et al. (2009), Gandhi et al. (2010) and Melotto-Passarin et al. (2011) observed similar results for Solanaceae, Brassicaceae and Poaceae families, respectively. We also found that penta- and hexanucleotide microsatellites occur predominantly in non-coding regions and these two motifs were the least frequent types of microsatellites for all studied Eucalyptus species. Penta- and hexanucleotide cpSSRs were observed only in E. baxteri and E. diversicolor, respectively (Table 2), all in non-coding regions. Similar to the results presented herein, penta- and hexanucleotide repeats occurred at a low frequency in Bambuseae species, but in coding sequences (Vieira et al. 2015). For cpDNA of R. piresii, Vieira et al. (2016a) found no penta- or hexanucleotide cpSSRs, which is consistent with the results found for most of the studied Eucalyptus species.
Sixty primer pairs were developed and only 20% (12) of them occur in genic regions. All primers virtually amplified products around the expected size and displayed interspecific polymorphism. According to Sumathi and Yasodha (2014), until 2014 only 35 cpSSRs had been identified for different species of Eucalyptus. Therefore, our results indicate a wide range of regions in cpDNA that have yet to be explored (Doorduin et al. 2011). As cpDNA markers are nonrecombinant, uniparentally inherited, and present low rates of mutation in the plastome sequence (Provan et al. 2001, Vieira et al. 2016b), they can be used in phylogenetic analyses and breeding and conservation programs for this genus.
Several studies have indicated that chloroplast genomes and cpSSRs have enormous potential (Zhang et al. 2017) and they are an important tool for plant biologists and breeders in assessing genetic diversity (Doorduin et al. 2011), spatial genetic structure (Islam et al. 2015), evolutionary history, and hybridization in native and improved species (Ebert and Peakall 2009, Doorduin et al. 2011). Thus, our analysis provides a rich database for researchers that can be used to study these Eucalyptus species across a wide range of scientific disciplines.
ACKNOWLEDGEMENTS
The authors wish to thank CNPq, Fundação Arauácia and CAPES for scholarships attributed to the first, second and third author, respectively. The authors are very grateful to the anonymous referees for their important suggestions in the data analysis. We also thank Dr. Evelyn R. Nimmo for editing the English of the manuscript.
DATA ARCHIVING STATEMENT
All cpSSRs are provided as additional files when requested to the Correspondent Author.
REFERENCES
- Clegg MT, Gaut BS, Learn Jr BH and Morton BR (1994) Rates and patterns of 643 chloroplast DNA evolution. Proceedings of the National Academy of Sciences 91: 6795-6801.
- Delgado P, Salas-Lizana R, Vázquez-Lobo A, Wegier A, Anzidei M, Alvarez Buylla ER, Vendramin GG and Piñero D (2007) Introgressive hybridization inPinus montezumaeLamb andPinus pseudostrobusLindl. (Pinaceae): Morphological and molecular (cpSSR) evidence. Journal of Plant Science 168: 861-875.
- Doorduin L, Gravendeel B, Lammers Y, Ariyurek Y, Chin-A-Woeng T and Vrieling K (2011) The complete chloroplast genome of 17 individuals of pest speciesJacobaea vulgaris: SNPs, microsatellites and barcoding markers for population and phylogenetic studies. DNA Research 18: 93-105.
- Ebert D and Peakall R (2009) Chloroplast simple sequence repeats (cpSSRs): technical resources and recommendations for expanding cpSSR discovery and applications to a wide array of plant species. Molecular Ecology Resources 9: 673-690.
- Eldridge K, Davidson J, Harwood C and Van Wyk G (1993) Eucalypt domestication and breeding. Claredon Press, Oxford, 312p.
- FAO (2016) Global forest resources assessment 2015. 2nd edn, FAO, Rome , 54p.
- Freeman JS, Jackson HD, Steane DA, McKinnon GE, Dutkowski GW, Potts BM and Vaillancourt RE (2001) Chloroplast DNA phylogeography of Eucalyptus globulus Australian Journal of Botany 49: 585-596.
- Fuchs MCP, Tambarussi EV, Lourenção JC, Nogueira LM, Bortoloto TM, González ER, Oda S and Marino CL (2015) Molecular marker associated with a deleterious recessive anomaly in Eucalyptus grandis seedlings. Annals of Forest Science 72: 1043-1052.
- Gandhi SG, Awasthi P and Bedi YS (2010) Analysis of SSR dynamics in chloroplast genomes of Brassicaceae family. Bioinformation 5: 16-20.
- George B, Bhatt BS, Awasthi M, George B and Singh AK (2015) Comparative analysis of microsatellites in chloroplast genomes of lower and higher plants. Current Genetics 61: 665-677.
- Islam MS, Lian C, Kameyama N and Hogetsu T (2015) Analysis of the mating system, reproductive characteristics, and spatial genetic structure in a natural mangrove tree (Bruguiera gymnorrhiza) population at its northern biogeographic limit in the southern Japanese archipelago. Journal of Forest Research 20: 293-300.
-
Kalendar R, Lee D and Schulman AH (2016) FastPCR software for PCR primer and probe design and repeat search. Genes, Genomes and Genomics 3. Available at <Available at www.biocenter.helsinki.fi/bi/Programs/fastpcr.htm
>. Accessed on Feb 16, 2016.
» www.biocenter.helsinki.fi/bi/Programs/fastpcr.htm - Kirst M, Cordeiro CM, Rezende GDSP and Grattapalglia D (2005) Power of microsatellite markers for fingerprinting and parentage analysis in Eucalyptus grandis breeding population. Journal of Heredity 96: 1-6.
- Latouche-Hallé C, Ramboer A, Bandou E, Caron H and Kremer A (2003) Nuclear and chloroplast genetic structure indicate fine-scale spatial dynamics in a Neotropical tree population. Heredity 91: 181-190.
- Mahajan R and Gupta P (2012) Molecular markers: their use in tree improvement. Journal of Forest Science 58: 137-144.
- McFadden GI (2001) Chloroplast origin and integration. Plant Physiology 125: 50-53.
- McKinnon GE, Smith JJ and Potts BM (2010) Recurrent nuclear DNA introgression accompanies chloroplast DNA exchange between two eucalypt species. Molecular Ecology 19: 1367-1380.
- McKinnon GE, Vaillancourt RE, Jackson HD and Potts BM (2001) Chloroplast sharing in the Tasmanian eucalypts. Evolution 55: 703-711.
- Melotto-Passarin DM, Tambarussi EV, Dressano K, De Martin VF and Carrer H (2011) Characterization of chloroplast DNA microsatellites from Saccharum spp. and related species. Genetics and Molecular Research 10: 2024-2033.
- Myburg AA, Grattapaglia D, Tuskan GA, Hellsten U, Hayes RD, Jane Grimwood J, Jenkins J, Lindquist E, Tice H, Bauer D, Goodstein DM, Dubchak I, Poliakov A, Mizrachi E, Kullan ARK,Hussey SG, Pinard D, Merwe K, Singh P, Jaarsveld I,Silva-Junior OB, Togawa RC, Pappas MR, Faria DA, Sansaloni CP, Petroli CD, Yang X, Ranjan P, Tschaplinski TJ, Ye C-Y, Li T, Sterck L, Vanneste K, Murat F, Soler M, Clemente HS, Saidi N, Cassan-Wang H, Dunand C, Hefer CA, Bornberg-Bauer E, Kersting AR, Vining K, Amarasinghe V, Ranik M, Naithani S, Elser J, Boyd AE, Liston A, Spatafora JW, Dharmwardhana P, Raja R, Sullivan C, Romanel E, Alves-Ferreira M, Külheim C, Foley W, Carocha V, Paiva J, Kudrna D, Brommonschenkel SH, Pasquali G, Byrne M, Rigault P, Tibbits J, Spokevicius A, Jones RC, Steane DA, Vaillancourt RE, Potts BM, Joubert F, Barry K, Pappas GJ, Strauss SH, Jaiswal P, Grima-Pettenati J, Salse J, Peer YV, Rokhsar DS and Schmutz J (2014) The genome of Eucalyptus grandis Nature 510: 356-362.
- Provan J, Powell W and Hollingsworth PM (2001) Chloroplast microsatellites: new tools for studies in plant ecology and evolution. Trends in Ecology & Evolution 16: 142-147.
- Raven JA and Allen JF (2003) Genomics and chloroplast evolution: what did cyanobacteria do for plants? Genome Biology 4: 209.
- Sands R (2005) Forestry in a global context. CABI Publishing, United Kingdom, 233p.
- Shaw J, Lickey EB, Beck JT, Farmer SB, Liu W, Miller J, Siripun KC, Winder CT, Schilling EE and Small RL (2005) The torroise and the hare II: relative utility of 21 noncoding chloroplast DNA sequences for phylogenetic analysis. American Journal of Botany 92: 142-166.
- Shaw J, Lickey EB, Schilling EE and Small RL (2007) Comparison of whole chloroplast genome sequences to choose noncoding regions for phylogenetic studies in angiosperms: the tortoise and the hare III analysis. American Journal of Botany 94: 275-288.
- Smith S, Hughes J and Wardell-Johnson G (2003) High population differentiation and extensive clonality in a rare mallee eucalypt: Eucalyptus curtisii Conservation Genetics 4: 289-300.
- Steane DA, Byrne M, Vaillancourt RE and Potts BM (1998) Chloroplast DNA polymorphism signals complex interspecific interactions in Eucalyptus (Myrtaceae). Australian Systematic Botany 11: 25-40.
- Steane DA, Jones RC and Vaillancourt RE (2005) A set of chloroplast microsatellite primers for Eucalyptus (Myrtaceae). Molecular Ecology Notes 5: 538-541.
- Steane DA, West AK, Potts BM, Ovenden JR and Reid JB (1991) Restriction fragment length polymorphisms in chloroplast DNA from six species of Eucalyptus Australian Journal of Botany 39: 399-414.
- Stothard P (2000) The Sequence Manipulation Suite: JavaScript programs for analyzing and formatting protein and DNA sequences. Biotechniques 28: 1102-1104.
- Sumathi M and Yasodha R (2014) Micorsatellite resources of Eucalyptus: current status and future perspectives. Botanical Studies 55: 73.
- Tambarussi EV, Melotto-Passarin DM, Dressano K, Guidetti-Gonzalez S, Brigati JB, Jesus FA, Barbosa AL and Carrer H (2009) In silico analysis of Simple Sequence Repeats from chloroplast genomes of Solanaceae species. Crop Breeding and Applied Biotechnology 9: 342-352.
- Untergasser A, Cutcutache I, Koressaar T, Ye J, Faircloth BC, Remm M and Rozen SG (2012) Primer3 - new capabilities and interfaces. Nucleic Acids Research 40: e 115.
- Vieira LN, Anjos KG, Faoro H, Fraga HPF, Greco TM, Pedrosa FO, Souza EM, Rogalski M, Souza RF and Guerra MP (2015) Phylogenetic inference and SSR characterization of tropical woody bamboos tribe Bambuseae (Poaceae: Bambusoideae) based on complete plastid genome sequences. Current Genetics 62: 443-453.
- Vieira LN, Rogalski M, Faoro H, Fraga HPF, Anjos KG, Picchi GFA, Nodari RO, Pedrosa FO, Souza EM and Guerra MP (2016a) The plastome sequence of the endemic Amazonian conifer, Retrophyllum piresii (Silba) C.N.Page, reveals different recombination events and plastome isoforms. Tree Genetics & Genomes 12: 10.
- Vieira MLCV, Santini L, Diniz AL and Munhoz CF (2016b) Microsatellite markers: what they mean and why they are so useful. Genetics & Molecular Biology 39: 312-328.
- Weising K and Gardner RC (1999) A set of conserved PCR primers for the analysis of simple sequence repeat polymorphisms in chloroplast genomes of dicotyledonous angiosperms. Genome 42: 9-19.
- Whitock S, Steane SA, Vaillancourt RE and Potts BM (2003) Molecular evidence shows that the tropical boxes (Eucalyptus subgenus Minutifructus) are over-ranked. Transactions of Royal Society of South Australian 127: 27-32.
- Zhang H, Hall N, McElroy JS, Lowe EK and Goertzen LR (2017) Complete plastid genome sequence of goosegrass (Eleusine indica) and comparison with other Poaceae. Gene 5: 36-43.
