Publications

Wheat genetic resources hold the diversity required to mitigate agricultural challenges from climate change and reduced inputs. Using DArTseq, we genotype 461 wheat landraces and cultivars and evaluate them for powdery mildew resistance. By developing a k-mer-based GWAS approach with fully assembled genomes of Triticum aestivum and its progenitors, we uncover 25% more resistance-associated k-mers than single-reference methods, outperforming SNP-based GWAS in both loci detection and mapping precision. In total, we detect 34 powdery mildew resistance loci, including 27 potentially novel regions. Our approach underscores the importance of integrating multiple reference genomes to unlock the potential of wheat germplasm.

Much of what we know about eukaryotic transcription stems from animals and yeast; however, plants evolved separately for over a billion years, leaving ample time for divergence in transcriptional regulation. Here we set out to elucidate fundamental properties of cis-regulatory sequences in plants. Using massively parallel reporter assays across four plant species, we demonstrate the central role of sequences downstream of the transcription start site (TSS) in transcriptional regulation. Unlike animal enhancers that are position independent, plant regulatory elements depend on their position, as altering their location relative to the TSS significantly affects transcription. We highlight the importance of the region downstream of the TSS in regulating transcription by identifying a DNA motif that is conserved across vascular plants and is sufficient to enhance gene expression in a dose-dependent manner. The identification of a large number of position-dependent enhancers points to fundamental differences in gene regulation between plants and animals.

Fertilization in Arabidopsis (Arabidopsis thaliana) is a highly coordinated process that begins with a pollen tube delivering the 2 sperm cells into the embryo sac. Each sperm cell can then fertilize either the egg or the central cell to initiate embryo or endosperm development, respectively. The success of this double fertilization process requires a tight cell cycle synchrony between the male and female gametes to allow karyogamy (nuclei fusion). However, the cell cycle status of the male and female gametes during fertilization remains elusive as DNA quantification and DNA replication assays have given conflicting results. Here, to reconcile these results, we quantified the DNA replication state by DNA sequencing and performed microscopic analyses of fluorescent markers covering all phases of the cell cycle. We show that male and female Arabidopsis gametes are both arrested prior to DNA replication at maturity and initiate their DNA replication only during fertilization.

Repeated herbicide applications in agricultural fields exert strong selection on weeds such as blackgrass (Alopecurus myosuroides), which is a major threat for temperate climate cereal crops. This inadvertent selection pressure provides an opportunity for investigating the underlying genetic mechanisms and evolutionary processes of rapid adaptation, which can occur both through mutations in the direct targets of herbicides and through changes in other, often metabolic, pathways, known as non-target-site resistance. How much target-site resistance (TSR) relies on de novo mutations vs. standing variation is important for developing strategies to manage herbicide resistance. We first generated a chromosome-level reference genome for A. myosuroides for population genomic studies of herbicide resistance and genome-wide diversity across Europe in this species. Next, through empirical data in the form of highly accurate long-read amplicons of alleles encoding acetyl-CoA carboxylase (ACCase) and acetolactate synthase (ALS) variants, we showed that most populations with resistance due to TSR mutations-23 out of 27 and six out of nine populations for ACCase and ALS, respectively-contained at least two TSR haplotypes, indicating that soft sweeps are the norm. Finally, through forward-in-time simulations, we inferred that TSR is likely to mainly result from standing genetic variation, with only a minor role for de novo mutations.

Structural variants and presence/absence polymorphisms are common in plant genomes, yet they are routinely overlooked in genome-wide association studies (GWAS). Here, we expand the type of genetic variants detected in GWAS to include major deletions, insertions and rearrangements. We first use raw sequencing data directly to derive short sequences, k-mers, that mark a broad range of polymorphisms independently of a reference genome. We then link k-mers associated with phenotypes to specific genomic regions. Using this approach, we reanalyzed 2,000 traits in Arabidopsis thaliana, tomato and maize populations. Associations identified with k-mers recapitulate those found with SNPs, but with stronger statistical support. Importantly, we discovered new associations with structural variants and with regions missing from reference genomes. Our results demonstrate the power of performing GWAS before linking sequence reads to specific genomic regions, which allows the detection of a wider range of genetic variants responsible for phenotypic variation.

