FastQFS – A tool for evaluating and filtering paired-end sequencing data generated from high throughput sequencing

Cox2 community barcoding at Prince Edward Island reveals long-distance dispersal of a downy mildew species and potentially marine members of the Saprolegniaceae

Marine oomycetes are highly diverse, globally distributed, and play key roles in marine food webs as decomposers, food source, and parasites. Despite their potential importance in global ocean ecosystems, marine oomycetes are comparatively little studied. Here, we tested if the primer pair cox2F_Hud and cox2-RC4, which is already well-established for phylogenetic investigations of terrestrial oomycetes, can also be used for high-throughput community barcoding. Community barcoding of a plankton sample from Brudenell River (Prince Edward Island, Canada), revealed six distinct oomycete OTU clusters. Two of these clusters corresponded to members of the Peronosporaceae—one could be assigned to Peronospora verna, an obligate biotrophic pathogen of the terrestrial plant Veronica serpyllifolia and related species, the other was closely related to Globisporangium rostratum. While the detection of the former in the sample is likely due to long-distance dispersal from the island, the latter might be a bona fide marine species, as several cultivable species of the Peronosporaceae are known to withstand high salt concentrations. Two OTU lineages could be assigned to the Saprolegniaceae. While these might represent marine species of the otherwise terrestrial genus, it is also conceivable that they were introduced on detritus from the island. Two additional OTU clusters were grouped with the early-diverging oomycete lineages but could not be assigned to a specific family. This reflects the current underrepresentation of cox2 sequence data which will hopefully improve with the increasing interest in marine oomycetes.

The Genome of Microthlaspi erraticum (Brassicaceae) Provides Insights Into the Adaptation to Highly Calcareous Soils

Microthlaspi erraticum is widely distributed in temperate Eurasia, but restricted to Ca²⁺-rich habitats, predominantly on white Jurassic limestone, which is made up by calcium carbonate, with little other minerals. Thus, naturally occurring Microthlaspi erraticum individuals are confronted with a high concentration of Ca²⁺ ions while Mg²⁺ ion concentration is relatively low. As there is a competitive uptake between these two ions, adaptation to the soil condition can be expected. In this study, it was the aim to explore the genomic consequences of this adaptation by sequencing and analysing the genome of Microthlaspi erraticum. Its genome size is comparable with other diploid Brassicaceae, while more genes were predicted. Two Mg²⁺ transporters known to be expressed in roots were duplicated and one showed a significant degree of positive selection. It is speculated that this evolved due to the pressure to take up Mg²⁺ ions efficiently in the presence of an overwhelming amount of Ca²⁺ ions. Future studies on plants specialized on similar soils and affinity tests of the transporters are needed to provide unequivocal evidence for this hypothesis. If verified, the transporters found in this study might be useful for breeding Brassicaceae crops for higher yield on Ca²⁺-rich and Mg²⁺ -poor soils.

The Genome of Peronospora belbahrii Reveals High Heterozygosity, a Low Number of Canonical Effectors, and TC-Rich Promoters

Along with Plasmopara destructor, Peronosopora belbahrii has arguably been the economically most important newly emerging downy mildew pathogen of the past two decades. Originating from Africa, it has started devastating basil production throughout the world, most likely due to the distribution of infested seed material. Here, we present the genome of this pathogen and results from comparisons of its genomic features to other oomycetes. The assembly of the nuclear genome was around 35.4 Mbp in length, with an N₅₀ scaffold length of around 248 kbp and an L₅₀ scaffold count of 46. The circular mitochondrial genome consisted of around 40.1 kbp. From the repeat-masked genome, 9,049 protein-coding genes were predicted, out of which 335 were predicted to have extracellular functions, representing the smallest secretome so far found in peronosporalean oomycetes. About 16% of the genome consists of repetitive sequences, and, based on simple sequence repeat regions, we provide a set of microsatellites that could be used for population genetic studies of P. belbahrii. P. belbahrii has undergone a high degree of convergent evolution with other obligate parasitic pathogen groups, reflecting its obligate biotrophic lifestyle. Features of its secretome, signaling networks, and promoters are presented, and some patterns are hypothesized to reflect the high degree of host specificity in Peronospora species. In addition, we suggest the presence of additional virulence factors apart from classical effector classes that are promising candidates for future functional studies.

Community barcoding reveals little effect of ocean acidification on the composition of coastal plankton communities: Evidence from a long-term mesocosm study in the Gullmar Fjord, Skagerrak

The acidification of the oceans could potentially alter marine plankton communities with consequences for ecosystem functioning. While several studies have investigated effects of ocean acidification on communities using traditional methods, few have used genetic analyses. Here, we use community barcoding to assess the impact of ocean acidification on the composition of a coastal plankton community in a large scale, in situ, long-term mesocosm experiment. High-throughput sequencing resulted in the identification of a wide range of planktonic taxa (Alveolata, Cryptophyta, Haptophyceae, Fungi, Metazoa, Hydrozoa, Rhizaria, Straminipila, Chlorophyta). Analyses based on predicted operational taxonomical units as well as taxonomical compositions revealed no differences between communities in high CO2 mesocosms (~ 760 μatm) and those exposed to present-day CO2 conditions. Observed shifts in the planktonic community composition were mainly related to seasonal changes in temperature and nutrients. Furthermore, based on our investigations, the elevated CO2 did not affect the intraspecific diversity of the most common mesozooplankter, the calanoid copepod Pseudocalanus acuspes. Nevertheless, accompanying studies found temporary effects attributed to a raise in CO2. Differences in taxa composition between the CO2 treatments could, however, only be observed in a specific period of the experiment. Based on our genetic investigations, no compositional long-term shifts of the plankton communities exposed to elevated CO2 conditions were observed. Thus, we conclude that the compositions of planktonic communities, especially those in coastal areas, remain rather unaffected by increased CO2.

