G3: Genes|Genomes|Genetics Early Online, published on May 25, 2017 as doi:10.1534/g3.117.041509

Assessment of gene flow between hirsutum and G. herbaceum: evidence of unreduced gametes in the diploid progenitor.

Montes E *1, Coriton O †, Eber F †, Huteau V †, Lacape JM ‡, Reinhardt C §2, Marais

D §, Hofs JL **, Chèvre AM † 2, Pannetier C * ‡ 2

*: Institut Jean-Pierre Bourgin, INRA, AgroParisTech, CNRS, Université Paris-Saclay,

RD10, F-78026 Versailles Cedex, France

†: Institut de Génétique, Environnement et Protection des Plantes, INRA, Agrocampus

Ouest, Université de Rennes I., BP35327, 35653 Le Rheu, France

‡: CIRAD, UMR AGAP, Amélioration Génétique et Adaptation des Plantes méditerranéennes et tropicales, 34398 Montpellier, France.

§: Department of and Soil Sciences, University of Pretoria, Pretoria 0001, South Africa

** : CIRAD, UR AIDA, Agro-écologie et Intensification Durable des cultures Annuelles,

34398 Montpellier, France.

1 present address: UMR 1349, IGEPP, INRA BP35327, 35653 Le Rheu, France

2 corresponding authors

Running title

Unreduced gametes in diploid

1

© The Author(s) 2013. Published by the Genetics Society of America.

Key words

Gene flow, natural hybridization, unreduced gamete, ,

Corresponding authors: UMR 1349, Institut de Génétique, Environnement et Protection des

Plantes, Institut National de la Recherche Agronomique (INRA), BP35327, F-35653 Le Rheu,

France. Email [email protected] and UMR1318, Institut Jean-Pierre Bourgin, INRA

F-78026 Versailles, France. Email [email protected].

2 Abstract

In the framework of a gene flow assessment, we investigated the natural hybridization rate between Gossypium hirsutum (AADD genome) and G. herbaceum (AA genome). The latter species, a diploid progenitor of G. hirsutum, is spontaneously present in South Africa.

Reciprocal crosses were performed without emasculation between G. herbaceum and G. hirsutum. Neither examination of the morphological characteristics nor flow cytometry analysis of the 335 resulting from the G. hirsutum x G. herbaceum cross showed any hybrid features. Of the 148 plants produced from the G. herbaceum x G. hirsutum cross, three showed a hybrid phenotype, and their hybrid status was confirmed by SSR markers. Analysis of DNA content by flow cytometry and morphological traits clearly showed that two of these plants were triploid (AAD). The third plant had a flow cytometry DNA content slightly higher than G. hirsutum. In addition, its morphological characteristics (plant architecture, presence and size of petal spots, leaf shape) led us to conclude that this plant was AAAD thus resulting from fertilization with an unreduced AA gamete of the female G. herbaceum parent.

Fluorescent In Situ Hybridization (FISH) and meiotic behavior confirmed this hypothesis. To the best of our knowledge, this is the first description of such gametes in G. herbaceum, and it opens new avenues in breeding programs. Furthermore, this plant material could provide a useful tool for studying the expression of genes duplicated in the A and D cotton genome.

3 INTRODUCTION

Assessment of gene flow from a crop to a wild relative is particularly relevant in the case of allopolyploid species, which are formed by hybridization between related species, when one of its progenitors is present in cultivated areas. In certain geographical areas, this may be the case for the cultivated tetraploid cotton (Gossypium hirsutum, AADD, 2n=4x=52) and a wild diploid species (Gossypium herbaceum, AA, 2n=2x=26). Cultivated cotton, the most important source of natural fibre, has been the long-standing subject of a number of taxonomical studies (Watt, 1907; Hutchinson et al. 1947; Saunders, 1961; cited in Wendel et al. 2010; Fryxell, 1979, 1992). The cotton genus (Gossypium L.) includes 50 diploid species

(2n=2x=26) that can be differentiated cytogenetically into eight genome groups (A-G & K), and seven allotetraploid species (2n=4x=52) (Fryxell, 1979; Fryxell, 1992; Stewart 1995,

Wendel and Croon, 2003, Grover et al. 2015, Gallagher et al. 2017). Remarkably, among the four domesticated species, there are two allotetraploids (the New World tetraploids G. hirsutum and G. barbadense) and two diploids (the Old World diploids G. arboreum and G. herbaceum). The tetraploid species originated from a single polyploidization event between an A-genome diploid species and a D-genome diploid species (Wendel and Cronn, 2003). The two African-Asian A-genome , G. arboreum and G. herbaceum, are the equivalent of the maternal genome donor to polyploid cotton (Wendel, 1989), although some evidence suggests that G. herbaceum may more closely resemble the A-genome donor than G. arboreum (Wendel et al. 2010). Recent results have suggested an independent domestication of these two Old World cotton species (Renny-Byfield et al. 2016). Grover et al (2015) demonstrated that this combination of a D-genome progenitor, resembling modern G. raimondii, and an A-genome progenitor, much like modern G. arboreum and G. herbaceum, occurred as a single hybridization event (genomic merger and doubling) that gave rise to all polyploid cottons. Interspecific hybridizations have been used extensively to improve yield,

4 fibre quality or other agronomic traits, and numerous studies have focused on the introgression of traits of interest from wild diploid cotton into cultivated cotton G. hirsutum

