Pangenome graph layout by Path-Guided Stochastic Gradient Descent.

Motivation: The increasing availability of complete genomes demands for models to study genomic variability within entire populations. Pangenome graphs capture the full genomic similarity and diversity between multiple genomes. In order to understand them, we need to see them. For visualization, we need a human readable graph layout: A graph embedding in low (e.g. two) dimensional depictions. Due to a pangenome graph's potential excessive size, this is a significant challenge. Results: In response, we introduce a novel graph layout algorithm: the Path-Guided Stochastic Gradient Descent (PG-SGD). PG-SGD uses the genomes, represented in the pangenome graph as paths, as an embedded positional system to sample genomic distances between pairs of nodes. This avoids the quadratic cost seen in previous versions of graph drawing by Stochastic Gradient Descent (SGD). We show that our implementation efficiently computes the low dimensional layouts of gigabase-scale pangenome graphs, unveiling their biological features. Availability: We integrated PG-SGD in ODGI which is released as free software under the MIT open source license. Source code is available at https://github.com/pangenome/odgi.

Building pangenome graphs.

Pangenome graphs can represent all variation between multiple genomes, but existing methods for constructing them are biased due to reference-guided approaches. In response, we have developed PanGenome Graph Builder (PGGB), a reference-free pipeline for constructing unbi-ased pangenome graphs. PGGB uses all-to-all whole-genome alignments and learned graph embeddings to build and iteratively refine a model in which we can identify variation, measure conservation, detect recombination events, and infer phylogenetic relationships.

Impaired KIN10 function restores developmental defects in the Arabidopsis trehalose 6-phosphate synthase1 (tps1) mutant.

Sensing carbohydrate availability is essential for plants to coordinate their growth and development. In Arabidopsis thaliana, TREHALOSE 6-PHOSPHATE SYNTHASE 1 (TPS1) and its product, trehalose 6-phosphate (T6P), are important for the metabolic control of development. tps1 mutants are embryo-lethal and unable to flower when embryogenesis is rescued. T6P regulates development in part through inhibition of SUCROSE NON-FERMENTING1 RELATED KINASE1 (SnRK1). Here, we explored the role of SnRK1 in T6P-mediated plant growth and development using a combination of a mutant suppressor screen and genetic, cellular and transcriptomic approaches. We report nonsynonymous amino acid substitutions in the catalytic KIN10 and regulatory SNF4 subunits of SnRK1 that can restore both embryogenesis and flowering of tps1 mutant plants. The identified SNF4 point mutations disrupt the interaction with the catalytic subunit KIN10. Contrary to the common view that the two A. thaliana SnRK1 catalytic subunits act redundantly, we found that loss-of-function mutations in KIN11 are unable to restore embryogenesis and flowering, highlighting the important role of KIN10 in T6P signalling.

MethylScore, a pipeline for accurate and context-aware identification of differentially methylated regions from population-scale plant whole-genome bisulfite sequencing data.

Whole-genome bisulfite sequencing (WGBS) is the standard method for profiling DNA methylation at single-nucleotide resolution. Different tools have been developed to extract differentially methylated regions (DMRs), often built upon assumptions from mammalian data. Here, we present MethylScore, a pipeline to analyse WGBS data and to account for the substantially more complex and variable nature of plant DNA methylation. MethylScore uses an unsupervised machine learning approach to segment the genome by classification into states of high and low methylation. It processes data from genomic alignments to DMR output and is designed to be usable by novice and expert users alike. We show how MethylScore can identify DMRs from hundreds of samples and how its data-driven approach can stratify associated samples without prior information. We identify DMRs in the A. thaliana 1,001 Genomes dataset to unveil known and unknown genotype-epigenotype associations .

Methylation of CpG 5962 in L1 of the human papillomavirus 16 genome as a potential predictive marker for viral persistence: A prospective large cohort study using cervical swab samples.

