African Journal of
Agricultural Research

  • Abbreviation: Afr. J. Agric. Res.
  • Language: English
  • ISSN: 1991-637X
  • DOI: 10.5897/AJAR
  • Start Year: 2006
  • Published Articles: 6894

Full Length Research Paper

In silico identification of putative expressed sequence tag (EST)-simple sequence repeats (SSRs) markers of resistance to Meloidogyne spp. in common bean

Lucas Donizetti Vieira
  • Lucas Donizetti Vieira
  • Department of Genetics and Molecular Biology, Universidade Federal de Goiás, Goiânia, Brazil.
  • Google Scholar
Juliana Oliveira da Silva
  • Juliana Oliveira da Silva
  • Department of Agronomy, IF Goiano, Urutaí, Goiás, Brazil.
  • Google Scholar
Caio César de Oliveira Pereira
  • Caio César de Oliveira Pereira
  • Department of Agronomy, IF Goiano, Urutaí, Goiás, Brazil.
  • Google Scholar
Solange Aline de Carvalho
  • Solange Aline de Carvalho
  • Department of Biology, IF Goiano, Urutaí, Goiás, Brazil.
  • Google Scholar
Ricardo Diogenes Dias Silveira
  • Ricardo Diogenes Dias Silveira
  • Department of Biology, IF Goiano, Urutaí, Goiás, Brazil.
  • Google Scholar
Guilherme Malafaia
  • Guilherme Malafaia
  • Department of Biology, IF Goiano, Urutaí, Goiás, Brazil.
  • Google Scholar
Ivandilson Pessoa Pinto de Menezes
  • Ivandilson Pessoa Pinto de Menezes
  • Department of Biology, IF Goiano, Urutaí, Goiás, Brazil.
  • Google Scholar


  •  Received: 10 April 2016
  •  Accepted: 16 May 2016
  •  Published: 09 June 2016

 ABSTRACT

Expressed sequences are important sources in the development of functional heterologous microsatellite markers in phylogenetic related groups, that is, soybean and commom bean. The objective of this work was to identify and characterize expressed sequence tag (EST)-simple sequence repeats (SSRs) in silico candidates of resistance to Meloidogyne spp. in common bean. Seven DNA sequences from soybean associated with genetic resistance have been identified and obtained in the NCBI database. Its homology in common bean genome was verified using the BLAST tool. The cellular processes involved were also checked using the Blast2GO program. The identification of microsatellite markers and design of the primer pairs was performed using the SSRLocator and Primer3 programs, respectively. The transferability rate of common beans the target sequences identified was 86%, demonstrating the power of success of this method. All the cellular processes involved in the original DNA sequences were verified from EST on beans, with E-value between 0 and 2.9×10-168. Fifteen EST-SSRs candidates for common bean resistance were identified, which have proved to be suitable for their amplification by PCR. The transferability analysis of ESTs related to resistance to Meloidogyne spp., especially among soybeans and common beans is efficient. Based on this study, 15 EST-SSRs candidates are available for validation and later use.

Key words: simple sequence repeats (SSRs) markers, assisted selection, ontology, pre-breeding.


 INTRODUCTION

The common bean (Phaseolus vulgaris) is a leguminous vegetable cultivated in 113 tropical and subtropical countries, which has a high economic value and market acceptance due  to  its  balanced  chemical  composition, including various proteins, complex carbohydrates, minerals and vitamin B complex (Broughton et al., 2003).
 
Brazil is the world leader in the production and consumption of beans, although it is not yet self-sufficient in the production of the same (Wander, 2005). The bean is cultivated in most Brazilian states, in a variety of soils, climates, growing season and farming systems (Silva and Wander, 2013). The adversities which the beans have been exposed contribute to high diversity found in the culture, enabling the selection, that is, genetic variation for breeding (Burle et al., 2010).
 
This availability of variability generates demand of detailed information about their genetic resources for use, which will be required to meet a process of expansion and growing interest in improved cultivars increasingly productive and adapted to adverse conditions of development. Another point required and of greatest interest to breeders is the introduction of desired genetic variation in the Active Germplasm Bank (AGB) on which lineages are selected, reducing the segregation of characters already improved. The assisted selection allied to the DNA markers is a useful method for this purpose.
 
