Genotypic Frequencies at Equilibrium for Polysomic Inheritance Under Double-Reduction

A locus-dependent mixed inheritance in the segmental allohexaploid sweetpotato (Ipomoea batatas [L.] Lam)

Two interrelated aspects of the sweetpotato genome, its polyploid origin and inheritance type, remain uncertain. We recently proposed a segmental allohexaploid sweetpotato and thus sought to clarify its inheritance type by direct analyses of homoeolog segregations at selected single-copy loci. For such analyses, we developed a digital quantitative PCR genotyping method using one nondiscriminatory and three discriminatory probes for each selected locus to discriminate and quantify three homoeolog-differentiating variation types (homoeolog-types) in genomic DNA samples for genotype fitting and constructed a F2 population for segregation analyses. We confirmed inter-subgenomic distinctions of three identified homoeolog-types at each of five selected loci by their interspecific differentiations among 14 species in Ipomoea section batatas and genotyped the loci in 549 F2 lines, selected F1 progenies, and their founding parents. Segregation and genotype analyses revealed a locus-dependent mixed inheritance (disomic, polysomic, and intermediate types) of the homoeolog-types at 4 loci in the F2 population, displaying estimated disomic-inheritance frequencies of 0, 2.72%, 14.52%, and 36.92%, and probably in the F1 population too. There were also low-frequency non-hexaploid F1 and F2 genotypes that were probably derived from double-reduction recombination or partially unreduced gametes, and F2 genotypes of apparent aneuploids/dysploids with neopolyploid-like frequencies. Additional analyses of homoeolog-type genotypes at the 5 loci in 46 lines from various regions revealed locus-dependent selection biases, favoring genotypes having more of one homoeolog-type, i.e. more of di- or homogenized homoeolog-type composition, and one-direction ploidy trending among apparent aneuploids/dysploids. These inheritance features pointed to an evolving segmental allohexaploid sweetpotato impacted by selection biases.

GenAPoPop 1.0: A user-friendly software to analyse genetic diversity and structure from partially clonal and selfed autopolyploid organisms

Autopolyploidy is quite common in most clades of eukaryotes. The emergence of sequence-based genotyping methods with individual and marker tags now enables confident allele dosage, overcoming the main obstacle to the democratization of the population genetic approaches when studying ecology and evolution of autopolyploid populations and species. Reproductive modes, including clonality, selfing and allogamy, have deep consequences on the ecology and evolution of population and species. Analysing genetic diversity and its dynamics over generations is one efficient way to infer the relative importance of clonality, selfing and allogamy in populations. GenAPoPop is a user-friendly solution to compute the specific corpus of population genetic indices, including indices about genotypic diversity, needed to analyse partially clonal, selfed and allogamous polysomic populations genotyped with confident allele dosage. It also easily provides the posterior probabilities of quantitative reproductive modes in autopolyploid populations genotyped at two-time steps and a graphical representation of the minimum spanning trees of the genetic distances between polyploid individuals, facilitating the interpretation of the genetic coancestry between individuals in hierarchically structured populations. GenAPoPop complements the previously existing solutions, including SPAGEDI and POLYGENE, to use genotypings to study the ecology and evolution of autopolyploid populations. It was specially developed with a simple graphical interface and workflow, and comes with a simulator to facilitate practical courses and teaching of population genetics for autopolyploid populations.

Bayesian tests for random mating in polyploids

Hardy-Weinberg proportions (HWP) are often explored to evaluate the assumption of random mating. However, in autopolyploids, organisms with more than two sets of homologous chromosomes, HWP and random mating are different hypotheses that require different statistical testing approaches. Currently, the only available methods to test for random mating in autopolyploids (i) heavily rely on asymptotic approximations and (ii) assume genotypes are known, ignoring genotype uncertainty. Furthermore, these approaches are all frequentist, and so do not carry the benefits of Bayesian analysis, including ease of interpretability, incorporation of prior information, and consistency under the null. Here, we present Bayesian approaches to test for random mating, bringing the benefits of Bayesian analysis to this problem. Our Bayesian methods also (i) do not rely on asymptotic approximations, being appropriate for small sample sizes, and (ii) optionally account for genotype uncertainty via genotype likelihoods. We validate our methods in simulations and demonstrate on two real datasets how testing for random mating is more useful for detecting genotyping errors than testing for HWP (in a natural population) and testing for Mendelian segregation (in an experimental S1 population). Our methods are implemented in Version 2.0.2 of the hwep R package on the Comprehensive R Archive Network https://cran.r-project.org/package=hwep.