Several studies have demonstrated that the viral genome can be methylated by the host cell during progression from persistent infection to cervical cancer. The aim of this study was to investigate whether methylation at a specific site could predict the development of viral persistence and whether viral load shows a correlation with specific methylation patterns. HPV16-positive samples from women aged 20-29 years (n = 99) with a follow-up time of 13 years, were included from a Danish cohort comprising 11 088 women. Viral load was measured by real-time PCR and methylation status was determined for 39 CpG sites in the upstream regulatory region (URR), E6/E7, and L1 region of HPV16 by next-generation sequencing. Participants were divided into two groups according to whether they were persistently (≥ 24 months) or transiently HPV16 infected. The general methylation status was significantly different between women with a persistent and women with a transient infection outcome (P = .025). One site located in L1 (nt. 5962) was statistically significantly (P = .00048) different in the methylation status after correction using the Holm-Sidak method (alpha = 0.05). Correlation analyses of samples from HPV16 persistently infected women suggest that methylation is higher although viral load is lower. This study indicates that methylation at position 5962 of the HPV16 genome within the L1 gene might be a predictive marker for the development of a persistent HPV16 infection.

RPW8/HR repeats control NLR activation in Arabidopsis thaliana.

In many plant species, conflicts between divergent elements of the immune system, especially nucleotide-binding oligomerization domain-like receptors (NLR), can lead to hybrid necrosis. Here, we report deleterious allele-specific interactions between an NLR and a non-NLR gene cluster, resulting in not one, but multiple hybrid necrosis cases in Arabidopsis thaliana. The NLR cluster is RESISTANCE TO PERONOSPORA PARASITICA 7 (RPP7), which can confer strain-specific resistance to oomycetes. The non-NLR cluster is RESISTANCE TO POWDERY MILDEW 8 (RPW8) / HOMOLOG OF RPW8 (HR), which can confer broad-spectrum resistance to both fungi and oomycetes. RPW8/HR proteins contain at the N-terminus a potential transmembrane domain, followed by a specific coiled-coil (CC) domain that is similar to a domain found in pore-forming toxins MLKL and HET-S from mammals and fungi. C-terminal to the CC domain is a variable number of 21- or 14-amino acid repeats, reminiscent of regulatory 21-amino acid repeats in fungal HET-S. The number of repeats in different RPW8/HR proteins along with the sequence of a short C-terminal tail predicts their ability to activate immunity in combination with specific RPP7 partners. Whether a larger or smaller number of repeats is more dangerous depends on the specific RPW8/HR autoimmune risk variant.

Long-term balancing selection drives evolution of immunity genes in Capsella.

Genetic drift is expected to remove polymorphism from populations over long periods of time, with the rate of polymorphism loss being accelerated when species experience strong reductions in population size. Adaptive forces that maintain genetic variation in populations, or balancing selection, might counteract this process. To understand the extent to which natural selection can drive the retention of genetic diversity, we document genomic variability after two parallel species-wide bottlenecks in the genus Capsella. We find that ancestral variation preferentially persists at immunity related loci, and that the same collection of alleles has been maintained in different lineages that have been separated for several million years. By reconstructing the evolution of the disease-related locus MLO2b, we find that divergence between ancient haplotypes can be obscured by referenced based re-sequencing methods, and that trans-specific alleles can encode substantially diverged protein sequences. Our data point to long-term balancing selection as an important factor shaping the genetics of immune systems in plants and as the predominant driver of genomic variability after a population bottleneck.

The rate and potential relevance of new mutations in a colonizing plant lineage.

By following the evolution of populations that are initially genetically homogeneous, much can be learned about core biological principles. For example, it allows for detailed studies of the rate of emergence of de novo mutations and their change in frequency due to drift and selection. Unfortunately, in multicellular organisms with generation times of months or years, it is difficult to set up and carry out such experiments over many generations. An alternative is provided by "natural evolution experiments" that started from colonizations or invasions of new habitats by selfing lineages. With limited or missing gene flow from other lineages, new mutations and their effects can be easily detected. North America has been colonized in historic times by the plant Arabidopsis thaliana, and although multiple intercrossing lineages are found today, many of the individuals belong to a single lineage, HPG1. To determine in this lineage the rate of substitutions-the subset of mutations that survived natural selection and drift-, we have sequenced genomes from plants collected between 1863 and 2006. We identified 73 modern and 27 herbarium specimens that belonged to HPG1. Using the estimated substitution rate, we infer that the last common HPG1 ancestor lived in the early 17th century, when it was most likely introduced by chance from Europe. Mutations in coding regions are depleted in frequency compared to those in other portions of the genome, consistent with purifying selection. Nevertheless, a handful of mutations is found at high frequency in present-day populations. We link these to detectable phenotypic variance in traits of known ecological importance, life history and growth, which could reflect their adaptive value. Our work showcases how, by applying genomics methods to a combination of modern and historic samples from colonizing lineages, we can directly study new mutations and their potential evolutionary relevance.

