A two-phase approach for detecting recombination in nucleotide sequences

Genetic recombination can produce heterogeneous phylogenetic histories within a set of homologous genes. Delineating recombination events is important in the study of molecular evolution, as inference of such events provides a clearer picture of the phylogenetic relationships among different gene sequences or genomes. Nevertheless, detecting recombination events can be a daunting task, as the performance of different recombination-detecting approaches can vary, depending on evolutionary events that take place after recombination. We recently evaluated the effects of post-recombination events on the prediction accuracy of recombination-detecting approaches using simulated nucleotide sequence data. The main conclusion, supported by other studies, is that one should not depend on a single method when searching for recombination events. In this paper, we introduce a two-phase strategy, applying three statistical measures to detect the occurrence of recombination events, and a Bayesian phylogenetic approach in delineating breakpoints of such events in nucleotide sequences. We evaluate the performance of these approaches using simulated data, and demonstrate the applicability of this strategy to empirical data. The two-phase strategy proves to be time-efficient when applied to large datasets, and yields high-confidence results.

Molecular phylogeny of Wolbachia endosymbionts in southeast asian mosquitoes (Diptera: Culicidae) based on wsp gene sequences

Wolbachia are maternally inherited intracellular bacteria that infect a wide range of arthropods and nematodes and are associated with various reproductive abnormalities in their hosts. Insect-associated Wolbachia form a monophyletic clade in the α-Proteobacteria and recently have been separated into two supergroups (A and B) and 19 groups. Our recent polymerase chain reaction (PCR) survey using wsp specific primers indicated that various strains of Wolbachia were present in mosquitoes collected from Southeast Asia. Here, we report the phylogenetic relationship of the Wolbachia strains found in these mosquitoes using wsp gene sequences. Our phylogenetic analysis revealed eight new Wolbachia strains, five in the A supergroup and three in the B supergroup. Most of the Wolbachia strains present in Southeast Asian mosquitoes belong to the established Mors, Con, and Pip groups.

wsp Gene sequences from the Wolbachia of filarial nematodes

Wolbachia endosymbiotic bacteria are widespread in arthropods and are also present in filarial nematodes. Almost all filarial species so far examined have been found to harbor these endosymbionts. The sequences of only three genes have been published for nematode Wolbachia (i.e., the genes coding for the proteins FtsZ and catalase and for 16S rRNA). Here we present the sequences of the genes coding for the Wolbachia surface protein (WSP) from the endosymbionts of eight species of filaria. Complete gene sequences were obtained from the endosymbionts of two different species, Dirofilaria immitis and Brugia malayi. These sequences allowed us to design general primers for amplification of the wsp gene from the Wolbachia of all filarial species examined. For these species, partial WSP sequences (about 600 base pairs) were obtained with these primers. Phylogenetic analysis groups these nematode wsp sequences into a coherent cluster. Within the nematode cluster, wsp-based Wolbachia phylogeny matches a previous phylogeny obtained with ftsZ gene sequences, with a good consistency of the phylogeny of hosts (nematodes) and symbionts (Wolbachia). In addition, different individuals of the same host species (Dirofilaria immitis and Wuchereria bancrofti) show identical wsp gene sequences.

A two-phase strategy for detecting recombination in nucleotide sequences

Genetic recombination can produce heterogeneous phylogenetic histories within a set of homologous genes. Delineating recombination events is important in the study of molecular evolution, as inference of such events provides a clearer picture of the phylogenetic relationships among different gene sequences or genomes. Nevertheless, detecting recombination events can be a daunting task, as the performance of different recombination-detecting approaches can vary, depending on evolutionary events that take place after recombination. We previously evaluated the effects of post-recombination events on the prediction accuracy of recombination-detecting approaches using simulated nucleotide sequence data. The main conclusion, supported by other studies, is that one should not depend on a single method when searching for recombination events. In this paper, we introduce a two-phase strategy, applying three statistical measures to detect the occurrence of recombination events, and a Bayesian phylogenetic approach to delineate breakpoints of such events in nucleotide sequences. We evaluate the performance of these approaches using simulated data, and demonstrate the applicability of this strategy to empirical data. The two-phase strategy proves to be time-efficient when applied to large datasets, and yields high-confidence results.

Bellerophon: a program to detect chimeric sequences in multiple sequence alignments

