Genome-wide association analysis identified both RNA-seq and DNA variants associated to paratuberculosis in Canadian Holstein cattle 'in vitro' experimentally infected macrophages

The role of interleukin-10 receptor alpha (IL10Rα) in Mycobacterium avium subsp. paratuberculosis infection of a mammary epithelial cell line

BACKGROUND

Johne's disease is a chronic wasting disease caused by the bacterium Mycobacterium avium subspecies paratuberculosis (MAP). Johne's disease is highly contagious and MAP infection in dairy cattle can eventually lead to death. With no available treatment for Johne's disease, genetic selection and improvements in management practices could help reduce its prevalence. In a previous study, the gene coding interleukin-10 receptor subunit alpha (IL10Rα) was associated with Johne's disease in dairy cattle. Our objective was to determine how IL10Rα affects the pathogenesis of MAP by examining the effect of a live MAP challenge on a mammary epithelial cell line (MAC-T) that had IL10Rα knocked out using CRISPR/cas9. The wild type and the IL10Rα knockout MAC-T cell lines were exposed to live MAP bacteria for 72 h. Thereafter, mRNA was extracted from infected and uninfected cells. Differentially expressed genes were compared between the wild type and the IL10Rα knockout cell lines. Gene ontology was performed based on the differentially expressed genes to determine which biological pathways were involved.

RESULTS

Immune system processes pathways were targeted to determine the effect of IL10Rα on the response to MAP infection. There was a difference in immune response between the wild type and IL10Rα knockout MAC-T cell lines, and less difference in immune response between infected and not infected IL10Rα knockout MAC-T cells, indicating IL10Rα plays an important role in the progression of MAP infection. Additionally, these comparisons allowed us to identify other genes involved in inflammation-mediated chemokine and cytokine signalling, interleukin signalling and toll-like receptor pathways.

CONCLUSIONS

Identifying differentially expressed genes in wild type and ILR10α knockout MAC-T cells infected with live MAP bacteria provided further evidence that IL10Rα contributes to mounting an immune response to MAP infection and allowed us to identify additional potential candidate genes involved in this process. We found there was a complex immune response during MAP infection that is controlled by many genes.

MicroRNAs modulate immunological and inflammatory responses in Holstein cattle naturally infected with Mycobacterium avium subsp. paratuberculosis

MicroRNAs (miRNAs) regulate the post-transcriptional expression of genes by binding to their target mRNAs. In this study, whole miRNA sequencing was used to compare the expression of miRNAs in ileocecal valve (ICV) and peripheral blood (PB) samples of cows with focal or diffuse paratuberculosis (PTB)-associated lesions in gut tissues versus (vs) control cows without lesions. Among the eight miRNAs differentially expressed in the PB samples from cows with diffuse lesions vs controls, three (miR-19a, miR-144, miR32) were also down-regulated in cows with diffuse vs focal lesions. In the ICV samples, we identified a total of 4, 5, and 18 miRNAs differentially expressed in cows with focal lesions vs controls, diffuse lesions vs controls, and diffuse vs focal lesions, respectively. The differential expression of five microRNAs (miR-19a, miR-144, miR-2425-3p, miR-139, miR-101) was confirmed by RT-qPCR. Next, mRNA target prediction was performed for each differentially expressed miRNA. A functional analysis using the predicted gene targets revealed a significant enrichment of the RNA polymerase and MAPK signaling pathways in the comparison of cows with focal vs no lesions and with diffuse vs focal lesions, respectively. The identified miRNAs could be used for the development of novel diagnostic and therapeutical tools for PTB control.

Modeling Paratuberculosis in Laboratory Animals, Cells, or Tissues: A Focus on Their Applications for Pathogenesis, Diagnosis, Vaccines, and Therapy Studies

Paratuberculosis is a chronic granulomatous enteritis caused by Mycobacterium avium subsp. Paratuberculosis that affects a wide variety of domestic and wild animals. It is considered one of the diseases with the highest economic impact on the ruminant industry. Despite many efforts and intensive research, paratuberculosis control still remains controversial, and the existing diagnostic and immunoprophylactic tools have great limitations. Thus, models play a crucial role in understanding the pathogenesis of infection and disease, and in testing novel vaccine candidates. Ruminant animal models can be restricted by several reasons, related to space requirements, the cost of the animals, and the maintenance of the facilities. Therefore, we review the potential and limitations of the different experimental approaches currently used in paratuberculosis research, focusing on laboratory animals and cell-based models. The aim of this review is to offer a vision of the models that have been used, and what has been achieved or discovered with each one, so that the reader can choose the best model to answer their scientific questions and prove their hypotheses. Also, we bring forward new approaches that we consider worth exploring in the near future.

Complete Genome Sequence of Lacticaseibacillus paracasei Strain VHProbi O44, Isolated from Feces from a Healthy Baby

