Ultra-fast sequence clustering from similarity networks with <Emphasis FontCategory="NonProportional">SiLiX</Emphasis>

Background  

The number of gene sequences that are available for comparative genomics approaches is increasing extremely quickly. A current challenge is to be able to handle this huge amount of sequences in order to build families of homologous sequences in a reasonable time.     

A fast algorithm for determining the best combination of local alignments to a query sequence

Background  

Existing sequence alignment algorithms assume that similarities between DNA or amino acid sequences are linearly ordered. That is, stretches of similar nucleotides or amino acids are in the same order in both sequences. Recombination perturbs this order. An algorithm that can reconstruct sequence similarity despite rearrangement would be helpful for reconstructing the evolutionary history of recombined sequences.     

Hepatitis C virus genotype 3a with phylogenetically distinct origin is circulating in Pakistan

Background

Hepatitis C virus (HCV) is one of the leading causes of viral hepatitis worldwide and its genotype 3a is predominant in vast areas of Pakistan.

Findings

The present study reports the first full sequence of HCV 3a isolate PK-1 from Pakistan. This nucleotide sequence was compared with six other HCV genotype 3a full length sequences from different regions of the world by using statistical methods of phylogenetic analysis.

Conclusion

The nucleotide difference of these seven sequences shows that HCV genotype 3a of phylogenetically distinct origin is circulating in Pakistan.     

Phenotypic effect of mutations in evolving populations of RNA molecules

Background  

The secondary structure of folded RNA sequences is a good model to map phenotype onto genotype, as represented by the RNA sequence. Computational studies of the evolution of ensembles of RNA molecules towards target secondary structures yield valuable clues to the mechanisms behind adaptation of complex populations. The relationship between the space of sequences and structures, the organization of RNA ensembles at mutation-selection equilibrium, the time of adaptation as a function of the population parameters, the presence of collective effects in quasispecies, or the optimal mutation rates to promote adaptation all are issues that can be explored within this framework.     

NBLAST: a cluster variant of BLAST for NxN comparisons

Background  

The BLAST algorithm compares biological sequences to one another in order to determine shared motifs and common ancestry. However, the comparison of all non-redundant (NR) sequences against all other NR sequences is a computationally intensive task. We developed NBLAST as a cluster computer implementation of the BLAST family of sequence comparison programs for the purpose of generating pre-computed BLAST alignments and neighbour lists of NR sequences.     

Grammar-based distance in progressive multiple sequence alignment

Background  

We propose a multiple sequence alignment (MSA) algorithm and compare the alignment-quality and execution-time of the proposed algorithm with that of existing algorithms. The proposed progressive alignment algorithm uses a grammar-based distance metric to determine the order in which biological sequences are to be pairwise aligned. The progressive alignment occurs via pairwise aligning new sequences with an ensemble of the sequences previously aligned.     

Contribution of Genome-Wide HCV Genetic Differences to Outcome of Interferon-Based Therapy in Caucasian American and African American Patients

Background

Hepatitis C virus (HCV) has six major genotypes, and patients infected with genotype 1 respond less well to interferon-based therapy than other genotypes. African American patients respond to interferon α-based therapy at about half the rate of Caucasian Americans. The effect of HCV''s genetic variation on treatment outcome in both racial groups is poorly understood.

Methodology

We determined the near full-length pre-therapy consensus sequences from 94 patients infected with HCV genotype 1a or 1b undergoing treatment with peginterferon α-2a and ribavirin through the Virahep-C study. The sequences were stratified by genotype, race and treatment outcome to identify HCV genetic differences associated with treatment efficacy.

Principal Findings

HCV sequences from patients who achieved sustained viral response were more diverse than sequences from non-responders. These inter-patient diversity differences were found primarily in the NS5A gene in genotype 1a and in core and NS2 in genotype 1b. These differences could not be explained by host selection pressures. Genotype 1b but not 1a African American patients had viral genetic differences that correlated with treatment outcome.

Conclusions & Significance

Higher inter-patient viral genetic diversity correlated with successful treatment, implying that there are HCV genotype 1 strains with intrinsic differences in sensitivity to therapy. Core, NS3 and NS5A have interferon-suppressive activities detectable through in vitro assays, and hence these activities also appear to function in human patients. Both preferential infection with relatively resistant HCV variants and host-specific factors appear to contribute to the unusually poor response to therapy in African American patients.     

A multistep bioinformatic approach detects putative regulatory elements in gene promoters

