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.

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 .

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.

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.

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.

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.

Mutagenesis of a Quintuple Mutant Impaired in Environmental Responses Reveals Roles for CHROMATIN REMODELING4 in the Arabidopsis Floral Transition.

Several pathways conferring environmental flowering responses in Arabidopsis (Arabidopsis thaliana) converge on developmental processes that mediate the floral transition in the shoot apical meristem. Many characterized mutations disrupt these environmental responses, but downstream developmental processes have been more refractory to mutagenesis. Here, we constructed a quintuple mutant impaired in several environmental pathways and showed that it possesses severely reduced flowering responses to changes in photoperiod and ambient temperature. RNA-sequencing (RNA-seq) analysis of the quintuple mutant showed that the expression of genes encoding gibberellin biosynthesis enzymes and transcription factors involved in the age pathway correlates with flowering. Mutagenesis of the quintuple mutant generated two late-flowering mutants, quintuple ems1 (qem1) and qem2 The mutated genes were identified by isogenic mapping and transgenic complementation. The qem1 mutant is an allele of the gibberellin 20-oxidase gene ga20ox2, confirming the importance of gibberellin for flowering in the absence of environmental responses. By contrast, qem2 is impaired in CHROMATIN REMODELING4 (CHR4), which has not been genetically implicated in floral induction. Using co-immunoprecipitation, RNA-seq, and chromatin immunoprecipitation sequencing, we show that CHR4 interacts with transcription factors involved in floral meristem identity and affects the expression of key floral regulators. Therefore, CHR4 mediates the response to endogenous flowering pathways in the inflorescence meristem to promote floral identity.

Bacterial Infection Disrupts Clock Gene Expression to Attenuate Immune Responses.

The circadian clock modulates immune responses in plants and animals; however, it is unclear how host-pathogen interactions affect the clock. Here we analyzed clock function in Arabidopsis thaliana mutants with defective immune responses and found that enhanced disease susceptibility 4 (eds4) displays alterations in several circadian rhythms. Mapping by sequencing revealed that EDS4 encodes the ortholog of NUCLEOPORIN 205, a core component of the inner ring of the nuclear pore complex (NPC). Consistent with the idea that the NPC specifically modulates clock function, we found a strong enrichment in core clock genes, as well as an increased nuclear to total mRNA accumulation, among genes that were differentially expressed in eds4 mutants. Interestingly, infection with Pseudomonas syringae in wild-type (WT) plants downregulated the expression of several morning core clock genes as early as 1 h post-infection, including all members of the NIGHT LIGHT-INDUCIBLE AND CLOCK-REGULATED (LNK) gene family, and this effect was attenuated in eds4. Furthermore, lnk mutants were more susceptible than the WT to P. syringae infection. These results indicate that bacterial infection, acting in part through the NPC, alters core clock gene expression and/or mRNA accumulation in a way that favors bacterial growth and disease susceptibility.

Mutations in the miR396 binding site of the growth-regulating factor gene VvGRF4 modulate inflorescence architecture in grapevine.

Bunch rot caused by Botrytis cinerea infections is a notorious problem in grapevine cultivation. To produce high quality fruits, grapevine plants are treated with fungicides, which is cost intensive and harmful to the environment. Conversely, loose cluster bunches show a considerably enhanced physical resilience to bunch diseases. With the aim to identify genetic determinants that modulate the development of bunch architecture, we have compared loose and compact 'Pinot noir' clones. Loose cluster architecture was found to be correlated with increased berry size, elongated rachis and elongated pedicels. Using transcriptome analysis in combination with whole genome sequencing, we have identified a growth-regulating factor gene, VvGRF4, upregulated and harbours heterozygous mutations in the loose cluster clones. At late stages of inflorescence development, the mRNA pools of loose cluster clones contain predominantly mRNAs derived from the mutated alleles, which are resistant to miR396 degradation. Expression of the VvGRF4 gene and its mutated variants in Arabidopsis demonstrates that it promotes pedicel elongation. Taken together, VvGRF4 modulates bunch architecture in grapevine 'Pinot noir' clones. This trait can be introduced into other cultivars using marker-assisted breeding or CRISPR-Cas9 technology. Related growth-regulating factors or other genes of the same pathway may have similar functions.

SyRI: finding genomic rearrangements and local sequence differences from whole-genome assemblies.

Genomic differences range from single nucleotide differences to complex structural variations. Current methods typically annotate sequence differences ranging from SNPs to large indels accurately but do not unravel the full complexity of structural rearrangements, including inversions, translocations, and duplications, where highly similar sequence changes in location, orientation, or copy number. Here, we present SyRI, a pairwise whole-genome comparison tool for chromosome-level assemblies. SyRI starts by finding rearranged regions and then searches for differences in the sequences, which are distinguished for residing in syntenic or rearranged regions. This distinction is important as rearranged regions are inherited differently compared to syntenic regions.

