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.

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.

Analysis and application of an equilibrium model for in vitro bioassay systems with three components: Receptor, hormone and hormone-binding-protein

A simple theoretical framework is presented for bioassay studies using three component in vitro systems. An equilibrium model is used to derive equations useful for predicting changes in biological response after addition of hormone-binding-protein or as a consequence of increased hormone affinity. Sets of possible solutions for receptor occupancy and binding protein occupancy are found for typical values of receptor and binding protein affinity constants. Unique equilibrium solutions are dictated by the initial condition of total hormone concentration. According to the occupancy theory of drug action, increasing the affinity of a hormone for its receptor will result in a proportional increase in biological potency. However, the three component model predicts that the magnitude of increase in biological potency will be a small fraction of the proportional increase in affinity. With typical initial conditions a two-fold increase in hormone affinity for its receptor is predicted to result in only a 33% increase in biological response. Under the same conditions an Ii-fold increase in hormone affinity for receptor would be needed to produce a two-fold increase in biological potency. Some currently used bioassay systems may be unrecognized three component systems and gross errors in biopotency estimates will result if the effect of binding protein is not calculated. An algorithm derived from the three component model is used to predict changes in biological response after addition of binding protein to in vitro systems. The algorithm is tested by application to a published data set from an experimental study in an in vitro system (Lim et al., 1990, Endocrinology 127, 1287-1291). Predicted changes show good agreement (within 8%) with experimental observations. (C) 1998 Academic Press Limited.

Analysis of the substrate specificity of human sulfotransferases SULT1A1 and SULT1A3: Site-directed mutagenesis and kinetic studies

Sulfonation is an important metabolic process involved in the excretion and in some cases activation of various endogenous compounds and xenobiotics. This reaction is catalyzed by a family of enzymes named sulfotransferases. The cytosolic human sulfotransferases SULT1A1 and SULT1A3 have overlapping yet distinct substrate specificities. SULT1A1 favors simple phenolic substrates such as p-nitrophenol, whereas SULT1A3 prefers monoamine substrates such as dopamine. In this study we have used a variety of phenolic substrates to functionally characterize the role of the amino acid at position 146 in SULT1A1 and SULT1A3. First, the mutation A146E in SULT1A1 yielded a SULT1A3-like protein with respect to the Michaelis constant for simple phenols. The mutation E146A in SULT1A3 resulted in a SULT1A1-like protein with respect to the Michaelis constant for both simple phenols and monoamine compounds. When comparing the specificity of SULT1A3 toward tyramine with that for p-ethylphenol (which differs from tyramine in having no amine group on the carbon side chain), we saw a 200-fold preference for tyramine. The kinetic data obtained with the E146A mutant of SULT1A3 for these two substrates clearly showed that this protein preferred substrates without an amine group attached. Second, changing the glutamic acid at position 146 of SULT1A3 to a glutamine, thereby neutralizing the negative charge at this position, resulted in a 360-fold decrease in the specificity constant for dopamine. The results provide strong evidence that residue 146 is crucial in determining the substrate specificity of both SULT1A1 and SULT1A3 and suggest that there is a direct interaction between glutamic acid 146 in SULT1A3 and monoamine substrates.

A new chemolithoautotrophic arsenite-oxidizing bacterium isolated from a gold mine: Phylogenetic, physiological, and preliminary biochemical studies

A previously unknown chemolithoautotrophic arsenite-oxidizing bacterium has been isolated from a gold mine in the Northern Territory of Australia. The organism, designated NT-26, was found to be a gram-negative motile rod with two subterminal flagella. In a minimal medium containing only arsenite as the electron donor (5 mM), oxygen as the electron acceptor, and carbon dioxide-bicarbonate as the carbon source, the doubling time for chemolithoautotrophic growth was 7.6 h. Arsenite oxidation was found to be catalyzed by a periplasmic arsenite oxidase (optimum pH, 5.5). Based upon 16S rDNA phylogenetic sequence analysis, NT-26 belongs to the Agrobacterium/Rhizbium branch of the alpha-Proteobacteria and may represent a new species. This recently discovered organism is the most rapidly growing chemolithoautotrophic arsenite oxidizer known.

Characterization of the acyl carrier protein gene and the fab gene locus in Xanthomonas albilineans

