A k-mer-based bulked segregant analysis approach to map seed traits in unphased heterozygous potato genomes.

High-throughput sequencing-based methods for bulked segregant analysis (BSA) allow for the rapid identification of genetic markers associated with traits of interest. BSA studies have successfully identified qualitative (binary) and quantitative trait loci (QTLs) using QTL mapping. However, most require population structures that fit the models available and a reference genome. Instead, high-throughput short-read sequencing can be combined with BSA of k-mers (BSA-k-mer) to map traits that appear refractory to standard approaches. This method can be applied to any organism and is particularly useful for species with genomes diverged from the closest sequenced genome. It is also instrumental when dealing with highly heterozygous and potentially polyploid genomes without phased haplotype assemblies and for which a single haplotype can control a trait. Finally, it is flexible in terms of population structure. Here, we apply the BSA-k-mer method for the rapid identification of candidate regions related to seed spot and seed size in diploid potato. Using a mixture of F1 and F2 individuals from a cross between 2 highly heterozygous parents, candidate sequences were identified for each trait using the BSA-k-mer approach. Using parental reads, we were able to determine the parental origin of the loci. Finally, we mapped the identified k-mers to a closely related potato genome to validate the method and determine the genomic loci underlying these sequences. The location identified for the seed spot matches with previously identified loci associated with pigmentation in potato. The loci associated with seed size are novel. Both loci are relevant in future breeding toward true seeds in potato.

LD-CNV: rapid and simple discovery of chromosomal translocations using linkage disequilibrium between copy number variable loci.

Large-scale structural variations, such as chromosomal translocations, can have profound effects on fitness and phenotype, but are difficult to identify and characterize. Here, we describe a simple and effective method aimed at identifying translocations using only the dosage of sequence reads mapped on the reference genome. We binned reads on genomic segments sized according to sequencing coverage and identified instances when copy number segregated in populations. For each dosage-polymorphic 1 Mb bin, we tested independence, effectively an apparent linkage disequilibrium (LD), with other variable bins. In nine potato (Solanum tuberosum) dihaploid families translocations affecting pericentromeric regions were common and in two cases were due to genomic misassembly. In two populations, we found evidence for translocation affecting euchromatic arms. In cv. PI 310467, a nonreciprocal translocation between chromosomes (chr.) 7 and 8 resulted in a 5-3 copy number change affecting several Mb at the respective chromosome tips. In cv. "Alca Tarma," the terminal arm of chr. 4 translocated to the tip of chr. 1. Using oligonucleotide-based fluorescent in situ hybridization painting probes (oligo-FISH), we tested and confirmed the predicted arrangement in PI 310467. In 192 natural accessions of Arabidopsis thaliana, dosage haplotypes tended to vary continuously and resulted in higher noise, while apparent LD between pericentromeric regions suggested the effect of repeats. This method, LD-CNV, should be useful in species where translocations are suspected because it tests linkage without the need for genotyping.

Rare instances of haploid inducer DNA in potato dihaploids and ploidy-dependent genome instability.

In cultivated tetraploid potato (Solanum tuberosum), reduction to diploidy (dihaploidy) allows for hybridization to diploids and introgression breeding and may facilitate the production of inbreds. Pollination with haploid inducers (HIs) yields maternal dihaploids, as well as triploid and tetraploid hybrids. Dihaploids may result from parthenogenesis, entailing the development of embryos from unfertilized eggs, or genome elimination, entailing missegregation and the loss of paternal chromosomes. A sign of genome elimination is the occasional persistence of HI DNA in some dihaploids. We characterized the genomes of 919 putative dihaploids and 134 hybrids produced by pollinating tetraploid clones with three HIs: IVP35, IVP101, and PL-4. Whole-chromosome or segmental aneuploidy was observed in 76 dihaploids, with karyotypes ranging from 2n = 2x - 1 = 23 to 2n = 2x + 3 = 27. Of the additional chromosomes in 74 aneuploids, 66 were from the non-inducer parent and 8 from the inducer parent. Overall, we detected full or partial chromosomes from the HI parent in 0.87% of the dihaploids, irrespective of parental genotypes. Chromosomal breaks commonly affected the paternal genome in the dihaploid and tetraploid progeny, but not in the triploid progeny, correlating instability to sperm ploidy and to haploid induction. The residual HI DNA discovered in the progeny is consistent with genome elimination as the mechanism of haploid induction.

Genomic Outcomes of Haploid Induction Crosses in Potato (Solanum tuberosum L.).