Double reduction estimation and equilibrium tests in natural autopolyploid populations

Many bioinformatics pipelines include tests for equilibrium. Tests for diploids are well studied and widely available, but extending these approaches to autopolyploids is hampered by the presence of double reduction, the comigration of sister chromatid segments into the same gamete during meiosis. Though a hindrance for equilibrium tests, double reduction rates are quantities of interest in their own right, as they provide insights about the meiotic behavior of autopolyploid organisms. Here, we develop procedures to (i) test for equilibrium while accounting for double reduction, and (ii) estimate the double reduction rate given equilibrium. To do so, we take two approaches: a likelihood approach, and a novel U-statistic minimization approach that we show generalizes the classical equilibrium χ2 test in diploids. For small sample sizes and uncertain genotypes, we further develop a bootstrap procedure based on our U-statistic to test for equilibrium. We validate our methods on both simulated and real data.

vcfpop: Performing population genetics analyses for autopolyploids and aneuploids based on next-generation sequencing data sets

Polyploids are cells or organisms with a genome consisting of more than two sets of homologous chromosomes. Polyploid plants have important traits that facilitate speciation and are thus often model systems for evolutionary, molecular ecology and agricultural studies. However, due to their unusual mode of inheritance and double-reduction, diploid models of population genetic analysis cannot properly be applied to autopolyploids. To overcome this problem, we developed a software package entitled vcfpop to perform a variety of population genetic analyses for autopolyploids, such as parentage analysis, analysis of molecular variance, principal coordinates analysis, hierarchical clustering analysis and Bayesian clustering. We used three data sets to evaluate the capability of vcfpop to analyse large data sets on a desktop computer. This software is freely available at http://github.com/huangkang1987/vcfpop.

The recombination landscape and multiple QTL mapping in a Solanum tuberosum cv. 'Atlantic'-derived F1 population

There are many challenges involved with the genetic analyses of autopolyploid species, such as the tetraploid potato, Solanum tuberosum (2n = 4x = 48). The development of new analytical methods has made it valuable to re-analyze an F₁ population (n = 156) derived from a cross involving 'Atlantic', a widely grown chipping variety in the USA. A fully integrated genetic map with 4285 single nucleotide polymorphisms, spanning 1630 cM, was constructed with MAPpoly software. We observed that bivalent configurations were the most abundant ones (51.0~72.4% depending on parent and linkage group), though multivalent configurations were also observed (2.2~39.2%). Seven traits were evaluated over four years (2006-8 and 2014) and quantitative trait loci (QTL) mapping was carried out using QTLpoly software. Based on a multiple-QTL model approach, we detected 21 QTL for 15 out of 27 trait-year combination phenotypes. A hotspot on linkage group 5 was identified with co-located QTL for maturity, plant yield, specific gravity, and internal heat necrosis resistance evaluated over different years. Additional QTL for specific gravity and dry matter were detected with maturity-corrected phenotypes. Among the genes around QTL peaks, we found those on chromosome 5 that have been previously implicated in maturity (StCDF1) and tuber formation (POTH1). These analyses have the potential to provide insights into the biology and breeding of tetraploid potato and other autopolyploid species.

Performing parentage analysis for polysomic inheritances based on allelic phenotypes

