Genetic association of polymorphisms in bovine TLR2 and TLR4 genes with Mycobacterium avium subspecies paratuberculosis infection in Indian cattle population

Genome-Wide Association Study Reveals Quantitative Trait Loci and Candidate Genes Associated with High Interferon-gamma Production in Holstein Cattle Naturally Infected with Mycobacterium Bovis

Mycobacterium bovis (Mb) is the causative agent of bovine tuberculosis (bTb). Genetic selection aiming to identify less susceptible animals has been proposed as a complementary measure in ongoing programs toward controlling Mb infection. However, individual animal phenotypes for bTb based on interferon-gamma (IFNɣ) and its use in bovine selective breeding programs have not been explored. In the current study, IFNɣ production was measured using a specific IFNɣ ELISA kit in bovine purified protein derivative (bPPD)-stimulated blood samples collected from Holstein cattle. DNA isolated from the peripheral blood samples collected from the animals included in the study was genotyped with the EuroG Medium Density bead Chip, and the genotypes were imputed to whole-genome sequences. A genome-wide association analysis (GWAS) revealed that the IFNɣ in response to bPPD was associated with a specific genetic profile (heritability = 0.23) and allowed the identification of 163 SNPs, 72 quantitative trait loci (QTLs), 197 candidate genes, and 8 microRNAs (miRNAs) associated with this phenotype. No negative correlations between this phenotype and other phenotypes and traits included in the Spanish breeding program were observed. Taken together, our results define a heritable and distinct immunogenetic profile associated with strong production of IFNɣ in response to Mb.

Elucidation of novel SNPs affecting immune response to classical swine fever vaccination in pigs using immunogenomics approach

The host genetic makeup plays a significant role in causing the within-breed variation among individuals after vaccination. The present study was undertaken to elucidate the genetic basis of differential immune response between high and low responder Landlly (Landrace X Ghurrah) piglets vis-à-vis CSF vaccination. For the purpose, E2 antibody response against CSF vaccination was estimated in sampled animals on the day of vaccination and 21-day post-vaccination as a measure of humoral immune response. Double-digestion restriction associated DNA (ddRAD) sequencing was undertaken on 96 randomly chosen Landlly piglets using Illumina HiSeq platform. SNP markers were called using standard methodology. Genome-wide association study (GWAS) was undertaken in PLINK program to identify the informative SNP markers significantly associated with differential immune response. The results revealed significant SNPs associated with E2 antibody response against CSF vaccination. The genome-wide informative SNPs for the humoral immune response against CSF vaccination were located on SSC10, SSC17, SSC9, SSC2, SSC3 and SSC6. The overlapping and flanking genes (500Kb upstream and downstream) of significant SNPs were CYB5R1, PCMTD2, WT1, IL9R, CD101, TMEM64, TLR6, PIGG, ADIPOR1, PRSS37, EIF3M, and DNAJC24. Functional enrichment and annotation analysis were undertaken for these genes in order to gain maximum insights into the association of these genes with immune system functionality in pigs. The genetic makeup was associated with differential immune response against CSF vaccination in Landlly piglets while the identified informative SNPs may be used as suitable markers for determining variation in host immune response against CSF vaccination in pigs.

Characterization of the complete mitochondrial genome and identification of signature sequence of Indian wild pig