Background  

Searching for approximate patterns in large promoter sequences frequently produces an exceedingly high numbers of results. Our aim was to exploit biological knowledge for definition of a sheltered search space and of appropriate search parameters, in order to develop a method for identification of a tractable number of sequence motifs.     

SeqTrim: a high-throughput pipeline for pre-processing any type of sequence read

Background  

High-throughput automated sequencing has enabled an exponential growth rate of sequencing data. This requires increasing sequence quality and reliability in order to avoid database contamination with artefactual sequences. The arrival of pyrosequencing enhances this problem and necessitates customisable pre-processing algorithms.     

Occult Hepatitis B Virus Infection in Anti-HBs-Positive Infants Born to HBsAg-Positive Mothers in China

Objective

To investigate the prevalence of occult HBV infection (OBI) among children and to characterize virology of occult HBV, we conducted an epidemiological survey.

Methods

186 HB-vaccinated infants born to HBsAg-positive mothers were included in the study. Serological tests for HBV markers were performed using commercial ELISA kits. Real-time quantitative PCR and nested PCR were used to detect HBV DNA. PCR products of the C and pre-S/S regions were sequenced and analyzed.

Results

1.61% (3/186) infants were HBsAg positive, and 4.92% (9/183) infants were considered as occult infection. The viral load of mothers was associated with occult infection (P = 0.020). Incomplete three-dose injections of HB vaccine was associated with HBV infection (P = 0.022). Six OBI infants were positive for anti-HBs, but their titers were not greater than 100 mIU/mL. Seven isolated HBV pre-S/S sequences were obtained from nine OBI infants. Three of the sequences were genotype C, and four of the sequences were genotype C/D. Escape mutation S143L was found in the four sequences of genotype C/D. All seven sequences lacked G145R and other escape mutation in S region.

Conclusions

Occult HBV infection was detected in anti-HBs positive infants born to HBsAg-positive mothers in China. Occult infection was associated with absent anti-HBs or with low anti-HBs level, high maternal viral loads and escape mutations in the S gene.     

New insight into HCV E1/E2 region of genotype 4a

Introduction

Hepatitis C virus (HCV) genome contains two envelope proteins (E1 and E2) responsible for the virus entry into the cell. There is a substantial lack of sequences covering the full length of E1/E2 region for genotype 4. Our study aims at providing new sequences as well as characterizing the genetic divergence of the E1/E2 region of HCV 4a using our new sequences along with all publicly available datasets.

Methods

The genomic segments covering the whole E1/E2 region were isolated from Egyptian HCV patients and sequenced. The resulting 36 sequences 36 were analyzed using sequence analysis techniques to study variability within and among hosts in the same time point. Furthermore, previously published HCV E1/E2 sequence datasets for genotype 4a were retrieved and categorized according to the geographical location and date of isolation and were used for further analysis of variability among Egyptian over a period of 15 years, also compared with non-Egyptian sequences to figure out region-specific variability.

Results

Phylogenetic analysis of the new sequences has shown variability within the host and among different individuals in the same time point. Analysis of the 36 sequences along with the Egyptian sequences (254 sequences in E1 in the period from 1997 to 2010 and 8 E2 sequences in the period from 2006 to 2010) has shown temporal change over time. Analysis of the new HCV sequences with the non-Egyptian sequences (182 sequences in E1 and 155 sequences in the E2) has shown region specific variability. The molecular clock rate of E1 was estimated to be 5E-3 per site per year for Egyptian and 5.38E-3 for non-Egyptian. The clock rate of E2 was estimated to be 8.48E per site per year for Egyptian and 6.3E-3 for non-Egyptian.

Conclusion

The results of this study support the high rate of evolution of the Egyptian HCV genotype 4a. It has also revealed significant level of genetic variability among sequences from different regions in the world.

     

Genotyping and inflated type I error rate in genome-wide association case/control studies

Background  

One common goal of a case/control genome wide association study (GWAS) is to find SNPs associated with a disease. Traditionally, the first step in such studies is to assign a genotype to each SNP in each subject, based on a statistic summarizing fluorescence measurements. When the distributions of the summary statistics are not well separated by genotype, the act of genotype assignment can lead to more potential problems than acknowledged by the literature.     

Distinct Distribution Pattern of Hepatitis B Virus Genotype C and D in Liver Tissue and Serum of Dual Genotype Infected Liver Cirrhosis and Hepatocellular Carcinoma Patients

