Open-access Mating system analysis of Açaí-do-Amazonas (Euterpe precatoria Mart.) using molecular markers

Abstract

Euterpe precatoria (Açaí-do-Amazonas) produces fruits of which the fresh pulp is consumed. It is almost exclusively collected by extractivist farmers, because no selected genotypes are available for the establishment of plantations. For the domestication and breeding of the species, mating system studies are needed for strategy formulation. This study evaluated the mating system of a natural population of E. precatoria. Thirteen progenies were genotyped with 13 microsatellite loci by capillary electrophoresis in an automated DNA sequencer. Estimates of single-locus and multilocus outcrossing rates were 1.0, and paternity correlation was low (r^p(m)>= 0.293). Euterpe precatoria families consist mainly of half-sibs and the reproductive strategy of the species is allogamy.

Keywords: Arecaceae; Euterpe precatoria; molecular marker.

INTRODUCTION

The palm Euterpe precatoria Mart. (acaí-do-Amazonas) belongs to the family Arecaceae, occurring in the western and central Brazilian Amazon and to within borders of the Amazon of Peru, Brazil, Colombia (Kahn 1991), and Bolivia (Bussmann and Zambrana 2012). The species is exploited by an extractive production chain, maintained by the local population since several decades and currently intensified due to the interest of different industries. Fruit for consumption in the form of beverage is the main product of E. precatoria (Noda 2012).

When correlating the total production of açaí fruits exploited in the different states of Brazil (IBGE 2018) with the information of geographic distribution and habitat (Lorenzi et al. 2010), it becomes clear that açaí from the states of Amazonas, Acre and Rondônia is extracted from fruits of E. precatoria. However, in the states of Pará, Maranhão and Amapá, it is extracted from Euterpe oleracea, another açaí species. Euterpe precatoria is single-stemmed, while E. oleracea has tillers (multi-stemmed). One of the advantages of E. precatoria is that the fruit pulp has better anti-oxidant and anti-inflammatory properties than E. oleracea (Kang et al. 2012). Moreover, the evolutionary adaptation to nutrient-poor and well-drained lowland, Latosolos and Argisolos is noteworthy (FAO 1987). The adaptation of E. precatoria to the Amazonian ecosystem indicates it for exploitation in the recovery of degraded areas, although no seeds of improved varieties recommended for the production of seedlings and planting are available, since to date, E. precatoria has been neglected in research.

The inflorescences of E. precatoria have numerous male (4.5 x 2.7 mm) and female (3.2 x 2.6 mm) flowers. The male flowers open and release pollen before the female flowers are receptive (Küchmeister et al. 1997, Lorenzi et al. 2010). However, this observation of protandry is no guarantee that the species is autogamous, since nonetheless some crosses may occur. Thus, research must be carried out to fill gaps in the knowledge about the gene flow and mating system to interpret its behavior within tropical forest ecosystems in in situ or ex situ conservation areas.

Euterpe precatoria is a highly important species for the Amazon region. As the genetic structure of tree populations is partially determined by the mating system and strongly by the gene flow among populations (Ramos et al. 2016a), this study will provide valuable information for the domestication and rational management of the species (Ramos et al. 2011). Studies of the mating system can be based on microsatellite markers or simple sequence repeats (SSRs), which are an adequate tool for this purpose, due to the high polymorphism in terms of number of alleles, co-dominant inheritance and low cost of the method (Wadt et al. 2015). Simple sequence repeat-based studies to determine the mating system were carried out for different Amazonian and tropical tree species (Ramos et al. 2011, Medina-Macedo et al. 2015, Picanço-Rodrigues et al. 2015; Moraes et al. 2018).

The objective of this study was to investigate the mating system of a natural E. precatoria population, to identify the levels of outcrossing, mating among relatives and correlated matings to better understand the genetic structure of open-pollinated progenies for ex situ conservation and breeding plans. Moreover, we tried to determine the coancestry coefficient, effective population size within progeny and the estimated number of trees to ensure enough seeds for conservation.

MATERIAL AND METHODS

Study area and sampling

