DNA sequences of Alu elements indicate a recent replacement of the human autosomal genetic complement.

DNA sequences of neutral nuclear autosomal loci, compared across diverse human populations, provide a previously untapped perspective into the mode and tempo of the emergence of modern humans and a critical comparison with published clonally inherited mitochondrial DNA and Y chromosome measurements of human diversity. We obtained over 55 kilobases of sequence from three autosomal loci encompassing Alu repeats for representatives of diverse human populations as well as orthologous sequences for other hominoid species at one of these loci. Nucleotide diversity was exceedingly low. Most individuals and populations were identical. Only a single nucleotide difference distinguished presumed ancestral alleles from descendants. These results differ from those expected if alleles from divergent archaic populations were maintained through multiregional continuity. The observed virtual lack of sequence polymorphism is the signature of a recent single origin for modern humans, with general replacement of archaic populations.

Contrasting histories of avian and mammalian Mhc genes revealed by class II B sequences from songbirds.

To explore the evolutionary dynamics of genes in the major histocompatibility complex (Mhc) in nonmammalian vertebrates, we have amplified complete sequences of the polymorphic second (beta1) and third (beta2) exons of class II beta chain genes of songbirds. The pattern of nucleotide substitution in the antigen-binding site of sequences cloned from three behaviorally and phylogenetically divergent songbirds [scrub jays Aphelocoma coerulescens), red-winged blackbirds (Agelaius phoeniceus), and house finches (Carpodacus mexicanus) reveals that class II B genes of songbirds are subject to the same types of diversifying forces as those observed at mammalian class II loci. By contrast, the tree of avian class II B genes reveals that orthologous relationships have not been retained as in placental mammals and that, unlike class II genes in mammals, genes in songbirds and chickens have had very recent common ancestors within their respective groups. Thus, whereas the selective forces diversifying class II B genes of birds are likely similar to those in mammals, their long-term evolutionary dynamics appear to be characterized by much higher rates of concerted evolution.

Documentation of reticulate evolution in peonies (Paeonia) using internal transcribed spacer sequences of nuclear ribosomal DNA: implications for biogeography and concerted evolution.

The internal transcribed spacers (ITS) of nuclear ribosomal DNA of 33 species of genus Paeonia (Paeoniaceae) were sequenced. In section Paeonia, different patterns of nucleotide additivity were detected in 14 diploid and tetraploid species at sites that are variable in the other 12 species of the section, suggesting that reticulate evolution has occurred. Phylogenetic relationships of species that do not show additivity, and thus ostensibly were not derived through hybridization, were reconstructed by parsimony analysis. The taxa presumably derived through reticulate evolution were then added to the phylogenetic tree according to additivity from putative parents. The study provides an example of successfully using ITS sequences to reconstruct reticulate evolution in plants and further demonstrates that the sequence data could be highly informative and accurate for detecting hybridization. Maintenance of parental sequences in the species of hybrid origin is likely due to slowing of concerted evolution caused by the long generation time of peonies. The partial and uneven homogenization of parental sequences displayed in nine species of putative hybrid origin may have resulted from gradients of gene conversion. The documented hybridizations may have occurred since the Pleistocene glaciations. The species of hybrid origin and their putative parents are now distantly allopatric. Reconstruction of reticulate evolution with sequence data, therefore, provides gene records for distributional histories of some of the parental species.

Operational RNA code for amino acids: species-specific aminoacylation of minihelices switched by a single nucleotide.

The genetic code is based on aminoacylation reactions where specific amino acids are attached to tRNAs bearing anticodon trinucleotides. However, the anticodon-independent specific aminoacylation of RNA minihelix substrates by bacterial and yeast tRNA synthetases suggested an operational RNA code for amino acids whereby specific RNA sequences/structures in tRNA acceptor stems correspond to specific amino acids. Because of the possible significance of the operational RNA code for the development of the genetic code, we investigated aminoacylation of synthetic RNA minihelices with a human enzyme to understand the sequences needed for that aminoacylation compared with those needed for a microbial system. We show here that the species-specific aminoacylation of glycine tRNAs is recapitulated by a species-specific aminoacylation of minihelices. Although the mammalian and Escherichia coli minihelices differ at 6 of 12 base pairs, two of the three nucleotides essential for aminoacylation by the E. coli enzyme are conserved in the mammalian minihelix. The two conserved nucleotides were shown to be also important for aminoacylation of the mammalian minihelix by the human enzyme. A simple interchange of the differing nucleotide enabled the human enzyme to now charge the bacterial substrate and not the mammalian minihelix. Conversely, this interchange made the bacterial enzyme specific for the mammalian substrate. Thus, the positional locations (if not the actual nucleotides) for the operational RNA code for glycine appear conserved from bacteria to mammals.

