Prakash G. Koringa

Metagenomic analysis of virulence-associated and antibiotic resistance genes of microbes in rumen of Indian buffalo (Bubalus bubalis)

A major research goal in rumen microbial ecology is to understand the relationship between community composition and its function, particularly involved in fermentation process is of a potential interest. The buffalo rumen microbiota impacts human food safety as well as animal health. Although the bacteria of bovine rumen have been well characterized, techniques have been lacking to correlate total community structure with gene function. We applied 454 next generations sequencing technology to characterize general microbial diversity present in buffalo rumen metagenome and also identified the repertoire of microbial genes present, including genes associated with antibiotic resistance and bacterial virulence. Results suggest that over six percent (6.44%) of the sequences from our buffalo rumen pool sample could be categorized as virulence genes and genes associated with resistance to antibiotic and toxic compounds (RATC), which is a higher proportion of virulence genes reported from metagenome samples of chicken cecum (5.39%), cow rumen (4.43%) and Sargasso sea (2.95%). However, it was lower than the proportion found in cow milk (11.33%) cattle faeces (8.4%), Antarctic marine derived lake (8.45%), human fecal (7.7%) and farm soil (7.79%). The dynamic nature of metagenomic data, together with the large number of RATC classes observed in samples from widely different ecologies indicates that metagenomic data can be used to track potential targets and relative amounts of antibiotic resistance genes in individual animals. In addition, these data can be also used to generate antibiotic resistance gene profiles to facilitate an understanding of the ecology of the microbial communities in each habitat as well as the epidemiology of antibiotic resistant gene transport between and among habitats.

Genome Sequence of Pasteurella multocida subsp. gallicida Anand1_poultry

We report the finished and annotated genome sequence of Pasteurella multocida gallicida strain Anand1_poultry, which was isolated from the liver of a diseased adult female chicken. The strain causes a disease called fowl cholera, which is a contagious disease in birds. We compared it with the published genome sequence of Pasteurella multocida Pm70.

Identification of novel transcripts deregulated in buccal cancer by RNA-seq

The differential transcriptome analysis provides better understanding of molecular pathways leading to cancer, which in turn allows designing the effective strategies for diagnosis, therapeutic intervention and prediction of therapeutic outcome. This study describes the transcriptome analysis of buccal cancer and normal tissue by CLC Genomics Workbench from the data generated by Roches 454 sequencing platform, which identified total of 1797 and 2655 genes uniquely expressed in normal and cancer tissues, respectively with 2466 genes expressed in both tissues. Among the genes expressed in both tissues, 1842 were up-regulated whereas 624 were down-regulated in cancer tissue. Besides transcripts known to be involved in cancer, this study led to the identification of novel transcripts, with significantly altered expression in buccal cancer tissue, providing potential targets for diagnosis and cancer therapeutics. The functional categorization by the KEGG pathway and gene ontology analysis revealed enrichment of differentially expressed transcripts to various pathways leading to cancer, including the p53 signaling pathway. Moreover, the gene ontology analysis unfolded suppression of transcripts involved in actin mediated cell contraction process. The down-regulation of four of these transcripts MYL1, ACTA1, TCAP and DESMIN in buccal cancer were further supported by quantitative PCR signifying its possible implication in the cancer progression.

RNA-Seq reveals differentially expressed isoforms and novel splice variants in buccal mucosal cancer.

Buccal mucosal cancer (BMC) is a multifactorial disease with poorly defined genetic profile and prognosis due to late detection stage and unavailability of reliable prognostic markers. To identify aberrant transcriptional events, we employed high throughput RNA-Seq analysis of BMC and normal tissue. Comparative transcriptome analysis with Cufflinks revealed 260 up and 328 down regulated genes whereas, 350 up and 397 down regulated isoforms by at least two folds over buccal normal in BMC. Study revealed 46 splice variants in normal and 106 in cancer, out of which 10 variants were validated with end point RT-PCR. Expression of two isoforms of CD74 was validated using RT-qPCR and found in accordance with RNA-Seq. Further extensive follow up analysis of modulator genes, isoforms and splice variants found in this study, might be useful in deep understanding of pathological changes in BMC and development of prospective intervention strategies.

Transcriptome analysis and SNP Identification in SCC of Horn in (Bos indicus) Indian cattle

Single Nucleotide Polymorphisms (SNPs) have become the marker of choice for genome wide association studies. In order to provide the best genome coverage for the analysis of disease, production and performance traits, a large number of relatively evenly distributed SNPs are needed. The main objective of present work was to identify large numbers of gene-associated SNPs using high-throughput sequencing in squamous cell carcinoma of horn. RNA-seq analysis was conducted on 2 tissues viz. Horn Cancer (HC) and Horn Normal (HN) in Kankrej breed of cattle. A total of 909,362 reads with average read length of 405 bp for HC and 583,491 reads with average read length of 411 bp for HN were obtained. We found 9532 and 7065 SNPs as well as 1771 and 1172 Indels in HC and HN, respectively, from which, 7889 SNPs and 1736 Indels were uniquely present in HC, 5886 SNPs and 1146 Indels were uniquely present in HN and reported first time in Bos indicus, whereas the rest are already reported in Bos taurus dbSNP database. The gene-associated SNPs and Indels were high in upregulated genes of HC as compared to HN. Analysis of differentially expressed genes was identified, these genes are involved in regulation of cell proliferation, apoptosis, gene transcription, cell survival and metabolism through various metabolic pathways. The result of transcriptome expression profiling was validated using Real Time quantitative PCR in nine randomly selected genes. We identified numbers aberrant signaling pathways responsible for carcinogenesis in HC which are also commonly altered in squamous cell carcinoma (SCC) of lung in human being. We conclude that a large number of altered genes and dysfunction of multiple pathways are involved in the development of Horn Cancer. The present findings contribute to theoretical information for further screening of genes and identification of markers for early diagnosis of HC as well as SNPs identified in this report provide a much needed resource for genetic studies in B. indicus and shall contribute to the development of a high density SNP array. Validation and testing of these SNPs using SNP arrays will form the material basis for gene associated SNPs in HC.