The study investigated a natural E. precatoria population in the community of Nossa Senhora de Fatima do Açaí, in the rural area of Vila Amazônia (lat 02º 36' 52.09" S, long 56º 33' 29.13" W, and alt around 40 m), in Parintins, state of Amazonas, Brazil, where the climate is tropical monsoon (Am) (Peel et al. 2007). All 13 trees bearing fruit within 10 hectares of climax forest were used for the study. Plant material (leaflets) and 100 fruits per mother tree were collected, i.e., a total sample of 1300 fruits. The fruits were transported to the seed laboratory of the experimental field Caldeirão of Embrapa Western Amazon, located at Rodovia Manoel Urbano, Km 13, Estrada do Caldeirão, Iranduba-AM, for seed germination, emergence and seedling production. To accelerate germination, fruits with pulp and seeds were immersed in water at 60 °C and removed when the water reached room temperature (Nogueira et al. 1995). After three months, seedlings of all progenies with at least two bifid leaves were sampled for DNA extraction and microsatellite analysis, resulting in 3 to 25 seedlings per progeny, totaling 227 seedlings (Table 2). One leaflet of each seedling (227 samples) and mother tree (13 samples) was collected and packed separately in a previously labelled zip lock plastic bag containing silica gel. These samples were transported to the Laboratory of Genetics and Plant Breeding of the Department of Agronomy of the Federal University of Amazonas (UFAM) for storage at -20 °C for further genomic DNA isolation.

Table 2
Mating system at the mother tree level

Amplification of microsatellite markers

The DNA was extracted based on the Cationic Hexadecyl Trimethyl Ammonium Bromide (CTAB) method described by Doyle and Doyle (1990) and then quantified (Ramos et al. 2011). A total of 13 microsatellite primers (Epr01, Epr02, Epr04, Epr05, Epr13, Epr14, Epr15, Epr15, Epr18, Epr19, Epr21, Epr22, Epr31, and Epr32) were used in the study (Ramos et al. 2016b). These microsatellite loci were amplified by the polymerase chain reaction (PCR) with a Veriti thermal cycler (Applied Biosystems) in a total volume of 10 μL per reaction (containing 10 ng of genomic DNA, 1X buffer, 210 μM of each dNTP, 1.5mM MgCl2, 0.16 µM forward primer and M13 label primers (FAM or NED) broth (Schuelke 2000), 0.32 μM reverse primer, 1.05 U Taq DNA polymerase (Invitrogen, Carlsbad, California, USA), and 3.49 μl ultra-pure water. These PCR amplifications consisted of two steps, of which the first was primer-specific and the second with M13 binding (Ramos et al. 2016b). The amplification products were checked by electrophoresis on 1.5% agarose gels stained with GelRed (Biotium) in 1x TBE buffer (pH 8.0). The amplified PCR products were subjected to an automated DNA analyzer by capillary electrophoresis (ABI 3130XL, Genetic Analyzer, Applied Biosystems). A standard size GeneScan™-500 ROX® (Life Technologies of Brazil Ltda.) was used to determine the allele size. The amplified fragments were observed and analyzed with software GeneMapper v4.0 (Applied Biosystems).

Data analysis

The mating system was analyzed using mixed mating and correlated mating models with software MLTR 3.4 (Ritland 2004). Analyses were based on probabilities of maximum expectations "ME". The following indices were estimated: multilocus outcrossing rate (t^m>), single-locus outcrossing rate (t^s>), biparental inbreeding or mating between relatives (t^m-t^s>), multilocus paternity correlation (r^p(m)>), selfing correlation (r^s>) and maternal fixation index (F^m>). The t m of each family was estimated by the Moment method. The 95% confidence interval (CI) of each index was calculated from 1000 bootstrap replications, where the sampling units were represented by plants within progenies for the individual analysis and by progenies for population analysis. The mean effective number of pollen donors was estimated by N^ep=1/rp(m)> (Ritland 1989) and the average proportions of pairwise self-sibs (P^ss>), half-sibs (P^hs>), full-sibs (P^fs>), and self-half-sibs (P^shs>) within families were estimated as: P^ss=s^2>;P^hs=t^m2(1-r^pm)>; P^fs=t^m2r^pm)> and P^shs=2s^t^m> (Sebbenn 2006, Wadt et al. 2015). The coancestry coefficient within family (Θ) was estimated by Θ=0.1251+F^m4s^+t^m2+t^ms^r^s1+r^p(m)>, where s is the selfing rate (s^=1-t^m>) and variance effective size within family, assuming an idealized reference population (infinite size, random mating, without selection, mutation or migration) : N^e=0.5/Θn-1/n+1+F^o/2n>, where n is the number of analyzed offspring within families (we used the average of 227 seedlings, n= 17.46) and F^o> is the level of inbreeding within families (Sebbenn 2006). The number of mother trees (m) for seed collection, aiming to retain the reference effective population size of 150 was calculated by m=Ne(reference)/Ne=150/Ne m=N^ereference/N^e=150/N^e>> (Sebbenn 2006), based on three assumptions: i) non-relatedness of the mother trees; ii) no intermating of the sampled mother trees; iii) no overlapping of the pollen pools and no pollen from same fathers received by the mother trees (Sebbenn 2006).

