Molecular, genetic and evolutionary analysis of a paracentric inversion in Arabidopsis thaliana

Epigenetic state and gene expression remain stable after CRISPR/Cas-mediated chromosomal inversions

The epigenetic state of chromatin, gene activity and chromosomal positions are interrelated in plants. In Arabidopsis thaliana, chromosome arms are DNA-hypomethylated and enriched with the euchromatin-specific histone mark H3K4me3, while pericentromeric regions are DNA-hypermethylated and enriched with the heterochromatin-specific mark H3K9me2. We aimed to investigate how the chromosomal location affects epigenetic stability and gene expression by chromosome engineering. Two chromosomal inversions of different sizes were induced using CRISPR/Cas9 to move heterochromatic, pericentric sequences into euchromatic regions. The epigenetic status of these lines was investigated using whole-genome bisulfite sequencing and chromatin immunoprecipitation. Gene expression changes following the induction of the chromosomal inversions were studied via transcriptome analysis. Both inversions had a minimal impact on the global distribution of histone marks and DNA methylation patterns, although minor epigenetic changes were observed across the genome. Notably, the inverted chromosomal regions and their borders retained their original epigenetic profiles. Gene expression analysis showed that only 0.5-1% of genes were differentially expressed genome-wide following the induction of the inversions. CRISPR/Cas-induced chromosomal inversions minimally affect epigenetic landscape and gene expression, preserving their profiles in subsequent generations.

Chromosome engineering points to the cis-acting mechanism of chromosome arm-specific telomere length setting and robustness of plant phenotype, chromatin structure and gene expression

The study investigates the impact of targeted chromosome engineering on telomere dynamics, chromatin structure, gene expression, and phenotypic stability in Arabidopsis thaliana. Using precise CRISPR/Cas-based engineering, reciprocal translocations of chromosome arms were introduced between non-homologous chromosomes. The subsequent homozygous generations of plants were assessed for phenotype, transcriptomic changes and chromatin modifications near translocation breakpoints, and telomere length maintenance. Phenotypically, translocated lines were indistinguishable from wild-type plants, as confirmed through morphological assessments and principal component analysis. Gene expression profiling detected minimal differential expression, with affected genes dispersed across the genome, indicating negligible transcriptional impact. Similarly, ChIPseq analysis showed no substantial alterations in the enrichment of key histone marks (H3K27me3, H3K4me1, H3K56ac) near junction sites or across the genome. Finally, bulk and arm-specific telomere lengths remained stable across multiple generations, except for minor variations in one translocation line. These findings highlight the remarkable genomic and phenotypic robustness of A. thaliana despite large-scale chromosomal rearrangements. The study offers insights into the cis-acting mechanisms underlying chromosome arm-specific telomere length setting and establishes the feasibility of chromosome engineering for studies of plant genome evolution and crop improvement strategies.

A catalogue of recombination coldspots in interspecific tomato hybrids

Increasing natural resistance and resilience in plants is key for ensuring food security within a changing climate. Breeders improve these traits by crossing cultivars with their wild relatives and introgressing specific alleles through meiotic recombination. However, some genomic regions are devoid of recombination especially in crosses between divergent genomes, limiting the combinations of desirable alleles. Here, we used pooled-pollen sequencing to build a map of recombinant and non-recombinant regions between tomato and five wild relatives commonly used for introgressive tomato breeding. We detected hybrid-specific recombination coldspots that underscore the role of structural variations in modifying recombination patterns and maintaining genetic linkage in interspecific crosses. Crossover regions and coldspots show strong association with specific TE superfamilies exhibiting differentially accessible chromatin between somatic and meiotic cells. About two-thirds of the genome are conserved coldspots, located mostly in the pericentromeres and enriched with retrotransposons. The coldspots also harbor genes associated with agronomic traits and stress resistance, revealing undesired consequences of linkage drag and possible barriers to breeding. We presented examples of linkage drag that can potentially be resolved by pairing tomato with other wild species. Overall, this catalogue will help breeders better understand crossover localization and make informed decisions on generating new tomato varieties.

Maintenance and dynamic reprogramming of chromatin organization during development

Controlled transcription of genes is critical for cell differentiation and development. Gene expression regulation therefore involves a multilayered control from nucleosome composition in histone variants and their post-translational modifications to higher-order folding of chromatin fibers and chromatin interactions in nuclear space. Recent technological advances have allowed gaining insight into these mechanisms, the interplay between local and higher-order chromatin organization, and the dynamic changes that occur during stress response and developmental transitions. In this review, we will discuss chromatin organization from the nucleosome to its three-dimensional structure in the nucleus, and consider how these different layers of organization are maintained during the cell cycle or rapidly reprogrammed during development.

A pan-genome of 69 Arabidopsis thaliana accessions reveals a conserved genome structure throughout the global species range

Although originally primarily a system for functional biology, Arabidopsis thaliana has, owing to its broad geographical distribution and adaptation to diverse environments, developed into a powerful model in population genomics. Here we present chromosome-level genome assemblies of 69 accessions from a global species range. We found that genomic colinearity is very conserved, even among geographically and genetically distant accessions. Along chromosome arms, megabase-scale rearrangements are rare and typically present only in a single accession. This indicates that the karyotype is quasi-fixed and that rearrangements in chromosome arms are counter-selected. Centromeric regions display higher structural dynamics, and divergences in core centromeres account for most of the genome size variations. Pan-genome analyses uncovered 32,986 distinct gene families, 60% being present in all accessions and 40% appearing to be dispensable, including 18% private to a single accession, indicating unexplored genic diversity. These 69 new Arabidopsis thaliana genome assemblies will empower future genetic research.

Two large inversions seriously suppress recombination and are essential for key genotype fixation in cabbage (Brassica oleracea L. var. capitata)

Chromosomal inversion is an important structural variation that usually suppresses recombination and is critical for key genotype fixation. In a previous study, an 11.47 Mb recombination suppression region was identified in the yellow-green leaf locus BoYgl-1 on chromosome 1, but the cause of recombination suppression is still unclear. In this study, chlorophyll and carotenoid contents were found to be significantly decreased in the yellow-green leaf mutant YL-1. Genome assembly and comparative analysis revealed that two large inversions in YL-1 were responsible for the severe recombination suppression in the BoYgl-1 locus. Analyses with inversion-specific markers revealed that the inversions were present in 44 (including all wild cabbage; INV1 and INV2) of 195 cabbage inbred lines and 15 (INV1) ornamental kale inbred lines, indicating that these species with INV1 or INV2 may have evolved much earlier than other types of cabbage. Analyses with inversion-correlated markers revealed that the genotypes of CoINV1, CoINV2 and CoINV3 were highly correlated with INV1 and INV2, indicating that INVs could fix the key genotypes of the involved region. In addition, a 5.87 Mb assembly inversion was identified at the BoYgl-1 locus in the TO1000 genome by genome comparative analysis. This study provides new insight into the recombination suppression mechanism of chromosomal inversion and the application of genome fragment fixation in cabbage breeding.

Plant chromosome engineering - past, present and future

Spontaneous chromosomal rearrangements (CRs) play an essential role in speciation, genome evolution and crop domestication. To be able to use the potential of CRs for breeding, plant chromosome engineering was initiated by fragmenting chromosomes by X-ray irradiation. With the rise of the CRISPR/Cas system, it became possible to induce double-strand breaks (DSBs) in a highly efficient manner at will at any chromosomal position. This has enabled a completely new level of predesigned chromosome engineering. The genetic linkage between specific genes can be broken by inducing chromosomal translocations. Natural inversions, which suppress genetic exchange, can be reverted for breeding. In addition, various approaches for constructing minichromosomes by downsizing regular standard A or supernumerary B chromosomes, which could serve as future vectors in plant biotechnology, have been developed. Recently, a functional synthetic centromere could be constructed. Also, different ways of genome haploidization have been set up, some based on centromere manipulations. In the future, we expect to see even more complex rearrangements, which can be combined with previously developed engineering technologies such as recombinases. Chromosome engineering might help to redefine genetic linkage groups, change the number of chromosomes, stack beneficial genes on mini cargo chromosomes, or set up genetic isolation to avoid outcrossing.

Comparative analysis of repetitive DNA in dysploid and non-dysploid Phaseolus beans

