Single-pollen genotyping to estimate mode of unreduced pollen formation in Citrus tamurana cv. Nishiuchi Konatsu

Consequences of whole genome duplication for 2n pollen performance

The vegetative cell of the angiosperm male gametophyte (pollen) functions as a free-living, single-celled organism that both produces and transports sperm to egg. Whole-genome duplication (WGD) should have strong effects on pollen because of the haploid to diploid transition and because of both genetic and epigenetic effects on cell-level phenotypes. To disentangle historical effects of WGD on pollen performance, studies can compare 1n pollen from diploids to neo-2n pollen from diploids and synthetic autotetraploids to older 2n pollen from established neo-autotetraploids. WGD doubles both gene number and bulk nuclear DNA mass, and a substantial proportion of diploid and autotetraploid heterozygosity can be transmitted to 2n pollen. Relative to 1n pollen, 2n pollen can exhibit heterosis due to higher gene dosage, higher heterozygosity and new allelic interactions. Doubled genome size also has consequences for gene regulation and expression as well as epigenetic effects on cell architecture. Pollen volume doubling is a universal effect of WGD, whereas an increase in aperture number is common among taxa with simultaneous microsporogenesis and pored apertures, mostly eudicots. WGD instantly affects numerous evolved compromises among mature pollen functional traits and these are rapidly shaped by highly diverse tissue interactions and pollen competitive environments in the early post-WGD generations. 2n pollen phenotypes generally incur higher performance costs, and the degree to which these are met or evolve by scaling up provisioning and metabolic vigor needs further study.

High-throughput measuring of meiotic recombination rates in barley pollen nuclei using Crystal Digital PCRTM

Breeding exploits novel allelic combinations assured by meiotic recombination. Barley (Hordeum vulgare) single pollen nucleus genotyping enables measurement of meiotic recombination rates in gametes before fertilization without the need for segregating populations. However, so far, established methods rely on whole-genome amplification of every single pollen nucleus due to their limited DNA content, thus restricting the number of analyzed samples. In this study, high-throughput measurements of meiotic recombination rates in barley pollen nuclei without whole-genome amplification were performed through a Crystal Digital PCR^TM -based genotyping assay. Meiotic recombination rates within two centromeric and two distal chromosomal intervals were measured in hybrid plants by genotyping a total of >42 000 individual pollen nuclei (up to 4900 nuclei analyzed per plant). Determined recombination frequencies in pollen nuclei were similar to frequencies in segregating populations. We improved the efficiency of the genotyping by pretreating the pollen nuclei with a thermostable restriction enzyme. Additional opportunities for a higher sample throughput and a further increase of the genotyping efficiency are presented and discussed. Taken together, single barley pollen nucleus genotyping based on Crystal Digital PCR^TM enables reliable, rapid and high-throughput meiotic recombination measurements within defined chromosomal intervals of intraspecific hybrid plants. The successful encapsulation of nuclei from a range of species with different nuclear and genome sizes suggests that the proposed method is broadly applicable to genotyping single nuclei.

Localization and characterization of Citrus centromeres by combining half-tetrad analysis and CenH3-associated sequence profiling

The physical locations of citrus centromere are revealed by combining genetic and immunological assays for the first time and nine citrus centromere-specific markers for cytogenetics are mined. Centromere localization is challenging, because highly redundant repetitive sequences in centromeric regions make sequence assembly difficult. Although several citrus genomes have been released, the centromeric regions and their characteristics remain to be elucidated. Here, we mapped citrus centromeres through half-tetrad analysis (HTA) that included the genotyping of 54 tetraploid hybrids derived from 2n megagametophytes of Nadorcott tangor with 212 single nucleotide polymorphism (SNP) markers. The sizes of centromeric regions, which estimated based on the heterozygosity restitution rate pattern along the chromosomes, ranged from 1.12 to 18.19 Mb. We also profiled the binding sequences with the centromere-specific histone variant CenH3 by chromatin immunoprecipitation sequencing (ChIP-seq). Based on the positions of the top ten CenH3-enriched contigs, the sizes of centromeric regions were estimated to range from 0.01 to 7.60 Mb and were either adjacent to or included in the centromeric regions identified by HTA. We used DNA probes from two repeats selected from the centromeric regions and seven CenH3-binding centromeric repeats to verify centromeric locations by fluorescence in situ hybridization (FISH). Centromere localization in citrus will contribute to the mining of centromeric/pericentromeric markers, thus to facilitate the rapid identification of mechanisms underlying 2n gamete formation and serve the polyploidy breeding.

Polyploidy in Fruit Tree Crops of the Genus Annona (Annonaceae)

