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.

The Arabidopsis Hop1 homolog ASY1 mediates cross-over assurance and interference.

The chromosome axis plays a crucial role in meiotic recombination. Here, we study the function of ASY1, the Arabidopsis homolog of the yeast chromosome axis-associated component Hop1. Specifically, we characterized cross-over (CO) distribution in female and male meiosis by deep sequencing of the progeny of an allelic series of asy1 mutants. Combining data from nearly 1,000 individual plants, we find that reduced ASY1 functionality leads to genomic instability and sometimes drastic genomic rearrangements. We further observed that COs are less frequent and appear in more distal chromosomal regions in plants with no or reduced ASY1 functionality, consistent with previous analyses. However, our sequencing approach revealed that the reduction in CO number is not as dramatic as suggested by cytological analyses. Analysis of double mutants of asy1 with mutants with three other CO factors, MUS81, MSH4, and MSH5, as well as the determination of foci number of the CO regulator MLH1 demonstrates that the majority of the COs in asy1, similar to the situation in the wildtype (WT), largely belong to the class I, which are subject to interference. However, these COs are redistributed in asy1 mutants and typically appear much closer than in the WT. Hence, ASY1 plays a key role in CO interference that spaces COs along a chromosome. Conversely, since a large proportion of chromosomes do not receive any CO, we conclude that CO assurance, the process that ensures the obligatory assignment of one CO per chromosome, is also affected in asy1 mutants.

Chromoanagenesis in plants: triggers, mechanisms, and potential impact.

Chromoanagenesis is a single catastrophic event that involves, in most cases, localized chromosomal shattering and reorganization, resulting in a dramatically restructured chromosome. First discovered in cancer cells, it has since been observed in various other systems, including plants. In this review, we discuss the origin, characteristics, and potential mechanisms underlying chromoanagenesis in plants. We report that multiple processes, including mutagenesis and genetic engineering, can trigger chromoanagenesis via a variety of mechanisms such as micronucleation, breakage-fusion-bridge (BFB) cycles, or chain-like translocations. The resulting rearranged chromosomes can be preserved during subsequent plant growth, and sometimes inherited to the next generation. Because of their high tolerance to genome restructuring, plants offer a unique system for investigating the evolutionary consequences and potential practical applications of chromoanagenesis.

Chromoanagenesis in the asy1 meiotic mutant of Arabidopsis.

Chromoanagenesis is a catastrophic event that involves localized chromosomal shattering and reorganization. In this study, we report a case of chromoanagenesis resulting from defective meiosis in the MEIOTIC ASYNAPTIC MUTANT 1 (asy1) background in Arabidopsis thaliana. We provide a detailed characterization of the genomic structure of this individual with a severely shattered segment of chromosome 1. We identified 260 novel DNA junctions in the affected region, most of which affect gene sequence on 1 or both sides of the junction. Our results confirm that asy1-related defective meiosis is a potential trigger for chromoanagenesis. This is the first example of chromoanagenesis associated with female meiosis and indicates the potential for genome evolution during oogenesis. PLAIN LANGUAGE SUMMARY: Chromoanagenesis is a complex and catastrophic event that results in severely restructured chromosomes. It has been identified in cancer cells and in some plant samples, after specific triggering events. Here, we identified this kind of genome restructuring in a mutant that exhibits defective meiosis in the model plant system Arabidopsis thaliana.

Diploid mint (M. longifolia) can produce spearmint type oil with a high yield potential.

Mint oil is a key source of natural flavors with wide industrial applications. Two unbalanced polyploid cultivars named Native (Mentha Spicata L) and Scotch (M. × gracilis Sole) are the main producers of spearmint type oil, which is characterized by high levels of the monoterpenes (-)-carvone and (-)-limonene. These cultivars have been the backbone of spearmint oil production for decades, while breeding and improvement remained largely unexplored, in part, due to sterility in cultivated lines. Here we show that sexual breeding at the diploid level can be leveraged to develop new varieties that produce spearmint type oil, along with the improvement of other important traits. Using field trials and GC-FID oil analysis we characterized plant materials from a public germplasm repository and identified a diploid accession that exhibited 89.5% increase in oil yield, compared to the industry standard, and another that produces spearmint type oil. Spearmint-type oil was present at high frequency in a segregating F2 population (32/160) produced from these two accessions. Field-testing of ten of these F2 lines showed segregation for oil yield and confirmed the production of spearmint-type oil profiles. Two of these lines combined high yield and spearmint-type oil with acceptable analytic and sensory profiles. These results demonstrate that spearmint-type oil can be produced in a diploid background with high yield potential, providing a simpler genetic system for the development of improved spearmint varieties.

