Deleterious phenotypes in wild Arabidopsis arenosa populations are common and linked to runs of homozygosity

In this study, we aimed to systematically assess the frequency at which potentially deleterious phenotypes appear in natural populations of the outcrossing model plant Arabidopsis arenosa, and to establish their underlying genetics. For this purpose, we collected seeds from wild A. arenosa populations and screened over 2,500 plants for unusual phenotypes in the greenhouse. We repeatedly found plants with obvious phenotypic defects, such as small stature and necrotic or chlorotic leaves, among first-generation progeny of wild A. arenosa plants. Such abnormal plants were present in about 10% of maternal sibships, with multiple plants with similar phenotypes in each of these sibships, pointing to a genetic basis of the observed defects. A combination of transcriptome profiling, linkage mapping and genome-wide runs of homozygosity patterns using a newly assembled reference genome indicated a range of underlying genetic architectures associated with phenotypic abnormalities. This included evidence for homozygosity of certain genomic regions, consistent with alleles that are identical by descent being responsible for these defects. Our observations suggest that deleterious alleles with different genetic architectures are segregating at appreciable frequencies in wild A. arenosa populations.

GWAS reveals a rapidly evolving candidate avirulence effector in the Cercospora leaf spot pathogen

The major resistance gene BvCR4 recently bred into sugar beet hybrids provides a high level of resistance to Cercospora leaf spot caused by the fungal pathogen Cercospora beticola. The occurrence of pathogen strains that overcome BvCR4 was studied using field trials in Switzerland conducted under natural disease pressure. Virulence of a subset of these strains was evaluated in a field trial conducted under elevated artificial disease pressure. We created a new C. beticola reference genome and mapped whole genome sequences of 256 isolates collected in Switzerland and Germany. These were combined with virulence phenotypes to conduct three separate genome-wide association studies (GWAS) to identify candidate avirulence genes. We identified a locus associated with avirulence containing a putative avirulence effector gene named AvrCR4. All virulent isolates either lacked AvrCR4 or had nonsynonymous mutations within the gene. AvrCR4 was present in all 74 isolates from non-BvCR4 hybrids, whereas 33 of 89 isolates from BvCR4 hybrids carried a deletion. We also mapped genomic data from 190 publicly available US isolates to our new reference genome. The AvrCR4 deletion was found in only one of 95 unique isolates from non-BvCR4 hybrids in the United States. AvrCR4 presents a unique example of an avirulence effector in which virulent alleles have only recently emerged. Most likely these were selected out of standing genetic variation after deployment of BvCR4. Identification of AvrCR4 will enable real-time screening of C. beticola populations for the emergence and spread of virulent isolates.

Improved Reference Genome Uncovers Novel Sex-Linked Regions in the Guppy (Poecilia reticulata)

Theory predicts that the sexes can achieve greater fitness if loci with sexually antagonistic polymorphisms become linked to the sex determining loci, and this can favor the spread of reduced recombination around sex determining regions. Given that sex-linked regions are frequently repetitive and highly heterozygous, few complete Y chromosome assemblies are available to test these ideas. The guppy system (Poecilia reticulata) has long been invoked as an example of sex chromosome formation resulting from sexual conflict. Early genetics studies revealed that male color patterning genes are mostly but not entirely Y-linked, and that X-linkage may be most common in low-predation populations. More recent population genomic studies of guppies have reached varying conclusions about the size and placement of the Y-linked region. However, this previous work used a reference genome assembled from short-read sequences from a female guppy. Here, we present a new guppy reference genome assembly from a male, using long-read PacBio single-molecule real-time sequencing and chromosome contact information. Our new assembly sequences across repeat- and GC-rich regions and thus closes gaps and corrects mis-assemblies found in the short-read female-derived guppy genome. Using this improved reference genome, we then employed broad population sampling to detect sex differences across the genome. We identified two small regions that showed consistent male-specific signals. Moreover, our results help reconcile the contradictory conclusions put forth by past population genomic studies of the guppy sex chromosome. Our results are consistent with a small Y-specific region and rare recombination in male guppies.

A Species-Wide Inventory of NLR Genes and Alleles in Arabidopsis thaliana

