In Silico and Fluorescence In Situ Hybridization Mapping Reveals Collinearity between the Pennisetum squamulatum Apomixis Carrier-Chromosome and Chromosome 2 of Sorghum and Foxtail Millet

Sirjan Sapkota; Joann A. Conner; Wayne W. Hanna; Bindu Simon; Kevin Fengler; Stéphane Deschamps; Mark Cigan; Peggy Ozias-Akins

doi:10.1371/journal.pone.0152411

Abstract

Apomixis, or clonal propagation through seed, is a trait identified within multiple species of the grass family (Poaceae). The genetic locus controlling apomixis in Pennisetum squamulatum (syn Cenchrus squamulatus) and Cenchrus ciliaris (syn Pennisetum ciliare, buffelgrass) is the apospory-specific genomic region (ASGR). Previously, the ASGR was shown to be highly conserved but inverted in marker order between P. squamulatum and C. ciliaris based on fluorescence in situ hybridization (FISH) and varied in both karyotype and position of the ASGR on the ASGR-carrier chromosome among other apomictic Cenchrus/Pennisetum species. Using in silico transcript mapping and verification of physical positions of some of the transcripts via FISH, we discovered that the ASGR-carrier chromosome from P. squamulatum is collinear with chromosome 2 of foxtail millet and sorghum outside of the ASGR. The in silico ordering of the ASGR-carrier chromosome markers, previously unmapped in P. squamulatum, allowed for the identification of a backcross line with structural changes to the P. squamulatum ASGR-carrier chromosome derived from gamma irradiated pollen.

Citation: Sapkota S, Conner JA, Hanna WW, Simon B, Fengler K, Deschamps S, et al. (2016) In Silico and Fluorescence In Situ Hybridization Mapping Reveals Collinearity between the Pennisetum squamulatum Apomixis Carrier-Chromosome and Chromosome 2 of Sorghum and Foxtail Millet. PLoS ONE 11(3): e0152411. https://doi.org/10.1371/journal.pone.0152411

Editor: Emidio Albertini, University of Perugia, ITALY

Received: January 25, 2016; Accepted: March 14, 2016; Published: March 31, 2016

Copyright: © 2016 Sapkota et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Data Availability: All relevant data are within the paper and its Supporting Information files.

Funding: This study was partially funded by a Sponsored Research Agreement between Peggy Ozias-Akins/UGA and Pioneer Hi-Bred. Pioneer Hi-Bred provided financial support for research materials, in kind support, and had intellectual input into the design, analysis and interpretation of the results.

Competing interests: The financial support of Pioneer Hi-Bred does not alter the authors' adherence to PLOS ONE policies on sharing data and materials.

Introduction

Most current grass species are found to be derived from a common ancestor that lived about 50–80 million years ago. Despite the relatively recent and monophyletic origin of the grass genomes, there is considerable divergence in genome size and chromosome number [1, 2]. Yet sequenced members of the Poaceae clade have shown conservation of gene order (collinearity) among species such as rice, sorghum, maize, and foxtail millet [3]. Pearl millet (Pennisetum glaucum) shares a common ancestor with foxtail millet (Setaria italica) ~8.3 million years ago and with maize and sorghum ~26 million years ago [4].

There are several apomictic species in the genus Cenchrus/Pennisetum. Apomixis has been defined as asexual reproduction through seed [5]. Apomixis in P. squamulatum (syn Cenchrus squamulatus) and the closely related species C. ciliaris (syn Pennisetum ciliare) was found to be controlled by a dominant and hemizygous genetic locus named the apospory-specific genomic region (ASGR) [6, 7]. The ASGR in P. squamulatum is identified as a large (>50 Mb in size), heterochromatic chromosomal block localized near the telomere on the short arm of the ASGR-carrier chromosome by fluorescence in situ hybridization (FISH). The ASGR contains a region of low copy DNA flanked by regions of high copy DNA [8]. Physical mapping of the ASGR between P. squamulatum and C. ciliaris identified an inversion but conservation of ASGR-BAC order between the two species [9]. Comparison of the ASGR across other apomictic Cenchrus/Pennisetum species revealed variation in the ASGR chromosomal position and ASGR-carrier chromosome structure, suggesting movement and divergence of the ASGR block within the Cenchrus/Pennisetum lineage [10]. Within the ASGR, partial sequencing of ASGR-linked BAC clones showed the presence of multiple regions of small-scale, but not large-scale, collinearity with the rice and sorghum genomes [11].

An apomictic backcross (06-A-58) of P. glaucum, originating from a cross between P. squamulatum and tetraploid P. glaucum, was identified by FISH to carry one alien chromosome, the ASGR-carrier chromosome from P. squamulatum [12]. Forty-nine contigs, generated from the assembly of 454 sequences derived from dissected apomictic ovules, were mapped to the ASGR-carrier chromosome via SCAR, CAPS or SSCP markers [13]. Contig PS26_c9369 demonstrated tight linkage to the ASGR. Three other contigs (PS26_c5080, PS26_c33813, and PS26_c2552) mapped as unlinked to the ASGR using CAPS markers. The remaining contig SCAR markers failed to identify polymorphisms for mapping within a segregating F₁ population. However, in silico mapping of the 49 ASGR-carrier chromosome contigs to the sorghum reference genome identified 21 sequences (BlastN, e value ≤e ^-20) with similarity to the sorghum genome of which 17 had a unique or highest similarity to sorghum chromosome 2.

As the 454 transcriptome data was 3ʹ biased due to T7 amplification of the ovule RNA, we generated additional transcriptome assemblies using RNA-seq data from an apomictic backcross and screened these to extend the length of the 454 ASGR-carrier chromosome contigs for additional comparison to the sorghum and foxtail millet genomes. The predicted in silico positions of 7 contigs from the P. squamulatum ASGR-carrier chromosome, based on hits to the sorghum and foxtail millet genome, were verified by cytogenetic mapping of BAC clones containing the SCAR marker for the ASGR-carrier chromosome contigs. Our analysis demonstrates that the ASGR-carrier chromosome from P. squamulatum, outside the ASGR boundary, is collinear with chromosome 2 of sorghum and foxtail millet. Using the established ASGR-carrier chromosome SCAR markers and physical mapping results, a screen of gamma irradiated offspring was tested to identify lines with structural changes to the ASGR-carrier chromosome. This screen identified a sexual line which has lost the ASGR, but retained most of the long arm of the ASGR-carrier chromosome.

