SciELO - Scientific Electronic Library Online

 
vol.5 número13Diversidad genética en cultivares de manzano establecidos en la Sierra de Arteaga Coahuila, MéxicoDosis letal media y reducción media del crecimiento por radiación gamma en pasto africano (Eragrostis lehmanniana Ness) índice de autoresíndice de assuntospesquisa de artigos
Home Pagelista alfabética de periódicos  

Serviços Personalizados

Journal

Artigo

Indicadores

Links relacionados

  • Não possue artigos similaresSimilares em SciELO

Compartilhar


Ecosistemas y recursos agropecuarios

versão On-line ISSN 2007-901Xversão impressa ISSN 2007-9028

Ecosistemas y recur. agropecuarios vol.5 no.13 Villahermosa Jan./Abr. 2018

https://doi.org/10.19136/era.a5nl3.1180 

Notas científicas

Differentiation of apomictic and sexual genotypes of Brachiaria spp., using molecularmarkers

Diferenciación de genotipos apomícticos y sexuales de Brachiaria spp., con marcadores Moleculares

Jaime Poblete-Vargas1 

Ernestina Valadez-Moctezuma2  * 

Gabino García de los Santos1 

Carlos Martínez-Flores3 

Armando Peralta-Martínez4 

1Programa de Semillas, Colegio de Postgraduados Campus Montecillo. Km 36.5 Carretera México- Texcoco, CP. 56230, Montecillo, Texcoco, Estado de México, México.

2Laboratorio de Biología Molecular, Departamento de Fitotecnia. Universidad Autónoma Chapingo. Carretera México-Texcoco Km 38.5 Chapingo, CP. 56230, Estado de México, México.

3Programa Global de Maíz. Centro Internacional de Mejoramiento de Maíz y Trigo. Carretera México-Veracruz. Km. 45, Colonia El Batán, CP. 56237, Texcoco, Estado de México, México.

4Agroproductos de Iguala, S.A. de C.V. Hidalgo No. 23, Col. Centro, CP. 40400, Iguala, Guerrero, México.


ABSTRACT.

Some species of Brachiaria are cultivated in the tropics because of their high productivity and drought resistance; their apomictic hybrids are of interest because of the almost nuil segregation and fixation of the hybrid vigor in seeds. In this research, Amplified Fragment Length Polymorphism (AFLP) and Cleaved Amplified Polymorphic Sequence (CAPS) markers were used to differentiate apomictic tetraploid and sexual diploid Brachiaria parents and their progeny. Polymorphism detected with AFLP was 91.3% and with CAPS 79.7%. Phenograms differentiated each genotype and the clustering of parents and their progeny was by their degree of genetic relatedness. AFLP did not separate the apomictic genotypes from the sexual ones, but CAPS markers did so through the C15-8 marker. Some apomictic genotypes of the progeny inherited it from their B. ruziziensis sexual mother, which may be related to some segregating character of the parental line.

Key words: AFLP; CAPS; DNA; forage; polymorphism

RESUMEN.

Algunas especies de Brachiaria se cultivan en los trópicos por su alta productividad y resistencia a sequía; sus híbridos apomicticos, son de interés por la casi nula segregación y fijación del vigor híbrido. En este trabajo se utilizaron marcadores Amplified Fragment Length Polymorphism (AFLPs) y Cleaved Amplified Polymorphic Sequences (CAPS) para diferenciar progenitores apomicticos tetraploides y sexuales diploides de Brachiaria y su progenie. El polimorfismo detectado con AFLPs fue del 91.3% y con CAPS del 79.7%. Los fenogramas diferenciaron a cada genotipo y agrupación de progenitores con la progenie por su grado de parentesco genético. Los AFLPs no separaron a los genotipos apomicticos de los sexuales, pero los CAPS si lo lograron mediante el marcador C15-8; algunos genotipos apomicticos de la progenie, lo heredaron de la madre sexual B. ruziziensis, lo que puede estar relacionado con algún carácter segregante de la línea progenitora.

INTRODUCTION

