Search

In the present work, two complete genome sequences of SARS-CoV-2 were obtained from nasal swab samples of Tunisian SARS-CoV-2 PCR-positive patients using nanopore sequencing. The virus genomes of two of the patients examined, a Tunisian soldier returning from a mission in Morocco and a member of another Tunisian family, showed significant differences in analyses of the total genome and single nucleotide polymorphisms (SNPs). Phylogenetic relationships with known SARS-CoV-2 genomes in the African region, some European and Middle Eastern countries and initial epidemiological conclusions indicate that the introduction of SARS-CoV-2 into Tunisia from two independent sources was travel-related.

The anthrax pathogen Bacillus anthracis poses a significant threat to human health. Identification of B. anthracis is challenging because of the bacterium’s close genetic relationship to other Bacillus cereus group species. Thus, molecular detection is founded on species-specific PCR targeting single-copy genes. Here, we validated a previously recognized multi-copy target, a species-specific single nucleotide polymorphism (SNP) present in 2–5 copies in every B. anthracis genome analyzed. For this, a hydrolysis probe-based real-time PCR assay was developed and rigorously tested. The assay was specific as only B. anthracis DNA yielded positive results, was linear over 9 log10 units, and was sensitive with a limit of detection (LoD) of 2.9 copies/reaction. Though not exhibiting a lower LoD than established single-copy PCR targets (dhp61 or PL3), the higher copy number of the B. anthracis–specific 16S rRNA gene alleles afforded ≤2 unit lower threshold (Ct) values. To push the detection limit even further, the assay was adapted for reverse transcription PCR on 16S rRNA transcripts. This RT-PCR assay was also linear over 9 log10 units and was sensitive with an LoD of 6.3 copies/reaction. In a dilution series of experiments, the 16S RT-PCR assay achieved a thousand-fold higher sensitivity than the DNA-targeting assays. For molecular diagnostics, we recommend a real-time RT-PCR assay variant in which both DNA and RNA serve as templates (thus, no requirement for DNase treatment). This can at least provide results equaling the DNA-based implementation if no RNA is present but is superior even at the lowest residual rRNA concentrations.

The anthrax pathogen Bacillus anthracis poses a significant threat to human health. Identification of B. anthracis is challenging because of the bacterium’s close genetic relationship to other Bacillus cereus group species. Thus, molecular detection is founded on species-specific PCR targeting single-copy genes. Here, we validated a previously recognized multi-copy target, a species-specific SNP present in 2-5 copies in every B. anthracis genome analyzed. For this, a hydrolysis probe-based real time PCR assay was developed and rigorously tested. The assay was specific as only B. anthracis DNA yielded positive results, was linear over 9 log10 units and was sensitive with a limit of detection (LoD) of 2.9 copies/reaction. Though not exhibiting a lower LoD than established single copy PCR targets (dhp61 or PL3), the higher copy number of the B. anthracis–specific 16S rRNA gene allele afforded ≤2 unit lower threshold (Ct) values. To push the detection limit even further, the assay was adapted for reverse transcription PCR on 16S rRNA transcripts. This RT-PCR assay was also linear over 9 log10 units and was sensitive with a LoD of 6.3 copies/reaction.In a dilution-series of experiments, the 16S RT-PCR assay achieved a thousand-fold higher sensitivity than the DNA-targeting assays. For molecular diagnostics, we recommend a real time RT-PCR assay variant in which both DNA and RNA serve as templates (thus, no requirement for DNase treatment). This will at least provide results equaling the DNA-based implementation if no RNA is present but will be superior even at the lowest residual rRNA concentrations.

Analysis of 16S rRNA (rRNA) genes provides a central means of taxonomic classification of bacterial species. Based on presumed sequence identity among species of the Bacillus cereus sensu lato group, the 16S rRNA genes of B. anthracis have been considered unsuitable for diagnosis of the anthrax pathogen. With the recent identification of a single nucleotide polymorphism in some 16S rRNA gene copies, specific identification of B. anthracis becomes feasible. Here, we designed and evaluated a set of in situ, in vitro, and in silico assays to assess the unknown 16S state of B. anthracis from different perspectives. Using a combination of digital PCR, fluorescence in situ hybridization, long-read genome sequencing, and bioinformatics, we were able to detect and quantify a unique 16S rRNA gene allele of B. anthracis (16S-BA-allele). This allele was found in all available B. anthracis genomes and may facilitate differentiation of the pathogen from any close relative. Bioinformatics analysis of 959 B. anthracis SRA data sets inferred that abundances and genomic arrangements of the 16S-BA-allele and the entire rRNA operon copy numbers differ considerably between strains. Expression ratios of 16S-BA-alleles were proportional to the respective genomic allele copy numbers. The findings and experimental tools presented here provide detailed insights into the intra- and intergenomic diversity of 16S rRNA genes and may pave the way for improved identification of B. anthracis and other pathogens with diverse rRNA operons. IMPORTANCE For severe infectious diseases, precise pathogen detection is crucial for antibiotic therapy and patient survival. Identification of Bacillus anthracis, the causative agent of the zoonosis anthrax, can be challenging when querying specific nucleotide sequences such as in small subunit rRNA (16S rRNA) genes, which are commonly used for typing of bacteria. This study analyzed on a broad genomic scale a cryptic and hitherto underappreciated allelic variant of the bacterium’s 16S rRNA genes and their transcripts using a set of in situ, in vitro, and in silico assays and found significant intra- and intergenomic heterogeneity in the distribution of the allele and overall rRNA operon copy numbers. This allelic variation was uniquely species specific, which enabled sensitive pathogen detection on both DNA and transcript levels. The methodology used here is likely also applicable to other pathogens that are otherwise difficult to discriminate from their less harmful relatives.