A Role for the F-Box Protein HAWAIIAN SKIRT in Plant microRNA Function.

As regulators of gene expression in multicellular organisms, microRNAs (miRNAs) are crucial for growth and development. Although a plethora of factors involved in their biogenesis and action in Arabidopsis (Arabidopsis thaliana) has been described, these processes and their fine-tuning are not fully understood. Here, we used plants expressing an artificial miRNA target mimic (MIM) to screen for negative regulators of miR156. We identified a new mutant allele of the F-box gene HAWAIIAN SKIRT (HWS; At3G61590), hws-5, as a suppressor of the MIM156-induced developmental and molecular phenotypes. In hws plants, levels of some endogenous miRNAs are increased and their mRNA targets decreased. Plants constitutively expressing full-length HWS-but not a truncated version lacking the F-box domain-display morphological and molecular phenotypes resembling those of mutants defective in miRNA biogenesis and activity. In combination with such mutants, hws loses its delayed floral organ abscission ("skirt") phenotype, suggesting epistasis. Also, the hws transcriptome profile partially resembles those of well-known miRNA mutants hyl1-2, se-3, and ago1-27, pointing to a role in a common pathway. We thus propose HWS as a novel, F-box dependent factor involved in miRNA function.

Arabidopsis BTB/POZ protein-dependent PENETRATION3 trafficking and disease susceptibility.

The outermost cell layer of plant roots (epidermis) constantly encounters environmental challenges. The epidermal outer plasma membrane domain harbours the PENETRATION3 (PEN3)/ABCG36/PDR8 ATP-binding cassette transporter that confers non-host resistance to several pathogens. Here, we show that the Arabidopsis ENDOPLASMIC RETICULUM-ARRESTED PEN3 (EAP3) BTB/POZ-domain protein specifically mediates PEN3 exit from the endoplasmic reticulum and confers resistance to a root-penetrating fungus, providing prime evidence for BTB/POZ-domain protein-dependent membrane trafficking underlying disease resistance.

Assessing Distribution and Variation of Genome-Wide DNA Methylation Using Short-Read Sequencing.

DNA methylation is one of the most prominent epigenetic marks and is particularly complex in plant genomes. Whole-genome bisulfite sequencing (WGBS) using short reads has become the standard tool to study genome-wide patterns of DNA methylation. The goal of the present protocol is to enable readers to perform WGBS on both the wet lab and the computational side. We briefly outline important steps in bisulfite library preparation, then focus on the different aspects of DNA methylation analysis, from read mapping to identifying biologically relevant differential methylation between different samples.

On the Origin of De Novo Genes in Arabidopsis thaliana Populations.

De novo genes, which originate from ancestral nongenic sequences, are one of the most important sources of protein-coding genes. This origination process is crucial for the adaptation of organisms. However, how de novo genes arise and become fixed in a population or species remains largely unknown. Here, we identified 782 de novo genes from the model plant Arabidopsis thaliana and divided them into three types based on the availability of translational evidence, transcriptional evidence, and neither transcriptional nor translational evidence for their origin. Importantly, by integrating multiple types of omics data, including data from genomes, epigenomes, transcriptomes, and translatomes, we found that epigenetic modifications (DNA methylation and histone modification) play an important role in the origination process of de novo genes. Intriguingly, using the transcriptomes and methylomes from the same population of 84 accessions, we found that de novo genes that are transcribed in approximately half of the total accessions within the population are highly methylated, with lower levels of transcription than those transcribed at other frequencies within the population. We hypothesized that, during the origin of de novo gene alleles, those neutralized to low expression states via DNA methylation have relatively high probabilities of spreading and becoming fixed in a population. Our results highlight the process underlying the origin of de novo genes at the population level, as well as the importance of DNA methylation in this process.

