Epigenetic and transcriptional consequences in the endosperm of chemically induced transposon mobilization in Arabidopsis.

Genomic imprinting, an epigenetic phenomenon leading to parent-of-origin-specific gene expression, has independently evolved in the endosperm of flowering plants and the placenta of mammals-tissues crucial for nurturing embryos. While transposable elements (TEs) frequently colocalize with imprinted genes and are implicated in imprinting establishment, direct investigations of the impact of de novo TE transposition on genomic imprinting remain scarce. In this study, we explored the effects of chemically induced transposition of the Copia element ONSEN on genomic imprinting in Arabidopsis thaliana. Through the combination of chemical TE mobilization and doubled haploid induction, we generated a line with 40 new ONSEN copies. Our findings reveal a preferential targeting of maternally expressed genes (MEGs) for transposition, aligning with the colocalization of H2A.Z and H3K27me3 in MEGs-both previously identified as promoters of ONSEN insertions. Additionally, we demonstrate that chemically-induced DNA hypomethylation induces global transcriptional deregulation in the endosperm, leading to the breakdown of MEG imprinting. This study provides insights into the consequences of chemically induced TE remobilization in the endosperm, revealing that chemically-induced epigenome changes can have long-term consequences on imprinted gene expression.

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.

Meiotic recombination dynamics in plants with repeat-based holocentromeres shed light on the primary drivers of crossover patterning.

Centromeres strongly affect (epi)genomic architecture and meiotic recombination dynamics, influencing the overall distribution and frequency of crossovers. Here we show how recombination is regulated and distributed in the holocentric plant Rhynchospora breviuscula, a species with diffused centromeres. Combining immunocytochemistry, chromatin analysis and high-throughput single-pollen sequencing, we discovered that crossover frequency is distally biased, in sharp contrast to the diffused distribution of hundreds of centromeric units and (epi)genomic features. Remarkably, we found that crossovers were abolished inside centromeric units but not in their proximity, indicating the absence of a canonical centromere effect. We further propose that telomere-led synapsis of homologues is the feature that best explains the observed recombination landscape. Our results hint at the primary influence of mechanistic features of meiotic pairing and synapsis rather than (epi)genomic features and centromere organization in determining the distally biased crossover distribution in R. breviuscula, whereas centromeres and (epi)genetic properties only affect crossover positioning locally.

Reliable genotyping of recombinant genomes using a robust hidden Markov model.

Meiotic recombination is an essential mechanism during sexual reproduction and includes the exchange of chromosome segments between homologous chromosomes. New allelic combinations are transmitted to the new generation, introducing novel genetic variation in the offspring genomes. With the improvement of high-throughput whole-genome sequencing technologies, large numbers of recombinant individuals can now be sequenced with low sequencing depth at low costs, necessitating computational methods for reconstructing their haplotypes. The main challenge is the uncertainty in haplotype calling that arises from the low information content of a single genomic position. Straightforward sliding window-based approaches are difficult to tune and fail to place recombination breakpoints precisely. Hidden Markov model (HMM)-based approaches, on the other hand, tend to over-segment the genome. Here, we present RTIGER, an HMM-based model that exploits in a mathematically precise way the fact that true chromosome segments typically have a certain minimum length. We further separate the task of identifying the correct haplotype sequence from the accurate placement of haplotype borders, thereby maximizing the accuracy of border positions. By comparing segmentations based on simulated data with known underlying haplotypes, we highlight the reasons for RTIGER outperforming traditional segmentation approaches. We then analyze the meiotic recombination pattern of segregants of 2 Arabidopsis (Arabidopsis thaliana) accessions and a previously described hyper-recombining mutant. RTIGER is available as an R package with an efficient Julia implementation of the core algorithm.

Gamete Binning to Achieve Haplotype-Resolved Genome Assembly.

Haplotype-resolved genome assemblies remain a challenge in practice. Here, we provide a step-by-step guide on gamete binning, a method to generate haplotype-resolved genome assemblies for diploid species. The protocol starts by phasing heterozygous variants to individual haplotypes of specific chromosomes using the genome information of individual haploid gametes of the focal individual. Using phased variants, the whole-genome sequencing reads from the diploid genome can be genotyped and assigned into groups, which represent the individual haplotypes of each of the chromosomes. Finally, haplotype-specific chromosomes can be assembled independently using standard assembly tools. First applications of gamete binning revealed a haplotyping accuracy over 99%, which outperformed sequence-only or Hi-C-based haplotype-resolved genome assemblies.Availability: github.com/schneebergerlab/GameteBinning_prac .

