Establishing the precise evolutionary history of a gene improves prediction of disease-causing missense mutations

PURPOSE: Predicting the phenotypic effects of mutations has become an important application in clinical genetic diagnostics. Computational tools evaluate the behavior of the variant over evolutionary time and assume that variations seen during the course of evolution are probably benign in humans. However, current tools do not take into account orthologous/paralogous relationships. Paralogs have dramatically different roles in Mendelian diseases. For example, whereas inactivating mutations in the NPC1 gene cause the neurodegenerative disorder Niemann-Pick C, inactivating mutations in its paralog NPC1L1 are not disease-causing and, moreover, are implicated in protection from coronary heart disease. METHODS: We identified major events in NPC1 evolution and revealed and compared orthologs and paralogs of the human NPC1 gene through phylogenetic and protein sequence analyses. We predicted whether an amino acid substitution affects protein function by reducing the organism’s fitness. RESULTS: Removing the paralogs and distant homologs improved the overall performance of categorizing disease-causing and benign amino acid substitutions. CONCLUSION: The results show that a thorough evolutionary analysis followed by identification of orthologs improves the accuracy in predicting disease-causing missense mutations. We anticipate that this approach will be used as a reference in the interpretation of variants in other genetic diseases as well. Genet Med 18 10, 1029–1036

A thyroid hormone regulated asymmetric responsive centre is correlated with eye migration during flatfish metamorphosis

