RNA-seq transcriptome profiling of porcine lung from two pig breeds in response to Mycoplasma hyopneumoniae infection

Revelation of genetic diversity and genomic footprints of adaptation in Indian pig breeds

In the present study, the genetic diversity measures among four Indian domestic breeds of pig namely Agonda Goan, Ghurrah, Ghungroo, and Nicobari, of different agro-climatic regions of country were explored and compared with European commercial breeds, European wild boar and Chinese domestic breeds. The double digest restriction site-associated DNA sequencing (ddRADseq) data of Indian pigs (102) and Landrace (10 animals) were generated and whole genome sequencing data of exotic pigs (60 animals) from public data repository were used in the study. The principal component analysis (PCA), admixture analysis and phylogenetic analysis revealed that Indian breeds were closer in ancestry to Chinese breeds than European breeds. European breeds exhibited highest genetic diversity measures among all the considered breeds. Among Indian breeds, Agonda Goan and Ghurrah were found to be more genetically diverse than Nicobari and Ghungroo. The selection signature regions in Indian pigs were explored using iHS and XP-EHH, and during iHS analysis, it was observed that genes related to growth, reproduction, health, meat quality, sensory perception and behavior were found to be under selection pressure in Indian pig breeds. Strong selection signatures were recorded in 24.25-25.25 Mb region of SSC18, 123.25-124 Mb region of SSC15 and 118.75-119.5 Mb region of SSC2 in most of the Indian breeds upon pairwise comparison with European commercial breeds using XP-EHH. These regions were harboring some important genes such as EPHA4 for thermotolerance, TAS2R16, FEZF1, CADPS2 and PTPRZ1 for adaptability to scavenging system of rearing, TRIM36 and PGGT1B for disease resistance and CCDC112, PIAS1, FEM1B and ITGA11 for reproduction.

The Relationship between Complement Components C1R and C5 Gene Polymorphism and the Values of Blood Indices in Suckling Piglets

The main mechanism of innate immunity is the complement system. Its components include the protein products of the C1R and C5 genes, which are involved in the classical activation pathway as well as the inflammatory and cytolytic immune responses, respectively. The aim of this study was to determine the relationship between PCR-restriction fragment length polymorphism in C1R (726T > C) and C5 (1044A > C) genes, and the values of hematological and biochemical blood indices in suckling crossbred (Polish Large White × Polish Landrace × Duroc × Pietrain) piglets (n = 473), considering their age (younger, 21 ± 3 days, n = 274; older, 35 ± 3 days, n = 199) and health status. The frequencies of the C5 genotypes deviated from the Hardy-Weinberg expectations. Younger piglets, healthy piglets, piglets that deviated from physiological norms and older piglets with the C1R TT genotype all had lower white and red blood cell indices. In piglets with the C5 CC genotype, younger piglets, piglets that deviated from physiological norms and older piglets, a greater number and/or percentage of monocytes were recorded in the blood. Older piglets also showed an increase in the number of leukocytes and granulocytes, along with a tendency for a decrease in the percentage of lymphocytes in their blood. We concluded that a polymorphism in the C1R gene may exhibit a functional association or genetic linkage with other genes involved in the process of erythropoiesis. Furthermore the relationship between the C5 gene polymorphism and the number and/or percentage of monocytes in the blood may modify the body's defense abilities. Piglets with the CC genotype, having an increased number/proportion of these cells in their blood, probably display a weakened immune response to pathogens or a chronic stimulation of the immune system.

Regulatory functional role of NLRP3 inflammasome during Mycoplasma hyopneumoniae infection in swine