Structural karyotype changes result from ectopic recombination events frequently associated with repetitive DNA. Although most Phaseolus species present relatively stable karyotypes with 2n = 22 chromosomes, the karyotypes of species of the Leptostachyus group show high rates of structural rearrangements, including a nested chromosome fusion that led to the dysploid chromosome number of the group (2n = 20). We examined the roles of repetitive landscapes in the rearrangements of species of the Leptostachyus group using genome-skimming data to characterize the repeatome in a range of Phaseolus species and compared them to species of that group (P. leptostachyus and P. macvaughii). LTR retrotransposons, especially the Ty3/gypsy lineage Chromovirus, were the most abundant elements in the genomes. Differences in the abundance of Tekay, Retand, and SIRE elements between P. macvaughii and P. leptostachyus were reflected in their total amounts of Ty3/gypsy and Ty1/copia. The satellite DNA fraction was the most divergent among the species, varying both in abundance and distribution, even between P. leptostachyus and P. macvaughii. The rapid turnover of repeats in the Leptostachyus group may be associated with the several rearrangements observed.

Transcriptome-wide association study coupled with eQTL analysis reveals the genetic connection between gene expression and flowering time in Arabidopsis

Genome-wide association study (GWAS) has improved our understanding of complex traits, but challenges exist in distinguishing causation versus association caused by linkage disequilibrium. Instead, transcriptome-wide association studies (TWAS) detect direct associations between expression levels and phenotypic variations, providing an opportunity to better prioritize candidate genes. To assess the feasibility of TWAS, we investigated the association between transcriptomes, genomes, and various traits in Arabidopsis, including flowering time. The associated genes formerly known to regulate growth allometry or metabolite production were first identified by TWAS. Next, for flowering time, six TWAS-newly identified genes were functionally validated. Analysis of the expression quantitative trait locus (eQTL) further revealed a trans-regulatory hotspot affecting the expression of several TWAS-identified genes. The hotspot covers the FRIGIDA (FRI) gene body, which possesses multiple haplotypes differentially affecting the expression of downstream genes, such as FLOWERING LOCUS C (FLC) and SUPPRESSOR OF OVEREXPRESSION OF CO 1 (SOC1). We also revealed multiple independent paths towards the loss of function of FRI in natural accessions. Altogether, this study demonstrates the potential of combining TWAS with eQTL analysis to identify important regulatory modules of FRI-FLC-SOC1 for quantitative traits in natural populations.

High rates of structural rearrangements have shaped the chromosome evolution in dysploid Phaseolus beans

KEY MESSAGE

Karyotypes evolve through numerical and structural chromosome rearrangements. We show that Phaseolus leptostachyus, a wild bean, underwent a rapid genome reshuffling associated with the reduction from 11 to 10 chromosome pairs, but without whole genome duplication, the highest chromosome evolution rate known for plants. Plant karyotypes evolve through structural rearrangements often associated with polyploidy or dysploidy. The genus Phaseolus comprises ~ 90 species, five of them domesticated due to their nutritional relevance. Most of the species have 2n = 22 karyotypes and are highly syntenic, except for three dysploid karyotypes of species from the Leptostachyus group (2n = 20) that have accumulated several rearrangements. Here, we investigated the degrees of structural rearrangements among Leptostachyus and other Phaseolus groups by estimating their chromosomal evolution rates (CER). For this, we combined our oligo-FISH barcode system for beans and chromosome-specific painting probes for chromosomes 2 and 3, with rDNA and a centromeric probe to establish chromosome orthologies and identify structural rearrangements across nine Phaseolus species. We also integrated the detected rearrangements with a phylogenomic approach to estimate the CERs for each Phaseolus lineage. Our data allowed us to identify translocations, inversions, duplications and deletions, mostly in species belonging to the Leptostachyus group. Phaseolus leptostachyus showed the highest CER (12.31 rearrangements/My), a tenfold increase in contrast to the 2n = 22 species analysed. This is the highest rate known yet for plants, making it a model species for investigating the mechanisms behind rapid genome reshuffling in early species diversification.

Genomics of plant speciation

Studies of plants have been instrumental for revealing how new species originate. For several decades, botanical research has complemented and, in some cases, challenged concepts on speciation developed via the study of other organisms while also revealing additional ways in which species can form. Now, the ability to sequence genomes at an unprecedented pace and scale has allowed biologists to settle decades-long debates and tackle other emerging challenges in speciation research. Here, we review these recent genome-enabled developments in plant speciation. We discuss complications related to identification of reproductive isolation (RI) loci using analyses of the landscape of genomic divergence and highlight the important role that structural variants have in speciation, as increasingly revealed by new sequencing technologies. Further, we review how genomics has advanced what we know of some routes to new species formation, like hybridization or whole-genome duplication, while casting doubt on others, like population bottlenecks and genetic drift. While genomics can fast-track identification of genes and mutations that confer RI, we emphasize that follow-up molecular and field experiments remain critical. Nonetheless, genomics has clarified the outsized role of ancient variants rather than new mutations, particularly early during speciation. We conclude by highlighting promising avenues of future study. These include expanding what we know so far about the role of epigenetic and structural changes during speciation, broadening the scope and taxonomic breadth of plant speciation genomics studies, and synthesizing information from extensive genomic data that have already been generated by the plant speciation community.

Five Regions of the Pea Genome Co-Control Partial Resistance to D. pinodes, Tolerance to Frost, and Some Architectural or Phenological Traits

Evidence for reciprocal links between plant responses to biotic or abiotic stresses and architectural and developmental traits has been raised using approaches based on epidemiology, physiology, or genetics. Winter pea has been selected for years for many agronomic traits contributing to yield, taking into account architectural or phenological traits such as height or flowering date. It remains nevertheless particularly susceptible to biotic and abiotic stresses, among which Didymella pinodes and frost are leading examples. The purpose of this study was to identify and resize QTL localizations that control partial resistance to D. pinodes, tolerance to frost, and architectural or phenological traits on pea dense genetic maps, considering how QTL colocalizations may impact future winter pea breeding. QTL analysis revealed five metaQTLs distributed over three linkage groups contributing to both D. pinodes disease severity and frost tolerance. At these loci, the haplotypes of alleles increasing both partial resistance to D. pinodes and frost tolerance also delayed the flowering date, increased the number of branches, and/or decreased the stipule length. These results question both the underlying mechanisms of the joint control of biotic stress resistance, abiotic stress tolerance, and plant architecture and phenology and the methods of marker-assisted selection optimizing stress control and productivity in winter pea breeding.

Comparison and Characterization of Phenotypic and Genomic Mutations Induced by a Carbon-Ion Beam and Gamma-ray Irradiation in Soybean (Glycine max (L.) Merr.)

Soybean (Glycine max (L.) Merr.) is a nutritious crop that can provide both oil and protein. A variety of mutagenesis methods have been proposed to obtain better soybean germplasm resources. Among the different types of physical mutagens, carbon-ion beams are considered to be highly efficient with high linear energy transfer (LET), and gamma rays have also been widely used for mutation breeding. However, systematic knowledge of the mutagenic effects of these two mutagens during development and on phenotypic and genomic mutations has not yet been elucidated in soybean. To this end, dry seeds of Williams 82 soybean were irradiated with a carbon-ion beam and gamma rays. The biological effects of the M₁ generation included changes in survival rate, yield and fertility. Compared with gamma rays, the relative biological effectiveness (RBE) of the carbon-ion beams was between 2.5 and 3.0. Furthermore, the optimal dose for soybean was determined to be 101 Gy to 115 Gy when using the carbon-ion beam, and it was 263 Gy to 343 Gy when using gamma rays. A total of 325 screened mutant families were detected from out of 2000 M₂ families using the carbon-ion beam, and 336 screened mutant families were found using gamma rays. Regarding the screened phenotypic M₂ mutations, the proportion of low-frequency phenotypic mutations was 23.4% when using a carbon ion beam, and the proportion was 9.8% when using gamma rays. Low-frequency phenotypic mutations were easily obtained with the carbon-ion beam. After screening the mutations from the M₂ generation, their stability was verified, and the genome mutation spectrum of M₃ was systemically profiled. A variety of mutations, including single-base substitutions (SBSs), insertion-deletion mutations (INDELs), multinucleotide variants (MNVs) and structural variants (SVs) were detected with both carbon-ion beam irradiation and gamma-ray irradiation. Overall, 1988 homozygous mutations and 9695 homozygous + heterozygous genotype mutations were detected when using the carbon-ion beam. Additionally, 5279 homozygous mutations and 14,243 homozygous + heterozygous genotype mutations were detected when using gamma rays. The carbon-ion beam, which resulted in low levels of background mutations, has the potential to alleviate the problems caused by linkage drag in soybean mutation breeding. Regarding the genomic mutations, when using the carbon-ion beam, the proportion of homozygous-genotype SVs was 0.45%, and that of homozygous + heterozygous-genotype SVs was 6.27%; meanwhile, the proportions were 0.04% and 4.04% when using gamma rays. A higher proportion of SVs were detected when using the carbon ion beam. The gene effects of missense mutations were greater under carbon-ion beam irradiation, and the gene effects of nonsense mutations were greater under gamma-ray irradiation, which meant that the changes in the amino acid sequences were different between the carbon-ion beam and gamma rays. Taken together, our results demonstrate that both carbon-ion beam and gamma rays are effective techniques for rapid mutation breeding in soybean. If one would like to obtain mutations with a low-frequency phenotype, low levels of background genomic mutations and mutations with a higher proportion of SVs, carbon-ion beams are the best choice.