Lacticaseibacillus paracasei strain VHProbi O44 is a Chinese commercial lactic acid bacterium with several probiotic functions. The whole genome contains a chromosome and three plasmids.

Control of Mycobacterium avium subsp. paratuberculosis load within infected bovine monocyte-derived macrophages is associated with host genetics

The genetic loci influencing individual resistance to Mycobacterium avium subsp. paratuberculosis (MAP) infection are still largely unknown. In the current study, we searched for genetic loci associated with resistance to MAP infection by evaluating the performance of monocyte-derived macrophages (MDMs) isolated from the peripheral blood of 75 healthy Holsteins cows and infected ex vivo with MAP. Bacterial load (log colony-forming units, log CFUs) within MDMs was quantified at 2 h and 7 days p. i. using a BACTEC MGIT 960 instrument. In addition, the expression levels of some genes with important roles in the innate immune response including epiregulin (EREG), complement component C3 (C3), galectin-9 (Gal9), and nitric oxide (NO^-) were measured in the supernatant of the infected cells. DNA from peripheral blood samples of the animals included in the study was isolated and genotyped with the EuroG MD bead Chip (44,779 single nucleotide-polymorphisms, SNPs). Linear mixed models were used to calculate the heritability (h² ) estimates for each indicator of MDM performance, MAP load within MDMs and EREG, C3, Gal9, and NO^-expression. After performing a genome-wide association study, the only phenotypes that showed SNPs with a significant association were the bacterial load within MDMs at 2 h (h² = 0. 87) and 7 days (h² = 0.83) p.i. A total of 6 SNPs, 5 candidate genes, and one microRNA on the Bos taurus chromosomes BTA2, BTA17, BTA18, and BTA21 were associated with MAP load at 2 h p.i. Overlap was seen in two SNPs associated with the log CFUs at 2 h and 7 d p.i. The identified SNPs had negative regression coefficients, and were, therefore, associated with a low bacterial load within MDMs. Some of the identified SNPs were located within QTLs previously associated with longevity, reproductive, and udder health traits. Some of the identified candidate genes; Oxysterol Binding Protein Like 6, Cysteine and Serine Rich Nuclear Protein 3, and the Coiled-Coil Domain Containing 92 regulate cellular cholesterol trafficking and efflux, apoptosis, and interferon production, respectively. Taken together, our results define a heritable and distinct immunogenetic profile in MAP-infected macrophages designed to limit bacterial load early after infection.

Variant Calling from RNA-seq Data Using the GATK Joint Genotyping Workflow

The Genome Analysis Toolkit (GATK) developed at the Broad Institute provides state-of-the-art pipelines for germline and somatic variant discovery and genotyping. Unfortunately, the fully validated GATK pipeline for calling variant on RNAseq data is a Per-sample workflow that does not include the recent improvements seen in modern workflows, especially the possibility to perform joint genotyping analysis. Here, we describe how modern GATK commands from distinct workflows can be combined to call variants on RNAseq samples. We provide a detailed tutorial that starts with raw RNAseq reads and ends with filtered variants, of which some were shown to be associated with bovine paratuberculosis.

Whole Genome Methylation Analysis Reveals Role of DNA Methylation in Cow's Ileal and Ileal Lymph Node Responses to Mycobacterium avium subsp. paratuberculosis Infection

Johne's Disease (JD), caused by Mycobacterium avium subsp paratuberculosis (MAP), is an incurable disease of ruminants and other animal species and is characterized by an imbalance of gut immunity. The role of MAP infection on the epigenetic modeling of gut immunity during the progression of JD is still unknown. This study investigated the DNA methylation patterns in ileal (IL) and ileal lymph node (ILLN) tissues from cows diagnosed with persistent subclinical MAP infection over a one to 4 years period. DNA samples from IL and ILLN tissues from cows negative (MAPneg) (n = 3) or positive for MAP infection (MAPinf) (n = 4) were subjected to whole genome bisulfite sequencing. A total of 11,263 and 62,459 differentially methylated cytosines (DMCs), and 1259 and 8086 differentially methylated regions (DMRs) (FDR<0.1) were found between MAPinf and MAPneg IL and ILLN tissues, respectively. The DMRs were found on 394 genes (denoted DMR genes) in the IL and on 1305 genes in the ILLN. DMR genes with hypermethylated promoters/5'UTR [3 (IL) and 88 (ILLN)] or hypomethylated promoters/5'UTR [10 (IL) and 25 (ILLN)] and having multiple functions including response to stimulus/immune response (BLK, BTC, CCL21, AVPR1A, CHRNG, GABRA4, TDGF1), cellular processes (H2AC20, TEX101, GLA, NCKAP5L, RBM27, SLC18A1, H2AC20BARHL2, NLGN3, SUV39H1, GABRA4, PPA1, UBE2D2) and metabolic processes (GSTO2, H2AC20, SUV39H1, PPA1, UBE2D2) are potential DNA methylation candidate genes of MAP infection. The ILLN DMR genes were enriched for more biological process (BP) gene ontology (GO) terms (n = 374), most of which were related to cellular processes (27.6%), biological regulation (16.6%), metabolic processes (15.4%) and response to stimulus/immune response (8.2%) compared to 75 BP GO terms (related to cellular processes, metabolic processes and transport, and system development) enriched for IL DMR genes. ILLN DMR genes were enriched for more pathways (n = 47) including 13 disease pathways compared with 36 enriched pathways, including 7 disease/immune pathways for IL DMR genes. In conclusion, the results show tissue specific responses to MAP infection with more epigenetic changes (DMCs and DMRs) in the ILLN than in the IL tissue, suggesting that the ILLN and immune processes were more responsive to regulation by methylation of DNA relative to IL tissue. Our data is the first to demonstrate a potential role for DNA methylation in the pathogenesis of MAP infection in dairy cattle.

