Evolution of the multidomain protein wheat germ agglutinin

Summary We compared the homologous amino acid sequences of hevein and each of the four domains (A, B, C, and D) of wheat germ agglutinin and used them to construct a pseudophylogenetic tree relating these sequences to a hypothetical common ancestor sequence. In the crystal structure of the wheat germ agglutinin dimer, six pseudo-twofold rotational symmetry axes have previously been located in addition to the true twofold axis. Four of these relate two nonidentical domains to each other in each of the four possible pairs constituting the sugar-binding sites (A₁D₂, A₂D₁, B₁C₂, and B₂C₁). The remaining two relate contiguous unique pairs of sugar-binding sites to each other (A₁D₂ to B₁C₂, and A₂D₁ to B₂C₁). These latter two sets of pairs are related to each other by the true twofold axis. Side chains that mediate sugar binding in the interfaces of each of the four pairs were found to be largely conserved. The sequence homology, taken together with these pseudo-symmetry elements in the dimer structure, suggests a pathway for the evolution of the four-domain molecule from a single-domain dimer that can be correlated with simultaneous development of the saccharide-binding sites.     

Molecular analysis of the split cox1 gene from the Basidiomycota Agrocybe aegerita: relationship of its introns with homologous Ascomycota introns and divergence levels from common ancestral copies

The Basidiomycota Agrocybe aegerita (Aa) mitochondrial cox1 gene (6790 nucleotides), encoding a protein of 527 aa (58 377 Da), is split by four large subgroup IB introns possessing site-specific endonucleases assumed to be involved in intron mobility. When compared to other fungal COX1 proteins, the Aa protein is closely related to the COX1 one of the Basidiomycota Schizophyllum commune (Sc). This clade reveals a relationship with the studied Ascomycota ones, with the exception of Schizosaccharomyces pombe (Sp) which ranges in an out-group position compared with both higher fungi divisions. When comparison is extended to other kingdoms, fungal COX1 sequences are found to be more related to algae and plant ones (more than 57.5% aa similarity) than to animal sequences (53.6% aa similarity), contrasting with the previously established close relationship between fungi and animals, based on comparisons of nuclear genes. The four Aa cox1 introns are homologous to Ascomycota or algae cox1 introns sharing the same location within the exonic sequences. The percentages of identity of the intronic nucleotide sequences suggest a possible acquisition by lateral transfers of ancestral copies or of their derived sequences. These identities extend over the whole intronic sequences, arguing in favor of a transfer of the complete intron rather than a transfer limited to the encoded ORF. The intron i4 shares 74% of identity, at the nucleotidic level, with the Podospora anserina (Pa) intron i14, and up to 90.5% of aa similarity between the encoded proteins, i.e. the highest values reported to date between introns of two phylogenetically distant species. This low divergence argues for a recent lateral transfer between the two species. On the contrary, the low sequence identities (below 36%) observed between Aa i1 and the homologous Sp i1 or Prototheca wickeramii (Pw) i1 suggest a long evolution time after the separation of these sequences. The introns i2 and i3 possessed intermediate percentages of identity with their homologous Ascomycota introns. This is the first report of the complete nucleotide sequence and molecular organization of a mitochondrial cox1 gene of any member of the Basidiomycota division.     

A maximum-likelihood approach to analyzing nonoverlapping and overlapping reading frames

A model is presented for sequence evolution on the basis of which one can analyze combinations of noncoding, singly coding, and multiply coding regions of aligned homologous DNA sequences. It is a generalization of Kimura's (J. Mol. Evol. 16:111–120, 1980) and Li et al.'s (J. Mol. Evol. 36:96–99, 1985) transition-transversion models with selection on replacement substitutions.Based on a hierarchy of hypotheses, one will be able to estimate selection factors and transition and transversion distances for different combinations of regions ranging from many regions, each with their private set of parameters, to one set of parameters for all regions.The method is demonstrated on two aligned HIV I retroviruses. Correspondence to: J. Hein     

In silico investigation of PHD-3 specific HIF1-α proline 567 hydroxylation: A new player in the VHL/HIF-1α interaction pathway?