Rapid large-scale genomic introgression in Arabidopsis suecica via an autoallohexaploid bridge.

Gene flow between species in the genus Arabidopsis occurs in significant amounts, but how exactly gene flow is achieved is not well understood. Polyploidization may be one avenue to explain gene flow between species. One problem, however, with polyploidization as a satisfying explanation is the occurrence of lethal genomic instabilities in neopolyploids as a result of genomic exchange, erratic meiotic behavior, and genomic shock. We have created an autoallohexaploid by pollinating naturally co-occurring diploid Arabidopsis thaliana with allotetraploid Arabidopsis suecica (an allotetraploid composed of A. thaliana and Arabidopsis arenosa). Its triploid offspring underwent spontaneous genome duplication and was used to generate a multigenerational pedigree. Using genome resequencing, we show that 2 major mechanisms promote stable genomic exchange in this population. Legitimate meiotic recombination and chromosome segregation between the autopolyploid chromosomes of the 2 A. thaliana genomes occur without any obvious bias for the parental origin and combine the A. thaliana haplotypes from the A. thaliana parent with the A. thaliana haplotypes from A. suecica similar to purely autopolyploid plants. In addition, we repeatedly observed that occasional exchanges between regions of the homoeologous chromosomes are tolerated. The combination of these mechanisms may result in gene flow leading to stable introgression in natural populations. Unlike the previously reported resynthesized neoallotetraploid A. suecica, this population of autoallohexaploids contains mostly vigorous, and genetically, cytotypically, and phenotypically variable individuals. We propose that naturally formed autoallohexaploid populations might serve as an intermediate bridge between diploid and polyploid species, which can facilitate gene flow rapidly and efficiently.

Barley Ror1 encodes a class XI myosin required for mlo-based broad-spectrum resistance to the fungal powdery mildew pathogen.

Loss-of-function alleles of plant MLO genes confer broad-spectrum resistance to powdery mildews in many eudicot and monocot species. Although barley (Hordeum vulgare) mlo mutants have been used in agriculture for more than 40 years, understanding of the molecular principles underlying this type of disease resistance remains fragmentary. Forward genetic screens in barley have revealed mutations in two Required for mlo resistance (Ror) genes that partially impair immunity conferred by mlo mutants. While Ror2 encodes a soluble N-ethylmaleimide-sensitive factor-attached protein receptor (SNARE), the identity of Ror1, located at the pericentromeric region of barley chromosome 1H, remained elusive. We report the identification of Ror1 based on combined barley genomic sequence information and transcriptomic data from ror1 mutant plants. Ror1 encodes the barley class XI myosin Myo11A (HORVU.MOREX.r3.1HG0046420). Single amino acid substitutions of this myosin, deduced from non-functional ror1 mutant alleles, map to the nucleotide-binding region and the interface between the relay-helix and the converter domain of the motor protein. Ror1 myosin accumulates transiently in the course of powdery mildew infection. Functional fluorophore-labeled Ror1 variants associate with mobile intracellular compartments that partially colocalize with peroxisomes. Single-cell expression of the Ror1 tail region causes a dominant-negative effect that phenocopies ror1 loss-of-function mutants. We define a myosin motor for the establishment of mlo-mediated resistance, suggesting that motor protein-driven intracellular transport processes are critical for extracellular immunity, possibly through the targeted transfer of antifungal and/or cell wall cargoes to pathogen contact sites.

The megabase-scale crossover landscape is largely independent of sequence divergence.

Meiotic recombination frequency varies along chromosomes and strongly correlates with sequence divergence. However, the causal relationship between recombination landscapes and polymorphisms is unclear. Here, we characterize the genome-wide recombination landscape in the quasi-absence of polymorphisms, using Arabidopsis thaliana homozygous inbred lines in which a few hundred genetic markers were introduced through mutagenesis. We find that megabase-scale recombination landscapes in inbred lines are strikingly similar to the recombination landscapes in hybrids, with the notable exception of heterozygous large rearrangements where recombination is prevented locally. In addition, the megabase-scale recombination landscape can be largely explained by chromatin features. Our results show that polymorphisms are not a major determinant of the shape of the megabase-scale recombination landscape but rather favour alternative models in which recombination and chromatin shape sequence divergence across the genome.