A genomic region containing the fatty acid biosynthetic (fab) genes was isolated from the sugarcane leaf-scald pathogen Xanthomonasalbilineans. The order and predicted products of fabG (beta -ketoacyl reductase), acpP (acyl carrier protein), fabF(ketoacyl synthase II) and downstream genes in X. albilineans are very similar to those in Escherichia coli, with one exception. Sequence analysis, confirmed by insertional knockout and specific substrate feeding experiments, shows that the position occupied by pabC (encoding aminodeoxychorismate lyase) in other bacteria is occupied instead by pabB (encoding aminodeoxychorismate synthase component I) in X. albilineans. Downstream of pabB, X. albilineans resumes the arrangement common to characterized Gram-negative bacteria, with three transcriptionally coupled genes, encoding an ORF340 protein of undefined function, thymidylate kinase and delta' subunit of DNA polymerase III holoenzyme (HolB). Different species may obtain a common advantage from coordinated regulation of the same biosynthetic pathways using different genes in this region. (C) 2000 Federation of European Microbiological Societies. Published by Elsevier Science B.V. All rights reserved.

Expression analysis of putative vitellogenin and lipophorin receptors in honey bee (Apis mellifera L.) queens and workers

Two members of the low density lipoprotein receptor (LDLR) family were identified as putative orthologs for a vitellogenin receptor (Amvgr) and a lipophorin receptor (Amlpr) in the Apis mellifera genome. Both receptor sequences have the structural motifs characteristic of LDLR family members and show a high degree of similarity with sequences of other insects. RT-PCR analysis of Amvgr and Amlpr expression detected the presence of both transcripts in different tissues of adult female (ovary, fat body, midgut, head and specifically hypopharyngeal gland), as well as in embryos. In the head RNA samples we found two variant forms of AmLpR: a full length one and a shorter one lacking 29 amino acids in the O-linked sugar domain. In ovaries the expression levels of the two honey bee LDLR members showed opposing trends: whereas Amvgr expression was upregulated as the ovaries became activated, Amlpr transcript levels gradually declined. In situ hybridization analysis performed on ovaries detected Amvgr mRNA exclusively in germ line cells and corroborated the qPCR results showing an increase in Amvgr gene expression concomitant with follicle growth. (C) 2008 Elsevier Ltd. All rights reserved.

Identification and characterization of a G protein-coupled receptor homolog encoded by murine cytomegalovirus

This report describes the identification of a murine cytomegalovirus (MCMV) G protein-coupled receptor (GCR) homolog. This open reading frame (M33) is most closely related to, and collinear with, human cytomegalovirus UL33, and homologs are also present in human herpesvirus 6 and 7 (U12 for both viruses). Conserved counterparts in the sequenced alpha- or gammaherpesviruses have not been identified to date, suggesting that these genes encode proteins which are important for the biological characteristics of betaherpesviruses. We have detected transcripts for both UL33 and M33 as early as 3 or 4 h postinfection, and these reappear at late times. In addition, we have identified N-terminal splicing for both the UL33 and M33 RNA transcripts. For both open reading frames, splicing results in the introduction of amino acids which are highly conserved among known GCRs. To characterise the function of the M33 in the natural host, two independent MCMV recombinant viruses were prepared, each of which possesses an M33 open reading frame which has been disrupted with the beta-galactosidase gene. While the recombinant M33 null viruses showed no phenotypic differences in replication from wild-type MCMV in primary mouse embryo fibroblasts in vitro, they showed severely restricted growth in the salivary glands of infected mice. These data suggest that M33 plays an important role in vivo, in particular in the dissemination to or replication in the salivary gland, and provide the first evidence for the function of a viral GCR homolog in vivo.

Recessive COL6A2 C-globular Missense Mutations in Ullrich Congenital Muscular Dystrophy ROLE OF THE C2a SPLICE VARIANT

Ullrich congenital muscular dystrophy (UCMD) is a disabling and life-threatening disorder resulting from either recessive or dominant mutations in genes encoding collagen VI. Although the majority of the recessive UCMD cases have frameshift or nonsense mutations in COL6A1, COL6A2, or COL6A3, recessive structural mutations in the COL6A2 C-globular region are emerging also. However, the underlying molecular mechanisms have remained elusive. Here we identified a homozygous COL6A2 E624K mutation (C1 subdomain) and a homozygous COL6A2 R876S mutation (C2 subdomain) in two UCMD patients. The consequences of the mutations were investigated using fibroblasts from patients and cells stably transfected with the mutant constructs. In contrast to expectations based on the clinical severity of these two patients, secretion and assembly of collagen VI were moderately affected by the E624K mutation but severely impaired by the R876S substitution. The E624K substitution altered the electrostatic potential of the region surrounding the metal ion-dependent adhesion site, resulting in a collagen VI network containing thick fibrils and spots with densely packed microfibrils. The R876S mutation prevented the chain from assembling into triple-helical collagen VI molecules. The minute amount of collagen VI secreted by the R876S fibroblasts was solely composed of a faster migrating chain corresponding to the C2a splice variant with an alternative C2 subdomain. In transfected cells, the C2a splice variant was able to assemble into short microfibrils. Together, the results suggest that the C2a splice variant may functionally compensate for the loss of the normal COL6A2 chain when mutations occur in the C2 subdomain.