(Meyer, 1974; Stewart, 1995; Mergeai, 2004, Chee et al. 2016). Colchicine doubling is an essential step in producing an improved through interspecific crossing, as the recovered triploid hybrids are sterile. Unreduced gametes have proven useful in enabling crosses between plants of different ploidy levels which otherwise often fail due to unbalanced parental contributions in the developing (Barcaccia et al. 2003; Carputo et al. 2003). If the plant of the lower ploidy level can be induced to produce unreduced gametes, such limitations can be overcome. This strategy has been successfully used in various species such as potato, alfalfa and manioc (reviewed in Brownfield and Kohler, 2011; Dewitte et al. 2012).

The large-scale development of transgenic insect-resistant and/or herbicide-tolerant cotton has promoted studies about the ability of native diploid A or D species to produce fertile hybrids with tetraploid cotton (Stewart, 1995; Brubaker et al. 1999). In South Africa, and particularly in the KwaZulu Natal Province, where the commercialization of transgenic Bt cotton began in

1998, a wild species (G. herbaceum), one of the progenitors of G. hirsutum, is found in areas neighbouring the cultivated cotton fields. Consequently, we have undertaken a study to evaluate the likelihood of gene flow between the allotetraploid cultivated cotton (G. hirsutum) and its wild diploid relative (G. herbaceum). Transgene escape depends on the generation of fertile hybrids, and the possible selective advantage that would result in the development of weedy derivatives. In this paper, we present our analysis of natural cotton gene flow in South

Africa, which furthermore provides evidence for unreduced gametes in the species G. herbaceum.

5

MATERIALS AND METHODS

Plant material

Reciprocal crosses were performed without emasculation between G. herbaceum Indian

accession Boumi Aria (A1A1, 2n=26) and G. hirsutum (AtAtDtDt, 2n=52; lower case letter “t” referring to the “tetraploid” genome) Delta Opal RR from Deltapine, South Africa. Forty plants from each accession were maintained in individual pots in a greenhouse at Pretoria

University and divided into two batches: donor plants for pollen harvest and recipient plants for hybrid production. To simulate both pollination occurring in the field and regular visits by pollinators to plants, pollen harvested from newly opened was spread on the recipient flowers without emasculation two or three times between 9 and 11 am. All the subsequently produced were sown (one seed per pot) in greenhouses at the Versailles Centre of the

French National Institute of Agricultural Research (INRA) to screen for hybrids.

Seeds with a thick seed coat harvested from G. herbaceum were scraped with sandpaper to facilitate germination, and then placed on wet paper at 24°C. Germinating seeds were planted in soil in a greenhouse following radicle emergence.

Morphological examinations

Morphological examinations were made both at the vegetative stage (six leaves) and at the reproductive stage. Details of morphology were recorded. In plants that exhibited hybrid characteristics, particular attention was paid to flower size and color, and especially to the presence and intensity of red petal spots, as G. hirsutum has no spot and G. herbaceum shows large, bright red spots.

6 Pollen viability

Pollen viability was assessed after Alexander dye staining (Alexander, 1969) with a 3-hour incubation period at 50°C. Several samples were taken in the greenhouse between 10 and 11 a.m. on several distinct days.

Estimation of nuclear DNA amount by flow cytometry

Flow cytometry analyses (Partec ploidy analyser, Munster, Germany) were performed for all progeny on leaf tissue to assess their nuclear DNA amount, as described by Eber et al. (1997).

Barley leaves served as the internal reference with set value of 100. Preliminary experiments confirmed that the measure was not affected neither by leaf age or conservation delay (5 to 20 minutes) of shredded leaves in the buffer (Kit Partec, Munster, Germany). Three measurements with independent leaf samples were performed for each control and potentially hybrid plant. Twenty plants were analysed for the G. hirsutum and G. herbaceum parents. The two-fold size difference between the two genomes A and D was taken into account to determine chromosome numbers.

Meiotic behavior and chromosome counting

Young squares (flower buds) were sampled in the greenhouse from the two parental lines and for deemed hybrid plants according to flow cytometric analysis and morphological observations. They were fixed in Carnoy’s solution (ethanol–chloroform–acetic acid, 6: 3: 1) for 24 hours and preserved in 70% ethanol until characterized. For the observation of pollen mother cells (PMCs), anthers were squashed in 1% acetocarmin and cells inspected with a light microscope. PMCs were observed at metaphase I to establish the meiotic behavior of the hybrids.

SSR marker analysis

Microsatellite (SSR) markers were used to verify the hybrid status of progenies obtained from the different hybridizations. Six SSRs (BNL 3264, BNL 3563, BNL 409, BNL 4049, NAU

7 0915, NAU 1246) that demonstrated clear polymorphisms between G. hirsutum and G. herbaceum were selected, based on prior screens of diploid and tetraploid accessions

(unpublished data). Their primer sequences are available from CottonGen

(https://www.cottongen.org/). Sixty-six DNAs were analyzed, including 64 progenies from the 2 reciprocal crosses (46 for G hirsutum x G herbaceum, and 18 for G. herbaceum x G. hirsutum), and from control DNAs of each parental species. Polymorphic markers between the two species were kept for hybrid characterization. DNA extraction, PCR amplification and electrophoresis-based allele separation were performed according to Lacape et al (2007). The progenies were assigned an allelic “hybrid” coding when alleles of both G. hirsutum (could be

1 or 2 amplicons) and G. herbaceum (an additional unique amplicon) were detected.