In July 2018, brucellosis was diagnosed in a German patient without a travel history to regions endemic for Brucella. Microbiological analysis, including whole-genome sequencing, revealed Brucella suis biovar 1 as the etiologic agent. Core-genome-based multilocus sequence-typing analysis placed the isolate in close proximity to strains originating from Argentina. Notably, despite a strong IgM response, the patient did not develop Brucella-specific IgG antibodies during infection. Here, we describe the clinical course of infection, the extensive epidemiological investigations, and discuss possible routes of transmission.

A number of RT-qPCR assays for the detection of SARS-CoV-2 have been published and are listed by the WHO as recommended assays. Furthermore, numerous commercial assays with undisclosed primer and probe sequences are on the market. As the SARS-CoV-2 pandemic progresses, the virus accrues mutations, which in some cases – as seen with the B.1.1.7 variant – can outperform and push back other strains of SARS-CoV-2. If mutations occur in primer or probe binding sites, this can impact RT-qPCR results and impede SARS-CoV-2 diagnostics. Here we tested the effect of primer mismatches on RT-qPCR performance in vitro using synthetic mismatch in vitro transcripts. The effects of the mismatches ranged from a shift in ct values from −0.13 to +7.61. Crucially, we found that a mismatch in the forward primer has a more detrimental effect for PCR performance than a mismatch in the reverse primer. Furthermore, we compared the performance of the original Charité RdRP primer set, which has several ambiguities, with a primer version without ambiguities and found that without ambiguities the ct values are ca. 3 ct lower. Finally, we investigated the shift in ct values observed with the Seegene Allplex kit with the B.1.1.7 SARS-CoV-2 variant and found a three-nucleotide mismatch in the forward primer of the N target.

Cases of pox-like lesions in horses and donkeys have been associated with poxviruses belonging to different genera of the family Poxviridae. These include the orthopoxviruses vaccinia virus (VACV), horsepoxvirus (HPXV) and cowpoxvirus (CPXV), as well as a potentially novel parapoxvirus and molluscum contagiosum virus (MOCV). However, with the exception of VACV, HPXV and CPXV, the genomic characterization of the causative agents remains largely elusive with only single short genome fragments available. Here we present the first full-length genome sequence of an equine molluscum contagiosum-like virus (EMCLV) directly determined from skin biopsies of a horse with generalized papular dermatitis. Histopathological analysis of the lesions revealed severe epidermal hyperplasia with numerous eosinophilic inclusion bodies within keratinocytes. Virions were detected in the lesions in embedded tissue by transmission electron microscopy. The genome sequence determined by next- and third-generation sequencing comprises 166 843 nt with inverted terminal repeats (ITRs) of 3473 nt. Overall, 20 of the predicted 159 ORFs have no equivalents in other poxviruses. Intriguingly, two of these ORFs were identified to encode homologues of mammalian proteins involved in immune signalling pathways, namely secreted and transmembrane protein 1 (SECTM1) and insulin growth factor-like family receptor 1 (IGFLR1), that were not described in any virus family so far. Phylogenetic analysis with all relevant representatives of the Poxviridae suggests that EMCLV should be nominated as a new species within the genus Molluscipoxvirus.