Hyperosmotic stress memory in Arabidopsis is mediated by distinct epigenetically labile sites in the genome and is restricted in the male germline by DNA glycosylase activity.

Inducible epigenetic changes in eukaryotes are believed to enable rapid adaptation to environmental fluctuations. We have found distinct regions of the Arabidopsis genome that are susceptible to DNA (de)methylation in response to hyperosmotic stress. The stress-induced epigenetic changes are associated with conditionally heritable adaptive phenotypic stress responses. However, these stress responses are primarily transmitted to the next generation through the female lineage due to widespread DNA glycosylase activity in the male germline, and extensively reset in the absence of stress. Using the CNI1/ATL31 locus as an example, we demonstrate that epigenetically targeted sequences function as distantly-acting control elements of antisense long non-coding RNAs, which in turn regulate targeted gene expression in response to stress. Collectively, our findings reveal that plants use a highly dynamic maternal 'short-term stress memory' with which to respond to adverse external conditions. This transient memory relies on the DNA methylation machinery and associated transcriptional changes to extend the phenotypic plasticity accessible to the immediate offspring.

Targeting Metabolic Symbiosis to Overcome Resistance to Anti-angiogenic Therapy.

Despite the approval of several anti-angiogenic therapies, clinical results remain unsatisfactory, and transient benefits are followed by rapid tumor recurrence. Here, we demonstrate potent anti-angiogenic efficacy of the multi-kinase inhibitors nintedanib and sunitinib in a mouse model of breast cancer. However, after an initial regression, tumors resume growth in the absence of active tumor angiogenesis. Gene expression profiling of tumor cells reveals metabolic reprogramming toward anaerobic glycolysis. Indeed, combinatorial treatment with a glycolysis inhibitor (3PO) efficiently inhibits tumor growth. Moreover, tumors establish metabolic symbiosis, illustrated by the differential expression of MCT1 and MCT4, monocarboxylate transporters active in lactate exchange in glycolytic tumors. Accordingly, genetic ablation of MCT4 expression overcomes adaptive resistance against anti-angiogenic therapy. Hence, targeting metabolic symbiosis may be an attractive avenue to avoid resistance development to anti-angiogenic therapy in patients.

KH domain protein RCF3 is a tissue-biased regulator of the plant miRNA biogenesis cofactor HYL1.

The biogenesis of microRNAs (miRNAs), which regulate mRNA abundance through posttranscriptional silencing, comprises multiple well-orchestrated processing steps. We have identified the Arabidopsis thaliana K homology (KH) domain protein REGULATOR OF CBF GENE EXPRESSION 3 (RCF3) as a cofactor affecting miRNA biogenesis in specific plant tissues. MiRNA and miRNA-target levels were reduced in apex-enriched samples of rcf3 mutants, but not in other tissues. Mechanistically, RCF3 affects miRNA biogenesis through nuclear interactions with the phosphatases C-TERMINAL DOMAIN PHOSPHATASE-LIKE1 and 2 (CPL1 and CPL2). These interactions are essential to regulate the phosphorylation status, and thus the activity, of the double-stranded RNA binding protein and DICER-LIKE1 (DCL1) cofactor HYPONASTIC LEAVES1 (HYL1).

Modulation of Ambient Temperature-Dependent Flowering in Arabidopsis thaliana by Natural Variation of FLOWERING LOCUS M.

