Amino Acid Sequence Alignment Research Articles

A Novel Major Pilin Subunit Protein FimM Is Involved in Adhesion of Bifidobacterium longum BBMN68 to Intestinal Epithelial Cells.

Adhesion to the gastrointestinal tract is considered to be important for bifidobacteria to colonize the human gut and exert their probiotic effects. Some cell surface proteins of bifidobacteria, known as adhesins, play critical roles in the binding to host cells or the extracellular matrix (ECM). To elucidate the mechanisms associated with the adhesion of Bifidobacterium longum BBMN68, a centenarian originated potential probiotic, PSORTdb was employed to identify putative extracellular localized proteins in the B. longum BBMN68. Of the 560 predicted extracellular proteins, 21 were further identified as putative adhesion proteins using the conserved domain database of NCBI, and four were successfully overexpressed in the heterologous host, Lactococcus lactis NZ9000. Notably, a recombinant strain expressing FimM showed a significantly increased adhesive affinity for both HT-29 and mucus-secreting LS174T goblet cells (2.2- and 5.4-fold higher than that of the control strain, respectively). Amino acid sequence alignment showed that FimM is a major pilin subunit protein containing a Cna-B type domain and a C-terminal LPKTG sequence. However, in silico analysis of the fimM-coding cluster revealed that BBMN68_RS10200, encoding a pilus-specific class C sortase, was a pseudogene, indicating that FimM may function as a surface adhesin that cannot polymerize into a pili-like structure. Immunogold electron microscopy results further confirmed that FimM localized to the surface of L. lactis NZfimM and B. longum BBMN68 but did not assemble into pilus filaments. Moreover, the adhesive affinity of L. lactis NZfimM to fibronectin, fibrinogen, and mucin were 3.8-, 2.1-, and 3.1-fold higher than that of the control. The affinity of FimM for its attachment receptors was further verified through an inhibition assay using anti-FimM antibodies. In addition, homologs of FimM were found in Bifidobacterium bifidum 85B, Bifidobacterium gallinarum CACC 514, and 23 other B. longum strains by sequence similarity analysis using BLASTP. Our results suggested that FimM is a novel surface adhesin that is mainly present in B. longum strains.

Caloric restriction upregulates the expression of DNMT3.1, lacking the conserved catalytic domain, in Daphnia magna.

DNA methylation plays an important role in many aspects of biology, including development, disease, and phenotypic plasticity. In the branchiopod crustacean, Daphnia, de novo DNA methylation has been detected in specific environmental contexts. However, fundamental information on de novo DNA methyltransferase DNMT3 orthologs, including domain organization, developmental expression, and response to environmental stimuli, is lacking. In this study, we examined two DNMT3 orthologs in Daphnia magna, DapmaDNMT3.1 and DapmaDNMT3.2. Amino acid sequence alignment revealed that DapmaDNMT3.1 and DapmaDNMT3.2 lack the conserved methyltransferase motifs of the catalytic domain and the PWWP domain, respectively. We profiled the expression of the two orthologs during embryogenesis and under various feeding levels. During embryogenesis, in contrast to the low DapmaDNMT3.1 expression, DapmaDNTM3.2 was highly expressed at specific stages, that is, in the one cell-stage and at 48 hr post ovulation. In nutrient-rich condition, both genes were lowly expressed, whereas DapmaDNMT3.1 was upregulated at the lower food levels, suggesting a potential role of DapmaDNMT3.1 in gene regulation in response to caloric restriction. These findings provide a basis for understanding the developmental stage- and stress-dependent function of DNMT3 orthologs in D. magna.

Molecular identification and functional analysis of Niemann-Pick type C2 protein in Macrocentrus cingulum Brischke (Hymenoptera: Braconidae)

