SSR markers to detect gene flow from upland to mocó cotton

Marcadores SSR para detecção de fluxo gênico de algodoeiro herbáceo para algodoeiro mocó

FITOTECNIA

SSR markers to detect gene flow from upland to mocó cotton¹

Marcadores SSR para detecção de fluxo gênico de algodoeiro herbáceo para algodoeiro mocó

Ana Carolina Assis Dantas^I; Paulo Augusto Vianna Barroso^II; Lúcia Vieira HoffmannII, ^*; Milena Ferreira Alves^III; Francisco Pereira de Andrade^IV

^IBolsista da CAPES, Programa de Pós Graduação em Fitotecnia, Universidade Federal do Semi-árido, Mossoró-RN, Brasil, [e-mail]

^IINúcleo Cerrado- Embrapa Algodão, Rodovia GO-462, km 12 Zona Rural Caixa Postal 179, Santo Antônio de Goiás-GO, Brasil, 75.375-000, [e-mail], [e-mail]

^IIIBolsista da CAPES, Programa de Pós Graduação em Genética e Melhoramento de Plantas, Universidade Federal de Goiás, Goiânia-GO, Brasil, [e-mail]

^IVEmbrapa Algodão - Campina Grande-PB, Brasil, [e-mail]

ABSTRACT

Mocó cotton was the most important crop in Northeastern Brazil. Local varieties had been developed before governmental breeding programs, which should be identified and preserved as a genetic resource. The aim of this work was to select SSR markers which can be used to monitor gene flow from herbaceous to mocó cotton plants, using non transgenic upland cotton as a pollen source. Unique mocó cotton feral populations were identified in the state of Rio Grande do Norte, and mother plants and their seeds were collected. Allele frequencies were studied among fifteen plants from these populations using 64 SSR primer pairs. Ten SSR markers, with alleles exclusive to mocó or upland genotypes, were selected to monitor gene flow by paternity analysis (BNL3627; 2960; 2572; 3261; 3398; 3948; 3502; 3646 and CIR 094; 097). The collected seeds had none of the upland cotton alleles. The absence of exclusive alleles from upland cotton in the offspring showed that gene flow to mocó cotton was absent or lower than the analysis could detect. Among the mother plants, the number of homozygous individuals at four polymorphic loci was greater than expected, showing that reproduction occurs preferentially by self-fertilization or crossing among related individuals. Paternity analysis is an accurate measure of gene flow, and the methodology can be used for monitor these populations or others. Preservation of the plants is spoiled mainly by drought.

Key words - Biosafety. Molecular markers. Genetic resources.

RESUMO

O algodoeiro mocó foi a planta cultivada mais importante do Nordeste do Brasil. Variedades locais desenvolvidas anteriormente aos programas de melhoramento governamentais devem ser identificadas e preservadas como recurso genético. O objetivo deste trabalho foi selecionar marcadores SSR para monitorar o fluxo de genes do algodoeiro herbáceo para o mocó, usando plantas de algodoeiro herbáceo não transgênicas como fonte de pólen. Populações ferais de algodoeiro mocó raras foram identificadas no estado Rio Grande do Norte, e plantas mães e suas sementes foram coletadas. As freqüências alélicas foram estudadas nas plantas mãe utilizando 64 marcadores SSR. Dez marcadores SSR com alelos exclusivos dos algodoeiros mocós ou dos algodoeiros herbáceos foram usados para monitorar o fluxo gênico pela análise de paternidade (BNL3627; 2960; 2572; 3261; 3398; 3948; 3502; 3646 and CIR 094; 097). As sementes coletadas não possuíam nenhum alelo de algodoeiro herbáceo. A ausência de alelos exclusivos de algodoeiro herbáceo na progênie mostrou que o fluxo gênico para o algodoeiro mocó era ausente ou menor do que a análise pode detectar. Entre as plantas mãe, o número de indivíduos homozigotos em quatro locos polimórficos era maior que o esperado, mostrando que a reprodução ocorre preferencialmente por autofecundação ou cruzamento entre indivíduos aparentados. O teste de paternidade mede com acurácia o fluxo de genes, e pode ser utilizada para monitorar esta ou outras populações. A preservação das populações é ameaçada principalmente pela seca.