Controlling and monitoring the still ongoing severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic regarding geographical distribution, evolution, and emergence of new mutations of the SARS-CoV-2 virus is only possible due to continuous next-generation sequencing (NGS) and sharing sequence data worldwide. Efficient sequencing strategies enable the retrieval of increasing numbers of high-quality, full-length genomes and are, hence, indispensable. Two opposed enrichment methods, tiling multiplex PCR and sequence hybridization by bait capture, have been established for SARS-CoV-2 sequencing and are both frequently used, depending on the quality of the patient sample and the question at hand. Here, we focused on the evaluation of the sequence hybridization method by studying five commercially available sequence capture bait panels with regard to sensitivity and capture efficiency. We discovered the SARS-CoV-2-specific panel of Twist Bioscience to be the most efficient panel, followed by two respiratory panels from Twist Bioscience and Illumina, respectively. Our results provide on the one hand a decision basis for the sequencing community including a computation for using the full capacity of the flow cell and on the other hand potential improvements for the manufacturers. IMPORTANCE Sequencing the genomes of the circulating SARS-CoV-2 strains is the only way to monitor the viral spread and evolution of the virus. Two different approaches, namely, tiling multiplex PCR and sequence hybridization by bait capture, are commonly used to fulfill this task. This study describes for the first time a combined approach of droplet digital PCR (ddPCR) and NGS to evaluate five commercially available sequence capture panels targeting SARS-CoV-2. In doing so, we were able to determine the most sensitive and efficient capture panel, distinguish the mode of action of the various bait panels, and compute the number of read pairs needed to recover a high-quality full-length genome. By calculating the minimum number of read pairs needed, we are providing optimized flow cell loading conditions for all sequencing laboratories worldwide that are striving for maximizing sequencing output and simultaneously minimizing time, costs, and sequencing resources.

Summary The Flexible Taxonomy Database framework provides a method for modification and merging official and custom taxonomic databases to create improved databases. Using such databases will increase accuracy and precision of existing methods to classify sequence reads. Availability and implementation Source code is freely available at https://github.com/FOI-Bioinformatics/flextaxd and installable through Bioconda.

Q (query) fever is an infectious zoonotic disease, caused by the Gram-negative bacteria Coxiella burnetii, that sometimes is transmitted to humans from small ruminants like sheep, goat and cattle. Although the disease has been studied since decades, it still represents a threat due to sporadic outbreaks across farms in Europe. The reason for this has been linked to the interaction of several dynamic factors including reservoir type and vector diversity. One important control measure we have identified is a central platform for Coxiella typing data management. This is particularly important in the case of an outbreak where the nature of the pathogen and type would need to be rapidly identified and compared to existing isolates as well as further documented and made available for researcher to aid future investigations. The existing platforms are focused on MLVA (multiple locus VNTR analysis) genotyping. We have designed and implemented an online, open, web-based platform called CoxBase (https://coxbase.q-gaps.de), that is capable of in silico genotyping of completely or minimally assembled Coxiella sequences using five different typing methods, included with a database that holds genotyping information of more than 400 Coxiella isolates with their metadata such as host type, source and year of isolation together with further metadata information. Also, it includes a query and submission interface for interrogating existing isolates and depositing new isolates. Attractive visualization features includes maps showing the geographical source of the isolates and plots that can be used to summarize isolates metadata on a country level. We tested our in silico typing method on 50 Coxiella genomes downloaded from the RefSeq database and we could type almost all the genomes except for cases where the sequences are poorly assembled. We identified new spacer sequences using our MST in silico typing methods, and could categorize adaA gene phenotypes for all 50 genomes as well as their plasmid types.

We are currently facing a pandemic of COVID-19, caused by a spillover from an animal-originating coronavirus to humans occuring in the Wuhan region, China, in December 2019. From China the virus has spread to 188 countries and regions worldwide, reaching the Sahel region on the 2nd of March 2020. Since whole genome sequencing (WGS) data is very crucial to understand the spreading dynamics of the ongoing pandemic, but only limited sequence data is available from the Sahel region to date, we have focused our efforts on generating the first Malian sequencing data available. Screening of 217 Malian patient samples for the presence of SARS-CoV-2 resulted in 38 positive isolates from which 21 whole genome sequences were generated. Our analysis shows that both, the early A (19B) and the fast evolving B (20A/C) clade, are present in Mali indicating multiple and independent introductions of the SARS-CoV-2 to the Sahel region.

Summary Background In December, 2019, the newly identified severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) emerged in Wuhan, China, causing COVID-19, a respiratory disease presenting with fever, cough, and often pneumonia. WHO has set the strategic objective to interrupt spread of SARS-CoV-2 worldwide. An outbreak in Bavaria, Germany, starting at the end of January, 2020, provided the opportunity to study transmission events, incubation period, and secondary attack rates. Methods A case was defined as a person with SARS-CoV-2 infection confirmed by RT-PCR. Case interviews were done to describe timing of onset and nature of symptoms and to identify and classify contacts as high risk (had cumulative face-to-face contact with a confirmed case for ≥15 min, direct contact with secretions or body fluids of a patient with confirmed COVID-19, or, in the case of health-care workers, had worked within 2 m of a patient with confirmed COVID-19 without personal protective equipment) or low risk (all other contacts). High-risk contacts were ordered to stay at home in quarantine for 14 days and were actively followed up and monitored for symptoms, and low-risk contacts were tested upon self-reporting of symptoms. We defined fever and cough as specific symptoms, and defined a prodromal phase as the presence of non-specific symptoms for at least 1 day before the onset of specific symptoms. Whole genome sequencing was used to confirm epidemiological links and clarify transmission events where contact histories were ambiguous; integration with epidemiological data enabled precise reconstruction of exposure events and incubation periods. Secondary attack rates were calculated as the number of cases divided by the number of contacts, using Fisher's exact test for the 95% CIs. Findings Patient 0 was a Chinese resident who visited Germany for professional reasons. 16 subsequent cases, often with mild and non-specific symptoms, emerged in four transmission generations. Signature mutations in the viral genome occurred upon foundation of generation 2, as well as in one case pertaining to generation 4. The median incubation period was 4·0 days (IQR 2·3–4·3) and the median serial interval was 4·0 days (3·0–5·0). Transmission events were likely to have occurred presymptomatically for one case (possibly five more), at the day of symptom onset for four cases (possibly five more), and the remainder after the day of symptom onset or unknown. One or two cases resulted from contact with a case during the prodromal phase. Secondary attack rates were 75·0% (95% CI 19·0–99·0; three of four people) among members of a household cluster in common isolation, 10·0% (1·2–32·0; two of 20) among household contacts only together until isolation of the patient, and 5·1% (2·6–8·9; 11 of 217) among non-household, high-risk contacts. Interpretation Although patients in our study presented with predominately mild, non-specific symptoms, infectiousness before or on the day of symptom onset was substantial. Additionally, the incubation period was often very short and false-negative tests occurred. These results suggest that although the outbreak was controlled, successful long-term and global containment of COVID-19 could be difficult to achieve. Funding All authors are employed and all expenses covered by governmental, federal state, or other publicly funded institutions.