Over the past decades, bean culture is going through several changes in the form of management, production systems and breeding. Coupled with the landscape changes for planting, the various plant health problems arise, making diseases, caused by plant pathogens, to cause a strong deterrent to obtain better yields in the country in diverse cultures, like the beans (Vieira et al., 1998). These changes contributed to the introduction of emerging pathogens, and old diseases that were not considered problems for the culture started to be highlighted. An important example is the root-knot nematodes, considered bean pests only in some producing regions of Brazil. Though the diseases caused by nematodes are not widely distributed in the Brazilian states, their damage lead to the Brazilian production not meeting the supply demand of the domestic market (Yokoyama, 2007).
 
An important form of management of nematodes is the use of resistant varieties, which, despite having been found for Meloidogyne spp. some find themselves still not very exploited in bean crop (Carneiro et al., 1992; Walber et al., 2003), a fact that underscores the importance in the development of tools that help produce detailed information for using these varieties, in addition define preventive recommendations. In this context, DNA markers are biotechnological tool of quick evaluation that can help design strategies of use of available genetic resources.
 
The DNA markers, such as Simple Sequence Repeats (SSR) and Single Nucleotide Polymorphism (SNP), have often been used in bean breeding programs (Müller et al., 2015; Müller et al., 2014). Its use has involved studies of quantification and determination of gene diversity (Cardoso et al., 2014); characterization of Value for Cultivation and Use (VCU) (Cardoso et al., 2013), assisted selection (Alzate-Marian et al., 2005), among others. The SSRs have been widely used in recent years and had selected population studies  and  use  of  genetic resources in plants for breeding purposes because of their hypervariability and ease polymorphism detection (Li et al., 2015a; Li et al., 2015b; Kaur et al., 2015; Zhai et al., 2015).
 
Obtaining selective microsatellite markers developed from databases, get easier, faster and economic compared to the ones developed by conventional methods (Buso et al., 2006). Thus, the recent full genomic sequencing advances and the increase in banks expressed sequence tag (EST) of different species, has facilitated studies of genetic similarity between phylogenetically related species (Schmutz et al., 2014; Wang et al., 2014). With this study we propose to identify and characterize in silico microsatellite markers candidates in P. vulgaris associated with genes that confer resistance to root-knot nematode (Meloidogyne spp.). The product generated from this study could be used as a preventive polymorphism identification tool related to the phenotypic behavior of bean caused by parasitism.


 MATERIALS AND METHODS

Initially this study used secondary source materials to find articles covering genes or DNA sequences identified that confer a phenotypic behavior of resistance to nematodes belonging to Meloidogyne spp. in soybean (Glycine max), within the following databases: NCBI (National Center for Biotechnology Information), Scopus, Science Direct and CAPES Journals. The search was conducted from 09/2014 to 03/2015 and used the following indexing terms in various combinations: (1) Resistance; (2) Glycine; (3) Meloidogyne; (4) Gene; (5) QTL; (6) Marker; all descriptors were searched in Portuguese and English languages. Articles, published between 1990 and 2014, were found to be related to the subject. From the first articles read other potential were identified according to the methodology "snowball" (Bernard, 2011).
 
In the selected articles, the #id sequence or genes related to a phenotypic behavior of resistance to root-knot nematode were transcribed. Based on these #id or genes sought to access their DNA target sequences in their entirety using the NCBI database. Then the "FASTA" file format of each target sequence was used separately for the identification of homologous regions in the common bean genome by the "BLAST" tool (Basic Local Alignment Search Tool) available on the online program PHYTOZOME 10.3. Homologous sequences in common beans were chosen by the value E-value <1e-50 with sizes above 1000 bp. Then, the functional annotation of target sequences identified in beans was performed using the program BLAST2GO (Conesa et al., 2005).
 
To identify putative markers SSRs in the identified target sequences in common bean was used the program SSRLocator (Maia et al., 2008). Only with microsatellite repeats larger than three motifs were used for the design of the foward (F) and reverse (R) primers using the program Primer 3 (Rozen and Skaletsky, 2000). The designed primer pairs for the localized SSRs were used on an in silico PCR test model by primer- BLAST tool available in NCBI. This tool is used to check specificity of the pair primers against a database of interest, which at the time was the genome of P. vulgaris.
 