• Niemann-Pick type C2 gene (McinNPC2) was cloned from Macrocentrus cingulum . • McinNPC2 is likely to be OBPs and the nearest evolutionary relation genes as OBPs were clustered. • McinNPC2 may be related protein as a new carrier of aroma compounds. • McinNPC2 may play roles in perception of plant volatiles. Niemann-Pick type C2 (NPC2) proteins in arthropods have been extensively differentiated and possibly duplicated according to environmental conditions and are probable to have different functions. The participation of NPC2 proteins in chemical communication in arthropods brings new objectives in environmental-friendly strategies for pest population control. In this study, NPC2 gene in Macrocentrus cingulum (McinNPC2) was newly identified by rapid amplification cDNA ends (RACE) technology. McinNPC2 amino acid sequence alignment with other representative NPC2 annotates to evaluate the highly conserved consensus amino acids, but with odorant binding proteins in M. cingulum show that only one consensus amino acid. Primary six-cysteine structures that are same to odorant binding proteins in M. cingulum were observed in McinNPC2. Phylogenetic analysis of McinNPC2 indicated that the nearest monophyletic group forming one clade with high posterior probability values clusters as Cyphomyrmex costatus (CcosNPC2) whereas the nearest evolutionary relation group as some odorant binding proteins. Moreover, quantitative real-time PCR (qPCR) measurements show that the McinNPC2 gene expression level in various tissues of the female is significantly and ubiquitously higher than in male, whereas the highest expression level in female antennae. We further explore the binding characterization of recombinant McinNPC2 to candidate odor molecules and did the modeling and docking simulations. The results showed ligands binding specificity and docking tests results indicate that β-ionone, an aroma compound commonly found in essential oils, can strongly bind with McinNPC2 . In conclusion, we proposed that McinNPC2 may be involved in chemical communication and play roles in perception of plant volatiles.

A sequence embedding method for enzyme optimal condition analysis

BackgroundAn enzyme activity is influenced by the external environment. It is important to have an enzyme remain high activity in a specific condition. A usual way is to first determine the optimal condition of an enzyme by either the gradient test or by tertiary structure, and then to use protein engineering to mutate a wild type enzyme for a higher activity in an expected condition.ResultsIn this paper, we investigate the optimal condition of an enzyme by directly analyzing the sequence. We propose an embedding method to represent the amino acids and the structural information as vectors in the latent space. These vectors contain information about the correlations between amino acids and sites in the aligned amino acid sequences, as well as the correlation with the optimal condition. We crawled and processed the amino acid sequences in the glycoside hydrolase GH11 family, and got 125 amino acid sequences with optimal pH condition. We used probabilistic approximation method to implement the embedding learning method on these samples. Based on these embedding vectors, we design a computational score to determine which one has a better optimal condition for two given amino acid sequences and achieves the accuracy 80% on the test proteins in the same family. We also give the mutation suggestion such that it has a higher activity in an expected environment, which is consistent with the previously professional wet experiments and analysis.ConclusionA new computational method is proposed for the sequence based on the enzyme optimal condition analysis. Compared with the traditional process that involves a lot of wet experiments and requires multiple mutations, this method can give recommendations on the direction and location of amino acid substitution with reference significance for an expected condition in an efficient and effective way.

Characterization of Cu/Zn-SODs in sympatric species: A comparison of zebra and quagga mussels

Zebra mussels (Dreissena polymorpha) are commonly used as sentinel organisms in freshwater biomonitoring. The displacement of this species by quagga mussels (Dreissena rostriformis bugensis) has raised the issue of their discrimination after field sampling. Even if those two species exhibit specific morphological characteristics, confusion can occur. However, behind the visible phenotypes lie important physiological and metabolic differences that may affect the interpretation of biomarkers in risk assessment. In this study, we present a characterization of Cu/Zn-SODs (superoxide dismutases) in Dreissena spp. with a comparison of activity profiles from individuals originating from North America and Europe. A phylogenetic tree was constructed and dreissenid SOD sequences appear to be closely related to clam species sequences. A specific isoform at an isoelectric point value of 4.95 was observed only for quagga mussels independently of their origin, while zebra mussels exhibited several polymorphic profiles, even for a single sampling site. A comparison of transcript sequences was conducted and the observed differences were highlighted. Alignment of amino acid sequences shows a high level of identity between the two species associated with minor changes in the predicted secondary structures.