The Striking Flower-in-Flower Phenotype of Arabidopsis thaliana Nossen (No-0) is Caused by a Novel LEAFY Allele.

The transition to reproduction is a crucial step in the life cycle of any organism. In Arabidopsis thaliana the establishment of reproductive growth can be divided into two phases: Firstly, cauline leaves with axillary meristems are formed and internode elongation begins. Secondly, lateral meristems develop into flowers with defined organs. Floral shoots are usually determinate and suppress the development of lateral shoots. Here, we describe a transposon insertion mutant in the Nossen accession with defects in floral development and growth. Most strikingly is the outgrowth of stems from the axillary bracts of the primary flower carrying secondary flowers. Therefore, we named this mutant flower-in-flower (fif). However, the transposon insertion in the annotated gene is not the cause for the fif phenotype. By means of classical and genome sequencing-based mapping, the mutation responsible for the fif phenotype was found to be in the LEAFY gene. The mutation, a G-to-A exchange in the second exon of LEAFY, creates a novel lfy allele and results in a cysteine-to-tyrosine exchange in the α1-helix of LEAFY's DNA-binding domain. This exchange abolishes target DNA-binding, whereas subcellular localization and homomerization are not affected. To explain the strong fif phenotype against these molecular findings, several hypotheses are discussed.

Linked-read sequencing of gametes allows efficient genome-wide analysis of meiotic recombination.

Meiotic crossovers (COs) ensure proper chromosome segregation and redistribute the genetic variation that is transmitted to the next generation. Large populations and the demand for genome-wide, fine-scale resolution challenge existing methods for CO identification. Taking advantage of linked-read sequencing, we develop a highly efficient method for genome-wide identification of COs at kilobase resolution in pooled recombinants. We first test this method using a pool of Arabidopsis F2 recombinants, and recapitulate results obtained from the same plants using individual whole-genome sequencing. By applying this method to a pool of pollen DNA from an F1 plant, we establish a highly accurate CO landscape without generating or sequencing a single recombinant plant. The simplicity of this approach enables the simultaneous generation and analysis of multiple CO landscapes, accelerating the pace at which mechanisms for the regulation of recombination can be elucidated through efficient comparisons of genotypic and environmental effects on recombination.

Interspecies association mapping links reduced CG to TG substitution rates to the loss of gene-body methylation.

Comparative genomics can unravel the genetic basis of species differences; however, successful reports on quantitative traits are still scarce. Here we present genome assemblies of 31 so-far unassembled Brassicaceae plant species and combine them with 16 previously published assemblies to establish the Brassicaceae Diversity Panel. Using a new interspecies association strategy for quantitative traits, we found a so-far unknown association between the unexpectedly high variation in CG to TG substitution rates in genes and the absence of CHROMOMETHYLASE3 (CMT3) orthologues. Low substitution rates were associated with the loss of CMT3, while species with conserved CMT3 orthologues showed high substitution rates. Species without CMT3 also lacked gene-body methylation (gbM), suggesting an evolutionary trade-off between the unknown function of gbM and low substitution rates in Brassicaceae, possibly due to low mutability of non-methylated cytosines.

findGSE: estimating genome size variation within human and Arabidopsis using k-mer frequencies.

Motivation: Analyzing k-mer frequencies in whole-genome sequencing data is becoming a common method for estimating genome size (GS). However, it remains uninvestigated how accurate the method is, especially if it can capture intra-species GS variation. Results: We present findGSE, which fits skew normal distributions to k-mer frequencies to estimate GS. findGSE outperformed existing tools in an extensive simulation study. Estimating GSs of 89 Arabidopsis thaliana accessions, findGSE showed the highest capability in capturing GS variations. In an application with 71 female and 71 male human individuals, findGSE delivered an average of 3039 Mb as haploid human GS, while female genomes were on average 41 Mb larger than male genomes, in astonishing agreement with size difference of the X and Y chromosomes. Further analysis showed that human GS variations link to geographical patterns and significant differences between populations, which can be explained by variable abundances of LINE-1 retrotransposons. Availability and implementation: R package of findGSE is freely available at https://github.com/schneebergerlab/findGSE and supported on linux and Mac systems. Contact: schneeberger@mpipz.mpg.de. Supplementary information: Supplementary data are available at Bioinformatics online.

Brassinosteroids Dominate Hormonal Regulation of Plant Thermomorphogenesis via BZR1.

