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.

A systems genetics approach to deciphering the effect of dosage variation on leaf morphology in Populus.

Gene copy number variation is frequent in plant genomes of various species, but the impact of such gene dosage variation on morphological traits is poorly understood. We used a large population of Populus carrying genomically characterized insertions and deletions across the genome to systematically assay the effect of gene dosage variation on a suite of leaf morphology traits. A systems genetics approach was used to integrate insertion and deletion locations, leaf morphology phenotypes, gene expression, and transcriptional network data, to provide an overview of how gene dosage influences morphology. Dosage-sensitive genomic regions were identified that influenced individual or pleiotropic morphological traits. We also identified cis-expression quantitative trait loci (QTL) within these dosage QTL regions, a subset of which modulated trans-expression QTL as well. Integration of data types within a gene co-expression framework identified co-expressed gene modules that are dosage sensitive, enriched for dosage expression QTL, and associated with morphological traits. Functional description of these modules linked dosage-sensitive morphological variation to specific cellular processes, as well as candidate regulatory genes. Together, these results show that gene dosage variation can influence morphological variation through complex changes in gene expression, and suggest that frequently occurring gene dosage variation has the potential to likewise influence quantitative traits in nature.

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.

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.

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

[This corrects the article DOI: 10.1371/journal.pgen.1008566.].

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.

A Y-Encoded Suppressor of Feminization Arose via Lineage-Specific Duplication of a Cytokinin Response Regulator in Kiwifruit.

Dioecy, the presence of male and female flowers on distinct individuals, has evolved independently in multiple plant lineages, and the genes involved in this differential development are just starting to be uncovered in a few species. Here, we used genomic approaches to investigate this pathway in kiwifruits (genus Actinidia). Genome-wide cataloging of male-specific subsequences, combined with transcriptome analysis, led to the identification of a type-C cytokinin response regulator as a potential sex determinant gene in this genus. Functional transgenic analyses in two model systems, Arabidopsis thaliana and Nicotiana tabacum, indicated that this gene acts as a dominant suppressor of carpel development, prompting us to name it Shy Girl (SyGI). Evolutionary analyses in a panel of Actinidia species revealed that SyGI is located in the Y-specific region of the genome and probably arose from a lineage-specific gene duplication. Comparisons with the duplicated autosomal counterpart, and with orthologs from other angiosperms, suggest that the SyGI-specific duplication and subsequent evolution of cis-elements may have played a key role in the acquisition of separate sexes in this species.

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.

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.

Rapid identification of lettuce seed germination mutants by bulked segregant analysis and whole genome sequencing.

Lettuce (Lactuca sativa) seeds exhibit thermoinhibition, or failure to complete germination when imbibed at warm temperatures. Chemical mutagenesis was employed to develop lettuce lines that exhibit germination thermotolerance. Two independent thermotolerant lettuce seed mutant lines, TG01 and TG10, were generated through ethyl methanesulfonate mutagenesis. Genetic and physiological analyses indicated that these two mutations were allelic and recessive. To identify the causal gene(s), we applied bulked segregant analysis by whole genome sequencing. For each mutant, bulked DNA samples of segregating thermotolerant (mutant) seeds were sequenced and analyzed for homozygous single-nucleotide polymorphisms. Two independent candidate mutations were identified at different physical positions in the zeaxanthin epoxidase gene (ABSCISIC ACID DEFICIENT 1/ZEAXANTHIN EPOXIDASE, or ABA1/ZEP) in TG01 and TG10. The mutation in TG01 caused an amino acid replacement, whereas the mutation in TG10 resulted in alternative mRNA splicing. Endogenous abscisic acid contents were reduced in both mutants, and expression of the ABA1 gene from wild-type lettuce under its own promoter fully complemented the TG01 mutant. Conventional genetic mapping confirmed that the causal mutations were located near the ZEP/ABA1 gene, but the bulked segregant whole genome sequencing approach more efficiently identified the specific gene responsible for the phenotype.

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.

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.

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.