Fast turnaround times are of utmost importance for biomedical reconnaissance, particularly regarding dangerous pathogens. Recent advances in sequencing technology and its devices allow sequencing within a short time frame outside stationary laboratories close to the epicenter of the outbreak. In our study, we evaluated the portable sequencing device MinION as part of a rapidly deployable laboratory specialized in identification of highly pathogenic agents. We tested the device in the course of a NATO live agent exercise in a deployable field laboratory in hot climate conditions. The samples were obtained from bio-terroristic scenarios that formed part of the exercise and contained unknown bacterial agents. To simulate conditions of a resource-limited remote deployment site, we operated the sequencer without internet access. Using a metagenomic approach, we were able to identify the causative agent in the analyzed samples. Furthermore, depending on the obtained data, we were able to perform molecular typing down to strain level. In our study we challenged the device and discuss advances as well as remaining limitations for sequencing biological samples outside of stationary laboratories. Nevertheless, massive parallel sequencing as a non-selective methodology yields important information and is able to support outbreak investigation − even in the field.

The autosomal recessive genetic disorder Nijmegen breakage syndrome (NBS) is characterized by a defect in DNA double-strand break repair protein nibrin and chromosome instability associated with a high predisposition to cancer. Here we hypothesized that impaired nibrin/MRE11/RAD50 telomere maintenance complex may also affect telomere length and modulate the cancer phenotype.Telomere length was studied in blood from 38 homozygous and 27 heterozygous individuals, in one homozygous fetus, and in sex NBS lymphoblastoid cell lines (all with the founder mutation c.657_661del5), and in three humanized Nbs mice, using qPCR, TRF and Q-FISH.Telomere lengths were markedly but uniformly reduced to 20-40% of healthy controls. There was no correlation between telomere length and severity of clinical phenotype or age of death. By contrast, individual patients with very short telomeres displayed long survival times after cancer manifestation. Mildly accelerated telomere attrition was found in older NBS heterozygotes. In the NBS-fetus, the spinal cord, brain and heart had the longest telomeres, skin the shortest. Humanized Nbs mice (with much longer telo-meres than those in human beings) did not show accelerated telomere attrition.Our data clearly show that NBS is a secondary telomeropathy with unique features. Te- lomere attrition in NBS may cause genetic instability and contribute to the high cancer incidence in NBS. On the other hand, short telomeres may prevent an even worse pheno-type when a tumor has developed. These data may help to understand the high cancer rate in NBS and also the bifunctional role of telomere shortening in cancerogenesis.Author SummaryDNA damage is harmful because it leads to mutations in genes that initiate or accelerate cancerogenesis. The devastating consequences of DNA damage are manifested in diseases with non-functional repair pathways such as Nijmegen breakage syndrome (NBS). A common feature of these diseases is a high tumor incidence. However, cancer incidence varies and is not clear why it is highest for NBS. In a previous study, we have shown that the underlying nebrin mutation not only leads to defective DNA repair but also to higher degree of oxidative stress that generates further DNA lesions. Nibrin may play also an important role in protecting chromosome ends, the telomeres, from inap-propriate DNA repair. Therefore we examined the telomere length in NBS and show markedly reduced values in affected patients but not in NBC mice (with much milder phenotype and longer telomeres). Telomere attrition contributes to genetic instability and may thus contribute to the high cancer incidence in NBS. Individual patients with very short telomeres, however, displayed long survival times after cancer manifestation. Thus, short telomeres may also prevent an even worse phenotype when a tumor has developed. These data are fundamental to understanding the high cancer rate in NBS and also the bifunctional role of telomere shortening in cancer.

