Lightweight genome viewer: portable software for browsing genomics data in its chromosomal context

An information system on genomic elements and predicted protein structures of buffalo (Bubalus bubalis)

Among the livestock species, buffalo remained as an integral part of the Indian rural economy. With the advent of genome sequencing technologies, it became possible to sequence the whole genome of Murrah buffalo. Also, significant amount of information on different genomic elements of buffalo is available at National Centre for Biotechnology Information (NCBI). However, the positions of these elements on the genome are not fully known. In addition, the 3D structures of buffalo proteins are not available and also there exist no browser to visualize important genic elements on buffalo genome. Hence, a study was taken up to develop a web-based information system having information on genomic elements, protein 3-D structures and genome browser. Initially, information on nucleotide and protein sequences were retrieved from NCBI and parsed suitably. Later, the protein structures were predicted, validated, refined and stabilized in silico. An Information System on Buffalo Genome (ISBG) with 3-tier architecture was developed containing the sequence and structural information. ISBG contains complete coding sequences (CDS), Mitochondrial DNAs, 1k upstream regions and Untranslated Regions (UTRs) of buffalo genome. The buffalo genes were also mapped onto the genome. The results revealed that maximum number of genes were found distributed on chromosome 4 followed by chromosome 18, which can also be visualized from the developed genome browser. ISBG can be accessed at http://cabgrid.res.in:8080/bgis. The proposed information system helps animal breeders and biotechnologist in animal improvement.

GinMicrosatDb: a genome-wide microsatellite markers database for sesame (Sesamum indicum L.)

Molecular breeding in sesame is still at infancy due to limited number of microsatellite markers available and the low level of polymorphism exhibited by them. Therefore, whole genome sequencing was used for development of microsatellite markers so as to ensure availability of substantial number of polymorphic markers for use in marker assisted breeding programs. Whole genome sequencing of sesame variety 'Swetha' was done using Illumina paired-end sequencing and Roche 454 shotgun sequencing technologies (GCA_000975565.1 in GenBank). 'GinMicrosatDb', a genome-wide microsatellite marker database has been developed using the whole genome sequence data of sesame variety 'Swetha'. The database consists of microsatellites localized on both linkage groups and scaffolds with their genomic co-ordinates. It provides five sets of forward and reverse primers for each of the microsatellite loci along with the flanking sequences, primer GC content, product size and melting temperature etc. The distribution of microsatellites can be viewed and selected through a genome browser as well as through a physical map. The newly identified microsatellite markers are expected to help sesame breeders in developing marker tags for traits of economic importance thereby bringing about greater efficiency in marker-assisted selection programs.

Pronounced strain-specific chemosensory receptor gene expression in the mouse vomeronasal organ

Background

The chemosensory system plays an important role in orchestrating sexual behaviors in mammals. Pheromones trigger sexually dimorphic behaviors and different mouse strains exhibit differential responses to pheromone stimuli. It has been speculated that differential gene expression in the sensory organs that detect pheromones may underlie sexually-dimorphic and strain-specific responses to pheromone cues.

Results

We have performed transcriptome analyses of the mouse vomeronasal organ, a sensory organ recognizing pheromones and interspecies cues. We find little evidence of sexual dimorphism in gene expression except for Xist, an essential gene for X-linked gene inactivation. Variations in gene expression are found mainly among strains, with genes from immune response and chemosensory receptor classes dominating the list. Differentially expressed genes are concentrated in genomic hotspots enriched in these families of genes. Some chemosensory receptors show exclusive patterns of expression in different strains. We find high levels of single nucleotide polymorphism in chemosensory receptor pseudogenes, some of which lead to functionalized receptors. Moreover, we identify a number of differentially expressed long noncoding RNA species showing strong correlation or anti-correlation with chemoreceptor genes.

Conclusions