Infectious disease is both a major force of selection in nature and a prime cause of yield loss in agriculture. In plants, disease resistance is often conferred by nucleotide-binding leucine-rich repeat (NLR) proteins, intracellular immune receptors that recognize pathogen proteins and their effects on the host. Consistent with extensive balancing and positive selection, NLRs are encoded by one of the most variable gene families in plants, but the true extent of intraspecific NLR diversity has been unclear. Here, we define a nearly complete species-wide pan-NLRome in Arabidopsis thaliana based on sequence enrichment and long-read sequencing. The pan-NLRome largely saturates with approximately 40 well-chosen wild strains, with half of the pan-NLRome being present in most accessions. We chart NLR architectural diversity, identify new architectures, and quantify selective forces that act on specific NLRs and NLR domains. Our study provides a blueprint for defining pan-NLRomes.

•

Species-wide NLR diversity is high but not unlimited

•

A large fraction of NLR diversity is recovered with 40–50 accessions

•

Presence/absence variation in NLRs is widespread, resulting in a mosaic population

•

A high diversity of NLR-integrated domains favor known virulence targets

A Species-Wide Inventory of NLR Genes and Alleles in Arabidopsis thaliana

Infectious disease is both a major force of selection in nature and a prime cause of yield loss in agriculture. In plants, disease resistance is often conferred by nucleotide-binding leucine-rich repeat (NLR) proteins, intracellular immune receptors that recognize pathogen proteins and their effects on the host. Consistent with extensive balancing and positive selection, NLRs are encoded by one of the most variable gene families in plants, but the true extent of intraspecific NLR diversity has been unclear. Here, we define a nearly complete species-wide pan-NLRome in Arabidopsis thaliana based on sequence enrichment and long-read sequencing. The pan-NLRome largely saturates with approximately 40 well-chosen wild strains, with half of the pan-NLRome being present in most accessions. We chart NLR architectural diversity, identify new architectures, and quantify selective forces that act on specific NLRs and NLR domains. Our study provides a blueprint for defining pan-NLRomes.

Complex Relationships between Chromatin Accessibility, Sequence Divergence, and Gene Expression in Arabidopsis thaliana

Variation in regulatory DNA is thought to drive phenotypic variation, evolution, and disease. Prior studies of regulatory DNA and transcription factors across animal species highlighted a fundamental conundrum: Transcription factor binding domains and cognate binding sites are conserved, while regulatory DNA sequences are not. It remains unclear how conserved transcription factors and dynamic regulatory sites produce conserved expression patterns across species. Here, we explore regulatory DNA variation and its functional consequences within Arabidopsis thaliana, using chromatin accessibility to delineate regulatory DNA genome-wide. Unlike in previous cross-species comparisons, the positional homology of regulatory DNA is maintained among A. thaliana ecotypes and less nucleotide divergence has occurred. Of the ∼50,000 regulatory sites in A. thaliana, we found that 15% varied in accessibility among ecotypes. Some of these accessibility differences were associated with extensive, previously unannotated sequence variation, encompassing many deletions and ancient hypervariable alleles. Unexpectedly, for the majority of such regulatory sites, nearby gene expression was unaffected. Nevertheless, regulatory sites with high levels of sequence variation and differential chromatin accessibility were the most likely to be associated with differential gene expression. Finally, and most surprising, we found that the vast majority of differentially accessible sites show no underlying sequence variation. We argue that these surprising results highlight the necessity to consider higher-order regulatory context in evaluating regulatory variation and predicting its phenotypic consequences.

Venus flytrap microbiotas withstand harsh conditions during prey digestion

The carnivorous Venus flytrap (Dionaea muscipula) overcomes environmental nutrient limitation by capturing small animals. Such prey is digested with an acidic enzyme-containing mucilage that is secreted into the closed trap. However, surprisingly little is known about associations with microorganisms. Therefore, we assessed microbiotas of traps and petioles for the Venus flytrap by 16S amplicon meta-barcoding. We also performed time-series assessments of dynamics during digestion in traps and experimental acidification of petioles. We found that the traps hosted distinct microbiotas that differed from adjacent petioles. Further, they showed a significant taxonomic turnover during digestion. Following successful catches, prey-associated bacteria had strong effects on overall composition. With proceeding digestion, however, microbiotas were restored to compositions resembling pre-digestion stages. A comparable, yet less extensive shift was found when stimulating digestion with coronatine. Artificial acidification of petioles did not induce changes towards trap-like communities. Our results show that trap microbiota were maintained during digestion despite harsh conditions and recovered after short-term disturbances through prey microbiota. This indicates trap-specific and resilient associations. By mapping to known genomes, we predicted putative adaptations and functional implications for the system, yet direct mechanisms and quantification of host benefits, like the involvement in digestion, remain to be addressed.

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.

Modulation of ACD6 dependent hyperimmunity by natural alleles of an Arabidopsis thaliana NLR resistance gene