The zoonosis anthrax caused by the bacterium Bacillus anthracis has a broad geographical distribution. Active enzootic areas are typically located away from central and northern Europe where cases of the disease occur only sporadically and in limited numbers. In contrast, a few out of the 64 districts of Bangladesh are hyper-endemic for anthrax and there the disease causes major losses in live-stock. In this study we genotyped eight strains of B. anthracis collected from the districts of Sirajganj and Tangail in 2013. All these strains belonged to canSNP group A.Br.001/002 Sterne differing only in a few of 31 tandem-repeat (MLVA)-markers. Whole genome sequences were obtained from five of these strains and compared with genomic information of B. anthracis strains originating from various geographical locations. Characteristic signatures were detected defining two “Bangladesh” clusters potentially useful for rapid molecular epidemiology. From this data high-resolution PCR assays were developed and subsequently tested on additional isolates from Bangladesh and Central Europe. Remarkably, this comparative genomic analysis focusing on SNP-discovery revealed a close genetic relationship between these strains from Bangladesh and historic strains collected between 1991 and 2008 in The Netherlands and Germany, respectively. Possible explanations for these phylogenetic relationships are discussed.

Coxiella burnetii is a pathogen causing Q fever in domestic animals and humans. Seabirds have been implicated as possible reservoirs of this bacterium in the Arabian Gulf and in the Western Indian Ocean. Recently, Coxiella species closely related to C. burnetii was detected from ticks collected from oil rigs used as roosting areas by Socotra Cormorants (Phalacrocorax nigrogularis) in the western Arabian Gulf. We collected ticks from the largest breeding colony of Socotra Cormorants in the United Arab Emirates on the eastern extreme of the species' breeding range to determine the prevalence of C. burnetii and evaluate its role as a wild reservoir. All ticks were identified as Ornithodoros muesebecki and genomic DNA was extracted from larval and nymph/adult tick pools. Multiplex PCR tests were performed targeting three C. burnetii specific genes. C. burnetii was not detected although a Coxiella-like endosymbiont was identified that was closely related to Coxiella symbionts from Ornithodoros capensis ticks. Because domestic and wild ungulates are the primary source of C. burnetii, we suggest that the presence of free-ranging, native and non-native ungulates in some off-shore islands in the Arabian Gulf could disseminate C. burnetii to seabirds. More comprehensive studies on seabird colonies are needed to better understand the diversity and prevalence of Coxiella symbionts and to establish if C. burnetii is endemic on some of these islands.

Armillaria species are both devastating forest pathogens and some of the largest terrestrial organisms on Earth. They forage for hosts and achieve immense colony sizes via rhizomorphs, root-like multicellular structures of clonal dispersal. Here, we sequenced and analysed the genomes of four Armillaria species and performed RNA sequencing and quantitative proteomic analysis on the invasive and reproductive developmental stages of A. ostoyae. Comparison with 22 related fungi revealed a significant genome expansion in Armillaria, affecting several pathogenicity-related genes, lignocellulose-degrading enzymes and lineage-specific genes expressed during rhizomorph development. Rhizomorphs express an evolutionarily young transcriptome that shares features with the transcriptomes of both fruiting bodies and vegetative mycelia. Several genes show concomitant upregulation in rhizomorphs and fruiting bodies and share cis-regulatory signatures in their promoters, providing genetic and regulatory insights into complex multicellularity in fungi. Our results suggest that the evolution of the unique dispersal and pathogenicity mechanisms of Armillaria might have drawn upon ancestral genetic toolkits for wood-decay, morphogenesis and complex multicellularity., Nature Ecology & Evolution, 1, ISSN:2397-334X

