Marcus Lechner

Orthology detection combining clustering and synteny for very large datasets.

The elucidation of orthology relationships is an important step both in gene function prediction as well as towards understanding patterns of sequence evolution. Orthology assignments are usually derived directly from sequence similarities for large data because more exact approaches exhibit too high computational costs. Here we present PoFF, an extension for the standalone tool Proteinortho, which enhances orthology detection by combining clustering, sequence similarity, and synteny. In the course of this work, FFAdj-MCS, a heuristic that assesses pairwise gene order using adjacencies (a similarity measure related to the breakpoint distance) was adapted to support multiple linear chromosomes and extended to detect duplicated regions. PoFF largely reduces the number of false positives and enables more fine-grained predictions than purely similarity-based approaches. The extension maintains the low memory requirements and the efficient concurrency options of its basis Proteinortho, making the software applicable to very large datasets.

Phylogenomics with paralogs

Significance We demonstrate that the distribution of paralogs in large gene families contains in itself sufficient phylogenetic signal to infer fully resolved species phylogenies. This source of phylogenetic information is independent of information contained in orthologous sequences and is resilient against horizontal gene transfer. An important consequence is that phylogenomics data sets need not be restricted to 1:1 orthologs. Phylogenomics heavily relies on well-curated sequence data sets that comprise, for each gene, exclusively 1:1 orthologos. Paralogs are treated as a dangerous nuisance that has to be detected and removed. We show here that this severe restriction of the data sets is not necessary. Building upon recent advances in mathematical phylogenetics, we demonstrate that gene duplications convey meaningful phylogenetic information and allow the inference of plausible phylogenetic trees, provided orthologs and paralogs can be distinguished with a degree of certainty. Starting from tree-free estimates of orthology, cograph editing can sufficiently reduce the noise to find correct event-annotated gene trees. The information of gene trees can then directly be translated into constraints on the species trees. Although the resolution is very poor for individual gene families, we show that genome-wide data sets are sufficient to generate fully resolved phylogenetic trees, even in the presence of horizontal gene transfer.

Genetic characterization of Tribec virus and Kemerovo virus, two tick-transmitted human-pathogenic Orbiviruses

We determined the complete genome sequences of Tribeč virus (TRBV) and Kemerovo virus (KEMV), two tick-transmitted Orbiviruses that can cause diseases of the central nervous system and that are currently classified into the Great Island virus serogroup. VP2 proteins of TRBV and KEMV show very low sequence similarity to the homologous VP4 protein of tick-transmitted Great Island virus (GIV). The new sequence data support previous serological classification of these Orbiviruses into the Kemerovo serogroup, which is different from the Great Island virus serogroup. Genome segment 9 of TRBV and KEMV encodes several overlapping ORFs in the +1 reading frame relative to VP6(Hel). A co-phylogenetic analysis indicates a host switch from insect-borne Orbiviruses toward Ixodes species, which is in disagreement with previously published data.

The correlation of genome size and DNA methylation rate in metazoans

Total DNA methylation rates are well known to vary widely between different metazoans. The phylogenetic distribution of this variation, however, has not been investigated systematically. We combine here publicly available data on methylcytosine content with the analysis of nucleotide compositions of genomes and transcriptomes of 78 metazoan species to trace the evolution of abundance and distribution of DNA methylation. The depletion of CpG and the associated enrichment of TpG and CpA dinucleotides are used to infer the intensity and localization of germline CpG methylation and to estimate its evolutionary dynamics. We observe a positive correlation of the relative methylation of CpG motifs with genome size. We tested this trend successfully by measuring total DNA methylation with LC/MS in orthopteran insects with very different genome sizes: house crickets, migratory locusts and meadow grasshoppers. We hypothesize that the observed correlation between methylation rate and genome size is due to a dependence of both variables from long-term effective population size and is driven by the accumulation of repetitive sequences that are typically methylated during periods of small population sizes. This process may result in generally methylated, large genomes such as those of jawed vertebrates. In this case, the emergence of a novel demethylation pathway and of novel reader proteins for methylcytosine may have enabled the usage of cytosine methylation for promoter-based gene regulation. On the other hand, persistently large populations may lead to a compression of the genome and to the loss of the DNA methylation machinery, as observed, e.g., in nematodes.

Mechanistic comparison of Bacillus subtilis 6S-1 and 6S-2 RNAs—commonalities and differences