Identification of structural variations related to drought tolerance in wheat (Triticum aestivum L.)

Structural variations are common in plant genomes, affecting meiotic recombination and distorted segregation in wheat. And presence/absence variations can significantly affect drought tolerance in wheat. Drought is a major abiotic stress limiting wheat production. Common wheat has a complex genome with three sub-genomes, which host large numbers of structural variations (SVs). SVs play critical roles in understanding the genetic contributions of plant domestication and phenotypic plasticity, but little is known about their genomic characteristics and their effects on drought tolerance. In the present study, high-resolution karyotypes of 180 doubled haploids (DHs) were developed. Signal polymorphisms between the parents involved with 8 presence-absence variations (PAVs) of tandem repeats (TR) distributed on the 7 (2A, 4A, 5A, 7A, 3B, 7B, and 2D) of 21 chromosomes. Among them, PAV on chromosome 2D showed distorted segregation, others transmit normal conforming to a 1:1 segregation ration in the population; and a PAVs recombination occurred on chromosome 2A. Association analysis of PAV and phenotypic traits under different water regimes, we found PAVs on chromosomes 4A, 5A, and 7B showed negative effect on grain length (GL) and grain width (GW); PAV.7A had opposite effect on grain thickness (GT) and spike length (SL), with the effect on traits differing under different water regimes. PAVs on linkage group 2A, 4A, 7A, 2D, and 7B associated with the drought tolerance coefficients (DTCs), and significant negative effect on drought resistance values (D values) were detected in PAV.7B. Additionally, quantitative trait loci (QTL) associated with phenotypic traits using the 90 K SNP array showed QTL for DTCs and grain-related traits in chromosomes 4A, and 5A, 3B were co-localized in differential regions of PAVs. These PAVs can cause the differentiation of the target region of SNP and could be used for genetic improvement of agronomic traits under drought stress via marker-assisted selection (MAS) breeding.

Genome-wide association study reveals structural chromosome variations with phenotypic effects in wheat (Triticum aestivum L.)

Structural chromosome variations (SCVs) are large-scale genomic variations that can be detected by fluorescence in situ hybridization (FISH). SCVs have played important roles in the genome evolution of wheat (Triticum aestivum L.), but little is known about their genetic effects. In this study, a total of 543 wheat accessions from the Chinese wheat mini-core collection and the Shanxi Province wheat collection were used for chromosome analysis using oligonucleotide probe multiplex FISH. A total of 139 SCVs including translocations, pericentric inversions, presence/absence variations (PAVs), and copy number variations (CNVs) in heterochromatin were identified at 230 loci. The distribution frequency of SCVs varied between ecological regions and between landraces and modern cultivars. Structural analysis using SCVs as markers clearly divided the landraces and modern cultivars into different groups. There are very clear instances illustrating alien introgression and wide application of foreign germplasms improved the chromosome diversity of Chinese modern wheat cultivars. A genome-wide association study (GWAS) identified 29 SCVs associated with 12 phenotypic traits, and five (RT4AS•4AL-1DS/1DL•1DS-4AL, Mg2A-3, Mr3B-10, Mr7B-13, and Mr4A-7) of them were further validated using a doubled haploid population and advanced sib-lines, implying the potential value of these SCVs. Importantly, the number of favored SCVs that were associated with agronomic trait improvement was significantly higher in modern cultivars compared to landraces, indicating positive selection in wheat breeding. This study demonstrates the significant effects of SCVs during wheat breeding and provides an efficient method of mining favored SCVs in wheat and other crops.

Massive crossover suppression by CRISPR-Cas-mediated plant chromosome engineering

Recent studies have demonstrated that not only genes but also entire chromosomes can be engineered using clustered regularly interspaced short palindromic repeats (CRISPR)-CRISPER-associated protein 9 (Cas9)^1-5. A major objective of applying chromosome restructuring in plant breeding is the manipulation of genetic exchange⁶. Here we show that meiotic recombination can be suppressed in nearly the entire chromosome using chromosome restructuring. We were able to induce a heritable inversion of a >17 Mb-long chromosome fragment that contained the centromere and covered most of chromosome 2 of the Arabidopsis ecotype Col-0. Only the 2 and 0.5 Mb-long telomeric ends remained in their original orientation. In single-nucleotide polymorphism marker analysis of the offspring of crosses with the ecotype Ler-1, we detected a massive reduction of crossovers within the inverted chromosome region, coupled with a shift of crossovers to the telomeric ends. The few genetic exchanges detected within the inversion all originated from double crossovers. This not only indicates that heritable genetic exchange can occur by interstitial chromosome pairing, but also that it is restricted to the production of viable progeny.

Comparative repeatome analysis reveals new evidence on genome evolution in wild diploid Arachis (Fabaceae) species

Opposing changes in the abundance of satellite DNA and long terminal repeat (LTR) retroelements are the main contributors to the variation in genome size and heterochromatin amount in Arachis diploids. The South American genus Arachis (Fabaceae) comprises 83 species organized in nine taxonomic sections. Among them, section Arachis is characterized by species with a wide genome and karyotype diversity. Such diversity is determined mainly by the amount and composition of repetitive DNA. Here we performed computational analysis on low coverage genome sequencing to infer the dynamics of changes in major repeat families that led to the differentiation of genomes in diploid species (x = 10) of genus Arachis, focusing on section Arachis. Estimated repeat content ranged from 62.50 to 71.68% of the genomes. Species with different genome composition tended to have different landscapes of repeated sequences. Athila family retrotransposons were the most abundant and variable lineage among Arachis repeatomes, with peaks of transpositional activity inferred at different times in the evolution of the species. Satellite DNAs (satDNAs) were less abundant, but differentially represented among species. High rates of evolution of an AT-rich superfamily of satDNAs led to the differential accumulation of heterochromatin in Arachis genomes. The relationship between genome size variation and the repetitive content is complex. However, largest genomes presented a higher accumulation of LTR elements and lower contents of satDNAs. In contrast, species with lowest genome sizes tended to accumulate satDNAs in detriment of LTR elements. Phylogenetic analysis based on repetitive DNA supported the genome arrangement of section Arachis. Altogether, our results provide the most comprehensive picture on the repeatome dynamics that led to the genome differentiation of Arachis species.

Epigenomic and structural events preclude recombination in Brassica napus

Meiotic recombination is a major evolutionary process generating genetic diversity at each generation in sexual organisms. However, this process is highly regulated, with the majority of crossovers lying in the distal chromosomal regions that harbor low DNA methylation levels. Even in these regions, some islands without recombination remain, for which we investigated the underlying causes. Genetic maps were established in two Brassica napus hybrids to detect the presence of such large nonrecombinant islands. The role played by DNA methylation and structural variations in this local absence of recombination was determined by performing bisulfite sequencing and whole genome comparisons. Inferred structural variations were validated using either optical mapping or oligo fluorescence in situ hybridization. Hypermethylated or inverted regions between Brassica genomes were associated with the absence of recombination. Pairwise comparisons of nine B. napus genome assemblies revealed that such inversions occur frequently and may contain key agronomic genes such as resistance to biotic stresses. We conclude that such islands without recombination can have different origins, such as DNA methylation or structural variations in B. napus. It is thus essential to take into account these features in breeding programs as they may hamper the efficient combination of favorable alleles in elite varieties.