Mitochondrial DNA (mtDNA) serves as a valuable molecular marker for constructing matrilineal genealogies and tracing the evolutionary history of animals. This study aimed to characterize the complete mitochondrial genome of the Indian wild pig (IWB) (Sus scrofa cristatus) and identify IWB-specific DNA sequences that could be used as genomic signatures to differentiate IWB from domestic Indian pigs (IDP) in forensic cases. For the purpose, three wild IWB from a rescue centre were used for the characterization of the mitochondrial genome of the IWB. The mitochondrial genome was sequenced by the primer walking technique using 30 overlapping primers. The mitochondrial genome of the IWB was found to be 16,689 bp long containing 37 genes coding for 2 rRNAs, 22 tRNAs, 13 protein coding genes, and 1 D-loop region similar to the mitogenome of other pigs. Sequence analysis of the D-loop of IWB with other IDP indicated some signature sequence for IWB like duplication and transition event from 1090th to 1099th position, deletion of a 10 bp sequence at the 755th position, insertion of (CA) at the 137th position, and substitution of AT to GA at the 638th position. These variations specially the duplication along with transition event causes creation of unique signature sequence (-ACACAAACCT-) in the IWB that could serve as signature sequences for the IWB and be used as markers for differentiation of IWB from IDP breeds in academic as well as forensic or vetero-legal cases. Overall, a total of 36 polymorphic positions were identified in the IWB, with 29 sites being unique to the IWB only and seven being common to the Doom and HDK75 pig breeds. None of the common polymorphic sites were identified in prevailing domestic pig populations. Phylogenetic analysis of the mitochondrial genome revealed the distinct separation of the IWB from IDP. The results of genetic distance evaluation showed that the Doom pig breed was the closest to the IWB. This study provides valuable insights into the mitogenome characterisation, signature sequence and genetic distance analysis of the IWB and establishes a foundation for future studies on the conservation of this protected species.

Comparative gene expression profile in circulating PBMCs of Bos indicus and crossbred cattle to understand disease tolerance mechanism

The present investigation was performed to compare the global gene expression profile in peripheral blood mononuclear cells (PBMCs) of Bos indicus and crossbred (Bos taurus × B. indicus) cattle. Previously, several studies revealed the disease tolerance potential of B. indicus cattle but underlying genetic mechanism is still not fully explored. The PBMCs model was used for this investigation as it plays crucial role in the immune system regulation. Transcriptomic analysis revealed total 6767 significantly differentially expressed transcripts (fold change (absolute) >2.0, p < .05). In addition, 4149 transcripts were upregulated, 2618 transcripts were downregulated and fold change (absolute) of differentially expressed transcript varied from -223.32 to 213.63. Functional annotation analysis of differentially expressed genes confirmed their role in various molecular pathways viz. innate immune response, antigen processing and presentation, MHC protein complex, defense response to bacterium, regulation of immune response, positive regulation of JAK-STAT cascade, cytoskeletal protein binding, etc. Protein-protein interaction network analysis provided understanding of inter-relationship of immune genes with differentially expressed genes. In conclusion, this study could provide comprehensive information about the dysregulated genes and biological pathways in PBMCs which might be responsible for disease tolerance in B. indicus cattle.

Effect of selected single nucleotide polymorphisms in SLC11A1, ANKRA2, IFNG and PGLYRP1 genes on host susceptibility to Mycobacterium avium subspecies paratuberculosis infection in Indian cattle

Paratuberculosis (PTB) is a chronic infectious enteritis of ruminants, caused by Mycobacterium avium subspecies paratuberculosis (MAP) that brings huge economic loss to the dairy farmers. The study was conducted to explore the association of selected SNPs in IFNG, SLC11A1, ANKRA2 and PGLYRP1 genes with resistance to PTB disease in Indian cattle population. A case-control resource population was established based on the results of diagnostic tests used for detection of MAP infection status viz. ELISA, Johnin PPD test, faecal microscopy and IS900 blood PCR. The PCR-RFLP method was used for genotyping of SNPs. SNPs rs109453173 in SLC11A1, rs110853455 in IFNG and rs41933863 in ANKRA2 genes were significantly (P<0.05) associated with resistance to MAP infection. For SNP rs109453173, GG genotype and G allele was found to be associated with resistance against MAP infection than CC and CG genotypes and C allele, respectively. For SNP rs110853455, AG genotype was found to be associated with susceptibility to MAP infection than AA and GG genotype. For SNP rs41933863, the AG genotype provided three and six times more resistance against MAP infection than GG and AA genotype. The results of this study are suggestive of SNPs rs109453173, rs110853455 and rs41933863 as potential markers for screening MAP resistant cattle and a breeding programme favouring GG genotype and G allele for rs109453173, AG genotype for rs41933863 and against AG genotype for rs110853455 might confer resistance against MAP infection in Indian cattle. However, investigation of these SNPs in an independent and larger population will warrant the strength of association for resistance against MAP infection in cattle.