Fluorescence in situ hybridization (FISH)

Mitotic chromosomal squashes were prepared with standard techniques. FISH experiments were carried out with two DNA probes: the 45S rDNA probe (pTa 71) is a cloned 9-kb EcoRI fragment of a rDNA repeat unit (18S-5.8S-26S genes and spacers) isolated from Triticum aestivum (Gerlach and Bedbrook 1979) and the pTa 794 probe is a 410 bp BamHI fragment of the 5S rDNA from wheat (Gerlach and Dyer 1980). The two probes were labelled by random priming with Alexa-488 dUTP (Invitrogen, Life Technologies) and Alexa 594 (Invitrogen,

Life Technologies) respectively. Genomic in situ hybridization (GISH) was performed with probe DNAs from the progenitors of the allotetraploids. Total genomic DNA from G. herbaceum (A-genome) was labeled by nick translation with biotin-16-dUTP (Roche

Diagnostics, Indianapolis, IN, USA). Total genomic DNAs from G. raimondi (D-genome) were autoclaved to yield fragments of 100-300 bp and used as a blocking DNA. A series of test assessments determined the amount of blocking DNA that was required to discriminate between the two genomes. Optimum results were obtained with a ratio of 1:50 for A genome probe and D genome blocking DNA in the hybridization mixture. Biotinylated probes were

8 immunodetected by Texas-red-conjugated with avidin antibodies (Vector Laboratories,

Burlingame, CA). The chromosomes were mounted and counterstained in Vectashield

(Vector Laboratories) containing 2.5µg/mL 4,6-diamidino-2-phenylindole (DAPI).

Fluorescence images were captured using a CoolSnap HQ camera (Photometrics, Tucson,

Arizona) on an Axioplan 2 microscope (Zeiss, Oberkochen, Germany) and analyzed with the

MetaVueTM software (Universal Imaging Corporation, Downington, PA).

RESULTS

Interspecific F1 hybrid Production

G. hirsutum (AtAtDtDt) x G. herbaceum (A1A1)

After pollination, moderate boll shedding was observed, as 60 mature bolls were obtained from 68 pollinated flowers. A total of 335 plants were obtained from 504 sown seeds. Flow cytometry analysis revealed that all of them were identical to the female allotetraploid parent,

G. hirsutum. This observation was confirmed by their morphological characteristics including plant architecture, leaf shape as well as flower size, shape and color. No interspecific hybrid was obtained.

G. herbaceum (A1A1) x G. hirsutum (AtAtDtDt)

Significant boll abscission was observed from the hybridizations involving G. herbaceum as the female (29 mature bolls out of 146 pollinated flowers). Similarly, a large number of aborted seeds were harvested (163 aborted seeds out of a total 498 harvested). Out of 335 seeds sown, 148 developed into plants: among these, three interspecific hybrids were detected. Flow cytometric analysis showed that two of them (namely H1 and H2) had an intermediate value possibly corresponding to AAD hybrids. The third hybrid (H3) exhibited a slightly higher DNA content than G. hirsutum.

9 Our results indicate a hybridization rate of 0.02 hybrid seed per pollinated flower, compared to ~30 seeds per boll observed after selfing of G. hirsutum. Among the 148 plants observed, the interspecific hybridization rate detected was 0.02.

Characterization of interspecific F1 hybrids

Morphological traits

Among the 148 plants obtained with G. herbaceum as the female parent, observations 20 days after sowing (at the three-leaf stage) led to the morphological identification of the three hybrids. Already at a young stage (three-leaf, 20 days after sowing), three plants exhibited general morphological traits similar to G. hirsutum. Several weeks later, two types of plants could be distinguished. Among the three hybrids, H1 and H2 displayed identical morphological characteristics and exhibited the phenotype typical of a triploid hybrid (AAD), with traits intermediate between the two parents (Fig 1). For the third hybrid (H3), the overall plant phenotype was intermediate between G. hirsutum and the triploid (Fig 1), with individual plant parts varying between the triploid and G. herbaceum. Hybrid H3 leaves had less marked serrations but more lobes and were intermediate between G. herbaceum and the triploid (Fig 2). In its first stages of development, the hybrid H3 displayed density of stem hairs and gossypol glands similar to those of G. herbaceum. The two types of hybrid phenotypes were also differentiated by flower characteristics: corolla size, petal red spot intensity and size and stamen color and size (Fig 3). The corollas of the H1 and H2 hybrid flowers had the same size as G. hirsutum, but they exhibited light red spots (Fig 3). H3 had larger, brighter red spots, intermediate between AAD and G. herbaceum (Fig 3), and its corolla size was slightly smaller than those of H1 and H2. These observations can support the hypothesis that the H3 hybrid carries more genome A than the G. herbaceum parent.

Moreover it has been reported that spot presence is controlled by a single gene carried by the

A genome (Stephens, 1974).

10

Fig 1: Plants at 70 days of culture. a: left to right G. hirsutum , G herbaceum, b: left to right, hybrid H1, hybrid H2, hybrid H3.

Fig 2. Leaf shapes. Left to right G. herbaceum, hybrid H1, hybrid H3, G. hirsutum. Bar=1cm