In silico analysis and molecular characterization of Influenza A (H1N1) pdm09 virus circulating and causing major outbreaks in central India, 2009-2019.

Background and Objectives:Influenza A/H1N1pdm09 causes respiratory illness and remains a concern for public health. Since its first emergence in 2009, the virus has been continuously circulating in the form of its genetic variants. Influenza A/H1N1pdm09 surveillance is essential for uncovering emerging variants of epidemiologic and vaccine efficacy. The present study attempts in silico analysis and molecular characterization of Influenza A (H1N1) pdm09 virus circulating and causing major outbreaks in central India during 2009–2019.Materials and Methods:We have investigated the antigenic drift analysis of 96 isolates’ hemagglutinin (HA) gene sequences (59 central Indian and 37 local Indian and 28 global reference HA gene sequences) of Influenza A/H1N1pdm09 viruses from 2009 to 2019. The study includes mutational (Multiple sequence Alignment), phylogenetic (Maximum Likelihood Method), and statistical analysis (Covariance and correlation) of HA sequences submitted in NCBI, IRD and GISAID from central India.Results:Phylogenetic analysis indicated maximum clustering of central Indian HA gene sequences in genogroup 6B. Analysis of amino acid sequence alignment revealed changes in receptor binding site (RBS). The frequency of S220T amino acid substitution was found to be high followed by S202T, K300E A273T, K180Q. The Karl Pearson correlation coefficient (r) and covariance between the number of mutations and the death toll was found +0.246 and +100.3 respectively.Conclusion:The study identifies the continuous genetic variations in the HA gene sequences of circulating Influenza A/H1N1pdm09 in central India from the year 2009 to 2019. Further suggesting importance of monitoring the gradual evolution of the virus with regards to an increase in virulence, pathogenicity and vaccine efficacy timely.

Genome analysis and phylogenetic characterization of two deformed wing virus strains from Apis cerana in Vietnam.

BackgroundDeformed wing virus (DWV) is a virulent virus that causes honeybee disease. DWV can exist as a latent infection in honeybees, outbreak into epidemics, and cause serious damage to beekeeping cross the world, including Vietnam.MethodsThe two DWV strains circulating in Vietnamese honeybee, Apis cerana, were first isolated from adult honeybees in North Vietnam (DWV-NVN) and South Vietnam (DWV-SVN). Their complete nucleotide sequences were determined, aligned, and compared with other DWV strains.ResultsThe two Vietnamese DWV strains comprised 10,113 bp and contained a large single open reading frame (ORF) of 2,893 amino acids, initiating at nucleotide 1,130 and terminating at nucleotide 9,812. Multiple nucleotide sequence alignment between these two DWV-VN strains and DWV strains in A. mellifera was performed. The DWV-VN strains showed a low genetic identity (from 91.4% to 92.0%) with almost of these strains, but lower identities (89.2% and 89.4%) with UK2 and (89.6%) with the China2 strain. Low identities (91.7% and 91.9%) were also observed between the China3 strain (in A. cerana) and the DWV-VN strains, respectively. The deduced amino acid sequence alignment showed high genetic similarities (97.0%–97.9%) when the USA1, Chile, Italy1, France, UK1, UK2, Japan, Korea2, China1, China2 and China3 strains were compared to the DWV-VN strains. This ratio was 96.7% and 96.8% when the Korea1 strain was compared to the DWV-SVN and DWV-NVN strains, respectively. Numerous amino acid substitutions were identified in the L, VP3, and RdRp sequences. Notably, we observed six substitutions positioned at amino acids 27 (E > I), 98 (S > T), 120 (A > V), 153 (M > T), 170 (D > F), and 174 (Y > F) in the L protein, two amino acid changes at positions 980 (S > A) and 1032 (E > T) in VP3, and one amino acid change at position 2627 (R > C) unique to the DWV-VN strains. Phylogenetic analysis based on complete genome sequences, RdRp sequences and Simplot analysis indicated that there was a significant difference between DWV-VN strains in A. cerana and DWV strains in A. mellifera. The results suggested that the genetic variations of the DWV-VN strains in A. cerana help them to adapt geographical conditions and may lead to change the viral pathogenicity of DWV-VN strains.