Polyploidy poses several problems for parentage analysis. We present a new polysomic inheritance model for parentage analysis based on genotypes or allelic phenotypes to solve these problems. The effects of five factors are simultaneously accommodated in this model: (1) double-reduction, (2) null alleles, (3) negative amplification, (4) genotyping errors and (5) self-fertilization. To solve genotyping ambiguity (unknown allele dosage), we developed a new method to establish the likelihood formulas for allelic phenotype data and to simultaneously include the effects of our five chosen factors. We then evaluated and compared the performance of our new method with three established methods by using both simulated data and empirical data from the cultivated blueberry (Vaccinium corymbosum). We also developed and compared the performance of two additional estimators to estimate the genotyping error rate and the sample rate. We make our new methods freely available in the software package polygene, at http://github.com/huangkang1987/polygene.

Chromosome inheritance and meiotic stability in allopolyploid Brassica napus

Homoeologous recombination, aneuploidy, and other genetic changes are common in resynthesized allopolyploid Brassica napus. In contrast, the chromosomes of cultivars have long been considered to be meiotically stable. To gain a better understanding of the underlying mechanisms leading to stabilization in the allopolyploid, the behavior of chromosomes during meiosis can be compared by unambiguous chromosome identification between resynthesized and natural B. napus. Compared with natural B. napus, resynthesized lines show high rates of nonhomologous centromere association, homoeologous recombination leading to translocation, homoeologous chromosome replacement, and association and breakage of 45S rDNA loci. In both natural and resynthesized B. napus, we observed low rates of univalents, A–C bivalents, and early sister chromatid separations. Reciprocal homoeologous chromosome exchanges and double reductions were photographed for the first time in meiotic telophase I. Meiotic errors were non-uniformly distributed across the genome in resynthesized B. napus, and in particular homoeologs sharing synteny along their entire length exhibited multivalents at diakinesis and polysomic inheritance at telophase I. Natural B. napus appeared to resolve meiotic errors mainly by suppressing homoeologous pairing, resolving nonhomologous centromere associations and 45S rDNA associations before diakinesis, and reducing homoeologous cross-overs.

Isolation and characterization of twelve polymorphic microsatellite markers in the endangered Hopea hainanensis (Dipterocarpaceae)

Microsatellite markers were isolated and characterized for Hopea hainanensis Merrill & Chun, an endangered tree species with scattered distribution in Hainan Island and northern Vietnam. Twenty-six microsatellite markers were developed based on next-generation sequencing data and were genotyped by capillary electrophoresis on an ABI 3730xl DNA Analyzer. Twelve markers were found to be polymorphic in H. hainanensis. GENODIVE analyses indicated that the number of alleles ranged from 2 to 6 per locus, and the observed and expected heterozygosity varied from 0 to 0.755 and from 0.259 to 0.779, respectively. Primer transferability was tested with Hopea chinensis Hand.-Mazz. and Hopea reticulata Tardieu, in which 3 and 7 microsatellite markers were found to be polymorphic, separately. The results showed that H. reticulata and H. hainanensis had similar levels of genetic diversity. A neighbor joining dendrogram clustered all individuals into two major groups, one of which was exclusively constituted by H. hainanensis, while the other consisted of two subgroups, corresponding to H. reticulata and H. chinensis, respectively. The 12 polymorphic microsatellite markers could be applied to study genetic diversity, population differentiation, mating system, and fine-scale spatial genetic structures of H. hainanensis as well as its close relatives, facilitating the conservation and restoration of these endangered but valuable Hopea species.

Genetic diversity and population structure in the endangered tree Hopea hainanensis (Dipterocarpaceae) on Hainan Island, China