Splicing expands, reshapes, and regulates the transcriptome of eukaryotic organisms. Despite its importance, key questions remain unanswered, including the following: Can splicing evolve when organisms adapt to new challenges? How does evolution optimize inefficiency of introns' splicing and of the splicing machinery? To explore these questions, we evolved yeast cells that were engineered to contain an inefficiently spliced intron inside a gene whose protein product was under selection for an increased expression level. We identified a combination of mutations in Cis (within the gene of interest) and in Trans (in mRNA-maturation machinery). Surprisingly, the mutations in Cis resided outside of known intronic functional sites and improved the intron's splicing efficiency potentially by easing tight mRNA structures. One of these mutations hampered a protein's domain that was not under selection, demonstrating the evolutionary flexibility of multi-domain proteins as one domain functionality was improved at the expense of the other domain. The Trans adaptations resided in two proteins, Npl3 and Gbp2, that bind pre-mRNAs and are central to their maturation. Interestingly, these mutations either increased or decreased the affinity of these proteins to mRNA, presumably allowing faster spliceosome recruitment or increased time before degradation of the pre-mRNAs, respectively. Altogether, our work reveals various mechanistic pathways toward optimizations of intron splicing to ultimately adapt gene expression patterns to novel demands.

DNA replication introduces a dosage imbalance between early and late replicating genes. In budding yeast, buffering gene expression against this imbalance depends on marking replicated DNA by H3K56 acetylation (H3K56ac). Whether additional processes are required for suppressing transcription from H3K56ac-labeled DNA remains unknown. Here, using a database-guided candidate screen, we find that COMPASS, the H3K4 methyltransferase, and its upstream effector, PAF1C, act downstream of H3K56ac to buffer expression. Replicated genes show reduced abundance of the transcription activating mark H3K4me3 and accumulate the transcription inhibitory mark H3K4me2 near transcription start sites. Notably, in hydroxyurea-exposed cells, the S phase checkpoint stabilizes H3K56ac and becomes essential for buffering. We suggest that H3K56ac suppresses transcription of replicated genes by interfering with post-replication recovery of epigenetic marks and assign a new function for the S phase checkpoint in stabilizing this mechanism during persistent dosage imbalance.

Sphingolipids (SLs) are essential components of cell membranes and are broad-range bioactive signaling molecules. SL levels must be tightly regulated as imbalances affect cellular function and contribute to pathologies ranging from neurodegenerative and metabolic disorders to cancer and aging. Deciphering how SL homeostasis is maintained and uncovering new regulators is required for understanding lipid biology and for identifying new targets for therapeutic interventions. Here we combine omics technologies to identify the changes of the transcriptome, proteome, and phosphoproteome in the yeast Saccharomyces cerevisiae upon SL depletion induced by myriocin. Surprisingly, while SL depletion triggers important changes in the expression of regulatory proteins involved in SL homeostasis, the most dramatic regulation occurs at the level of the phosphoproteome, suggesting that maintaining SL homeostasis demands rapid responses. To discover which of the phosphoproteomic changes are required for the cell's first-line response to SL depletion, we overlaid our omics results with systematic growth screens for genes required during growth in myriocin. By following the rate of SL biosynthesis in those candidates that are both affecting growth and are phosphorylated in response to the drug, we uncovered Atg9, Stp4, and Gvp36 as putative new regulators of SL homeostasis.

DNA replication perturbs the dosage balance between genes that replicate early during S phase and those that replicate late. If propagated to influence protein content, this dosage imbalance could influence cellular functions. In bacteria, mechanisms have evolved to use this imbalance to tune certain processes with the rate of cell growth. By contrast, eukaryotes buffer this dosage imbalance to ensure gene expression homeostasis also during S phase. Here, we outline classical and more recent studies describing how different organisms deal with this replication-dependent dosage imbalance, and describe recent results linking the eukaryotic buffering mechanism to replication-dependent histone acetylation. Finally, we discuss the possible implications of this buffering mechanism and speculate why it is specific to eukaryote cells.