Mycoplasma hyopneumoniae causes enzootic pneumonia, a highly contagious respiratory disease in swine that causes significant economic losses worldwide. It is unknown whether the nucleotide oligomerization domain-like receptor (NLR) family pyrin domain containing 3 (NLRP3) inflammasome regulates the immune response in swine during M. hyopneumoniae infection. The current study utilized an in vivo swine model of M. hyopneumoniae infection to investigate the regulatory functional role of the NLRP3 inflammasome during M. hyopneumoniae infection. Notable histopathological alterations were observed in M. hyopneumoniae-infected swine tissues, which were associated with an inflammatory response and disease progression. Swine M. hyopneumoniae infection was associated with an increase in the expression of the NLRP3 inflammasome, which stimulated pro-inflammatory cytokines such as tumor necrosis factor-alpha, interleukin 18, and interleukin 1 beta (IL-1β). The impact of the NLRP3 inhibitor, MCC950 on NLRP3 and pro-inflammatory cytokines in M. hyopneumoniae-infected swine was examined to investigate the relationship between the NLRP3 inflammasome and M. hyopneumoniae infection. Taken together, our findings provide strong evidence that the NLRP3 inflammasome plays a critical regulatory functional role in M. hyopneumoniae infection in swine.

Integrated Proteotranscriptomics Reveals Differences in Molecular Immunity between Min and Large White Pig Breeds

Long-term selection or evolution is an important factor governing the development of disease resistance in pigs. To better clarify the molecular mechanisms underlying different levels of disease resistance, we used transcriptomics and proteomics analysis to characterize differences in the immunities between six resistant (Min pig) and six susceptible (Large White, LW) pigs which were raised in the same environment. A total of 135 proteins and 791 genes were identified as being differentially expressed between the Large White and Min pig groups. Protein expression clustering and functional analysis revealed that proteins related to immune system process, humoral immune response, the B cell receptor signaling pathway, lymphocyte-mediated immunity, and innate immune responses were more highly expressed in Min pigs. Transcriptome gene set enrichment analysis was used to reveal that pathways of cell adhesion molecules and antigen processing and presentation are significantly enriched in Min pigs. Integrated proteomics and transcriptomics data analysis identified 16 genes that are differentially expressed at both the mRNA and protein levels. In addition, 13 out of these 16 genes were related to the quantitative trait loci of immune diseases, including neural EGFL-like 2 (NELL2) and lactate dehydrogenase B (LDHB), which are involved in innate immunity. Correlation analysis between the genes/proteins and cytokines shows upregulated proteins in LW pigs in association with immunosuppressive/pro-inflammatory cytokines, such as interleukin (IL) 10, IL6, and tumor necrosis factor alpha. This was further validated using parallel reaction monitoring analysis. In summary, we discovered several potential candidate pathways and key genes/proteins involved in determining differences in disease resistance between the two studied pig breeds, which could provide new insights into the breeding of pigs for disease resistance.

Post-transcriptional regulation through alternative splicing in the lungs of Tibetan pigs under hypoxia

In multicellular organisms, alternative splicing (AS) is central to the regulation of multiple biological processes. To further elucidate the adaptive strategy of AS in the lungs of Tibetan pigs in response to hypoxia, we identified and analyzed five basic AS types and 59,930 AS events in 18,179 genes. We found that approximately 65.10% of the total expressed genes underwent AS in the lungs of Tibetan pigs at a high altitude (TH). The frequencies of AS events were similar among the different groups (5.06-5.30 events in each gene on average). Skipped exons (SEs) were the predominant type of AS event, followed by mutually exclusive exons (MXEs), alternative 3' splice sites (A3SSs) and alternative 5' splice sites (A5SSs). Retained introns (RIs), the remaining type of AS event, showed lower frequencies. Further comparison analysis of differentially expressed genes (DEGs) and differentially spliced genes (DSGs) identified 2,209 differential splicing events in the above 18,000 expressed genes, including 918 increased and 1,291 decreased splicing events between the TH and Tibetan pigs at a low altitude (TL) groups. We identified 227 hypoxia-related genes involved in lung development that were differentially regulated through AS. Gene Ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis clearly identified many DEGs and DSGs at high or low altitude. Seven pathways in the top 20 enriched KEGG terms overlapped for the DEGs and DSGs, including the chemokine signaling pathway, B cell receptor signaling pathway, and cytokine-cytokine receptor interaction, which exert many immunoregulatory and inflammatory actions critical to the lung under hypoxia. Twelve pathways overlapped in hypoxic DEGs and DSGs and included antigen processing, presentation and biosynthesis. ​GO analysis of the DEGs and DSGs among the four groups showed that numerous GO terms were enriched in the biological category, and the proportion of genes with downregulated expression was greater among 227 hypoxic genes than that of all genes. The results suggest that AS plays an essential role in the regulation of gene expression during hypoxia and that numerous genes involved in lung development are differentially regulated through AS.