Hopea hainanensis Merrill & Chun (Dipterocarpaceae) is an endangered tree species restricted to Hainan Island, China and a small part of Northern Vietnam. On Hainan Island, it is an important indicator species for tropical forests. However, because of its highly valued timber, H. hainanensis has suffered from overexploitation, leading to a sharp population decline. To facilitate the conservation of this species, genetic diversity and population structure were assessed using 12 SSR markers for 10 populations sampled across Hainan Island. Compared to non-threatened Hopea species, H. hainanensis exhibited reduced overall genetic diversity and increased population differentiation (AMOVA: F_ST = 0.23). Bayesian model-based clustering and principal coordinate analysis consistently assigned H. hainanensis individuals into three genetic groups, which were found to be widespread and overlapping geographically. A Mantel test found no correlation between genetic and geographical distances (r = 0.040, p = 0.418). The observed genetic structure suggests that long-distance gene flow occurred among H. hainanensis populations prior to habitat fragmentation. A recent population bottleneck was revealed, which may cause rapid loss of genetic diversity and increased differentiation across populations. Based on these findings, appropriate strategies for the long-term conservation of the endangered species H. hainanensis are proposed.

A generalized framework for AMOVA with multiple hierarchies and ploidies

The analysis of molecular variance (AMOVA) is a widely used statistical method in population genetics and molecular ecology. The classic framework of AMOVA only supports haploid and diploid data, in which the number of hierarchies ranges from two to four. In practice, natural populations can be classified into more hierarchies, and polyploidy is frequently observed in extant species. The ploidy level may even vary within the same species, and/or within the same individual. We generalized the framework of AMOVA such that it can be used for any number of hierarchies and any level of ploidy. Based on this framework, we present four methods to account for data that are multilocus genotypic and allelic phenotypic (with unknown allele dosage). We use simulated datasets and an empirical dataset to evaluate the performance of our framework. We make freely available our methods in a new software package, polygene, which is freely available at https://github.com/huangkang1987/polygene.

High-Resolution Linkage Map With Allele Dosage Allows the Identification of Regions Governing Complex Traits and Apospory in Guinea Grass (Megathyrsus maximus)

Forage grasses are mainly used in animal feed to fatten cattle and dairy herds, and guinea grass (Megathyrsus maximus) is considered one of the most productive of the tropical forage crops that reproduce by seeds. Due to the recent process of domestication, this species has several genomic complexities, such as autotetraploidy and aposporous apomixis. Consequently, approaches that relate phenotypic and genotypic data are incipient. In this context, we built a linkage map with allele dosage and generated novel information of the genetic architecture of traits that are important for the breeding of M. maximus. From a full-sib progeny, a linkage map containing 858 single nucleotide polymorphism (SNP) markers with allele dosage information expected for an autotetraploid was obtained. The high genetic variability of the progeny allowed us to map 10 quantitative trait loci (QTLs) related to agronomic traits, such as regrowth capacity and total dry matter, and 36 QTLs related to nutritional quality, which were distributed among all homology groups (HGs). Various overlapping regions associated with the quantitative traits suggested QTL hotspots. In addition, we were able to map one locus that controls apospory (apo-locus) in HG II. A total of 55 different gene families involved in cellular metabolism and plant growth were identified from markers adjacent to the QTLs and APOSPORY locus using the Panicum virgatum genome as a reference in comparisons with the genomes of Arabidopsis thaliana and Oryza sativa. Our results provide a better understanding of the genetic basis of reproduction by apomixis and traits important for breeding programs that considerably influence animal productivity as well as the quality of meat and milk.

Development of simple sequence repeat markers for Chamerion angustifolium (Onagraceae)

PREMISE

Rosebay willowherb, or fireweed (Chamerion angustifolium: Onagraceae), has diploid, tetraploid, and hexaploid cytotypes. There are known physiological and ecological differences among the three cytotypes, but genetic differences remain undetermined. We developed simple sequence repeat (SSR) markers for this species.

METHODS AND RESULTS

Leaf samples were collected from three hexaploid C. angustifolium populations. We successfully amplified 16 SSR loci, which were found to be highly polymorphic. The number of alleles, the observed heterozygosity levels, and the expected heterozygosity levels ranged from four to 13, 0.286-0.899, and 0.372-0.871, respectively. Most primers could also be amplified successfully in C. conspersum and the closely related species Epilobium palustre.

CONCLUSIONS

The 16 polymorphic SSR markers developed here will be useful for genetic studies in C. angustifolium and related species.