To verify the presence of dimers between the forward and reverse sequences of each primer was used AutoDimer, a program that makes use of a similar slide algorithm capable of comparing two overlapping strands of DNA (Vallone and Butler, 2004). Then the  pairs  designed EST-SSR primers were named "IFRTXX" being IFRT relative to ‘Federal Institute Goiano - Campus Urutaí’ and XX the SSR primer number (Table 1).


 RESULTS AND DISCUSSION

Seven genes described in the literature which confer resistance to root-knot nematode (Meloidogyne spp.) in Soybean (Glycine max) were found namely: Rhg1, Rhg4 (Concibido et al., 2004), Rmi1 (Luzzi et al., 1994), Extensin 1, Extensin 2, Pectin Esterase 1 and EREBP (Pham et al., 2013) (Table 1). From the sequences of the selected genes were identified regions of about 70% homology in common beans, except for the Pectin Esterase 1 gene, generating a sequence transferability rate of 87%. This success of genomics synteny between the two species is already well reported in the literature, which is due to the fact that they are phylogenetically sister groups, favoring the high occurrence of orthologous genes between them (Schmutz et al., 2014). A result that confirms the story of divergence between species from a complete duplication of the genome event approximately 56.5 million years ago (Lavin et al., 2005).
 
Using homologous target sequences in common bean was possible to recover the functional processes checked in selected genes in soybeans, with the exception of Ext1 gene, averaging 82.5% similarity and E-value values ​​ranging from 0 to 2.9 and e-168. In general, the functional processes involved five categories: Cell wall structure, organization of the cell wall, integral membrane component, protein kinase activity and transcriptional regulation (Table 2). This comparison confirms the functional completeness homology, obtained using the transferability of soybean ESTs in common beans.
 
The identification of heterologous EST through in silico approach on comparative studies of intraspecific genes evolutionarily has become an easy and low-cost strategy, not practiced before the Post-Genomic Era (Zane et al., 2002; Pandey and Sharma, 2012). This understanding has allowed advances in the development of efficient functional SSRs markers (Sterky et al., 1998; Kaliswamy et al., 2015), as this work and others (Gupta et al., 2010; Victoria et al., 2011; Wang et al., 2014). Demonstrating that ESTs are good sources of sequences for prospecting functional microsatellite markers (Andersen and Lubberstedt, 2003), in this study, and associated with resistance to root-knot nematodes in bean.
 
Twenty-six different SSRs were identified for the ESTs found in common beans from genes of genetic resistance to root-knot nematodes in soybean. Of those twenty-six only fifteen were considered because present a number of repetitions larger than three. For each of the six genes have been identified SSRs with at least two different repeat motifs, with the exception of Rhg1, Rhg4 and Rmi1 genes identified for which only one SSR (Table 1). The length of the pairs of designed primers was suitable, ranging from 18 to 26 base pairs (bp), whereas the approximate length of 20 bp defined as optimal (Wang  et al., 2014). The expected size of amplification product by PCR ranged between 182 and 357 pb, with a range of melting temperature (Tm) from 59 to 61ºC, which allows to estimate an average value approximately to 60°C among the sequences (forward and reverse) and primer pairs.
 
The small average difference in Tm between the sequences facilitates adjustment of the annealing temperature (Ta) of the F and R primers defined by "Ta = Tm ± 5", which should be similar. Defining Ta as the temperature at which pairing occurs between the primer and the flanking region of the target in the DNA of interest. The wide temperature range and length of the sequences of SSR primers >18 bp gives to the designed markers a high stringency (Bornet and Branchard, 2001). This increased reliability and reproducibility provides the certainty to be evaluating specific loci, and to facilitate the comparison of data by different researchers.
 
We expect that the fifteen EST-SSRs show a high amplification rate, since there was no competition between F and R sequences of each primer, either the formation of self-dimers. It was still possible to recover ESTs from 87% of the designed primers using the PrimerBLAST tool, allowing us to perform better thorough pre-selection, seeking to increase this amplification rate (Table 1). The success of this rate transferability studies using ESTs has been high, reported in literature values ​​from 60 to 90% in different species (Wang et al., 2014; Luro et al., 2008; Saha et al., 2004).
 