Livestock production is a key sector in the economies of Latin American countries, occupying a large area of land with productive potential (FAO 2009). Among the forage plants used, the genus Brachiaria, which is characterized by having cultivated forage pastures, stands out (Do Valle and Miles et al. 2001). The success of these pastures is due to the fact that many species have a wide natural variation in agronomic traits, such as biomass production, nutrient quality, drought tolerance, tolerance to flooded soils and high seed yields (Miles et al. 2004). The qualities of these forage species have been the reason for developing hybrids with resistance to cuckoo spit, the frothed-up sap caused by the spittlebug (Aeneolamia postica), which limits the passage of water and nutrients in the plant (De la Cruz-Llanas et al. 2005), and to foliar fungi such as Rhizotocnia; hybrids have also been developed that are adapted to acid soils with high leveis of aluminum (Rao et al. 2006) and that have higher production of high digestibility forage (Lascano 2002). Most commercially cultivated Brachiaria species are tetraploid apomicts; only B. ruzizienzis reproduces sexually and its recombination with apomictic genotypes can generate hybrids with desirable characteristics (Miles et al. 2004).

Apomixis is present in some Brachiaria species, which implies that the progeny are genetically identical to the mother plant, so they do not present segregation, a phenomenon that is often of importance in plant breeding (Grimanelli et al. 2001). This character is considered relevant in agriculture, due to the possibility of transforming crops of worldwide importance from sexual reproduction to apomictic clones (Bicknell and Koltunow 2004). Apomixis has great agricultural interest, as its fixa-tion would allow maintaining the genetic stability of hybrids and varieties for an indefinite period, while the progeny of sexual reproduction maintain their genetic variation (Koltunow et al. 1995). B. decumbens, B. brizantha and B. ruziziensis are considered the most important species in this regard, since they belong to the same agomictic complex, with the first two being apomictic and the last one sexual (Do Valle and Miles 2001, Risso-Pascotto et al. 2005).

Traditionally in Brachiaria, morphological descriptors and embryo sac analysis with pistil thinning techniques are used to differentiate progenyvor sexual (Savidan 2000), while genetic-molecular marker techniques have been used successfully since the 1990s, with the advantage that they provide results in a short time from any plant tissue and any plant age; they are efficient because they are not affected by the environment. In most cases, the markers complement the information from the morphological markers (Azofeifa-Delgado 2006). In particular, the Amplified Fragment Length Polymorphism (AFLP) and Cleaved Amplified Polymorphic Sequence (CAPS) techniques are used to sample the entire genome (Zorzatto et al. 2010). There-fore, the objective of the present study was to apply the AFLP and CAPS techniques to differentiate interspecific and sexual apomictic genotypes of Brachiaria spp, selected as apomictic in the field.

MATERIALS AND METHODS

Plant material

We used young seedling leaves from 30 hybrid genotypes from three artificial hybrid tetraploid (2n + 2n= 4n = 36) male (♂) apomictic Brachiaria briza ntha x B. ruziziensis parents; five tetraploid female (♀) Braploiaria ruziziensis (4x = 36) parents and 22 hybrids of the progeny resulting from 10 crosses selected in the field with the apomixis character (Table 1). The parent ♂ 18 (BR/NO 1873) is registered in México as cultivar Mulato I, characterized by having high genetic uniformity; parent ♂71 (BR/NO 1371) is facultative, only produces viable pollen and is identified as self-compatible (CIAT 2002); it was used as a pollinator because it is a cuckoo spit-resistant apomict. Parent ♂ 94 (AIG 2094) is a seed-producing, drought-resistant hybrid. The progeny tests for identifying the reproductive mode of the mother plants and the 22 apomictic hybrids were carried out in an experimental lot in Tuxpan, Iguala county, Guerrero state, México, located at 18° 21' NL, 99° 29' WL.

DNA extraction

From 1 g of lyophilized tissue obtained from the mixture of the genotypes listed in Table 1, the DNA was obtained with the cetyltrimethylammonium bromide (CTAB) method (Dellaporta et al. 1983, Zhang and Stewart 2000).

Table 1 Genetic material and keys used in molecular characterization. 

*ln the hybrids (9-30) the number after the hyphen corresponds to the number of plants resulting from the cross with apomictic characteristics considered in the study. ** In the hybrids (9-30) the numbers correspond to the parents involved in the crossing.

Molecular evaluation

DNA quantification was performed with spec-trophotometry at 260 nm and its quality was esti-mated on 1% agarose gels with IX TAE buffer (Sambrook et al. 1989).

AFLP