RESULTS AND DISCUSSION

The maternal fixation index (F^m>) was not different from zero, indicating no inbreeding level of the mother trees. A similar result was found for Malpighia emarginata (Lopes et al. 2002). The multilocus outcrossing rate (t^m>) and single-locus outcrossing rate (t^s>) were equal to 1 (Table 1), indicating that the progenies were originated by outcrossing and that the species is allogamous. A similar finding was reported for the palm Astrocaryum aculeatum (Ramos et al. 2011). The difference between the multilocus and the single-locus outcrossing rate measured as the mating rate among related individuals (t^m-t^s>), was zero (0.001), evidencing the absence of mating between relatives of the sampled seedlings.

Table 1
Mating system at the population level

The estimated selfing correlation (r^s>) among progenies was low (0.11), indicating low variation in t^m> among mother trees. Accordingly, the t^m> among seed trees ranged from 0.79 to 1, in other words, it was significantly lower than 1.0 in 11 of the 13 families (Table 2), showing that some seedlings were originated by inbreeding. This result also suggests that the species is self-incompatible (Moraes et al. 2018). The individual variation in the outcrossing rate may be attributed to flowering asynchrony.

The paternity correlation (r^p(m)>) was moderate (0.293) and indicated that about 3.4 effective pollen donors (N^ep>) fertilized the seed trees in the investigated reproductive event (Table 1). These results show that the families are mainly composed of half-sibs (71%) and that mating was not random. The highest N^ep> was reported for the palm A. aculeatum (5.7, Ramos et al. 2011). Correlated matings can be caused by the behavior of pollinators that systematically visit nearby trees (Spoladore et al. 2017). Non-random matings were also evidenced by the variation among families for mating among relatives (t^m-t^s>: between 0.05 - 0.35), paternity correlation rates (r^p(m)>: 0.08 - 0.58) and effective number of pollen-donor trees (N^ep>: 1.7 - 11.9, Table 2).

The mean coancestry coefficient within progenies (Θ = 0.162) was higher than expected in half-sib progenies (0.125). Thus, estimates of addictive genetic variance and heritability must be calculated using a relatedness coefficient (Sobierajski et al. 2006) of 0.324 (Θ 2) instead of 0.25. Knowledge about the coancestry coefficient is also important when estimating the variance effective size (N^e>), which was lower (N^e>= 2.76) than expected in the random mating populations (4, Furlani et al. 2005). In infinite samples of progeny structures, the variance effective size varies from ​​1 to 4, where value 1 indicates selfed progenies, 2 full-sibs and 4 half-sib progenies (Sebbenn 2006). Due to the estimated effective population size, the seeds must be collected from at least 54 trees to retain the effective reference population size of 150 in progeny array samples. These analyses are important for estimating sample sizes in breeding, genetic conservation and seed collection programs addressing environmental recovery, as well as for the monitoring of genetic diversity in manipulated populations (Moraes et al. 2018).

CONCLUSION

The mating system indices estimated for E. precatoria in this study indicated that the species is allogamous but self-compatible. The studied progenies were mainly represented by half-sibs, but matings were not random due to the occurrence of some correlated mating, resulting in a few full-sibs within progenies.

ACKNOWLEDGEMENTS

