Genome sequencing of cells that live inside glass cages reveals their past history

One-quarter of the total primary production on earth is contributed by diatoms1. These are photosynthetic, unicellular algae with ornamented silica shells found in all aquatic and moist environments. They form the base of energy-efficient food webs that support all aquatic life forms. More than 250 genera of living diatoms, with as many as 100,000 species are known2. Fossil diatoms are known as early as the Cretaceous, 144–65 m.y. ago3. In India, deposits of diatoms occur in Rajasthan and are known as ‘multani mitti’. Multani mitti or Indian Fuller’s earth or diatomaceous earth as it is called in the West, is applied as a paste on the surface of the skin for 15–20 min and then washed-off. This leaves the skin feeling smooth, soft, moist and rejuvenated. Diatomaceous earth is now being used in the formulation of soaps, cleansing products, face powders and skincare preparations. Diatomaceous earth is a mineral material consisting mainly of siliceous fragments of various species of fossilized remains of diatoms.

Relative stability of DNA as a generic criterion for promoter prediction: whole genome annotation of microbial genomes with varying nucleotide base composition

The rapid increase in genome sequence information has necessitated the annotation of their functional elements, particularly those occurring in the non-coding regions, in the genomic context. Promoter region is the key regulatory region, which enables the gene to be transcribed or repressed, but it is difficult to determine experimentally. Hence an in silico identification of promoters is crucial in order to guide experimental work and to pin point the key region that controls the transcription initiation of a gene. In this analysis, we demonstrate that while the promoter regions are in general less stable than the flanking regions, their average free energy varies depending on the GC composition of the flanking genomic sequence. We have therefore obtained a set of free energy threshold values, for genomic DNA with varying GC content and used them as generic criteria for predicting promoter regions in several microbial genomes, using an in-house developed tool `PromPredict'. On applying it to predict promoter regions corresponding to the 1144 and 612 experimentally validated TSSs in E. coli (50.8% GC) and B. subtilis (43.5% GC) sensitivity of 99% and 95% and precision values of 58% and 60%, respectively, were achieved. For the limited data set of 81 TSSs available for M. tuberculosis (65.6% GC) a sensitivity of 100% and precision of 49% was obtained.

Suite of tools for statistical N-gram language modeling for pattern mining in whole genome sequences

Genome sequences contain a number of patterns that have biomedical significance. Repetitive sequences of various kinds are a primary component of most of the genomic sequence patterns. We extended the suffix-array based Biological Language Modeling Toolkit to compute n-gram frequencies as well as n-gram language-model based perplexity in windows over the whole genome sequence to find biologically relevant patterns. We present the suite of tools and their application for analysis on whole human genome sequence.

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.

Whole exome sequencing identifies a novel splice-site mutation in ADAMTS17 in an Indian family with Weill-Marchesani syndrome

Purpose: Weill-Marchesani syndrome (WMS) is a rare connective tissue disorder, characterized by short stature, micro-spherophakic lens, and stubby hands and feet (brachydactyly). WMS is caused by mutations in the FBN1, ADAMTS10, and LTBP2 genes. Mutations in the LTBP2 and ADAMTS17 genes cause a WMS-like syndrome, in which the affected individuals show major features of WMS but do not display brachydactyly and joint stiffness. The main purpose of our study was to determine the genetic cause of WMS in an Indian family. Methods: Whole exome sequencing (WES) was used to identify the genetic cause of WMS in the family. The cosegregation of the mutation was determined with Sanger sequencing. Reverse transcription (RT)-PCR analysis was used to assess the effect of a splice-site mutation on splicing of the ADAMTS17 transcript. Results: The WES analysis identified a homozygous novel splice-site mutation c.873+1G>T in a known WMS-like syndrome gene, ADAMTS17, in the family. RT-PCR analysis in the patient showed that exon 5 was skipped, which resulted in the deletion of 28 amino acids in the ADAMTS17 protein. Conclusions: The mutation in the WMS-like syndrome gene ADAMTS17 also causes WMS in an Indian family. The present study will be helpful in genetic diagnosis of this family and increases the number of mutations of this gene to six.

Streptococcus pneumoniae Genome Database (SPGDB): A database for strain specific comparative analysis of Streptococcus pneumoniae genes and proteins

Streptococcus pneumoniae causes pneumonia, septicemia and meningitis. S. pneumoniae is responsible for significant mortality both in children and in the elderly. In recent years, the whole genome sequencing of various S. pneumoniae strains have increased manifold and there is an urgent need to provide organism specific annotations to the scientific community. This prompted us to develop the Streptococcus pneumoniae Genome Database (SPGDB) to integrate and analyze the completely sequenced and available S. pneumoniae genome sequences. Further, links to several tools are provided to compare the pool of gene and protein sequences, and proteins structure across different strains of S. pneumoniae. SPGDB aids in the analysis of phenotypic variations as well as to perform extensive genomics and evolutionary studies with reference to S. pneumoniae. (C) 2014 Elsevier Inc. All rights reserved.

Whole exome sequencing in an Indian family links Coats plus syndrome and dextrocardia with a homozygous novel CTC1 and a rare HES7 variation

Background: Coats plus syndrome is an autosomal recessive, pleiotropic, multisystem disorder characterized by retinal telangiectasia and exudates, intracranial calcification with leukoencephalopathy and brain cysts, osteopenia with predisposition to fractures, bone marrow suppression, gastrointestinal bleeding and portal hypertension. It is caused by compound heterozygous mutations in the CTC1 gene. Case presentation: We encountered a case of an eight-year old boy from an Indian family with manifestations of Coats plus syndrome along with an unusual occurrence of dextrocardia and situs inversus. Targeted resequencing of the CTC1 gene as well as whole exome sequencing (WES) were conducted in this family to identify the causal variations. The identified candidate variations were screened in ethnicity matched healthy controls. The effect of CTC1 variation on telomere length was assessed using Southern blot. A novel homozygous missense mutation c.1451A > C (p.H484P) in exon 9 of the CTC1 gene and a rare 3'UTR known dbSNP variation (c.*556 T > C) in HES7 were identified as the plausible candidates associated with this complex phenotype of Coats plus and dextrocardia. This CTC1 variation was absent in the controls and we also observed a reduced telomere length in the affected individual's DNA, suggesting its likely pathogenic nature. The reported p.H484P mutation is located in the N-terminal 700 amino acid regionthat is important for the binding of CTC1 to ssDNA through its two OB domains. WES data also showed a rare homozygous missense variation in the TEK gene in the affected individual. Both HES7 and TEK are targets of the Notch signaling pathway. Conclusions: This is the first report of a genetically confirmed case of Coats plus syndrome from India. By means of WES, the genetic variations in this family with unique and rare complex phenotype could be traced effectively. We speculate the important role of Notch signaling in this complex phenotypic presentation of Coats plus syndrome and dextrocardia. The present finding will be useful for genetic diagnosis and carrier detection in the family and for other patients with similar disease manifestations.

Novel Rearrangements in the Staphylococcal Cassette Chromosome Mec Type V Elements of Indian ST772 and ST672 Methicillin Resistant Staphylococcus aureus Strains

Staphylococcus aureus is a commensal gram positive bacteria which causes severe and non severe infections in humans and livestock. In India, ST772 is a dominant and ST672 is an emerging clone of Staphylococcus aureus. Both cause serious human diseases, and carry type V SCCmec elements. The objective of this study was to characterize SCCmec type V elements of ST772 and ST672 because the usual PCR methods did not amplify all primers specific to the type. Whole genome sequencing analysis of seven ST772 and one ST672 S. aureus isolates revealed that the SCCmec elements of six of the ST772 isolates were the smallest of the extant type V elements and in addition have several other novel features. Only one ST772 isolate and the ST672 isolate carried bigger SCCmec cassettes which were composites carrying multiple ccrC genes. These cassettes had some similarities to type V SCCmec element from M013 isolate (ST59) from Taiwan in certain aspects. SCCmec elements of all Indian isolates had an inversion of the mec complex, similar to the bovine SCCmec type X. This study reveals that six out of seven ST772 S. aureus isolates have a novel type V (5C2) SCCmec element while one each of ST772 and ST672 isolates have a composite SCCmec type V element (5C2&5) formed by the integration of type V SCCmec into a MSSA carrying a SCC element, in addition to the mec gene complex inversions and extensive recombinations.

Genome analysis: A new approach for visualization of sequence organization in genomes

In this article we describe and demonstrate the versatility of a computer program, GENOME MAPPING, that uses interactive graphics and runs on an IRIS workstation. The program helps to visualize as well as analyse global and local patterns of genomic DNA sequences. It was developed keeping in mind the requirements of the human genome sequencing programme, which requires rapid analysis of the data. Using GENOME MAPPING one can discern signature patterns of different kinds of sequences and analyse such patterns for repetitive as well as rare sequence strings. Further, one can visualize the extent of global homology between different genomic sequences. An application of our method to the published yeast mitochondrial genome data shows similar sequence organizations in the entire sequence and in smaller subsequences.

Haemophilus influenzae Genome Database (HIGDB): A single point web resource for Haemophilus influenzae

Background: Haemophilus influenzae (H. Influenzae) is the causative agent of pneumonia, bacteraemia and meningitis. The organism is responsible for large number of deaths in both developed and developing countries. Even-though the first bacterial genome to be sequenced was that of H. Influenzae, there is no exclusive database dedicated for H. Influenzae. This prompted us to develop the Haemophilus influenzae Genome Database (HIGDB). Methods: All data of HIGDB are stored and managed in MySQL database. The HIGDB is hosted on Solaris server and developed using PERL modules. Ajax and JavaScript are used for the interface development. Results: The HIGDB contains detailed information on 42,741 proteins, 18,077 genes including 10 whole genome sequences and also 284 three dimensional structures of proteins of H. influenzae. In addition, the database provides ``Motif search'' and ``GBrowse''. The HIGDB is freely accessible through the URL:http://bioserverl.physicslisc.ernetin/HIGDB/. Discussion: The HIGDB will be a single point access for bacteriological, clinical, genomic and proteomic information of H. influenzae. The database can also be used to identify DNA motifs within H. influenzae genomes and to compare gene or protein sequences of a particular strain with other strains of H. influenzae. (C) 2014 Elsevier Ltd. All rights reserved.

Evolution of bacterial transcription factors: how proteins take on new tasks, but do not always stop doing the old ones

Many bacterial transcription factors do not behave as per the textbook operon model. We draw on whole genome work, as well as reported diversity across different bacteria, to argue that transcription factors may have evolved from nucleoid-associated proteins. This view would explain a large amount of recent data gleaned from high-throughput sequencing and bioinformatic analyses.

Diversity of DNA methyltransferases that recognize asymmetric target sequences

DNA methyltransferases (MTases) are a group of enzymes that catalyze the methyl group transfer from S-adenosyl-L-methionine in a sequence-specific manner. Orthodox Type II DNA MTases usually recognize palindromic DNA sequences and add a methyl group to the target base (either adenine or cytosine) on both strands. However, there are a number of MTases that recognize asymmetric target sequences and differ in their subunit organization. In a bacterial cell, after each round of replication, the substrate for any MTase is hemimethylated DNA, and it therefore needs only a single methylation event to restore the fully methylated state. This is in consistent with the fact that most of the DNA MTases studied exist as monomers in solution. Multiple lines of evidence suggest that some DNA MTases function as dimers. Further, functional analysis of many restriction-modification systems showed the presence of more than one or fused MTase genes. It was proposed that presence of two MTases responsible for the recognition and methylation of asymmetric sequences would protect the nascent strands generated during DNA replication from cognate restriction endonuclease. In this review, MTases recognizing asymmetric sequences have been grouped into different subgroups based on their unique properties. Detailed characterization of these unusual MTases would help in better understanding of their specific biological roles and mechanisms of action. The rapid progress made by the genome sequencing of bacteria and archaea may accelerate the identification and study of species- and strain-specific MTases of host-adapted bacteria and their roles in pathogenic mechanisms.

Diversity of DNA methyltransferases that recognize asymmetric target sequences

DNA methyltransferases (MTases) are a group of enzymes that catalyze the methyl group transfer from S-adenosyl-L-methionine in a sequence-specific manner. Orthodox Type II DNA MTases usually recognize palindromic DNA sequences and add a methyl group to the target base (either adenine or cytosine) on both strands. However, there are a number of MTases that recognize asymmetric target sequences and differ in their subunit organization. In a bacterial cell, after each round of replication, the substrate for any MTase is hemimethylated DNA, and it therefore needs only a single methylation event to restore the fully methylated state. This is in consistent with the fact that most of the DNA MTases studied exist as monomers in solution. Multiple lines of evidence suggest that some DNA MTases function as dimers. Further, functional analysis of many restriction-modification systems showed the presence of more than one or fused MTase genes. It was proposed that presence of two MTases responsible for the recognition and methylation of asymmetric sequences would protect the nascent strands generated during DNA replication from cognate restriction endonuclease. In this review, MTases recognizing asymmetric sequences have been grouped into different subgroups based on their unique properties. Detailed characterization of these unusual MTases would help in better understanding of their specific biological roles and mechanisms of action. The rapid progress made by the genome sequencing of bacteria and archaea may accelerate the identification and study of species- and strain-specific MTases of host-adapted bacteria and their roles in pathogenic mechanisms.

A homozygous mutation in LTBP2 causes isolated microspherophakia

Microspherophakia is an autosomal-recessive congenital disorder characterized by small spherical lens. It may be isolated or occur as part of a hereditary systemic disorder, such as Marfan syndrome, autosomal dominant and recessive forms of Weill-Marchesani syndrome, autosomal dominant glaucoma–lens ectopia–microspherophakia–stiVness– shortness syndrome, autosomal dominant microspherophakia with hernia, and microspherophakia-metaphyseal dysplasia. The purpose of this study was to map and identify the gene for isolated microspherophakia in two consanguineous Indian families. Using a whole-genome linkage scan in one family, we identiWed a likely locus for microspherophakia (MSP1) on chromosome 14q24.1–q32.12 between markers D14S588 and D14S1050 in a physical distance of 22.76 Mb. The maximum multi-point lod score was 2.91 between markers D14S1020 and D14S606. The MSP1 candidate region harbors 110 reference genes. DNA sequence analysis of one of the genes, LTBP2, detected a homozygous duplication (insertion) mutation, c.5446dupC, in the last exon (exon 36) in aVected family members. This homozygous mutation is predicted to elongate the LTBP2 protein by replacing the last 6 amino acids with 27 novel amino acids. Microspherophakia in the second family did not map to this locus, suggesting genetic heterogeneity. The present study suggests a role for LTBP2 in the structural stability of ciliary zonules, and growth and development of lens.

A homozygous mutation in LTBP2 causes isolated microspherophakia

Microspherophakia is an autosomal-recessive congenital disorder characterized by small spherical lens. It may be isolated or occur as part of a hereditary systemic disorder, such as Marfan syndrome, autosomal dominant and recessive forms of Weill-Marchesani syndrome, autosomal dominant glaucoma-lens ectopia-microspherophakia-stiffness-shortness syndrome, autosomal dominant microspherophakia with hernia, and microspherophakia-metaphyseal dysplasia. The purpose of this study was to map and identify the gene for isolated microspherophakia in two consanguineous Indian families. Using a whole-genome linkage scan in one family, we identified a likely locus for microspherophakia (MSP1) on chromosome 14q24.1-q32.12 between markers D14S588 and D14S1050 in a physical distance of 22.76 Mb. The maximum multi-point lod score was 2.91 between markers D14S1020 and D14S606. The MSP1 candidate region harbors 110 reference genes. DNA sequence analysis of one of the genes, LTBP2, detected a homozygous duplication (insertion) mutation, c.5446dupC, in the last exon (exon 36) in affected family members. This homozygous mutation is predicted to elongate the LTBP2 protein by replacing the last 6 amino acids with 27 novel amino acids. Microspherophakia in the second family did not map to this locus, suggesting genetic heterogeneity. The present study suggests a role for LTBP2 in the structural stability of ciliary zonules, and growth and development of lens.