Palavras-chave - Biossegurança. Marcadores moleculares. Recursos genéticos.

Introduction

Cotton was the main agricultural economical activity for Brazilian semiarid region, at the North-east region of the country, during the twentieth century (IBGE, 1979). The cotton type used for cultivation is known as mocó (Gossypium hirsutum r. marie galante Hutch), which cultivation started around 1860, using seeds developed by farmers (MOREIRA et al., 1982).

The municipality of Acari, in the semi-arid-region, is probably the center of origin of at least part of local varieties of the mocó cotton planted in the past in Northeastern Brazil, since, after the drought of 1877 to 1879, in a small property next to a pond at Acari, called Olho d'Água da Seriema - where a farmer, noticing the drought resistance of a cotton plant, started to preserve it and obtain its seeds, distributing them to the neighborhood (MOREIRA et al,1989). Furthermore, the first report about mocó, written in 1860, described it growing like wild-type plants at Acari (MOREIRA et al, 1994).

Mocó cotton genotypes are adapted to the cultivation in this semi arid region. Cultivation in drought conditions is viable, since the plant maintain vigorous and productive for until eight years.

Local varieties are important genetic resources (CAMPBELL et al, 2003) and collections were intensified due to the accelerated decline of the in farm maintenance. Various local varieties have been planted in Northeast Brazil, which coexisted varied developed by governmental breeding programs, developed by Superintendência do Desenvolvimento do Nordeste (Sudene), Centro Nacional de Pesquisa da Algodão (CNPA, currently Embrapa Algodão) and Universidade Federal do Ceará (MOREIRA et al., 1972).

The release of transgenic cotton in Brazil was done considering that gene flow for the traditional mocó cotton varieties could possibly occur (BARROSO et al, 2005), but current gene flow levels between commercial and local varieties has not been determined. The main cotton type which is sympatric to mocó cotton is the upland cotton, G. hirsutum r. latifolium. G. mustelinum is also sympatric, although in Rio Grande do Norte the population is very small (BARROSO et al., 2010). G. barbadense is frequently present in gardens as an ornamental or medicinal plant (ALMEIDA et al, 2009).

Given that high variations can be found when gene flow is measured for particular populations, even considering the same species (ELLSTRAND , 2003), monitoring gene flow is proposed as a tool to prevent possible undesirable environmental impacts. Gene flow from crops to local varieties may just increase slightly the variability in the target population, or otherwise it may harm population conservation, especially if genes or genomes hybridization increases or decreases adaptability or reproduction taxes. Hybrid fitness higher than the fitness of both parents may lead to progenitor extinction by replacement (HEDGE et al, 2006), while if hybrids are less fertile than parents, populations may shrink (HAYGOOD et al., 2003). Paternity analysis may provide accurate information about gene flow (NEFF et al, 2000). When transgenic plants are not available or are not allowed, conventional non transgenic cultivars of the same species of the potential transgenic cultivars may be a confident model to measure gene flow, since plant transformation per se would change the amount of transgene flow. SSR markers, also known as microsatellite markers, are largely used in paternity analysis, and are well established for cotton genetic studies (OLIVEIRA et al., 2010). This report aimed to assess the existence of gene flow towards feral mocó cotton populations thought paternity analysis of offspring.

Materials and methods

Plant collections

With the aim to find mocó plants which are certainly traditional cultivars expeditions were organized to Acari (06°26'08'' S - 36°38'20''W), to the hills next to Seridó river Valey (MOREIRA et al, 1989), where this kind of plants have been previously found (FREIRE; MOREIRA, 1991), and the neighbor municipality Santana do Seridó (06°46'14'' S - 36°44'00'' W), in expeditions during the year 2003.

Shoot buds, seeds and leaves were collected from 32 different plants, 26 from Acari and 6 from Santana do Seridó (TAB. 1).

Plants were identified individually, and seeds and shoots of the same plants received the same designations. Shoots were collected and maintained in humidified plastic bags during transportation. Then, they were treated with indol acetic acid 1 g L^-1 for ten minutes, just before planting them in pots in greenhouses. Seeds were planted in autoclaved sand and, after germination, transplanted to pots.