Genetic characterization of structural and open reading Fram-8 proteins of SARS-CoV-2 isolates from different countries

Since December 2019, a severe pandemic of pneumonia, COVID-19 associated with a novel coronavirus (SARS-CoV-2), have emerged in Wuhan, China and spreading throughout the world. As RNA viruses have a high mutation rate therefore we wanted to identify whether this virus is also prone to mutations. For this reason we selected four major structural (Spike protein (S), Envelope protein (E), Membrane glycoprotein (M), Nucleocapsid phosphoprotein (N)) and ORF8 protein of 100 different SARS-CoV-2 isolates of fifteen countries from NCBI database and compared these to the reference sequence, Wuhan NC_045512.2, which was the first isolate of SARS-CoV-2 that was sequenced. By multiple sequence alignment of amino acids, we observed substitutions and deletion in S protein at 13 different sites in the isolates of five countries (China, USA, Finland, India and Australia) as compared to the reference sequence. Similarly, alignment of N protein revealed substitutions at three different sites in isolates of China, Spain and Japan. M protein exhibits substitution only in one isolates from USA, however, no mutation was observed in E protein of any isolate. Interestingly, in ORF8 substitution of Leucine, a nonpolar to Serine a polar amino acid at same position (aa84 L to S) in 23 isolates of five countries i.e. China, USA, Spain, Taiwan and India were observed, which may affect the conformation of peptides. Thus, we observed several mutations in the isolates thereafter the first sequencing of SARS-CoV-2 isolate, NC_045512.2, which suggested that this virus might be a threat to the whole world and therefore further studies are needed to characterize how these mutations in different proteins affect the functionality and pathogenesis of SARS-CoV-2.

LvCSN5 is involved in WSSV infection via interaction with wsv006

As an extremely virulent pathogen, white spot syndrome virus (WSSV) greatly threatens shrimp aquaculture worldwide. The interaction between virus and host is important for viral infection. In the present study, a yeast two-hybrid (Y2H) library was constructed to clarify the functions of wsv006, and the interaction between wsv006 and shrimp Litopenaeus vannamei (L. vannamei) was analyzed. Furthermore, we explored the role of the wsv006-interacting molecule L. vannamei COP9 constitutive photomorphogenic-like protein subunit 5 (LvCSN5) in WSSV infection. Y2H assay showed that wsv006 interacted with LvCSN5, and co-immunoprecipitation (Co-IP) assay confirmed such interaction. Multiple alignments of amino acid sequences with other species revealed that the LvCSN5 had high identity with Penaeusmonodon CSN5 (PmCSN5). LvCSN5 was mainly expressed in intestine, eye and hepatopancreas. In addition, the relative expression of LvCSN5 was significantly up-regulated both in intestine and hepatopancreas following the WSSV challenge. Besides, the relative expressions of IE1 and VP28, as well as the viral copy numbers were significantly increased in the LvCSN5-silenced shrimp. Our findings suggested that LvCSN5 was involved in WSSV infection by interacting with wsv006.

Identification of 22 Novel Motifs of the Cell Entry Fusion Glycoprotein B of Oncolytic Herpes Simplex Viruses: Sequence Analysis and Literature Review