Bellerophon is a program for detecting chimeric sequences in multiple sequence datasets by an adaption of partial treeing analysis. Bellerophon was specifically developed to detect 16S rRNA gene chimeras in PCR-clone libraries of environmental samples but can be applied to other nucleotide sequence alignments.

PCR bias toward the wild-type k-ras and p53 sequences: Implications for PCR detection of mutations and cancer diagnosis

PCR-based cancer diagnosis requires detection of rare mutations in k-ras, p53 or other genes. The assumption has been that mutant and wild-type sequences amplify with near equal efficiency, so that they are eventually present in proportions representative of the starting material. Work factor IX suggests that this assumption is invalid for one case of near-sequence identity To test the generality of this phenomenon and its relevance to cancer diagnosis, primers distant from point mutations in p53 and k-ras were used to amplify, wild-type and mutant sequences from these genes. A substantial bias against PCR amplification of mutants was observed for two regions of the p53 gene and one region of k-ras. For kras and p53, bias was observed when the wild-type and mutant sequences were amplified separately or when mixed in equal proportions before PCR. Bias was present with proofreading and non-proofreading polymerases. Mutant and wild-type segments of the factor V cystic fibrosis transmembrane conductance regulator and prothrombin genes were amplified and did not exhibit PCR bias. Therefore, the assumption of equal PCR efficiency for point mutant and wild-type sequences is invalid in several systems. Quantitative or diagnostic PCR will require validation for each locus, and enrichment strategies may be needed to optimize detection of mutants.

Structural basis of recognition of monopartite and bipartite nuclear localization sequences by mammalian importin-alpha

Importin-alpha is the nuclear import receptor that recognizes cargo proteins which contain classical monopartite and bipartite nuclear localization sequences (NLSs), and facilitates their transport into the nucleus. To determine the structural basis of the recognition of the two classes of NLSs by mammalian importin-alpha, we co-crystallized an N-terminally truncated mouse receptor protein with peptides corresponding to the monopartite NLS from the simian virus 40 (SV40) large T-antigen, and the bipartite NLS from nucleoplasmin. We show that the monopartite SV40 large T-antigen NLS binds to two binding sites on the receptor, similar to what was observed in yeast importin-alpha. The nucleoplasmin NLS-importin-alpha complex shows, for the first time, the mode of binding of bipartite NLSs to the receptor. The two basic clusters in the NLS occupy the two binding sites used by the monopartite NLS, while the sequence linking the two basic clusters is poorly ordered, consistent with its tolerance to mutations. The structures explain the structural basis for binding of diverse NLSs to the sole receptor protein. (C) 2000 Academic Press.

Analysis of Mouse Keratin 6a Regulatory Sequences in Transgenic Mice Reveals Constitutive, Tissue-Specific Expression by a Keratin 6a Minigene

The analysis of keratin 6 expression is complicated by the presence of multiple isoforms that are expressed constitutively in a number of internal stratified epithelia, in palmoplantar epidermis, and in the companion cell layer of the hair follicle. In addition, keratin 6 expression is inducible in interfollicular epidermis and the outer root sheath of the follicle, in response to wounding stimuli, phorbol esters, or retinoic acid. In order to establish the critical regions involved in the regulation of keratin 6a (the dominant isoform in mice), we generated transgenic mice with two different-sized mouse keratin 6a constructs containing either 1.3 kb or 0.12 kb of 5' flanking sequence linked to the lacZ reporter gene. Both constructs also contained the first intron and the 3' flanking sequence of mouse keratin 6a. Ectopic expression of either transgene was not observed. Double-label immunofluorescence analyses demonstrated expression of the reporter gene in keratin 6 expressing tissues, including the hair follicle, tongue, footpad, and nail bed, showing that both transgenes retained keratinocyte-specific expression. Quantitative analysis of beta -galactosidase activity verified that both the 1.3 and 0.12 kb keratin 6a promoter constructs produced similar levels of the reporter. Notably, both constructs were constitutively expressed in the outer root sheath and interfollicular epidermis in the absence of any activating stimulus, suggesting that they lack the regulatory elements that normally silence transcription in these cells. This study has revealed that a keratin 6a minigene contains critical cis elements that mediate tissue-specific expression and that the elements regulating keratin 6 induction lie distal to the 1.3 kb promoter region.

Chimeric 16S rDNA sequences of diverse origin are accumulating in the public databases