[Opportunities and risks of the use of genetic resistances to infectious diseases in pigs - an overview]

Demands for health, performance and welfare in pigs, as well as the desire for consumer protection and reduced antibiotic use, require optimal measures in advance of disease development. This includes, in principle, the use of genetically more resistant lines and breeding animals, whose existence has been proven for a wide range of pathogen-host interactions. In addition, attempts are being made to identify the gene variants responsible for disease resistance in order to force the selection of suitable populations, also using modern biotechnical technics. The present work is intended to provide an overview of the research status achieved in this context and to highlight opportunities and risks for the future.The evaluation of the international literature shows that genetic disease resistance exist in many areas of swine diseases. However, polygenic inheritance, lack of animal models and the influence of environmental factors during evaluation render their implementation in practical breeding programs demanding. This is where modern molecular genetic methods, such as Gene Editing, come into play. Both approaches possess their pros and cons, which are discussed in this paper. The most important infectious diseases in pigs, including general diseases and epizootics, diseases of the respiratory and digestive tract and diseases of the immune system are taken into account.

Blood transcriptome analysis revealed the immune changes and immunological adaptation of wildness training giant pandas

The giant panda (Ailuropoda melanoleuca) is a global flagship species for biodiversity conservation. As the time for captive giant pandas to be released into the wild matures, wildness training is provided to allow adaptation to their natural environment. It is assumed that changes in the immune system would be integral in this adaptation from captive to wild, where many more pathogens would be encountered in their natural habitats. Therefore, this study aims to determine the expression changes of immune-related genes and their potential as immunoassay markers for adaptation monitoring in wildness training giant pandas, and then to understand the adaptation strategy of wildness training giant pandas to the wild environment, thereby improving the success rate of panda reintroduction. We obtained 300 differentially expressed genes (DEGs) by RNA-seq, with 239 up-regulated and 61 down-regulated DEGs in wildness training giant pandas compared to captive pandas. Functional enrichment analysis indicated that up-regulated DEGs were enriched in several immune-related terms and pathways. There were 21 immune-related DEGs, in which most of them were up-regulated in wildness training giant pandas, including several critical innate and cellular immune genes. IL1R2 was the most significantly up-regulated gene and is a signature of homeostasis within the immune system. In the protein-protein interaction (PPI) analysis, CXCL8, CXCL10, and CCL5 were identified as the hub immune genes. Our results suggested that wildness training giant pandas have stronger innate and cellular immunity than captive giant pandas, and we proposed that a gene set of CXCL8, CXCL10, CCL5, CD3D, NFKBIA, TBX21, IL12RB2, and IL1R2 may serve as potential immunoassay markers to monitor and assess the immune status of wildness training giant pandas. Our study offers the first insight into immune alterations of wildness training giant pandas, paving the way for monitoring and evaluating the immune status of giant pandas when reintroducing them into the wild.

Comparative Transcriptome Analysis Reveals the Molecular Immunopathogenesis of Chinese Soft-Shelled Turtle (Trionyx sinensis) Infected with Aeromonas hydrophila

Simple Summary