We report the draft genome sequence of clindamycin-resistant Bacillus safensis strain Ingolstadt isolated from a patient with bacterial colonization after heart surgery. The draft genome comprises 3.75 Mbp and harbors 3,793 predicted protein-encoding genes and a small plasmid.

We report the draft genome sequence of clindamycin-resistant Bacillus safensis strain Ingolstadt isolated from a patient with bacterial colonization after heart surgery. The draft genome comprises 3.75 Mbp and harbors 3,793 predicted protein-encoding genes and a small plasmid.

The causative agent of Q fever, Coxiella burnetii, is a query agent occurring naturally all over the world. We studied 104 German Coxiella burnetii strains/DNA samples obtained between 1969 and 2011 using a 14 microsatellite marker Multiple-locus variable-number of tandem repeat (VNTR) analysis (MLVA) technique. We were able to divide our collection into 32 different genotypes clustered into four major groups (A-D). Two of these (A and C) formed predominant clonal complexes that covered 97% of all studied samples. Group C consisted exclusively of cattle-associated isolates/DNA specimens, while group A comprised all other affected species including all sheep-derived strains/DNA samples. Within this second cluster, two major genotypes (A1, A2) were identified. Genotype A2 occurred in strains isolated from ewes in northern and central Germany, whereas genotype A1 was found in most areas of Germany. MLVA analysis of C. burnetii strains from neighbouring countries revealed a close relationship to German strains. We thus hypothesize that there is a western and central European cluster of C. burnetii. We identified predominant genotypes related to relevant host species and geographic regions which is in line with findings of the Dutch Q fever outbreak (2007-2010). Furthermore three of our analyzed German strains are closely related to the Dutch outbreak clone. These findings support the theory of predominant genotypes in the context of regional outbreaks. Our results show that a combination of 8 MLVA markers provides the highest discriminatory power for attributing C. burnetii isolates to genotypes. For future epidemiological studies we propose the use of three MLVA markers for easy and rapid classification of C. burnetii into 4 main clusters.

eggNOG is a public resource that provides Orthologous Groups (OGs) of proteins at different taxonomic levels, each with integrated and summarized functional annotations. Developments since the latest public release include changes to the algorithm for creating OGs across taxonomic levels, making nested groups hierarchically consistent. This allows for a better propagation of functional terms across nested OGs and led to the novel annotation of 95 890 previously uncharacterized OGs, increasing overall annotation coverage from 67% to 72%. The functional annotations of OGs have been expanded to also provide Gene Ontology terms, KEGG pathways and SMART/Pfam domains for each group. Moreover, eggNOG now provides pairwise orthology relationships within OGs based on analysis of phylogenetic trees. We have also incorporated a framework for quickly mapping novel sequences to OGs based on precomputed HMM profiles. Finally, eggNOG version 4.5 incorporates a novel data set spanning 2605 viral OGs, covering 5228 proteins from 352 viral proteomes. All data are accessible for bulk downloading, as a web-service, and through a completely redesigned web interface. The new access points provide faster searches and a number of new browsing and visualization capabilities, facilitating the needs of both experts and less experienced users. eggNOG v4.5 is available at http://eggnog.embl.de.