Plants defend themselves against pathogens by activating an array of immune responses. Unfortunately, immunity programs may also cause unintended collateral damage to the plant itself. The quantitative disease resistance gene ACCELERATED CELL DEATH 6 (ACD6) serves to balance growth and pathogen resistance in natural populations of Arabidopsis thaliana. An autoimmune allele, ACD6-Est, which strongly reduces growth under specific laboratory conditions, is found in over 10% of wild strains. There is, however, extensive variation in the strength of the autoimmune phenotype expressed by strains with an ACD6-Est allele, indicative of genetic modifiers. Quantitative genetic analysis suggests that ACD6 activity can be modulated in diverse ways, with different strains often carrying different large-effect modifiers. One modifier is SUPPRESSOR OF NPR1-1, CONSTITUTIVE 1 (SNC1), located in a highly polymorphic cluster of nucleotide-binding domain and leucine-rich repeat (NLR) immune receptor genes, which are prototypes for qualitative disease resistance genes. Allelic variation at SNC1 correlates with ACD6-Est activity in multiple accessions, and a common structural variant affecting the NL linker sequence can explain differences in SNC1 activity. Taken together, we find that an NLR gene can mask the activity of an ACD6 autoimmune allele in natural A. thaliana populations, thereby linking different arms of the plant immune system.

Plants defend themselves against pathogens by activating immune responses. Unfortunately, these can cause unintended collateral damage to the plant itself. Nevertheless, some wild plants have genetic variants that confer a low threshold for the activation of immunity. While these enable a plant to respond particularly quickly to pathogen attack, such variants might be potentially dangerous. We are investigating one such variant of the immune gene ACCELERATED CELL DEATH 6 (ACD6) in the plant Arabidopsis thaliana. We discovered that there are variants at other genetic loci that can mask the effects of an overly active ACD6 gene. One of these genes, SUPPRESSOR OF NPR1-1, CONSTITUTIVE 1 (SNC1), codes for a known immune receptor. The SNC1 variant that attenuates ACD6 activity is rather common in A. thaliana populations, suggesting that new combinations of the hyperactive ACD6 variant and this antagonistic SNC1 variant will often arise by natural crosses. Similarly, because the two genes are unlinked, outcrossing will often lead to the hyperactive ACD6 variants being unmasked again. We propose that allelic diversity at SNC1 contributes to the maintenance of the hyperactive ACD6 variant in natural A. thaliana populations.

High V-PPase activity is beneficial under high salt loads, but detrimental without salinity

The membrane-bound proton-pumping pyrophosphatase (V-PPase), together with the V-type H+ -ATPase, generates the proton motive force that drives vacuolar membrane solute transport. Transgenic plants constitutively overexpressing V-PPases were shown to have improved salinity tolerance, but the relative impact of increasing PPi hydrolysis and proton-pumping functions has yet to be dissected. For a better understanding of the molecular processes underlying V-PPase-dependent salt tolerance, we transiently overexpressed the pyrophosphate-driven proton pump (NbVHP) in Nicotiana benthamiana leaves and studied its functional properties in relation to salt treatment by primarily using patch-clamp, impalement electrodes and pH imaging. NbVHP overexpression led to higher vacuolar proton currents and vacuolar acidification. After 3 d in salt-untreated conditions, V-PPase-overexpressing leaves showed a drop in photosynthetic capacity, plasma membrane depolarization and eventual leaf necrosis. Salt, however, rescued NbVHP-hyperactive cells from cell death. Furthermore, a salt-induced rise in V-PPase but not of V-ATPase pump currents was detected in nontransformed plants. The results indicate that under normal growth conditions, plants need to regulate the V-PPase pump activity to avoid hyperactivity and its negative feedback on cell viability. Nonetheless, V-PPase proton pump function becomes increasingly important under salt stress for generating the pH gradient necessary for vacuolar proton-coupled Na+ sequestration.

Single-Molecule Sequencing Reveals the Chromosome-Scale Genomic Architecture of the Nematode Model Organism Pristionchus pacificus

The nematode Pristionchus pacificus is an established model for integrative evolutionary biology and comparative studies with Caenorhabditis elegans. While an existing genome draft facilitated the identification of several genes controlling various developmental processes, its high degree of fragmentation complicated virtually all genomic analyses. Here, we present a de novo genome assembly from single-molecule, long-read sequencing data consisting of 135 P. pacificus contigs. When combined with a genetic linkage map, 99% of the assembly could be ordered and oriented into six chromosomes. This allowed us to robustly characterize chromosomal patterns of gene density, repeat content, nucleotide diversity, linkage disequilibrium, and macrosynteny in P. pacificus. Despite widespread conservation of synteny between P. pacificus and C. elegans, we identified one major translocation from an autosome to the sex chromosome in the lineage leading to C. elegans. This highlights the potential of the chromosome-scale assembly for future genomic studies of P. pacificus.