Conservation of hepatitis C virus nonstructural protein 3 amino acid sequence in viral isolates during liver transplantation

As a consequence of selective pressure exerted by the immune response during hepatitis C virus (HCV) infection, a high rate of nucleotide mutations in the viral genome is observed which leads to the emergence of viral escape mutants. The aim of this study was to evaluate the evolution of the amino acid (aa) sequence of the HCV nonstructural protein 3 (NS3) in viral isolates after liver transplantation. Six patients with HCV-induced liver disease undergoing liver transplantation (LT) were followed up for sequence analysis. Hepatitis C recurrence was observed in all patients after LT. The rate of synonymous (dS) nucleotide substitutions was much higher than that of nonsynonymous (dN) ones in the NS3 encoding region. The high values of the dS/dN ratios suggest no sustained adaptive evolution selection pressure and, therefore, absence of specific NS3 viral populations. Clinical genotype assignments were supported by phylogenetic analysis. Serial samples from each patient showed lower mean nucleotide genetic distance when compared with samples of the same HCV genotype and subtype. The NS3 samples studied had an N-terminal aa sequence with several differences as compared with reference ones, mainly in genotype 1b-infected patients. After LT, as compared with the sequences before, a few reverted aa substitutions and several established aa substitutions were observed at the N-terminal of NS3. Sites described to be involved in important functions of NS3, notably those of the catalytic triad and zinc binding, remained unaltered in terms of aa sequence. Rare or frequent aa substitutions occurred indiscriminately in different positions. Several cytotoxic T lymphocyte epitopes described for HCV were present in our 1b samples. Nevertheless, the deduced secondary structure of the NS3 protease showed a few alterations in samples from genotype 3a patients, but none were seen in 1b cases. Our data, obtained from patients under important selective pressure during LT, show that the NS3 protease remains well conserved, mainly in HCV 3a patients. It reinforces its potential use as an antigenic candidate for further studies aiming at the development of a protective immune response.

Analysis of thermodynamic non-ideality in terms of protein solvation

The effects of thermodynamic non-ideality on the forms of sedimentation equilibrium distributions for several isoelectric proteins have been analysed on the statistical-mechanical basis of excluded volume to obtain an estimate of the extent of protein solvation. Values of the effective solvation. parameter delta are reported for ellipsoidal as well as spherical models of the proteins, taken to be rigid, impenetrable macromolecular structures. The dependence of the effective solvated radius upon protein molecular mass exhibits reasonable agreement with the relationship calculated for a model in which the unsolvated protein molecule is surrounded by a 0.52-nm solvation shell. Although the observation that this shell thickness corresponds to a double layer of water molecules may be of questionable relevance to mechanistic interpretation of protein hydration, it augurs well for the assignment of magnitudes to the second virial coefficients of putative complexes in the quantitative characterization of protein-protein interactions under conditions where effects of thermodynamic non-ideality cannot justifiably be neglected. (C) 2001 Elsevier Science B.V. All rights reserved.

Analysis of the TGF beta functional pathway in epithelial ovarian carcinoma

Epithelial ovarian carcinoma is often diagnosed at an advanced stage of disease and is the leading cause of death from gynaecological neoplasia. The genetic changes that occur during the development of this carcinoma are poorly understood. It has been proposed that IGFIIR, TGF beta1 and TGF beta RII act as a functional unit in the TGF beta growth inhibitory pathway, and that somatic loss-of-function mutations in any one of these genes could lead to disruption of the pathway and subsequent loss of cell cycle control. We have examined these 3 genes in 25 epithelial ovarian carcinomas using single-stranded conformational polymorphism analysis and DNA sequence analysis. A total of 3 somatic missense mutations were found in the TGF beta RII gene, but none in IGFRII or TGF beta1. An association was found between TGF beta RII mutations and histology, with 2 out of 3 clear cell carcinomas having TGF beta RII mutations. This data supports other evidence from mutational analysis of the PTEN and beta -catenin genes that there are distinct developmental pathways responsible for the progression of different epithelial ovarian cancer histologic subtypes. (C) 2001 Cancer Research Campaign.