Chromoanagenesis from radiation-induced genome damage in Populus.

Chromoanagenesis is a genomic catastrophe that results in chromosomal shattering and reassembly. These extreme single chromosome events were first identified in cancer, and have since been observed in other systems, but have so far only been formally documented in plants in the context of haploid induction crosses. The frequency, origins, consequences, and evolutionary impact of such major chromosomal remodeling in other situations remain obscure. Here, we demonstrate the occurrence of chromoanagenesis in poplar (Populus sp.) trees produced from gamma-irradiated pollen. Specifically, in this population of siblings carrying indel mutations, two individuals exhibited highly frequent copy number variation (CNV) clustered on a single chromosome, one of the hallmarks of chromoanagenesis. Using short-read sequencing, we confirmed the presence of clustered segmental rearrangement. Independently, we identified and validated novel DNA junctions and confirmed that they were clustered and corresponded to these rearrangements. Our reconstruction of the novel sequences suggests that the chromosomal segments have reorganized randomly to produce a novel rearranged chromosome but that two different mechanisms might be at play. Our results indicate that gamma irradiation can trigger chromoanagenesis, suggesting that this may also occur when natural or induced mutagens cause DNA breaks. We further demonstrate that such events can be tolerated in poplar, and even replicated clonally, providing an attractive system for more in-depth investigations of their consequences.

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.

The persimmon genome reveals clues to the evolution of a lineage-specific sex determination system in plants.

Most angiosperms bear hermaphroditic flowers, but a few species have evolved outcrossing strategies, such as dioecy, the presence of separate male and female individuals. We previously investigated the mechanisms underlying dioecy in diploid persimmon (D. lotus) and found that male flowers are specified by repression of the autosomal gene MeGI by its paralog, the Y-encoded pseudo-gene OGI. This mechanism is thought to be lineage-specific, but its evolutionary path remains unknown. Here, we developed a full draft of the diploid persimmon genome (D. lotus), which revealed a lineage-specific whole-genome duplication event and provided information on the architecture of the Y chromosome. We also identified three paralogs, MeGI, OGI and newly identified Sister of MeGI (SiMeGI). Evolutionary analysis suggested that MeGI underwent adaptive evolution after the whole-genome duplication event. Transformation of tobacco plants with MeGI and SiMeGI revealed that MeGI specifically acquired a new function as a repressor of male organ development, while SiMeGI presumably maintained the original function. Later, a segmental duplication event spawned MeGI's regulator OGI on the Y-chromosome, completing the path leading to dioecy, and probably initiating the formation of the Y-chromosome. These findings exemplify how duplication events can provide flexible genetic material available to help respond to varying environments and provide interesting parallels for our understanding of the mechanisms underlying the transition into dieocy in plants.

A comprehensive genomic scan reveals gene dosage balance impacts on quantitative traits in Populus trees.

Gene dosage variation and the associated changes in gene expression influence a wide variety of traits, ranging from cancer in humans to yield in plants. It is also expected to affect important traits of ecological and agronomic importance in forest trees, but this variation has not been systematically characterized or exploited. Here we performed a comprehensive scan of the Populus genome for dosage-sensitive loci affecting quantitative trait variation for spring and fall phenology and biomass production. The study population was a large collection of clonally propagated F1 hybrid lines of Populus that saturate the genome 10-fold with deletions and insertions (indels) of known sizes and positions. As a group, the phenotypic means of the indel lines consistently differed from control nonindel lines, with an overall negative effect of both insertions and deletions on all biomass-related traits but more diverse effects and an overall wider phenotypic distribution of the indel lines for the phenology-related traits. We also investigated the correlation between gene dosage at specific chromosomal locations and phenotype, to identify dosage quantitative trait loci (dQTL). Such dQTL were detected for most phenotypes examined, but stronger effect dQTL were identified for the phenology-related traits than for the biomass traits. Our genome-wide screen for dosage sensitivity in a higher eukaryote demonstrates the importance of global genomic balance and the impact of dosage on life history traits.