The varieties that could have possibly acted as a source of pollen are upland cotton (G. hirsutum var latifolium), which have been planted in Northeastern Brazil semiarid region. Therefore two varieties known to be planted in this region were used for comparison of allele size and frequency: CNPA 8H and Precoce 3. These varieties both developed by Brazilian Agriculture research corporation (Embrapa). A third upland cotton Guazuncho, developed in Argentina, was chosen as a control because the allele sizes of this cultivar are known, since it has been used for molecular mapping (NGUYEN, 2004).

Allele frequencies and marker selection

The calculation of allele frequencies among these populations was made from DNA collected from leaves of 15 adult plants. These leaves, directly collected in the fields where the plants had been found, were placed in 50 ml tubes containing TE (Tris EDTA pH 8.0; 0.5 M) and maintained in boxes in ice for transportation to the laboratory, where the entire tubes containing the lives were frozen to -20ºC. For DNA extraction, leaf discs (50 to 100 mg from leaves of each plant) were placed in microtubes, frozen in liquid nitrogen and ground.

Buffer containing CTAB 2%; NaCl 1.4 M; EDTA 0.2 M; Tris HCl 0.1 M pH 8.0; PVP 2% and 2-mercaptoetanol 0.2% was added (600 µL). The suspension was incubated at 65ºC for 30 minutes, then 600 µL of a mixture of chloroform, isoamilic alcohol 24:1 (CIA). Was added, mixed by tube inversion and centrifuged at 12000 g for 10 minutes. DNA was recovered from supernatant and transferred to a new tube. After a second extraction by CIA, the collected DNA was precipitated with isopropanol, sequentially washed with 70% and 100% ethanol, and suspended in 100 µL Tris EDTA buffer. The concentration was estimated on agarose gels, and adjusted to 10 ng µL^-1.

The PCR reactions to obtain Simple Sequence Repeats (SSR) markers contained 20 ng of genomic DNA, 0.2 µM of each primer, 4 mM MgCl₂, PCR reaction buffer (10 mM Tris-HCl pH 8.3 and 50 mM KCl), 0.2 mM of dNTP and 1 U of Taq DNA polymerase.

Amplification reactions were achieved using a touchdown program, with an initial denaturation step at 94 ºC for 12 minutes, followed by 10 cycles with denaturation at 94 ºC for 15 s, in which annealing temperature decreased 1 ºC each cycle, from 65ºC to 55 ºC, for 30 s; followed by 72 ºC by 1 minute. Additional 35 cycles were at 94 ºC for 15 s, 55 ºC for 30 s and 72 ºC for 1 minute. A final elongation step at 72 °C for 5 minutes was added.

Reaction products were mixed to a loading buffer containing formaldehyde 95%; EDTA 0.5 M - pH 8.0, bromophenol blue and xylene cyanol 1 mg.mL^-1, and heated to 95 °C for 5 minutes prior to loading the samples to 6% acrylamide gels under denaturing conditions. The gels were run at 80 watts for 2 hours and 30 minutes and then stained with silver nitrate (CRESTE et al., 2001). The band patterns led two measure allele frequencies and the number of exclusive alleles among herbaceous cotton or mocó cotton genotypes. The expected (He) and observed (Ho) heterozygosities were calculated using the softwares GDA-Genetic Data Analysis version 1.1 and Popgen 1.31 (EXCOFFIER; HECKEL, 2006).

Gene flow

Gene flow was estimated by paternity analysis, comparing the DNA from the mother plants to the DNA obtained from seedlings germinated from seeds collected from the same mother plant. The DNAs were extracted from these 12 families, each constituted by the mother plant and 15 offspring plants.

Seeds received the same designation of the mother plant when they were collected. They were germinated in the laboratory and transplanted to field. The DNA was extracted the same way as described for the mother plants.

For this comparison just a selected set of primers was used. The selected primer pairs were those which revealed alleles that were present exclusively at mocó cotton plants, and absent at upland control plants of the varieties CNPA 8H, Precoce 3 e Guazuncho. The PCR amplification procedure and gel analysis were the same as described for the primer pair selection. The analyzed loci are presented in Table 2.

Results and discussion

Plants remained from abandoned fields of mocó cotton were localized in the municipalities of Acari and Santana do Seridó, state of Rio Grande do Norte. These two fields, for at least twenty years, had not been used for new plantations. Old and young cotton plants could be seen, and native vegetation grew among them.