Due to their economic relevance, the study of plant pathogen interactions is of importance. However, elucidating these interactions and their underlying molecular mechanisms remains challenging since both host and pathogen need to be fully genetically accessible organisms. Here we present milestones in the establishment of a new biotrophic model pathosystem: Ustilago bromivora and Brachypodium sp. We provide a complete toolset, including an annotated fungal genome and methods for genetic manipulation of the fungus and its host plant. This toolset will enable researchers to easily study biotrophic interactions at the molecular level on both the pathogen and the host side. Moreover, our research on the fungal life cycle revealed a mating type bias phenomenon. U. bromivora harbors a haplo-lethal allele that is linked to one mating type region. As a result, the identified mating type bias strongly promotes inbreeding, which we consider to be a potential speciation driver. DOI: http://dx.doi.org/10.7554/eLife.20522.001, eLife digest Fungi cause many diseases in plants, and reduce the yield of important crops like wheat, corn and rice – all of which belong to the family of grasses. Much research into how disease-causing fungi infect plants will look at a given fungus that infects a specific plant in order to understand plant diseases in general. Over the years, scientists have generated suites of research tools to study these pairs of fungi and plants. However, many of these organism pairs (often called “model pathosystems”) have drawbacks when it comes to research in the laboratory, either on the side of the fungus or the side of plant. Brachypodium is a small grass that grows quickly and, unlike crop plants, it grows well in the laboratory. These characteristics make Brachypodium a promising model organism for studying many aspects of plant biology. Recently, a fungus called Ustilago bromivora – which is related to a fungus that infects corn – was reported to infect Brachypodium. This raised the question: could this fungus and this small grass become a new model pathosystem? Rabe, Bosch et al. set out to answer this question and now provide a toolkit that will help to establish U. bromivora and Brachypodium as a new model pathosystem. In all of U. bromivora’s close relatives, two compatible strains must meet and mate before the fungus can infect the plant; first Rabe, Bosch et al. confirmed that this is also the case for U. bromivora. Studying the life cycle of the U. bromivora fungus also unexpectedly revealed that while both mating partners are needed to infect the plant, only one of the strains survives outside of the plant after the infection. This phenomenon, referred to as a “mating type bias”, has been described for a few other related fungi. Next, Rabe, Bosch et al. conducted a genetic screen and identified two compatible strains that can grow without the plant as yeast-like cells. This means that these cells can be manipulated genetically, and indeed protocols to grow and genetically engineer the fungus and plant to address different research questions are included in the toolkit as well. Other new tools include the complete genetic sequence of the fungus with all its genes annotated, and a dataset of which genes are active in U. bromivora growing yeast-like in liquid culture versus those active when the fungus grows as a pathogen inside the plant. Together these new tools and datasets will provide a foundation to study different aspects of the interactions between grasses and disease-causing fungi. This in turn may lead to new methods to reduce fungal growth and reduce yield losses caused by fungal diseases in crop plants. Finally, the discovery that U. bromivora shows a mating type bias could provide a starting point for future studies into sexual reproduction in fungi and how new species arise. DOI: http://dx.doi.org/10.7554/eLife.20522.002

It is generally accepted that functionally important RNA structure is more conserved than sequence due to compensatory mutations that may alter the sequence without disrupting the structure. For small RNA molecules sequence-structure relationships are relatively well understood. However, structural bioinformatics of mRNAs is still in its infancy due to a virtual absence of experimental data. This report presents the first quantitative assessment of sequence-structure divergence in the coding regions of mRNA molecules based on recently published transcriptome-wide experimental determination of their base paring patterns. Structural resemblance in paralogous mRNA pairs quickly drops as sequence identity decreases from 100% to 85-90%. Structures of mRNAs sharing sequence identity below roughly 85% are essentially uncorrelated. This outcome is in dramatic contrast to small functional non-coding RNAs where sequence and structure divergence are correlated at very low levels of sequence similarity. The fact that very similar mRNA sequences can have vastly different secondary structures may imply that the particular global shape of base paired elements in coding regions does not play a major role in modulating gene expression and translation efficiency. Apparently, the need to maintain stable three-dimensional structures of encoded proteins places a much higher evolutionary pressure on mRNA sequences than on their RNA structures.

In the version of this Article originally published, it was incorrectly stated that "16,687 protein-coding genes were inferred for the most recent common ancestor (MRCA) of Armillaria"; the value was incorrect and it should have read "15,787". This has now been corrected.

Motivation: Accurate automatic assignment of protein functions remains a challenge for genome annotation. We have developed and compared the automatic annotation of four bacterial genomes employing a 5-fold cross-validation procedure and several machine learning methods. Results: The analyzed genomes were manually annotated with FunCat categories in MIPS providing a gold standard. Features describing a pair of sequences rather than each sequence alone were used. The descriptors were derived from sequence alignment scores, InterPro domains, synteny information, sequence length and calculated protein properties. Following training we scored all pairs from the validation sets, selected a pair with the highest predicted score and annotated the target protein with functional categories of the prototype protein. The data integration using machine-learning methods provided significantly higher annotation accuracy compared to the use of individual descriptors alone. The neural network approach showed the best performance. The descriptors derived from the InterPro domains and sequence similarity provided the highest contribution to the method performance. The predicted annotation scores allow differentiation of reliable versus non-reliable annotations. The developed approach was applied to annotate the protein sequences from 180 complete bacterial genomes. Availability: The FUNcat Annotation Tool (FUNAT) is available on-line as Web Services at http://mips.gsf.de/proj/funat Contact: i.tetko@gsf.de Supplementary information: Supplementary data are available at Bioinformatics online.

Coxiella burnetii strain AuQ01 was isolated from the serum of an Australian acute Q fever patient and represents the first whole genome from this historical Q fever country. This new genome shows distinct differences from existing genomic data and will enhance the understanding of this query pathogen.