Objective: Herpes simplex viruses (HSVs) are widely spread throughout the world, causing infections from oral, and genital mucous membrane ulcerations to severe viral encephalitis. Glycoprotein B (gB) was the first HSV envelope glycoprotein identified to induce cell fusion. This glycoprotein initiates viral entry and thereby determines the infectivity of HSV, as well as oncolytic HSV (oHSV). Clarifying its molecular characterization and enlarging its motif reservoir will help to engineer oHSV and in cancer treatment applications. Only in recent years has the importance of gB been acknowledged in HSV infection and oHSV engineering. Although gB-modified oHSVs have been developed, the detailed molecular biology of gB needs to be illustrated more clearly in order to construct more effective oHSVs.Method: Here, we performed a systematic comparative sequence analysis of gBs from the 9 HSV-1 and 2 HSV-2 strains, including HSV-1-LXMW, which was isolated by our lab. Online software was implemented to predict gB secondary structure and motifs. Based on extensive literature reviews, a functional analysis of the predicted motifs was performed.Results: Here, we reported the DNA and predicted amino acid sequences of our recently isolated HSV-1-LXMW and found that the strain was evolutionarily close to HSV-1 strains F, H129, and SC16 based on gB analysis. The 22 novel motifs of HSV gB were identified for the first time. An amino acid sequence alignment of the 11 HSV strains showed that the gB motifs are conserved among HSV strains, suggesting that they are functional in vivo. Additionally, we found that certain amino acids within the 13 motifs out of the 22 were reported to be functional in vivo. Furthermore, the gB mutants and gB-engineered oHSVs were also summarized.Conclusion: Our identification of the 22 novel motifs shed light on HSV gB biology and provide new options for gB engineering to improve the efficiency and safety of oHSVs.

Unveiling the immunomodulatory properties of Haemonchus contortus adhesion regulating molecule 1 interacting with goat T cells

BackgroundGastrointestinal nematodes could release excretory-secretory (ES) proteins into the host environment to ensure their survival. These ES proteins act as immunomodulators to suppress or subvert the host immune response via the impairment of immune cell functions, especially in chronic infections. In our preliminary study, Haemonchus contortus adhesion-regulating molecule 1 (HcADRM1) was identified from H. contortus ES proteins (HcESPs) that interacted with host T cells via liquid chromatography-tandem mass spectrometry analysis. However, little is known about HcADRM1 as an ES protein which may play a pivotal role at the parasite-host interface.MethodsBased on bioinformatics approaches, multiple amino acid sequence alignment was conducted and the evolutionary relationship of HcADRM1 with ADRM1 orthologues was extrapolated. Employing RT-qPCR and immunohistochemistry assays, temporal transcriptional and spatial expression profiles of HcADRM1 were investigated. Using immunostaining approaches integrated with immunological bioassays, the immunomodulatory potentials of HcADRM1 on goat T cells were assessed.ResultsWe hereby demonstrated that HcADRM1 with immunodiagnostic utility was a mammalian ADRM1 orthologue abundantly expressed at all developmental stages of H. contortus. Given the implications of ADRM1 proteins in cell growth, survival and development, we further investigated the immunomodulatory property of HcADRM1 as an individual ES protein acting at the parasite-host interface. The rHcADRM1 stimuli notably suppressed T cell viability, promoted intrinsic and extrinsic T cell apoptosis, inhibited T cell proliferation and induced cell cycle arrest at G1 phase. Simultaneously, rHcADRM1 stimuli exerted critical controls on T cell cytokine secretion profiles, predominantly by restraining the secretions of interleukin (IL)-4, IL-10 and interferon-gamma.ConclusionsImportantly, HcADRM1 protein may have prophylactic potential for anti-H. contortus vaccine development. Together, these findings may contribute to the clarification of molecular and immunomodulatory traits of ES proteins, as well as improvement of our understanding of parasite immune evasion mechanism in H. contortus-host biology.

Phylogenetic Analysis of South African Bovine Leukaemia Virus (BLV) Isolates.

Bovine leukaemia virus (BLV) causes chronic lymphoproliferative disorder and fatal lymphosarcoma in cattle, leading to significant economic losses in the beef and dairy industries. BLV is endemic globally and eleven genotypes have been identified. To date, only Zambian isolates have been genotyped from Africa. Although high BLV prevalence has been reported in South Africa, there has been no molecular characterisation of South African BLV isolates. To characterise BLV isolates in South Africa for the first time, we investigated the phylogenetic relationships and compared the genetic variability of eight South African BLV isolates with BLV isolates representing the eleven known genotypes from different geographical regions worldwide. Phylogenetic analyses based on full-length and partial env sequences as well as full-length gag sequences revealed that at least two genotypes, genotypes 1 (G1) and 4 (G4), are present in cattle in South Africa, which is consistent with studies from Zambia. However, our analysis revealed that the G1 South African isolate is more similar to other G1 isolates than the G1 Zambian isolates whereas, the G4 South African isolates are more divergent from other G4 isolates but closely related to the G4 Zambian isolate. Lastly, amino acid sequence alignment identified genotype-specific as well as novel amino acid substitutions in the South African isolates. The detection of two genotypes (G1 and G4) in southern Africa highlights the urgent need for disease management and the development of an efficacious vaccine against local strains.