The challenges of breeding autotetraploid potato (Solanum tuberosum) have motivated the development of alternative breeding strategies. A common approach is to obtain uniparental dihaploids from a tetraploid of interest through pollination with S. tuberosum Andigenum Group (formerly S. phureja) cultivars. The mechanism underlying haploid formation of these crosses is unclear, and questions regarding the frequency of paternal DNA transmission remain. Previous reports have described aneuploid and euploid progeny that, in some cases, displayed genetic markers from the haploid inducer (HI). Here, we surveyed a population of 167 presumed dihaploids for large-scale structural variation that would underlie chromosomal addition from the HI, and for small-scale introgression of genetic markers. In 19 progeny, we detected 10 of the 12 possible trisomies and, in all cases, demonstrated the noninducer parent origin of the additional chromosome. Deep sequencing indicated that occasional, short-tract signals appearing to be of HI origin were better explained as technical artifacts. Leveraging recurring copy number variation patterns, we documented subchromosomal dosage variation indicating segregation of polymorphic maternal haplotypes. Collectively, 52% of the assayed chromosomal loci were classified as dosage variable. Our findings help elucidate the genomic consequences of potato haploid induction and suggest that most potato dihaploids will be free of residual pollinator DNA.

Regeneration of Solanum tuberosum Plants from Protoplasts Induces Widespread Genome Instability.

Nontransgenic genome editing in regenerable protoplasts, plant cells free of their cell wall, could revolutionize crop improvement because it reduces regulatory and technical complexity. However, plant tissue culture is known to engender frequent unwanted variation, termed somaclonal variation. To evaluate the contribution of large-scale genome instability to this phenomenon, we analyzed potatoes (Solanum tuberosum) regenerated from either protoplasts or stem explants for copy number changes by comparison of Illumina read depth. Whereas a control set of eight plants that had been propagated by cuttings displayed no changes, all 15 protoplast regenerants tested were affected by aneuploidy or structural chromosomal changes. Certain chromosomes displayed segmental deletions and duplications ranging from one to many. Resampling different leaves of the same plant found differences in three regenerants, indicating frequent persistence of instability. By comparison, 33 regenerants from stem explants used for Agrobacterium-mediated transformation displayed less frequent but still considerable (18%) large-scale copy number changes. Repetition of certain instability patterns suggested greater susceptibility in specific genomic sites. These results indicate that tissue culture, depending on the protocol used, can induce genomic instability resulting in large-scale changes likely to compromise final plant phenotype.

Point Mutations in Centromeric Histone Induce Post-zygotic Incompatibility and Uniparental Inheritance.

The centromeric histone 3 variant (CENH3, aka CENP-A) is essential for the segregation of sister chromatids during mitosis and meiosis. To better define CENH3 functional constraints, we complemented a null allele in Arabidopsis with a variety of mutant alleles, each inducing a single amino acid change in conserved residues of the histone fold domain. Many of these transgenic missense lines displayed wild-type growth and fertility on self-pollination, but exhibited frequent post-zygotic death and uniparental inheritance when crossed with wild-type plants. The failure of centromeres marked by these missense mutation in the histone fold domain of CENH3 reproduces the genome elimination syndromes described with chimeric CENH3 and CENH3 from diverged species. Additionally, evidence that a single point mutation is sufficient to generate a haploid inducer provide a simple one-step method for the identification of non-transgenic haploid inducers in existing mutagenized collections of crop species. As proof of the extreme simplicity of this approach to create haploid-inducing lines, we performed an in silico search for previously identified point mutations in CENH3 and identified an Arabidopsis line carrying the A86V substitution within the histone fold domain. This A87V non-transgenic line, while fully fertile on self-pollination, produced postzygotic death and uniparental haploids when crossed to wild type.

A System for Dosage-Based Functional Genomics in Poplar.

Altering gene dosage through variation in gene copy number is a powerful approach to addressing questions regarding gene regulation, quantitative trait loci, and heterosis, but one that is not easily applied to sexually transmitted species. Elite poplar (Populus spp) varieties are created through interspecific hybridization, followed by clonal propagation. Altered gene dosage relationships are believed to contribute to hybrid performance. Clonal propagation allows for replication and maintenance of meiotically unstable ploidy or structural variants and provides an alternative approach to investigating gene dosage effects not possible in sexually propagated species. Here, we built a genome-wide structural variation system for dosage-based functional genomics and breeding of poplar. We pollinated Populus deltoides with gamma-irradiated Populus nigra pollen to produce >500 F1 seedlings containing dosage lesions in the form of deletions and insertions of chromosomal segments (indel mutations). Using high-precision dosage analysis, we detected indel mutations in ∼55% of the progeny. These indels varied in length, position, and number per individual, cumulatively tiling >99% of the genome, with an average of 10 indels per gene. Combined with future phenotype and transcriptome data, this population will provide an excellent resource for creating and characterizing dosage-based variation in poplar, including the contribution of dosage to quantitative traits and heterosis.

Two callose synthases, GSL1 and GSL5, play an essential and redundant role in plant and pollen development and in fertility.