High contiguity Arabidopsis thaliana genome assembly with a single nanopore flow cell

The handheld Oxford Nanopore MinION sequencer generates ultra-long reads with minimal cost and time requirements, which makes sequencing genomes at the bench feasible. Here, we sequence the gold standard Arabidopsis thaliana genome (KBS-Mac-74 accession) on the bench with the MinION sequencer, and assemble the genome using typical consumer computing hardware (4 Cores, 16 Gb RAM) into chromosome arms (62 contigs with an N50 length of 12.3 Mb). We validate the contiguity and quality of the assembly with two independent single-molecule technologies, Bionano optical genome maps and Pacific Biosciences Sequel sequencing. The new A. thaliana KBS-Mac-74 genome enables resolution of a quantitative trait locus that had previously been recalcitrant to a Sanger-based BAC sequencing approach. In summary, we demonstrate that even when the purpose is to understand complex structural variation at a single region of the genome, complete genome assembly is becoming the simplest way to achieve this goal.

Epigenomic Diversity in a Global Collection of Arabidopsis thaliana Accessions

The epigenome orchestrates genome accessibility, functionality, and three-dimensional structure. Because epigenetic variation can impact transcription and thus phenotypes, it may contribute to adaptation. Here, we report 1,107 high-quality single-base resolution methylomes and 1,203 transcriptomes from the 1001 Genomes collection of Arabidopsis thaliana. Although the genetic basis of methylation variation is highly complex, geographic origin is a major predictor of genome-wide DNA methylation levels and of altered gene expression caused by epialleles. Comparison to cistrome and epicistrome datasets identifies associations between transcription factor binding sites, methylation, nucleotide variation, and co-expression modules. Physical maps for nine of the most diverse genomes reveal how transposons and other structural variants shape the epigenome, with dramatic effects on immunity genes. The 1001 Epigenomes Project provides a comprehensive resource for understanding how variation in DNA methylation contributes to molecular and non-molecular phenotypes in natural populations of the most studied model plant.

TBro: visualization and management of de novo transcriptomes

RNA sequencing (RNA-seq) has become a powerful tool to understand molecular mechanisms and/or developmental programs. It provides a fast, reliable and cost-effective method to access sets of expressed elements in a qualitative and quantitative manner. Especially for non-model organisms and in absence of a reference genome, RNA-seq data is used to reconstruct and quantify transcriptomes at the same time. Even SNPs, InDels, and alternative splicing events are predicted directly from the data without having a reference genome at hand. A key challenge, especially for non-computational personnal, is the management of the resulting datasets, consisting of different data types and formats. Here, we present TBro, a flexible de novo transcriptome browser, tackling this challenge. TBro aggregates sequences, their annotation, expression levels as well as differential testing results. It provides an easy-to-use interface to mine the aggregated data and generate publication-ready visualizations. Additionally, it supports users with an intuitive cart system, that helps collecting and analysing biological meaningful sets of transcripts. TBro’s modular architecture allows easy extension of its functionalities in the future. Especially, the integration of new data types such as proteomic quantifications or array-based gene expression data is straightforward. Thus, TBro is a fully featured yet flexible transcriptome browser that supports approaching complex biological questions and enhances collaboration of numerous researchers.

Database URL: tbro.carnivorom.com

Venus flytrap carnivorous lifestyle builds on herbivore defense strategies

Although the concept of botanical carnivory has been known since Darwin's time, the molecular mechanisms that allow animal feeding remain unknown, primarily due to a complete lack of genomic information. Here, we show that the transcriptomic landscape of the Dionaea trap is dramatically shifted toward signal transduction and nutrient transport upon insect feeding, with touch hormone signaling and protein secretion prevailing. At the same time, a massive induction of general defense responses is accompanied by the repression of cell death-related genes/processes. We hypothesize that the carnivory syndrome of Dionaea evolved by exaptation of ancient defense pathways, replacing cell death with nutrient acquisition.

Integration of trap- and root-derived nitrogen nutrition of carnivorous Dionaea muscipula