The Chinese soft-shelled turtle (Trionyx sinensis) is an important cultured reptile in East Asia. Hemorrhagic sepsis caused by Aeromonas hydrophila infection is the dominant disease in the aquaculture of Chinese soft-shelled turtles, while the molecular pathology is far from clear due to the lag of research on turtle immunology. It has been reported in mammals and fish that the dysfunction of immune responses to pathogen infections causes host tissue hemorrhagic sepsis. In this study, two groups of turtles with different susceptibility to A. hydrophila infection are found. A comparative transcriptome strategy is adopted to examine the gene expression profiles in liver and spleen for these two phenotypes of turtles post A. hydrophila infection, for the first time revealing the full picture of immune mechanisms against A. hydrophila, which provides new insight into the molecular pathology during A. hydrophila infection in T. sinensis. The findings will promote further investigations on pathogenic mechanisms of hemorrhagic sepsis caused by A. hydrophila infection in T. sinensis, and also will benefit their culture industry.

Abstract

Although hemorrhagic sepsis caused by Aeromonas hydrophila infection is the dominant disease in the aquaculture of Chinese soft-shelled turtle, information on its molecular pathology is seriously limited. In this study, ninety turtles intraperitoneally injected with A. hydrophila exhibited two different phenotypes based on the pathological symptoms, referred to as active and inactive turtles. Comparative transcriptomes of liver and spleen from these two groups at 6, 24, and 72 h post-injection (hpi) were further analyzed. The results showed that cytokine–cytokine receptor interaction, PRRs mediated signaling pathway, apoptosis, and phagocytosis enriched in active and inactive turtles were significantly different. Pro-inflammatory cytokines, the TLR signaling pathway, NLR signaling pathway, and RLR signaling pathway mediating cytokine expression, and apoptosis-related genes, were significantly up-regulated in inactive turtles at the early stage (6 hpi). The significant up-regulation of phagocytosis-related genes occurred at 24 hpi in inactive turtles and relatively lagged behind those in active turtles. The anti-inflammatory cytokine, IL10, was significantly up-regulated during the tested periods (6, 24, and 72 hpi) in active turtles. These findings offer valuable information for the understanding of molecular immunopathogenesis after A. hydrophila infection, and facilitate further investigations on strategies against hemorrhagic sepsis in Chinese soft-shelled turtle T. sinensis.

The Expression Regulatory Network in the Lung Tissue of Tibetan Pigs Provides Insight Into Hypoxia-Sensitive Pathways in High-Altitude Hypoxia

To adapt to a low-oxygen environment, Tibetan pigs have developed a series of unique characteristics and can transport oxygen more effectively; however, the regulation of the associated processes in high-altitude animals remains elusive. We performed mRNA-seq and miRNA-seq, and we constructed coexpression regulatory networks of the lung tissues of Tibetan and Landrace pigs. HBB, AGT, COL1A2, and EPHX1 were identified as major regulators of hypoxia-induced genes that regulate blood pressure and circulation, and they were enriched in pathways related to signal transduction and angiogenesis, such as HIF-1, PI3K-Akt, mTOR, and AMPK. HBB may promote the combination of hemoglobin and oxygen as well as angiogenesis for high-altitude adaptation in Tibetan pigs. The expression of MMP2 showed a similar tendency of alveolar septum thickness among the four groups. These results indicated that MMP2 activity may lead to widening of the alveolar wall and septum, alveolar structure damage, and collapse of alveolar space with remarkable fibrosis. These findings provide a perspective on hypoxia-adaptive genes in the lungs in addition to insights into potential candidate genes in Tibetan pigs for further research in the field of high-altitude adaptation.

Pathogenicity & virulence of Mycoplasma hyopneumoniae

Mycoplasma hyopneumoniae: is the etiological agent of porcine enzootic pneumonia (EP), a disease that impacts the swine industry worldwide. Pathogen-induced damage, as well as the elicited host-response, contribute to disease. Here, we provide an overview of EP epidemiology, control and prevention, and a more in-depth review of M. hyopneumoniae pathogenicity determinants, highlighting some molecular mechanisms of pathogen-host interactions relevant for pathogenesis. Based on recent functional, immunological, and comparative “omics” results, we discuss the roles of many known or putative M. hyopneumoniae virulence factors, along with host molecules involved in EP. Moreover, the known molecular bases of pathogenicity mechanisms, including M. hyopneumoniae adhesion to host respiratory epithelium, protein secretion, cell damage, host microbicidal response and its modulation, and maintenance of M. hyopneumoniae homeostasis during infection are described. Recent findings regarding M. hyopneumoniae pathogenicity determinants also contribute to the development of novel diagnostic tests, vaccines, and treatments for EP.