FLOWERING REPRESSOR AAA+ ATPase 1 is a novel regulator of perennial flowering in Arabis alpina.

Arabis alpina is a polycarpic perennial, in which PERPETUAL FLOWERING1 (PEP1) regulates flowering and perennial traits in a vernalization-dependent manner. Mutagenesis screens of the pep1 mutant established the role of other flowering time regulators in PEP1-parallel pathways. Here we characterized three allelic enhancers of pep1 (eop002, 085 and 091) which flower early. We mapped the causal mutations and complemented mutants with the identified gene. Using quantitative reverse transcriptase PCR and reporter lines, we determined the protein spatiotemporal expression patterns and localization within the cell. We also characterized its role in Arabidopsis thaliana using CRISPR and in A. alpina by introgressing mutant alleles into a wild-type background. These mutants carried lesions in an AAA+ ATPase of unknown function, FLOWERING REPRESSOR AAA+ ATPase 1 (AaFRAT1). AaFRAT1 was detected in the vasculature of young leaf primordia and the rib zone of flowering shoot apical meristems. At the subcellular level, AaFRAT1 was localized at the interphase between the endoplasmic reticulum and peroxisomes. Introgression lines carrying Aafrat1 alleles required less vernalization to flower and reduced number of vegetative axillary branches. By contrast, A. thaliana CRISPR lines showed weak flowering phenotypes. AaFRAT1 contributes to flowering time regulation and the perennial growth habit of A. alpina.

plotsr: visualizing structural similarities and rearrangements between multiple genomes.