Nonhomologous end joining as key to CRISPR/Cas-mediated plant chromosome engineering

Although clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein (Cas)-mediated gene editing has revolutionized biology and plant breeding, large-scale, heritable restructuring of plant chromosomes is still in its infancy. Duplications and inversions within a chromosome, and also translocations between chromosomes, can now be achieved. Subsequently, genetic linkages can be broken or can be newly created. Also, the order of genes on a chromosome can be changed. While natural chromosomal recombination occurs by homologous recombination during meiosis, CRISPR/Cas-mediated chromosomal rearrangements can be obtained best by harnessing nonhomologous end joining (NHEJ) pathways in somatic cells. NHEJ can be subdivided into the classical (cNHEJ) and alternative NHEJ (aNHEJ) pathways, which partially operate antagonistically. The cNHEJ pathway not only protects broken DNA ends from degradation but also suppresses the joining of previously unlinked broken ends. Hence, in the absence of cNHEJ, more inversions or translocations can be obtained which can be ascribed to the unrestricted use of the aNHEJ pathway for double-strand break (DSB) repair. In contrast to inversions or translocations, short tandem duplications can be produced by paired single-strand breaks via a Cas9 nickase. Interestingly, the cNHEJ pathway is essential for these kinds of duplications, whereas aNHEJ is required for patch insertions that can also be formed during DSB repair. As chromosome engineering has not only been accomplished in the model plant Arabidopsis (Arabidopsis thaliana) but also in the crop maize (Zea mays), we expect that this technology will soon transform the breeding process.

Domestication shapes recombination patterns in tomato

Meiotic recombination is a biological process of key importance in breeding, to generate genetic diversity and develop novel or agronomically relevant haplotypes. In crop tomato, recombination is curtailed as manifested by linkage disequilibrium decay over a longer distance and reduced diversity compared with wild relatives. Here, we compared domesticated and wild populations of tomato and found an overall conserved recombination landscape, with local changes in effective recombination rate in specific genomic regions. We also studied the dynamics of recombination hotspots resulting from domestication and found that loss of such hotspots is associated with selective sweeps, most notably in the pericentromeric heterochromatin. We detected footprints of genetic changes and structural variants, among them associated with transposable elements, linked with hotspot divergence during domestication, likely causing fine-scale alterations to recombination patterns and resulting in linkage drag.

Transposon-induced inversions activate gene expression in the maize pericarp

Chromosomal inversions can have considerable biological and agronomic impacts including disrupted gene function, change in gene expression, and inhibited recombination. Here, we describe the molecular structure and functional impact of six inversions caused by Alternative Transpositions between p1 and p2 genes responsible for floral pigmentation in maize. In maize line p1-wwB54, the p1 gene is null and the p2 gene is expressed in anther and silk but not in pericarp, making the kernels white. By screening for kernels with red pericarp, we identified inversions in this region caused by transposition of Ac and fractured Ac (fAc) transposable elements. We hypothesize that these inversions place the p2 gene promoter near a p1 gene enhancer, thereby activating p2 expression in kernel pericarp. To our knowledge, this is the first report of multiple recurrent inversions that change the position of a gene promoter relative to an enhancer to induce ectopic expression in a eukaryote.

CRISPR-Cas-mediated chromosome engineering for crop improvement and synthetic biology

Plant breeding relies on the presence of genetic variation, as well as on the ability to break or stabilize genetic linkages between traits. The development of the genome-editing tool clustered regularly interspaced short palindromic repeats (CRISPR)-CRISPR-associated protein (Cas) has allowed breeders to induce genetic variability in a controlled and site-specific manner, and to improve traits with high efficiency. However, the presence of genetic linkages is a major obstacle to the transfer of desirable traits from wild species to their cultivated relatives. One way to address this issue is to create mutants with deficiencies in the meiotic recombination machinery, thereby enhancing global crossover frequencies between homologous parental chromosomes. Although this seemed to be a promising approach at first, thus far, no crossover frequencies could be enhanced in recombination-cold regions of the genome. Additionally, this approach can lead to unintended genomic instabilities due to DNA repair defects. Therefore, efforts have been undertaken to obtain predefined crossovers between homologues by inducing site-specific double-strand breaks (DSBs) in meiotic, as well as in somatic plant cells using CRISPR-Cas tools. However, this strategy has not been able to produce a substantial number of heritable homologous recombination-based crossovers. Most recently, heritable chromosomal rearrangements, such as inversions and translocations, have been obtained in a controlled way using CRISPR-Cas in plants. This approach unlocks a completely new way of manipulating genetic linkages, one in which the DSBs are induced in somatic cells, enabling the formation of chromosomal rearrangements in the megabase range, by DSB repair via non-homologous end-joining. This technology might also enable the restructuring of genomes more globally, resulting in not only the obtainment of synthetic plant chromosome, but also of novel plant species.

A New Catalog of Structural Variants in 1,301 A. thaliana Lines from Africa, Eurasia, and North America Reveals a Signature of Balancing Selection at Defense Response Genes

Genomic variation in the model plant Arabidopsis thaliana has been extensively used to understand evolutionary processes in natural populations, mainly focusing on single-nucleotide polymorphisms. Conversely, structural variation has been largely ignored in spite of its potential to dramatically affect phenotype. Here, we identify 155,440 indels and structural variants ranging in size from 1 bp to 10 kb, including presence/absence variants (PAVs), inversions, and tandem duplications in 1,301 A. thaliana natural accessions from Morocco, Madeira, Europe, Asia, and North America. We show evidence for strong purifying selection on PAVs in genes, in particular for housekeeping genes and homeobox genes, and we find that PAVs are concentrated in defense-related genes (R-genes, secondary metabolites) and F-box genes. This implies the presence of a "core" genome underlying basic cellular processes and a "flexible" genome that includes genes that may be important in spatially or temporally varying selection. Further, we find an excess of intermediate frequency PAVs in defense response genes in nearly all populations studied, consistent with a history of balancing selection on this class of genes. Finally, we find that PAVs in genes involved in the cold requirement for flowering (vernalization) and drought response are strongly associated with temperature at the sites of origin.

Detection of ploidy and chromosomal aberrations in commercial oil palm using high-throughput SNP markers

Karyotyping using high-density genome-wide SNP markers identified various chromosomal aberrations in oil palm (Elaeis guineensis Jacq.) with supporting evidence from the 2C DNA content measurements (determined using FCM) and chromosome counts. Oil palm produces a quarter of the world's total vegetable oil. In line with its global importance, an initiative to sequence the oil palm genome was carried out successfully, producing huge amounts of sequence information, allowing SNP discovery. High-capacity SNP genotyping platforms have been widely used for marker-trait association studies in oil palm. Besides genotyping, a SNP array is also an attractive tool for understanding aberrations in chromosome inheritance. Exploiting this, the present study utilized chromosome-wide SNP allelic distributions to determine the ploidy composition of over 1,000 oil palms from a commercial F₁ family, including 197 derived from twin-embryo seeds. Our method consisted of an inspection of the allelic intensity ratio using SNP markers. For palms with a shifted or abnormal distribution ratio, the SNP allelic frequencies were plotted along the pseudo-chromosomes. This method proved to be efficient in identifying whole genome duplication (triploids) and aneuploidy. We also detected several loss of heterozygosity regions which may indicate small chromosomal deletions and/or inheritance of identical by descent regions from both parents. The SNP analysis was validated by flow cytometry and chromosome counts. The triploids were all derived from twin-embryo seeds. This is the first report on the efficiency and reliability of SNP array data for karyotyping oil palm chromosomes, as an alternative to the conventional cytogenetic technique. Information on the ploidy composition and chromosomal structural variation can help to better understand the genetic makeup of samples and lead to a more robust interpretation of the genomic data in marker-trait association analyses.

CRISPR/Cas-mediated chromosome engineering: opening up a new avenue for plant breeding