Hypoxia inducible factor 1α (HIF-1α) regulates oxygen homeostasis in the cell through a sensing mechanism involving its hydroxylation and binding to the von Hippel–Lindau (VHL) tumor suppressor. This mechanism is mediated through hydroxylation of HIF-1α proline 564, although in vitro tests have previously shown an alternative hydroxylation at proline 567 by PHD-3. Here, molecular dynamics simulations were used to investigate the structural effect of this alternative hydroxylation. A specific hydrogen bond network rearrangement and improved electrostatic energy for hydroxylated P567 are compatible with an increase in HIF-1α binding affinity. Sequence analysis also confirms P567 to be vastly conserved during evolution, indicating a possible role for this alternative, PHD-3 driven, post translational modification in pVHL–HIF-1α complex formation.     

Classification of fungal and bacterial lytic polysaccharide monooxygenases

Background

Lytic polysaccharide monooxygenases are important enzymes for the decomposition of recalcitrant biological macromolecules such as plant cell wall and chitin polymers. These enzymes were originally designated glycoside hydrolase family 61 and carbohydrate-binding module family 33 but are now classified as auxiliary activities 9, 10 and 11 in the CAZy database. To obtain a systematic analysis of the divergent families of lytic polysaccharide monooxygenases we used Peptide Pattern Recognition to divide 5396 protein sequences resembling enzymes from families AA9 (1828 proteins), AA10 (2799 proteins) and AA11 (769 proteins) into subfamilies.

Results

The results showed that the lytic polysaccharide monooxygenases have two conserved regions identified by conserved peptides specific for each AA family. The peptides were used for in silico PCR discovery of the lytic polysaccharide monooxygenases in 79 fungal and 95 bacterial genomes. The bacterial genomes encoded 0 – 7 AA10s (average 0.6). No AA9 or AA11 were found in the bacteria. The fungal genomes encoded 0 – 40 AA9s (average 7) and 0 – 15 AA11s (average 2) and two of the fungi possessed a gene encoding a putative AA10. The AA9s were mainly found in plant cell wall-degrading asco- and basidiomycetes in agreement with the described role of AA9 enzymes. In contrast, the AA11 proteins were found in 36 of the 39 ascomycetes and in only two of the 32 basidiomycetes and their abundance did not correlate to the degradation of cellulose and hemicellulose.

Conclusions

These results provides an overview of the sequence characteristics and occurrence of the divergent AA9, AA10 and AA11 families and pave the way for systematic investigations of the of lytic polysaccharide monooxygenases and for structure-function studies of these enzymes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1601-6) contains supplementary material, which is available to authorized users.     

Novel Y-chromosome Short Tandem Repeat Variants Detected Through the Use of Massively Parallel Sequencing

Massively parallel sequencing(MPS) technology is capable of determining the sizes of short tandem repeat(STR) alleles as well as their individual nucleotide sequences. Thus, single nucleotide polymorphisms(SNPs) within the repeat regions of STRs and variations in the pattern of repeat units in a given repeat motif can be used to differentiate alleles of the same length. In this study, MPS was used to sequence 28 forensically-relevant Y-chromosome STRs in a set of 41 DNA samples from the 3 major U.S. population groups(African Americans, Caucasians, and Hispanics).The resulting sequence data, which were analyzed with STRait Razor v2.0, revealed 37 unique allele sequence variants that have not been previously reported. Of these, 19 sequences were variations of documented sequences resulting from the presence of intra-repeat SNPs or alternative repeat unit patterns. Despite a limited sampling, two of the most frequently-observed variants were found only in African American samples. The remaining 18 variants represented allele sequences for which there were no published data with which to compare. These findings illustrate the great potential of MPS with regard to increasing the resolving power of STR typing and emphasize the need for sample population characterization of STR alleles.     

Sequence Analysis and Domain Motifs in the Porcine Skin Decorin Glycosaminoglycan Chain

Decorin proteoglycan is comprised of a core protein containing a single O-linked dermatan sulfate/chondroitin sulfate glycosaminoglycan (GAG) chain. Although the sequence of the decorin core protein is determined by the gene encoding its structure, the structure of its GAG chain is determined in the Golgi. The recent application of modern MS to bikunin, a far simpler chondroitin sulfate proteoglycans, suggests that it has a single or small number of defined sequences. On this basis, a similar approach to sequence the decorin of porcine skin much larger and more structurally complex dermatan sulfate/chondroitin sulfate GAG chain was undertaken. This approach resulted in information on the consistency/variability of its linkage region at the reducing end of the GAG chain, its iduronic acid-rich domain, glucuronic acid-rich domain, and non-reducing end. A general motif for the porcine skin decorin GAG chain was established. A single small decorin GAG chain was sequenced using MS/MS analysis. The data obtained in the study suggest that the decorin GAG chain has a small or a limited number of sequences.     