Materials and Methods

Plant material

Apomictic P. glaucum backcross 8 (BC₈) (06-A-58) derived seedlings carrying the ASGR-carrier chromosome [12] were used for FISH and tissue collection for RNA extraction. Backcross 06-A-58 is a facultative apomict and therefore produces progeny derived through both modes of reproduction, apomictic and sexual. DNA from seedlings was extracted using a modified CTAB method [14] and screened for the ASGR using the ASGR-linked SCAR marker p787/788 [12].

Extension of ASGR-carrier chromosome transcript information

RNA was extracted from unfertilized ovaries of apomictic and sexual BC₈ (06-A-58) derived plants collected on the day of anthesis (anther exsertion). Approximately 50 ovaries were collected from individual plants based on their mode of reproduction. RNA was isolated from 12 apomictic and 12 sexual plants usingQiagen Plant RNAeasy kit. RNA from individuals was pooled based on mode of reproduction to make an apomictic and sexual RNA sample for sequencing. Sequencing libraries were constructed according to manufacturer’s instructions and Illumina sequenced to yield 2 x 76 paired-end reads. Sequences from both libraries were quality trimmed, separated into paired-end and single-end reads (~100 million reads for each library), and assembled together with the Velvet de novo [15] assembly algorithm. Two different assemblies were used (varying in K-mer value) to identify the longest ASGR-carrier chromosome transcript available. PS contig sequences used for in silico analysis are located in S1 File.

Identification of ASGR-carrier chromosome BACs

BACs linked to the ASGR-carrier chromosome were identified by screening the polyhaploid BAC library [16] with probes derived from the ASGR-carrier chromosome transcripts [13]. Hybridizing BAC clones were confirmed as linked to the ASGR-carrier chromosome via PCR amplification with the respective ASGR-carrier chromosome SCAR marker. Markers with isolated BACs are indicated in Table 1.

Download:

Table 1. Information for in silico and deletion mapping line.

https://doi.org/10.1371/journal.pone.0152411.t001

Southern blot hybridization of HindIII digested fragments of BAC clones with ³²P label genomic DNA from apomictic BC₈ (06-A-58) was used to assess the level of DNA repetitiveness within the BAC clones based on signal intensity. A centromeric BAC clone was identified from the polyhaploid BAC library using a 160 bp KpnI repeat probe [17].

FISH

FISH probes.

BAC DNA for nick translation was extracted using an alkaline lysis method http://www.protocolpedia.com/component/sobipro/?pid=69&sid=2209:BAC-DNA-Isolation-from-200-ml-Cultures-by-a-Cleared-Lysate-Method-Followed-by-Double-Acetate-Precipitation&Itemid=0 with the following modifications. The RNase treatment was done with 10 μl Ambion® RNase cocktail (Life Technologies, Grand Island, NY) consisting of 5 U of RNase A and 200 U of RNase T1. The BAC DNA was suspended in a final volume of 50 μl Buffer EB (QIAGEN Inc., Valencia, CA, USA). 1–2 μg of BAC DNA was labeled with biotin (bio)-11-dUTP (Roche, Indianapolis, IN) or digoxigenin (dig)-11-dUTP (Roche), using the nick translation kit (Roche) according to manufacturer’s instructions.

PCR centromere probes were prepared by labeling with biotin-11-dUTP using primers 5’-GGTACCCCGTAATAGTGCATTC-3’ and 5’-GAAAATGGTTTCGCAACAAAAG-3’ designed from the 160 bp KpnI repeat family sequence [17].

Chromosome preparation.

Root tips from apomictic BC₈ (06-A-58) derived seedlings were collected, washed, placed in a 0.5ml Eppendorf tube with a hole in the lid in 300 μl distilled water and treated with nitrous oxide at 1 to 1.5 Mpa for 3 to 4 hours at room temperature in a Nitrous Oxide gas chamber [18] prior to fixing in 3:1 (V:V) ethanol to acetic acid solution. Root caps were removed and 2–3 mm of the meristematic region was incubated in an enzyme mix containing 2% (w/v) cellulose RS (Karlan Research, Santa Rosa, CA), 1% (w/v), pectolyase Y23 (Karlan Research, Santa Rosa, CA), 1% (w/v) macerozyme R 10 (Desert Biologicals, Phoenix, AZ) in citrate buffer (10 mM sodium citrate, 10 mM sodium EDTA, pH 5.5) [19] for 90 minutes at 37°C. Slide preparation for chromosome spreads after digestion was done either through air-drying [20] or a “SteamDrop” method [21].

Fluorescence in situ hybridization and detection.