Analysis on genetic characteristics of H9N2 avian influenza virus isolated from human infection and external environment in Gansu province

Objective: To analyze the genomic characteristics of human infection with H9N2 avian influenza virus in Gansu province. Methods: The etiological analysis was conducted for human infection with H9N2 avian influenza virus detected in influenza like illness cases in northwestern China in 2016. Molecular bioinformatics Mega 7.0 software was used to analyze the full genomic sequences of the viral isolate. Results: The gene fragments of HA, NA, MP, NP, NS, PA, PB1 and PB2 of the isolate were highly similar (>90%) to those of H9N2 avian influenza virus strain isolated in external environment in Gansu from 2014 to 2019. The HA gene belonged to BJ/94-like branch, PB2 and MP belonged to G1/97-like branch, and the PB1, PA, NS, and NP genes belonged to F/98-like branch. MP and PB2 were closely related to H7N9, H10N8 and H5N6 viruses. Amino acid sequence alignment showed that the HA cleavage site was arranged in PSRSSR ↓ GLF, H183N and Q226L mutated which included 7 HA glycosylated sites; 62-64 sites of NA absented 3 amino acids (ITE); and M2-31N, NS1-42S, PA-356R, and PA-409N mutated. Conclusions: Apparently, this case of human infection with human infection with H9N2 avian influenza virus was an incidental. However, the isolates of H9N2 influenza virus in external environment of Gansu had a series of mammalian adaptive molecular markers, suggesting that the risk of human infection is higher. It is necessary to strengthen the surveillance by multi departments to deal with influenza pandemic.

Identification of the prion-related protein gene (PRNT) sequences in various species of the Cervidae family.

Chronic wasting disease (CWD) is caused by abnormal deleterious prion protein (PrPSc), and transmissible spongiform encephalopathy occurs in the Cervidae family. In recent studies, the susceptibility of prion disease has been affected by polymorphisms of the prion gene family. However, the study of the prion-related protein gene (PRNT) is rare, and the DNA sequence of this gene was not fully reported in all Cervidae families. In the present study, we amplified and first identified PRNT DNA sequences in the Cervidae family, including red deer, elk, sika deer and Korean water deer, using polymerase chain reaction (PCR). We aligned nucleotide sequences of the PRNT gene and the amino acid sequences of prion-related protein (Prt) protein among several species. In addition, we performed phylogenetic analysis to measure the evolutionary relationships of the PRNT gene in the Cervidae family. Furthermore, we performed homology modeling of the Prt protein using SWISS-MODEL and compared the structure of Prt protein between sheep and the Cervidae family using the Swiss-PdbViewer program. We obtained much longer PRNT sequences of red deer compared to the PRNT gene sequence registered in GenBank. Korean water deer denoted more close evolutionary distances with goats and cattle than the Cervidae family. We found 6 Cervidae family-specific amino acids by the alignment of Prt amino acid sequences. There are significantly different distributions of hydrogen bonds and the atomic distance of the N-terminal tail and C-terminal tail between sheep and the Cervidae family. We also detected the mRNA expression of PRNT gene in 3 tissues investigated. To our knowledge, this report is the first genetic study of the PRNT gene in the Cervidae family.

Identification and characterization of novel thermostable α-amylase from Geobacillus sp. GS33