Our analyses provide little evidence supporting sexually dimorphic gene expression in the vomeronasal organ that may underlie dimorphic pheromone responses. In contrast, we find pronounced variations in the expression of immune response related genes, vomeronasal and G-protein coupled receptor genes among different mouse strains. These findings raised the possibility that diverse strains of mouse perceive pheromone cues differently and behavioral difference among strains in response to pheromone may first arise from differential detection of pheromones. On the other hand, sexually dimorphic responses to pheromones more likely originate from dimorphic neural circuits in the brain than from differential detection. Moreover, noncoding RNA may offer a potential regulatory mechanism controlling the differential expression patterns.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4364-4) contains supplementary material, which is available to authorized users.

Isolation and Characterization of a Novel Gammaherpesvirus from a Microbat Cell Line

While employing deep sequencing and de novo assembly to characterize the mRNA transcript profile of a cell line derived from the microbat Myotis velifer incautus, we serendipitously identified mRNAs encoding proteins with a high level of identity to herpesviruses. A majority were closely related to proteins of equine herpesvirus 2 (EHV-2), a horse gammaherpesvirus. We demonstrated by electron microscopy the presence of herpesvirus-like particles in the microbat cells. Passage of supernatants from microbat cells to Vero cells resulted in syncytium formation, and expression of viral genes and amplification of viral DNA were demonstrated by quantitative PCR. Susceptibility of human cell lines to productive infection was also demonstrated. Next-generation sequencing and de novo assembly of the viral genome from supernatants from Vero cells yielded a single contig of approximately 130 kb with at least 77 open reading frames (ORFs), predicted microRNAs (miRNAs), and a gammaherpesvirus genomic organization. Phylogenic analysis of the envelope glycoprotein (gB) and DNA polymerase (POLD1) revealed similarity to multiple gammaherpesviruses, including those from as-yet-uncultured viruses of the Rhadinovirus genus that were obtained by deep sequencing of bat tissues. Moreover, the assembled genome revealed ORFs that share little or no homology to known ORFs in EHV-2 but are similar to accessory proteins of other gammaherpesviruses. Some also have striking homology to predicted Myotis bat proteins. Cumulatively, this study provides the first isolation and characterization of a replication-competent bat gammaherpesvirus. IMPORTANCE Bats are of significant interest as reservoirs for zoonotic viral pathogens; however, tools to dissect bat-virus interactions are limited in availability. This study serendipitously identified, in an established bat cell line, a fully replication-competent gammaherpesvirus; determined the complete genome sequence of the virus; and generated a viral transcript map. This virus can replicate in select human and nonhuman primate cell lines. However, analyses of viral sequences support a bat origin for this virus; we therefore refer to the virus as bat gammaherpesvirus 8 (BGHV8). The viral genome contains unique open reading frames that likely encode modulators of bat innate and adaptive immune signaling pathways and expresses viral miRNAs. The virus and its gene products should provide a unique tool to dissect both bat and gammaherpesvirus biology.

GenPlay Multi-Genome, a tool to compare and analyze multiple human genomes in a graphical interface

SUMMARY

Parallel visualization of multiple individual human genomes is a complex endeavor that is rapidly gaining importance with the increasing number of personal, phased and cancer genomes that are being generated. It requires the display of variants such as SNPs, indels and structural variants that are unique to specific genomes and the introduction of multiple overlapping gaps in the reference sequence. Here, we describe GenPlay Multi-Genome, an application specifically written to visualize and analyze multiple human genomes in parallel. GenPlay Multi-Genome is ideally suited for the comparison of allele-specific expression and functional genomic data obtained from multiple phased genomes in a graphical interface with access to multiple-track operation. It also allows the analysis of data that have been aligned to custom genomes rather than to a standard reference and can be used as a variant calling format file browser and as a tool to compare different genome assembly, such as hg19 and hg38.

AVAILABILITY AND IMPLEMENTATION

GenPlay is available under the GNU public license (GPL-3) from http://genplay.einstein.yu.edu. The source code is available at https://github.com/JulienLajugie/GenPlay.

The shared antibiotic resistome of soil bacteria and human pathogens