A genome-wide analysis of Cas9 binding specificity using ChIP-seq and targeted sequence capture.

Clustered regularly interspaced short palindromic repeat (CRISPR) RNA-guided nucleases have gathered considerable excitement as a tool for genome engineering. However, questions remain about the specificity of target site recognition. Cleavage specificity is typically evaluated by low throughput assays (T7 endonuclease I assay, target amplification followed by high-throughput sequencing), which are limited to a subset of potential off-target sites. Here, we used ChIP-seq to examine genome-wide CRISPR binding specificity at gRNA-specific and gRNA-independent sites for two guide RNAs. RNA-guided Cas9 binding was highly specific to the target site while off-target binding occurred at much lower intensities. Cas9-bound regions were highly enriched in NGG sites, a sequence required for target site recognition by Streptococcus pyogenes Cas9. To determine the relationship between Cas9 binding and endonuclease activity, we applied targeted sequence capture, which allowed us to survey 1200 genomic loci simultaneously including potential off-target sites identified by ChIP-seq and by computational prediction. A high frequency of indels was observed at both target sites and one off-target site, while no cleavage activity could be detected at other ChIP-bound regions. Our results confirm the high-specificity of CRISPR endonucleases and demonstrate that sequence capture can be used as a high-throughput genome-wide approach to identify off-target activity.

Insights into the Prunus-Specific S-RNase-Based Self-Incompatibility System from a Genome-Wide Analysis of the Evolutionary Radiation of S Locus-Related F-box Genes.

Self-incompatibility (SI) is an important plant reproduction mechanism that facilitates the maintenance of genetic diversity within species. Three plant families, the Solanaceae, Rosaceae and Plantaginaceae, share an S-RNase-based gametophytic SI (GSI) system that involves a single S-RNase as the pistil S determinant and several F-box genes as pollen S determinants that act via non-self-recognition. Previous evidence has suggested a specific self-recognition mechanism in Prunus (Rosaceae), raising questions about the generality of the S-RNase-based GSI system. We investigated the evolution of the pollen S determinant by comparing the sequences of the Prunus S haplotype-specific F-box gene (SFB) with those of its orthologs in other angiosperm genomes. Our results indicate that the Prunus SFB does not cluster with the pollen S of other plants and diverged early after the establishment of the Eudicots. Our results further indicate multiple F-box gene duplication events, specifically in the Rosaceae family, and suggest that the Prunus SFB gene originated in a recent Prunus-specific gene duplication event. Transcriptomic and evolutionary analyses of the Prunus S paralogs are consistent with the establishment of a Prunus-specific SI system, and the possibility of subfunctionalization differentiating the newly generated SFB from the original pollen S determinant.

Highly active zinc-finger nucleases by extended modular assembly.

Zinc-finger nucleases (ZFNs) are important tools for genome engineering. Despite intense interest by many academic groups, the lack of robust noncommercial methods has hindered their widespread use. The modular assembly (MA) of ZFNs from publicly available one-finger archives provides a rapid method to create proteins that can recognize a very broad spectrum of DNA sequences. However, three- and four-finger arrays often fail to produce active nucleases. Efforts to improve the specificity of the one-finger archives have not increased the success rate above 25%, suggesting that the MA method might be inherently inefficient due to its insensitivity to context-dependent effects. Here we present the first systematic study on the effect of array length on ZFN activity. ZFNs composed of six-finger MA arrays produced mutations at 15 of 21 (71%) targeted loci in human and mouse cells. A novel drop-out linker scheme was used to rapidly assess three- to six-finger combinations, demonstrating that shorter arrays could improve activity in some cases. Analysis of 268 array variants revealed that half of MA ZFNs of any array composition that exceed an ab initio B-score cutoff of 15 were active. These results suggest that, when used appropriately, MA ZFNs are able to target more DNA sequences with higher success rates than other current methods.