Bacterial 6S RNAs bind to the housekeeping RNA polymerase (σ(A)-RNAP in Bacillus subtilis) to regulate transcription in a growth phase-dependent manner. B. subtilis expresses two 6S RNAs, 6S-1 and 6S-2 RNA, with different expression profiles. We show in vitro that 6S-2 RNA shares hallmark features with 6S-1 RNA: Both (1) are able to serve as templates for pRNA transcription; (2) bind with comparable affinity to σ(A)-RNAP; (3) are able to specifically inhibit transcription from DNA promoters, and (4) can form stable 6S RNA:pRNA hybrid structures that (5) abolish binding to σ(A)-RNAP. However, pRNAs of equal length dissociate faster from 6S-2 than 6S-1 RNA, owing to the higher A,U-content of 6S-2 pRNAs. This could have two mechanistic implications: (1) Short 6S-2 pRNAs (<10 nt) dissociate faster instead of being elongated to longer pRNAs, which could make it more difficult for 6S-2 RNA-stalled RNAP molecules to escape from the sequestration; and (2) relative to 6S-1 RNA, 6S-2 pRNAs of equal length will dissociate more rapidly from 6S-2 RNA after RNAP release, which could affect pRNA turnover or the kinetics of 6S-2 RNA binding to a new RNAP molecule. As 6S-2 pRNAs have not yet been detected in vivo, we considered that cellular RNAP release from 6S-2 RNA might occur via 6S-1 RNA displacing 6S-2 RNA from the enzyme, either in the absence of pRNA transcription or upon synthesis of very short 6S-2 pRNAs (∼ 5-mers, which would escape detection by deep sequencing). However, binding competition experiments argued against these possibilities.

Genomewide comparison and novel ncRNAs of Aquificales

BackgroundThe Aquificales are a diverse group of thermophilic bacteria that thrive in terrestrial and marine hydrothermal environments. They can be divided into the families Aquificaceae, Desulfurobacteriaceae and Hydrogenothermaceae. Although eleven fully sequenced and assembled genomes are available, only little is known about this taxonomic order in terms of RNA metabolism.ResultsIn this work, we compare the available genomes, extend their protein annotation, identify regulatory sequences, annotate non-coding RNAs (ncRNAs) of known function, predict novel ncRNA candidates, show idiosyncrasies of the genetic decoding machinery, present two different types of transfer-messenger RNAs and variations of the CRISPR systems. Furthermore, we performed a phylogenetic analysis of the Aquificales based on entire genome sequences, and extended this by a classification among all bacteria using 16S rRNA sequences and a set of orthologous proteins.Combining several in silico features (e.g. conserved and stable secondary structures, GC-content, comparison based on multiple genome alignments) with an in vivo dRNA-seq transcriptome analysis of Aquifex aeolicus, we predict roughly 100 novel ncRNA candidates in this bacterium.ConclusionsWe have here re-analyzed the Aquificales, a group of bacteria thriving in extreme environments, sharing the feature of a small, compact genome with a reduced number of protein and ncRNA genes. We present several classical ncRNAs and riboswitch candidates. By combining in silico analysis with dRNA-seq data of A. aeolicus we predict nearly 100 novel ncRNA candidates.

A genomic view of food-related and probiotic Enterococcus strains

Abstract The study of enterococcal genomes has grown considerably in recent years. While special attention is paid to comparative genomic analysis among clinical relevant isolates, in this study we performed an exhaustive comparative analysis of enterococcal genomes of food origin and/or with potential to be used as probiotics. Beyond common genetic features, we especially aimed to identify those that are specific to enterococcal strains isolated from a certain food-related source as well as features present in a species-specific manner. Thus, the genome sequences of 25 Enterococcus strains, from 7 different species, were examined and compared. Their phylogenetic relationship was reconstructed based on orthologous proteins and whole genomes. Likewise, markers associated with a successful colonization (bacteriocin genes and genomic islands) and genome plasticity (phages and clustered regularly interspaced short palindromic repeats) were investigated for lifestyle specific genetic features. At the same time, a search for antibiotic resistance genes was carried out, since they are of big concern in the food industry. Finally, it was possible to locate 1617 FIGfam families as a core proteome universally present among the genera and to determine that most of the accessory genes code for hypothetical proteins, providing reasonable hints to support their functional characterization.

Identification of specific residues in avian influenza A virus NS1 that enhance viral replication and pathogenicity in mammalian systems