FISH was performed according to Zhong [22] with modifications. Slides with chromosomal spreads were treated with 5 μg/ml pepsin in 0.01M HCl for 5 to 10 minutes, fixed in 1% formaldehyde with 50 mM MgCl₂ in 2x SSC and dehydrated in a series of 70, 90 and 100% (v/v) ethanol. The hybridization mix consisted of 1–5 ng/μl of each probe, 50% formamide, 10% dextran sulfate, 75–85 ng/μl salmon sperm DNA and 2× SSC in a final volume of 18–20 μl. If necessary, P. squamulatum blocking DNA (10–50 ng/μl) was added to the hybridization mix to block signal from minor repetitive sequence within BAC clones. Hybridization mixtures were denatured at 80°C for 5 minutes, snap chilled on ice, applied to the chromosome spread, and covered with a 22 × 30 mm coverslip. Slides were placed on an 80°C heat block for 3 minutes then incubated in a moist chamber at 60°C for 90 minutes followed by a 37°C incubation for 64–67 hours. Two post-hybridization washes were done in 50% formamide in 2× SSC at 37°C for 10 minutes each. Slides were blocked in TNB (100 mM Tris-HCl, pH 7.5; 150 mM NaCl; 1× DIG blocking solution, Roche) for 30 minutes at 37°C, and blocked again in 5% (w/v) IgG-free bovine serum albumin (Sigma) in TN (100 mM Tris-HCl, pH 7.5; 150 mM NaCl) for 30 minutes at 37°C. Two-color detection was carried out according to Zhong et al. (1996) with modifications. The biotin-labeled probes were detected with Texas red using a three step amplification and DIG-labeled probes were detected with FITC with a two-step amplification. Antibodies were diluted in TNB. Preparations were counterstained by mounting in Vectashield (Vector Laboratories) containing 1.5 μg/ml DAPI. Slides were examined with a Zeiss Axioskop 2 plus fluorescence microscope. Fluorescent signals were detected for DAPI (λ_ex = 360 nm, λ_em = 420 nm), FITC (λ_ex = 480 nm, λ_em = 515 nm), and Texas red (λ_ex = 560 nm, λ_em = 645 nm). Monochrome digital images were captured using a charge-coupled device AxioCam camera and stored using AxioVision, version 4.8 for Windows. Composite images were constructed using Adobe® Photoshop CS2 version 9.0.

Deletion Study

Pollen irradiation.

Pollen was collected in glassine bags from individual heads each morning between 10:30 and 11:00 am from greenhouse-grown plants. Individual plants were derived from six backcross 8 lines and one backcross 7 line, genotyped as ASGR positive. Pollen was irradiated within 30 minutes of collection in the glassine bags using a J.L. Shepherd Model 109-GR-12 self-contained Cobalt-60 irradiator set for either 2 or 3 Kr of Cobalt-60 gamma radiation. After irradiation, the pollen was immediately used to pollinate inflorescences of sexual tetraploid pearl millet which had stigmas exserted, but at least a day before pollen shed.

Molecular screen of plant lines derived from irradiated pollen crosses.

DNA was extracted from equal amounts of tissue from 4 green plants within a line using a modified CTAB method [14]. DNA was then amplified with ASGR specific primers 787/788 and ASGR-carrier chromosome CAPS marker p1670/71 and SCAR marker 1656/1657 [13].

Results and Discussion

Longer sequences were identified for 44 of the 49 ASGR-carrier chromosome 454 contigs within the Velvet assemblies (Table 1, S1 File). The Velvet-assembled contigs had to share at least 90% or greater sequence identity to the ASGR-carrier chromosome 454 contigs [13]. The longest available ASGR-carrier chromosome contig was used for BlastN (cutoff of e^-20) in silico analysis against the foxtail millet (NW_004675962.1) and sorghum (NC_012877.1) genomes at the National Center for Biotechnology Information (NCBI). With the additional contig lengths, 47 and 44 of the ASGR-carrier chromosome contigs showed similarity to the foxtail millet and sorghum genomes, respectively. Forty-three (88%) and 37 (82%) of the ASGR-carrier chromosome contigs showed a unique or highest similarity to foxtail millet and sorghum chromosome 2 with individual contig hits distributed along the length of chromosome 2 of both species (Table 1). Contig PS26_c9369, tightly linked to the ASGR [13], had similarity to chromosome 1 in foxtail millet and chromosomes 3, 4, 6, and 10 in sorghum. PS_c194, PS_c28392 and PS_c2807 had similarity to foxtail millet chromosome 2 but did not have corresponding BlastN hits to the sorghum genome. PS_c9776, PS_c9993 and PS_c13655 had hits to foxtail millet chromosome 2 but identified more significant similarity to genes on sorghum chromosomes 6 or 7. These noted PS contigs did not tightly cluster in a particular area of the foxtail millet chromosome. PS26_c283, aligned on sorghum chromosome 2 but to foxtail millet chromosome 9, although it was identified on a scaffold mapped to chromosome 2 in the Beijing Genomics Institute (BGI) foxtail millet genome assembly (http://foxtailmillet.genomics.org.cn). The identification of large-scale collinearity between sorghum and foxtail millet for chromosome 2 was expected based on whole-genome dot plot comparisons which show that chromosome 2 in sorghum and foxtail millet share large degrees of similarity except at the centromeric region [4]. While the pearl millet genome is not yet available, comparative mapping revealed that pearl millet linkage group 7 is homoeolgous to foxtail millet 2 [23] and is likely the homoeolgous chromosome for the P. squamulatum ASGR-carrier chromosome.

A range of one to five BAC clones were isolated from the polyhaploid BAC library [15] for 25 of the 49 ASGR-carrier chromosome transcripts (Table 1). The relative amount of repetitive DNA within each BAC clone was assayed by the signal strength and number of restriction fragments of the BAC DNA hybridizing to labeled 06-A-58 total DNA when compared to ASGR BAC clones p109 and p800 which were used as low and high copy controls, respectively [8]. BAC clones with moderately (m) or highly (h) repetitive DNA are noted in Table 1. The ASGR-carrier chromosome BAC clones were selected for physical mapping based on the in silico mapping of the contig to the foxtail millet and sorghum genomes and their level of repetitive DNA. Seven BAC clones were physically mapped to the ASGR-carrier chromosome in apomictic 06-A-58. The ASGR on the ASGR-carrier chromosome was detected using either a high copy ASGR-BAC (red pseudo color signal denoted by the yellow ASGR arrow) which hybridizes to the high copy regions of the ASGR [8] (Fig 1B, 1C and 1F) or a combination of a high and low copy ASGR-BAC (green pseudo color signal denoted by the yellow ASGR arrow Fig 1A, 1D, 1E and 1G) near the telomere of the ASGR-carrier chromosome. A centromere probe was used in most FISH experiments (red pseudo color signal denoted by the red arrow). The position of the ASGR-carrier chromosome BACs are green pseudo color and denoted by a green arrow (Fig 1A–1G). Two of the seven ASGR-carrier chromosome BAC clones (Ps26_c2838 and Ps26_c583) mapped to the short arm of the ASGR-carrier chromosome while the remaining BAC clones mapped to the long arm. This result agrees with the in silico prediction based on centromere positions on chromosome 2 of sorghum and foxtail millet at ~30–35 Mb and ~17–20 Mb, respectively (http://ensembl.gramene.org). Using dual labeling of BAC probes, the order of PS_c3993/PS_c10535, PS26_5080 and PS26_c5851 was verified as linear from the centromere to telomere on the long arm of the ASGR-carrier chromosome. The signal order of PS_c3993 and PS_c10535 to each other was not verified. Our attempt to physically map the PS26_c1279 BAC, which in silico maps closer to the centromere on the long arm of the chromosome, was unsuccessful, even with blocking, due to repetitive DNA, which hybridized as a large signal on both sides of the centromere on the ASGR carrier chromosome.

Download:

Fig 1. FISH images of BAC probes mapped to the ASGR-carrier chromosome.

Example of a mitotic chromosome spread of apomictic BC₈ (06-A-58) with the ASGR identified with a yellow arrow and the centromere of the ASGR-carrier chromosome identified by a red arrow. A-G insets are mapped ASGR-carrier chromosome BACs. The yellow arrow denotes the ASGR signal. The ASGR signal was identified either through a high copy ASGR-BAC clone (red pseudo-color) or a combination of a high (red pseudo-color) and low (green pseudo-color) ASGR-BAC clone. The red arrow denotes the centromere signal and the green arrow denotes the mapped ASGR carrier chromosome BAC signal from the following a) p285J18/PS_c2838, b) p220A02/PS_c583 and p236E19/PS_c10535 (*), c) p036L06/PS_c2448, d) p258L05/PS_c3993, e) p236E19/PS_c10535, f) p057M05/PS_c5080, and g) p142D19/PS_c5851.