Considering the SSRs classification based on the motifs identified, was verified the occurrence almost exclusively of trinucleotide markers (n = 13) and only two (02) dinucleotide. As the number of tandem repeats found a range from four to seventeen repetitions. Of the 15 SSRs, eight were classified as perfect, five compounds and one discontinued. The perfect SSRs are those with uninterrupted sequence motif; the interrupted show a sequence of discontinuous motif for a small and different sequence; and the compounds are those formed with at least two distinct motifs (Oliveira et al., 2006).
 
The high abundance of di- and trinucleotide and small amplitude in the length of motifs in the SSR target structure are related to the sequence of EST and its mutational dynamic (Schlötterer, 2000). In expressed plant sequences, the prevalence of these repetitions has been widely reported (Gupta et al., 2010; Wang et al., 2014; Hong et al., 2015; Boccacci et al., 2015). The regions of ESTs are transcriptional targets rich in dimers in the untranslated regions (UTR) and in the predominance of trimers in the open reading frame (ORF) (Ellegren, 2004). These regions exhibit less mutational rate compared to the others repeating regions of the genome not associated with genes (Li et al., 2002), thus supporting a lower amplitude of motifs, since abrupt changes in ORFs were being selected negatively (Metzgar et al., 2002). This indicates that the reduced size is due to then action selection designed to preserve the structure of the transcript because abrupt changes easily affect the corresponding gene product.
 
 


 CONCLUSION

With this early work, fifteen (15) EST-SSRs candidate markers in six different genes that give soybean root-knot nematode resistance are available for    validation in common bean. Moreover, it was possible to demonstrate that transferability between different species, but close evolutionarily following the example soybean and common bean, can increase the efficiency in the transfer of genomic information by homology ESTs gene of interest. Finally, it is noteworthy the promising results of this herein study, once validated the genetic link by phenotypic evaluation essays of hospitability to root-knot nematode in beans and association of polymorphism amplified by suggested markers studies, they can help programs improvement, outlining their crosses directed to the phenotype of interest, in this case the resistance to root-knot nematode. Yet it is noteworthy that these markers will facilitate the development of superior genotypes with multiple sources of resistance.
 


 CONFLICT OF INTERESTS

The authors have not declared any conflict of interests.



 REFERENCES

Alzate-Marin AL, Cervigni GD, Moreira MA, Barros EG (2005). Seleção assistida por marcadores moleculares visando ao desenvolvimento de plantas resistentes a doenças, com ênfase em feijoeiro e soja. Fitopatol. Bras. 30(4):333-342.
Crossref

 

Andersen JR, Lübberstedt T (2003). Functional markers in plants. Trends Plant. Sci. 8(11):554-560.
Crossref

 
 

Bernard HR (2011). Research methods in anthropology: Qualitative and quantitative approaches. Rowman Altamira.

 
 

Boccacci P, Beltramo C, Prando MS,Lembo A, Sartor C, Mehlenbacher S, Botta R, Marinoni DT (2015). In silico mining, characterization and cross-species transferability of EST-SSR markers for European hazelnut (Corylus avellana L.). Mol. Breed. 35(1):1-14.
Crossref

 
 

Bornet B, Branchard M (2001). Nonanchored inter simple sequence repeat (ISSR) markers: reproducible and specific tools for genome fingerprinting. Plant Mol. Biol. Rep. 19(3):209-215.
Crossref

 
 

Broughton WJ, Hernandez G, Blair M, Beebe S, Gepts P, Vanderleyden J (2003). Beans (Phaseolus spp.)–model food legumes. Plant Soil 252(1):55-128.
Crossref

 
 

Burle ML, Fonseca JR, Kami JA, Gepts P (2010). Microsatellite diversity and genetic structure among common bean (Phaseolus vulgaris L.) landraces in Brazil, a secondary center of diversity. Theor. Appl. Genet. 121(5):801-813.
Crossref

 
 

Buso G,Amaral Z, Brondani R, Ferreira M (2006). Microsatellite markers for the common bean Phaseolus vulgaris. Mol. Ecol. Notes 6(1):252-254.
Crossref

 
 

Cardoso PC, Veiga MM, Menezes IPP; Valdisser PAMR, Borba TC, Melo LC, Del Peloso MJ, Brondani C, Vianello RP (2013). Molecular characterization of high performance inbred lines of Brazilian common beans. Genet. Mol. Res. 12(4):5467-5484.
Crossref

 
 