One Hundred Ways to Invent the Sexes: Theoretical and Observed Paths to Dioecy in Plants.

Dioecy, the presence of male and female flowers on separate individuals, is both widespread and uncommon within flowering plants, with only a few percent of dioecious species spread across most major phylogenetic taxa. It is therefore safe to assume that dioecy evolved independently in these different groups, which allows us to ask questions regarding the molecular and developmental mechanisms underlying these independent transitions to dioecy. We start this review by examining the problem from the standpoint of a genetic engineer trying to develop dioecy, discuss various potential solutions, and compare them to models proposed in the past and based on genetic and evolutionary considerations. Next, we present recent information regarding candidate sex determinants in three species, acquired using newly established genomic approaches. Although such specific information is still scarce, it is slowly becoming apparent that various genes or pathways can be altered to evolve dioecy.

Significant enhancement of fatty acid composition in seeds of the allohexaploid, Camelina sativa, using CRISPR/Cas9 gene editing.

The CRISPR/Cas9 nuclease system is a powerful and flexible tool for genome editing, and novel applications of this system are being developed rapidly. Here, we used CRISPR/Cas9 to target the FAD2 gene in Arabidopsis thaliana and in the closely related emerging oil seed plant, Camelina sativa, with the goal of improving seed oil composition. We successfully obtained Camelina seeds in which oleic acid content was increased from 16% to over 50% of the fatty acid composition. These increases were associated with significant decreases in the less desirable polyunsaturated fatty acids, linoleic acid (i.e. a decrease from ~16% to <4%) and linolenic acid (a decrease from ~35% to <10%). These changes result in oils that are superior on multiple levels: they are healthier, more oxidatively stable and better suited for production of certain commercial chemicals, including biofuels. As expected, A. thaliana T2 and T3 generation seeds exhibiting these types of altered fatty acid profiles were homozygous for disrupted FAD2 alleles. In the allohexaploid, Camelina, guide RNAs were designed that simultaneously targeted all three homoeologous FAD2 genes. This strategy that significantly enhanced oil composition in T3 and T4 generation Camelina seeds was associated with a combination of germ-line mutations and somatic cell mutations in FAD2 genes in each of the three Camelina subgenomes.

Epigenetic Regulation of the Sex Determination Gene MeGI in Polyploid Persimmon.

Epigenetic regulation can add a flexible layer to genetic variation, potentially enabling long-term but reversible cis-regulatory changes to an allele while maintaining its DNA sequence. Here, we present a case in which alternative epigenetic states lead to reversible sex determination in the hexaploid persimmon Diospyros kaki Previously, we elucidated the molecular mechanism of sex determination in diploid persimmon and demonstrated the action of a Y-encoded sex determinant pseudogene called OGI, which produces small RNAs targeting the autosomal gene MeGI, resulting in separate male and female individuals (dioecy). We contrast these findings with the discovery, in hexaploid persimmon, of an additional layer of regulation in the form of DNA methylation of the MeGI promoter associated with the production of both male and female flowers in genetically male trees. Consistent with this model, developing male buds exhibited higher methylation levels across the MeGI promoter than developing female flowers from either monoecious or female trees. Additionally, a DNA methylation inhibitor induced developing male buds to form feminized flowers. Concurrently, in Y-chromosome-carrying trees, the expression of OGI is silenced by the presence of a SINE (short interspersed nuclear element)-like insertion in the OGI promoter. Our findings provide an example of an adaptive scenario involving epigenetic plasticity.

Catastrophic chromosomal restructuring during genome elimination in plants.

Genome instability is associated with mitotic errors and cancer. This phenomenon can lead to deleterious rearrangements, but also genetic novelty, and many questions regarding its genesis, fate and evolutionary role remain unanswered. Here, we describe extreme chromosomal restructuring during genome elimination, a process resulting from hybridization of Arabidopsis plants expressing different centromere histones H3. Shattered chromosomes are formed from the genome of the haploid inducer, consistent with genomic catastrophes affecting a single, laggard chromosome compartmentalized within a micronucleus. Analysis of breakpoint junctions implicates breaks followed by repair through non-homologous end joining (NHEJ) or stalled fork repair. Furthermore, mutation of required NHEJ factor DNA Ligase 4 results in enhanced haploid recovery. Lastly, heritability and stability of a rearranged chromosome suggest a potential for enduring genomic novelty. These findings provide a tractable, natural system towards investigating the causes and mechanisms of complex genomic rearrangements similar to those associated with several human disorders.