https://doi.org/10.1371/journal.pone.0152411.g001

Positions determined with cytogenetic mapping along the ASGR-carrier chromosome corresponded to the predicted in silico locations based on the sorghum/foxtail millet genomes as shown in Figs 1 and 2. Based on the in silico and cytogenetic mapping data, the collinearity of the P. squamulatum ASGR-carrier chromosome outside the ASGR and foxtail millet starts from ~5.5 Mb of foxtail millet and continues across the entire chromosome. Whether the synteny extends to the very beginning of chromosome 2 to right outside the ASGR remains unknown until additional BACs can be identified. Of the seven P. squamulatum ASGR-carrier chromosome BACs used for physical mapping, none showed additional FISH signals on the pearl millet chromosomes. The lack of signal on the pearl millet genome could be due to FISH hybridization conditions or to the divergence of the non-genic regions between the P. glaucum and P. squamulatum chromosomes. Additional research will be required to determine if the ASGR-carrier chromosome markers will be useful to identify similar regions in other apomictic Pennisetum species. If the ASGR-carrier chromosome BACs show signal on other Pennisetum species, this may allow us to more fully understand the evolution of the apomixis locus within the Pennisetum species.

Download:

Fig 2. Diagram of collinearity identified using FISH signals on the P. squamulatum ASGR-carrier chromosome and the in silico positions in sorghum and foxtail millet chromosome 2.

Chromosome lengths are not drawn to scale. Black box denotes BAC clones where linear order was confirmed. *Order of Ps26_c3993 and Ps26_c10535 BACs to each other was not verified.

https://doi.org/10.1371/journal.pone.0152411.g002

A published RFLP mapping study placed the apospory locus on linkage group 7b in buffelgrass (C. ciliaris) [24]. Sequence information could be found for 3 RFLP markers within the 6 marker linkage group covering 78.8 cM. Two markers, HHU27 (gb|H54993.1) and pPAP3A07 (gb|BM084123.1) flanked the aposporous locus at 10.7 cM and 1.4 cM, respectively. The third marker, pPAP08H05 (gb|BM084577.1) resided 43.2 cM from HHU27. All 3 markers showed highest similarity in a linear order to chromosome 7 in foxtail millet and chromosome 6 in sorghum. In foxtail millet, pPAP3A07 was located at ~14 Mb, HHU27 at ~24 Mb and pPAP08H05 at ~31 Mb. In sorghum, pPAP3A07 was located at ~19.1 Mb, HHU27 at ~50.8 Mb and pPAP08H05 at ~58.4 Mb. While limited in data points, the marker sequence comparison and karyotype differences identified between the ASGR-carrier chromosomes in C. ciliaris and P. squamulatum [9, 10] further supports the idea of translocation of the ASGR to different chromosomes between apomictic Pennisetum/Cenchrus species. Changes at the apomixis-controlling-locus (ACL) in the Paspalum genus have also been identified. Markers from the telomeric portion of the long arm of rice chromosome 12 flank the ACL in mapping studies with P.simplex [25] and P. malacophyllum [26]. However, P. notatum showed both rice chromosome 2 and chromosome 12 markers flanking the ACL [26, 27]. Markers identified as ACL linked in P. simplex were not linked to the ACL in P. procurrens [28]. Sorghum chromosome 2 shows synteny blocks with rice chromosome 3, 7, 8 and 9, but not with rice chromosome 2 or 12 [29].

The DNA content of the P. squamulatum chromosome was roughly estimated at ~200 Mbp [8] of which ~50 Mbp is the ASGR. Therefore the P. squamulatum ASGR-carrier chromosome has expanded roughly 2 to 3 times when compared to the corresponding sorghum and foxtail millet chromosome 2. It is likely that much of that expansion is caused by transposable elements as has been shown when comparing the number of predicted genes within a genome to genome size in many plant species [30].