The age of the plants was estimated as approximately 20 to 50 years, according to morphological characteristics and interviews with inhabitants in the nearby area, who could remember approximated planting dates. Young plants could be rarely observed.

The plants were harmed by drought when collected. Seeds were frequently damaged, what may have been caused by rodents, since excrements of mocó, a rodent whose excrements resembles cotton seeds, and from which the name of the cotton race can be derived (STEPHENS, 1967), could be observed at the time of our collection.

Considering visits to this region, which have been made approximately each five years during the last twenty years, the number of plants is diminishing.

The analysis for the presence of different alleles revealed that, among sixty four primer pairs, forty four (68.75%) presented the same alleles for all genotypes. Twenty loci were polymorphic between these 15 mother plants or the three upland cotton varieties used as controls.

The loci CIR 251, CIR170, CIR043 and BNL 1434 were polymorphic inside the mocó population (TAB. 3), and presented a medium number of alleles per locus of 3.4. The polymorphic mocó alleles, their frequencies, the expected and observed heteregozyties and the intrapopulation fixation indexes are presented in Table 3. A homozygous excess is observed in most of the loci.

The polymorphic loci CIR 364, CIR 121, CIR 249, BNL 3902, BNL 3816, CIR 251, CIR 170, CIR 049, CIR 043, BNL 1434 presented at least one common allele between mocó and herbaceous cotton. Therefore, they were not selected to paternity analyses studies.

Private alleles of mocó population were identified for ten primer pairs, which lead to the distinction between mocó and upland cotton genotypes. These alleles are presented in Table 4. Among them, eight primer pairs, BNL3627; 2960; 2572; 3261; 3398; 3948; 3502 and 3646, amplified only one locus each. Two primer pairs, CIR 094 and 097, amplified two loci each, both of them showing alleles which led two distinguish between the two cotton types.

The seeds that were collected from the two feral populations presented low germination rates. The plants that had not at least 15 germinated seeds were not considered for this study. Nevertheless, 12 plants (mothers) had at least 15 offspring plants. Among these, 11 were from Acari (61; 45; 46; 27; 28; 30; 32; 48; 49; 50; 59) and only one from Santana do Seridó (36).

These twelve loci were used to compare mother plants with their offspring. All the fifteen offspring plants, at all the twelve families and twelve studied loci, always presented only the mother plants alleles. None heterozygous plants were found. These suggest the absence of gene flow from upland to mocó cotton, for these loci, in this generation and in this population.

The populations used in these studies are unique among others found trough Northeast Brazil, mainly due to the relatively height number of plants and age. Although impediments for regeneration could be demonstrated, young plants were observed, showing that population had the ability of self propagation; hence it can be considered feral. Nevertheless, it is known that regeneration is rare, and populations are declining, as noticed by periodical visits to these places.

The observed growth of native vegetation among old fields is relatively common at this region, and in some occasions regenerated plants are used as pasture by caprines or bovine cattle. They are called capoeiras. In these two fields caprines were not found, probably due to the considerable declivity at the place of collection. Therefore plants were relatively well preserved. Another cotton species, G. mustelinum, which is endemic of Northeast Brazil, is also frequently used by cattle as a pasture, and grows much better when protected from them (BARROSO et al., 2010).

Cotton has a mixed matting system. The low numbers of heterozygote and the positive inbreeding coefficients suggest that reproduction occurs preferentially through selfing or between related individuals. Crosspollination may have been restricted by the distance among the plants or absence of efficient pollinators. The same has been found in other mocó populations in Northeast Brazil (MENEZES et al., 2010).

Gene flow by pollen could be achieved by insect pollination. The most frequent floral visitors of cotton belongs to the genus Apis sp., responsible for about 50% of cotton floral visits in Brazil (SILVEIRA, 2003). Pollen transportation by honey bees can achieve distances of 10km (BARROSO; FREIRE, 2003).

Even though upland and mocó cotton are sympatric, thus gene flow is possible, upland cotton alleles were not found among mocó cotton plants, showing that they are reproductively isolated. Paternity analysis can lead to accurate detection of gene flow (ELLSTRAND , 2003) and the data presented lead to the conclusion that pollen gene flow could be noticed in the studied generation.