Plants integrate seasonal cues such as temperature and day length to optimally adjust their flowering time to the environment. Compared to the control of flowering before and after winter by the vernalization and day length pathways, mechanisms that delay or promote flowering during a transient cool or warm period, especially during spring, are less well understood. Due to global warming, understanding this ambient temperature pathway has gained increasing importance. In Arabidopsis thaliana, FLOWERING LOCUS M (FLM) is a critical flowering regulator of the ambient temperature pathway. FLM is alternatively spliced in a temperature-dependent manner and the two predominant splice variants, FLM-ß and FLM-δ, can repress and activate flowering in the genetic background of the A. thaliana reference accession Columbia-0. The relevance of this regulatory mechanism for the environmental adaptation across the entire range of the species is, however, unknown. Here, we identify insertion polymorphisms in the first intron of FLM as causative for accelerated flowering in many natural A. thaliana accessions, especially in cool (15°C) temperatures. We present evidence for a potential adaptive role of this structural variation and link it specifically to changes in the abundance of FLM-ß. Our results may allow predicting flowering in response to ambient temperatures in the Brassicaceae.

FYVE1 is essential for vacuole biogenesis and intracellular trafficking in Arabidopsis.

The plant vacuole is a central organelle that is involved in various biological processes throughout the plant life cycle. Elucidating the mechanism of vacuole biogenesis and maintenance is thus the basis for our understanding of these processes. Proper formation of the vacuole has been shown to depend on the intracellular membrane trafficking pathway. Although several mutants with altered vacuole morphology have been characterized in the past, the molecular basis for plant vacuole biogenesis has yet to be fully elucidated. With the aim to identify key factors that are essential for vacuole biogenesis, we performed a forward genetics screen in Arabidopsis (Arabidopsis thaliana) and isolated mutants with altered vacuole morphology. The vacuolar fusion defective1 (vfd1) mutant shows seedling lethality and defects in central vacuole formation. VFD1 encodes a Fab1, YOTB, Vac1, and EEA1 (FYVE) domain-containing protein, FYVE1, that has been implicated in intracellular trafficking. FYVE1 localizes on late endosomes and interacts with Src homology-3 domain-containing proteins. Mutants of FYVE1 are defective in ubiquitin-mediated protein degradation, vacuolar transport, and autophagy. Altogether, our results show that FYVE1 is essential for plant growth and development and place FYVE1 as a key regulator of intracellular trafficking and vacuole biogenesis.

Century-scale methylome stability in a recently diverged Arabidopsis thaliana lineage.

There has been much excitement about the possibility that exposure to specific environments can induce an ecological memory in the form of whole-sale, genome-wide epigenetic changes that are maintained over many generations. In the model plant Arabidopsis thaliana, numerous heritable DNA methylation differences have been identified in greenhouse-grown isogenic lines, but it remains unknown how natural, highly variable environments affect the rate and spectrum of such changes. Here we present detailed methylome analyses in a geographically dispersed A. thaliana population that constitutes a collection of near-isogenic lines, diverged for at least a century from a common ancestor. Methylome variation largely reflected genetic distance, and was in many aspects similar to that of lines raised in uniform conditions. Thus, even when plants are grown in varying and diverse natural sites, genome-wide epigenetic variation accumulates mostly in a clock-like manner, and epigenetic divergence thus parallels the pattern of genome-wide DNA sequence divergence.

Genome expansion of Arabis alpina linked with retrotransposition and reduced symmetric DNA methylation.

Despite evolutionary conserved mechanisms to silence transposable element activity, there are drastic differences in the abundance of transposable elements even among closely related plant species. We conducted a de novo assembly for the 375 Mb genome of the perennial model plant, Arabis alpina. Analysing this genome revealed long-lasting and recent transposable element activity predominately driven by Gypsy long terminal repeat retrotransposons, which extended the low-recombining pericentromeres and transformed large formerly euchromatic regions into repeat-rich pericentromeric regions. This reduced capacity for long terminal repeat retrotransposon silencing and removal in A. alpina co-occurs with unexpectedly low levels of DNA methylation. Most remarkably, the striking reduction of symmetrical CG and CHG methylation suggests weakened DNA methylation maintenance in A. alpina compared with Arabidopsis thaliana. Phylogenetic analyses indicate a highly dynamic evolution of some components of methylation maintenance machinery that might be related to the unique methylation in A. alpina.