Carnivorous Dionaea muscipula operates active snap traps for nutrient acquisition from prey; so what is the role of D. muscipula's reduced root system? We studied the capacity for nitrogen (N) acquisition via traps, and its effect on plant allometry; the capacity of roots to absorb NO₃(-), NH₄(+) and glutamine from the soil solution; and the fate and interaction of foliar- and root-acquired N. Feeding D. muscipula snap traps with insects had little effect on the root : shoot ratio, but promoted petiole relative to trap growth. Large amounts of NH₄(+) and glutamine were absorbed upon root feeding. The high capacity for root N uptake was maintained upon feeding traps with glutamine. High root acquisition of NH₄(+) was mediated by 2.5-fold higher expression of the NH₄(+) transporter DmAMT1 in the roots compared with the traps. Electrophysiological studies confirmed a high constitutive capacity for NH₄(+) uptake by roots. Glutamine feeding of traps inhibited the influx of (15)N from root-absorbed (15)N/(13)C-glutamine into these traps, but not that of (13)C. Apparently, fed traps turned into carbon sinks that even acquired organic carbon from roots. N acquisition at the whole-plant level is fundamentally different in D. muscipula compared with noncarnivorous species, where foliar N influx down-regulates N uptake by roots.

Family 34 glycosyltransferase (GT34) genes and proteins in Pinus radiata (radiata pine) and Pinus taeda (loblolly pine)

Using a functional genomics approach, four candidate genes (PtGT34A, PtGT34B, PtGT34C and PtGT34D) were identified in Pinus taeda. These genes encode CAZy family GT34 glycosyltransferases that are involved in the synthesis of cell-wall xyloglucans and heteromannans. The full-length coding sequences of three orthologs (PrGT34A, B and C) were isolated from a xylem-specific cDNA library from the closely related Pinus radiata. PrGT34B is the ortholog of XXT1 and XXT2, the two main xyloglucan (1→6)-α-xylosyltransferases in Arabidopsis thaliana. PrGT34C is the ortholog of XXT5 in A. thaliana, which is also involved in the xylosylation of xyloglucans. PrGT34A is an ortholog of a galactosyltransferase from fenugreek (Trigonella foenum-graecum) that is involved in galactomannan synthesis. Truncated coding sequences of the genes were cloned into plasmid vectors and expressed in a Sf9 insect cell-culture system. The heterologous proteins were purified, and in vitro assays showed that, when incubated with UDP-xylose and cellotetraose, cellopentaose or cellohexaose, PrGT34B showed xylosyltransferase activity, and, when incubated with UDP-galactose and the same cello-oligosaccharides, PrGT34B showed some galactosyltransferase activity. The ratio of xylosyltransferase to galactosyltransferase activity was 434:1. Hydrolysis of the galactosyltransferase reaction products using galactosidases showed the linkages formed were α-linkages. Analysis of the products of PrGT34B by MALDI-TOF MS showed that up to three xylosyl residues were transferred from UDP-xylose to cellohexaose. The heterologous proteins PrGT34A and PrGT34C showed no detectable enzymatic activity.

The protein composition of the digestive fluid from the venus flytrap sheds light on prey digestion mechanisms

The Venus flytrap (Dionaea muscipula) is one of the most well-known carnivorous plants because of its unique ability to capture small animals, usually insects or spiders, through a unique snap-trapping mechanism. The animals are subsequently killed and digested so that the plants can assimilate nutrients, as they grow in mineral-deficient soils. We deep sequenced the cDNA from Dionaea traps to obtain transcript libraries, which were used in the mass spectrometry-based identification of the proteins secreted during digestion. The identified proteins consisted of peroxidases, nucleases, phosphatases, phospholipases, a glucanase, chitinases, and proteolytic enzymes, including four cysteine proteases, two aspartic proteases, and a serine carboxypeptidase. The majority of the most abundant proteins were categorized as pathogenesis-related proteins, suggesting that the plant's digestive system evolved from defense-related processes. This in-depth characterization of a highly specialized secreted fluid from a carnivorous plant provides new information about the plant's prey digestion mechanism and the evolutionary processes driving its defense pathways and nutrient acquisition.

A kinome of 2600 in the ciliate Paramecium tetraurelia

Protein kinases play a crucial role in the regulation of cellular processes. Most eukaryotes reserve about 2.5% of their genes for protein kinases. We analysed the genome of the single-celled ciliate Paramecium tetraurelia and identified 2606 kinases, about 6.6% of its genes, representing the largest kinome to date. A gene tree combined with human kinases revealed a massive expansion of the calcium calmodulin regulated subfamily, underlining the importance of calcium in the physiology of P. tetraurelia. The kinases are embedded in only 40 domain architectures, contrasting 134 in human. This might indicate different mechanisms to achieve target specificity.