The advent of powerful site-specific nucleases, particularly the clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein (Cas) system, which enables precise genome manipulation, has revolutionized plant breeding. Until recently, the main focus of researchers has been to simply knock-in or knock-out single genes, or to induce single base changes, but constant improvements of this technology have enabled more ambitious applications that aim to improve plant productivity or other desirable traits. One long-standing aim has been the induction of targeted chromosomal rearrangements (crossovers, inversions, or translocations). The feasibility of this technique has the potential to transform plant breeding, because natural rearrangements, like inversions, for example, typically present obstacles to the breeding process. In this way, genetic linkages between traits could be altered to combine or separate favorable and deleterious genes, respectively. In this review, we discuss recent breakthroughs in the field of chromosome engineering in plants and their potential applications in the field of plant breeding. In the future, these approaches might be applicable in shaping plant chromosomes in a directed manner, based on plant breeding needs.

Mutagenic Effect of Proton Beams Characterized by Phenotypic Analysis and Whole Genome Sequencing in Arabidopsis

Protons may have contributed to the evolution of plants as a major component of cosmic-rays and also have been used for mutagenesis in plants. Although the mutagenic effect of protons has been well-characterized in animals, no comprehensive phenotypic and genomic analyses has been reported in plants. Here, we investigated the phenotypes and whole genome sequences of Arabidopsis M₂ lines derived by irradiation with proton beams and gamma-rays, to determine unique characteristics of proton beams in mutagenesis. We found that mutation frequency was dependent on the irradiation doses of both proton beams and gamma-rays. On the basis of the relationship between survival and mutation rates, we hypothesized that there may be a mutation rate threshold for survived individuals after irradiation. There were no significant differences between the total mutation rates in groups derived using proton beam or gamma-ray irradiation at doses that had similar impacts on survival rate. However, proton beam irradiation resulted in a broader mutant phenotype spectrum than gamma-ray irradiation, and proton beams generated more DNA structural variations (SVs) than gamma-rays. The most frequent SV was inversion. Most of the inversion junctions contained sequences with microhomology and were associated with the deletion of only a few nucleotides, which implies that preferential use of microhomology in non-homologous end joining was likely to be responsible for the SVs. These results show that protons, as particles with low linear energy transfer (LET), have unique characteristics in mutagenesis that partially overlap with those of low-LET gamma-rays and high-LET heavy ions in different respects.

Ultra-High-Density QTL Marker Mapping for Seedling Photomorphogenesis Mediating Arabidopsis Establishment in Southern Patagonia

Arabidopsis thaliana shows a wide range of genetic and trait variation among wild accessions. Because of its unparalleled biological and genomic resources, Arabidopsis has a high potential for the identification of genes underlying ecologically important complex traits, thus providing new insights on genome evolution. Previous research suggested that distinct light responses were crucial for Arabidopsis establishment in a peculiar ecological niche of southern Patagonia. The aim of this study was to explore the genetic basis of contrasting light-associated physiological traits that may have mediated the rapid adaptation to this new environment. From a biparental cross between the photomorphogenic contrasting accessions Patagonia (Pat) and Columbia (Col-0), we generated a novel recombinant inbred line (RIL) population, which was entirely next-generation sequenced to achieve ultra-high-density saturating molecular markers resulting in supreme mapping sensitivity. We validated the quality of the RIL population by quantitative trait loci (QTL) mapping for seedling de-etiolation, finding seven QTLs for hypocotyl length in the dark and continuous blue light (Bc), continuous red light (Rc), and continuous far-red light (FRc). The most relevant QTLs, Rc1 and Bc1, were mapped close together to chromosome V; the former for Rc and Rc/dark, and the latter for Bc, FRc, and dark treatments. The additive effects of both QTLs were confirmed by independent heterogeneous inbred families (HIFs), and we explored TZP and ABA1 as potential candidate genes for Rc1 and Bc1QTLs, respectively. We conclude that the Pat × Col-0 RIL population is a valuable novel genetic resource to explore other adaptive traits in Arabidopsis.

Inside the supergene of the bird with four sexes

The white-throated sparrow (Zonotrichia albicollis) offers unique opportunities to understand the adaptive value of supergenes, particularly their role in alternative phenotypes. In this species, alternative plumage morphs segregate with a nonrecombining segment of chromosome 2, which has been called a 'supergene'. The species mates disassortatively with respect to the supergene; that is, each breeding pair consists of one individual with it and one without it. This species has therefore been called the "bird with four sexes". The supergene segregates with a behavioral phenotype; birds with it are more aggressive and less parental than birds without it. Here, we review our efforts to identify the genes inside the supergene that are responsible for the behavioral polymorphism. The gene ESR1, which encodes estrogen receptor α, differs between the morphs and predicts both territorial and parental behavior. Variation in the regulatory regions of ESR1 causes an imbalance in expression of the two alleles, and the degree to which this imbalance favors the supergene allele predicts territorial singing. In heterozygotes, knockdown of ESR1 causes a phenotypic switch, from more aggressive to less aggressive. We recently showed that another gene important for social behavior, vasoactive intestinal peptide (VIP), is differentially expressed between the morphs and predicts territorial singing. We hypothesize that ESR1 and VIP contribute to behavior in a coordinated way and could represent co-adapted alleles. Because the supergene contains more than 1000 individual genes, this species provides rich possibilities for discovering alleles that work together to mediate life-history trade-offs and maximize the fitness of alternative complex phenotypes.

Ribosomal DNA localization on Lathyrus species chromosomes by FISH

BACKGROUND

Fluorescence In Situ Hybridization (FISH) played an essential role to locate the ribosomal RNA genes on the chromosomes that offered a new tool to study the chromosome structure and evolution in plant. The 45S and 5S rRNA genes are independent and localized at one or more loci per the chromosome complement, their positions along chromosomes offer useful markers for chromosome discriminations. In the current study FISH has been performed to locate 45S and 5S rRNA genes on the chromosomes of nine Lathyrus species belong to five different sections, all have chromosome number 2n=14, Lathyrus gorgoni Parl, Lathyrus hirsutus L., Lathyrus amphicarpos L., Lathyrus odoratus L., Lathyrus sphaericus Retz, Lathyrus incospicuus L, Lathyrus paranensis Burkart, Lathyrus nissolia L., and Lathyrus articulates L.

RESULTS

The revealed loci of 45S and 5S rDNA by FISH on metaphase chromosomes of the examined species were as follow: all of the studied species have one 45S rDNA locus and one 5S rDNA locus except L. odoratus L., L. amphicarpos L. and L. sphaericus Retz L. have two loci of 5S rDNA. Three out of the nine examined species have the loci of 45S and 5S rRNA genes on the opposite arms of the same chromosome (L. nissolia L., L. amphicarpos L., and L. incospicuus L.), while L. hirsutus L. has both loci on the same chromosome arm. The other five species showed the loci of the two types of rDNA on different chromosomes.

CONCLUSION

The detected 5S and 45S rDNA loci in Lathyrus could be used as chromosomal markers to discriminate the chromosome pairs of the examined species. FISH could discriminate only one chromosome pair out of the seven pairs in three species, in L. hirsutus L., L. nissolia L. and L. incospicuus L., and two chromosome pairs in five species, in L. paranensis Burkart, L. odoratus L., L. amphicarpos L., L. gorgoni Parl. and L. articulatus L., while it could discriminate three chromosome pairs in L. sphaericus Retz. these results could contribute into the physical genome mapping of Lathyrus species and the evolution of rDNA patterns by FISH in the coming studies in future.

Changing local recombination patterns in Arabidopsis by CRISPR/Cas mediated chromosome engineering

Chromosomal inversions are recurrent rearrangements that occur between different plant isolates or cultivars. Such inversions may underlie reproductive isolation in evolution and represent a major obstacle for classical breeding as no crossovers can be observed between inverted sequences on homologous chromosomes. The heterochromatic knob (hk4S) on chromosome 4 is the most well-known inversion of Arabidopsis. If a knob carrying accession such as Col-0 is crossed with a knob-less accession such as Ler-1, crossovers cannot be recovered within the inverted region. Our work shows that by egg-cell specific expression of the Cas9 nuclease from Staphylococcus aureus, a targeted reversal of the 1.1 Mb long hk4S-inversion can be achieved. By crossing Col-0 harbouring the rearranged chromosome 4 with Ler-1, meiotic crossovers can be restored into a region with previously no detectable genetic exchange. The strategy of somatic chromosome engineering for breaking genetic linkage has huge potential for application in plant breeding.

Changing local recombination patterns in Arabidopsis by CRISPR/Cas mediated chromosome engineering