Sequence analysis and phylogenetic study of some toxin proteins of snakes and related non-toxin proteins of chordates

Snakes are equipped with their venomic armory to tackle different prey and predators in adverse natural world. The venomic composition of snakes is a mix of biologically active proteins and polypeptides. Among different components snake venom cytotoxins and short neurotoxin are non-enzymatic polypeptide candidates with in the venom. These two components structurally resembled to three-finger protein superfamily specific scaffold. Different non-toxin family members of three-finger protein superfamily are involved in different biological roles. In the present study we analyzed the snake venom cytotoxins, short neurotoxins and related non-toxin proteins of different chordates in terms of amino acid sequence level diversification profile, polarity profile of amino acid sequences, conserved pattern of amino acids and phylogenetic relationship of these toxin and nontoxin protein sequences. Sequence alignment analysis demonstrates the polarity specific molecular enrichment strategy for better system adaptivity. Occurrence of amino acid substitution is high in number in toxin sequences. In non-toxin body proteins there are less amino acid substitutions. With the help of conserved residues these proteins maintain the three-finger protein scaffold. Due to system specific adaptation toxin and non-toxin proteins exhibit a varied type of amino acid residue distribution in sequence stretch. Understanding of Natural invention scheme (recruitment of venom proteins from normal body proteins) may help us to develop futuristic engineered bio-molecules with remedial properties.     

Phylogenetic analysis of cubilin (CUBN) gene

Cubilin, (CUBN; also known as intrinsic factor-cobalamin receptor [Homo sapiens Entrez Pubmed ref {"type":"entrez-nucleotide","attrs":{"text":"NM_001081.3","term_id":"126091151","term_text":"NM_001081.3"}}NM_001081.3; {"type":"entrez-nucleotide","attrs":{"text":"NG_008967.1","term_id":"212549718","term_text":"NG_008967.1"}}NG_008967.1; GI: 119606627]), located in the epithelium of intestine and kidney acts as a receptor for intrinsic factor – vitamin B12 complexes. Mutations in CUBN may play a role in autosomal recessive megaloblastic anemia. The current study investigated the possible role of CUBN in evolution using phylogenetic testing. A total of 588 BLAST hits were found for the cubilin query sequence and these hits showed putative conserved domain, CUB superfamily (as on 27^th Nov 2012). A first-pass phylogenetic tree was constructed to identify the taxa which most often contained the CUBN sequences. Following this, we narrowed down the search by manually deleting sequences which were not CUBN. A repeat phylogenetic analysis of 25 taxa was performed using PhyML, RAxML and TreeDyn softwares to confirm that CUBN is a conserved protein emphasizing its importance as an extracellular domain and being present in proteins mostly known to be involved in development in many chordate taxa but not found in prokaryotes, plants and yeast.. No horizontal gene transfers have been found between different taxa.     

The impact of post-alignment processing procedures on whole-exome sequencing data

The use of post-alignment procedures has been suggested to prevent the identification of false-positives in massive DNA sequencing data. Insertions and deletions are most likely to be misinterpreted by variant calling algorithms. Using known genetic variants as references for post-processing pipelines can minimize mismatches. They allow reads to be correctly realigned and recalibrated, resulting in more parsimonious variant calling. In this work, we aim to investigate the impact of using different sets of common variants as references to facilitate variant calling from whole-exome sequencing data. We selected reference variants from common insertions and deletions available within the 1K Genomes project data and from databases from the Latin American Database of Genetic Variation (LatinGen). We used the Genome Analysis Toolkit to perform post-processing procedures like local realignment, quality recalibration procedures, and variant calling in whole exome samples. We identified an increased number of variants from the call set for all groups when no post-processing procedure was performed. We found that there was a higher concordance rate between variants called using 1K Genomes and LatinGen. Therefore, we believe that the increased number of rare variants identified in the analysis without realignment or quality recalibration indicated that they were likely false-positives.