Genome duplication or polyploidy is one of the main factors of speciation in plants. It is especially frequent in hybrids and very valuable in many crops. The genus Annona belongs to the Annonaceae, a family that includes several fruit tree crops, such as cherimoya (Annona cherimola), sugar apple (Annona squamosa), their hybrid atemoya (A. cherimola × A. squamosa) or pawpaw (Asimina triloba). In this work, genome content was evaluated in several Annona species, A. triloba and atemoya. Surprisingly, while the hybrid atemoya has been reported as diploid, flow cytometry analysis of a progeny obtained from an interspecific cross between A. cherimola and A. squamosa showed an unusual ploidy variability that was also confirmed karyotype analysis. While the progeny from intraspecific crosses of A. cherimola showed polyploid genotypes that ranged from 2.5 to 33%, the hybrid atemoyas from the interspecific cross showed 35% of triploids from a total of 186 genotypes analyzed. With the aim of understanding the possible implications of the production of non-reduced gametes, pollen performance, pollen size and frequency distribution of pollen grains was quantified in the progeny of this cross and the parents. A large polymorphism in pollen grain size was found within the interspecific progeny with higher production of unreduced pollen in triploids (38%) than in diploids (29%). Moreover, using PCR amplification of selected microsatellite loci, while 13.7% of the pollen grains from the diploids showed two alleles, 41.28% of the grains from the triploids amplified two alleles and 5.63% showed up to three alleles. This suggests that the larger pollen grains could correspond to diploid and, in a lower frequency, to triploid pollen. Pollen performance was also affected with lower pollen germination in the hybrid triploids than in both diploid parents. The results confirm a higher percentage of polyploids in the interspecific cross, affecting pollen grain size and pollen performance. The occurrence of unreduced gametes in A. cherimola, A. squamosa and their interspecific progeny that may result in abnormalities of ploidy such as the triploids and tetraploids observed in this study, opens an interesting opportunity to study polyploidy in Annonaceae.

Assessing Ploidy Level Analysis and Single Pollen Genotyping of Diploid and Euploid Citrus Genotypes by Fluorescence-Activated Cell Sorting and Whole-Genome Amplification

Flow cytometry is widely used to determine genome size and ploidy level in plants. This technique, when coupled with fluorescence-activated cell sorting (FACS), whole genome amplification and genotyping (WGA), opens up new opportunities for genetic studies of individualized nuclei. This strategy was used to analyze the genetic composition of single pollen nuclei of different citrus species. The flow cytometry and microscope observations allowed us to differentiate the populations of pollen nuclei present in the diploid and euploid genotypes analyzed, showing that citrus has binuclear pollen. We have identified in the "CSO" tangor an additional nuclei population composed by the vegetative plus generative nuclei. Genotyping of this nuclei population revealed that vegetative and generative nuclei show the same genetic configuration. In addition, we have demonstrated the presence of unreduced gametes in the diploid genotype "Mexican lime." Genomic amplification is a robust method for haploid nuclei genotyping with several molecular markers, whereas in diploid nuclei using heterozygous markers showed a bias towards one of the two alleles, limiting the use of this tool in this type of nuclei. We further discuss the importance and applications of single pollen genotyping in citrus genetic studies.

High-quality genetic mapping with ddRADseq in the non-model tree Quercus rubra

Background

Restriction site associated DNA sequencing (RADseq) has the potential to be a broadly applicable, low-cost approach for high-quality genetic linkage mapping in forest trees lacking a reference genome. The statistical inference of linear order must be as accurate as possible for the correct ordering of sequence scaffolds and contigs to chromosomal locations. Accurate maps also facilitate the discovery of chromosome segments containing allelic variants conferring resistance to the biotic and abiotic stresses that threaten forest trees worldwide. We used ddRADseq for genetic mapping in the tree Quercus rubra, with an approach optimized to produce a high-quality map. Our study design also enabled us to model the results we would have obtained with less depth of coverage.

Results

Our sequencing design produced a high sequencing depth in the parents (248×) and a moderate sequencing depth (15×) in the progeny. The digital normalization method of generating a de novo reference and the SAMtools SNP variant caller yielded the most SNP calls (78,725). The major drivers of map inflation were multiple SNPs located within the same sequence (77% of SNPs called). The highest quality map was generated with a low level of missing data (5%) and a genome-wide threshold of 0.025 for deviation from Mendelian expectation. The final map included 849 SNP markers (1.8% of the 78,725 SNPs called). Downsampling the individual FASTQ files to model lower depth of coverage revealed that sequencing the progeny using 96 samples per lane would have yielded too few SNP markers to generate a map, even if we had sequenced the parents at depth 248×.

Conclusions

The ddRADseq technology produced enough high-quality SNP markers to make a moderately dense, high-quality map. The success of this project was due to high depth of coverage of the parents, moderate depth of coverage of the progeny, a good framework map, an optimized bioinformatics pipeline, and rigorous premapping filters. The ddRADseq approach is useful for the construction of high-quality genetic maps in organisms lacking a reference genome if the parents and progeny are sequenced at sufficient depth. Technical improvements in reduced representation sequencing (RRS) approaches are needed to reduce the amount of missing data.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3765-8) contains supplementary material, which is available to authorized users.