With the location of the ASGR-carrier chromosome markers identified, a screen to identify backcross lines with potential structural changes to the ASGR-carrier chromosome was undertaken. We sought to identify a line where a functional ASGR locus had been moved from the P. squamulatum ASGR-carrier chromosome to a pearl millet chromosome via irradiated pollen (S1 Fig). For the study, irradiated pollen from offspring derived from six different apomictic backcross 8 and one apomictic backcross 7 tetraploid pearl millet lines was used to pollinate 71 sexual tetraploid IA4X heads. As the ASGR is a single-dose locus, approximately half of the pollen used in the crosses would not carry the ASGR. To help identify lines which reproduced sexually, heads from 1962 plants derived from the irradiated pollen by IA4X cross were pollinated with Red IA4X pollen and seed collected. Red IA4X pollen contains the dominant Rp1 allele which confers a red midrib color [31]. Plants producing only red progeny, indicating obligate sexuality, were not screened by DNA markers. Seventy-eight lines were initially tested for structural changes to the ASGR-carrier chromosome using ASGR SCAR marker p787/788 and ASGR-carrier chromosome SCAR markers from PS_c6373 (in silico mapped approximately half way between the centromere and telomere on the long arm of foxtail millet and sorghum) and CAPS marker from PS26_c2552 (in silico mapped close to the telomere of the long arm in foxtail millet and sorghum) based on their production of green, and therefore potentially apomictic, progeny. The markers chosen for screening are co-dominant and therefore would eliminate false negative results in the PCR- based screen. Fifty-one lines did not carry either the ASGR or ASGR-carrier chromosome markers tested. These lines producing green progeny were generated by self-pollination after unsuccessful crossing with Red IA4X pollen. Twenty-six lines carried all three markers and therefore did not contain large structural changes to the ASGR-carrier chromosome. One line, 312, contained both ASGR-carrier chromosome markers but not the ASGR marker. Line 312 was derived from 2 Kr irradiated pollen from the BC₇ line. Plant 312 was screened for reproductive phenotype by ovary clearing [32]. As expected, plant 312 formed mature sexual embryo sacs. Further screening using the ASGR-carrier chromosome SCAR markers showed that the 312 plant inherited most of the long arm of the ASGR-carrier chromosome from P. squamulatum (Table 1).

As shown in our preliminary screen of gamma irradiated apomictic pollen offspring, the mapped ASGR-carrier chromosome markers can be used to identify structural changes in the ASGR-carrier chromosome as found in plant 312. Additional screening of more plants may allow us to identify apomictic lines with large deletions in the P. squamulatum ASGR-carrier chromosome. If found, these plants could be subjected to both genomic and transcriptional sequencing which could help as a process of elimination to identify the genes controlling apomixis in P. squamulatum.

Supporting Information

S1 Fig. Graphical Overview of Deletion Study Screen.

Graphical overview detailing the steps of the deletion study screen.

https://doi.org/10.1371/journal.pone.0152411.s001

(TIF)

S1 File. PS_Contigs.txt.

Fasta file of PS_contigs.

https://doi.org/10.1371/journal.pone.0152411.s002

(TXT)

Acknowledgments

We would like to thank Tracey Vellidis for help in generating the figures.

Author Contributions

Conceived and designed the experiments: JAC PO-A WWH KF SD MC. Performed the experiments: SS JAC WWH BS KF SD. Analyzed the data: SS JAC WWH SD KF. Wrote the paper: JAC SS PO-A.

References