A preliminary sketch of horn cancer transcriptome in Indian zebu cattle

Horn cancer, a type of squamous cell carcinoma, in zebu cattle is an expensive affair in Indian agriculture sector, which accounts for 83.34% of total tumors found. In general, cancer tissue confirms considerably different expression patterns when compared to a normal stage. This includes not only up/down regulation, but also, the aberrant gene expression, the presence of different non-coding RNAs (ncRNAs), pseudogenes expression and genes involved in unusual pathways. We employed Roche 454 next generation sequencing platform to sequence Bos indicus cancerous and normal horn tissue transcripts. This resulted into a total of 909,345 high-confidence deep sequencing reads and detected a range of unusual transcriptional events including tumor associated genes. We also validated expression of two of the four tested genes in five other similar tissue samples by RT-qPCR. Further, seven cancer specific non-coding transcripts were accessed and a few of them have been suggested as cancer specific markers. This study for the first time provides primary transcriptome sketch of Bos indicus horn cancer tissue, and also demonstrates the suitability of the 454 sequencer for transcriptome analysis, which supports the concept of varied gene expression in cancerous condition.

Comprehensive transcriptome profiling of squamous cell carcinoma of horn in Bos indicus.

Squamous cell carcinoma (SCC) of horn is frequently observed in Bos indicus affecting 1% of cattle population and accounting 83.34% of total tumours found. The transcriptome profile of horn cancer (HC) tissue and the matched normal (HN) tissue were analysed by RNA-seq using Roche 454 sequencing. A total of 1u2009504u2009900 reads comprising of 612 MB data were used to identify differentially expressed genes using CLC Genomic Workbench. These include up-regulation of KRT6A, KRT6B, KRT6C, KRT14, SFN, KRT84, PI3, COL17A1, ANLN, SERPINB5 and down-regulation of BOLA, SCGB1A1, CXCL17, KRT19, BPIFB1, NR4A1 and TFF3 in HC, which are involved in regulation of gene transcription, cell proliferation, apoptosis, cell survival and metabolic pathways. The qPCR analysis of several targets suggested concordance of gene expression profile with RNA-seq analysis. The present findings would provide basis for further screening of genes and identification of markers for early diagnosis and therapeutic intervention of HC.

Identification of novel splice variants in horn cancer by RNA-Seq analysis in Zebu cattle.

Horn cancer accounts for nearly 83% of total tumors found in Indian Zebu cattle, which results in chronic suffering and causes heavy economic losses. Alternative splicing has been frequently implicated in the various types of cancer progression. Utilizing the transcriptome sequence generated by next generation sequencing, we analyzed the transcript data for the presence of alternative splicing using BLAT program and identified 27 alternatively spliced genes, of which 12 spliced variants appeared to be the novel spliced candidates. Protein prediction of these novel spliced variants revealed that splice variation has caused either truncation of protein, insertion/deletion of stretch of amino acids or formation of unique carboxy terminus. The RT-PCR analysis confirmed the expression of 8 of the 12 novel spliced variants observed by transcriptome sequencing. Additionally, altered splicing/expression of these novel candidates between cancer and normal tissues revealed by qPCR suggests their potential involvement in the development of horn cancer.

In silico analysis of consequences of non-synonymous SNPs of Slc11a2 gene in Indian bovines.

The aim of our study was to analyze the consequences of non-synonymous SNPs in Slc11a2 gene using bioinformatic tools. There is a current need of efficient bioinformatic tools for in-depth analysis of data generated by the next generation sequencing technologies. SNPs are known to play an imperative role in understanding the genetic basis of many genetic diseases. Slc11a2 is one of the major metal transporter families in mammals and plays a critical role in host defenses. In this study, we performed a comprehensive analysis of the impact of all non-synonymous SNPs in this gene using multiple tools like SIFT, PROVEAN, I-Mutant and PANTHER. Among the total 124 SNPs obtained from amplicon sequencing of Slc11a2 gene by Ion Torrent PGM involving 10 individuals of Gir cattle and Murrah buffalo each, we found 22 non-synonymous. Comparing the prediction of these 4 methods, 5 nsSNPs (G369R, Y374C, A377V, Q385H and N492S) were identified as deleterious. In addition, while tested out for polar interactions with other amino acids in the protein, from above 5, Y374C, Q385H and N492S showed a change in interaction pattern and further confirmed by an increase in total energy after energy minimizations in case of mutant protein compared to the native.

Draft Genome Sequence of Pasteurella multocida subsp. multocida B:2 Strain VTCCBAA264 Isolated from Bubalus bubalis in North India

ABSTRACT The Pasteurella multocida subsp. multocida B:2 serotype causes hemorrhagic septicemia in bubalines with high morbidity and mortality in the Indian subcontinent. We report the draft genome sequence of Pasteurella multocida strain VTCCBAA264 isolated from the small-intestine of a buffalo calf that died of high fever.