Cardoso PC, Brondani C, Menezes IPP,Valdisser PAMR, Borba TC, Del Peloso MJ,Vianello RP (2014). Discrimination of common bean cultivars using multiplexed microsatellite markers. Genet. Mol. Res. 13(1):1964-1978.
Crossref

 
 

Carneiro RG, Ferraz S, Regazzi AJ (1992). Estudo de mecanismo de resistência a Meloidogyne incognita raça 3 em variedades de feijoeiro. Nematol. Bras. 16(1/2):41-52.

 
 

Concibido VC, Diers BW, Arelli PR (2004). A decade of QTL mapping for cyst nematode resistance in soybean. Crop Sci. 44(4):1121-1131.
Crossref

 
 

Conesa A, Götz S, García-Gómez JM, Terol J, Talón M, Robles M (2005). Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics 21(18):3674-3676.
Crossref

 
 

Ellegren H (2004). Microsatellites: simple sequences with complex evolution. Nat. Rev. Genet. 5(6):435-445.
Crossref

 
 

Gupta S, Shukla R, Roy S, Sen N, Sharma A (2010). In silico SSR and FDM analysis through EST sequences in Ocimum basilicum. Plant Omics 3(4): 121-128.

 
 

Hong JH, Kwon YS, Mishra RK, Kim DH (2015). Construction of EST-SSR databases for effective cultivar identification and their applicability to complement for Lettuce (Lactuca sativa L.) Distinctness Test. Am. J. Plant Sci. 6(01):113.
Crossref

 
 

Kaliswamy P, Vellingiri S, Nathan B, Selvaraj S (2015). Microsatellite analysis in the genome of Acanthaceae: An in silico approach. Pharmacogn. Mag. 11(41):152-156.
Crossref

 
 

Kaur S, Panesar PS, Bera MB, Kaur V (2015). Simple sequence repeat markers in genetic divergence and marker-assisted selection of rice cultivars: A review. Crit. Rev. Food. Sci. Nutr. 55(1):41-49.
Crossref

 
 

Lavin M, Herendeen PS, Wojciechowski MF (2005). Evolutionary rates analysis of Leguminosae implicates a rapid diversification of lineages during the Tertiary. Syst. Biol. 54(4):575-594.
Crossref

 
 

Li OJ, Chen XM, Xia PX, Pei ZY, Wang Y, Lan QK, Zhang RW (2015a). EST-SSR marker-based assay for purity identification of melon "Green Angle". Adv. Appl. Biotechnol. 333:637-642.
Crossref

 
 

Li X, Hu X, Hu T, Li G, Ru Z, Zhang L, Lang Y (2015b). Identification of a novel wheat-Thinopyrum ponticumaddition line revealed with cytology, SSR, EST-SSR, EST-STS and PLUG markers. Cereal Res. Commun. 43(4):1-10.
Crossref

 
 

Li YC, Korol AB, Fahima T, Beiles A, Nevo E (2002). Microsatellites: genomic distribution, putative functions and mutational mechanisms: a review. Mol. Ecol. 11(12):2453-2465.
Crossref

 
 

Luro FL, Costantino G, Terol J, Argou X, Allario T, Wincker P, Talon M, Ollitrault P, Morillon R (2008). Transferability of the EST-SSRs developed on Nules clementine (Citrus clementina Hort ex Tan) to other Citrus species and their effectiveness for genetic mapping. BMC Genome 9:1-13.
Crossref

 
 

Luzzi B, Boerma H, Hussey R (1994). A gene for resistance to the southern root-knot nematode in soybean. J. Hered. 85(6):484-486.

 
 

Maia LC, Palmieri DA, Souza VQ, Kopp MM, Carvalho FIF, Oliveira AC (2008). SSR Locator: Tool for simple sequence repeat discovery integrated with primer design and PCR simulation. Int. J. Plant. Genome 1-9
Crossref

 
 

Metzgar D, Liu L, Hansen C, Dybvig K, Wills C (2002). Domain-level differences in microsatellite distribution and content result from different relative rates of insertion and deletion mutations. Genome Res.12(3):408-413.
Crossref

 
 

Müller BS, Pappas GJ Jr, Valdisser PA, Coelho GR, Menezes IP, Abreu AG, Borba TC, Sakamoto T, Brondani C, Barros EG (2015). An Operational SNP Panel Integrated to SSR Marker for the Assessment of Genetic Diversity and Population Structure of the Common Bean. Plant Mol. Biol. Rep. 33(6):1697-1711.
Crossref

 
 