1. Bennetzen JL. Patterns in grass genome evolution. Curr Opin Plant Biol. 2007; 10: 176–181. pmid:17291821
- View Article
- PubMed/NCBI
- Google Scholar
2. Gaut BS. Evolutionary dynamics of grass genomes. New Phytol. 2002; 154: 15–28.
- View Article
- Google Scholar
3. Wang X, Wang J, Jin D, Guo H, Lee T-H, Liu T, et al. (2015) Genome alignment spanning major Poaceae lineages reveals heterogeneous evolutionary rates and alters inferred dates for key evolutionary events, Molec Plant 2015; 8:885–898. http://dx.doi.org/10.1016/j.molp.2015.04.004.
- View Article
- Google Scholar
4. Bennetzen JL, Schmutz J, Wang H, Percifield R, Hawkins J, Pontaroli AC, et al. Reference genome sequence of the model plant Setaria. Nature Biotech. 2012; 30: 555–561.
- View Article
- Google Scholar
5. Nogler GA. Gametophytic apomixis. In: Johri BM (ed) Embryology of Angiosperms, Springer Verlag, pp 475–518; 1984.
6. Ozias-Akins P, Roche D, Hanna WW. 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. Proc Natl Acad Sci U S A. 1998; 95: 5127–5132. pmid:9560240
- View Article
- PubMed/NCBI
- Google Scholar
7. Roche D, Cong P, Chen Z, Hanna WW, Gustine DL, Sherwood RT, et al. An apospory-specific genomic region is conserved between Buffelgrass (Cenchrus ciliaris L.) and Pennisetum squamulatum Fresen. Plant J. 1999; 19: 203–208. pmid:10476067
- View Article
- PubMed/NCBI
- Google Scholar
8. Akiyama Y, Conner JA, Goel S, Morishige DT, Mullet JE, Hanna WW, et al. High-resolution physical mapping in Pennisetum squamulatum reveals extensive chromosomal heteromorphism of the genomic region associated with apomixis. Plant Physiol. 2004; 134: 1733–1741. pmid:15064383
- View Article
- PubMed/NCBI
- Google Scholar
9. Goel S, Chen Z, Akiyama Y, Conner JA, Basu M, Gualtieri G, et al. Comparative physical mapping of the apospory-specific genomic region in two apomictic grasses: Pennisetum squamulatum and Cenchrus ciliaris. Genetics. 2006; 173: 389–400. pmid:16547108
- View Article
- PubMed/NCBI
- Google Scholar
10. Akiyama Y, Goel S, Conner JA, Hanna WW, Yamada-Akiyama H, Ozias-Akins P. Evolution of the apomixis transmitting chromosome in Pennisetum. BMC Evol Biol. 2011; 11: 289–289. pmid:21975191
- View Article
- PubMed/NCBI
- Google Scholar
11. Conner JA, Goel S, Gunawan G, Cordonnier-Pratt MM, Johnson VE, Liang C, et al. Sequence analysis of bacterial artificial chromosome clones from the apospory-specific genomic region of Pennisetum and Cenchrus. Plant Physiol. 2008; 147: 1396–1411. pmid:18508959
- View Article
- PubMed/NCBI
- Google Scholar
12. Singh M, Conner JA, Zeng YJ, Hanna WW, Johnson VE, Ozias-Akins P. Characterization of apomictic BC₇ and BC₈ pearl millet: meiotic chromosome behavior and construction of an ASGR-carrier chromosome-specific library. Crop Sci. 2010; 50: 892–902.
- View Article
- Google Scholar
13. Zeng Y, Conner J, Ozias-Akins P. Identification of ovule transcripts from the Apospory-Specific Genomic Region (ASGR)-carrier chromosome. BMC Genomics. 2011; 12: 206–206. pmid:21521529
- View Article
- PubMed/NCBI
- Google Scholar
14. Conner JA, Gunawan G, Ozias-Akins P. Recombination within the apospory specific genomic region leads to the uncoupling of apomixis components in Cenchrus ciliaris. Planta. 2013; 238: 51–63. pmid:23553451
- View Article
- PubMed/NCBI
- Google Scholar
15. Zerbino DR. Using the Velvet de novo assembler for short-read sequencing technologies. Current Protocols In Bioinformatics / Editoral Board, Andreas D. Baxevanis … [et al.] Chapter 11: 11.15–11.15; 2010.
16. Roche D, Conner JA, Budiman A, Frisch D, Wing R, Hanna W, Ozias-Akins P. Construction of BAC libraries from two apomictic grasses to study the microcolinearity of their apospory-specific genomic regions. Theor Appl Genet. 2002; 104: 804–812. pmid:12582640
- View Article
- PubMed/NCBI
- Google Scholar
17. Ingham LD, Hanna WW, Baier JW, Hannah LC. Origin of the main class of repetitive DNA within selected Pennisetum species. Mol Gen Genet. 1993; 238: 350–356. pmid:8492802
- View Article
- PubMed/NCBI
- Google Scholar
18. Andres RJ and Kuraparthy V. Development of an improved method of mitotic metaphase chromosome preparation compatible for fluorescence in situ hybridization in cotton. J Cotton Sci. 2013; 17:149–156.
- View Article
- Google Scholar
19. Iwata A, Tek AL, Richard MM, Abernathy B, Fonsêca A, Schmutz J, et al. Identification and characterization of functional centromeres of the common bean. Plant J. 2013; 76: 47–60. pmid:23795942
- View Article
- PubMed/NCBI
- Google Scholar
20. Fukui K. Plant chromosome at mitosis. In: Fukui K and Nagayama S (eds) Plant chromosomes laboratory methods, edited by FUKUI K. and NAGAYAMA S.. CRC Press, Boca Raton, FL, pp 1–17; 1996.
21. Kirov I, Divashuk M, Van Laere K, Soloviev A, Khrustaleva L. An easy "SteamDrop" method for high quality plant chromosome preparation. Mol Cytogenet. 2014; 7: 21 pmid:24602284
- View Article
- PubMed/NCBI
- Google Scholar
22. Zhong X, Fransz PF, van Eden JW, Zabel P, van Kammen Ab, de Jong JH. High-resolution mapping on pachytene chromosomes and extended DNA fibres by fluorescence in-situ hybridisation. Plant Molec Biol Rep. 1996; 14: 232–242.
- View Article
- Google Scholar
23. Devos KM, Pittaway TS, Reynolds A, Gale MD. Comparative mapping reveals complex relationship between the pearl millet genome and those of foxtail millet and rice. Theor Appl Genet. 2000; 100:190–198.
- View Article
- Google Scholar
24. Jessup RW, Burson BL, Burow GB, Wang YWCC, Li Z, Paterson AH, Hussey MA. Disomic inheritance, suppressed recombination, and allelic interactions govern apospory in buffelgrass as revealed by genome mapping. Crop Sci. 2002; 42: 1688–1694.
- View Article
- Google Scholar
25. Pupilli F, Labombarda P, Caceres ME, Quarin CL, Arcioni S. The chromosome segment related to apomixis in Paspalum simplex is homoeologous to the telomeric region of the long arm of rice chromosome 12. Mol Breeding 2001, 8:53–61.
- View Article
- Google Scholar
26. Pupilli F, Martinez EJ, Busti A, Calderini O, Quarin CL, Arcioni S. Comparative mapping reveals partial conservation of synteny at the apomixis locus in Paspalum spp. Molecular Genetics and Genomics 2004, 270:539–548. pmid:14648202
- View Article
- PubMed/NCBI
- Google Scholar
27. Podio M. Rodriuez MP, Felitti S, Stein J, Martinez EJ, Siena LA, et al. Sequence characterization, in silico mapping and cytosine methylation analysis of markers linked to apospory in Paspalum notatum Genetics and Molecular Biology 2012, 35 (4): 827–837. pmid:23271945
- View Article
- PubMed/NCBI
- Google Scholar
28. Hojsgaard DH, Martinez EJ, Acuna CA, Quarin CL, Pupilli F. A molecular map of the apomixis-control locus in Paspalum procurrens and its comparative analysis with other species of Paspalum. Theor Appl Genet 2011, 123:959–971. pmid:21713535
- View Article
- PubMed/NCBI
- Google Scholar
29. Paterson AH, Bowers JE, Bruggmann R, Dubchak I, Grimwood J, Gundlach H, et al. The Sorghum bicolor genome and the diversification of grasses. Nature 2009, 457:551–556. pmid:19189423
- View Article
- PubMed/NCBI
- Google Scholar
30. Paterson AH, Wang X, Tang H, Lee TH. Synteny and Genomic Rearrangements. In: Wendel JF, Greilhuber J, Dolezel J, Leitch IJ, editors. Plant Genome Diversity Volume 1: Springer Verlag Wien; pp195–207; 2012.
31. Hanna WW, Burton GW. Genetics of red and purple plant color in pearl millet. J Hered. 1992; 83, 386–388.
- View Article
- Google Scholar
32. Young BA, Sherwood RT, Bashaw EC. Cleared-pistil and thick-sectioning techniques for detecting aposporous apomixis in grasses. Can J Bot. 1979; 57: 1668–1672.
- View Article
- Google Scholar