CRISPR-Cas9-mediated knockout of TLR4 modulates Mycobacterium avium ssp. paratuberculosis cell lysate-induced inflammation in bovine mammary epithelial cells

Toll-like receptor 4 (TLR4) is a pattern-recognition receptor involved in the recognition of microbial pathogens and host alarmins. Ligation to TLR4 initiates a signaling cascade that leads to inflammation. Polymorphisms in bovine TLR4 have been associated with Mycobacterium avium ssp. paratuberculosis (MAP) susceptibility and resistance, the cause of Johne's disease, and milk somatic cell score, a biomarker of mastitis. Although the contribution of TLR4 to recognition of bacterial lipopolysaccharide (LPS) has been well characterized, its role in MAP recognition is less certain. Clustered regularly interspaced short palindromic repeats-Cas9 mediated gene editing was performed to generate TLR4 knockout (KO) mammary epithelial cells to determine if TLR4 expression is involved in the initiation of the host inflammatory response to MAP cell lysate (5 and 10 µg/mL) and Escherichia coli LPS (5 µg/mL). The absence of TLR4 in KO cells resulted in enhanced expression of key inflammatory genes (TNFA and IL6), anti-inflammatory genes (IL10 and SOCS3), and supernatant cytokine and chemokine levels (TNF-α, IL-6, IL-10, CCL3) in response to the MAP cell lysate (10 µg/mL). However, in response to LPS, the KO cells showed reduced expression of key inflammatory genes (TNFA, IL1A, IL1B, and IL6) and supernatant cytokine levels (TNF-α, IL-6, CCL2, IL-8) as compared with unedited cells. Overall, these results confirm that TLR4 is essential for eliciting inflammation in response to LPS; however, exacerbated gene and protein expression in TLR4 KO cells in response to MAP cell lysate suggests a different mechanism of infection and host response for MAP, at least in terms of how it interacts with TLR4. These novel findings show potential divergent roles for TLR4 in mycobacterial infections, and this may have important consequences for the therapeutic control of inflammation in cattle.

Association of genetic variability in CD209 gene with bovine paratuberculosis disease: a case-control study in the Indian cattle population

This study aimed to explore the association of single nucleotide polymorphisms (SNPs) in CD209 gene with the occurrence of bovine paratuberculosis (PTB) disease caused by Mycobacterium avium subspecies paratuberculosis (MAP) in Indian cattle. A total of 213 animals were preliminarily selected on the basis of physical body condition score, which was then screened by a panel of diagnostic tests viz. Johnin, ELISA, fecal microscopy, and fecal culture, for the establishment of a case-control resource population. A total of four SNPs viz. rs208222804, rs211654540, rs208814257, and rs210748127 in CD209 gene were genotyped by PCR-RFLP. All SNPs, except rs210748127, were polymorphic in our population. Genotypic-phenotypic associations were assessed by the PROCLOGISTIC procedure of SAS 9.3. The SNP rs208814257 yielded three genotypes viz. CC, CG, and GG, which were significantly (p < 0.05) different in case as compared to the control population. The odds of CC and CG in comparison to GG genotype were 1.21 and 0.40, respectively. The CG genotype was significantly higher in control population, indicating that this genotype may provide resistance against PTB in our resource population. Upon validation in an independent, larger test population and following biological characterization, SNP rs208814257 can be incorporated in marker panel for selection of animals with greater resistance to MAP infection.

Systematic Assessment of Mycobacterium avium Subspecies Paratuberculosis Infections from 1911-2019: A Growth Analysis of Association with Human Autoimmune Diseases

