SEARCH · Search PubMed
Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Explore the source record for details and available documents.
Members of the Emx homeobox gene class are expressed during embryogenesis in the brain and/or other head structures of phylogenetically diverse phyla. Here, we describe sequence, genomic structure, and molecular phylogenetic analysis of a cephalochordate (amphioxus) Emx class gene termed AmphiEmxA. The genomic structure of AmphiEmxA is very similar to that of vertebrate Emx genes, with two conserved intron sites. The Drosophila homolog empty spiracles (ems) has just one intron, which may be shared with chordates; the other has been secondarily lost in this Drosophila gene and in a cnidarian Emx-related gene. We identify a highly conserved peptide motif close to the amino terminus of Emx proteins, demonstrate its similarity to a sequence found in a variety of transcription factors, and argue that it arose through convergent evolution in homeobox and forkhead genes. Finally, our molecular phylogenetic analysis strongly supports the presence of a single Emx gene in the ancestor of chordates and gene duplication along the vertebrate lineage.
We have correlated myelin membrane structure with biochemical composition in the CNS and PNS of a phylogenetic series of animals, including elasmobranchs, teleosts, amphibians, and mammals. X-ray diffraction patterns were recorded from freshly dissected, unfixed tissue and used to determine the thicknesses of the liquid bilayer and the widths of the spaces between membranes at their cytoplasmic and extracellular appositions. The lipid and protein compositions of myelinated tissue from selected animals were determined by TLC and sodium dodecyl sulfate-polyacrylamide gel electrophoresis/immunoblotting, respectively. We found that (1) there were considerable differences in lipid (particularly glycolipid) composition, but no apparent phylogenetic trends; (2) the lipid composition did not seem to affect either the bilayer thickness, which was relatively constant, or the membrane separation; (3) the CNS of elasmobranch and teleost and the PNS of all four classes contained polypeptides that were recognized by antibodies against myelin P0 glycoprotein; (4) antibodies against proteolipid protein (PLP) were recognized only by amphibian and mammalian CNS; (5) wide extracellular spaces (ranging from 36 to 48 A) always correlated with the presence of P0-immunoreactive protein; (6) the narrowest extracellular spaces (approximately 31 A) were observed only in PLP-containing myelin; (7) the cytoplasmic space in PLP-containing myelin (approximately 31 A) averaged approximately 5 A less than that in P0-containing myelin; (8) even narrower cytoplasmic spaces (approximately 24 A) were measured when both P0 and 11-13-kilodalton basic protein were detected; (9) proteins immunoreactive to antibodies against myelin P2 basic protein were present in elasmobranch and teleost CNS and/or PNS, and in mammalian PNS, but not in amphibian tissues; and (10) among mammalian PNS myelins, the major difference in structure was a variation in membrane separation at the cytoplasmic apposition. These findings demonstrate which features of myelin structure have remained constant and which have become specifically altered as myelin composition changed during evolutionary development.
The biting midge Culicoides imicola Kieffer (Diptera: Ceratopogonidae) is the most important Old World vector of African horse sickness (AHS) and bluetongue (BT). Recent increases of BT incidence in the Mediterranean basin are attributed to its increased abundance and distribution. The phylogenetic status and genetic structure of C. imicola in this region are unknown, despite the importance of these aspects for BT epidemiology in the North American BT vector. In this study, analyses of partial mitochondrial cytochrome oxidase subunit I gene (COI) sequences were used to infer phylogenetic relationships among 50 C. imicola from Portugal, Rhodes, Israel, and South Africa and four other species of the Imicola Complex from southern Africa, and to estimate levels of matrilineal subdivision in C. imicola between Portugal and Israel. Eleven haplotypes were detected in C. imicola, and these formed one well-supported clade in maximum likelihood and Bayesian trees implying that the C. imicola samples comprise one phylogenetic species. Molecular variance was distributed mainly between Portugal and Israel, with no haplotypes shared between these countries, suggesting that female-mediated gene flow at this scale has been either limited or non-existent. Our results provide phylogenetic evidence that C. imicola in the study areas are potentially competent AHS and BT vectors. The geographical structure of the C. imicola COI haplotypes was concordant with that of BT virus serotypes in recent BT outbreaks in the Mediterranean basin, suggesting that population subdivision in its vector can impose spatial constraints on BT virus transmission.
A central task in the study of molecular evolution is the reconstruction of a phylogenetic tree from sequences of current-day taxa. The most established approach to tree reconstruction is maximum likelihood (ML) analysis. Unfortunately, searching for the maximum likelihood phylogenetic tree is computationally prohibitive for large data sets. In this paper, we describe a new algorithm that uses Structural Expectation Maximization (EM) for learning maximum likelihood phylogenetic trees. This algorithm is similar to the standard EM method for edge-length estimation, except that during iterations of the Structural EM algorithm the topology is improved as well as the edge length. Our algorithm performs iterations of two steps. In the E-step, we use the current tree topology and edge lengths to compute expected sufficient statistics, which summarize the data. In the M-Step, we search for a topology that maximizes the likelihood with respect to these expected sufficient statistics. We show that searching for better topologies inside the M-step can be done efficiently, as opposed to standard methods for topology search. We prove that each iteration of this procedure increases the likelihood of the topology, and thus the procedure must converge. This convergence point, however, can be a suboptimal one. To escape from such "local optima," we further enhance our basic EM procedure by incorporating moves in the flavor of simulated annealing. We evaluate these new algorithms on both synthetic and real sequence data and show that for protein sequences even our basic algorithm finds more plausible trees than existing methods for searching maximum likelihood phylogenies. Furthermore, our algorithms are dramatically faster than such methods, enabling, for the first time, phylogenetic analysis of large protein data sets in the maximum likelihood framework.
Very little is known about the structure of microbial communities, despite their abundance and importance to ecosystem processes. Recent work suggests that bacterial biodiversity might exhibit patterns similar to those of plants and animals. However, relative to our knowledge about the diversity of macro-organisms, we know little about patterns of relatedness in free-living bacterial communities, and relatively few studies have quantitatively examined community structure in a phylogenetic framework. Here we apply phylogenetic tools to bacterial diversity data to determine whether bacterial communities are phylogenetically structured. We find that bacterial communities tend to contain lower taxonomic diversity and are more likely to be phylogenetically clustered than expected by chance. Such phylogenetic clustering may indicate the importance of habitat filtering (where a group of closely related species shares a trait, or suite of traits, that allow them to persist in a given habitat) in the assembly of bacterial communities. Microbial communities are especially accessible for phylogenetic analysis and thus have the potential to figure prominently in the integration of evolutionary biology and community ecology.
The 18S-26S nuclear rDNA external transcribed spacer (ETS) has recently gained attention as a region that is valuable in phylogenetic analyses of angiosperms primarily because it can supplement nucleotide variation from the widely used and generally shorter internal transcribed spacers (ITS-1 and ITS-2) and thereby improve phylogenetic resolution and clade support in rDNA trees. Subrepeated ETS sequences (often occurring in the 5(') region) can, however, create a challenge for systematists interested in using ETS sequence data for phylogeny reconstruction. We sequenced the 5(')ETS for members of Lessingia (Compositae, Astereae) and close relatives (26 taxa total) to characterize the subrepeat variation across a group of closely related plant lineages and to gain improved understanding of the structure, molecular evolution, and phylogenetic utility of the region. The 5(')ETS region of Lessingia and relatives varied in length from approximately 245 to 1009 bp due to the presence of a variable number of subrepeats (one to eight). We assessed homology of the subrepeats using phylogenetic analysis and concluded that only two of the subrepeats and a portion of a third ( approximately 282 bp in total) were orthologous across Lessingia and could be aligned with confidence and included in further analyses. When the partial 5(')ETS data were combined with 3(')ETS and ITS data in phylogenetic analyses, no additional resolution of relationships among taxa was obtained beyond that found from analysis of 3(')ETS + ITS sequences. Inferred patterns of concerted evolution indicate that homogenization is occurring at a faster rate in the 3(')ETS and ITS regions than in the 5(')ETS region. Additionally, homogenization appears to be acting within but not among subrepeats of the same rDNA array. We conclude that challenges in assessing subrepeat orthology across taxa greatly limit the utility of the 5(')ETS region for phylogenetic analyses among species of Lessingia.
We developed a theoretical framework based on phylogenetic comparative methods to integrate phylogeny into three measures of biodiversity: species variability, richness, and evenness. These metrics can be used in conjunction with permutation procedures to test for phylogenetic community structure. As an illustration, we analyzed data on the composition of 58 lake fish communities in Wisconsin. The fish communities showed phylogenetic underdispersion, with communities more likely to contain closely related species. Using information about differences in environmental characteristics among lakes, we demonstrated that phylogenetic underdispersion in fish communities was associated with environmental factors. For example, lakes with low pH were more likely to contain species in the same clade of acid-tolerant species. Our metrics differ from existing metrics used to calculate phylogenetic community structure, such as net relatedness index and Faith's phylogenetic diversity. Our metrics have the advantage of providing an integrated and easy-to-understand package of phylogenetic measures of species variability, richness, and evenness with well-defined statistical properties. Furthermore, they allow the easy evaluation of contributions of individual species to different aspects of the phylogenetic organization of communities. Therefore, these metrics should aid with the incorporation of phylogenetic information into strategies for understanding biodiversity and its conservation.
The ITS1, ITS2, and 5.8S gene sequences of nuclear ribosomal DNA from 40 taxa of the family Heteroderidae (including the genera Afenestrata, Cactodera, Heterodera, Globodera, Punctodera, Meloidodera, Cryphodera, and Thecavermiculatus) were sequenced and analyzed. The ITS regions displayed high levels of sequence divergence within Heteroderinae and compared to outgroup taxa. Unlike recent findings in root knot nematodes, ITS sequence polymorphism does not appear to complicate phylogenetic analysis of cyst nematodes. Phylogenetic analyses with maximum-parsimony, minimum-evolution, and maximum-likelihood methods were performed with a range of computer alignments, including elision and culled alignments. All multiple alignments and phylogenetic methods yielded similar basic structure for phylogenetic relationships of Heteroderidae. The cyst-forming nematodes are represented by six main clades corresponding to morphological characters and host specialization, with certain clades assuming different positions depending on alignment procedure and/or method of phylogenetic inference. Hypotheses of monophyly of Punctoderinae and Heteroderinae are, respectively, strongly and moderately supported by the ITS data across most alignments. Close relationships were revealed between the Avenae and the Sacchari groups and between the Humuli group and the species H. salixophila within Heteroderinae. The Goettingiana group occupies a basal position within this subfamily. The validity of the genera Afenestrata and Bidera was tested and is discussed based on molecular data. We conclude that ITS sequence data are appropriate for studies of relationships within the different species groups and less so for recovery of more ancient speciations within Heteroderidae.
Circular DNA, derived from lymphocytes of juvenile channel catfish, was used to construct lambda libraries that were screened to identify the products of immunoglobulin DH-JH excision events. Clones were characterized that contained DH to JH recombination signal joints. The signal joints represented 23-bp recombination signal sequences (RSS) identical to germline JH segments that were adjacent to DH 12-bp RSS elements. DH flanking regions within the clones were used to probe a genomic library. Three germline DH gene segments containing 11-19 bp coding regions flanked by 12-bp RSS elements with conserved heptamers and nonamers were identified. The DH locus is closely linked to the JH locus, and Southern blots indicate that the DH segments represent different single member gene families. Analysis of H chain cDNA shows that each germline DH segment was expressed in functional VDJ recombination events involving different JH segments and members of different VH families. Several aspects of CDR3 junctional diversity were evident, including deletion of coding region nucleotides, N- and P-region nucleotide additions, alternate DH reading frame utilization, and point mutations. Coding region motifs of catfish DH segments are phylogenetically conserved in some DH segments of higher vertebrates. These studies indicate that the structure, genomic organization, and recombination patterns of DH segments typically associated with higher vertebrates evolved early in vertebrate phylogeny at the level of the bony fish.
A significant fraction of the bases in a folded, structured RNA molecule participate in noncanonical base pairing interactions, often in the context of internal loops or multi-helix junction loops. The appearance of each new high-resolution RNA structure provides welcome data to guide efforts to understand and predict RNA 3D structure, especially when the RNA in question is a functionally conserved molecule. The recent publication of the crystal structure of the "Loop E" region of bacterial 5S ribosomal RNA is such an event [Correll CC, Freeborn B, Moore PB, Steitz TA, 1997, Cell 91:705-712]. In addition to providing more examples of already established noncanonical base pairs, such as purine-purine sheared pairings, trans-Hoogsteen UA, and GU wobble pairs, the structure provides the first high-resolution views of two new purine-purine pairings and a new GU pairing. The goal of the present analysis is to expand the capabilities of both chemical probing and phylogenetic analysis to predict with greater accuracy the structures of RNA molecules. First, in light of existing chemical probing data, we investigate what lessons could be learned regarding the interpretation of this widely used method of RNA structure probing. Then we analyze the 3D structure with reference to molecular phylogeny data (assuming conservation of function) to discover what alternative base pairings are geometrically compatible with the structure. The comparisons between previous modeling efforts and crystal structures show that the intricate involvements of ions and water molecules in the maintenance of non-Watson-Crick pairs render the process of correctly identifying the interacting sites in such pairs treacherous, except in cases of trans-Hoogsteen A/U or sheared A/G pairs for the adenine N1 site. The phylogenetic analysis identifies A/A, A/C, A/U and C/A, C/C, and C/U pairings isosteric with sheared A/G, as well as A/A and A/C pairings isosteric with both G/U and G/G bifurcated pairings. Thus, each non-Watson-Crick pair could be characterized by a phylogenetic signature of variations between isosteric-like pairings. In addition to the conservative changes, which form a dictionary of pairings isosterically compatible with those observed in the crystal structure, concerted changes involving several base pairs also occur. The latter covariations may indicate transitions between related but distinctive motifs within the loop E of 5S ribosomal RNA.
The study among domestic fowls (Gallus domesticus L.) of the growth allometry of various brain structures in comparison with the body weight during the embryonic and post natal period leads to two different conclusions: Every structure shows a growth alteration which occurs about hatching time and is characterized by an abrupt decrease of the slope of the reducted major axis. The comparison of the allometric characteristics of every structure with those of the entire brain allows to distinguish between three kinds of growth; regressive, steady or progressive. A slow and early growth is typical of a phylogenetically ancient structure; a dynamic and late growth shows a phylogenetically recent structure. This corroborates the corticoid trend among some striated complexes (Neostriatum, ventral Hyperstriatum) that several authors had thought of.
BACKGROUND: Endogenous retroviral sequences (ERVs) are integral parts of most eukaryotic genomes and vastly outnumber exogenous retroviruses (XRVs). ERVs with a relatively complete structure were retrieved from the genetic archives of humans and chickens, diametrically opposite representatives of vertebrate retroviruses (over 3300 proviruses), and analyzed, using a bioinformatic program, RetroTector, developed by us. This rich source of proviral information, accumulated in a local database, and a collection of XRV sequences from the literature, allowed the reconstruction of a Pol based phylogenetic tree, more extensive than previously possible. The aim was to find traits useful for classification and evolutionary studies of retroviruses. Some of these traits have been used by others, but they are here tested in a wider context than before. RESULTS: In the ERV collection we found sequences similar to the XRV-based genera: alpha-, beta-, gamma-, epsilon- and spumaretroviruses. However, the occurrence of intermediates between them indicated an evolutionary continuum and suggested that taxonomic changes eventually will be necessary. No delta or lentivirus representatives were found among ERVs. Classification based on Pol similarity is congruent with a number of structural traits. Acquisition of dUTPase occurred three times in retroviral evolution. Loss of one or two NC zinc fingers appears to have occurred several times during evolution. Nucleotide biases have been described earlier for lenti-, delta- and betaretroviruses and were here confirmed in a larger context. CONCLUSION: Pol similarities and other structural traits contribute to a better understanding of retroviral phylogeny. "Global" genomic properties useful in phylogenies are i.) translational strategy, ii.) number of Gag NC zinc finger motifs, iii.) presence of Pro N-terminal dUTPase (dUTPasePro), iv.) presence of Pro C-terminal G-patch and v.) presence of a GPY/F motif in the Pol integrase (IN) C-terminal domain. "Local" retroviral genomic properties useful for delineation of lower level taxa are i.) host species range, ii.) nucleotide compositional bias and iii.) LTR lengths.
A gene for the Alternaria major allergen, Alt a 1, was amplified from 52 species of Alternaria and related genera, and sequence information was used for phylogenetic study. Alt a 1 gene sequences evolved 3.8 times faster and contained 3.5 times more parsimony-informative sites than glyceraldehyde-3-phosphate dehydrogenase (gpd) sequences. Analyses of Alt a 1 gene and gpd exon sequences strongly supported grouping of Alternaria spp. and related taxa into several species-groups described in previous studies, especially the infectoria, alternata, porri, brassicicola, and radicina species-groups and the Embellisia group. The sonchi species-group was newly suggested in this study. Monophyly of the Nimbya group was moderately supported, and monophyly of the Ulocladium group was weakly supported. Relationships among species-groups and among closely related species of the same species-group were not fully resolved. However, higher resolution could be obtained using Alt a 1 sequences or a combined dataset than using gpd sequences alone. Despite high levels of variation in amino acid sequences, results of in silico prediction of protein secondary structure for Alt a 1 demonstrated a high degree of structural similarity for most of the species suggesting a conservation of function.
Multiple Alignment is a new interface for performing and analyzing multiple protein structure alignments. It enables viewing levels of sequence and structure similarity on the aligned structures and performing a variety of evolutionary and bioinformatic tasks, including the construction of structure-based phylogenetic trees and minimal basis sets of structures that best represent the topology of the phylogenetic tree. It is implemented as a plugin for VMD (Visual Molecular Dynamics), which is distributed by the NIH Resource for Macromolecular Modeling and Bioinformatics at the University of Illinois.
Protein glutamine methylation at GGQ sites of protein chain release factors plays a pivotal role in the termination of translation. We report here the crystal structure of the Escherichia coli HemK protein (N5)-glutamine methyltransferase (MTase) in a binary complex with the methyl-donor product S-adenosyl-L-homocysteine (AdoHcy). HemK contains two domains: a putative substrate binding domain at the N terminus consisting of a five helix bundle and a seven-stranded catalytic domain at the C terminus that harbors the binding site for AdoHcy. The two domains are linked by a beta-hairpin. Structure-guided sequence analysis of the HemK family revealed 11 invariant residues functioning in methyl-donor binding and catalysis of methyl transfer. The putative substrate-binding domains of HemK from E.coli and Thermotoga maritima are structurally similar, despite the fact that they share very little sequence similarity. When the two proteins are aligned structurally, the helical N-terminal domain is subject to approximately 10 degrees of hinge movement relative to the C-terminal domain. The apparent hinge mobility of the two domains may reflect functional importance during the reaction cycle. Comparative phylogenetic analysis of the hemK gene and its frequent neighbor gene, prfA, which encodes a major substrate, provides evidence for several examples of lateral gene transfer.
Although it has been known for half a century that unique structures evolved in the cerebellum of anthropoid apes and became greatly enlarged in the human brain, the function of these structures still remains unknown. In an attempt to explain their function, a new concept of cerebellar capabilities is proposed, which is based both on neural evidence and on information-processing theory. The phylogenetically newest structures of the cerebellum may contribute to mental skills in much the same way that the phylogenetically older structures contribute to motor skills. In both cases, the cerebellum can send signals from the dentate nucleus to the cerebral frontal cortex via the thalamus. Signals from the older part of the dentate nucleus certainly help the frontal motor cortex to effect the skilled manipulation of muscles, and signals from the newest part of the dentate nucleus may help the frontal association cortex to effect the skilled manipulation of information or ideas. How such mental skills could have evolved in higher primates in the course of phylogenetic and ontogenetic development is shown. The validity of this new concept of cerebellar function can be tested on humans by means of tomographic brain scans.
While advances in molecular epidemiology and computational modeling have enhanced our capacity to track pathogen evolution, the accurate reconstruction of spatiotemporal transmission dynamics remains essential for developing epidemic preparedness frameworks and implementing outbreak response measures. Structured coalescent models offer a phylogeographic framework by restricting lineage coalescence events to geographically proximate host populations. Although the Bayesian structured coalescent approximation (BASTA) provides a tractable approach, contemporary phylogeographic analyses involving dozens of geographic localities and hundreds to thousands of viral genomes substantially exceed the computational capacity of existing implementations. The BASTA likelihood scales cubically with deme count and quadratically with sequence count due to matrix exponentiation and pairwise coalescent probability calculations. Here, we introduce a comprehensive algorithmic restructuring of the structured coalescent likelihood that eliminates redundancies, optimizes memory access, and exposes parallelization opportunities. Our approach reorganizes computations along three dimensions: (i) independent calculation of deme-transition probability matrices across time intervals; (ii) simultaneous evaluation of partial likelihood vectors within temporal slices; and (iii) concurrent aggregation of coalescent probabilities. Algorithmic restructuring cuts average coalescent likelihood computation by 7-8 fold, and parallelization further boosts performance to 10-26 fold, enabling joint phylogeographic analyses of dengue virus across 10 South American countries and H5N1 avian influenza across 20 Eurasian regions to finish in a fraction of prior time. This computational efficiency also enables comparison between backward-in-time structured coalescent approximations and forward-in-time phylogeographic methods, revealing that the former provides appropriately conservative posterior estimates, particularly at intermediate phylogenetic depths. We integrate our implementation into the popular BEAST X and BEAGLE software packages, with an accompanying interface in BEAUti X to easily set up the analyses, providing researchers with an accessible and scalable tool for real-time phylogeographic surveillance of rapidly evolving pathogens.