Reassortment of their segmented genomes allows influenza A viruses (IAV) to gain new characteristics, which potentially enable them to cross the species barrier and infect new hosts. Improved replication was observed for reassortants of the strictly avian IAV A/FPV/Rostock/34 (FPV, H7N1) containing the NS segment from A/Goose/Guangdong/1/1996 (GD, H5N1), but not for reassortants containing the NS segment of A/Mallard/NL/12/2000 (MA, H7N3). The NS1 of GD and MA differ only in 8 aa positions. Here, we show that efficient replication of FPV-NSMA-derived mutants was linked to the presence of a single substitution (D74N) and more prominently to a triple substitution (P3S+R41K+D74N) in the NS1MA protein. The substitution(s) led to (i) increased virus titres, (ii) larger plaque sizes and (iii) increased levels and faster kinetics of viral mRNA and protein accumulation in mammalian cells. Interestingly, the NS1 substitutions did not affect viral growth characteristics in avian cells. Furthermore, we show that an FPV mutant with N74 in the NS1 (already possessing S3+K41) is able to replicate and cause disease in mice, demonstrating a key role of NS1 in the adaptation of avian IAV to mammalian hosts. Our data suggest that (i) adaptation to mammalian hosts does not necessarily compromise replication in the natural (avian) host and (ii) very few genetic changes may pave the way for zoonotic transmission. The study reinforces the need for close surveillance and characterization of circulating avian IAV to identify genetic signatures that indicate a potential risk for efficient transmission of avian strains to mammalian hosts.

Minimal and RNA-free RNase P in Aquifex aeolicus

Significance RNase P is a tRNA-processing enzyme of unique architectural diversity: either a catalytic RNA plus one or more (up to 10) proteins, or one (or three) unrelated proteins only. We identified yet another enzyme form in the bacterium Aquifex aeolicus, a 23-kDa protein and the smallest known form of RNase P. Apparently, it was acquired by horizontal gene transfer from Archaea. In some other bacteria and many archaea, it is simultaneously present with the presumably more ancient RNA-based enzyme form. Bacteria with both activities may represent the missing link of RNase P evolution, a transition state that had also been once traversed by the Aquificaceae, which, however, later lost their RNA-based RNase P. RNase P is an essential tRNA-processing enzyme in all domains of life. We identified an unknown type of protein-only RNase P in the hyperthermophilic bacterium Aquifex aeolicus: Without an RNA subunit and the smallest of its kind, the 23-kDa polypeptide comprises a metallonuclease domain only. The protein has RNase P activity in vitro and rescued the growth of Escherichia coli and Saccharomyces cerevisiae strains with inactivations of their more complex and larger endogenous ribonucleoprotein RNase P. Homologs of Aquifex RNase P (HARP) were identified in many Archaea and some Bacteria, of which all Archaea and most Bacteria also encode an RNA-based RNase P; activity of both RNase P forms from the same bacterium or archaeon could be verified in two selected cases. Bioinformatic analyses suggest that A. aeolicus and related Aquificaceae likely acquired HARP by horizontal gene transfer from an archaeon.

Protein-only RNase P function in Escherichia coli: viability, processing defects and differences between PRORP isoenzymes

Abstract The RNase P family comprises structurally diverse endoribonucleases ranging from complex ribonucleoproteins to single polypeptides. We show that the organellar (AtPRORP1) and the two nuclear (AtPRORP2,3) single-polypeptide RNase P isoenzymes from Arabidopsis thaliana confer viability to Escherichia coli cells with a lethal knockdown of its endogenous RNA-based RNase P. RNA-Seq revealed that AtPRORP1, compared with bacterial RNase P or AtPRORP3, cleaves several precursor tRNAs (pre-tRNAs) aberrantly in E. coli. Aberrant cleavage by AtPRORP1 was mainly observed for pre-tRNAs that can form short acceptor-stem extensions involving G:C base pairs, including tRNAAsp(GUC), tRNASer(CGA) and tRNAHis. However, both AtPRORP1 and 3 were defective in processing of E. coli pre-tRNASec carrying an acceptor stem expanded by three G:C base pairs. Instead, pre-tRNASec was degraded, suggesting that tRNASec is dispensable for E. coli under laboratory conditions. AtPRORP1, 2 and 3 are also essentially unable to process the primary transcript of 4.5S RNA, a hairpin-like non-tRNA substrate processed by E. coli RNase P, indicating that PRORP enzymes have a narrower, more tRNA-centric substrate spectrum than bacterial RNA-based RNase P enzymes. The cells’ viability also suggests that the essential function of the signal recognition particle can be maintained with a 5΄-extended 4.5S RNA.