Chromosomal inversions are recurrent rearrangements that occur between different plant isolates or cultivars. Such inversions may underlie reproductive isolation in evolution and represent a major obstacle for classical breeding as no crossovers can be observed between inverted sequences on homologous chromosomes. The heterochromatic knob (hk4S) on chromosome 4 is the most well-known inversion of Arabidopsis. If a knob carrying accession such as Col-0 is crossed with a knob-less accession such as Ler-1, crossovers cannot be recovered within the inverted region. Our work shows that by egg-cell specific expression of the Cas9 nuclease from Staphylococcus aureus, a targeted reversal of the 1.1 Mb long hk4S-inversion can be achieved. By crossing Col-0 harbouring the rearranged chromosome 4 with Ler-1, meiotic crossovers can be restored into a region with previously no detectable genetic exchange. The strategy of somatic chromosome engineering for breaking genetic linkage has huge potential for application in plant breeding.

QTL for horticulturally important traits associated with pleiotropic andromonoecy and carpel number loci, and a paracentric inversion in cucumber

The legendary cucumber inbred line WI2757 possesses a rare combination of resistances against nine pathogens, which is an important germplasm for cucumber breeding. However, WI2757 flowers late and does not perform well under field conditions. The genetic basis for horticulturally important traits other than disease resistances in WI2757 is largely unknown. In this study, we conducted QTL mapping using F₂ and recombinant inbred line (RIL) populations from the WI2757 × True Lemon cross that were segregating for multiple traits. Phenotypic data were collected in replicated field trials across multiple years for seven traits including fruit carpel number (CN) and sex expression. A high-density SNP-based genetic map was developed with genotyping by sequencing of the RIL population, which revealed a region on chromosome 1 with strong recombination suppression. The reduced recombination in this region was due to a ~ 10-Mbp paracentric inversion in WI2757 that was confirmed with additional segregation and cytological (FISH) analyses. Thirty-six QTL were detected for flowering time, fruit length (FL), fruit diameter (FD), fruit shape (LD), fruit number (FN), CN, and powdery mildew resistance. Five moderate- or major-effect QTL for FL, FD, LD, and FN inside the inversion are likely the pleiotropic effects of the andromonoecy (m), or the cn locus. The major-effect flowering time QTL ft1.1 was also mapped inside the inversion, which seems to be different from the previously assigned delayed flowering in WI2757. Implications of these findings on the use of WI2757 in cucumber breeding are discussed.

Comparative genomics of four strains of the edible brown alga, Cladosiphon okamuranus

BACKGROUND

The brown alga, Cladosiphon okamuranus (Okinawa mozuku), is one of the most important edible seaweeds, and it is cultivated for market primarily in Okinawa, Japan. Four strains, denominated S, K, O, and C, with distinctively different morphologies, have been cultivated commercially since the early 2000s. We previously reported a draft genome of the S-strain. To facilitate studies of seaweed biology for future aquaculture, we here decoded and analyzed genomes of the other three strains (K, O, and C).

RESULTS

Here we improved the genome of the S-strain (ver. 2, 130 Mbp, 12,999 genes), and decoded the K-strain (135 Mbp, 12,511 genes), the O-strain (140 Mbp, 12,548 genes), and the C-strain (143 Mbp, 12,182 genes). Molecular phylogenies, using mitochondrial and nuclear genes, showed that the S-strain diverged first, followed by the K-strain, and most recently the C- and O-strains. Comparisons of genome architecture among the four strains document the frequent occurrence of inversions. In addition to gene acquisitions and losses, the S-, K-, O-, and C-strains possess 457, 344, 367, and 262 gene families unique to each strain, respectively. Comprehensive Blast searches showed that most genes have no sequence similarity to any entries in the non-redundant protein sequence database, although GO annotation suggested that they likely function in relation to molecular and biological processes and cellular components.

CONCLUSIONS

Our study compares the genomes of four strains of C. okamuranus and examines their phylogenetic relationships. Due to global environmental changes, including temperature increases, acidification, and pollution, brown algal aquaculture is facing critical challenges. Genomic and phylogenetic information reported by the present research provides useful tools for isolation of novel strains.

Chromosomal Evolution and Apomixis in the Cruciferous Tribe Boechereae

The mustard family (Brassicaceae) comprises several dozen monophyletic clades usually ranked as tribes. The tribe Boechereae plays a prominent role in plant research due to the incidence of apomixis and its close relationship to Arabidopsis. This tribe, largely confined to western North America, harbors nine genera and c. 130 species, with >90% of species belonging to the genus Boechera. Hundreds of apomictic diploid and triploid Boechera hybrids have spurred interest in this genus, but the remaining Boechereae genomes remain virtually unstudied. Here we report on comparative genome structure of six genera (Borodinia, Cusickiella, Phoenicaulis, Polyctenium, Nevada, and Sandbergia) and three Boechera species as revealed by comparative chromosome painting (CCP). All analyzed taxa shared the same seven-chromosome genome structure. Comparisons with the sister Halimolobeae tribe (n = 8) showed that the ancestral Boechereae genome (n = 7) was derived from an older n = 8 genome by descending dysploidy followed by the divergence of extant Boechereae taxa. As tribal divergence post-dated the origin of four tribe-specific chromosomes, it is proposed that these chromosomal rearrangements were a key evolutionary innovation underlaying the origin and diversification of the Boechereae in North America. Although most Boechereae genera exhibit genomic conservatism, intra-tribal cladogenesis has occasionally been accompanied by chromosomal rearrangements (particularly inversions). Recently, apomixis was reported in the Boechereae genera Borodinia and Phoenicaulis. Here, we report sexual reproduction in diploid Nevada, diploid Sandbergia, and tetraploid Cusickiella and aposporous apomixis in tetraploids of Polyctenium and Sandbergia. In sum, apomixis is now known to occur in five of the nine Boechereae genera.

Frequency, Origins, and Evolutionary Role of Chromosomal Inversions in Plants

Chromosomal inversions have the potential to play an important role in evolution by reducing recombination between favorable combinations of alleles. Until recently, however, most evidence for their likely importance derived from dipteran flies, whose giant larval salivary chromosomes aided early cytogenetic studies. The widespread application of new genomic technologies has revealed that inversions are ubiquitous across much of the plant and animal kingdoms. Here we review the rapidly accumulating literature on inversions in the plant kingdom and discuss what we have learned about their establishment and likely evolutionary role. We show that inversions are prevalent across a wide range of plant groups. We find that inversions are often associated with locally favored traits, as well as with traits that contribute to assortative mating, suggesting that they may be key to adaptation and speciation in the face of gene flow. We also discuss the role of inversions in sex chromosome formation, and explore possible parallels with inversion establishment on autosomes. The identification of inversion origins, as well as the causal variants within them, will advance our understanding of chromosomal evolution in plants.

Frequency, Origins, and Evolutionary Role of Chromosomal Inversions in Plants

Chromosomal inversions have the potential to play an important role in evolution by reducing recombination between favorable combinations of alleles. Until recently, however, most evidence for their likely importance derived from dipteran flies, whose giant larval salivary chromosomes aided early cytogenetic studies. The widespread application of new genomic technologies has revealed that inversions are ubiquitous across much of the plant and animal kingdoms. Here we review the rapidly accumulating literature on inversions in the plant kingdom and discuss what we have learned about their establishment and likely evolutionary role. We show that inversions are prevalent across a wide range of plant groups. We find that inversions are often associated with locally favored traits, as well as with traits that contribute to assortative mating, suggesting that they may be key to adaptation and speciation in the face of gene flow. We also discuss the role of inversions in sex chromosome formation, and explore possible parallels with inversion establishment on autosomes. The identification of inversion origins, as well as the causal variants within them, will advance our understanding of chromosomal evolution in plants.

Chromosomal features revealed by comparison of genetic maps of Glycine max and Glycine soja

Recombination is a crucial component of evolution and breeding. New combinations of variation on chromosomes are shaped by recombination. Recombination is also involved in chromosomal rearrangements. However, recombination rates vary tremendously among chromosome segments. Genome-wide genetic maps are one of the best tools to study variation of recombination. Here, we describe high density genetic maps of Glycine max and Glycine soja constructed from four segregating populations. The maps were used to identify chromosomal rearrangements and find the highly predictable pattern of cross-overs on the broad scale in soybean. Markers on these genetic maps were used to evaluate assembly quality of the current soybean reference genome sequence. We find a strong inversion candidate larger than 3 Mb based on patterns of cross-overs. We also identify quantitative trait loci (QTL) that control number of cross-overs. This study provides fundamental insights relevant to practical strategy for breeding programs and for pan-genome researches.