Phylogenetic revision of Agapophytinae subf.n. (Diptera : Therevidae) based on molecular and morphological evidence

Agapophytinae subf.n. is a highly diverse lineage of Australasian Therevidae, comprising eight described and two new genera: Agapophytus Guerin-Meneville, Acupalpa Krober, Acraspisa Krober, Belonalys Krober, Bonjeania Irwin & Lyneborg, Parapsilocephala Krober, Acatopygia Krober, Laxotela Winterton & Irwin, Pipinnipons gen.n. and Patanothrix gen.n. A genus-level cladistic analysis of the subfamily was undertaken using sixty-eight adult morphological characters and c. 1000 base pairs of the elongation factor-1 alpha (EF-1 alpha) protein coding gene. The morphological data partition produced three most parsimonious cladograms, whereas the molecular data partition gave a single most parsimonious cladogram, which did not match any of the cladograms found in the morphological analysis. The level of congruence between the data partitions was determined using the partition homogeneity test (HTF) and Wilcoxon signed ranks rest. Despite being significantly incongruent in at least one of the incongruence tests, the partitions were combined in a simultaneous analysis. The combined data yielded a single cladogram that was better supported than that of the individual partitions analysed separately. The relative contributions of the data partitions to support for individual nodes on the combined cladogram were investigated using Partitioned Bremer Support. The level of support for many nodes on the combined cladogram was non-additive and often greater than the sum of support for the respective nodes on individual partitions. This synergistic interaction between incongruent data partitions indicates a common phylogenetic signal in both partitions. It also suggests that criteria for partition combination based solely on incongruence may be misleading. The phylogenetic relationships of the genera are discussed using the combined data. A key to genera of Agapophytinae is presented, with genera diagnosed and figured. Two new genera are described: Patanothrix with a new species (Pat. skevingtoni) and Pat. wilsoni (Mann) transferred from Parapsilocephala, and Pipinnipons with a new species (Pip. kroeberi). Pipinnipons fascipennis (Krober) is transferred from Squamopygin Krober and Pip. imitans (Mann) is transferred from Agapophytus. Agapophytus bicolor (Krober) is transferred from Parapsilocephala. Agapophytus varipennis Mann is synonymised with Aga, queenslandi Krober and Aga. flavicornis Mann is synonymised with Aga. pallidicornis (Krober).

A dileucine motif targets E-cadherin to the basolateral cell surface in Madin-Darby canine kidney and LLC-PK1 epithelial cells

E-cadherin is a major adherens junction protein of epithelial cells, with a central role in cell-cell adhesion and cell polarity. Newly synthesized E-cadherin is targeted to the basolateral cell surface, We analyzed targeting information in the cytoplasmic tail of E-cadherin by utilizing chimeras of E-cadherin fused to the ectodo- main of the interleukin-2 alpha (IL-2 alpha) receptor expressed in Madin-Darby canine kidney and LLC-PK1 epithelial cells, Chimeras containing the full-length or membrane-proximal half of the E-cadherin cytoplasmic tail were correctly targeted to the basolateral domain. Sequence analysis of the membrane-proximal tail region revealed the presence of a highly conserved dileucine motif, which was analyzed as a putative targeting signal by mutagenesis. Elimination of this motif resulted in the loss of Tac/E-cadherin basolateral localization, pinpointing this dileucine signal as being both necessary and sufficient for basolateral targeting of E-cadherin, Truncation mutants unable to bind beta -catenin were correctly targeted, showing, contrary to current understanding, that beta -catenin is not required for basolateral trafficking. Our results also provide evidence that dileucine mediated targeting is maintained in UC-PK, cells despite the altered polarity of basolateral proteins with tyrosine-based signals in this cell line, These results provide the first direct insights into how E-cadherin is targeted to the basolateral membrane.

Isolation of Gemmata-like and Isosphaera-like planctomycete bacteria from soil and freshwater

New cultured strains of the planctomycete division (order Planctomycetales) of the domain Bacteria related to species in the genera Gemmata and Isosphaera were isolated from soil, freshwater, and a laboratory ampicillin solution. Phylogenetic analysis of the 16S rRNA gene from eight representative isolates showed that all the isolates were members of the planctomycete division. Six isolates clustered with Gemmata obscuriglobus and related strains, while two isolates clustered with Isosphaera pallida. A double-membrane-bounded nucleoid was observed in Gemmata-related isolates but not in Isosphaera-related isolates, consistent with the ultrastructures of existing species of each genus. Two isolates from this study represent the first planctomycetes successfully cultivated from soil.