GIV: A Tool for Genomic Islands Visualization

A Genomic Islands (GI) is a chunk of DNA sequence in a genome whose origin can be traced back to other organisms or viruses. The detection of GIs plays an indispensable role in biomedical research, due to the fact that GIs are highly related to special functionalities such as disease-causing GIs - pathogenicity islands. It is also very important to visualize genomic islands, as well as the supporting features corresponding to the genomic islands in the genome. We have developed a program, Genomic Island Visualization (GIV), which displays the locations of genomic islands in a genome, as well as the corresponding supportive feature information for GIs. GIV was implemented in C++, and was compiled and executed on Linux/Unix operating systems.

Availability

GIV is freely available for non-commercial use at http://www5.esu.edu/cpsc/bioinfo/software/GIV     

Insect Barcode Information System

Insect Barcode Information System called as Insect Barcode Informática (IBIn) is an online database resource developed by the National Bureau of Agriculturally Important Insects, Bangalore. This database provides acquisition, storage, analysis and publication of DNA barcode records of agriculturally important insects, for researchers specifically in India and other countries. It bridges a gap in bioinformatics by integrating molecular, morphological and distribution details of agriculturally important insects. IBIn was developed using PHP/My SQL by using relational database management concept. This database is based on the client– server architecture, where many clients can access data simultaneously. IBIn is freely available on-line and is user-friendly. IBIn allows the registered users to input new information, search and view information related to DNA barcode of agriculturally important insects.This paper provides a current status of insect barcode in India and brief introduction about the database IBIn.

Availability

http://www.nabg-nbaii.res.in/barcode     

FMiR: A Curated Resource of Mitochondrial DNA Information for Fish

Mitochondrial genome sequences have been widely used for evolutionary and phylogenetic studies. Among vertebrates, fish are an important, diverse group, and their mitogenome sequences are growing rapidly in public repositories. To facilitate mitochondrial genome analysis and to explore the valuable genetic information, we developed the Fish Mitogenome Resource (FMiR) database to provide a workbench for mitogenome annotation, species identification and microsatellite marker mining. The microsatellites are also known as simple sequence repeats (SSRs) and used as molecular markers in studies on population genetics, gene duplication and marker assisted selection. Here, easy-to-use tools have been implemented for mining SSRs and for designing primers to identify species/habitat specific markers. In addition, FMiR can analyze complete or partial mitochondrial genome sequence to identify species and to deduce relational distances among sequences across species. The database presently contains curated mitochondrial genomes from 1302 fish species belonging to 297 families and 47 orders reported from saltwater and freshwater ecosystems. In addition, the database covers information on fish species such as conservation status, ecosystem, family, distribution and occurrence downloaded from the FishBase and IUCN Red List databases. Those fish information have been used to browse mitogenome information for the species belonging to a particular category. The database is scalable in terms of content and inclusion of other analytical modules. The FMiR is running under Linux operating platform on high performance server accessible at URL http://mail.nbfgr.res.in/fmir.     

NABIC marker database: A molecular markers information network of agricultural crops

In 2013, National Agricultural Biotechnology Information Center (NABIC) reconstructs a molecular marker database for useful genetic resources. The web-based marker database consists of three major functional categories: map viewer, RSN marker and gene annotation. It provides 7250 marker locations, 3301 RSN marker property, 3280 molecular marker annotation information in agricultural plants. The individual molecular marker provides information such as marker name, expressed sequence tag number, gene definition and general marker information. This updated marker-based database provides useful information through a user-friendly web interface that assisted in tracing any new structures of the chromosomes and gene positional functions using specific molecular markers.

Availability

The database is available for free at http://nabic.rda.go.kr/gere/rice/molecularMarkers/     

Carbohydrate Metabolism Differences between Subgroup A1 and B2 Strains of Bacillus anthracis as Assessed by Comparative Genomics and Functional Genetics