Digestión, adapter ligation, pre-amplification reactions and selective amplification were per formed as indicated by GIBCO-BRL Life Technolo gies AFLP™. For the selective amplification, the following primer combinations were evaluated: E-AAC + M-CAC, E-AAC + M-CTG, E-ACC + M-CAC, E-ACC + M-CTG, E-AGG + M-CTC, E-AGG + M-CTA, E-ACT + M-CTC y E-ACT + M-CTA and those which produced clear differences between parents and progeny were selected (Table 2). The fragments obtained in the selective amplification were separated on denaturing 6% polyacrylamide gels (19:1), 7.5 M urea and IX TBE buffer in a BioRAD® vertical sequencing chamber and an E-C Apparatus Corporati ® model 3000PW power source. The reagents required for electrophoresis and staining were prepared according to Sambrook et al. (1989). Prior to running the polymerase chain reaction (PCR) products, the gel was pre-run for 30 min at 55 °C. Subsequently, 10 µL of the samples previously mixed and heated with a run buffer at Electrophoresis was run at 80 watts (1000 volts) for 3 h; the gel was then stained in 0.2% silver nitrate solution for 30 min and developed with sodium carbonate solution, made with 30 g of sodium carbonate (Na2CO3) per L of distilled water, which formaldehyde and 400 µL of sodium thiosulfate (10 mg mL-1) were added. When the first bands were developed, after approximately 4 min, the solution was discarded and replaced with fresh solution to continue the development for another 15 min. Finally, the gels were photodocumented for analysis.

Table 2 Primers, number of fragments and level of polymorphism detected with the Amplified Fragment Length Polymorphism (AFLP) and Cleaved Amplified Polymorphic Sequences (CAPS) techniques in parents and progeny of the genus Brachiaria

CAPS

It was performed with a modification, which consisted in digesting the DNA with the Eco Rl enzyme prior to PCR, while in the original technique the amplified fragments are digested (Konieczny and Ausubel 1993). The modification provided reliable, informative and reproducible results between experiments, because the conformational change in the molecule used in PCR (fragments) exhibits a less complicated topography. PCR amplifications were perfomed on an Applied Biosystems® 9700 thermocycler according to the methodology of Williams et al. (1990). The selected RAPD primers were C-06, C-08, C-15 and C-19, of 25 primers previously evaluated by the ROTH Company, Germany (Table 2). The thermocycling consisted of one minute at 94o C for pre-denaturation, 38 cycles [94oC, 20 s; 40oC, 15 sec; 72oC, 60 s] and a 5-min final extension cycle at 72 oC. The products were separated in acrylamide (29: 1) double-layer gels where the run-off gel was prepared at 16%. The run buffer was IX Tris-Glycine and it was run for 4 h with 40 V/cm in a C. B. S. Scientific CO® model MVG-216-33 vertical chamber, with an EC-105 source. The gels were stained with 0.2% silver nitrate (AgNO3), visualized and documented with a white light transilluminator (Sambrook et al. 1989).

Statistical analysis

A binary matrix was constructed with the products amplified in each technique. Genetic similarities were calculated with Nei and Li/Dice distance (Nei and Li 1979) and the method for clustering was the Unweighted Pair Group Method with Arithmetic mean (UPGMA) with a 500-replicate bootstrap analysis. The program used was Free-Tree v 0.9.1.50 and the trees were drawn with the TreeView program (Win32) 1.6.6. A simple correspondence factor analysis was also carried out with the molecular data obtained with both tech niques, to determine the percent contribution of the markers to the differentiation of the genotypes with respect to apomictic or sexual character.

RESULTS AND DISCUSSION

