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.

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.

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.

Multiple modes of convergent adaptation in the spread of glyphosate-resistant Amaranthus tuberculatus.

The selection pressure exerted by herbicides has led to the repeated evolution of herbicide resistance in weeds. The evolution of herbicide resistance on contemporary timescales in turn provides an outstanding opportunity to investigate key questions about the genetics of adaptation, in particular the relative importance of adaptation from new mutations, standing genetic variation, or geographic spread of adaptive alleles through gene flow. Glyphosate-resistant Amaranthus tuberculatus poses one of the most significant threats to crop yields in the Midwestern United States, with both agricultural populations and herbicide resistance only recently emerging in Canada. To understand the evolutionary mechanisms driving the spread of resistance, we sequenced and assembled the A. tuberculatus genome and investigated the origins and population genomics of 163 resequenced glyphosate-resistant and susceptible individuals from Canada and the United States. In Canada, we discovered multiple modes of convergent evolution: in one locality, resistance appears to have evolved through introductions of preadapted US genotypes, while in another, there is evidence for the independent evolution of resistance on genomic backgrounds that are historically nonagricultural. Moreover, resistance on these local, nonagricultural backgrounds appears to have occurred predominantly through the partial sweep of a single haplotype. In contrast, resistant haplotypes arising from the Midwestern United States show multiple amplification haplotypes segregating both between and within populations. Therefore, while the remarkable species-wide diversity of A. tuberculatus has facilitated geographic parallel adaptation of glyphosate resistance, more recently established agricultural populations are limited to adaptation in a more mutation-limited framework.

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.

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.

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.

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.

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.

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.

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.

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.

Different but overlapping populations of Strongyloides stercoralis in dogs and humans-Dogs as a possible source for zoonotic strongyloidiasis.

Strongyloidiasis is a much-neglected soil born helminthiasis caused by the nematode Strongyloides stercoralis. Human derived S. stercoralis can be maintained in dogs in the laboratory and this parasite has been reported to also occur in dogs in the wild. Some authors have considered strongyloidiasis a zoonotic disease while others have argued that the two hosts carry host specialized populations of S. stercoralis and that dogs play a minor role, if any, as a reservoir for zoonotic S. stercoralis infections of humans. We isolated S. stercoralis from humans and their dogs in rural villages in northern Cambodia, a region with a high incidence of strongyloidiasis, and compared the worms derived from these two host species using nuclear and mitochondrial DNA sequence polymorphisms. We found that in dogs there exist two populations of S. stercoralis, which are clearly separated from each other genetically based on the nuclear 18S rDNA, the mitochondrial cox1 locus and whole genome sequence. One population, to which the majority of the worms belong, appears to be restricted to dogs. The other population is indistinguishable from the population of S. stercoralis isolated from humans. Consistent with earlier studies, we found multiple sequence variants of the hypervariable region I of the 18 S rDNA in S. stercoralis from humans. However, comparison of mitochondrial sequences and whole genome analysis suggest that these different 18S variants do not represent multiple genetically isolated subpopulations among the worms isolated from humans. We also investigated the mode of reproduction of the free-living generations of laboratory and wild isolates of S. stercoralis. Contrary to earlier literature on S. stercoralis but similar to other species of Strongyloides, we found clear evidence of sexual reproduction. Overall, our results show that dogs carry two populations, possibly different species of Strongyloides. One population appears to be dog specific but the other one is shared with humans. This argues for the strong potential of dogs as reservoirs for zoonotic transmission of S. stercoralis to humans and suggests that in order to reduce the exposure of humans to infective S. stercoralis larvae, dogs should be treated for the infection along with their owners.

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.

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.

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.

Calcium sensor kinase activates potassium uptake systems in gland cells of Venus flytraps.

The Darwin plant Dionaea muscipula is able to grow on mineral-poor soil, because it gains essential nutrients from captured animal prey. Given that no nutrients remain in the trap when it opens after the consumption of an animal meal, we here asked the question of how Dionaea sequesters prey-derived potassium. We show that prey capture triggers expression of a K(+) uptake system in the Venus flytrap. In search of K(+) transporters endowed with adequate properties for this role, we screened a Dionaea expressed sequence tag (EST) database and identified DmKT1 and DmHAK5 as candidates. On insect and touch hormone stimulation, the number of transcripts of these transporters increased in flytraps. After cRNA injection of K(+)-transporter genes into Xenopus oocytes, however, both putative K(+) transporters remained silent. Assuming that calcium sensor kinases are regulating Arabidopsis K(+) transporter 1 (AKT1), we coexpressed the putative K(+) transporters with a large set of kinases and identified the CBL9-CIPK23 pair as the major activating complex for both transporters in Dionaea K(+) uptake. DmKT1 was found to be a K(+)-selective channel of voltage-dependent high capacity and low affinity, whereas DmHAK5 was identified as the first, to our knowledge, proton-driven, high-affinity potassium transporter with weak selectivity. When the Venus flytrap is processing its prey, the gland cell membrane potential is maintained around -120 mV, and the apoplast is acidified to pH 3. These conditions in the green stomach formed by the closed flytrap allow DmKT1 and DmHAK5 to acquire prey-derived K(+), reducing its concentration from millimolar levels down to trace levels.

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 NO3(-), NH4(+) 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 NH4(+) 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 NH4(+) was mediated by 2.5-fold higher expression of the NH4(+) transporter DmAMT1 in the roots compared with the traps. Electrophysiological studies confirmed a high constitutive capacity for NH4(+) 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.