In France, Bacillus anthracis subgroup B2 strains do not metabolize starch or glycogen but can use gluconate, whereas subgroup A1 strains show the inverse pattern. Functional genetic analysis revealed that mutations in the amyS and gntK genes encoding an alpha-amylase and a gluconate kinase, respectively, were responsible for these phenotypes.Bacillus anthracis, the etiological agent of anthrax, is a gram-positive, aerobic soil bacterium. Multilocus variable-number tandem repeat analysis of a collection of French isolates shows that the main groups of B. anthracis groups A (subgroup A1) and B (subgroup B2) described worldwide are represented (1, 2). Subgroup B2 isolates are the most common isolates in France and are found particularly in southern mountain regions, but they are extremely rare elsewhere in the world. Biochemical characterization of French isolates indicates that subgroup A1 and B2 strains have different carbohydrate utilization patterns (P. Vaissaire, A. Fouet, K. L. Smith, C. Keys, C. Le Doujet, P. Sylvestre, M. Levy, P. Keim, and M. Mock, presented at the 5th International Conference on Anthrax and 3rd International Workshop on the Molecular Biology of Bacillus cereus, B. anthracis and B. thuringiensis, 30 March to 3 April 2003, Nice, France). French subgroup A1 strains metabolize starch and glycogen but not gluconate, and the inverse is true for subgroup B2 strains. The genomes of several B. anthracis strains are available on the NCBI website (http://www.ncbi.nlm.nih.gov/), and two of these strains, Ames and CNEVA, are representative of groups A and B, respectively. We compared the genomic sequences of Ames and CNEVA to identify mutations that may affect metabolic activities involved in the phenotypic differences.The Kegg pathway database (http://www.genome.jp/kegg/pathway.html) was used to select enzyme activities involved in the metabolic pathways for starch, glycogen, and gluconate. BLAST analysis of the corresponding open reading frame in the Ames (subgroup A3) and CNEVA (subgroup B2) genomes was then used to identify the selected genes that were interrupted or mutated. The functions and localizations of these open reading frames were then investigated with the Pfam (http://pfam.sanger.ac.uk/), CDD (http://www.ncbi.nlm.nih.gov/Structure/cdd/cdd.shtml), SMART (http://smart.embl-heidelberg.de/), SignalP (http://www.cbs.dtu.dk/services/SignalP/), and TMHMM (http://www.cbs.dtu.dk/services/TMHMM-2.0/) search programs. A single-base deletion in the amyS gene (BA3551) encoding an alpha-amylase linked to starch and glycogen metabolism was found in the CNEVA genome. The wild-type AmyS protein contains 513 amino acids, and its predicted molecular mass is 58.4 kDa. In subgroup B2, there is a frameshift due to deletion of an adenosine in the 7th position of the nucleotide sequence that leads to a premature stop codon in the 13th position. In the Ames genome, a single-base substitution was found in the gntK gene (BA0162) encoding a gluconate kinase linked to gluconate metabolism. The predicted wild-type GntK protein contains 511 amino acids, and its predicted molecular mass is 56.7 kDa. The mutation identified is a cytosine-to-adenosine substitution at position 530 of the nucleotide sequence that leads to a premature stop codon at amino acid position 176. We confirmed the presence of these two mutations in the other B. anthracis subgroup genomes accessible in the NCBI unfinished microbial genome database and sequenced 12 isolates with various genotypes belonging to subgroups A1 and B2 (6 isolates in each subgroup) originating from outbreaks that occurred in different regions of France over the last 15 years. These analyses revealed that the deletion in amyS is restricted to strains belonging to group B subgroups, whereas the substitution in gntK is restricted to strains belonging to group A subgroups. The mutations identified in amyS and gntK both result in premature stop codons that lead to a loss of the enzymatic activities and may thus account for the observed phenotypic differences between subgroup A1 and B2 strains. We therefore focused on these two genes and used French strains 9602R and RA3R belonging to subgroups A1 and B2, respectively, for further analysis.     

A database of natural products and chemical entities from marine habitat

Marine compound database consists of marine natural products and chemical entities, collected from various literature sources, which are known to possess bioactivity against human diseases. The database is constructed using html code. The 12 categories of 182 compounds are provided with the source, compound name, 2-dimensional structure, bioactivity and clinical trial information. The database is freely available online and can be accessed at http://www.progenebio.in/mcdb/index.htm     

EMu: probabilistic inference of mutational processes and their localization in the cancer genome

The spectrum of mutations discovered in cancer genomes can be explained by the activity of a few elementary mutational processes. We present a novel probabilistic method, EMu, to infer the mutational signatures of these processes from a collection of sequenced tumors. EMu naturally incorporates the tumor-specific opportunity for different mutation types according to sequence composition. Applying EMu to breast cancer data, we derive detailed maps of the activity of each process, both genome-wide and within specific local regions of the genome. Our work provides new opportunities to study the mutational processes underlying cancer development. EMu is available at http://www.sanger.ac.uk/resources/software/emu/.     

Development of eSSR-Markers in Setaria italica and Their Applicability in Studying Genetic Diversity,Cross-Transferability and Comparative Mapping in Millet and Non-Millet Species

Foxtail millet ( Setaria italica L.) is a tractable experimental model crop for studying functional genomics of millets and bioenergy grasses. But the limited availability of genomic resources, particularly expressed sequence-based genic markers is significantly impeding its genetic improvement. Considering this, we attempted to develop EST-derived-SSR (eSSR) markers and utilize them in germplasm characterization, cross-genera transferability and in silico comparative mapping. From 66,027 foxtail millet EST sequences 24,828 non-redundant ESTs were deduced, representing ~16 Mb, which revealed 534 (~2%) eSSRs in 495 SSR containing ESTs at a frequency of 1/30 kb. A total of 447 pp were successfully designed, of which 327 were mapped physically onto nine chromosomes. About 106 selected primer pairs representing the foxtail millet genome showed high-level of cross-genera amplification at an average of ~88% in eight millets and four non-millet species. Broad range of genetic diversity (0.02–0.65) obtained in constructed phylogenetic tree using 40 eSSR markers demonstrated its utility in germplasm characterizations and phylogenetics. Comparative mapping of physically mapped eSSR markers showed considerable proportion of sequence-based orthology and syntenic relationship between foxtail millet chromosomes and sorghum (~68%), maize (~61%) and rice (~42%) chromosomes. Synteny analysis of eSSRs of foxtail millet, rice, maize and sorghum suggested the nested chromosome fusion frequently observed in grass genomes. Thus, for the first time we had generated large-scale eSSR markers in foxtail millet and demonstrated their utility in germplasm characterization, transferability, phylogenetics and comparative mapping studies in millets and bioenergy grass species.     

Development of 5123 Intron-Length Polymorphic Markers for Large-Scale Genotyping Applications in Foxtail Millet

Generating genomic resources in terms of molecular markers is imperative in molecular breeding for crop improvement. Though development and application of microsatellite markers in large-scale was reported in the model crop foxtail millet, no such large-scale study was conducted for intron-length polymorphic (ILP) markers. Considering this, we developed 5123 ILP markers, of which 4049 were physically mapped onto 9 chromosomes of foxtail millet. BLAST analysis of 5123 expressed sequence tags (ESTs) suggested the function for ∼71.5% ESTs and grouped them into 5 different functional categories. About 440 selected primer pairs representing the foxtail millet genome and the different functional groups showed high-level of cross-genera amplification at an average of ∼85% in eight millets and five non-millet species. The efficacy of the ILP markers for distinguishing the foxtail millet is demonstrated by observed heterozygosity (0.20) and Nei''s average gene diversity (0.22). In silico comparative mapping of physically mapped ILP markers demonstrated substantial percentage of sequence-based orthology and syntenic relationship between foxtail millet chromosomes and sorghum (∼50%), maize (∼46%), rice (∼21%) and Brachypodium (∼21%) chromosomes. Hence, for the first time, we developed large-scale ILP markers in foxtail millet and demonstrated their utility in germplasm characterization, transferability, phylogenetics and comparative mapping studies in millets and bioenergy grass species.     

Remarkably Divergent Regions Punctuate the Genome Assembly of the Caenorhabditis elegans Hawaiian Strain CB4856

The Hawaiian strain (CB4856) of Caenorhabditis elegans is one of the most divergent from the canonical laboratory strain N2 and has been widely used in developmental, population, and evolutionary studies. To enhance the utility of the strain, we have generated a draft sequence of the CB4856 genome, exploiting a variety of resources and strategies. When compared against the N2 reference, the CB4856 genome has 327,050 single nucleotide variants (SNVs) and 79,529 insertion–deletion events that result in a total of 3.3 Mb of N2 sequence missing from CB4856 and 1.4 Mb of sequence present in CB4856 but not present in N2. As previously reported, the density of SNVs varies along the chromosomes, with the arms of chromosomes showing greater average variation than the centers. In addition, we find 61 regions totaling 2.8 Mb, distributed across all six chromosomes, which have a greatly elevated SNV density, ranging from 2 to 16% SNVs. A survey of other wild isolates show that the two alternative haplotypes for each region are widely distributed, suggesting they have been maintained by balancing selection over long evolutionary times. These divergent regions contain an abundance of genes from large rapidly evolving families encoding F-box, MATH, BATH, seven-transmembrane G-coupled receptors, and nuclear hormone receptors, suggesting that they provide selective advantages in natural environments. The draft sequence makes available a comprehensive catalog of sequence differences between the CB4856 and N2 strains that will facilitate the molecular dissection of their phenotypic differences. Our work also emphasizes the importance of going beyond simple alignment of reads to a reference genome when assessing differences between genomes.     

CytokineDB: a database collecting biological information

Cytokines are subdivided in 12 sub-families and are described as multi-functional molecules that play an important biological activity in host defense system against pathogens, in homeostasis, tissue repair, cell growth and development. CytokineDB is an annotated database that collects biological information regarding the cytokines family in human and will be periodically updated by including new biological information. This database is freely available online and can be accessed at the URL: http://www.cro-m.eu/CytokineDB/     

Replisome Function During Replicative Stress Is Modulated by Histone H3 Lysine 56 Acetylation Through Ctf4

Histone H3 lysine 56 acetylation in Saccharomyces cerevisiae is required for the maintenance of genome stability under normal conditions and upon DNA replication stress. Here we show that in the absence of H3 lysine 56 acetylation replisome components become deleterious when replication forks collapse at natural replication block sites. This lethality is not a direct consequence of chromatin assembly defects during replication fork progression. Rather, our genetic analyses suggest that in the presence of replicative stress H3 lysine 56 acetylation uncouples the Cdc45–Mcm2-7–GINS DNA helicase complex and DNA polymerases through the replisome component Ctf4. In addition, we discovered that the N-terminal domain of Ctf4, necessary for the interaction of Ctf4 with Mms22, an adaptor protein of the Rtt101-Mms1 E3 ubiquitin ligase, is required for the function of the H3 lysine 56 acetylation pathway, suggesting that replicative stress promotes the interaction between Ctf4 and Mms22. Taken together, our results indicate that Ctf4 is an essential member of the H3 lysine 56 acetylation pathway and provide novel mechanistic insights into understanding the role of H3 lysine 56 acetylation in maintaining genome stability upon replication stress.     

RADtyping: An Integrated Package for Accurate De Novo Codominant and Dominant RAD Genotyping in Mapping Populations

Genetic linkage maps are indispensable tools in genetic, genomic and breeding studies. As one of genotyping-by-sequencing methods, RAD-Seq (restriction-site associated DNA sequencing) has gained particular popularity for construction of high-density linkage maps. Current RAD analytical tools are being predominantly used for typing codominant markers. However, no genotyping algorithm has been developed for dominant markers (resulting from recognition site disruption). Given their abundance in eukaryotic genomes, utilization of dominant markers would greatly diminish the extensive sequencing effort required for large-scale marker development. In this study, we established, for the first time, a novel statistical framework for de novo dominant genotyping in mapping populations. An integrated package called RADtyping was developed by incorporating both de novo codominant and dominant genotyping algorithms. We demonstrated the superb performance of RADtyping in achieving remarkably high genotyping accuracy based on simulated and real mapping datasets. The RADtyping package is freely available at http://www2.ouc.edu.cn/mollusk/ detailen.asp?id=727.     

A low-latency,big database system and browser for storage,querying and visualization of 3D genomic data

Recent releases of genome three-dimensional (3D) structures have the potential to transform our understanding of genomes. Nonetheless, the storage technology and visualization tools need to evolve to offer to the scientific community fast and convenient access to these data. We introduce simultaneously a database system to store and query 3D genomic data (3DBG), and a 3D genome browser to visualize and explore 3D genome structures (3DGB). We benchmark 3DBG against state-of-the-art systems and demonstrate that it is faster than previous solutions, and importantly gracefully scales with the size of data. We also illustrate the usefulness of our 3D genome Web browser to explore human genome structures. The 3D genome browser is available at http://3dgb.cs.mcgill.ca/.     

Phyto-mellitus: A phyto-chemical database for diabetes

Herbs are the base used for treatment in Ayurveda. We describe a database named Phyto-Mellitus with information on plants traditionally used for diabetes with their chemical constituents. The active principles of these plants are antioxidant and free radical scavenging.

Availability

http://www.bicmlacw.org/bt/     

SomatiCA: Identifying,Characterizing and Quantifying Somatic Copy Number Aberrations from Cancer Genome Sequencing Data

Whole genome sequencing of matched tumor-normal sample pairs is becoming routine in cancer research. However, analysis of somatic copy-number changes from sequencing data is still challenging because of insufficient sequencing coverage, unknown tumor sample purity and subclonal heterogeneity. Here we describe a computational framework, named SomatiCA, which explicitly accounts for tumor purity and subclonality in the analysis of somatic copy-number profiles. Taking read depths (RD) and lesser allele frequencies (LAF) as input, SomatiCA will output 1) admixture rate for each tumor sample, 2) somatic allelic copy-number for each genomic segment, 3) fraction of tumor cells with subclonal change in each somatic copy number aberration (SCNA), and 4) a list of substantial genomic aberration events including gain, loss and LOH. SomatiCA is available as a Bioconductor R package at http://www.bioconductor.org/packages/2.13/bioc/html/SomatiCA.html.