Thermomorphogenesis is defined as the suite of morphological changes that together are likely to contribute to adaptive growth acclimation to usually elevated ambient temperature [1, 2]. While many details of warmth-induced signal transduction are still elusive, parallels to light signaling recently became obvious (reviewed in [3]). It involves photoreceptors that can also sense changes in ambient temperature [3-5] and act, for example, by repressing protein activity of the central integrator of temperature information PHYTOCHROME-INTERACTING FACTOR 4 (PIF4 [6]). In addition, PIF4 transcript accumulation is tightly controlled by the evening complex member EARLY FLOWERING 3 [7, 8]. According to the current understanding, PIF4 activates growth-promoting genes directly but also via inducing auxin biosynthesis and signaling, resulting in cell elongation. Based on a mutagenesis screen in the model plant Arabidopsis thaliana for mutants with defects in temperature-induced hypocotyl elongation, we show here that both PIF4 and auxin function depend on brassinosteroids. Genetic and pharmacological analyses place brassinosteroids downstream of PIF4 and auxin. We found that brassinosteroids act via the transcription factor BRASSINAZOLE RESISTANT 1 (BZR1), which accumulates in the nucleus at high temperature, where it induces expression of growth-promoting genes. Furthermore, we show that at elevated temperature BZR1 binds to the promoter of PIF4, inducing its expression. These findings suggest that BZR1 functions in an amplifying feedforward loop involved in PIF4 activation. Although numerous negative regulators of PIF4 have been described, we identify BZR1 here as a true temperature-dependent positive regulator of PIF4, acting as a major growth coordinator.

SHOREmap v3.0: fast and accurate identification of causal mutations from forward genetic screens.

Whole-genome resequencing of pools of recombinant mutant genomes allows direct linking of phenotypic traits to causal mutations. Such analysis, called mapping-by-sequencing, combines classical genetic mapping and next-generation sequencing by relying on selection-induced patterns within genome-wide allele frequency (AF) in pooled genomes. Mapping-by-sequencing can be performed with computational tools such as SHOREmap. Previous versions of SHOREmap, however, did not implement standardized analyses, but were specifically designed for particular experimental settings. Here, we introduce the usage of a novel and advanced implementation of SHOREmap (version 3.0), including several new features like file readers for commonly used file formats, SNP marker selection, and a stable calculation of mapping intervals. SHOREmap can be downloaded at shoremap.org.

Rapid identification of a natural knockout allele of ARMADILLO REPEAT-CONTAINING KINESIN1 that causes root hair branching by mapping-by-sequencing.

In Arabidopsis (Arabidopsis thaliana), branched root hairs are an indicator of defects in root hair tip growth. Among 62 accessions, one accession (Heiligkreuztal2 [HKT2.4]) displayed branched root hairs, suggesting that this accession carries a mutation in a gene of importance for tip growth. We determined 200- to 300-kb mapping intervals using a mapping-by-sequencing approach of F2 pools from crossings of HKT2.4 with three different accessions. The intersection of these mapping intervals was 80 kb in size featuring not more than 36 HKT2.4-specific single nucleotide polymorphisms, only two of which changed the coding potential of genes. Among them, we identified the causative single nucleotide polymorphism changing a splicing site in ARMADILLO REPEAT-CONTAINING KINESIN1. The applied strategies have the potential to complement statistical methods in high-throughput phenotyping studies using different natural accessions to identify causative genes for distinct phenotypes represented by only one or a few accessions.

Identification of a new export signal in Plasmodium yoelii: identification of a new exportome.

Development of the erythrocytic malaria parasite requires targeting of parasite proteins into multiple compartments located within and beyond the parasite confine. Beyond the PEXEL/VTS pathway and its characterized players, increasing amount of evidence has highlighted the existence of proteins exported using alternative export-signal(s)/pathway(s); hence, the exportomes currently predicted are incomplete. The nature of these exported proteins which could have a prominent role in most of the Plasmodium species remains elusive. Using P. yoelii variant proteins, we identified a signal associated to lipophilic region that mediates export of P. yoelii proteins. This non-PEXEL signal termed PLASMED is defined by semi-conserved residues and possibly a secondary structure. In vivo characterization of exported-proteins indicated that PLASMED is a bona fide export-signal that allowed us to identify an unseen P. yoelii exportome. The repertoire of the newly predicted exported proteins opens up perspectives for unravelling the remodelling of the host-cell by the parasite, against which new therapies could be elaborated.

Exhaled droplets due to talking and coughing.

Respiratory infections can be spread via 'contact' with droplets from expiratory activities such as talking, coughing and sneezing, and also from aerosol-generating clinical procedures. Droplet sizes predominately determine the times they can remain airborne, the possibility of spread of infectious diseases and thus the strategies for controlling the infections. While significant inconsistencies exist between the existing measured data on respiratory droplets generated during expiratory activities, a food dye was used in the mouth during measurements of large droplets, which made the expiratory activities 'unnatural'. We carried out a series of experiments using glass slides and a microscope as well as an aerosol spectrometer to measure the number and size of respiratory droplets produced from the mouth of healthy individuals during talking and coughing with and without a food dye. The total mass of respiratory droplets was measured using a mask, plastic bag with tissue and an electronic balance with a high precision. Considerable subject variability was observed and the average size of droplets captured using glass slides and microscope was about 50-100 microm. Smaller droplets were also detected by the aerosol spectrometer. More droplets seemed to be generated when a food dye was used.