Efficient Genome-Wide Detection and Cataloging of EMS-Induced Mutations Using Exome Capture and Next-Generation Sequencing.

Chemical mutagenesis efficiently generates phenotypic variation in otherwise homogeneous genetic backgrounds, enabling functional analysis of genes. Advances in mutation detection have brought the utility of induced mutant populations on par with those produced by insertional mutagenesis, but systematic cataloguing of mutations would further increase their utility. We examined the suitability of multiplexed global exome capture and sequencing coupled with custom-developed bioinformatics tools to identify mutations in well-characterized mutant populations of rice (Oryza sativa) and wheat (Triticum aestivum). In rice, we identified ∼18,000 induced mutations from 72 independent M2 individuals. Functional evaluation indicated the recovery of potentially deleterious mutations for >2600 genes. We further observed that specific sequence and cytosine methylation patterns surrounding the targeted guanine residues strongly affect their probability to be alkylated by ethyl methanesulfonate. Application of these methods to six independent M2 lines of tetraploid wheat demonstrated that our bioinformatics pipeline is applicable to polyploids. In conclusion, we provide a method for developing large-scale induced mutation resources with relatively small investments that is applicable to resource-poor organisms. Furthermore, our results demonstrate that large libraries of sequenced mutations can be readily generated, providing enhanced opportunities to study gene function and assess the effect of sequence and chromatin context on mutations.

The BOY NAMED SUE quantitative trait locus confers increased meiotic stability to an adapted natural allopolyploid of Arabidopsis.

Whole-genome duplication resulting from polyploidy is ubiquitous in the evolutionary history of plant species. Yet, polyploids must overcome the meiotic challenge of pairing, recombining, and segregating more than two sets of chromosomes. Using genomic sequencing of synthetic and natural allopolyploids of Arabidopsis thaliana and Arabidopsis arenosa, we determined that dosage variation and chromosomal translocations consistent with homoeologous pairing were more frequent in the synthetic allopolyploids. To test the role of structural chromosomal differentiation versus genetic regulation of meiotic pairing, we performed sequenced-based, high-density genetic mapping in F2 hybrids between synthetic and natural lines. This F2 population displayed frequent dosage variation and deleterious homoeologous recombination. The genetic map derived from this population provided no indication of structural evolution of the genome of the natural allopolyploid Arabidopsis suecica, compared with its predicted parents. The F2 population displayed variation in meiotic regularity and pollen viability that correlated with a single quantitative trait locus, which we named BOY NAMED SUE, and whose beneficial allele was contributed by A. suecica. This demonstrates that an additive, gain-of-function allele contributes to meiotic stability and fertility in a recently established allopolyploid and provides an Arabidopsis system to decipher evolutionary and molecular mechanisms of meiotic regularity in polyploids.

Production of a high-efficiency TILLING population through polyploidization.

Targeting Induced Local Lesions in Genomes (TILLING) provides a nontransgenic method for reverse genetics that is widely applicable, even in species where other functional resources are missing or expensive to build. The efficiency of TILLING, however, is greatly facilitated by high mutation density. Species vary in the number of mutations induced by comparable mutagenic treatments, suggesting that genetic background may affect the response. Allopolyploid species have often yielded higher mutation density than diploids. To examine the effect of ploidy, we autotetraploidized the Arabidopsis (Arabidopsis thaliana) ecotype Columbia, whose diploid has been used for TILLING extensively, and mutagenized it with 50 mm ethylmethane sulfonate. While the same treatment sterilized diploid Columbia, the tetraploid M1 plants produced good seed. To determine the mutation density, we searched 528 individuals for induced mutations in 15 genes for which few or no knockout alleles were previously available. We constructed tridimensional pools from the genomic DNA of M2 plants, amplified target DNA, and subjected them to Illumina sequencing. The results were analyzed with an improved version of the mutation detection software CAMBa that accepts any pooling scheme. This small population provided a rich resource with approximately 25 mutations per queried 1.5-kb fragment, including on average four severe missense and 1.3 truncation mutations. The overall mutation density of 19.4 mutations Mb(-1) is 4 times that achieved in the corresponding diploid accession, indicating that genomic redundancy engenders tolerance to high mutation density. Polyploidization of diploids will allow the production of small populations, such as less than 2,000, that provide allelic series from knockout to mild loss of function for virtually all genes.

