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.

Radiation of the Australian flora: What can comparisons of molecular phylogenies across multiple taxa tell us about the evolution of diversity in present-day communities?

The Australian fossil record shows that from ca. 25 Myr ago, the aseasonal-wet biome (rainforest and wet heath) gave way to the unique Australian sclerophyll biomes dominated by eucalypts, acacias and casuarinas. This transition coincided with tectonic isolation of Australia, leading to cooler, drier, more seasonal climates. From 3 Myr ago, aridification caused rapid opening of the central Australian and zone. Molecular phylogenies with dated nodes have provided new perspectives on how these events could have affected the evolution of the Australian flora. During the Mid-Cenozoic (25-10 Myr ago) period of climatic change, there were rapid radiations in sclerophyll taxa, such as Banksia, eucalypts, pea-flowered legumes and Allocasuarina. At the same time, taxa restricted to the aseasonal-wet biome (Nothofagus, Podocarpaceae and Araucariaceae) did not radiate or were depleted by extinction. During the Pliocene aridification, two Eremean biome taxa (Lepidium and Chenopodiaceae) radiated rapidly after dispersing into Australia from overseas. It is clear that the biomes have different histories. Lineages in the aseasonal-wet biome are species poor, with sister taxa that are species rich, either outside Australia or in the sclerophyll biomes. In conjunction with the fossil record, this indicates depletion of the Australian aseasonal-wet biome from the Mid-Cenozoic. In the sclerophyll biomes, there have been multiple exchanges between the southwest and southeast, rather than single large endemic radiations after a vicariance event. There is need for rigorous molecular phylogenetic studies so that additional questions can be addressed, such as how interactions between biomes may have driven the speciation process during radiations. New studies should include the hither-to neglected monsoonal tropics.

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.

Essential role of TM V and VI for binding the C-terminal sequences of Des-Arg-kinins

In the kallikrein-kinin and renin-angiotensin systems the main receptors, B-1 and B-2 (kinin receptors) and AT(1) and AT(2) (angiotensin receptors) respectively, are seven-transmembrane domain G-protein-coupled receptors. Considering that the B, agonists Des-Arg(9)-BK (Arg-Pro-Pro-Gly-Phe-Ser-Pro-Phe), Lys-desArg(9)-BK or Des-Arg(10)-KD (Lys-Arg-Pro-Pro-Gly-Phe-Ser-Pro-Phe) and the AT, agonist (Asp-Arg-Val-Tyr-lle-His-Pro-Phe) have the same two residues at the C-terminal region (i.e. Pro-Phe), we hypothesized that TM V and TM VI of the B-1 receptor could play an essential role in agonist binding and activity, being these regions receptor sites for binding the C-terminal sequences of Des-Arg-kinins similarly to that observed to AT, receptor. To investigate this hypothesis, we replaced Arg(212) for Ala at the top of the TM V and the sequence 274-282 (CPYHFFAFL) in TM VI of the rat kinin B, receptor by the 32 receptor homologous sequence, 289-297 (FPFQISTFL) and subsequently analyzed the consequences of these mutations by competition binding and functional assays. Despite correct expression, observed at the mRNA and protein level by RT-PCR and confocal microscopy, respectively, no agonist binding and function was verified for the mutated receptors. Therefore, our results suggest an important role for Arg(212) in the TM V and a region of TM VI of rat B, receptor in the interaction with the C-terminal residues of Des-Arg-kinins, similar to that observed with AngII. (c) 2007 Elsevier B.V. All rights reserved.

Phylogenetic analysis of the Monocotylidae (Monogenea) inferred from 28S rDNA sequences

The current classification of the Monocotylidae (Monogenea) is based on a phylogeny generated from morphological characters. The present study tests the morphological phylogenetic hypothesis using molecular methods. Sequences from domains C2 and D1 and the partial domains C1 and D2 from the 28S rDNA gene for 26 species of monocotylids from six of the seven subfamilies were used. Trees were generated using maximum parsimony, neighbour joining and maximum likelihood algorithms. The maximum parsimony tree, with branches showing less than 70% bootstrap support collapsed, had a topology identical to that obtained using the maximum likelihood analysis. The neighbour joining tree, with branches showing less than 70% support collapsed. differed only in its placement of Heterocotyle capricornensis as the sister group to the Decacotylinae clade. The molecular tree largely supports the subfamilies established using morphological characters. Differences are primarily how the subfamilies are related to each other. The monophyly of the Calicotylinae and Merizocotylinae and their sister group relationship is supported by high bootstrap values in all three methods, but relationships within the Merizocotylinae are unclear. Merizocotyle is paraphyletic and our data suggest that Mycteronastes and Thaumatocotyle, which were synonymized with Merizocotyle after the morphological cladistic analysis, should perhaps be resurrected as valid genera. The monophyly of the Monocotylinae and Decacotylinae is also supported by high bootstrap values. The Decacotylinae, which was considered previously to be the sister group to the Calicotylinae plus Merizocotylinae, is grouped in an unresolved polychotomy with the Monocotylinae and members of the Heterocotylinae. According to our molecular data, the Heterocotylinae is paraphyletic. Molecular data support a sister group relationship between Troglocephalus rhinobatidis and Neoheterocotyle rhinobatidis to the exclusion of the other species of Neoheterocotyle and recognition of Troglocephalus renders Neoheterocotyle,le paraphyletic. We propose Troglocephalus incertae sedis. An updated classification and full species list of the Monocotylidae is provided. (C) 2001 Australian Society for Parasitology Inc. Published by Elsevier Science Ltd. All rights reserved.