Fig 3. Flower characteristics: Left to right: G. herbaceum, hybrid H1, hybrid H3, G. hirsutum.

Bar=1cm

11 Microsatellite-based allelic status of the interspecific hybrids

One of the six SSRs tested was not exploited because of ambiguous allele assignment. For the five remaining SSRs (BNL3563, BNL409, BNL4049, NAU0915, and NAU1246), the hybrid status of hybrids H1, H2 and H3 was confirmed (Fig 4). The three hybrids exhibited the alleles from G. hirsutum as well as from G. herbaceum. However, the amplification profiles did not reveal any quantitative difference between these hybrids.

Fig 4: Molecular profiles of parents G. hirsutum, G. herbaceum and 20 progenies showing the hybrid profile of H1, H2 and H3. a :marker BNL409, b : marker NAU1246

Identification of the genomic structure of hybrids

Results of the FISH analysis at mitosis showed that G. hirsutum (AtAtDtDt) carried six 45S loci with 2 loci on At genome (Fig 5 a-b), in accordance with previous reports (Hanson et al

1996; Gan et al. 2013). However, six 45S loci were also observed on the A genome of G.

herbaceum (A1A1) (Fig 5c). The structure of H1 and H2 hybrids was confirmed with 39 chromosomes, and we found the expected number of 45S loci for an A1AtDt genomic structure, i.e. six in total with three from G. herbaceum and three from G. hirsutum (Fig 5d).

Analysis of H3 resulted in a count of 52 chromosomes carrying nine 45S loci six from G.

herbaceum plus three from G. hirsutum thus compatible with the genome structure A1A1AtDt

(Fig 5e).

12 Therefore, this hybrid was formed from a female unreduced gamete of G. herbaceum and a

reduced gamete of G. hirsutum (Table 1). According to the genomic composition (A1AtDt) of

H1 and H2 plants, the expected meiotic behavior is 13 univalents corresponding to the D genome and 13 bivalents between the homologous chromosomes of the A genomes. Thirty- five percent of the cells of 26 analyzed PMCs matched this configuration, but 57.7% revealed the presence of multivalents, trivalents or quadrivalents (Fig 6a, Table 2). In the third H3 hybrid, all the 15 PMCs observed indicated the presence of multivalents likely due to the presence of three homologous A genomes (Fig 6b, Table 2). Tetrads harbouring 3 to 7 microspores have been observed from young flowers of the triploid hybrids H1 and H2.

Fig 5: In situ hybridization analysis. (a) The GISH probe hybridized specifically to all A genome chromosomes (pink) in G. hirsutum (AADD); 5SrDNA in red, (b-e) FISH with 45S rDNA in green. FISH analyses of somatic metaphase chromosomes from G. hirsutum (b), G. herbaceum (c), hybrid H2 (d) and hybrid H3 (e). Chromosomes were counterstained with

DAPI (blue).

13

CHROMOSOMES 45SrDNA SIGNAL GENOME NUMBERS NUMBERS STRUCTURE

G. HIRSUTUM 52 At At Dt Dt 6 : 2 (At At),4 (Dt Dt)

HYBRID H1 39 6 : 1 At, 3 A1, 2 D A1 At D

HYBRID H2 39 6 : 1 At, 3 A1, 2 D A1 At D HYBRID H3 52 9 : 3 A1, 3 A1, 1 At, 2 D A1 A1 At D

Table 1: 45S rDNA signal numbers and genome structures

Fig 6: PMC in metaphase I. a from AAD hybrid H1, 13 univalents (1-13), 11 bivalents (I-

XI) and 1 quadrivalent; b from AAAD hybrid H3: 14 univalents (1-14), 16 bivalents (I-XVI) and 2 trivalents.

