An algorithm to find similar internal sequence repeats

In recent years, identification of sequence patterns has been given immense importance to understand better their significance with respect to genomic organization and evolutionary processes. To this end, an algorithm has been derived to identify all similar sequence repeats present in a protein sequence. The proposed algorithm is useful to correlate the three-dimensional structure of various similar sequence repeats available in the Protein Data Bank against the same sequence repeats present in other databases like SWISS-PROT, PIR and Genome databases.

Whole-genome sequencing identifies EN1 as a determinant of bone density and fracture

The extent to which low-frequency (minor allele frequency (MAF) between 1-5%) and rare (MAF genome sequencing (n = 2,882 from UK10K (ref. 10); a population-based genome sequencing consortium), whole-exome sequencing (n = 3,549), deep imputation of genotyped samples using a combined UK10K/1000 Genomes reference panel (n = 26,534), and de novo replication genotyping (n = 20,271). We identified a low-frequency non-coding variant near a novel locus, EN1, with an effect size fourfold larger than the mean of previously reported common variants for lumbar spine BMD (rs11692564(T), MAF = 1.6%, replication effect size = +0.20 s.d., Pmeta = 2 x 10(-14)), which was also associated with a decreased risk of fracture (odds ratio = 0.85; P = 2 x 10(-11); ncases = 98,742 and ncontrols = 409,511). Using an En1(cre/flox) mouse model, we observed that conditional loss of En1 results in low bone mass, probably as a consequence of high bone turnover. We also identified a novel low-frequency non-coding variant with large effects on BMD near WNT16 (rs148771817(T), MAF = 1.2%, replication effect size = +0.41 s.d., Pmeta = 1 x 10(-11)). In general, there was an excess of association signals arising from deleterious coding and conserved non-coding variants. These findings provide evidence that low-frequency non-coding variants have large effects on BMD and fracture, thereby providing rationale for whole-genome sequencing and improved imputation reference panels to study the genetic architecture of complex traits and disease in the general population.

Simple repetitive sequences in the genome: structure and functional significance

The current explosion of DNA sequence information has generated increasing evidence for the claim that noncoding repetitive DNA sequences present within and around different genes could play an important role in genetic control processes, although the precise role and mechanism by which these sequences function are poorly understood. Several of the simple repetitive sequences which occur in a large number of loci throughout the human and other eukaryotic genomes satisfy the sequence criteria for forming non-B DNA structures in vitro. We have summarized some of the features of three different types of simple repeats that highlight the importance of repetitive DNA in the control of gene expression and chromatin organization. (i) (TG/CA)n repeats are widespread and conserved in many loci. These sequences are associated with nucleosomes of varying linker length and may play a role in chromatin organization. These Z-potential sequences can help absorb superhelical stress during transcription and aid in recombination. (ii) Human telomeric repeat (TTAGGG)n adopts a novel quadruplex structure and exhibits unusual chromatin organization. This unusual structural motif could explain chromosome pairing and stability. (iii) Intragenic amplification of (CTG)n/(CAG)n trinucleotide repeat, which is now known to be associated with several genetic disorders, could down-regulate gene expression in vivo. The overall implications of these findings vis-à-vis repetitive sequences in the genome are summarized.

Unusual RNA plant virus integration in the soybean genome leads to the production of small RNAs

Horizontal gene transfer (HGT) is known to be a major force in genome evolution. The acquisition of genes from viruses by eukaryotic genomes is a well-studied example of HGT, including rare cases of non-retroviral RNA virus integration. The present study describes the integration of cucumber mosaic virus RNA-1 into soybean genome. After an initial metatranscriptomic analysis of small RNAs derived from soybean, the de novo assembly resulted a 3029-nt contig homologous to RNA-1. The integration of this sequence in the soybean genome was confirmed by DNA deep sequencing. The locus where the integration occurred harbors the full RNA-1 sequence followed by the partial sequence of an endogenous mRNA and another sequence of RNA-1 as an inverted repeat and allowing the formation of a hairpin structure. This region recombined into a retrotransposon located inside an exon of a soybean gene. The nucleotide similarity of the integrated sequence compared to other Cucumber mosaic virus sequences indicates that the integration event occurred recently. We described a rare event of non-retroviral RNA virus integration in soybean that leads to the production of a double-stranded RNA in a similar fashion to virus resistance RNAi plants.

Molecular methods for the epidemiological analysis of methicillin-resistant Staphylococcus aureus (MRSA) and Streptococcus pneumoniae

Methicillin-resistant Staphylococcus aureus (MRSA) and Streptococcus pneumoniae are major health problems worldwide, both found in symptomless carriage but also causing even life-threatening infections. The aim of this thesis was to characterise MRSA and S. pneumoniae in detail by using several molecular typing methods for various epidemiological purposes: clonality analysis, epidemiological surveillance, outbreak investigation, and virulence factor analysis. The characteristics of MRSA isolates from the strain collection of the Finnish National Infectious Disease Register (NIDR) and pneumococcal isolates collected from military recruits and children with acute otitis media (AOM) were analysed using various typing techniques. Antimicrobial susceptibility testing, pulsed-field gel electrophoresis (PFGE), multilocus sequence typing (MLST), spa typing, staphylococcal cassette chromosome mec (SCCmec) typing, and the detection of Panton-Valentine leukocidin (PVL) genes were performed for MRSA isolates. Pneumococcal isolates were analysed using antimicrobial susceptibility testing, serotyping, MLST, and by detecting pilus islet 1 (PI-1) and 2 (PI-2) genes. Several international community- and hospital-associated MRSA clones were recognised in Finland. The genetic diversity among MRSA FIN-4 isolates and among FIN-16 isolates was low. Overall, MRSA blood isolates from 1997 to 2006 were genetically diverse. spa typing was found to be a highly discriminatory, rapid and accurate typing method and it also qualifies as the primary typing method in countries with a long history of PFGE-based MRSA strain nomenclature. However, additional typing by another method, e.g. PFGE, is needed in certain situations to be able to provide adequate discrimination for epidemiological surveillance and outbreak investigation. An outbreak of pneumonia was associated with one pneumococcal strain among military recruits, previously healthy young men living in a crowded setting. The pneumococcal carriage rate after the outbreak was found to be exceptionally high. PI-1 genes were detected at a rather low prevalence among pneumococcal isolates from children with AOM. However, the study demonstrated that PI-1 has existed among pneumococcal isolates prior to pneumococcal conjugate vaccine and the increased antimicrobial resistance era. Moreover, PI-1 was found to associate with the serotype rather than the genotype. This study adds to our understanding of the molecular epidemiology of MRSA strains in Finland and the importance of an appropriate genotyping method to be able to perform high-level laboratory-based surveillance of MRSA. Epidemiological and molecular analyses of S. pneumoniae add to our knowledge of the characteristics of pneumococcal strains in Finland.

Stochastic dynamics modeling of the protein sequence length distribution in genomes: implications for microbial evolution

In this paper, we report an analysis of the protein sequence length distribution for 13 bacteria, four archaea and one eukaryote whose genomes have been completely sequenced, The frequency distribution of protein sequence length for all the 18 organisms are remarkably similar, independent of genome size and can be described in terms of a lognormal probability distribution function. A simple stochastic model based on multiplicative processes has been proposed to explain the sequence length distribution. The stochastic model supports the random-origin hypothesis of protein sequences in genomes. Distributions of large proteins deviate from the overall lognormal behavior. Their cumulative distribution follows a power-law analogous to Pareto's law used to describe the income distribution of the wealthy. The protein sequence length distribution in genomes of organisms has important implications for microbial evolution and applications. (C) 1999 Elsevier Science B.V. All rights reserved.

Crystallization and preliminary X-ray diffraction studies of sortase A from Streptococcus pneumoniae

Sortases are cell-membrane-anchored cysteine transpeptidases that are essential for the assembly and anchoring of cell-surface adhesins in Gram-positive bacteria. Thus, they play critical roles in virulence, infection and colonization by pathogens. Sortases have been classified into four types based on their primary sequence and the target-protein motifs that they recognize. All Gram-positive bacteria express a class A housekeeping sortase (SrtA). Sortase A from Streptococcus pneumoniae (NP_358691) has been crystallized in two crystal forms. Diamond-shaped crystals of Delta N(59)SrtA diffracted to 4.0 angstrom resolution and belonged to a tetragonal system with unit-cell parameters a = b = 122.8, c = 86.5 angstrom, alpha = beta = gamma = 90 degrees, while rod-shaped crystals of Delta N(81)SrtA diffracted to 2.91 angstrom resolution and belonged to the monoclinic space group P2(1) with unit-cell parameters a = 66.8, b = 103.47, c = 74.79 angstrom, alpha = gamma = 90, beta = 115.65 degrees. The Matthews coefficient (V(M) = 2.77 angstrom(3) Da(-1)) with similar to 56% solvent content suggested the presence of four molecules in the asymmetric unit for Delta N(81)SrtA. Also, a multi-copy search using a monomer as a probe in the molecular-replacement method resulted in the successful location of four sortase molecules in the asymmetric unit, with statistics R = 41.61, R(free) = 46.44, correlation coefficient (CC) = 64.31, CC(free) = 57.67.

Highly fluorescent GFPm2+-based genome integration-proficient promoter probe vector to study Mycobacterium tuberculosis promoters in infected macrophages

Study of activity of cloned promoters in slow-growing Mycobacterium tuberculosis during long-term growth conditions in vitro or inside macrophages, requires a genome-integration proficient promoter probe vector, which can be stably maintained even without antibiotics, carrying a substrate-independent, easily scorable and highly sensitive reporter gene. In order to meet this requirement, we constructed pAKMN2, which contains mycobacterial codon-optimized gfpm2+ gene, coding for GFPm2+ of highest fluorescence reported till date, mycobacteriophage L5 attP-int sequence for genome integration, and a multiple cloning site. pAKMN2 showed stable integration and expression of GFPm2+ from M. tuberculosis and M. smegmatis genome. Expression of GFPm2+, driven by the cloned minimal promoters of M. tuberculosis cell division gene, ftsZ (MtftsZ), could be detected in the M. tuberculosis/pAKMN2-promoter integrants, growing at exponential phase in defined medium in vitro and inside macrophages. Stable expression from genome-integrated format even without antibiotic, and high sensitivity of detection by flow cytometry and fluorescence imaging, in spite of single copy integration, make pAKMN2 useful for the study of cloned promoters of any mycobacterial species under long-term in vitro growth or stress conditions, or inside macrophages.

Sequence and structural basis for chromosomal fragility during translocations in cancer

Chromosomal aberration is considered to be one of the major characteristic features in many cancers. Chromosomal translocation, one type of genomic abnormality, can lead to deregulation of critical genes involved in regulating important physiological functions such as cell proliferation and DNA repair. Although chromosomal translocations were thought to be random events, recent findings suggest that certain regions in the human genome are more susceptible to breakage than others. The possibility of deviation from the usual B-DNA conformation in such fragile regions has been an active area of investigation. This review summarizes the factors that contribute towards the fragility of these regions in the chromosomes, such as DNA sequences and the role of different forms of DNA structures. Proteins responsible for chromosomal fragility, and their mechanism of action are also discussed. The effect of positioning of chromosomes within the nucleus favoring chromosomal translocations and the role of repair mechanisms are also addressed.

Genome sequencing unveils a novel sea enterotoxin-carrying PVL phage in staphylococcus aureus ST772 from India

Staphylococcus aureus is a major human pathogen, first recognized as a leading cause of hospital-acquired infections. Community-associated S. aureus (CA-SA) pose a greater threat due to increase in severity of infection and disease among children and healthy adults. CA-SA strains in India are genetically diverse, among which is the sequence type (ST) 772, which has now spread to Australia, Europe and Japan. Towards understanding the genetic characteristics of ST772, we obtained draft genome sequences of five relevant clinical isolates and studied the properties of their PVL-carrying prophages, whose presence is a defining hallmark of CA-SA. We show that this is a novel prophage, which carries the structural genes of the hlb-carrying prophage and includes the sea enterotoxin. This architecture probably emerged early within the ST772 lineage, at least in India. The sea gene, unique to ST772 PVL, despite having promoter sequence characteristics typical of low expression, appears to be highly expressed during early phase of growth in laboratory conditions. We speculate that this might be a consequence of its novel sequence context. The crippled nature of the hlb-converting prophage in ST772. suggests that widespread mobility of the sea enterotoxin might be a selective force behind its `transfer' to the PVL prophage. Wild type ST772 strains induced strong proliferative responses as well as high cytotoxic activity against neutrophils, likely mediated by superantigen SEA and the PVL toxin respectively. Both proliferation and cytotoxicity were markedly reduced in a cured ST772 strain indicating the impact of the phage on virulence. The presence of SEA alongside he genes for the immune system-modulating PVL toxin may contribute to the success and virulence of ST772.

Common recognition principles across diverse sequence and structural families of sialic acid binding proteins

Sialic acids form a large family of 9-carbon monosaccharides and are integral components of glycoconjugates. They are known to bind to a wide range of receptors belonging to diverse sequence families and fold classes and are key mediators in a plethora of cellular processes. Thus, it is of great interest to understand the features that give rise to such a recognition capability. Structural analyses using a non-redundant data set of known sialic acid binding proteins was carried out, which included exhaustive binding site comparisons and site alignments using in-house algorithms, followed by clustering and tree computation, which has led to derivation of sialic acid recognition principles. Although the proteins in the data set belong to several sequence and structure families, their binding sites could be grouped into only six types. Structural comparison of the binding sites indicates that all sites contain one or more different combinations of key structural features over a common scaffold. The six binding site types thus serve as structural motifs for recognizing sialic acid. Scanning the motifs against a non-redundant set of binding sites from PDB indicated the motifs to be specific for sialic acid recognition. Knowledge of determinants obtained from this study will be useful for detecting function in unknown proteins. As an example analysis, a genome-wide scan for the motifs in structures of Mycobacterium tuberculosis proteome identified 17 hits that contain combinations of the features, suggesting a possible function of sialic acid binding by these proteins.

REPERTOIRE OF PROTEIN KINASES ENCODED IN THE GENOME OF ZEBRAFISH SHOWS REMARKABLY LARGE POPULATION OF PIM KINASES

In recent times, zebrafish has garnered lot of popularity as model organism to study human cancers. Despite high evolutionary divergence from humans, zebrafish develops almost all types of human tumors when induced. However, mechanistic details of tumor formation have remained largely unknown. Present study is aimed at analysis of repertoire of kinases in zebrafish proteome to provide insights into various cellular components. Annotation using highly sensitive remote homology detection methods revealed ``substantial expansion'' of Ser/Thr/Tyr kinase family in zebrafish compared to humans, constituting over 3% of proteome. Subsequent classification of kinases into subfamilies revealed presence of large number of CAMK group of kinases, with massive representation of PIM kinases, important for cell cycle regulation and growth. Extensive sequence comparison between human and zebrafish PIM kinases revealed high conservation of functionally important residues with a few organism specific variations. There are about 300 PIM kinases in zebrafish kinome, while human genome codes for only about 500 kinases altogether. PIM kinases have been implicated in various human cancers and are currently being targeted to explore their therapeutic potentials. Hence, in depth analysis of PIM kinases in zebrafish has opened up new avenues of research to verify the model organism status of zebrafish.

Evolution of genome organizations of squirrels (Sciuridae) revealed by cross-species chromosome painting

With complete sets of chromosome-specific painting probes derived from flow-sorted chromosomes of human and grey squirrel (Sciurus carolinensis), the whole genome homologies between human and representatives of tree squirrels (Sciurus carolinensis, Callosciurus erythraeus), flying squirrels (Petaurista albiventer) and chipmunks (Tamias sibiricus) have been defined by cross-species chromosome painting. The results show that, unlike the highly rearranged karyotypes of mouse and rat, the karyotypes of squirrels are highly conserved. Two methods have been used to reconstruct the genome phylogeny of squirrels with the laboratory rabbit (Oryctolagus cuniculus) as the out-group: ( 1) phylogenetic analysis by parsimony using chromosomal characters identified by comparative cytogenetic approaches; ( 2) mapping the genome rearrangements onto recently published sequence-based molecular trees. Our chromosome painting results, in combination with molecular data, show that flying squirrels are phylogenetically close to New World tree squirrels. Chromosome painting and G-banding comparisons place chipmunks ( Tamias sibiricus), with a derived karyotype, outside the clade comprising tree and flying squirrels. The superorder Glires (order Rodentia + order Lagomorpha) is firmly supported by two conserved syntenic associations between human chromosomes 1 and 10p homologues, and between 9 and 11 homologues.

The mitochondrial genome organization of the rice frog, Fejervarya limnocharis (Amphibia : Anura): a new gene order in the vertebrate mtDNA

The mitochondrial DNA of the rice frog, Fejervarya limnocharis (Amphibia, Anura), was obtained using long-and-accurate polymerase chain reaction (LA-PCR) combining with subcloning method. The complete nucleotide sequence (17,717 bp) of mitochondrial genome was determined subsequently. This mitochondrial genome is characterized by four distinctive features: the translocation of ND5 gene, a cluster of rearranged tRNA genes (tRNA(Thr), tRNA(Pro), tRNA(Leu) ((CUN))) a tandem duplication of tRNA(Mer) gene, and eight large 89-bp tandem repeats in the control region, as well as three short noncoding regions containing two repeated motifs existing in the gene cluster of ND5/tRNA(Thr)/tRNA(Pro)/tRNA(Leu)/tRNA(Phe). The tandem duplication of gene regions followed by deletions of supernumerary genes can be invoked to explain the shuffling of tRNAM(Met) and a cluster of tRNA and ND5 genes, as observed in this study. Both ND5 gene translocation and tandem duplication of tRNA(Met) were first observed in the vertebrate mitochondrial genomes. (c) 2004 Elsevier B.V. All rights reserved.

DIFFERENT RATES OF MITOCHONDRIAL-DNA SEQUENCE EVOLUTION IN KIRK DIK-DIK (MADOQUA-KIRKII) POPULATIONS

We have investigated evolutionary rates of the mitochondrial genome among individuals of Madoqua kirkii using the relative rate test. Our results demonstrate that individuals of two chromosome races, East African cytotype A and Southwest African cytotype D, evolve about 2.3 times faster than East African cytotype B. Cytogenetic changes, DNA repair efficiency, mutagens, and more likely, hitherto unrecognized factors will account for the rate difference we have observed. Our results suggest additional caution when using molecular clocks in the estimation of divergence time, even within lineages of closely related taxa. Rate heterogeneity in microevolutionary timescales represents a potentially important aspect of basic evolutionary processes and may provide additional insights into factors which affect genome evolution. (C) 1995 Academic Press, Inc.