The nucleotide sequence of chromosome I from Saccharomyces cerevisiae.

Chromosome I from the yeast Saccharomyces cerevisiae contains a DNA molecule of approximately 231 kbp and is the smallest naturally occurring functional eukaryotic nuclear chromosome so far characterized. The nucleotide sequence of this chromosome has been determined as part of an international collaboration to sequence the entire yeast genome. The chromosome contains 89 open reading frames and 4 tRNA genes. The central 165 kbp of the chromosome resembles other large sequenced regions of the yeast genome in both its high density and distribution of genes. In contrast, the remaining sequences flanking this DNA that comprise the two ends of the chromosome and make up more than 25% of the DNA molecule have a much lower gene density, are largely not transcribed, contain no genes essential for vegetative growth, and contain several apparent pseudogenes and a 15-kbp redundant sequence. These terminally repetitive regions consist of a telomeric repeat called W', flanked by DNA closely related to the yeast FLO1 gene. The low gene density, presence of pseudogenes, and lack of expression are consistent with the idea that these terminal regions represent the yeast equivalent of heterochromatin. The occurrence of such a high proportion of DNA with so little information suggests that its presence gives this chromosome the critical length required for proper function.

Use of a novel DNA melting profile assay for the identification of PCR-amplified genomic sequences encoding for variant-specific surface proteins from the clonal GS/M-83-H7 line of Giardia lamblia.

During infections, Giardia lamblia undergoes a continuous change of its major surface antigens, the variant-specific surface proteins (VSPs). Many studies on antigenic variation have been performed using G. lamblia clone GS/M-83-H7, which expresses surface antigen VSP H7. The present study was focused on the identification and characterization of vsp gene sequences within the genome of the clonal G. lamblia GS/M-83-H7 line. For this purpose, we applied a PCR which specifically amplified truncated sequences from the 3'-terminal region of the vsp genes. Upon cloning, most of the vsp gene amplification products were shown to be approximately identical in size and thus could not be distinguished from each other by conventional gel electrophoresis. In order to pre-estimate the sequence complexity within the large panel of vsp clones isolated, we elaborated a novel concept which facilitated our large-scale genetic screening approach: PCR products from cloned DNA molecules were generated and then subjected to a DNA melting profile assay based on the use of the LightCycler Instrument. This high-throughput assay system proved to be well suited to monitor sequence differences between the amplification products from closely related vsp genes and thus could be used for the primary, sequence-related discrimination of the corresponding clones. After testing 50 candidates, vsp clones could be divided into five groups, each characterized by an individual DNA melting profile of the corresponding amplification products. Sequence analysis of some of these 50 candidates confirmed data from the aforementioned assay in that clones were demonstrated to be identical within, but different between, the distinct groups. The nucleotide and deduced amino acid sequences of five representative vsp clones showed high similarities both among each other and also with the corresponding gene segment of the variant-specific surface antigen (VSP H7) expressed by the original GS/M-83-H7 variant type. Furthermore, three of the genomic vsp sequences turned out to be identical to vsp sequences that represented previously characterized transcription products from in vivo- or in vitro-switched GS/M-83-H7 trophozoites. In conclusion, the DNA melting profile assay seems to be a versatile tool for the PCR-based genotyping of moderately or highly diversified sequence orthologues.

Deduced amino acid sequences surrounding the fusion glycoprotein cleavage site and of the carboxyl-terminus of haemagglutinin-neuraminidase protein of the avirulent thermostable vaccine strain I-2 of Newcastle disease virus