Chromatin is composed of DNA and histones, which provide a unified platform for regulating DNA-related processes, mostly through their post- Translational modification. During DNA replication, histone arrangement is perturbed, first to allow progression of DNA polymerase and then during repackaging of the replicated DNA. To study how DNA replication influences the pattern of histone modification, we followed the cell-cycle dynamics of 10 histone marks in budding yeast. We find that histones deposited on newly replicatedDNAare modified at different rates: While some marks appear immediately upon replication (e.g., H4K16ac, H3K4me1), others increase with transcription-dependent delays (e.g., H3K4me3, H3K36me3). Notably, H3K9ac was deposited as a wave preceding the replication fork by ∼5-6 kb. This replication-guided H3K9ac was fully dependent on the acetyltransferase Rtt109, while expression-guided H3K9ac was deposited by Gcn5. Further, topoisomerase depletion intensified H3K9ac in front of the replication fork and in sites where RNA polymerase II was trapped, suggesting supercoiling stresses trigger H3K9 acetylation. Our results assign complementary roles for DNA replication and gene expression in defining the pattern of histone modification.

Genome replication introduces a stepwise increase in the DNA template available for transcription. Genes replicated early in S phase experience this increase before latereplicating genes, raising the question of how expression levels are affected by DNA replication. We show that in budding yeast, messenger RNA (mRNA) synthesis rate is buffered against changes in gene dosage during S phase. This expression homeostasis depends on acetylation of H3 on its internal K56 site by Rtt109/Asf1. Deleting these factors, mutating H3K56 or up-regulating its deacetylation, increases gene expression in S phase in proportion to gene replication timing. Therefore, H3K56 acetylation on newly deposited histones reduces transcription efficiency from replicated DNA, complementing its role in guarding genome stability. Our study provides molecular insight into the mechanism maintaining expression homeostasis during DNA replication.

Chromatin can function as an integrator of DNA-related processes, allowing communication, for example, between DNA replication and gene transcription. Such communication is needed to overcome the gene-dosage imbalance introduced during DNA replication, when certain genes are replicated prior to others. Increased transcription of early replicating genes could alter regulatory balances. This does not occur, suggesting a mechanism that suppresses expression from newly replicated DNA. Critical to this buffering is Rtt109, which acetylates the internal K56 residue of newly synthesized histone H3 prior to incorporation onto DNA. H3K56ac distinguishes replicated from non-replicated DNA, communicating this information to the transcription machinery to ensure expression homeostasis during S phase.

4sUDRB-seq separately measures, on a genomic scale, the distinct contributions of transcription elongation speed and rate of RNA polymerase II (Pol II) transition into active elongation (TAE) to the overall mRNA production rate. It uses reversible inhibition of transcription elongation with 5,6-dichloro-1-D-ribofuranosylbenzimidazole (DRB), combined with a pulse of 4-thiouridine (4sU), to tag newly transcribed RNA. After DRB removal, cells are collected at several time points, and tagged RNA is biotinylated, captured on streptavidin beads and sequenced. 4sUDRB-seq enables the comparison of elongation speeds between different developmental stages or different cell types, and it allows the impact of specific transcription factors on transcription elongation speed versus TAE to be studied. RNA preparation takes ∼4 d to complete, with deep sequencing requiring an additional ∼4-11 d plus 1-3 d for bioinformatics analysis. The experimental protocol requires basic molecular biology skills, whereas data analysis requires knowledge in bioinformatics, particularly MATLAB and the Linux environment.

Various histone modifications decorate nucleosomes within transcribed genes. Among these, monoubiquitylation of histone H2B (H2Bub1) and methylation of histone H3 on lysines 36 (H3K36me2/3) and 79 (H3K79me2/3) correlate positively with gene expression. By measuring the progression of the transcriptional machinery along genes within live cells, we now report that H2B monoubiquitylation occurs cotranscriptionally and accurately reflects the advance of RNA polymerase II (Pol II). In contrast, H3K36me3 and H3K79me2 are less dynamic and represent Pol II movement less faithfully. High-resolution ChIP-seq reveals that H2Bub1 levels are selectively reduced at exons and decrease in an exon-dependent stepwise manner toward the 39 end of genes. Exonic depletion of H2Bub1 in gene bodies is highly correlated with Pol II pausing at exons, suggesting elongation rate changes associated with intron-exon structure. In support of this notion, H2Bub1 levels were found to be significantly correlatedwith transcription elongation rates measured in various cell lines.Overall, our data shed light on the organization of H2Bub1 within transcribed genes and single out H2Bub1 as a reliable marker for ongoing transcription elongation.