In this study, the heterologous expression and biochemical characterization of a thermostable α-amylase from Geobacillus sp. GS33 was investigated. The recombinant α-amylase was overexpressed in Escherichia coli BL21 (λDE) and purified via anion exchange and size-exclusion chromatography. The purified α-amylase had a molecular weight of about 60 kDa, and was active in a broad range of pH 3–10 and temperature (40–90 °C) with maximum activity at pH 7–8 and 60 °C. The enzyme retained 50% residual activity at 65 °C, but only 20% at 85 °C after 16 h. At pH 9 and pH 7, the residual activity at 65 °C was 50% and 30%, respectively. The enzyme was remarkably activated by Co2+, Ca2+, Mg2+, PMSF, DTT, and Triton X-100, but partially inhibited by Cu2+, methanol, hexane, ethanol, acetone, SDS, and Tween 20. A molecular phylogeny analysis showed that the enzyme's amino acid sequence had the closest connection with an α-amylase from Geobacillus thermoleovorans subsp. stromboliensis nov. 3D-structure-based amino acid sequence alignments revealed that the three catalytic residues (D217, E246, D314) and the four Ca2+ ion coordination residues (N143, E177, D186, H221) were conserved in α-amylase from Geobacillus sp. GS33. The temperature stability and neutral pH optimum suggest that the enzyme may be useful for industrial applications.

Cloning and expression of a novel α-galactosidase from Lactobacillus amylolyticus L6 with hydrolytic and transgalactosyl properties.

Lactobacillus amylolyticus L6, a gram-positive amylolytic bacterium isolated from naturally fermented tofu whey (NFTW), was able to hydrolyze raffinose and stachyose for the production of α-galactosidase. The cell-free extract of L. amylolyticus L6 was found to exhibit glycosyltransferase activity to synthesize α-galacto-oligosaccharides (GOS) with melibiose as substrate. The coding genes of α-galactosidase were identified in the genome of L. amylolyticus L6. The α-galactosidase (AglB) was placed into GH36 family by amino acid sequence alignments with other α-galactosidases from lactobacilli. The optimal reaction conditions of pH and temperature for AglB were pH 6.0 and 37°C, respectively. Besides, potassium ion was found to improve the activity of AglB while divalent mercury ion, copper ion and zinc ion displayed different degrees of inhibition effect. Under the optimum reaction condition, AglB could catalyze the synthesis of GOS with degree of polymerization (DP) ≥5 by using 300 mM melibiose concentration as substrate. The maximum yield of GOS with (DP) ≥3 could reach 31.56% (w/w). Transgalactosyl properties made AglB a potential candidate for application in the production of GOS.

Epidemiology, genetic diversity and evolution of canine astrovirus.

Astroviruses (AstVs) are major causative agents of gastroenteritis in children and have been reported in many species. Canine astrovirus (CaAstV), as an enteric pathogen, has been widely detected worldwide, but little is known about their genetic diversity and evolution, partly owing to a lack of genomic data. Here, we sequenced 12 nearly full-length CaAstV genomes to address the gap in knowledge. We found 14 (13.2%) and 7 (3.35%) CaAstV-positive samples from pet dogs with and without diarrhoea, respectively. Co-infections were with co-infection with Torque teno canis virus (TTCaV) reported for the first time. Phylogenetic analysis of the ORF2 gene revealed four major lineages. In particular, lineage 4 might have evolved from a recombinant virus from lineage 2 and lineage 3. The strains sequenced here clustered with lineages 2, 3 and 4 in contrast with other Chinese strains identified previously that clustered with lineages 2 and 4. Amino acid sequence alignment within lineage revealed intralineage amino acid diversity and that the type of epidemic strains within lineages changes over time. Three amino acids substitutions located in predicted B-cell epitopes, which might be involved escape of host immunity. Moreover, frequent inter-clade ORF2 gene recombinants were identified. The identification of individual recombination events and a recombinant lineage indicated that recombination plays a crucial role in CaAstV genetic evolution and diversity by generating divergent viruses. Moreover, phylogenetic analysis of ORF1b, the most conserved gene of astrovirus, revealed a close relationship between CaAstV and California sea lion astroviruses. Overall, we report detailed information on the genetic evolution and diversity of CaAstV, which indicates that CaAstV may pose challenges for diagnostics and future control strategies.