A significant number of chimeric 16S rDNA sequences of diverse origin were identified in the public databases by partial treeing analysis. This suggests that chimeric sequences, representing phylogenetically novel non-existent organisms, are routinely being overlooked in molecular phylogenetic surveys despite a general awareness of PCR-generated artefacts amongst researchers.

Intra- and interspecific divergence in the nuclear sequences of the clock gene period in species of the Drosophila buzzatii cluster

We have measured nucleotide variation in the CLOCK/CYCLE heterodimer inhibition domain (CCID) of the clock X-linked gene period in seven species belonging to the Drosophila buzzatii cluster, namely D. buzzatii, Drosophila koepferae, Drosophila antonietae, Drosophila serido, Drosophila gouveai, Drosophila seriema and Drosophila borborema. We detected that the purifying selection is the main force driving the sequence evolution in period, in agreement with the important role of CCID in clock machinery. Our survey revealed that period provides valuable phylogenetic information that allowed to resolve phylogenetic relationships among D. gouveai, D. borborema and D. seriema, which composed a polytomic clade in preliminary studies. The analysis of patterns of intraspecific variation revealed two different lineages of period in D. koepferae, probably reflecting introgressive hybridization from D. buzzatii, in concordance with previous molecular data.

Phylogenetic Analysis of Hepatitis B Virus Genotype F Complete Genome Sequences From Chilean Patients With Chronic Infection

Molecular epidemiological data concerning the hepatitis B virus (HBV) in Chile are not known completely. Since the HBV genotype F is the most prevalent in the country, the goal of this study was to obtain full HBV genome sequences from patients infected chronically in order to determine their subgenotypes and the occurrence of resistance-associated mutations. Twenty-one serum samples from antiviral drug-naive patients with chronic hepatitis B were subjected to full-length PCR amplification, and both strands of the whole genomes were fully sequenced. Phylogenetic analyses were performed along with reference sequences available from GenBank (n = 290). The sequences were aligned using Clustal X and edited in the SE-AL software. Bayesian phylogenetic analyses were conducted by Markov Chain Monte Carlo simulations (MCMC) for 10 million generations in order to obtain the substitution tree using BEAST. The sequences were also analyzed for the presence of primary drug resistance mutations using CodonCode Aligner Software. The phylogenetic analyses indicated that all sequences were found to be the HBV subgenotype F1b, clustered into four different groups, suggesting that diverse lineages of this subgenotype may be circulating within this population of Chilean patients. J. Med. Virol. 83: 1530-1536, 2011. (C) 2011 Wiley-Liss, Inc.

Prevalence and distribution of the GBV-C/HGV among HIV-1-infected patients under anti-retroviral therapy

Infection with GB virus C (GBV-C) or hepatitis G virus (HGV) is highly prevalent among HIV/AIDS patients. GBV-C/HGV viremia has not been associated with liver disease and seems to slow HIV disease progression. To study the GBV-C/HGV genotypes prevalence among HIV/AIDS patients and its association with HIV viral load (VL) and CD4+ lymphocyte counts. From February 2003 to February 2004, we analyzed 210 HIV-1-infected subjects who were on anti-retroviral therapy (ART). For 63 of them a PCR-nested to the non-coding 5` (5`NCR) region of the GBV-C/HGV was done, and for 49 a DNA direct sequencing was done. A phylogenetic analysis was performed by PHYLIP program. 63(30%) of the HIV-1-infected patients were co-infected with GBV-C/HGV. The phylogenetic analysis revealed the following genotypes (and respective relative frequencies): 1(10%), 2a (41%), 2b (43%), and 3 (6%). Co-infected patients presented lower HIV-1 VL and higher T CD4+ lymphocyte cells counts as compared with patients negative for GBV-C/HGV sequences (log = 4.52 vs. 4.71, p = 0.036), and T CD4+ lymphocyte counts (cells/mm(3) = 322.6 vs. 273.5, p = 0.081, respectively). T CD4+ cells counts equal to, or higher than, 200/mm(3) were significantly more common among co-infected patients than among HIV-infected-only patients (p = 0.042). The lowest T CD4+ cells counts were associated with genotype 1 and the highest with genotype 2b (p = 0.05). The GBV-C/HGV infection prevalence was 30% among HIV-1-infected subjects, and was associated with lower VL and higher CD4+ lymphocyte counts. GBV-C/HGV genotype 2b may be associated with better immunological response. Published by Elsevier B.V.