Various histone modifications are widely associated with gene expression, but their functional selectivity at individual genes remains to be characterized. Here, we identify widespread differences between genome-wide patterns of two prominent marks, H3K9ac and H3K4me3, in budding yeasts. As well as characteristic gene profiles, relative modification levels vary significantly amongst genes, irrespective of expression. Interestingly, we show that these differences couple to contrasting features: higher methylation to essential, periodically expressed, 'DPN' (Depleted Proximal Nucleosome) genes, and higher acetylation to non-essential, responsive, 'OPN' (Occupied Proximal Nucleosome) genes. Thus, H3K4me3 may generally associate with expression stability, and H3K9ac, with variability. To evaluate this notion, we examine their association with expression divergence between the closely related species, S. cerevisiae and S. paradoxus. Although individually well conserved at orthologous genes, changes between modifications are mostly uncorrelated, indicating largely non-overlapping regulatory mechanisms. Notably, we find that inter-species differences in methylation, but not acetylation, are well correlated with expression changes, thereby proposing H3K4me3 as a candidate regulator of expression divergence. Taken together, our results suggest distinct evolutionary roles for expression-linked modifications, wherein H3K4me3 may contribute to stabilize average expression, whilst H3K9ac associates with more indirect aspects such as responsiveness.

Cells adapt to environmental changes through genetic mutations that stabilize novel phenotypes. Often, this adaptation involves regulatory changes which modulate gene expression. In the budding yeast, ribosomal-related gene expression correlates with cell growth rate across different environments. To examine whether the same relationship between gene expression and growth rate is observed also across natural populations, we measured gene expression, growth rate and ethanol production of twenty-four wild type yeast strains originating from diverse habitats, grown on the pentose sugar xylulose. We found that expression of ribosome-related genes did not correlate with growth rate. Rather, growth rate was correlated with the expression of amino acid biosynthesis genes. Searching other databases, we observed a similar correlation between growth rate and amino-acid biosyntehsis genes in a library of gene deletions. We discuss the implications of our results for understanding how cells coordinate their translation capacity with available nutrient resources.

Although transcriptional elongation by RNA polymerase II is coupled with many RNA-related processes, genomewide elongation rates remain unknown. We describe a method, called 4sUDRB-seq, based on reversible inhibition of transcription elongation coupled with tagging newly transcribed RNA with 4-thiouridine and high throughput sequencing to measure simultaneously with high confidence genome-wide transcription elongation rates in cells. We find that most genes are transcribed at about 3.5 Kb/min, with elongation rates varying between 2 Kb/min and 6 Kb/min. 4sUDRB-seq can facilitate genomewide exploration of the involvement of specific elongation factors in transcription and the contribution of deregulated transcription elongation to various pathologies.

Gene expression shows a significant variation (noise) between genetically identical cells. Noise depends on the gene expression process regulated by the chromatin environment. We screened for chromatin factors that modulate noise in S. cerevisiae and analyzed the results using a theoretical model that infers regulatory mechanisms from the noise versus mean relationship. Distinct activities of the Rpd3(L) and Set3 histone deacetylase complexes were predicted. Both HDACs repressed expression. Yet, Rpd3(L)C decreased the frequency of transcriptional bursts, while Set3C decreased the burst size, as did H2B monoubiquitination (ubH2B). We mapped the acetylation of H3 lysine 9 (H3K9ac) upon deletion of multiple subunits of Set3C and Rpd3(L)C and of ubH2B effectors. ubH2B and Set3C appear to function in the same pathway to reduce the probability that an elongating PolII produces a functional transcript (PolII processivity), while Rpd3(L)C likely represses gene expression at a step preceding elongation.