From gene editing to genome engineering: restructuring plant chromosomes via CRISPR/Cas

In the last years, tremendous progress has been achieved in the field of gene editing in plants. By the induction of single site-specific double-strand breaks (DSBs), the knockout of genes by non-homologous end joining has become routine in many plant species. Recently, the efficiency of inducing pre-planned mutations by homologous recombination has also been improved considerably. However, very little effort has been undertaken until now to achieve more complex changes in plant genomes by the simultaneous induction of several DSBs. Several reports have been published on the efficient induction of deletions. However, the induction of intrachromosomal inversions and interchromosomal recombination by the use of CRISPR/Cas has only recently been reported. In this review, we want to sum up these results and put them into context with regards to what is known about natural chromosome rearrangements in plants. Moreover, we review the recent progress in CRISPR/Cas-based mammalian chromosomal rearrangements, which might be inspiring for plant biologists. In the long run, the controlled restructuring of plant genomes should enable us to link or break linkage of traits at will, thus defining a new area of plant breeding.

Genome Size Variation and Comparative Genomics Reveal Intraspecific Diversity in Brassica rapa

Traditionally, reference genomes in crop species rely on the assembly of one accession, thus occulting most of intraspecific diversity. However, rearrangements, gene duplications, and transposable element content may have a large impact on the genomic structure, which could generate new phenotypic traits. Comparing two Brassica rapa genomes recently sequenced and assembled using long-read technology and optical mapping, we investigated structural variants and repetitive content between the two accessions and genome size variation among a core collection. We explored the structural consequences of the presence of large repeated sequences in B. rapa 'Z1' genome vs. the B. rapa 'Chiifu' genome, using comparative genomics and cytogenetic approaches. First, we showed that large genomic variants on chromosomes A05, A06, A09, and A10 are due to large insertions and inversions when comparing B. rapa 'Z1' and B. rapa 'Chiifu' at the origin of important length differences in some chromosomes. For instance, lengths of 'Z1' and 'Chiifu' A06 chromosomes were estimated in silico to be 55 and 29 Mb, respectively. To validate these observations, we compared using fluorescent in situ hybridization (FISH) the two A06 chromosomes present in an F1 hybrid produced by crossing these two varieties. We confirmed a length difference of 17.6% between the A06 chromosomes of 'Z1' compared to 'Chiifu.' Alternatively, using a copy number variation approach, we were able to quantify the presence of a higher number of rDNA and gypsy elements in 'Z1' genome compared to 'Chiifu' on different chromosomes including A06. Using flow cytometry, the total genome size of 12 Brassica accessions corresponding to a B. rapa available core collection was estimated and revealed a genome size variation of up to 16% between these accessions as well as some shared inversions. This study revealed the contribution of long-read sequencing of new accessions belonging to different cultigroups of B. rapa and highlighted the potential impact of differential insertion of repeat elements and inversions of large genomic regions in genome size intraspecific variability.

Recombination suppression in heterozygotes for a pericentric inversion induces the interchromosomal effect on crossovers in Arabidopsis

During meiosis, recombination ensures allelic exchanges through crossovers (COs) between the homologous chromosomes. Advances in our understanding of the rules of COs have come from studies of mutations including structural chromosomal rearrangements that, when heterozygous, are known to impair COs in various organisms. In this work, we investigate the effect of a large heterozygous pericentric inversion on male and female recombination in Arabidopsis. The inversion was discovered in the Atmcc1 mutant background and was characterized through genetic and next-generation sequencing analysis. Reciprocal backcross populations, each consisting of over 400 individuals, obtained from the mutant and the wild type, both crossed with Landsberg erecta, were analyzed genome-wide by 143 single-nucleotide polymorphisms. The negative impact of inversion became evident in terms of CO loss in the rearranged chromosome in both male and female meiosis. No single-CO event was detected within the inversion, consistent with a post-meiotic selection operating against unbalanced gametes. Cytological analysis of chiasmata in F₁ plants confirmed that COs were reduced in male meiosis in the chromosome with inversion. Crossover suppression on the rearranged chromosome is associated with a significant increase of COs in the other chromosomes, thereby maintaining unchanged the number of COs per cell. The CO pattern observed in our study is consistent with the interchromosomal (IC) effect as first described in Drosophila. In contrast to male meiosis, in female meiosis no IC effect is visible. This may be related to the greater strength of interference that constrains the CO number in excess of the minimum value imposed by CO assurance in Arabidopsis female meiosis.

An Ultra High-Density Arabidopsis thaliana Crossover Map That Refines the Influences of Structural Variation and Epigenetic Features

Many environmental, genetic, and epigenetic factors are known to affect the frequency and positioning of meiotic crossovers (COs). Suppression of COs by large, cytologically visible inversions and translocations has long been recognized, but relatively little is known about how smaller structural variants (SVs) affect COs. To examine fine-scale determinants of the CO landscape, including SVs, we used a rapid, cost-effective method for high-throughput sequencing to generate a precise map of >17,000 COs between the Col-0 and Ler-0 accessions of Arabidopsis thaliana COs were generally suppressed in regions with SVs, but this effect did not depend on the size of the variant region, and was only marginally affected by the variant type. CO suppression did not extend far beyond the SV borders and CO rates were slightly elevated in the flanking regions. Disease resistance gene clusters, which often exist as SVs, exhibited high CO rates at some loci, but there was a tendency toward depressed CO rates at loci where large structural differences exist between the two parents. Our high-density map also revealed in fine detail how CO positioning relates to genetic (DNA motifs) and epigenetic (chromatin structure) features of the genome. We conclude that suppression of COs occurs over a narrow region spanning large- and small-scale SVs, representing an influence on the CO landscape in addition to sequence and epigenetic variation along chromosomes.

Population Genetics of the Highly Polymorphic RPP8 Gene Family

Plant nucleotide-binding domain and leucine-rich repeat containing (NLR) genes provide some of the most extreme examples of polymorphism in eukaryotic genomes, rivalling even the vertebrate major histocompatibility complex. Surprisingly, this is also true in Arabidopsis thaliana, a predominantly selfing species with low heterozygosity. Here, we investigate how gene duplication and intergenic exchange contribute to this extraordinary variation. RPP8 is a three-locus system that is configured chromosomally as either a direct-repeat tandem duplication or as a single copy locus, plus a locus 2 Mb distant. We sequenced 48 RPP8 alleles from 37 accessions of A. thaliana and 12 RPP8 alleles from Arabidopsis lyrata to investigate the patterns of interlocus shared variation. The tandem duplicates display fixed differences and share less variation with each other than either shares with the distant paralog. A high level of shared polymorphism among alleles at one of the tandem duplicates, the single-copy locus and the distal locus, must involve both classical crossing over and intergenic gene conversion. Despite these polymorphism-enhancing mechanisms, the observed nucleotide diversity could not be replicated under neutral forward-in-time simulations. Only by adding balancing selection to the simulations do they approach the level of polymorphism observed at RPP8. In this NLR gene triad, genetic architecture, gene function and selection all combine to generate diversity.

Genome invasion by a hypomethylated satellite repeat in Australian crucifer Ballantinia antipoda

Repetitive sequences are ubiquitous components of all eukaryotic genomes. They contribute to genome evolution and the regulation of gene transcription. However, the uncontrolled activity of repetitive sequences can negatively affect genome functions and stability. Therefore, repetitive DNAs are embedded in a highly repressive heterochromatic environment in plant cell nuclei. Here, we analyzed the sequence, composition and the epigenetic makeup of peculiar non-pericentromeric heterochromatic segments in the genome of the Australian crucifer Ballantinia antipoda. By the combination of high throughput sequencing, graph-based clustering and cytogenetics, we found that the heterochromatic segments consist of a mixture of unique sequences and an A-T-rich 174 bp satellite repeat (BaSAT1). BaSAT1 occupies about 10% of the B. antipoda nuclear genome in >250 000 copies. Unlike many other highly repetitive sequences, BaSAT1 repeats are hypomethylated; this contrasts with the normal patterns of DNA methylation in the B. antipoda genome. Detailed analysis of several copies revealed that these non-methylated BaSAT1 repeats were also devoid of heterochromatic histone H3K9me2 methylation. However, the factors decisive for the methylation status of BaSAT1 repeats remain currently unknown. In summary, we show that even highly repetitive sequences can exist as hypomethylated in the plant nuclear genome.