Quantitative analysis of the blood transcriptome of young healthy pigs and its relationship with subsequent disease resilience

Background

Disease resilience, which is the ability of an animal to maintain performance under disease, is important for pigs in commercial herds, where they are exposed to various pathogens. Our objective was to investigate population-level gene expression profiles in the blood of 912 healthy F1 barrows at ~ 27 days of age for associations with performance and health before and after their exposure to a natural polymicrobial disease challenge at ~ 43 days of age.

Results

Most significant (q < 0.20) associations of the level of expression of individual genes in blood of young healthy pigs were identified for concurrent growth rate and subjective health scores prior to the challenge, and for mortality, a combined mortality-treatment trait, and feed conversion rate after the challenge. Gene set enrichment analyses revealed three groups of gene ontology biological process terms that were related to disease resilience: 1) immune and stress response-related terms were enriched among genes whose increased expression was unfavorably associated with both pre- and post-challenge traits, 2) heme-related terms were enriched among genes that had favorable associations with both pre- and post-challenge traits, and 3) terms related to protein localization and viral gene expression were enriched among genes that were associated with reduced performance and health traits after but not before the challenge.

Conclusions

Gene expression profiles in blood from young healthy piglets provide insight into their performance when exposed to disease and other stressors. The expression of genes involved in stress response, heme metabolism, and baseline expression of host genes related to virus propagation were found to be associated with host response to disease.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07912-8.

Sex-Biased lncRNA Signature in Fetal Growth Restriction (FGR)

Impaired fetal growth is one of the most important causes of prematurity, stillbirth and infant mortality. The pathogenesis of idiopathic fetal growth restriction (FGR) is poorly understood but is thought to be multifactorial and comprise a range of genetic causes. This research aimed to investigate non-coding RNAs (lncRNAs) in the placentas of male and female fetuses affected by FGR. RNA-Seq data were analyzed to detect lncRNAs, their potential target genes and circular RNAs (circRNAs); a differential analysis was also performed. The multilevel bioinformatic analysis enabled the detection of 23,137 placental lncRNAs and 4263 of them were classified as novel. In FGR-affected female fetuses’ placentas (ff-FGR), among 19 transcriptionally active regions (TARs), five differentially expressed lncRNAs (DELs) and 12 differentially expressed protein-coding genes (DEGs) were identified. Within 232 differentially expressed TARs identified in male fetuses (mf-FGR), 33 encompassed novel and 176 known lncRNAs, and 52 DEGs were upregulated, while 180 revealed decreased expression. In ff-FGR ACTA2-AS1, lncRNA expression was significantly correlated with five DEGs, and in mf-FGR, 25 TARs were associated with DELs correlated with 157 unique DEGs. Backsplicing circRNA processes were detected in the range of H19 lncRNA, in both ff- and mf-FGR placentas. The performed global lncRNAs characteristics in terms of fetal sex showed dysregulation of DELs, DEGs and circRNAs that may affect fetus growth and pregnancy outcomes. In female placentas, DELs and DEGs were associated mainly with the vasculature, while in male placentas, disturbed expression predominantly affected immune processes.

Transcriptome Analysis Identifies Candidate Genes and Signaling Pathways Associated With Feed Efficiency in Xiayan Chicken