Callose, a beta-1,3-glucan that is widespread in plants, is synthesized by callose synthase. Arabidopsis thaliana contains a family of 12 putative callose synthase genes (GSL1-12). The role of callose and of the individual genes in plant development is still largely uncertain. We have now used TILLING and T-DNA insertion mutants (gsl1-1, gsl5-2 and gsl5-3) to study the role of two closely related and linked genes, GSL1 and GSL5, in sporophytic development and in reproduction. Both genes are expressed in all parts of the plant. Sporophytic development was nearly normal in gsl1-1 homozygotes and only moderately defective in homozygotes for either of the two gsl5 alleles. On the other hand, plants that were gsl1-1/+ gsl5/gsl5 were severely defective, with smaller leaves, shorter roots and bolts and smaller flowers. Plants were fertile when the sporophytes had either two wild-type GSL1 alleles, or one GSL5 allele in a gsl1-1 background, but gsl1-1/+ gsl5/gsl5 plants produced an extremely reduced number of viable seeds. A chromosome with mutations in both GSL1 and GSL5 rendered pollen infertile, although such a chromosome could be transmitted via the egg. As a result, it was not possible to obtain plants that were homozygous for mutations in both the GSL genes. Pollen grain development was severely affected in double mutant plants. Many pollen grains were collapsed and inviable in the gsl1-1/gsl1-1 gsl5/+ and gsl1-1/+ gsl5/gsl5 plants. In addition, gsl1-1/+ gsl5/gsl5 plants produced abnormally large pollen with unusual pore structures, and had problems with tetrad dissociation. In this particular genotype, while the callose wall formed around the pollen mother cells, no callose wall separated the resulting tetrads. We conclude that GSL1 and GSL5 play important, but at least partially redundant roles in both sporophytic development and in the development of pollen. They are responsible for the formation of the callose wall that separates the microspores of the tetrad, and also play a gametophytic role later in pollen grain maturation. Other GSL genes may control callose formation at different steps during pollen development.

Aneuploidy and genetic variation in the Arabidopsis thaliana triploid response.

Polyploidy, the inheritance of more than two genome copies per cell, has played a major role in the evolution of higher plants. Little is known about the transition from diploidy to polyploidy but in some species, triploids are thought to function as intermediates in this transition. In contrast, in other species triploidy is viewed as a block. We investigated the responses of Arabidopsis thaliana to triploidy. The role of genetic variability was tested by comparing triploids generated from crosses between Col-0, a diploid, and either a natural autotetraploid (Wa-1) or an induced tetraploid of Col-0. In this study, we demonstrate that triploids of A. thaliana are fertile, producing a swarm of different aneuploids. Propagation of the progeny of a triploid for a few generations resulted in diploid and tetraploid cohorts. This demonstrated that, in A. thaliana, triploids can readily form tetraploids and function as bridges between euploid types. Genetic analysis of recombinant inbred lines produced from a triploid identified a locus on chromosome I exhibiting allelic bias in the tetraploid lines but not in the diploid lines. Thus, genetic variation was subject to selection contingent on the final ploidy and possibly acting during the protracted aneuploid phase.

Mismatch cleavage by single-strand specific nucleases.

We have investigated the ability of single-strand specific (sss) nucleases from different sources to cleave single base pair mismatches in heteroduplex DNA templates used for mutation and single-nucleotide polymorphism analysis. The TILLING (Targeting Induced Local Lesions IN Genomes) mismatch cleavage protocol was used with the LI-COR gel detection system to assay cleavage of amplified heteroduplexes derived from a variety of induced mutations and naturally occurring polymorphisms. We found that purified nucleases derived from celery (CEL I), mung bean sprouts and Aspergillus (S1) were able to specifically cleave nearly all single base pair mismatches tested. Optimal nicking of heteroduplexes for mismatch detection was achieved using higher pH, temperature and divalent cation conditions than are routinely used for digestion of single-stranded DNA. Surprisingly, crude plant extracts performed as well as the highly purified preparations for this application. These observations suggest that diverse members of the S1 family of sss nucleases act similarly in cleaving non-specifically at bulges in heteroduplexes, and single-base mismatches are the least accessible because they present the smallest single-stranded region for enzyme binding. We conclude that a variety of sss nucleases and extracts can be effectively used for high-throughput mutation and polymorphism discovery.

Centromeric localization and adaptive evolution of an Arabidopsis histone H3 variant.

Centromeric H3-like histones, which replace histone H3 in the centromeric chromatin of animals and fungi, have not been reported in plants. We identified a histone H3 variant from Arabidopsis thaliana that encodes a centromere-identifying protein designated HTR12. By immunological detection, HTR12 localized at centromeres in both mitotic and meiotic cells. HTR12 signal revealed tissue- and stage-specific differences in centromere morphology, including a distended bead-like structure in interphase root tip cells. The anti-HTR12 antibody also detected spherical organelles in meiotic cells. Although the antibody does not label centromeres in the closely related species Arabidopsis arenosa, HTR12 signal was found on all centromeres in allopolyploids of these two species. Comparison of the HTR12 genes of A. thaliana and A. arenosa revealed striking adaptive evolution in the N-terminal tail of the protein, similar to the pattern seen in its counterpart in Drosophila. This finding suggests that the same evolutionary forces shape centromeric chromatin in both animals and plants.