Comparative Tyramide-FISH mapping of the genes controlling flavor and bulb color in Allium species revealed an altered gene order

Evolutionarily related species often share a common order of genes along homeologous chromosomes. Here we report the collinearity disruption of genes located on homeologous chromosome 4 in Allium species. Ultra-sensitive fluorescence in situ hybridization with tyramide signal amplification (tyr-FISH) allowed the visualization of the alliinase multigene family, chalcon synthase gene and EST markers on Allium cepa and Allium fistulosum chromosomes. In A. cepa, bulb alliinase, root alliinase (ALL1) and chalcon synthase (CHS-B) genes were located in the long arm but EST markers (API18 and ACM082) were located in the short arm. In A. fistulosum, all the visualized genes and markers were located in the short arm. Moreover, root alliinase genes (ALL1 and AOB249) showed contrast patterns in number of loci. We suppose that the altered order of the genes/markers is the result of a large pericentric inversion. To get insight into the evolution of the chromosome rearrangement, we mapped the bulb alliinase gene in phylogenetically close and distant species. In the taxonomic clade including A. fistulosum, A. altaicum, A. oschaninii and A. pskemense and in phylogenetically distant species A. roylei and A. nutans, the bulb alliinase gene was located on the short arm of chromosome 4 while, in A. cepa and A. schoenoprasum, the bulb alliinase gene was located on the long arm of chromosome 4. These results have encouraging implications for the further tracing of inverted regions in meiosis of interspecific hybrids and studding chromosome evolution. Also, this finding may have a practical benefit as closely related species are actively used for improving onion crop stock.

Natural Variation in TBP-ASSOCIATED FACTOR 4b Controls Meiotic Crossover and Germline Transcription in Arabidopsis

Meiotic crossover frequency varies within genomes, which influences genetic diversity and adaptation. In turn, genetic variation within populations can act to modify crossover frequency in cis and trans. To identify genetic variation that controls meiotic crossover frequency, we screened Arabidopsis accessions using fluorescent recombination reporters. We mapped a genetic modifier of crossover frequency in Col × Bur populations of Arabidopsis to a premature stop codon within TBP-ASSOCIATED FACTOR 4b (TAF4b), which encodes a subunit of the RNA polymerase II general transcription factor TFIID. The Arabidopsis taf4b mutation is a rare variant found in the British Isles, originating in South-West Ireland. Using genetics, genomics, and immunocytology, we demonstrate a genome-wide decrease in taf4b crossovers, with strongest reduction in the sub-telomeric regions. Using RNA sequencing (RNA-seq) from purified meiocytes, we show that TAF4b expression is meiocyte enriched, whereas its paralog TAF4 is broadly expressed. Consistent with the role of TFIID in promoting gene expression, RNA-seq of wild-type and taf4b meiocytes identified widespread transcriptional changes, including in genes that regulate the meiotic cell cycle and recombination. Therefore, TAF4b duplication is associated with acquisition of meiocyte-specific expression and promotion of germline transcription, which act directly or indirectly to elevate crossovers. This identifies a novel mode of meiotic recombination control via a general transcription factor.

Extensive intraspecific gene order and gene structural variations in upland cotton cultivars

Multiple cotton genomes (diploid and tetraploid) have been assembled. However, genomic variations between cultivars of allotetraploid upland cotton (Gossypium hirsutum L.), the most widely planted cotton species in the world, remain unexplored. Here, we use single-molecule long read and Hi-C sequencing technologies to assemble genomes of the two upland cotton cultivars TM-1 and zhongmiansuo24 (ZM24). Comparisons among TM-1 and ZM24 assemblies and the genomes of the diploid ancestors reveal a large amount of genetic variations. Among them, the top three longest structural variations are located on chromosome A08 of the tetraploid upland cotton, which account for ~30% total length of this chromosome. Haplotype analyses of the mapping population derived from these two cultivars and the germplasm panel show suppressed recombination rates in this region. This study provides additional genomic resources for the community, and the identified genetic variations, especially the reduced meiotic recombination on chromosome A08, will help future breeding.

Efficient induction of heritable inversions in plant genomes using the CRISPR/Cas system

During the evolution of plant genomes, sequence inversions occurred repeatedly making the respective regions inaccessible for meiotic recombination and thus for breeding. Therefore, it is important to develop technologies that allow the induction of inversions within chromosomes in a directed and efficient manner. Using the Cas9 nuclease from Staphylococcus aureus (SaCas9), we were able to obtain scarless heritable inversions with high efficiency in the model plant Arabidopsis thaliana. Via deep sequencing, we defined the patterns of junction formation in wild-type and in the non-homologous end-joining (NHEJ) mutant ku70-1. Surprisingly, in plants deficient of KU70, inversion induction is enhanced, indicating that KU70 is required for tethering the local broken ends together during repair. However, in contrast to wild-type, most junctions are formed by microhomology-mediated NHEJ and thus are imperfect with mainly deletions, making this approach unsuitable for practical applications. Using egg-cell-specific expression of Cas9, we were able to induce heritable inversions at different genomic loci and at intervals between 3 and 18 kb, in the percentage range, in the T1 generation. By screening individual lines, inversion frequencies of up to the 10% range were found in T2. Most of these inversions had scarless junctions and were without any sequence change within the inverted region, making the technology attractive for use in crop plants. Applying our approach, it should be possible to reverse natural inversions and induce artificial ones to break or fix linkages between traits at will.

On the postglacial spread of human commensal Arabidopsis thaliana: journey to the East

With more sequenced genomes, our understanding of the demographic history of Arabidopsis thaliana is rapidly expanding. However, no-one has yet compiled previous data to investigate patterns of genetic variation across Eurasia. While sub-Saharan accessions have been reported to be the most divergent group, in the nuclear genome we found accessions from Yunnan, China to be genetically closest to the sub-Saharan group. In chloroplast, several deeply diverged haplogroups exist only in Eurasia, and African populations have lower variation in many haplogroups that they share with the Eurasian populations. These patterns cannot be easily explained by a single out-of-Africa event suggested previously. For more recent demographic history, we dated the nonrelict expansion to 10 ka. In the Chinese Yangtze nonrelicts, we found clear traces of gene flow with local relicts, and genes under strong selection were enriched for traces of relict introgression, especially those related to biotic and immune responses. The results suggest the ability of nonrelicts to obtain locally adaptive alleles through admixture with relicts is important for the expansion across environmental gradients of Eurasia. Our re-analyses provide another model for the early history as well as elucidating factors contributing to the recent demographic turnover event of this species.

GOOGA: A platform to synthesize mapping experiments and identify genomic structural diversity

Understanding genomic structural variation such as inversions and translocations is a key challenge in evolutionary genetics. We develop a novel statistical approach to comparative genetic mapping to detect large-scale structural mutations from low-level sequencing data. The procedure, called Genome Order Optimization by Genetic Algorithm (GOOGA), couples a Hidden Markov Model with a Genetic Algorithm to analyze data from genetic mapping populations. We demonstrate the method using both simulated data (calibrated from experiments on Drosophila melanogaster) and real data from five distinct crosses within the flowering plant genus Mimulus. Application of GOOGA to the Mimulus data corrects numerous errors (misplaced sequences) in the M. guttatus reference genome and confirms or detects eight large inversions polymorphic within the species complex. Finally, we show how this method can be applied in genomic scans to improve the accuracy and resolution of Quantitative Trait Locus (QTL) mapping.

Evolutionary and ecological functional genomics, from lab to the wild

Plant phenotypes are the result of both genetic and environmental forces that act to modulate trait expression. Over the last few years, numerous approaches in functional genomics and systems biology have led to a greater understanding of plant phenotypic variation and plant responses to the environment. These approaches, and the questions that they can address, have been loosely termed evolutionary and ecological functional genomics (EEFG), and have been providing key insights on how plants adapt and evolve. In particular, by bringing these studies from the laboratory to the field, EEFG studies allow us to gain greater knowledge of how plants function in their natural contexts.