Feed efficiency is an important economic factor in poultry production, and the rate of feed efficiency is generally evaluated using residual feed intake (RFI). The molecular regulatory mechanisms of RFI remain unknown. Therefore, the objective of this study was to identify candidate genes and signaling pathways related to RFI using RNA-sequencing for low RFI (LRFI) and high RFI (HRFI) in the Xiayan chicken, a native chicken of the Guangxi province. Chickens were divided into four groups based on FE and sex: LRFI and HRFI for males and females, respectively. We identified a total of 1,015 and 742 differentially expressed genes associated with RFI in males and females, respectively. The 32 and 7 Gene Ontology (GO) enrichment terms, respectively, identified in males and females chiefly involved carbohydrate, amino acid, and energy metabolism. Additionally, Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis identified 11 and 5 significantly enriched signaling pathways, including those for nutrient metabolism, insulin signaling, and MAPK signaling, respectively. Protein-protein interaction (PPI) network analysis showed that the pathways involving CAT, ACSL1, ECI2, ABCD2, ACOX1, PCK1, HSPA2, and HSP90AA1 may have an effect on feed efficiency, and these genes are mainly involved in the biological processes of fat metabolism and heat stress. Gene set enrichment analysis indicated that the increased expression of genes in LRFI chickens was related to intestinal microvilli structure and function, and to the fat metabolism process in males. In females, the highly expressed set of genes in the LRFI group was primarily associated with nervous system and cell development. Our findings provide further insight into RFI regulation mechanisms in chickens.

Transcriptomic Analysis of Long-Term Protective Immunity Induced by Vaccination With Mycoplasma gallisepticum Strain ts-304

Live attenuated vaccines are commonly used to control Mycoplasma gallisepticum infections in chickens. M. gallisepticum ts-304 is a novel live attenuated vaccine strain that has been shown to be safe and effective. In this study, the transcriptional profiles of genes in the tracheal mucosa in chickens challenged with the M. gallisepticum wild-type strain Ap3AS at 57 weeks after vaccination with ts-304 were explored and compared with the profiles of unvaccinated chickens that had been challenged with strain Ap3AS, unvaccinated and unchallenged chickens, and vaccinated but unchallenged chickens. At two weeks after challenge, pair-wise comparisons of transcription in vaccinated-only, vaccinated-and-challenged and unvaccinated and unchallenged birds detected no differences. However, the challenged-only birds had significant up-regulation in the transcription of genes and enrichment of gene ontologies, pathways and protein classes involved in infiltration and proliferation of inflammatory cells and immune responses mediated through enhanced cytokine and chemokine production and signaling, while those predicted to be involved in formation and motor movement of cilia and formation of the cellular cytoskeleton were significantly down-regulated. The transcriptional changes associated with the inflammatory response were less severe in these mature birds than in the relatively young birds examined in a previous study. The findings of this study demonstrated that vaccination with the attenuated M. gallisepticum strain ts-304 protects against the transcriptional changes associated with the inflammatory response and pathological changes in the tracheal mucosa caused by infection with M. gallisepticum in chickens for at least 57 weeks after vaccination.

Improvement, identification, and target prediction for miRNAs in the porcine genome by using massive, public high-throughput sequencing data

Despite the broad variety of available microRNA (miRNA) research tools and methods, their application to the identification, annotation, and target prediction of miRNAs in nonmodel organisms is still limited. In this study, we collected nearly all public sRNA-seq data to improve the annotation for known miRNAs and identify novel miRNAs that have not been annotated in pigs (Sus scrofa). We newly annotated 210 mature sequences in known miRNAs and found that 43 of the known miRNA precursors were problematic due to redundant/missing annotations or incorrect sequences. We also predicted 811 novel miRNAs with high confidence, which was twice the current number of known miRNAs for pigs in miRBase. In addition, we proposed a correlation-based strategy to predict target genes for miRNAs by using a large amount of sRNA-seq and RNA-seq data. We found that the correlation-based strategy provided additional evidence of expression compared with traditional target prediction methods. The correlation-based strategy also identified the regulatory pairs that were controlled by nonbinding sites with a particular pattern, which provided abundant complementarity for studying the mechanism of miRNAs that regulate gene expression. In summary, our study improved the annotation of known miRNAs, identified a large number of novel miRNAs, and predicted target genes for all pig miRNAs by using massive public data. This large data-based strategy is also applicable for other nonmodel organisms with incomplete annotation information.