Microbial diversity and community composition of caecal microbiota in commercial and indigenous Indian chickens determined using 16s rDNA amplicon sequencing.

BACKGROUND: The caecal microbiota plays a key role in chicken health and performance, influencing digestion and absorption of nutrients, and contributing to defence against colonisation by invading pathogens. Measures of productivity and resistance to pathogen colonisation are directly influenced by chicken genotype, but host driven variation in microbiome structure is also likely to exert a considerable indirect influence. METHODS: Here, we define the caecal microbiome of indigenous Indian Aseel and Kadaknath chicken breeds and compare them with the global commercial broiler Cobb400 and Ross 308 lines using 16S rDNA V3-V4 hypervariable amplicon sequencing. RESULTS: Each caecal microbiome was dominated by the genera Bacteroides, unclassified bacteria, unclassified Clostridiales, Clostridium, Alistipes, Faecalibacterium, Eubacterium and Blautia. Geographic location (a measure recognised to include variation in environmental and climatic factors, but also likely to feature varied management practices) and chicken line/breed were both found to exert significant impacts (p < 0.05) on caecal microbiome composition. Linear discriminant analysis effect size (LEfSe) revealed 42 breed-specific biomarkers in the chicken lines reared under controlled conditions at two different locations. CONCLUSION: Chicken breed-specific variation in bacterial occurrence, correlation between genera and clustering of operational taxonomic units indicate scope for quantitative genetic analysis and the possibility of selective breeding of chickens for defined enteric microbiota.

A genome-wide approach to screen for genetic variants in broilers (Gallus gallus) with divergent feed conversion ratio.

Feed conversion ratio (FCR) is an economically important trait in broilers and feed accounts for a significant proportion of the costs involved in broiler production. To explore the contribution of functional variants to FCR trait, we analyzed coding and non-coding single-nucleotide variants (SNVs) across the genome by exome sequencing in seven pairs of full-sibs broilers with divergent FCR and with a sequence coverage at an average depth of fourfold. We identified 192,119 high-quality SNVs, including 30,380 coding SNVs (cSNVs) in the experimental population. We discovered missense SNVs in PGM2, NOX4, TGFBR3, and TMX4, and synonymous SNVs in TSNAX, ITA, HSP90B1, and COL18A1 associated with FCR. Haplotype analyses of genome-wide significant SNVs in PGM2, PHKG1, DGKZ, and SOD2 were also observed with suggestive evidence of haplotype association with FCR. Single-variant and FCR QTL-related genes-based association analyses of SNVs identified newly associated genes for FCR in the regions subjected to targeted exome sequencing. The top seven SNVs were next evaluated in independent replication data sets where SNV chr. 3: 13,990,160 (c. 961G>C) at TMX4 was replicated (p < 0.05). Collectively, we have detected SNVs associated with FCR in broiler as well as identification of SNVs in known FCR QTL region. These findings should facilitate the discovery of causative variants for FCR and contribute to marker-assisted selection.

VDAP-GUI: a user-friendly pipeline for variant discovery and annotation of raw next-generation sequencing data.

Even though next-generation sequencing (NGS) has become an invaluable tool in molecular biology, several laboratories with NGS facilities lack trained Bioinformaticians for data analysis. Here, focusing on the variant detection application of NGS analysis, we have developed a fully automated pipeline, namely Variant Discovery and Annotation Tool-Graphical User Interface (VDAP-GUI), which detects and annotates single nucleotide polymorphisms and insertions/deletions from raw sequence reads. VDAP-GUI consolidates several proven methods in each step such as quality control, trimming, mapping, variant detection and annotation. It supports multiple NGS platforms and has four methodological choices for variant detection. Further, it can re-analyze existing data with alternate thresholds and generates easily interpretable reports in html and tab-delimited formats. Using VDAP-GUI, we have analyzed a publically available human whole-exome sequence dataset. VDAP-GUI is developed using Perl/Tk programming, and is available for free download and use at http://sourceforge.net/projects/vdapgui/ .

Exploring genetic polymorphism in innate immune genes in Indian cattle (Bos indicus) and buffalo (Bubalus bubalis) using next generation sequencing technology.

Activation of innate immunity initiates various cascades of reactions that largely contribute to defense against physical, microbial or chemical damage, prompt for damage repair and removal of causative organisms as well as restoration of tissue homeostasis. Genetic polymorphism in innate immune genes plays prominent role in disease resistance capabilities in various breeds of cattle and buffalo. Here we studied single nucleotide variations (SNP/SNV) and haplotype structure in innate immune genes viz CHGA, CHGB, CHGC, NRAMP1, NRAMP2, DEFB1, BNBD4, BNBD5, TAP and LAP in Gir cattle and Murrah buffalo. Targeted sequencing of exonic regions of these genes was performed by Ion Torrent PGM sequencing platform. The sequence reads obtained corresponding to coding regions of these genes were mapped to reference genome of cattle BosTau7 by BWA program using genome analysis tool kit (GATK). Further variant analysis by Unified Genotyper revealed 54 and 224 SNPs in Gir and Murrah respectively and also 32 SNVs was identified. Among these SNPs 43, 36, 11,32,81,21 and 22 variations were in CHGA, CHGB, CHGC, NRAMP1, NRAMP2, DEFB1 and TAP genes respectively. Among these identified 278 SNPs, 24 were found to be reported in the dbSNP database. Variant analysis was followed by structure formation of haplotypes based on multiple SNPs using SAS software revealed a large number of haplotypes. The SNP discovery in innate immune genes in cattle and buffalo breeds of India would advance our understanding of role of these genes in determining the disease resistance/susceptibility in Indian breeds. The identified SNPs and haplotype data would also provide a wealth of sequence information for conservation studies, selective breeding and designing future strategies for identifying disease associations involving samples from distinct populations.

Isolation and characterization of H9N2 influenza virus isolates from poultry respiratory disease outbreak.

The present study reports isolation and characterization of H9N2 virus responsible for disease characterized by symptoms including difficulty in respiration, head swelling, nasal discharge, reduced feed intake, cyanotic comb, reduced egg production and mortality. Virus isolation from allantoic fluid inoculated with tracheal aspirates and whole genome sequencing of two isolates were performed on an Ion-Torrent sequencer. Phylogenetic analysis revealed that the two H9N2 isolates are reassortant viruses showing a G1-like lineage for HA, NA and NP, a Hok/49/98-like lineage for PB1 and PA, PK/UDL-01/05-like lineage for PB2, IL/90658/00-like lineage for NS and an unknown lineage for M gene. Analyses of the HA cleavage site showed a sequence of (333PARSSR↓GL340) indicating that these isolates are of low pathogenicity. Isolate 2 has leucine at amino acid position 226, a substitution which is associated with mammalian adaptation of avian influenza virus. Isolate 1 has the S31N substitution in the M2 gene that has been associated with drug resistance as well as R57Q and C241Y mutations in the NP gene which are associated with human adaptation. The result reported here gives deep insight in to H9N2 viruses circulating in domestic poultry of India and supports the policy of active efforts to control and manage H9N2 infections in Indian poultry.