Aims

The impact of co-infection of several hepatitis B virus (HBV) genotypes on the clinical outcome remains controversial. This study has for the first time investigated the distribution of HBV genotypes in the serum and in the intrahepatic tissue of liver cirrhotic (LC) and hepatocellular carcinoma (HCC) patients from India. In addition, the genotype-genotype interplay and plausible mechanism of development of HCC has also been explored.

Methods

The assessment of HBV genotypes was performed by nested PCR using either surface or HBx specific primers from both the circulating virus in the serum and replicative virus that includes covalently closed circular DNA (cccDNA) and relaxed circular DNA (rcDNA) of HBV from the intrahepatic tissue. The integrated virus within the host chromosome was genotyped by Alu-PCR method. Each PCR products were cloned and sequences of five randomly selected clones were subsequently analysed.

Results

HBV/genotype D was detected in the serum of all LC and HCC patients whereas the sequences of the replicative HBV DNA (cccDNA and rcDNA) from the intrahepatic tissue of the same patients revealed the presence of both HBV/genotype C and D. The sequences of the integrated viruses exhibited the solo presence of HBV/genotype C in the majority of LC and HCC tissues while both HBV/genotype C and D clones were found in few patients in which HBV/genotype C was predominated. Moreover, compared to HBV/genotype D, genotype C had higher propensity to generate double strand breaks, ER stress and reactive oxygen species and it had also showed higher cellular homologous-recombination efficiency that engendered more chromosomal rearrangements, which ultimately led to development of HCC.

Conclusions

Our study highlights the necessity of routine analysis of HBV genotype from the liver tissue of each chronic HBV infected patient in clinical practice to understand the disease prognosis and also to select therapeutic strategy.     

The role of genome and gene regulatory network canalization in the evolution of multi-trait polymorphisms and sympatric speciation

Background  

Sexual reproduction has classically been considered as a barrier to the buildup of discrete phenotypic differentiation. This notion has been confirmed by models of sympatric speciation in which a fixed genetic architecture and a linear genotype phenotype mapping were assumed. In this paper we study the influence of a flexible genetic architecture and non-linear genotype phenotype map on differentiation under sexual reproduction.     

Biological sequences as pictures – a generic two dimensional solution for iterated maps

Background  

Representing symbolic sequences graphically using iterated maps has enjoyed an enduring popularity since it was first proposed in Jeffrey 1990 as chaos game representation (CGR). The usefulness of this representation goes beyond the convenience of a scale independent representation, it provides a variable memory length representation of transition. This includes the representation of succession with non-integer order, which comes with the promise of generalizing Markovian formalisms. The original proposal targeted genomic sequences only but since then several generalizations have been proposed, many specifically designed to handle protein data.     

Angiotensin-converting enzyme gene insertion/deletion polymorphism in migraine patients

Background  

The main objective of this study was to investigate the angiotensin converting enzyme (ACE) genotype as a possible risk factor for migraine (both with and without aura) compared to controls. We also wanted to examine whether a clinical response to an ACE inhibitor, lisinopril, or an angiotensin II receptor blocker, candesartan, in migraine prophylaxis was related to ACE genotype.     

Comparing transformation methods for DNA microarray data

Background  

When DNA microarray data are used for gene clustering, genotype/phenotype correlation studies, or tissue classification the signal intensities are usually transformed and normalized in several steps in order to improve comparability and signal/noise ratio. These steps may include subtraction of an estimated background signal, subtracting the reference signal, smoothing (to account for nonlinear measurement effects), and more. Different authors use different approaches, and it is generally not clear to users which method they should prefer.     

Pattern statistics on Markov chains and sensitivity to parameter estimation

Background:  

In order to compute pattern statistics in computational biology a Markov model is commonly used to take into account the sequence composition. Usually its parameter must be estimated. The aim of this paper is to determine how sensitive these statistics are to parameter estimation, and what are the consequences of this variability on pattern studies (finding the most over-represented words in a genome, the most significant common words to a set of sequences,...).     

Genotype networks in metabolic reaction spaces

Background  

A metabolic genotype comprises all chemical reactions an organism can catalyze via enzymes encoded in its genome. A genotype is viable in a given environment if it is capable of producing all biomass components the organism needs to survive and reproduce. Previous work has focused on the properties of individual genotypes while little is known about how genome-scale metabolic networks with a given function can vary in their reaction content.