The authors wish to thank Professor Paulo Yoshio Kageyama (in memoriam) of the Escola Superior de Agricultura “Luiz de Queiroz”/Universidade de São Paulo (ESALQ/USP), Departamento de Ciências Florestais, for his contribution as coordinator and advisor to this research. The authors are also indebted to the São Paulo Research Foundation (FAPESP), for scholarships of SLFR and GD; National Council for Scientific and Technological Development - CNPq for scholarships of MTGL and AMS and a post-doctoral scholarship (grant # 150297/2018-1) of GD, and to FAPESP (Grant #2014/10947-8); and the Fundação de Amparo à Pesquisa do Estado do Amazonas - FAPEAM (process 062.00669/2015), for financial support.

REFERENCES

  • Bussmann RW and Zambrana NYP (2012) Facing global markets - usage changes in Western Amazonian plants: the example of Euterpe precatoria Mart. and E. oleracea Mart. Acta Societatis Botanicorum Poloniae 81: 257-261.
  • Doyle JJ and Doyle JL (1990) Isolation of plant DNA from fresh tissue. Focus 12: 13-15.
  • FAO (1987) Especies forestales productoras de frutas y otros alimentos, 3. Ejemplos de América Latina. Estudio FAO Montes 44/3. FAO, Roma, 265p. Available at <Available at http://www.fao.org/docrep/015/an785s/an785s00.pdf > Accessed in March, 2018.
    » http://www.fao.org/docrep/015/an785s/an785s00.pdf
  • Furlani RCM, Moraes CMB, Moraes MLT, Paiva JR and Sebbenn AM (2005) Mating system in a Hevea brasiliensis population by isozyme loci. Crop Breeding and Applied Biotechnology 5: 402-409.
  • IBGE - Instituto Brasileiro de Geografia e Estatística (2018) Banco de dados agregados: Sistema IBGE de recuperação automática (SIDRA) - ano. Available at < Available at http://www.sidra.ibge.gov.br > Accessed in March, 2018.
    » http://www.sidra.ibge.gov.br
  • Kahn F (1991) Palms as key swamp forest resources in Amazonia. Forest Ecology and Management 38: 133-142.
  • Kang J, Thakali KM, Xie C, Kondo M, Tong Y, Ou B, Gitte J, Medina MB, Schauss AG and Wu X (2012) Bioactivities of açaí (Euterpe precatoria Mart.) fruit pulp, superior antioxidant and anti-inflammatory properties to Euterpe oleracea Mart. Food Chemistry 133: 671-677.
  • Küchmeister H, Silberbauer-Gottsberger I and Gottsberger G (1997) Flowering, pollination, nectar standing crop, and nectaries of Euterpe precatoria (Arecaceae) an Amazonian rain forest palm. Plant Systematics and Evolution 206: 71-97.
  • Lopes R, Bruckner CH and Lopes MTG (2002) Estimação da taxa de cruzamento da aceroleira com base em dados isoenzimáticos. Pesquisa Agropecuária Brasileira 37: 321-327.
  • Lorenzi H, Noblick L, Kahn F and Ferreira E (2010) Flora brasileira Lorenzi: Arecaceae (palmeiras). Instituto Plantarum, Nova Odessa, 384p.
  • Medina-Macedo L, Lacerda AEB, Sebbenn AM, Ribeiro JZ, Soccol CR and Bittencourt JVM (2015) Using genetic diversity and mating system parameters estimated from genetic markers to determine strategies for the conservation of Araucaria angustifolia (Bert.) O. Kuntze (Araucariaceae). Conservation Genetics 17: 1-10.
  • Moraes MA, Kubota TYK, Rossini BC, Marino CL, Freitas MLM, Moraes MLT, Silva AM, Cambuim J and Sebbenn AM (2018) Long-distance pollen and seed dispersal and inbreeding depression in Hymenaea stigonocarpa (Fabaceae: Caesalpinioideae) in the Brazilian savannah. Ecology and Evolution 9: 1-17.
  • Noda H (2012) In situ breeding and conservation of Amazonian horticultural species. In Borém A, Lopes MTG, Clement CR and Noda H (eds) Domestication and breeding: Amazonian species. Universidade Federal de Viçosa, Viçosa, p. 170-208.
  • Nogueira OL, Carvalho CJR, Müller CH, Galvao EUP, Silva HM, Rodrigues JELF, Oliveira MSP, Carvalho JEU, Rocha-Neto OG, Nascimento WMO and Calzavara BBG (1995) A cultura do Açaí. Embrapa. Centro de pesquisa Agroflorestal da Amazônia Oriental, Brasília, 50p.
  • Peel MC, Finlayson BL and McMahon TA (2007) Updated world map of the Köppen-Geiger climate classification. Hydrology and Earth System Sciences 11: 1633-1644.
  • Picanço-Rodrigues D, Astolfi-Filho S, Lemes MR, Gribel R, Sebbenn AM and Clement CR (2015) Conservation implications of the mating system of the Pampa hermosa landrace of peach palm analyzed with microsatellite markers. Genetics and Molecular Biology 38: 59-66.
  • Ramos SLF, Dequigiovanni G, Sebbenn AM, Lopes MTG, Kageyama PY, Macêdo JLVD, Matias K and Veasey EA (2016a) Spatial genetic structure, genetic diversity and pollen dispersal in a harvested population ofAstrocaryum aculeatumin the Brazilian Amazon. BMC Genetics 17: 1-63.
  • Ramos SLF, Dequigiovanni G, Lopes MTG, Veasey EA, Macêdo JLVD, Batista JS, Formiga KM and Kageyama PY (2016b) Microsatellite records for volume 8, issue 1: Microsatellite markers for Euterpe precatoria Mart. (Arecaceae) a palm species used by extractive traditional farmers of Amazon. Conservation Genetics Resources 8: 43-81.
  • Ramos SLF, Lopes MTG, Lopes R, Cunha RNVD, Macêdo JLVD, Contim LAS, Clement CR, Rodrigues DP and Bernardes LG (2011) Determination of the mating system of Tucumã palm using microsatellite markers. Crop Breeding and Applied Biotechnology 11: 181-185.
  • Ritland K (1989) Correlated matings in the partial selfer Mimulus guttatus Evolution 43: 848-859.
  • Ritland K (2004) Multilocus mating system program - MLTR Version 3.0. Vancouver. Available at <Available at http://genetics.forestry.ubc.ca/ritland/programs.htlm >. Accessed in March, 2018.
    » http://genetics.forestry.ubc.ca/ritland/programs.htlm
  • Schuelke M (2000) An economic method for the fluorescent labeling of PCR fragments. Nature Biotechnology 18: 233-234.
  • Sebbenn AM (2006) Sistema de reprodução em espécies arbóreas tropicais e suas implicações para a seleção de árvores matrizes para reflorestamentos ambientais. In Higa AR and Silva LD (eds) Pomares de sementes de espécies nativas. FUPEF, Curitiba, p. 193-138.
  • Sobierajski GR, Kageyama PY and Sebbenn AM (2006) Estimates of genetic parameters in Mimosa scabrella populations by random and mixed reproduction models. Crop Breeding and Applied Biotechnology 6: 47-54.
  • Spoladore J, Mansano VF, Lemes MR, Freitas LCD and Sebbenn AM (2017) Genetic conservation of small populations of the endemic tree Swartzia glazioviana (Taub.) Glaz. (Leguminosae) in the Atlantic Forest. Conservation Genetics 18: 1-13.
  • Wadt LHO, Baldoni AB, Silva VS, Campos T, Martins K, Azevedo VCR, Mata LR, Botin AA, Hoogerheide ESS, Tonini H and Sebbenn AM (2015) Mating system variation among populations, individuals and within and among fruits in Bertholletia excelsa Silvae Genetica 64: 248-259.

Publication Dates

  • Publication in this collection
    11 Apr 2019
  • Date of issue
    Jan-Mar 2019

History

  • Received
    23 Apr 2018
  • Accepted
    13 Nov 2018
location_on
Crop Breeding and Applied Biotechnology Universidade Federal de Viçosa, Departamento de Fitotecnia, 36570-000 Viçosa - Minas Gerais/Brasil, Tel.: (55 31)3899-2611, Fax: (55 31)3899-2611 - Viçosa - MG - Brazil
E-mail: cbab@ufv.br
rss_feed Acompanhe os números deste periódico no seu leitor de RSS
Acessibilidade / Reportar erro