If an unidentified amount of gene flow occurs, it is low, so it is supposed to be not sufficient to harm mocó cotton genetic identity.

Twelve identified loci can be used to paternity inference, differentiating between mocó and upland cotton. These loci that were shown to distinguish between upland and mocó cotton plants can be used for monitoring gene flow for these populations, considering the interest of in situ preservation, as a genetic resource for cotton. The other four primer pairs shown to be polymorphic among mocó cotton plants can be further used to diversity studies. Other loci shown to be polymorphic among upland cotton and mocó cotton could also be evaluated (ALVES et al., 2009).

Conventional non transgenic plants were used to estimate gene flow, and it is expected that if existed any cotton transgenic plant growing at these places, the estimated gene flow would be equivalent to the one measured for the conventional analysis, since the transgenes that have been used for cotton, conferring resistance to insects or herbicide tolerance, bear pollen grains which have the same characteristics of the conventional genotypes, and these transgenes do not interfere with pollen transportation. The estimative of potential transgene flow using conventional genotypes as a model can be useful when transgenes are not available or in ex-ante biosafety experiments.

Genetic preservation of the traditional mocó cultivars could be harmed not only by transgenes, but also by gene flow from conventional non transgenic cotton cultivars. These genes, possibly acquired through gene flow, could also confer favorable or non favorable traits, which would be a potential harm to the genetic maintenance. In this sense, paternity analysis is a broader genetic analysis than transgenic tracking as a tool to monitor gene flow to populations.

In situ maintenance of these genotypes can be achieved, including the use of exclusion zones for transgenic cotton planting (BARROSO et al, 2005). Moreover, collections have been made for ex situ preservation. They survive and reproduce under harsh climate conditions, and since they are of the same species as upland cotton, they might constitute a relevant genetic resource.

Conclusion

Paternity analysis, an efficient methodology do monitor gene, was applied to two populations of mocó cotton, which historical evidences indicates to have unique genetic traits, showing that gene flow from upland cotton has been absent or low, therefore representing low danger or none danger for populations preservation. In situ observations showed that drought spoils the populations.

Acknowledgments

We are grateful to José Henrique Assunção for technical assistance. This work was conducted with Ibama authorization number 13973. The International Project on GMO Risk Assessment Methodologies (GMO ERA Project) coordinated by Eliana G. Fontes, David A. Andow and Angelika Hilbeck has provided useful information and discussions.

Recebido para publicação em 26/01/2011; aprovado em 26/07/2011