[ref1] 1. Bennetzen JL. Patterns in grass genome evolution. Curr Opin Plant Biol. 2007; 10: 176–181. pmid:17291821
View Article
PubMed/NCBI
Google Scholar

[2] View Article

[3] PubMed/NCBI

[4] Google Scholar

[ref2] 2. Gaut BS. Evolutionary dynamics of grass genomes. New Phytol. 2002; 154: 15–28.
View Article
Google Scholar

[6] View Article

[7] Google Scholar

[ref3] 3. Wang X, Wang J, Jin D, Guo H, Lee T-H, Liu T, et al. (2015) Genome alignment spanning major Poaceae lineages reveals heterogeneous evolutionary rates and alters inferred dates for key evolutionary events, Molec Plant 2015; 8:885–898. http://dx.doi.org/10.1016/j.molp.2015.04.004.
View Article
Google Scholar

[9] View Article

[10] Google Scholar

[ref4] 4. Bennetzen JL, Schmutz J, Wang H, Percifield R, Hawkins J, Pontaroli AC, et al. Reference genome sequence of the model plant Setaria. Nature Biotech. 2012; 30: 555–561.
View Article
Google Scholar

[12] View Article

[13] Google Scholar

[ref5] 5. Nogler GA. Gametophytic apomixis. In: Johri BM (ed) Embryology of Angiosperms, Springer Verlag, pp 475–518; 1984.

[ref6] 6. Ozias-Akins P, Roche D, Hanna WW. 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. Proc Natl Acad Sci U S A. 1998; 95: 5127–5132. pmid:9560240
View Article
PubMed/NCBI
Google Scholar

[16] View Article

[17] PubMed/NCBI

[18] Google Scholar

[ref7] 7. Roche D, Cong P, Chen Z, Hanna WW, Gustine DL, Sherwood RT, et al. An apospory-specific genomic region is conserved between Buffelgrass (Cenchrus ciliaris L.) and Pennisetum squamulatum Fresen. Plant J. 1999; 19: 203–208. pmid:10476067
View Article
PubMed/NCBI
Google Scholar

[20] View Article

[21] PubMed/NCBI

[22] Google Scholar

[ref8] 8. Akiyama Y, Conner JA, Goel S, Morishige DT, Mullet JE, Hanna WW, et al. High-resolution physical mapping in Pennisetum squamulatum reveals extensive chromosomal heteromorphism of the genomic region associated with apomixis. Plant Physiol. 2004; 134: 1733–1741. pmid:15064383
View Article
PubMed/NCBI
Google Scholar

[24] View Article

[25] PubMed/NCBI

[26] Google Scholar

[ref9] 9. Goel S, Chen Z, Akiyama Y, Conner JA, Basu M, Gualtieri G, et al. Comparative physical mapping of the apospory-specific genomic region in two apomictic grasses: Pennisetum squamulatum and Cenchrus ciliaris. Genetics. 2006; 173: 389–400. pmid:16547108
View Article
PubMed/NCBI
Google Scholar

[28] View Article

[29] PubMed/NCBI

[30] Google Scholar

[ref10] 10. Akiyama Y, Goel S, Conner JA, Hanna WW, Yamada-Akiyama H, Ozias-Akins P. Evolution of the apomixis transmitting chromosome in Pennisetum. BMC Evol Biol. 2011; 11: 289–289. pmid:21975191
View Article
PubMed/NCBI
Google Scholar

[32] View Article

[33] PubMed/NCBI

[34] Google Scholar

[ref11] 11. Conner JA, Goel S, Gunawan G, Cordonnier-Pratt MM, Johnson VE, Liang C, et al. Sequence analysis of bacterial artificial chromosome clones from the apospory-specific genomic region of Pennisetum and Cenchrus. Plant Physiol. 2008; 147: 1396–1411. pmid:18508959
View Article
PubMed/NCBI
Google Scholar

[36] View Article

[37] PubMed/NCBI

[38] Google Scholar

[ref12] 12. Singh M, Conner JA, Zeng YJ, Hanna WW, Johnson VE, Ozias-Akins P. Characterization of apomictic BC₇ and BC₈ pearl millet: meiotic chromosome behavior and construction of an ASGR-carrier chromosome-specific library. Crop Sci. 2010; 50: 892–902.
View Article
Google Scholar

[40] View Article

[41] Google Scholar

[ref13] 13. Zeng Y, Conner J, Ozias-Akins P. Identification of ovule transcripts from the Apospory-Specific Genomic Region (ASGR)-carrier chromosome. BMC Genomics. 2011; 12: 206–206. pmid:21521529
View Article
PubMed/NCBI
Google Scholar

[43] View Article

[44] PubMed/NCBI

[45] Google Scholar

[ref14] 14. Conner JA, Gunawan G, Ozias-Akins P. Recombination within the apospory specific genomic region leads to the uncoupling of apomixis components in Cenchrus ciliaris. Planta. 2013; 238: 51–63. pmid:23553451
View Article
PubMed/NCBI
Google Scholar

[47] View Article

[48] PubMed/NCBI

[49] Google Scholar

[ref15] 15. Zerbino DR. Using the Velvet de novo assembler for short-read sequencing technologies. Current Protocols In Bioinformatics / Editoral Board, Andreas D. Baxevanis … [et al.] Chapter 11: 11.15–11.15; 2010.

[ref16] 16. Roche D, Conner JA, Budiman A, Frisch D, Wing R, Hanna W, Ozias-Akins P. Construction of BAC libraries from two apomictic grasses to study the microcolinearity of their apospory-specific genomic regions. Theor Appl Genet. 2002; 104: 804–812. pmid:12582640
View Article
PubMed/NCBI
Google Scholar

[52] View Article

[53] PubMed/NCBI

[54] Google Scholar