Müller BS, Sakamoto T, Menezes IPP, Prado GS, Martins WS, Brondani C, Barros EG, Vianello RP (2014). Analysis of BAC-end sequences in common bean (Phaseolus vulgaris L.) towards the development and characterization of long motifs SSRs. Plant Mol. Biol. 86(45):455-470.
Crossref

 
 

Oliveira EJ, Pádua JG, Zucchi MI, Vencovsky R, Vieira MLC (2006). Origin, evolution and genome distribution of microsatellites. Genet. Mol. Biol. 29(2):294-307.
Crossref

 
 

Pandey M, Sharma J (2012). Efficiency of microsatellite isolation from orchids via next generation sequencing.

 
 

Pham AT, Mcnally K, Abdel-Haleem H, Boerma HR, Li Z (2013). Fine mapping and identification of candidate genes controlling the resistance to southern root-knot nematode in PI 96354. Theor. Appl. Genet. 126(7):1825-1838.
Crossref

 
 

Rozen S, Skaletsky H (2000). Primer3 on the WWW for general users and for biologist programmers. Meth. Mol. Biol. 132:365-386.

 
 

Saha MC, Mian MR, Eujayl I, Zwonitzer JC, Wang L, May GD (2004). Tall fescue EST-SSR markers with transferability across several grass species. Theor. Appl. Genet. 109(4):783-791.
Crossref

 
 

Schlötterer C (2000). Evolutionary dynamics of microsatellite DNA. Chromosoma 109(6):365-371.
Crossref

 
 

Schmutz J, Mcclean PE, Mamidi S, Wu GA, Cannon SB, Grimwood J, Jenkins J, Shu S, Song Q, Chavarro C (2014). A reference genome for common bean and genome-wide analysis of dual domestications. Nat. Genet. 46(7):707-713.
Crossref

 
 

Silva OF, Wander AE (2013). O feijão-comum no Brasil passado, presente e futuro.

 
 

Sterky F, Regan S, Karlsson J, Hertzberg M, Rohde A, Holmberg A, Amini B, Bhalerao R, Larsson M, Villarroel R (1998). Gene discovery in the wood-forming tissues of poplar: analysis of 5,692 expressed sequence tags. Proc. Natl. Acad. Sci. 95(22):13330-13335.
Crossref

 
 

Vallone PM, Butler JM (2004). AutoDimer: a screening tool for primer-dimer and hairpin structures. Biotechniques 37(2): 226-231.

 
 

Victoria FC, Maia LC, Oliveira AC (2011). In silico comparative analysis of SSR markers in plants. BMC Plant Biol. 11(1):15.
Crossref

 
 

Vieira C, Paula Júnior T, Borém A (1998). Adubação mineral e calagem. Feijão 2:115-142.

 
 

Walber R, Juliatti F, Santos M (2003). Avaliação de acessos de feijoeiro em relação aos nematóides das galhas. In: Congresso Brasileiro de Fitopatologia. Soc. Bras. Fitopatol. 293-294.

 
 

Wander AE (2005). Perspectivas de mercado interno e externo para o feijão. In: Congresso Nacional De Pesquisa De Feijão. pp. 892-895.

 
 

Wang B, Zhu P, Yuan Y, Wang C, Yu C, Zhang H, Zhu X, Wang W, Yao C, Zhuang Z (2014). Development of EST-SSR markers related to salt tolerance and their application in genetic diversity and evolution analysis in Gossypium. Genet. Mol. Res. 13(2):3732-3746.
Crossref

 
 

Yokoyama L (2007). Cultivo do feijoeiro comum: importância econômica. Sistemas de Produção. 2.

 
 

Zane L, Bargelloni L, Patarnello T (2002). Strategies for microsatellite isolation: a review. Mol. Ecol. 11(1):1-16.
Crossref

 
 

Zhai C, Xu P, Zhang X, Guo Q, Zhang X, Xu Z, Shen X (2015). Development of Gossypium anomalum-derived microsatellite markers and their use for genome-wide identification of recombination between the G. anomalum and G. hirsutum genomes. Theor. Appl. Genet. 128(8):1531-1540.
Crossref

 

 




          */?>