2n PMC UNIV. BIVAL. TRIVAL. QUADRIVAL. 12.61* 11.81 0.31 0.17 HYBRID H1 – H2 39 26 (9-14)** (7-14) (0-2) (0-2) 14.67 16.8 1.07 0.13 HYBRID H3 52 15 (12-17 (14-20 (0-4) (0-1)

Table 2: Meiotic behavior of hybrids H1, H2 and H3. *: average per cell, ** range

14 Fertility of interspecific hybrids

As expected, pollen grains from the parental lines had a high viability (Fig 7). In contrast, very few pollen grains (< 5%) from hybrids H1 and H2 were stained by Alexander dye (Fig

7). Pollen grains from the hybrid H3 exhibited a range of viability, up to 20%. Sampling was performed on all the plants at the same time, and the fertility of G. herbaceum was observed to be lower than that of G. hirsutum. This difference in the viability of the pollen grains from

G. hirsutum and G. herbaceum could be explained by the time of the sampling. Flowers were collected at the same time in the morning, which might not have been optimal for the sampling of G. herbaceum pollen.

None of the three hybrids produced any seed, either through open or self-pollination. In order to save the hybrid materials, the plants were reproduced through cutting when they became too large.

Fig 7: Pollen viability (percentage). Seven repetitions (40 to more than 1000 pollen grains per sampling) for G. hirsutum and H1, H2, and H3. One observation on 120 pollen grains for

G. herbaceum.

15 DISCUSSION

The results showed a relatively high hybridization rate (three hybrids among 148 progeny) between G. herbaceum and G. hirsutum without emasculation of the female parent and when the plant with the lower ploidy level was the female. This result is similar to those observed with Brassica napus, which is also an allotetraploid species; when one of its progenitors, B. rapa was used as female, the spontaneous interspecific hybridization ranged from 3 to 8%

(Warwick et al. 2003). When the diploid parent G. herbaceum was used as male, no hybrids were obtained. These results are not in accordance with those of Brubacker et al. (1999), working on the improvement of cultivated cotton by crosses with diploid Australian species

(genomes C, G and K). They found that it was more difficult to obtain triploid hybrids when the female parent was diploid. Similarly, Beasley (1940) concluded that the likelihood of successful hybridization is higher if the female parent has the higher ploidy level, reflecting the traditional recommendations concerning the critical role of ploidy ratio of endosperm and zygote for successful embryo development (Beasley, 1942; Stephens, 1942 cited in Ahoton et al. 2003). Similarly, it has been observed among Brassicaceae species that the production of interspecific hybrids using B. napus as female or male with different diploid species depends on the species combination (Chèvre et al. 2004). Our combined analyses based on cytometry, fluorescent in-situ hybridization and chromosome counting indicated that the three hybrids obtained were of two types. Two, H1 and H2, were triploids (AAD) and, remarkably, the third one (H3) was a tetraploid carrying three A genomes and a single D genome. The meiotic behavior of triploid hybrids is close to what can be expected with a majority of pollen mother cells with 13 bivalents corresponding to homologous pairing of the 13 A chromosomes, whereas the D chromosomes remained at univalent stage with 13 univalents. However, it is important to note that several cells showed the presence of multivalents, indicating the possibility of homoeologous pairing between A and D genomes, even if these exchanges were

16 difficult to detect after a long evolution in G. hirsutum (Page et al. 2016). In the third hybrid, the presence of more than two A genomes was shown by the morphological characteristics

(leaf shape, flower size and red spots on the corolla). Its genomic structure, A1A1AtDt, was confirmed by observation of metaphase in mitosis, revealing the presence of the nine 45S loci expected for this genomic structure. Additionally, the high frequency of multivalents was in agreement with the presence o f three A genomes.

Our primary objective was to estimate the risk of cross-pollination between allotetraploid cultivated cotton (G. hirsutum) and wild diploid cotton (G. herbaceum), since in the

Makhathini flats area of South Africa the diploid progenitor grows in proximity to the cultivated transgenic cotton fields. Most of published studies were concerned either with gene flow between transgenic and non-transgenic cultivated cotton in various areas of the world

(Zhang et al. 2005; Van Deynze 2005; Llewellyn et al. 2007; Heuberger et al. 2010; Loureiro et al. 2016), or between wild and cultivated tetraploid cotton (Wegier et al. 2011; Pereira et al. 2012). Our study on the hybridization rate between diploid and tetraploid cotton, showed that the three recovered hybrids produced no seeds and that the viability of their pollen was very low. The single seed was produced after crossing a triploid hybrid and G. hirsutum developed after hand-pollination on an emasculated flower of the triploid (data not shown). It is reasonable to conclude from the present study that the probability of obtaining offspring is very low. Moreover, only the species G. herbaceum var. africanum is endemic in South

Africa (van Wyk and van Wyk, 1997). Entomological investigations have shown that the beetle, Astylus atromaculatus could be an unexpected but efficient pollinator on G. hirsutum in South Africa (Pierre and Hofs, 2010). However, A. atromaculatus and domestic bees are not present on G. herbaceum (JL Hofs personal communication). In contrast, Anthophoridae are common on both G. hirsutum and G. herbaceum and have a great capacity for carrying pollen and may thus be a potential vector for pollen transfer between the two species (JL Hofs

17 personal communication). However, in the case of trials involving interspecific hybridization, reverting to a more fertile form required several generations of backcrosses with the cultivated species (Kammacher, 1966; Mergeai, 2004).

Our results have demonstrated the possibility of the formation of unreduced gametes. To our knowledge, this phenomenon as been poorly documented in cotton, while it is prevalent in diverse plant species and considered to be an important mechanism for polyploidization

(reviewed in Mason and Pires, 2015). The formation of unreduced gametes in Gossypium interspecific hybrids has been reported in crosses between G. hirsutum and Australian K- genome diploids G. exiguum or G. nobile (Brubaker et al. 1999). Observations of meiosis in these hybrids suggested that some may have derived from unreduced gametes arising from dysfunctions during M2, hence the duplication of the segregating chromosomes. Furthermore,

Skovsted (1934), working on the production of hybrids between diploid Asiatic and tetraploid

New World cotton, obtained a hybrid with 52 chromosomes and hypothesized the production of 2n gametes in the diploid parent. Likewise, in cytogenetic studies of cotton grown in Iran, the occurrence of large pollen grains was assumed to be potential unreduced gametes

(Sheidai, 2008; Noormohammadi et al. 2012), although other methods are considered to be necessary to confirm the correlation between size and ploidy in pollen (Dewitte et al. 2012).

For the same type of sterile hybrids (triploids 2n = 39, one tetraploid 2n = 52) obtained by crossing A-genome Asiatic cottons with tetraploid New World cottons, Amin (1940) reported that backcrossing to New World cottons gave occasional success in boll setting on the F1s.

Such bolls contained very few seeds, and the plants developed from them had varying degrees of fertility. In this context, by crossing one of our triploid hybrid (H1) as female with G. hirsutum (25 flowers emasculated and hand pollinated with G. hirsutum pollen) we obtained a single large seed. Its germination gave rise to a plant with 2n=65 chromosomes and a putative

A1AtAtDt structure (data not shown). The structure may be the result of unreduced gametes

18 from the triploid hybrid. This ability of Gossypium to produce such gametes from diploid parents and from F1 hybrids might point to an alternative explanation for the formation of the allotetraploid cultivated species, G. hirsutum. In fact, it has been hypothesized the formation of unreduced gametes offers the possibility to generate a triploid bridge which can be at the origin of polyploid species (reviewed in Mason and Pires 2015). Interspecific hybridization has often been carried out to provide new sources of germplasm for incorporation into breeding programs (Chee et al. 2016). However, it is not possible to obtain fertile hybrids from crosses between cultivated tetraploid G. hirsutum and wild diploids species without treating the floral bud at the time of fertilization with growth substances in order to prevent boll abscission (Meyer, 1974; Altman, 1988). Brubaker et al. (1999) have shown that, in the best case, the natural crossing of diploid wild Australian species with tetraploid cotton cultivars only resulted in the production of sterile hybrids. However, fertility can be restored by doubling the number of chromosomes using colchicine to produce hexaploid plants (2n=78 chromosomes). This approach has been widely used in interspecific hybridization to improve cotton plants, but such a doubling has never been observed occurring spontaneously in nature

(Mergeai, 2004, 2006). 2n gametes have already been exploited to create new cultivars at higher ploidy levels, but more interestingly, they have also been used in various species to create a bridge for the transfer of desirable genes from wild diploid species into the cultivated polyploid gene pool in various species (reviewed in Dewitte et al. 2012). This approach is particularly relevant for cultivated cotton, as the wild diploid species of Gossypium are rich in rare desirable attributes that are not available in the germplasm of cultivated species (Stewart,

1995; Mergeai, 2006). This methodology has been used to transfer genetic diversity from diploids through 2n gametes to polyploid crop varieties, as demonstrated for example in potato and alfalfa (Peloquin et al. 1999; Ramanna and Jacobsen, 2003). The production of autopolyploids to overcome the difficulty of transferring desirable traits from wild or

19 cultivated diploids to G. hirsutum was investigated by Mehetre et al. (2003). The use of unreduced gametes would facilitate the direct production of such AAAD hybrids either by direct crosses (H3 hybrid) or by crossing an AAD hybrid with G. hirsutum. Even though the production of unreduced gametes has rarely been described in cotton, various methods have been tested in other species to increase their occurrence (Brownfield and Kohler, 2011,

Dewitte et al. 2012, Younis et al. 2014). New strategies (Ravi et al. 2008; Chan, 2010; Ravi and Chan, 2010) developed on the model plant Arabidopsis thaliana (reviewed in Crismani et al. 2013) are considered by the authors to be possibly applicable to crops.

The possibility of obtaining an interspecific hybrid between tetraploid cultivated and diploid wild cotton through fertilization with an unreduced gamete raises the question of its evolution in natural populations. We have shown that this hybrid produced no seeds and its pollen has very low viability. In conclusion, the likelihood of natural transfer of chromosomal material from cultivated cotton to the diploid G. herbaceum is very low; even though one hybrid with a particular genomic structure arose via an unreduced gamete of the diploid parent, it did not produce seeds after selfing. However, the assessment of the hybrids' ability to produce a backcross progeny would require further experimentation, such as cultivating hybrids in the presence of G. herbaceum. Furthermore, the hybrid materials obtained in this study, with different genome dosages (triploids AAD and unreduced gamete-derived hybrid AAAD), could provide useful tools for the study of the expression of duplicated genes and the effects of allele dosage from A and D cotton genomes.

20 ACKNOWLEDGMENTS

The authors offer sincere thanks to Jonathan Wendel for reading the manuscript and valuable exchanges about their work. They acknowledge support for this project from the French

Ministry for Education and Research/ ANR (Agence Nationale de la Recherche), grant

MFD/AO OGM 0217. They are thankful to Mark Tepfer and Pierre Hilson for critical reading of the manuscript. The IJPB benefits from the support of the LabEx Saclay Plant Sciences-

SPS (ANR-10-LABX-0040-SPS).

REFERENCES

Ahoton L, Lacape JM, Baudoin JP, Mergeai G (2003) Introduction of Australian diploid cotton genetic variation into upland cotton. Crop Science 43:1999-2005

Alexander MP (1969) Differential staining of aborted and nonaborded pollen. Stain

Technology 44:117-122

Altman DW (1988) Exogenous hormone applications at pollination for in vitro and in vivo production of cotton interspecific hybrids. Plant Cell Reports 7:257-261

Amin KC (1940) Interspecific hybridization between Asiatic and New World cottons. Indian

J. Agric. Sci. 10 pp. 404- 413

Barcaccia G, Tavoletti S, Mariani A, Veronesi F (2003) Occurrence, inheritance and use of reproductive mutants in alfalfa improvement. Euphytica 133:37-56

21

Beasley JO (1940) The production of polyploids in Gossypium. J Hered 31: 39-48

Beasley JO (1942) Meiotic chromosome behavior in species, species hybrids, haploids and induced polyploids of Gossypium. Genetics 27: 25-54

Brownfield L, Köhler C (2011) Unreduced gamete formation in plants: mechanisms and prospects. Journal of Experimental Botany 62:1659-1668

Brubaker CL, Brown AHD, Stewart JMcD, Kilby MJ, Grace JP (1999) Production of fertile hybrid germplasm with diploid Australian Gossypium species for cotton improvement.

Euphytica 108:199-213

Carputo D, Frusciante L, Peloquin SJ (2003) The role of 2n gametes and endosperm balance number in the origin and evolution of polyploids in the tuber-bearing solanums. Genetics

163:287-294

Chan SWL (2010) Chromosome engineering: power tools for plant genetics. Trends in

Biotechnology 28, 12: 605-610

Chee PW, Paterson AH, Udall JA, Wendel JF (2016) Interspecific hybridization for upland cotton improvement. In Mason AS (Ed.), Polyploidy and interspecific hybridization for crop improvement CRC Press pp 1-20

22 Chèvre AM, Ammitzbøll H, Breckling B, Dietz-Pfeilstetter A, Eber F, Fargue A, Gomez-

Campo C, Jenczewski E, Jørgensen R, Lavigne C, Meier MS, Den Nijs H, Pascher K, Seguin-

Swartz G, Sweet J, Stewart Jr. CN, Warwick S (2004) A review on interspecific gene flow from oilseed rape to wild relatives. In: Introgression from Genetically Modified Plants into wild relatives. Den Nijs H.C.M, Bartsch D, Sweet J (Eds). CABI Publishing, Cambridge, pp

235-251

Crismani W, Girard C, Mercier R (2013) Tinkering with meiosis. Journal of Experimental

Botany 64:55-65

Dewitte A, Van Laere K, Van Huylenbroeck J (2012) Use of 2n Gametes in Plant Breeding.

Plant Breeding. Dr Ibrokhim Abdurakmonov (Ed.), InTech, pp 59-86

Eber F, Letanneur JC, Chèvre AM (1997) Chromosome number of oilseed rape (Brassica napus L.)-wild radish (Raphanus raphanistrum) spontaneous hybrids and of their progeny estimated by flow cytometry. Cruciferae Newsletter 19: 17-18

Fryxell, P.A. (1979) The Natural History of the Cotton Tribe. Texas A&M University Press,

College Station, TX, USA

Fryxell, P.A. (1992) A revised taxonomic interpretation of Gossypium L. ().

Rheedea 2, 2:108-165

Gallagher, J P.; Grover, C E.; Rex, K; Moran, M; Wendel, J F. (2017) A new species of cotton from wake atoll , Gossypium stephensii (Malvaceae). Systematic Botany, 42, 1, 115-

123

23 Gan YM, Liu F, Chen D, Wu Q, Qin Q, Wang CY, Li SH, Zhang XD, Wang YH, Kunbo

Wang XB (2013) Chromosomal Locations of 5S and 45S rDNA in Gossypium Genus and Its

Phylogenetic Implications Revealed by FISH. PLoS ONE 8(6): e68207

Gerlach WL, Bedbrook JR. (1979) Cloning and characterization of ribosomal RNA genes in wheat. Nucleic Acids Research 7:1869-1885

Gerlach WL and Dyer TA (1980) Sequence organisation of the repeating units in the nucleus of wheat which contain 5S rRNA genes. Nucl Acids Res. 8: 4851–4865

Grover CE, Gallagher JP, Jareczek JJ, Page JT, Udall JA, Gore MA, Wendel JF (2015) Re- evaluating the phylogeny of allopolyploid Gossypium L. Molecular Phylogenetics and

Evolution 92:45-52

Hanson RE, Islam-Faridi MN, Percival EA, Crane CF, Ji YF, McKnight TD, Stelly DM, Price

HJ (1996) Distribution of 5S and 18S-28S rDNA loci in a tetraploid cotton (Gossypium hirsutum L) and its putative diploid ancestors. Chromosoma 105:55-61

Heuberger S, Ellers-Kirk C, Tabashnik BE, Carriere Y (2010) Pollen- and Seed-Mediated

Transgene Flow in Commercial Cotton Seed Production Fields. PLoS ONE 5, 11: e14128

Kammacher P (1966) Etudes des relations génétiques et caryologiques entre génomes voisins du genre Gossypium. Coton et Fibres Tropicales 21: 263-290

24 Lacape JM, Dessauw D, Rajab M, Noyer JL, Hau B (2007) Microsatellite diversity in tetraploid Gossypium germplasm: assembling a highly informative genotyping set of cotton

SSRs. Molecular Breeding 19:45-58

Llewellyn D, Tyson C, Constable G, Duggan B, Beale S (2007) Containment of regulated genetically modified cotton in the field. Agr Ecosyst Environ 121: 419–429

Loureiro I, García-Ruiz E, Gutiérrez E, Gómez P, Escorial MC, Chueca MC (2016) Pollen- mediated gene flow in the cultivation of transgenic cotton under experimental field conditions in Spain. Industrial Crops and Products 85: 22–28

Mason AS, Pires JC (2015) Unreduced gametes: meiotic mishap or evolutionary mechanism?

Trends in Genetics 31, 1: 5-10

Mehetre SS, Aher AR, Gawande VL, Patil VR, Mokate AS (2003) Induced polyploidy in

Gossypium: A tool to overcome interspecific incompatibility of cultivated tetraploid and diploid cottons. Current Science 84:1510-1512

Mergeai, G (2004) Forty years of genetic improvement of cotton through interspecific hybridisation at Gembloux Agricultural University: Achievement and prospects. In: Proc.

World Cotton Research Conference-3. Cotton Production for the New Millennium, A

Swanepoel (Ed.). Cape Town, South Africa 9-13 March 2003, pp. 120–133

Mergeai G (2006) Cotton improvement through interspecific hybridization. Cahiers

Agricultures 15:135-143

25

Meyer VG (1974) Interspecific hydridization. Economic Botany 28: 56-60

Noormohammadi Z, Sheidai M, Shojaei F, Jeshvaghani FS, Alishah O (2012) Cytogenetic

Analysis of Mehr Cotton Cultivar and Its Crossing Progenies: A Search for Unreduced Pollen

Grains. Cytologia 77:107-112

Page JT, Liechty ZS, Alexander RH, Clemons K, Hulse-Kemp AM, Ashrafi H, Van Deynze

A, Stelly DM, Udall JA (2016) DNA Sequence Evolution and Rare Homoeologous

Conversion in Tetraploid Cotton. PloS Genetics 12, 5: e1006012

Peloquin SJ, Boiteux LS, Carputo D (1999) Meiotic mutants in potato: Valuable variants.

Genetics 153:1493-1499

Pereira GS, Sousa RL, Araújo RL, Hoffmann LV, Silva EF, Barroso PAV (2012) Selective fertilization in interspecific crosses of allotetraploid species of Gossypium. Botany 90:159-

166

Pierre J, Hofs JL (2010) Astylus atromaculatus (Coleoptera: Melyridae) : abundance and role in pollen dispersal in Bt and non-Bt cotton in South Africa. Environ Entomol (2010) 39 (5):

1523-1531

Ravi M, Chan SWL (2010) Haploid plants produced by centromere- mediated genome elimination. Nature 464: 615–618

Ravi M, Marimuthu MPA, Siddiqi I (2008) Gamete formation without meiosis in

26 Arabidopsis. Nature 451: 1121–1124

Ramanna MS, Jacobsen E (2003) Relevance of sexual polyploidization for crop improvement

- A review. Euphytica 133:3-18

Renny-Byfield S, Page JT, Udall JA, Sanders WS, Peterson DG, Arick II MA, Grover CE,

Wendel JF (2016) Independent domestication of two Old World cotton species. Genome

Biology and Evolution 8, 7: 1940-1947

Sheidai M (2008) Cytogenetic distinctiveness of sixty-six tetraploid cotton (Gossypium hirsutum L.) cultivars based on meiotic data. Acta Botanica Croatica 67:209-220

Skovsted A (1934) Cytological studies in Cotton. II. Two interspecific hybrids between

Asiatic and new world cottons. J. Genet. 28: 407–424

Stephens SG (1942) Colchicine-produced polyploids in Gossypium I. And autotetraploid asiatic cotton and certain of its hybrids with wild diploid species. J. Genet. 44: 272-295

Stephens SG (1974) Geographic and taxonomic distribution of anthocyanin genes in new world cottons. Journal series of the North Carolina State University Agricultural Expimental

Station, Raleigh, North Carolina. Paper 4061: 128-141

Stewart JMcD (1995) Potential for crop improvement with exotic germsplams and genetic engineering. In Challenging the Future - Proc. World. Cotton Research Conf. -1; Brisbane,

Australia 14-17 February 1994. Constable GA and Forrester NW (Eds). CSIRO, Melbourne pp 313-327

27 Van Deynze AE, Sundstrom FJ, Bradford KJ (2005) Pollen-mediated gene flow in California cotton depends on pollinator activity. Crop Sci 45: 1565–1570

van Wyk B and van Wyk P (1997). Field guide to trees of southern Africa. Struik Publishers,

Cape Town, South Africa. 536 pp

Warwick SI, Simard MJ, Legere A, Beckie HJ, Braun L, Zhu B, Mason P, Seguin-Swartz G,

Stewart CN (2003) Hybridization between transgenic Brassica napus L. and its wild relatives:

Brassica rapa L., Raphanus raphanistrum L., Sinapis arvensis L., and Erucastrum gallicum

(Willd.) O.E. Schulz. Theor Appl Genet 107:528–539

Wendel JF (1989) New world tetraploid cottons contained old world cytoplasm. Proc. Natl.

Acad. Sci. USA 86:4132-4136

Wendel JF, Cronn RC (2003) Polyploidy and the evolutionary history of cotton. Advances in

Agronomy 78:139-186

Wendel JF, Brubaker CL, Seelanan T (2010) The origin and evolution of Gossypium. In

JMcD Stewart, JM Oosterhuis, JJ Heitholt and JR Mauney (Eds.), Physiology of Cotton,

Springer, Netherlands pp1-18

Younis A, Hwang YJ, Lim KB (2014) Exploitation of induced 2n-gametes for plant breeding.

Plant Cell Reports 33:2015-223

28 Zhang BH, Pan XP, Guo TL, Wang QL, Anderson TA (2005) Measuring gene flow in the cultivation of transgenic cotton (Gossypium hirsutum L.). Mol Biotechnol 31:11–20

29 1

30