A single-tube RT-PCR technique generated a 387 bp or 300 bp cDNA amplicon covering the F-0 cleavage site or the carboxyl (C)-terminus of the HN gene, respectively, of Newcastle disease virus (NDV) strain 1-2. Sequence analysis was used to deduce the amino acid sequences of the cleavage site of F protein and the C-terminus of HN protein, which were then compared with sequences for other NDV strains. The cleavage site of NDV strain 1-2 had a sequence Motif of (112)RKQGRLIG(119), consistent with an avirulent phenotype. Nucleotide sequencing and deduction of amino acids at the C-terminus of HN revealed that strain 1-2 had a 7-amino-acid extension (VEILKDGVREARSSR). This differs from the virulent viruses that caused outbreaks of Newcastle disease in Australia in the 1930s and 1990s, which have HN extensions of 0 and 9 amino acids, respectively. Amino acid sequence analyses of the F and HN genes of strain 1-2 confirmed its avirulent nature and its Australian origin.

On the structural repertoire of pools of short, random RNA sequences

A detailed knowledge of the mapping between sequence and structure spaces in populations of RNA molecules is essential to better understand their present-day functional properties, to envisage a plausible early evolution of RNA in a prebiotic chemical environment and to improve the design of in vitro evolution experiments, among others. Analysis of natural RNAs, as well as in vitro and computational studies, show that certain RNA structural motifs are much more abundant than others, pointing out a complex relation between sequence and structure. Within this framework, we have investigated computationally the structural properties of a large pool (10 molecules) of single-stranded, 35 nt-long, random RNA sequences. The secondary structures obtained are ranked and classified into structure families. The number of structures in main families is analytically calculated and compared with the numerical results. This permits a quantification of the fraction of structure space covered by a large pool of sequences. We further show that the number of structural motifs and their frequency is highly unbalanced with respect to the nucleotide composition: simple structures such as stem-loops and hairpins arise from sequences depleted in G, while more complex structures require an enrichment of G. In general, we observe a strong correlation between subfamilies-characterized by a fixed number of paired nucleotides-and nucleotide composition. Our results are compared to the structural repertoire obtained in a second pool where isolated base pairs are prohibited. © 2008 Elsevier Ltd. All rights reserved.

Towards the prediction of mutations in genomic sequences

Bio-systems are inherently complex information processing systems. Furthermore, physiological complexities of biological systems limit the formation of a hypothesis in terms of behavior and the ability to test hypothesis. More importantly the identification and classification of mutation in patients are centric topics in today's cancer research. Next generation sequencing (NGS) technologies can provide genome-wide coverage at a single nucleotide resolution and at reasonable speed and cost. The unprecedented molecular characterization provided by NGS offers the potential for an individualized approach to treatment. These advances in cancer genomics have enabled scientists to interrogate cancer-specific genomic variants and compare them with the normal variants in the same patient. Analysis of this data provides a catalog of somatic variants, present in tumor genome but not in the normal tissue DNA. In this dissertation, we present a new computational framework to the problem of predicting the number of mutations on a chromosome for a certain patient, which is a fundamental problem in clinical and research fields. We begin this dissertation with the development of a framework system that is capable of utilizing published data from a longitudinal study of patients with acute myeloid leukemia (AML), who's DNA from both normal as well as malignant tissues was subjected to NGS analysis at various points in time. By processing the sequencing data at the time of cancer diagnosis using the components of our framework, we tested it by predicting the genomic regions to be mutated at the time of relapse and, later, by comparing our results with the actual regions that showed mutations (discovered at relapse time). We demonstrate that this coupling of the algorithm pipeline can drastically improve the predictive abilities of searching a reliable molecular signature. Arguably, the most important result of our research is its superior performance to other methods like Radial Basis Function Network, Sequential Minimal Optimization, and Gaussian Process. In the final part of this dissertation, we present a detailed significance, stability and statistical analysis of our model. A performance comparison of the results are presented. This work clearly lays a good foundation for future research for other types of cancer.^

Single nucleotide polymorphisms in the IS900 sequence of Mycobacterium avium subsp. paratuberculosis are strain type specific

Insertion sequence IS900 is used as a target for the identification of Mycobacterium avium subsp. paratuberculosis. Previous reports have revealed single nucleotide polymorphisms within IS900. This study, which analyzed the IS900 sequences of a panel of isolates representing M. avium subsp. paratuberculosis strain types I, II, and III, revealed conserved type-specific polymorphisms that could be utilized as a tool for diagnostic and epidemiological purposes.