Soil microbiota represent one of the ancient evolutionary origins of antibiotic resistance and have been proposed as a reservoir of resistance genes available for exchange with clinical pathogens. Using a high-throughput functional metagenomic approach in conjunction with a pipeline for the de novo assembly of short-read sequence data from functional selections (termed PARFuMS), we provide evidence for recent exchange of antibiotic resistance genes between environmental bacteria and clinical pathogens. We describe multidrug-resistant soil bacteria containing resistance cassettes against five classes of antibiotics (β-lactams, aminoglycosides, amphenicols, sulfonamides, and tetracyclines) that have perfect nucleotide identity to genes from diverse human pathogens. This identity encompasses noncoding regions as well as multiple mobilization sequences, offering not only evidence of lateral exchange but also a mechanism by which antibiotic resistance disseminates.

Geoseq: a tool for dissecting deep-sequencing datasets

Background

Datasets generated on deep-sequencing platforms have been deposited in various public repositories such as the Gene Expression Omnibus (GEO), Sequence Read Archive (SRA) hosted by the NCBI, or the DNA Data Bank of Japan (ddbj). Despite being rich data sources, they have not been used much due to the difficulty in locating and analyzing datasets of interest.

Results

Geoseq http://geoseq.mssm.edu provides a new method of analyzing short reads from deep sequencing experiments. Instead of mapping the reads to reference genomes or sequences, Geoseq maps a reference sequence against the sequencing data. It is web-based, and holds pre-computed data from public libraries. The analysis reduces the input sequence to tiles and measures the coverage of each tile in a sequence library through the use of suffix arrays. The user can upload custom target sequences or use gene/miRNA names for the search and get back results as plots and spreadsheet files. Geoseq organizes the public sequencing data using a controlled vocabulary, allowing identification of relevant libraries by organism, tissue and type of experiment.

Conclusions

Analysis of small sets of sequences against deep-sequencing datasets, as well as identification of public datasets of interest, is simplified by Geoseq. We applied Geoseq to, a) identify differential isoform expression in mRNA-seq datasets, b) identify miRNAs (microRNAs) in libraries, and identify mature and star sequences in miRNAS and c) to identify potentially mis-annotated miRNAs. The ease of using Geoseq for these analyses suggests its utility and uniqueness as an analysis tool.

Viruses in the faecal microbiota of monozygotic twins and their mothers

Viral diversity and life cycles are poorly understood in the human gut and other body habitats. Phages and their encoded functions may provide informative signatures of a human microbiota and of microbial community responses to various disturbances, and may indicate whether community health or dysfunction is manifest after apparent recovery from a disease or therapeutic intervention. Here we report sequencing of the viromes (metagenomes) of virus-like particles isolated from faecal samples collected from healthy adult female monozygotic twins and their mothers at three time points over a one-year period. We compared these data sets with data sets of sequenced bacterial 16S ribosomal RNA genes and total-faecal-community DNA. Co-twins and their mothers share a significantly greater degree of similarity in their faecal bacterial communities than do unrelated individuals. In contrast, viromes are unique to individuals regardless of their degree of genetic relatedness. Despite remarkable interpersonal variations in viromes and their encoded functions, intrapersonal diversity is very low, with >95% of virotypes retained over the period surveyed, and with viromes dominated by a few temperate phages that exhibit remarkable genetic stability. These results indicate that a predatory viral-microbial dynamic, manifest in a number of other characterized environmental ecosystems, is notably absent in the very distal intestine.

Many Microbe Microarrays Database: uniformly normalized Affymetrix compendia with structured experimental metadata

Many Microbe Microarrays Database (M^3D) is designed to facilitate the analysis and visualization of expression data in compendia compiled from multiple laboratories. M^3D contains over a thousand Affymetrix microarrays for Escherichia coli, Saccharomyces cerevisiae and Shewanella oneidensis. The expression data is uniformly normalized to make the data generated by different laboratories and researchers more comparable. To facilitate computational analyses, M^3D provides raw data (CEL file) and normalized data downloads of each compendium. In addition, web-based construction, visualization and download of custom datasets are provided to facilitate efficient interrogation of the compendium for more focused analyses. The experimental condition metadata in M^3D is human curated with each chemical and growth attribute stored as a structured and computable set of experimental features with consistent naming conventions and units. All versions of the normalized compendia constructed for each species are maintained and accessible in perpetuity to facilitate the future interpretation and comparison of results published on M^3D data. M^3D is accessible at http://m3d.bu.edu/.