Mycobacterium avium subsp. paratuberculosis (MAP) is an understudied pathogen worldwide with continuous implications in human autoimmune diseases (ADs). The awareness of MAP appears to be low in many places and its research is at infant stage in many countries. The lack of worldwide coverage of the MAP research landscape calls for urgent research attention and prioritization. This present study aimed to assess MAP global research productivity with an emphasis on its implications in ADs via bibliometric and growth analytic frameworks from authors, countries, institutions, international, disciplines and collaboration network perspectives. MAP primary articles were retrieved from the Scopus database and the Web of Science from 1911 to 2019 via title-specific algorithm. Analytic results of dataset yielded a total of 3889 articles from 581 journals and 20.65 average citations per documents. The annual growth rate of MAP research for the period was 6.31%. Based on a country’s productivity (articles (%), freq. of publication (%)), the USA (887 (22.81%), 26.72%), and Australia (236 (6.07%), 6.07%) ranked the top 2 countries but Egypt and Germany had the highest average growth rate (AGR, 170%) in the last 3 years. MAP studies are generally limited to Europe, Australia, Asia, South America and few nations in Africa. It had positive growth rate (30%–100%) in relation to type 1 diabetes mellitus and rheumatoid arthritis ADs; food science and technology, immunology, agriculture, pathology, and research and experimental medicine, wildlife, environments, virulence, disease resistance, meat and meat products, osteopontin, waste milk and slurry/sludge digestion subjects; but negative growth (−130% to −30%) in ulcerative colitis and Parkinson’s disease and no growth in multiple sclerosis, sarcoidosis, thyroid disorders, psoriasis, and lupus. The mapping revealed a gross lack of collaboration networking in terms of authorship, (intra- and inter-) nationally and institutionally with a generalized collaboration index of 1.82. In conclusion, inadequate resources-, knowledge- and scientific-networking hampered growth and awareness of MAP research globally. The study recommends further research to strengthen evidence of MAP’s epidemiologic prevalence in ADs and proffer practical solution(s) for drug development and point-of-care diagnostics amongst other extended themes.

Potential of TLR-gene diversity in Czech indigenous cattle for resistance breeding as revealed by hybrid sequencing

A production herd of Czech Simmental cattle (Czech Red Pied, CRP), the conserved subpopulation of this breed, and the ancient local breed Czech Red cattle (CR) were screened for diversity in the antibacterial toll-like receptors (TLRs), which are members of the innate immune system. Polymerase chain reaction (PCR) amplicons of TLR1, TLR2, TLR4, TLR5, and TLR6 from pooled DNA samples were sequenced with PacBio technology, with 3- 5 ×  coverage per gene per animal. To increase the reliability of variant detection, the gDNA pools were sequenced in parallel with the Illumina X-ten platform at low coverage ( 60 × per gene). The diversity in conserved CRP and CR was similar to the diversity in conserved and modern CRP, representing 76.4 % and 70.9 % of its variants, respectively. Sixty-eight (54.4 %) polymorphisms in the five TLR genes were shared by the two breeds, whereas 38 (30.4 %) were specific to the production herd of CRP; 4 (3.2 %) were specific to the broad CRP population; 7 (5.6 %) were present in both conserved populations; 5 (4.0 %) were present solely for the conserved CRP; and 3 (2.4 %) were restricted to CR. Consequently, gene pool erosion related to intensive breeding did not occur in Czech Simmental cattle. Similarly, no considerable consequences were found from known bottlenecks in the history of Czech Red cattle. On the other hand, the distinctness of the conserved populations and their potential for resistance breeding were only moderate. This relationship might be transferable to other non-abundant historical cattle breeds that are conserved as genetic resources. The estimates of polymorphism impact using Variant Effect Predictor and SIFT software tools allowed for the identification of candidate single-nucleotide polymorphisms (SNPs) for association studies related to infection resistance and targeted breeding. Knowledge of TLR-gene diversity present in Czech Simmental populations may aid in the potential transfer of variant characteristics from other breeds.