Sex-biased hatching sequences in the cooperatively breeding Noisy Miner

The Noisy Miner Manorina melanocephala (Meliphagidae) is a cooperatively breeding bird species in which sons often remain on their natal home ranges and help one or both of their parents. In a population of Noisy Miners in SE Queensland, Australia, a molecular technique was used to explore adult and offspring sex ratios. and also hatching sequences. Among the adult population, there were 2.31 males for every female, and roughly 99% of helping was performed by males. At hatching and fledging, the population sex ratio was even, with exactly 57 males and 57 females. However, in 17 out of 18 broods the first egg to hatch was male, First-hatched males were significantly larger and heavier than their sisters just prior to fledging. Through their helping behaviour, large healthy sons could clearly enhance the future reproductive success of parents. and benefit the entire group. Sex-biased hatching sequences could potentially provide cooperatively breeding birds with a subtle and precise way of varying investment in the helping sex.

Phylogenetic analysis of evolutionary relationships of the planctomycete division of the domain bacteria based on amino acid sequences of elongation factor Tu

Sequences from the tuf gene coding for the elongation factor EF-Tu were amplified and sequenced from the genomic DNA of Pirellula marina and Isosphaera pallida, two species of bacteria within the order Planctomycetales. A near-complete (1140-bp) sequence was obtained from Pi. marina and a partial (759-bp) sequence was obtained for I. pallida. Alignment of the deduced Pi. marina EF-Tu amino acid sequence against reference sequences demonstrated the presence of a unique Il-amino acid sequence motif not present in any other division of the domain Bacteria. Pi. marina shared the highest percentage amino acid sequence identity with I. pallida but showed only a low percentage identity with other members of the domain Bacteria. This is consistent with the concept of the planctomycetes as a unique division of the Bacteria. Neither primary sequence comparison of EF-Tu nor phylogenetic analysis supports any close relationship between planctomycetes and the chlamydiae, which has previously been postulated on the basis of 16S rRNA. Phylogenetic analysis of aligned EF-Tu amino acid sequences performed using distance, maximum-parsimony, and maximum likelihood approaches yielded contradictory results with respect to the position of planctomycetes relative to other bacteria, It is hypothesized that long-branch attraction effects due to unequal evolutionary rates and mutational saturation effects may account for some of the contradictions.

Essential role of cyclization sequences in Flavivirus RNA replication

A possible role in RNA replication for interactions between conserved complementary (cyclization) sequences in the 5'- and 3'-terminal regions of Flavivirus RNA was previously suggested but never tested in vivo. Using the M-fold program for RNA secondary-structure predictions, we examined for the first time the base-pairing interactions between the covalently linked 5' genomic region (first similar to 160 nucleotides) and the 3' untranslated region (last similar to 115 nucleotides) for a range of mosquito borne Flavivirus species. Base-pairing occurred as predicted for the previously proposed conserved cyclization sequences. In order to obtain experimental evidence of the predicted interactions, the putative cyclization sequences (5' or 3') in the replicon RNA of the mosquito-borne Kunjin virus,were mutated either separately, to destroy base-pairing, or simultaneously, to restore the complementarity. None of the RNAs with separate mutations in only the 5' or only the 3' cyclization sequences was able to replicate after transfection into BHK cells, while replicon RNA with simultaneous compensatory mutations in both cyclization sequences was replication competent. This was detected by immunofluorescence for expression of the major nonstructural protein NS3 and by Northern blot analysis for amplification and accumulation of replicon RNA. We then used the M-fold program to analyze RNA secondary structure of the covalently linked 5'- and 3'-terminal regions of three tick borne virus species and identified a previously undescribed additional pair of conserved complementary sequences in locations similar to those of the mosquito borne species. They base-paired with DeltaG values of approximately -20 kcal, equivalent or greater in stability than those calculated for the originally proposed cyclization sequences. The results show that the base-pairing between 5' and 3' complementary sequences, rather than the nucleotide sequence per se, is essential for the replication of mosquito-borne Kunjin virus RNA and that more than one pair of cyclization sequences might be involved in the replication of the tick-borne Flavivirus species.