The number of markers obtained was acceptable and informative; 147 corresponded to AFLP, of which 91.3% were polymorphic, while 87 were obtained with CAPS with 79.7% polymorphism. Both types of markers differentiated each of the genotypes studied (Table 2). In the consensus phenogram, constructed with the AFLP markers (Figure 1), three main groups are distinguished: group A comprised by the hybrid 143A-18-3; group B defined by the ♂ 94 and distant from the other two males; and group H with genotypes 4171-9, 11094-1 and 11094-16. Group B was related to the rest of the genotypes, which were clustered based on close genetic similarities. Groups C and G show groupings of full siblings. The C group concentrates genotypes with greater influence of ♂ 18 [10018-1, 10018-16 and 10018-6] and [4118-12, 4118-10 and 4118-6], which is a dominant apomict, and ♀ 110, 29 and 41; this group also includes the genotype 4171-20, which must share characters with the parent ♀ 41 and the indicated half-siblings. Group G contains full siblings originating from ♀ 143A and ♂ 94 [143A94-11, 143A94-3 and 143A94-4, in addition to genotype 143A18-23, which apparently with ♂18. On the other hand, the genotype 143A18-3 is found as an independent group (group A) when it was expected to be grouped with its genetic similarities. The discrepancy in this behavior may be due to mutations in the genome that generate polymorphisms in the profiles obtained. These alterations are common in plants of asexual/vegetative reproduction and it is relatively easy to detect them with molecular markers. Group D contains ♂ 71 and some of its progeny, group E has ♂ 18, and finally group F constains ♀ 100 and 143A. The fact that all the females were distributed throughout the phenogram suggests the genetic contribution of B. ruziziensis used as a female parent. The AFLP markers were able to differentiate all the genotypes and to show the close genomic relationship between some of them; however, they were insufficient to mark the locus or loci responsible for the apomixis (Hand and Koltunow 2014).

Figure 1 Consensus phenogram of apomictic and sexual genotypes of Brachiaria spp., constructed with Amplified Fragments Length Polymorphism (AFLP) genetic markers. Un-weighted Pair Group Method with Arithmetic Mean (UPGMA) clustering and random sam-pling with 500 replicates (bootstrapping). 

Figure 2 shows the grouping formed with CAPS data; ♂ 18 remained independent of the two main groupings, group I and group II, suggesting less genetic similarity according to these markers. Group I includes the parent ♂ 94, the genotypes with greater genetic influence of this parent, as well as some genotypes related to the parents ♂ 18 and ♂71. Group II comprises the females 29, 110 and 41, as well as the hybrids that shared more characters with them, but also ♂ 71, probably due to its facultative character. The other two females, 100 and 143A, were grouped independently. As in the phenogram with AFLP data, the groupings with CAPS were also made up of subgroups of full siblings; forexample, group I included 4: [10018-16,10018-6]; [143A94-3, 143A94-4]; [11094-16, 110941, 11094-7] and [2918-18, 2818-14], while group II comprised 2 subgroups, [4118-10 and 4118-12] and [4194-19, 4194-24], mainly. There are several matches between the two phenograms; there were groupings of full siblings, as well as females 100 and 143A or 29 and 41, which stayed together but in sepárate groups, which explains their close ge netic resemblance. On the other hand, in both phenograms it was possible to sepárate the females from the males; however, in the case of CAPS, ♂71 was not separated from the group of females 29, 110 and 41, perhaps due to their particular qualities. The distribution behavior observed in both phenograms can be useful, especially when considering the alternative of selecting any of the females used in the present study due to their degree of genetic similarity, or ♂18 and 94 that showed less relatedness between them in both marker systems. Another observation was that the progeny was mainly influenced by the genetic characters of the parents ♂18 and ♂94, and to a lesser degree by ♂71.

Figure 2 Consensus phenogram of apomictic and sexual genotypes of Brachiaria spp., constructed with Cleaved Amplified Polymorphic Sequences (CAPS) genetic markers, Un-weighted Pair Group Method with Arithmetic Mean (UPGMA) clustering and random sam-pling with 500 replicates (bootstrapping). 

The search for characters associated with apomixis or sexuality with techniques that detect anonymous markers such as those used in this study has proved difficult, unless sufficient molecular data are considered; this probability increases if precise information is available for the two loci that are assumed to control apomixis in plants (Noyes and Rieseberg 2000). DNA markers associated with apomixis (AFLPs and SCARs) have been re-ported for the genera Paspalum (Labombarda et al. 2002), Pennisetum (Ozias-Akins et al. 1998), Hieriacium (Catanach et al. 2006) and Trichloris crinite (Cavagnaro et al. 2006). In this study, the CAPS markers detected a DNA fragment called C15-8 that differentiated the sexual parents (♀) from the apomictic ones (♂) and that could be associated with the sexuality of B. ruziziensis; this may indicate the absence of this distinctive character in the apomictic genotypes used as males. This fragment was inherited from hybrids such as 9, 10, 12, 14, 15, 16, 18, 19, 22, 26 and 27 (Figure 3) that were selected at the field level as apomicts, one (♂) and that could be associated which may indicate that the B. ruziziensis line used as female is not the same as that of the tetraploid males. This observation can also be supported by separating the females in the two phenograms obtained. An alternative way of determining the effect of this marker is to observe the behavior at field level with respect to sexuality in different environments, to confirm it in a traditional way. The informative capacity of this CAPS technique has also been shown by authors such as Móhring et al. (2005) who detected self-incompatibility characters in Brassica species. The relative ease of differentiating sexual genotypes (Brachiaria ruziziensis) from interspecific apomicts (B. brizantha and B. ruziziensis) in the present study could be favored by the genomic contrast of the species types involved in the crosses.

Figure 3 Cleaved Amplified Polymorphic Sequences (CAPS) profiles of 30 genotypes of Brachiaria spp., obtained with the ROTH C-15 primer. Lanes 1-3 correspond to 4X males, 4-8 to 2X females and 9-30 to hybrids. The ellipses indicate a distinctive DNA fragment in the females (lanes 4-8) that were inherited by some hybrids and named C15-8. Lane M indicates the molecular weight marker øX174 DNA /Haelll. 

The simple correspondence factor analyses (Table 3) performed for the AFLP and CAPS data allowed weighting the value contributed by each component (amplicon) (Dore and Ojasso 2001). With both techniques, the first six principal components explain 56.47% for AFLP and 54.27% for CAPS, which suggests the effectiveness of the markers used to differentiate each of the genomes. The AFLP and CAPS techniques were appropriate to estimate the genetic relationship of parents and hybrids of the genus Brachiaria. The genetic contribution of male parents 18 and 94 was shown in the progeny, compared to that provided by parent ♂71. The CAPS technique was able to clearly differentiate the diploid genotypes of B. ruziziensis (female) from the tetraploid genotypes of B. brizantha X B. ruziziensis used as males. The CAPS factorial analysis allowed correlating the DNA fragment called C15-8 as a distinctive character of the female sexual parents, unlike that of the male apomictic parents.

Table 3 Eigenvalues of genotypes of the genus Brachiaria, product of simple correspondence 

LITERATURE CITED

Azofeita-Delgado A (2006) Uso de marcadores moleculares en plantas; aplicación en frutales del trópico. Agronomía Mesoamericana 17: 221-242. [ Links ]

Bicknell RA, Koltunow AM (2004) Understanding apomixis: recent advances and remaining conundrums. The Plant Cell 16: S228-S245. [ Links ]

Catanach AS, Erasmuson SK, Podivinsky E, Jordán BR, Bicknell R (2006) Deletion mapping of genetic regions associated with apomixis in Hieracium. Proceedings of the National Academy of Sciences of the United States of America 103: 18650-18655. [ Links ]

Cavagnaro PF, Cavagnaro JB, Lemes JL, Masuell RW, Passera CB (2006) Genetic diversity among varieties of the native forage grass Trichloris crinite based on AFLP markers, morphological characters, and quantitative agronomic traits. Genome 49: 906-918. [ Links ]

CIAT (2002) Annual Report: Part 1: Grass and legume genotypes with high forage attributes. CIAT. Cali, Colombia. 23p. [ Links ]

De la Cruz-Llanas JJ, Vera-Graziano J, López-Collado J, Pinto MV, Garza-García R (2005) Una técnica simple para el desarrollo de ninfas de Aeneolamia postica (Homoptera: Cercopidae). Folia Entomológica Mexicana 44: 91-93. [ Links ]

Dellaporta SL, Wood J, Hicks JB (1983) A plant DNA minipreparathion: Versión II. Plant Molecular Biology Reporter 1: 19-21. [ Links ]

Do Valle CB, Miles JW (2001) Breeding of apomictic species. In: CYMMYT (Ed.). The flowering of apomixis: from mechanisms to genetic engineering. CYMMYT/IRD. México, pp: 137-152. [ Links ]

Dore JC, Ojasoo T (2001) How to analyze publication time trends by correspondence factor analysis: Analysis of publications by 48 countries in 18 disciplines over 12 years. Journal of the Association for Information Science and Technology 52: 763-769. [ Links ]

FAO (2009) El estado mundial de la agricultura y la alimentación. Organización de las Naciones Unidas para la Agricultura y la Alimentación. Roma, Italia. 184p. [ Links ]

Grimanelli D, Leblanc O, Perotti E, Grossniklaus U (2001) Developmental genetics of gametophytic apomixis. Trends in Genetics 17: 597-604. [ Links ]

Hand ML, Koltunow AM (2014) The genetic control of apomixis: Asexual seed formation. Genetics 197: 441-450. [ Links ]

Koltunow AM, Bicknell RA, Chaudhury AM (1995) Apomixis: Molecular strategies for the generation of genetically identical seeds without fertilizaron. Plant Physiology 108: 1345-1352. [ Links ]

Konieczny A, Ausubel FM (1993) A procedure for mapping Arabidopsis mutations using co-dominant ecotype-specific PCR-based markers. The Plant Journal 4: 403-410. [ Links ]

Labombarda P, Busti A, Caceres MA, Pupilli F, Arcioni S (2002) An AFLP marker tightly linked to apomixis reveáis hemizygosity in a portion of the apomixis-controlling locus in Paspalum simplex. Genome 45: 513-519. [ Links ]

Lascano CE (2002) Caracterización de las pasturas para maximizar producción animal. Archivos Latinoamericanos de Producción Animal 10: 126-132. [ Links ]

Miles JW, Do Valle CB, Rao IM, Euclides VP (2004) Brachiaria grasses. In: Sollenberger LE, Moser LE, Burson BL (eds) Warm-season (C4) grasses. Agronomy monograph 45, American Society of Agronomy, Crop Science Society of America, Soil Science Society of America. Madison, Wl, USA. pp: 745-783. [ Links ]

Möhring S, Horstmann V, Esch E (2005) Development of a molecular CAPS marker for the self-incompatibility locus in Brassica napus and identificaron of different S alleles. Plant Breeding 124: 105-110. [ Links ]

Nei M, Li WH (1979) Mathematical model for studying genetic variation in terms of restriction endonucleases. Proceedings of the National Academy of Sciences United States of America 76: 5269-5273. [ Links ]

Noyes RD, Rieseberg LH (2000) Two independent loci control agamospermy (apomixis) in the triploid flowering plant Erigeron annuus. Genetics 155: 379-390. [ Links ]

Ozias-Akins P, Roche D, Hanna WW (1998) Tight clustering and hemizygosity of apomixis-linked molecular markers in Pennisetum squamulatum implies genetic control of apospory by a divergent locus that may have no allelic form in sexual genotypes. Proceedings of the National Academy of Sciences of the United States of America 95: 5127-5132. [ Links ]

Rao IM, Miles JW, García R, Ricaurte J (2006) Selección de híbridos de Brachiaria con resistencia a aluminio. Pasturas Tropicales 28: 12-15. [ Links ]

Risso-Pascotto C, Pagliarini MS, Valle CB (2005) Meiotic behavior in interspecific hybrids between Brachiaria ruziziensis and Brachiaria brizantha (Poaceae). Euphytica 145: 155-159. [ Links ]

Sambrook J, Fritsch EF, Maniatis T (1989) Molecular cloning: a laboratory manual. 2a. ed. Cold Spring Harbor Laboratory Press, NY, USA. pp: 18.51-18.57. [ Links ]

Savidan Y (2000) Apomixis: Genetics and breeding. In: Janick J (Ed.). Plant Breeding Reviews Volume 18. Wiley. New York, USA. pp: 13-86. [ Links ]

Williams JGK, Kubelik AR, Livak KJ, Rafalski JA, Tingey SV (1990) DNA polymorphisms amplified by arbitrary primers are useful as genetic markers. Nucleic Acids Research 18: 6531-6535. [ Links ]

Zhang J, Stewart JMD (2000) Economical and rapid method for extracting cotton genomic DNA. The Journal of Cotton Science 4: 193-201. [ Links ]

Zorzatto C, Chiari L, Araújo-Bitencourt G, Do Valle CB, De Campos-Leguizamón GO, Schuster I, et al. (2010) Identification of a molecular marker linked to apomixis in Brachiaria humidicola (Poaceae). Plant Breeding 129: 734-736. [ Links ]

Received: August 03, 2016; Accepted: May 11, 2017

*Corresponding author: evaladez@chapingo.mx

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License