Sex- and melanism-specific variations in the oxidative status of adult tawny owls in response to manipulated reproductive effort.

Oxidative stress, determined by the balance between the production of damaging reactive oxygen species (ROS) and antioxidant defences, is hypothesized to play an important role in shaping the cost of reproduction and life history trade-offs. To test this hypothesis, we manipulated reproductive effort in 94 breeding pairs of tawny owls (Strix aluco) to investigate the sex- and melanism-specific effects on markers of oxidative stress in red blood cells (RBCs). This colour polymorphic bird species shows sex-specific division of labour and melanism-specific history strategies. Brood sizes at hatching were experimentally enlarged or reduced to increase or decrease reproductive effort, respectively. We obtained an integrative measure of the oxidative balance by measuring ROS production by RBCs, intracellular antioxidant glutathione levels and membrane resistance to ROS. We found that light melanic males (the sex undertaking offspring food provisioning) produced more ROS than darker conspecifics, but only when rearing an enlarged brood. In both sexes, light melanic individuals had also a larger pool of intracellular antioxidant glutathione than darker owls under relaxed reproductive conditions (i.e. reduced brood), but not when investing substantial effort in current reproduction (enlarged brood). Finally, resistance to oxidative stress was differently affected by the brood size manipulation experiment in males and females independently of their plumage coloration. Altogether, our results support the hypothesis that reproductive effort can alter the oxidative balance in a sex- and colour-specific way. This further emphasizes the close link between melanin-based coloration and life history strategies.

Rapid creation of Arabidopsis doubled haploid lines for quantitative trait locus mapping.

Quantitative trait loci (QTL) mapping is a powerful tool for investigating the genetic basis of natural variation. QTL can be mapped using a number of different population designs, but recombinant inbred lines (RILs) are among the most effective. Unfortunately, homozygous RIL populations are time consuming to construct, typically requiring at least six generations of selfing starting from a heterozygous F(1). Haploid plants produced from an F(1) combine the two parental genomes and have only one allele at every locus. Converting these sterile haploids into fertile diploids (termed "doubled haploids," DHs) produces immortal homozygous lines in only two steps. Here we describe a unique technique for rapidly creating recombinant doubled haploid populations in Arabidopsis thaliana: centromere-mediated genome elimination. We generated a population of 238 doubled haploid lines that combine two parental genomes and genotyped them by reduced representation Illumina sequencing. The recombination rate and parental allele frequencies in our population are similar to those found in existing RIL sets. We phenotyped this population for traits related to flowering time and for petiole length and successfully mapped QTL controlling each trait. Our work demonstrates that doubled haploid populations offer a rapid, easy alternative to RILs for Arabidopsis genetic analysis.

Melanin-specific life-history strategies.

The maintenance of genetic variation is a long-standing issue because the adaptive value of life-history strategies associated with each genetic variant is usually unknown. However, evidence for the coexistence of alternative evolutionary fixed strategies at the population level remains scarce. Because in the tawny owl (Strix aluco) heritable melanin-based coloration shows different physiological and behavioral norms of reaction, we investigated whether coloration is associated with investment in maintenance and reproduction. Light melanic owls had lower adult survival compared to dark melanic conspecifics, and color variation was related to the trade-off between offspring number and quality. When we experimentally enlarged brood size, light melanic males produced more fledglings but in poorer condition, and they were less often recruited in the local breeding population than those of darker melanic conspecifics. Our results also suggest that dark melanic males allocate a constant effort to raise their brood independently of environmental conditions, whereas lighter melanic males finely adjust reproductive effort in relation to changes in environmental conditions. Color traits can therefore be associated with life-history strategies, and stochastic environmental perturbation can temporarily favor one phenotype over others. The existence of fixed strategies implies that some phenotypes can sometimes display a "maladapted" strategy. Long-term population monitoring is therefore vital for a full understanding of how different genotypes deal with trade-offs.