SPAligner: alignment of long diverged molecular sequences to assembly graphs

BackgroundGraph-based representation of genome assemblies has been recently used in different contexts — from improved reconstruction of plasmid sequences and refined analysis of metagenomic data to read error correction and reference-free haplotype reconstruction. While many of these applications heavily utilize the alignment of long nucleotide sequences to assembly graphs, first general-purpose software tools for finding such alignments have been released only recently and their deficiencies and limitations are yet to be discovered. Moreover, existing tools can not perform alignment of amino acid sequences, which could prove useful in various contexts — in particular the analysis of metagenomic sequencing data.ResultsIn this work we present a novel SPAligner (Saint-Petersburg Aligner) tool for aligning long diverged nucleotide and amino acid sequences to assembly graphs. We demonstrate that SPAligner is an efficient solution for mapping third generation sequencing reads onto assembly graphs of various complexity and also show how it can facilitate the identification of known genes in complex metagenomic datasets.ConclusionsOur work will facilitate accelerating the development of graph-based approaches in solving sequence to genome assembly alignment problem. SPAligner is implemented as a part of SPAdes tools library and is available on Github.

Improving the catalytic thermostability of Bacillus altitudinis W3 ω-transaminase by proline substitutions

As a green biocatalyst, transaminase with high thermostability can be better employed to synthesize many pharmaceutical intermediates in industry. To improve the thermostability of (R)-selective amine transaminase from Bacillus altitudinis W3, related mutation sites were determined by multiple amino acid sequence alignment between wild-type ω-transaminase and four potential thermophilic ω-transaminases, followed by replacement of the related amino acid residues with proline by site-directed mutagenesis. Three stabilized mutants (D192P, T237P, and D192P/T237P) showing the highest stability were obtained and used for further analysis. Comparison with the wild-type enzyme revealed that the double mutant D192P/T237P exhibited the largest shift in thermostability, with a 2.5-fold improvement of t 1/2 at 40°C, and a 6.3°C increase in T 50 15, and a 5°C higher optimal catalytic temperature. Additionally, this mutant exhibited an increase in catalytic efficiency (k cat/K m) relative to the wild-type enzyme. Modeling analysis indicated that the improved thermostability of the mutants could be associated with newly formed hydrophobic interactions and hydrogen bonds. This study shown that proline substitutions guided by sequence alignment to improve the thermostability of (R)-selective amine transaminase was effective and this method can also be used to engineering other enzymes.

Molecular characterization, catalytic, kinetic and thermodynamic properties of protease produced by a mutant of Bacillus cereus-S6-3

The proteolytic strain Bacillus cereus-S6-3 was subjected to mutagenic treatments viz. UV irradiations and methyl methane sulfonate (MMS). The obtained mutant strain, B. cereus-S6-3/UM90 showed 1.34 fold over the parent strain. Molecular characterization of proteases from the parent (PP/S6-3) and mutant (PM/UM90) strains indicated that they were consisted of two domains and binds a zinc ion and 4 calcium ions in the active site. Amino acid sequence alignment of PM/UM90 protease showed 19 amino acid residues were substituted compared to that of the wild-type enzyme. However, both proteases contained equal number of aromatic and hydrophobic amino acids. Protease from PM/UM90 showed an effective improvement in thermal properties in terms of reaction temperature, t1/2, the values of kd, activation energy (Ea), and decimal reduction time (D) within the temperature range from 60 to 80 °C. In addition, the kinetic and thermodynamic parameters for substrate hydrolysis (i.e., Km, Vmax, ΔH*, ΔG*, ΔS*, kcat, Vmax/Km, kcat/Km, ΔG*E–T and ΔG*E–S) showed a significant improvement of the catalytic efficiency for PM/UM90 protease. Furthermore, the correlation between thermodynamic properties and the patterns of amino acid substitution of wild-type enzyme to has been discussed.