Adaptive differentiation coincides with local bioclimatic conditions along an elevational cline in populations of a lichen-forming fungus

Background

Many fungal species occur across a variety of habitats. Particularly lichens, fungi forming symbioses with photosynthetic partners, have evolved remarkable tolerances for environmental extremes. Despite their ecological importance and ubiquity, little is known about the genetic basis of adaption in lichen populations. Here we studied patterns of genome-wide differentiation in the lichen-forming fungus Lasallia pustulata along an altitudinal gradient in the Mediterranean region. We resequenced six populations as pools and identified highly differentiated genomic regions. We then detected gene-environment correlations while controlling for shared population history and pooled sequencing bias, and performed ecophysiological experiments to assess fitness differences of individuals from different environments.

Results

We detected two strongly differentiated genetic clusters linked to Mediterranean and temperate-oceanic climate, and an admixture zone, which coincided with the transition between the two bioclimates. High altitude individuals showed ecophysiological adaptations to wetter and more shaded conditions. Highly differentiated genome regions contained a number of genes associated with stress response, local environmental adaptation, and sexual reproduction.

Conclusions

Taken together our results provide evidence for a complex interplay between demographic history and spatially varying selection acting on a number of key biological processes, suggesting a scenario of ecological speciation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-0929-8) contains supplementary material, which is available to authorized users.

Genome analyses of the sunflower pathogen Plasmopara halstedii provide insights into effector evolution in downy mildews and Phytophthora

BACKGROUND

Downy mildews are the most speciose group of oomycetes and affect crops of great economic importance. So far, there is only a single deeply-sequenced downy mildew genome available, from Hyaloperonospora arabidopsidis. Further genomic resources for downy mildews are required to study their evolution, including pathogenicity effector proteins, such as RxLR effectors. Plasmopara halstedii is a devastating pathogen of sunflower and a potential pathosystem model to study downy mildews, as several Avr-genes and R-genes have been predicted and unlike Arabidopsis downy mildew, large quantities of almost contamination-free material can be obtained easily.

RESULTS

Here a high-quality draft genome of Plasmopara halstedii is reported and analysed with respect to various aspects, including genome organisation, secondary metabolism, effector proteins and comparative genomics with other sequenced oomycetes. Interestingly, the present analyses revealed further variation of the RxLR motif, suggesting an important role of the conservation of the dEER-motif. Orthology analyses revealed the conservation of 28 RxLR-like core effectors among Phytophthora species. Only six putative RxLR-like effectors were shared by the two sequenced downy mildews, highlighting the fast and largely independent evolution of two of the three major downy mildew lineages. This is seemingly supported by phylogenomic results, in which downy mildews did not appear to be monophyletic.

CONCLUSIONS

The genome resource will be useful for developing markers for monitoring the pathogen population and might provide the basis for new approaches to fight Phytophthora and downy mildew pathogens by targeting core pathogenicity effectors.

Comparative Genomics Including the Early-Diverging Smut Fungus Ceraceosorus bombacis Reveals Signatures of Parallel Evolution within Plant and Animal Pathogens of Fungi and Oomycetes

Ceraceosorus bombacis is an early-diverging lineage of smut fungi and a pathogen of cotton trees (Bombax ceiba). To study the evolutionary genomics of smut fungi in comparison with other fungal and oomycete pathogens, the genome of C. bombacis was sequenced and comparative genomic analyses were performed. The genome of 26.09 Mb encodes for 8,024 proteins, of which 576 are putative-secreted effector proteins (PSEPs). Orthology analysis revealed 30 ortholog PSEPs among six Ustilaginomycotina genomes, the largest groups of which are lytic enzymes, such as aspartic peptidase and glycoside hydrolase. Positive selection analyses revealed the highest percentage of positively selected PSEPs in C. bombacis compared with other Ustilaginomycotina genomes. Metabolic pathway analyses revealed the absence of genes encoding for nitrite and nitrate reductase in the genome of the human skin pathogen Malassezia globosa, but these enzymes are present in the sequenced plant pathogens in smut fungi. Interestingly, these genes are also absent in cultivable oomycete animal pathogens, while nitrate reductase has been lost in cultivable oomycete plant pathogens. Similar patterns were also observed for obligate biotrophic and hemi-biotrophic fungal and oomycete pathogens. Furthermore, it was found that both fungal and oomycete animal pathogen genomes are lacking cutinases and pectinesterases. Overall, these findings highlight the parallel evolution of certain genomic traits, revealing potential common evolutionary trajectories among fungal and oomycete pathogens, shaping the pathogen genomes according to their lifestyle.