Flatfish metamorphosis is a unique post-embryonic developmental event in which thyroid hormones (THs) drive the development of symmetric pelagic larva into asymmetric benthic juveniles. One of the eyes migrates to join the other eye on the opposite side of the head. Developmental mechanisms at the basis of the acquisition of flatfish anatomical asymmetry remain an open question. Here we demonstrate that an TH responsive asymmetric centre, determined by deiodinase 2 expression, ventrally juxtaposed to the migrating eye in sole (Solea senegalensis) correlates with asymmetric cranial ossification that in turn drives eye migration. Besides skin pigmentation that is asymmetric between dorsal and ventral sides, only the most anterior head region delimited by the eyes becomes asymmetric whereas the remainder of the head and organs therein stay symmetric. Sub-ocular ossification is common to all flatfish analysed to date, so we propose that this newly discovered mechanism is universal and is associated with eye migration in all flatfish.Fundacao para a Ciencia e Tecnologia (FCT) [SFRH/BPD/66808/2009, IF/01274/2014]; FCT [SFRH/BPD/79105/2011, SFRH/BPD/89889/2012, PTDC/MAR/115005/2009, PEst-C/MAR/LA0015/2011, UID/Multi/04326/2013, Pest-OE/EQB/LA0023/2013, UID/BIM/04773/2013]; European Regional Development Fund through COMPETE; INIA; EU [RTA2013-00023-C02-01

Evolution of the mammalian lysozyme gene family

<p>Abstract</p> <p>Background</p> <p>Lysozyme <it>c </it>(chicken-type lysozyme) has an important role in host defense, and has been extensively studied as a model in molecular biology, enzymology, protein chemistry, and crystallography. Traditionally, lysozyme <it>c </it>has been considered to be part of a small family that includes genes for two other proteins, lactalbumin, which is found only in mammals, and calcium-binding lysozyme, which is found in only a few species of birds and mammals. More recently, additional testes-expressed members of this family have been identified in human and mouse, suggesting that the mammalian lysozyme gene family is larger than previously known.</p> <p>Results</p> <p>Here we characterize the extent and diversity of the lysozyme gene family in the genomes of phylogenetically diverse mammals, and show that this family contains at least eight different genes that likely duplicated prior to the diversification of extant mammals. These duplicated genes have largely been maintained, both in intron-exon structure and in genomic context, throughout mammalian evolution.</p> <p>Conclusions</p> <p>The mammalian lysozyme gene family is much larger than previously appreciated and consists of at least eight distinct genes scattered around the genome. Since the lysozyme <it>c </it>and lactalbumin proteins have acquired very different functions during evolution, it is likely that many of the other members of the lysozyme-like family will also have diverse and unexpected biological properties.</p

An Introductory Guide to Aligning Networks Using SANA, the Simulated Annealing Network Aligner.

Sequence alignment has had an enormous impact on our understanding of biology, evolution, and disease. The alignment of biological networks holds similar promise. Biological networks generally model interactions between biomolecules such as proteins, genes, metabolites, or mRNAs. There is strong evidence that the network topology-the "structure" of the network-is correlated with the functions performed, so that network topology can be used to help predict or understand function. However, unlike sequence comparison and alignment-which is an essentially solved problem-network comparison and alignment is an NP-complete problem for which heuristic algorithms must be used.Here we introduce SANA, the Simulated Annealing Network Aligner. SANA is one of many algorithms proposed for the arena of biological network alignment. In the context of global network alignment, SANA stands out for its speed, memory efficiency, ease-of-use, and flexibility in the arena of producing alignments between two or more networks. SANA produces better alignments in minutes on a laptop than most other algorithms can produce in hours or days of CPU time on large server-class machines. We walk the user through how to use SANA for several types of biomolecular networks

FLORA: a novel method to predict protein function from structure in diverse superfamilies

Predicting protein function from structure remains an active area of interest, particularly for the structural genomics initiatives where a substantial number of structures are initially solved with little or no functional characterisation. Although global structure comparison methods can be used to transfer functional annotations, the relationship between fold and function is complex, particularly in functionally diverse superfamilies that have evolved through different secondary structure embellishments to a common structural core. The majority of prediction algorithms employ local templates built on known or predicted functional residues. Here, we present a novel method (FLORA) that automatically generates structural motifs associated with different functional sub-families (FSGs) within functionally diverse domain superfamilies. Templates are created purely on the basis of their specificity for a given FSG, and the method makes no prior prediction of functional sites, nor assumes specific physico-chemical properties of residues. FLORA is able to accurately discriminate between homologous domains with different functions and substantially outperforms (a 2–3 fold increase in coverage at low error rates) popular structure comparison methods and a leading function prediction method. We benchmark FLORA on a large data set of enzyme superfamilies from all three major protein classes (α, β, αβ) and demonstrate the functional relevance of the motifs it identifies. We also provide novel predictions of enzymatic activity for a large number of structures solved by the Protein Structure Initiative. Overall, we show that FLORA is able to effectively detect functionally similar protein domain structures by purely using patterns of structural conservation of all residues

annot8r: GO, EC and KEGG annotation of EST datasets

<p>Abstract</p> <p>Background</p> <p>The expressed sequence tag (EST) methodology is an attractive option for the generation of sequence data for species for which no completely sequenced genome is available. The annotation and comparative analysis of such datasets poses a formidable challenge for research groups that do not have the bioinformatics infrastructure of major genome sequencing centres. Therefore, there is a need for user-friendly tools to facilitate the annotation of non-model species EST datasets with well-defined ontologies that enable meaningful cross-species comparisons. To address this, we have developed annot8r, a platform for the rapid annotation of EST datasets with GO-terms, EC-numbers and KEGG-pathways.</p> <p>Results</p> <p>annot8r automatically downloads all files relevant for the annotation process and generates a reference database that stores UniProt entries, their associated Gene Ontology (GO), Enzyme Commission (EC) and Kyoto Encyclopaedia of Genes and Genomes (KEGG) annotation and additional relevant data. For each of GO, EC and KEGG, annot8r extracts a specific sequence subset from the UniProt dataset based on the information stored in the reference database. These three subsets are then formatted for BLAST searches. The user provides the protein or nucleotide sequences to be annotated and annot8r runs BLAST searches against these three subsets. The BLAST results are parsed and the corresponding annotations retrieved from the reference database. The annotations are saved both as flat files and also in a relational postgreSQL results database to facilitate more advanced searches within the results. annot8r is integrated with the PartiGene suite of EST analysis tools.</p> <p>Conclusion</p> <p>annot8r is a tool that assigns GO, EC and KEGG annotations for data sets resulting from EST sequencing projects both rapidly and efficiently. The benefits of an underlying relational database, flexibility and the ease of use of the program make it ideally suited for non-model species EST-sequencing projects.</p

Application of COMPOCHIP Microarray to Investigate the Bacterial Communities of Different Composts

A microarray spotted with 369 different 16S rRNA gene probes specific to microorganisms involved in the degradation process of organic waste during composting was developed. The microarray was tested with pure cultures, and of the 30,258 individual probe-target hybridization reactions performed, there were only 188 false positive (0.62%) and 22 false negative signals (0.07%). Labeled target DNA was prepared by polymerase chain reaction amplification of 16S rRNA genes using a Cy5-labeled universal bacterial forward primer and a universal reverse primer. The COMPOCHIP microarray was applied to three different compost types (green compost, manure mix compost, and anaerobic digestate compost) of different maturity (2, 8, and 16 weeks), and differences in the microorganisms in the three compost types and maturity stages were observed. Multivariate analysis showed that the bacterial composition of the three composts was different at the beginning of the composting process and became more similar upon maturation. Certain probes (targeting Sphingobacterium, Actinomyces, Xylella/Xanthomonas/ Stenotrophomonas, Microbacterium, Verrucomicrobia, Planctomycetes, Low G + C and Alphaproteobacteria) were more influential in discriminating between different composts. Results from denaturing gradient gel electrophoresis supported those of microarray analysis. This study showed that the COMPOCHIP array is a suitable tool to study bacterial communities in composts

Anaerobic Carbon Monoxide Dehydrogenase Diversity in the Homoacetogenic Hindgut Microbial Communities of Lower Termites and the Wood Roach

Anaerobic carbon monoxide dehydrogenase (CODH) is a key enzyme in the Wood-Ljungdahl (acetyl-CoA) pathway for acetogenesis performed by homoacetogenic bacteria. Acetate generated by gut bacteria via the acetyl-CoA pathway provides considerable nutrition to wood-feeding dictyopteran insects making CODH important to the obligate mutualism occurring between termites and their hindgut microbiota. To investigate CODH diversity in insect gut communities, we developed the first degenerate primers designed to amplify cooS genes, which encode the catalytic (β) subunit of anaerobic CODH enzyme complexes. These primers target over 68 million combinations of potential forward and reverse cooS primer-binding sequences. We used the primers to identify cooS genes in bacterial isolates from the hindgut of a phylogenetically lower termite and to sample cooS diversity present in a variety of insect hindgut microbial communities including those of three phylogenetically-lower termites, Zootermopsis nevadensis, Reticulitermes hesperus, and Incisitermes minor, a wood-feeding cockroach, Cryptocercus punctulatus, and an omnivorous cockroach, Periplaneta americana. In total, we sequenced and analyzed 151 different cooS genes. These genes encode proteins that group within one of three highly divergent CODH phylogenetic clades. Each insect gut community contained CODH variants from all three of these clades. The patterns of CODH diversity in these communities likely reflect differences in enzyme or physiological function, and suggest that a diversity of microbial species participate in homoacetogenesis in these communities

Increased retention of functional fusions to toxic genes in new two-hybrid libraries of the E. coli strain MG1655 and B. subtilis strain 168 genomes, prepared without passaging through E. coli

BACKGROUND: Cloning of genes in expression libraries, such as the yeast two-hybrid system (Y2H), is based on the assumption that the loss of target genes is minimal, or at worst, managable. However, the expression of genes or gene fragments that are capable of interacting with E. coli or yeast gene products in these systems has been shown to be growth inhibitory, and therefore these clones are underrepresented (or completely lost) in the amplified library. RESULTS: Analysis of candidate genes as Y2H fusion constructs has shown that, while stable in E. coli and yeast for genetic studies, they are rapidly lost in growth conditions for genomic libraries. This includes the rapid loss of a fragment of the E. coli cell division gene ftsZ which encodes the binding site for ZipA and FtsA. Expression of this clone causes slower growth in E. coli. This clone is also rapidly lost in yeast, when expressed from a GAL1 promoter, relative to a vector control, but is stable when the promoter is repressed. We have demonstrated in this report that the construction of libraries for the E. coli and B. subtilis genomes without passaging through E. coli is practical, but the number of transformants is less than for libraries cloned using E. coli as a host. Analysis of several clones in the libraries that are strongly growth inhibitory in E. coli include genes for many essential cellular processes, such as transcription, translation, cell division, and transport. CONCLUSION: Expression of Y2H clones capable of interacting with E. coli and yeast targets are rapidly lost, causing a loss of complexity. The strategy for preparing Y2H libraries described here allows the retention of genes that are toxic when inappropriately expressed in E. coli, or yeast, including many genes that represent potential antibacterial targets. While these methods are generally applicable to the generation of Y2H libraries from any source, including mammalian and plant genomes, the potential of functional clones interacting with host proteins to inhibit growth would make this approach most relevant for the study of prokaryotic genomes