SUMMARY: Third-generation genome sequencing technologies have led to a sharp increase in the number of high-quality genome assemblies. This allows the comparison of multiple assembled genomes of individual species and demands new tools for visualizing their structural properties. Here, we present plotsr, an efficient tool to visualize structural similarities and rearrangements between genomes. It can be used to compare genomes on chromosome level or to zoom in on any selected region. In addition, plotsr can augment the visualization with regional identifiers (e.g. genes or genomic markers) or histogram tracks for continuous features (e.g. GC content or polymorphism density). AVAILABILITY AND IMPLEMENTATION: plotsr is implemented as a python package and uses the standard matplotlib library for plotting. It is freely available under the MIT license at GitHub (https://github.com/schneebergerlab/plotsr) and bioconda (https://anaconda.org/bioconda/plotsr). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Chromosome-scale and haplotype-resolved genome assembly of a tetraploid potato cultivar.

Potato is the most widely produced tuber crop worldwide. However, reconstructing the four haplotypes of its autotetraploid genome remained an unsolved challenge. Here, we report the 3.1 Gb haplotype-resolved (at 99.6% precision), chromosome-scale assembly of the potato cultivar 'Otava' based on high-quality long reads, single-cell sequencing of 717 pollen genomes and Hi-C data. Unexpectedly, ~50% of the genome was identical-by-descent due to recent inbreeding, which was contrasted by highly abundant structural rearrangements involving ~20% of the genome. Among 38,214 genes, only 54% were present in all four haplotypes with an average of 3.2 copies per gene. Taking the leaf transcriptome as an example, 11% of the genes were differently expressed in at least one haplotype, where 25% of them were likely regulated through allele-specific DNA methylation. Our work sheds light on the recent breeding history of potato, the functional organization of its tetraploid genome and has the potential to strengthen the future of genomics-assisted breeding.

Genetic and Molecular Analysis of Root Hair Development in Arabis alpina.

Root hair formation in Arabidopsis thaliana is a well-established model system for epidermal patterning and morphogenesis in plants. Over the last decades, many underlying regulatory genes and well-established networks have been identified by thorough genetic and molecular analysis. In this study, we used a forward genetic approach to identify genes involved in root hair development in Arabis alpina, a related crucifer species that diverged from A. thaliana approximately 26-40 million years ago. We found all root hair mutant classes known in A. thaliana and identified orthologous regulatory genes by whole-genome or candidate gene sequencing. Our findings indicate that the gene-phenotype relationships regulating root hair development are largely conserved between A. thaliana and A. alpina. Concordantly, a detailed analysis of one mutant with multiple hairs originating from one cell suggested that a mutation in the SUPERCENTIPEDE1 (SCN1) gene is causal for the phenotype and that AaSCN1 is fully functional in A. thaliana. Interestingly, we also found differences in the regulation of root hair differentiation and morphogenesis between the species, and a subset of root hair mutants could not be explained by mutations in orthologs of known genes from A. thaliana. This analysis provides insight into the conservation and divergence of root hair regulation in the Brassicaceae.

The genetic and epigenetic landscape of the Arabidopsis centromeres.

Centromeres attach chromosomes to spindle microtubules during cell division and, despite this conserved role, show paradoxically rapid evolution and are typified by complex repeats. We used long-read sequencing to generate the Col-CEN Arabidopsis thaliana genome assembly that resolves all five centromeres. The centromeres consist of megabase-scale tandemly repeated satellite arrays, which support CENTROMERE SPECIFIC HISTONE H3 (CENH3) occupancy and are densely DNA methylated, with satellite variants private to each chromosome. CENH3 preferentially occupies satellites that show the least amount of divergence and occur in higher-order repeats. The centromeres are invaded by ATHILA retrotransposons, which disrupt genetic and epigenetic organization. Centromeric crossover recombination is suppressed, yet low levels of meiotic DNA double-strand breaks occur that are regulated by DNA methylation. We propose that Arabidopsis centromeres are evolving through cycles of satellite homogenization and retrotransposon-driven diversification.

Transposition and duplication of MADS-domain transcription factor genes in annual and perennial Arabis species modulates flowering.

The timing of reproduction is an adaptive trait in many organisms. In plants, the timing, duration, and intensity of flowering differ between annual and perennial species. To identify interspecies variation in these traits, we studied introgression lines derived from hybridization of annual and perennial species, Arabis montbretiana and Arabis alpina, respectively. Recombination mapping identified two tandem A. montbretiana genes encoding MADS-domain transcription factors that confer extreme late flowering on A. alpina These genes are related to the MADS AFFECTING FLOWERING (MAF) cluster of floral repressors of other Brassicaceae species and were named A. montbretiana (Am) MAF-RELATED (MAR) genes. AmMAR1 but not AmMAR2 prevented floral induction at the shoot apex of A. alpina, strongly enhancing the effect of the MAF cluster, and MAR1 is absent from the genomes of all A. alpina accessions analyzed. Exposure of plants to cold (vernalization) represses AmMAR1 transcription and overcomes its inhibition of flowering. Assembly of the tandem arrays of MAR and MAF genes of six A. alpina accessions and three related species using PacBio long-sequence reads demonstrated that the MARs arose within the Arabis genus by interchromosomal transposition of a MAF1-like gene followed by tandem duplication. Time-resolved comparative RNA-sequencing (RNA-seq) suggested that AmMAR1 may be retained in A. montbretiana to enhance the effect of the AmMAF cluster and extend the duration of vernalization required for flowering. Our results demonstrate that MAF genes transposed independently in different Brassicaceae lineages and suggest that they were retained to modulate adaptive flowering responses that differ even among closely related species.

Using the longest run subsequence problem within homology-based scaffolding.

Genome assembly is one of the most important problems in computational genomics. Here, we suggest addressing an issue that arises in homology-based scaffolding, that is, when linking and ordering contigs to obtain larger pseudo-chromosomes by means of a second incomplete assembly of a related species. The idea is to use alignments of binned regions in one contig to find the most homologous contig in the other assembly. We show that ordering the contigs of the other assembly can be expressed by a new string problem, the longest run subsequence problem (LRS). We show that LRS is NP-hard and present reduction rules and two algorithmic approaches that, together, are able to solve large instances of LRS to provable optimality. All data used in the experiments as well as our source code are freely available. We demonstrate its usefulness within an existing larger scaffolding approach by solving realistic instances resulting from partial Arabidopsis thaliana assemblies in short computation time.

Anno genominis XX: 20 years of Arabidopsis genomics.

Twenty years ago, the Arabidopsis thaliana genome sequence was published. This was an important moment as it was the first sequenced plant genome and explicitly brought plant science into the genomics era. At the time, this was not only an outstanding technological achievement, but it was characterized by a superb global collaboration. The Arabidopsis genome was the seed for plant genomic research. Here, we review the development of numerous resources based on the genome that have enabled discoveries across plant species, which has enhanced our understanding of how plants function and interact with their environments.

Natural variation identifies a Pxy gene controlling vascular organisation and formation of nodules and lateral roots in Lotus japonicus.

Forward and reverse genetics using the model legumes Lotus japonicus and Medicago truncatula have been instrumental in identifying the essential genes governing legume-rhizobia symbiosis. However, little information is known about the effects of intraspecific variation on symbiotic signalling. Here, we use quantitative trait locus sequencing (QTL-seq) to investigate the genetic basis of the differentiated phenotypic responses shown by the Lotus accessions Gifu and MG20 to inoculation with the Mesorhizobium loti exoU mutant that produces truncated exopolysaccharides. We identified through genetic complementation the Pxy gene as a component of this differential exoU response. Lotus Pxy encodes a leucine-rich repeat receptor-like kinase similar to Arabidopsis thaliana PXY, which regulates stem vascular development. We show that Lotus pxy insertion mutants displayed defects in root and stem vascular organisation, as well as lateral root and nodule formation. Our work links Pxy to de novo organogenesis in the root, highlights the genetic overlap between regulation of lateral root and nodule formation, and demonstrates that natural variation in Pxy affects nodulation signalling.

A mis-regulated cyclic nucleotide-gated channel mediates cytosolic calcium elevation and activates immunity in Arabidopsis.

Calcium (Ca2+ ) is a second messenger for plant cell surface and intracellular receptors mediating pattern-triggered and effector-triggered immunity (respectively, PTI and ETI). Several CYCLIC NUCLEOTIDE-GATED CHANNELS (CNGCs) were shown to control transient cytosolic Ca2+ influx upon PTI activation. The contributions of specific CNGC members to PTI and ETI remain unclear. ENHANCED DISEASE SUSCEPTIBLITY1 (EDS1) regulates ETI signaling. In an Arabidopsis genetic screen for suppressors of eds1, we identify a recessive gain-of-function mutation in CNGC20, denoted cngc20-4, which partially restores disease resistance in eds1. cngc20-4 enhances PTI responses and ETI hypersensitive cell death. A cngc20-4 single mutant exhibits autoimmunity, which is dependent on genetically parallel EDS1 and salicylic acid (SA) pathways. CNGC20 self-associates, forms heteromeric complexes with CNGC19, and is phosphorylated and stabilized by BOTRYTIS INDUCED KINASE1 (BIK1). The cngc20-4 L371F exchange on a predicted transmembrane channel inward surface does not disrupt these interactions but leads to increased cytosolic Ca2+ accumulation, consistent with mis-regulation of CNGC20 Ca2+ -permeable channel activity. Our data show that ectopic Ca2+ influx caused by a mutant form of CNGC20 in cngc20-4 affects both PTI and ETI responses. We conclude that tight control of the CNGC20 Ca2+ ion channel is important for regulated immunity.

The Evolutionary Dynamics of Genetic Incompatibilities Introduced by Duplicated Genes in Arabidopsis thaliana.

Although gene duplications provide genetic backup and allow genomic changes under relaxed selection, they may potentially limit gene flow. When different copies of a duplicated gene are pseudofunctionalized in different genotypes, genetic incompatibilities can arise in their hybrid offspring. Although such cases have been reported after manual crosses, it remains unclear whether they occur in nature and how they affect natural populations. Here, we identified four duplicated-gene based incompatibilities including one previously not reported within an artificial Arabidopsis intercross population. Unexpectedly, however, for each of the genetic incompatibilities we also identified the incompatible alleles in natural populations based on the genomes of 1,135 Arabidopsis accessions published by the 1001 Genomes Project. Using the presence of incompatible allele combinations as phenotypes for GWAS, we mapped genomic regions that included additional gene copies which likely rescue the genetic incompatibility. Reconstructing the geographic origins and evolutionary trajectories of the individual alleles suggested that incompatible alleles frequently coexist, even in geographically closed regions, and that their effects can be overcome by additional gene copies collectively shaping the evolutionary dynamics of duplicated genes during population history.

Gamete binning: chromosome-level and haplotype-resolved genome assembly enabled by high-throughput single-cell sequencing of gamete genomes.

Generating chromosome-level, haplotype-resolved assemblies of heterozygous genomes remains challenging. To address this, we developed gamete binning, a method based on single-cell sequencing of haploid gametes enabling separation of the whole-genome sequencing reads into haplotype-specific reads sets. After assembling the reads of each haplotype, the contigs are scaffolded to chromosome level using a genetic map derived from the gametes. We assemble the two genomes of a diploid apricot tree based on whole-genome sequencing of 445 individual pollen grains. The two haplotype assemblies (N50: 25.5 and 25.8 Mb) feature a haplotyping precision of greater than 99% and are accurately scaffolded to chromosome-level.