International audience; Bacterial infectious diseases are the result of multifactorial processes affected by the interplay between virulence factors and host targets. The host-Pseudomonas and Coxiella interaction database (HoPaCI-DB) is a publicly available manually curated integrative database (http://mips.helmholtz-muenchen.de/HoPaCI/) of host-pathogen interaction data from Pseudomonas aeruginosa and Coxiella burnetii. The resource provides structured information on 3585 experimentally validated interactions between molecules, bioprocesses and cellular structures extracted from the scientific literature. Systematic annotation and interactive graphical representation of disease networks make HoPaCI-DB a versatile knowledge base for biologists and network biology approaches.

From Blight to Powdery Mildew Pathogenic effects of microbes on plants have widespread consequences. Witness, for example, the cultural upheavals driven by potato blight in the 1800s. A variety of microbial pathogens continue to afflict crop plants today, driving both loss of yield and incurring the increased costs of control mechanisms. Now, four reports analyze microbial genomes in order to understand better how plant pathogens function (see the Perspective by Dodds ). Raffaele et al. (p. 1540 ) describe how the genome of the potato blight pathogen accommodates transfer to different hosts. Spanu et al. (p. 1543 ) analyze what it takes to be an obligate biotroph in barley powdery mildew, and Baxter et al. (p. 1549 ) ask a similar question for a natural pathogen of Arabidopsis . Schirawski et al. (p. 1546 ) compared genomes of maize pathogens to identify virulence determinants. Better knowledge of what in a genome makes a pathogen efficient and deadly is likely to be useful for improving agricultural crop management and breeding.

The MIPS Fusarium graminearum Genome Database (FGDB) was established as a comprehensive genome database on one of the most devastating fungal plant pathogens of wheat, barley and maize. The current version of FGDB v3.1 provides information on the full manually revised gene set based on the Broad Institute assembly FG3 genome sequence. The results of gene prediction tools were integrated with the help of comparative data on related species to result in a set of 13.718 annotated protein coding genes. This rigorous approach involved adding or modifying gene models and represents a coding sequence gold standard for the genus Fusarium. The gene loci improvements results in 2461 genes which either are new or have different structures compared to the Broad Institute assembly 3 gene set. Moreover the database serves as a convenient entry point to explore expression data results and to obtain information on the Affymetrix GeneChip probe sets. The resource is accessible on http://mips.gsf.de/genre/proj/FGDB/.

The acute disease antigen A (adaA) gene is believed to be associated with Coxiella burnetii strains causing acute Q fever. The detailed analysis of the adaA genomic region of 23 human- and 86 animal-derived C. burnetii isolates presented in this study reveals a much more polymorphic appearance and distribution of the adaA gene, resulting in a classification of C. burnetii strains of better differentiation than previously anticipated. Three different genomic variants of the adaA gene were identified which could be detected in isolates from acute and chronic patients, rendering the association of adaA positive strains with acute Q fever disease disputable. In addition, all adaA positive strains in humans and animals showed the occurrence of the QpH1 plasmid. All adaA positive isolates of acute human patients except one showed a distinct SNP variation at position 431, also predominant in sheep strains, which correlates well with the observation that sheep are a major source of human infection. Furthermore, the phylogenetic analysis of the adaA gene revealed three deletion events and supported the hypothesis that strain Dugway 5J108-111 might be the ancestor of all known C. burnetii strains. Based on our findings, we could confirm the QpDV group and we were able to define a new genotypic cluster. The adaA gene polymorphisms shown here improve molecular typing of Q fever, and give new insights into microevolutionary adaption processes in C. burnetii.

The PEDANT genome database provides exhaustive annotation of nearly 3000 publicly available eukaryotic, eubacterial, archaeal and viral genomes with more than 4.5 million proteins by a broad set of bioinformatics algorithms. In particular, all completely sequenced genomes from the NCBI's Reference Sequence collection (RefSeq) are covered. The PEDANT processing pipeline has been sped up by an order of magnitude through the utilization of pre-calculated similarity information stored in the similarity matrix of proteins (SIMAP) database, making it possible to process newly sequenced genomes immediately as they become available. PEDANT is freely accessible to academic users at http://pedant.gsf.de. For programmatic access Web Services are available at http://pedant.gsf.de/webservices.jsp.

High volumes of serial computational tasks in bioinformatics, such as homology searches or profile matching, are often executed in distributed environments using classical batching systems like LSF or Sun Grid-Engine. Despite their simple usability they are limited to organizationally owned resources. In contrast to proprietary projects that implement large-scale grids we report on a grid-enabled solution for sequence homology and protein domain searches using BOINC, the Berkeley Open Infrastructure for Network Computing. We demonstrate that BOINC provides a powerful and versatile platform for public resource computing in bioinformatics that makes large-scale pre-calculation of sequence analyses feasible. The FASTA-/Smith-Waterman-and HMMer applications for BOINC are freely available from the authors upon request. Data from SIMAP is publicly available through Web-Services at http://mips.gsf.de/simap.