[ref17] 17. Ingham LD, Hanna WW, Baier JW, Hannah LC. Origin of the main class of repetitive DNA within selected Pennisetum species. Mol Gen Genet. 1993; 238: 350–356. pmid:8492802
View Article
PubMed/NCBI
Google Scholar

[56] View Article

[57] PubMed/NCBI

[58] Google Scholar

[ref18] 18. Andres RJ and Kuraparthy V. Development of an improved method of mitotic metaphase chromosome preparation compatible for fluorescence in situ hybridization in cotton. J Cotton Sci. 2013; 17:149–156.
View Article
Google Scholar

[60] View Article

[61] Google Scholar

[ref19] 19. Iwata A, Tek AL, Richard MM, Abernathy B, Fonsêca A, Schmutz J, et al. Identification and characterization of functional centromeres of the common bean. Plant J. 2013; 76: 47–60. pmid:23795942
View Article
PubMed/NCBI
Google Scholar

[63] View Article

[64] PubMed/NCBI

[65] Google Scholar

[ref20] 20. Fukui K. Plant chromosome at mitosis. In: Fukui K and Nagayama S (eds) Plant chromosomes laboratory methods, edited by FUKUI K. and NAGAYAMA S.. CRC Press, Boca Raton, FL, pp 1–17; 1996.

[ref21] 21. Kirov I, Divashuk M, Van Laere K, Soloviev A, Khrustaleva L. An easy "SteamDrop" method for high quality plant chromosome preparation. Mol Cytogenet. 2014; 7: 21 pmid:24602284
View Article
PubMed/NCBI
Google Scholar

[68] View Article

[69] PubMed/NCBI

[70] Google Scholar

[ref22] 22. Zhong X, Fransz PF, van Eden JW, Zabel P, van Kammen Ab, de Jong JH. High-resolution mapping on pachytene chromosomes and extended DNA fibres by fluorescence in-situ hybridisation. Plant Molec Biol Rep. 1996; 14: 232–242.
View Article
Google Scholar

[72] View Article

[73] Google Scholar

[ref23] 23. Devos KM, Pittaway TS, Reynolds A, Gale MD. Comparative mapping reveals complex relationship between the pearl millet genome and those of foxtail millet and rice. Theor Appl Genet. 2000; 100:190–198.
View Article
Google Scholar

[75] View Article

[76] Google Scholar

[ref24] 24. Jessup RW, Burson BL, Burow GB, Wang YWCC, Li Z, Paterson AH, Hussey MA. Disomic inheritance, suppressed recombination, and allelic interactions govern apospory in buffelgrass as revealed by genome mapping. Crop Sci. 2002; 42: 1688–1694.
View Article
Google Scholar

[78] View Article

[79] Google Scholar

[ref25] 25. Pupilli F, Labombarda P, Caceres ME, Quarin CL, Arcioni S. The chromosome segment related to apomixis in Paspalum simplex is homoeologous to the telomeric region of the long arm of rice chromosome 12. Mol Breeding 2001, 8:53–61.
View Article
Google Scholar

[81] View Article

[82] Google Scholar

[ref26] 26. Pupilli F, Martinez EJ, Busti A, Calderini O, Quarin CL, Arcioni S. Comparative mapping reveals partial conservation of synteny at the apomixis locus in Paspalum spp. Molecular Genetics and Genomics 2004, 270:539–548. pmid:14648202
View Article
PubMed/NCBI
Google Scholar

[84] View Article

[85] PubMed/NCBI

[86] Google Scholar

[ref27] 27. Podio M. Rodriuez MP, Felitti S, Stein J, Martinez EJ, Siena LA, et al. Sequence characterization, in silico mapping and cytosine methylation analysis of markers linked to apospory in Paspalum notatum Genetics and Molecular Biology 2012, 35 (4): 827–837. pmid:23271945
View Article
PubMed/NCBI
Google Scholar

[88] View Article

[89] PubMed/NCBI

[90] Google Scholar

[ref28] 28. Hojsgaard DH, Martinez EJ, Acuna CA, Quarin CL, Pupilli F. A molecular map of the apomixis-control locus in Paspalum procurrens and its comparative analysis with other species of Paspalum. Theor Appl Genet 2011, 123:959–971. pmid:21713535
View Article
PubMed/NCBI
Google Scholar

[92] View Article

[93] PubMed/NCBI

[94] Google Scholar

[ref29] 29. Paterson AH, Bowers JE, Bruggmann R, Dubchak I, Grimwood J, Gundlach H, et al. The Sorghum bicolor genome and the diversification of grasses. Nature 2009, 457:551–556. pmid:19189423
View Article
PubMed/NCBI
Google Scholar

[96] View Article

[97] PubMed/NCBI

[98] Google Scholar

[ref30] 30. Paterson AH, Wang X, Tang H, Lee TH. Synteny and Genomic Rearrangements. In: Wendel JF, Greilhuber J, Dolezel J, Leitch IJ, editors. Plant Genome Diversity Volume 1: Springer Verlag Wien; pp195–207; 2012.

[ref31] 31. Hanna WW, Burton GW. Genetics of red and purple plant color in pearl millet. J Hered. 1992; 83, 386–388.
View Article
Google Scholar

[101] View Article

[102] Google Scholar

[ref32] 32. Young BA, Sherwood RT, Bashaw EC. Cleared-pistil and thick-sectioning techniques for detecting aposporous apomixis in grasses. Can J Bot. 1979; 57: 1668–1672.
View Article
Google Scholar

[104] View Article

[105] Google Scholar

Figures

Abstract

Introduction

Materials and Methods

Plant material

Extension of ASGR-carrier chromosome transcript information

Identification of ASGR-carrier chromosome BACs

FISH

FISH probes.

Chromosome preparation.

Fluorescence in situ hybridization and detection.

Deletion Study

Pollen irradiation.

Molecular screen of plant lines derived from irradiated pollen crosses.

Results and Discussion

Supporting Information

S1 Fig. Graphical Overview of Deletion Study Screen.

S1 File. PS_Contigs.txt.

Acknowledgments

Author Contributions

References