• ALMEIDA, V. C. et al. In situ and genetic characterization of Gossypium barbadense L. from the States of Pará and Amapá, Brazil. Pesquisa Agropecuária Brasileira, v. 44, n. 07, p. 719-725, 2009.
• ALVES, M. F. et al. Marcadores moleculares polimórficos entre algodoeiros mocos e herbáceos. Revista Ciência Agronômica, v. 40, n. 03, p. 406-411, 2009.
• BARROSO, P. A. V. et al. Zonas de Exclusão de Algodoeiros Transgênicos para Preservação de Espécies de Gossypium Nativas ou Naturalizadas. Campina Grande: Embrapa Algodão, 2005. (Comunicado Técnico, 242).
• BARROSO, P. A. V. et al. In situ conservation and genetic diversity of three populations of Gossypium mustelinum Miers ex Watt. Genetic Resources and Crop Evolution. v. 57, n. 03, p. 343-349, 2010.
• BARROSO, P. A. V.; FREIRE, E. C. Fluxo gênico de algodão no Brasil. In: PIRES, C. S. S.; FONTES, E. M. G. e SUJII, E. R. Impacto ecológico de plantas geneticamente modificadas Brasília: Embrapa Recursos Genéticos e Biotecnologia, 2003. cap. 7, p. 163-194.
• CAMPBELL, B. T. et al. Status of the Global Cotton Germplasm Resources. Crop Science, v. 50, p. 1161-1179, 2010.
• CRESTE, S.; TULMANN, A.; FIGUEIRA, A. Detection of single sequence repeat polymorphisms in denaturing polyacrylamide sequencing gels by silver staining. Plant Molecular Biology Reporter, v. 19, n. 04, p. 299-306, 2001.
• ELLSTRAND, N. C. Current knowledge of gene flow in plants: implications for transgene flow. Philosophical Transactions of the Royal Society of London, v. 358, n. 1434, p. 1163-1170, 2003.
• EXCOFFIER, L.; HECKEL, G. Computer programs for population genetics data analysis: a survival guide. Nature Reviews Genetics, v. 07, p. 745-758, 2006.
• FREIRE, E. C.; MOREIRA, J. A. Relações genéticas entre o algodoeiro Mocó e diferentes. espécies e raças de algodoeiro. Revista Brasileira de Genética, v. 14, p. 393-411, 1991.
• HAYGOOD, R.; IVES, A. R.; ANDOW, D. A. Consequences of current gene flow from crops to wild relatives. Proceedings of the Royal Society of London, v. 270, p. 1879-1886, 2003.
• HEGDE, S. G. et al. The evolution of California's wild radish has results in the extinction of its progenitors. Evolution, v. 60, n. 06, p. 1187-1197, 2006.
• INSTITUTO BRASILEIRO DE GEOGRAFIA E ESTATÍSTICA (IBGE). Censos Econômicos de 1975 - Censo Agropecuário, Brasil. Rio de Janeiro: Fundação Instituto Brasileiro de Geografia e Estatística, 1979. 472 p.
• MENEZES, I. P. P. et al. Genetic diversity of moco´ cotton (Gossypium hirsutum race marie-galante) from the northeast of Brazil: implications for conservation. Botany, v. 88, n. 08, p. 765-773, 2010.
• MOREIRA, J. A. N.; BELTRÃO, N. E.; FREIRE, E. C. Organografia do algodoeiro mocó e sua relação com o crescimento e a produção. Campina Grande: Empresa Brasileira de Pesquisa Agropecuária, 1994. 66 p.
• MOREIRA, J. A. N.; SILVA; F. P.; ALBUQUERQUE, J. J. L. Melhoramento genético do algodão "mocó" no estado do Ceará-Brasil. Ciência Agronômica, v. 02, p. 19-23, 1972.
• MOREIRA, J. A. N. et al. Visão retrospectiva do melhoramento genético no algodoeiro mocó (Gossypium hirsutum r. marie galante Hutch) no nordeste do Brasil. Campina Grande: Centro Nacional de Pesquisa do Algodão, 1982. (Documento, 12).
• MOREIRA, J. A. N. et al. Algodoeiro Mocó: Uma Lavoura Ameaçada de Extinção. Campina Grande: Embrapa Algodão. 1989. (Documento, 36).
• NEFF, B. D.; REPKA, J.; GROSS, M. R. Statistical confidence in parentage analysis with incomplete sampling: how many loci and offspring are needed? Molecular Ecology, v. 09, n. 05, p. 529-539, 2000.
• NGUYEN, T. B. et al. Wide coverage of the tetraploid cotton genome using newly developed microsatellite markers. Theoretical and Applied Genetics, v. 109, n. 01, p. 167-175, 2004.
• OLIVEIRA, T. S. et al Validação de métodos laboratoriais aplicados a análises com marcadores microssatélites. Revista Ciencia Agronomica, v. 41, n. 02, p. 279-284, 2010.
• SILVEIRA, F. A. As abelhas e o algodão Bt no Brasil: uma avaliação preliminar. In: PIRES, C. S. S.; FONTES, E. M. G. E SUJII, E. R. Impacto ecológico de plantas geneticamente modificadas: o algodão resistente a insetos como estudo de caso. Brasília: Embrapa Recursos Genéticos e Biotecnologia, 2003. cap. 8, p. 195-215.
• STEPHENS, S. G. Evolution seeds domestication of the new world cottons (Gossypium spp.). Ciência e Cultura, v. 19, p. 118-134, 1967.

*

Autor para correspondência

1

Pesquisa financiamenta pela Embrapa (Rede Biossegurança) e Ana Carolina Assis Dantas que foi bolsista PIBIC/CNPq

Publication Dates

• Publication in this collection
  20 Oct 2011
• Date of issue
  Mar 2012

History

• Received
  26 Jan 2011
• Accepted
  26 July 2011