Global analysis of the small RNA transcriptome in different ploidies and genomic combinations of a vertebrate complex--the Squalius alburnoides.

The Squalius alburnoides complex (Steindachner) is one of the most intricate hybrid polyploid systems known in vertebrates. In this complex, the constant switch of the genome composition in consecutive generations, very frequently involving a change on the ploidy level, promotes repetitive situations of potential genomic shock. Previously in this complex, it was showed that in response to the increase in genome dosage, triploids hybrids could regulate gene expression to a diploid state. In this work we compared the small RNA profiles in the different genomic compositions interacting in the complex in order to explore the miRNA involvement in gene expression regulation of triploids. Using high-throughput arrays and sequencing technologies we were able to verify that diploid and triploid hybrids shared most of their sequences and their miRNA expression profiles were high correlated. However, an overall view indicates an up-regulation of several miRNAs in triploids and a global miRNA expression in triploids higher than the predicted from an additive model. Those results point to a participation of miRNAs in the cellular functional stability needed when the ploidy change.

Cis- and trans-regulatory divergence between progenitor species determines gene-expression novelty in Arabidopsis allopolyploids.

Gene-expression divergence between species shapes morphological evolution, but the molecular basis is largely unknown. Here we show cis- and trans-regulatory elements and chromatin modifications on gene-expression diversity in genetically tractable Arabidopsis allotetraploids. In Arabidopsis thaliana and Arabidopsis arenosa, both cis and trans with predominant cis-regulatory effects mediate gene-expression divergence. The majority of genes with both cis- and trans-effects are subjected to compensating interactions and stabilizing selection. Interestingly, cis- and trans-regulation is associated with chromatin modifications. In F1 allotetraploids, Arabidopsis arenosa trans factors predominately affect allelic expression divergence. Arabidopsis arenosa trans factors tend to upregulate Arabidopsis thaliana alleles, whereas Arabidopsis thaliana trans factors up- or down-regulate Arabidopsis arenosa alleles. In resynthesized and natural allotetraploids, trans effects drive expression of both homoeologous loci into the same direction. We provide evidence for natural selection and chromatin regulation in shaping gene-expression diversity during plant evolution and speciation.

Homoeolog-specific retention and use in allotetraploid Arabidopsis suecica depends on parent of origin and network partners.

BACKGROUND: Allotetraploids carry pairs of diverged homoeologs for most genes. With the genome doubled in size, the number of putative interactions is enormous. This poses challenges on how to coordinate the two disparate genomes, and creates opportunities by enhancing the phenotypic variation. New combinations of alleles co-adapt and respond to new environmental pressures. Three stages of the allopolyploidization process--parental species divergence, hybridization, and genome duplication--have been well analyzed. The last stage of evolutionary adjustments remains mysterious. RESULTS: Homoeolog-specific retention and use were analyzed in Arabidopsis suecica (As), a species derived from A. thaliana (At) and A. arenosa (Aa) in a single event 12,000 to 300,000 years ago. We used 405,466 diagnostic features on tiling microarrays to recognize At and Aa contributions to the As genome and transcriptome: 324 genes lacked Aa contributions and 614 genes lacked At contributions within As. In leaf tissues, 3,458 genes preferentially expressed At homoeologs while 4,150 favored Aa homoeologs. These patterns were validated with resequencing. Genes with preferential use of Aa homoeologs were enriched for expression functions, consistent with the dominance of Aa transcription. Heterologous networks--mixed from At and Aa transcripts--were underrepresented. CONCLUSIONS: Thousands of deleted and silenced homoeologs in the genome of As were identified. Since heterologous networks may be compromised by interspecies incompatibilities, these networks evolve co-biases, expressing either only Aa or only At homoeologs. This progressive change towards predominantly pure parental networks might contribute to phenotypic variability and plasticity, and enable the species to exploit a larger range of environments.