Exploring Biological Impacts of Prenatal Nutrition and Selection for Residual Feed Intake on Beef Cattle Using Omics Technologies: A Review

Approximately 70% of the cost of beef production is impacted by dietary intake. Maximizing production efficiency of beef cattle requires not only genetic selection to maximize feed efficiency (i.e., residual feed intake (RFI)), but also adequate nutrition throughout all stages of growth and development to maximize efficiency of growth and reproductive capacity, even during gestation. RFI as a measure of feed efficiency in cattle has been recently accepted and used in the beef industry, but the effect of selection for RFI upon the dynamics of gestation has not been extensively studied, especially in the context of fluctuating energy supply to the dam and fetus. Nutrient restriction during gestation has been shown to negatively affect postnatal growth and development as well as fertility of beef cattle offspring. This, when combined with the genetic potential for RFI, may significantly affect energy partitioning in the offspring and subsequently important performance traits. In this review, we discuss: 1) the importance of RFI as a measure of feed efficiency and how it can affect other economic traits in beef cattle; 2) the influence of prenatal nutrition on physiological phenotypes in calves; 3) the benefits of investigating the interaction of genetic selection for RFI and prenatal nutrition; 4) how metabolomics, transcriptomics, and epigenomics have been employed to investigate the underlying biology associated with prenatal nutrition, RFI, or their interactions in beef cattle; and 5) how the integration of omics information is adding a level of deeper understanding of the genetic architecture of phenotypic traits in cattle.

Response to Oxidative Burst-Induced Hypoxia Is Associated With Macrophage Inflammatory Profiles as Revealed by Cellular Genome-Wide Association

Background

In mammalian species, hypoxia is a prominent feature of inflammation. The role of hypoxia in regulating macrophage responses via alteration in metabolic pathways is well established. Recently, oxidative burst-induced hypoxia has been shown in murine macrophages after phagocytosis. Despite the available detailed information on the regulation of macrophage function at transcriptomic and epigenomic levels, the association of genetic polymorphism and macrophage function has been less explored. Previously, we have shown that host genetics controls approximately 80% of the variation in an oxidative burst as measured by nitric oxide (NO^-). Further studies revealed two clusters of transcription factors (hypoxia-related and inflammatory-related) are under the genetic control that shapes macrophages’ pro-inflammatory characteristics.

Material and Methods

In the current study, the association between 43,066 autosomal Single Nucleic Polymorphism (SNPs) and the ability of MDMs in production of NO^- in response to E. coli was evaluated in 58 Holstein cows. The positional candidate genes near significant SNPs were selected to perform functional analysis. In addition, the interaction between the positional candidate genes and differentially expressed genes from our previous study was investigated.

Results

Sixty SNPs on 22 chromosomes of the bovine genome were found to be significantly associated with NO^- production of macrophages. The functional genomic analysis showed a significant interaction between positional candidate genes and mitochondria-related differentially expressed genes from the previous study. Further examination showed 7 SNPs located in the vicinity of genes with roles in response to hypoxia, shaping approximately 73% of the observed individual variation in NO^- production by MDM. Regarding the normoxic condition of macrophage culture in this study, it was hypothesized that oxidative burst is responsible for causing hypoxia at the cellular level.

Conclusion

The results suggest that the genetic polymorphism via regulation of response to hypoxia is a candidate step that perhaps shapes macrophage functional characteristics in the pathway of phagocytosis leading to oxidative burst, hypoxia, cellular response to hypoxia and finally the pro-inflammatory responses. Since all cells in one individual carry the same alleles, the effect of genetic predisposition of sensitivity to hypoxia will likely be notable on the clinical outcome to a broad range of host-pathogen interactions.