Major histocompatibility complex genes exhibit a potential immunological role in mixed Eimeria-infected broiler cecum analyzed using RNA sequencing

OBJECTIVE

This study was conducted to investigate the differential expression of the major histocompatibility complex (MHC) gene region in Eimeria-infected broiler.

METHODS

We profiled gene expression of Eimeria-infected and uninfected ceca of broilers sampled at 4, 7, and 21 days post-infection (dpi) using RNA sequencing. Differentially expressed genes (DEGs) between two sample groups were identified at each time point. DEGs located on chicken chromosome 16 were used for further analysis. Kyoto encyclopedia of genes and genomes (KEGG) pathway analysis was conducted for the functional annotation of DEGs.

RESULTS

Fourteen significant (false discovery rate <0.1) DEGs were identified at 4 and 7 dpi and categorized into three groups: MHC-Y class I genes, MHC-B region genes, and non-MHC genes. In Eimeria-infected broilers, MHC-Y class I genes were upregulated at 4 dpi but downregulated at 7 dpi. This result implies that MHC-Y class I genes initially activated an immune response, which was then suppressed by Eimeria. Of the MHC-B region genes, the DMB1 gene was upregulated, and TAP-related genes significantly implemented antigen processing for MHC class I at 4 dpi, which was supported by KEGG pathway analysis.

CONCLUSION

This study is the first to investigate MHC gene responses to coccidia infection in chickens using RNA sequencing. MHC-B and MHC-Y genes showed their immune responses in reaction to Eimeria infection. These findings are valuable for understanding chicken MHC gene function.

Changes in expression levels of erythrocyte and immune-related genes are associated with high altitude polycythemia

BACKGROUND

As a chronic mountain sickness(CMS) with the highest incidence and the greatest harm, the pathogenesis of high altitude polycythemia (HAPC) is still not fully understood.

METHODS

37 HAPC patients and 42 healthy subjects were selected from plateau, and peripheral venous blood samples were collected for transcriptome sequencing on Illumina NovaSeq platform. The sequenced data were analyzed by bioinformatics and phenotypic association analysis.

RESULTS

The results showed significant differences in multiple clinical indicators including RBC and HGB et al. existed between HAPC and control. Based on the RNA-seq data, 550 genes with significant differential expression were identified in HAPC patients. GO and KEGG pathway enrichment analysis showed that the up-regulated genes were mainly enriched in processes such as erythrocyte differentiation and development and homeostasis of number of cells, while the down-regulated genes were mainly enriched in categories such as immunoglobulin production, classical pathway of complement activation and other biological processes. The coupling analysis of differential expression genes(DEGs) and pathological phenotypes revealed that 91 DEGs were in close correlation with in the phenotype of red blood cell volume distribution (width-CV and width-SD), and they were all up-regulated in HAPC and involved in the process of erythrocyte metabolism. Combined with the functional annotation of DEGs and literature survey, we found that the expression of several potential genes might be responsible for pathogenesis of HAPC. Besides, cell type deconvolution analysis result suggested that the changes in the number of some immune cell types was significantly lower in HAPC patients than control, implying the autoimmune level of HAPC patients was affected to a certain extent.

CONCLUSION

This study provides an important data source for understanding the pathogenesis and screening pathogenic genes of HAPC. We found for the first time that there was a significant correlation between HAPC and the pathological phenotype of width-CV and width-SD, wherein the enriched genes were all up-regulated expressed and involved in the process of erythrocyte metabolism. Although the role of these genes needs to be further studied, the candidate genes can provide a starting point for functionally pinning down the underlying mechanism of HAPC.

Characterizing Milk Production Related Genes in Holstein Using RNA-seq

Although the chemical, physical, and nutritional properties of bovine milk have been extensively studied, only a few studies have attempted to characterize milk-synthesizing genes using RNA-seq data. RNA-seq data was collected from 21 Holstein samples, along with group information about milk production ability; milk yield; and protein, fat, and solid contents. Meta-analysis was employed in order to generally characterize genes related to milk production. In addition, we attempted to investigate the relationship between milk related traits, parity, and lactation period. We observed that milk fat is highly correlated with lactation period; this result indicates that this effect should be considered in the model in order to accurately detect milk production related genes. By employing our developed model, 271 genes were significantly (false discovery rate [FDR] adjusted p-value<0.1) detected as milk production related differentially expressed genes. Of these genes, five (albumin, nitric oxide synthase 3, RNA-binding region (RNP1, RRM) containing 3, secreted and transmembrane 1, and serine palmitoyltransferase, small subunit B) were technically validated using quantitative real-time polymerase chain reaction (qRT-PCR) in order to check the accuracy of RNA-seq analysis. Finally, 83 gene ontology biological processes including several blood vessel and mammary gland development related terms, were significantly detected using DAVID gene-set enrichment analysis. From these results, we observed that detected milk production related genes are highly enriched in the circulation system process and mammary gland related biological functions. In addition, we observed that detected genes including caveolin 1, mammary serum amyloid A3.2, lingual antimicrobial peptide, cathelicidin 4 (CATHL4), cathelicidin 6 (CATHL6) have been reported in other species as milk production related gene. For this reason, we concluded that our detected 271 genes would be strong candidates for determining milk production.

Molecular Characterization and Expression Analysis of Equine Vascular Endothelial Growth Factor Alpha (VEGFα) Gene in Horse (Equus caballus)

The objective of this study was to determine the molecular characteristics of the horse vascular endothelial growth factor alpha gene (VEGFα) by constructing a phylogenetic tree, and to investigate gene expression profiles in tissues and blood leukocytes after exercise for development of suitable biomarkers. Using published amino acid sequences of other vertebrate species (human, chimpanzee, mouse, rat, cow, pig, chicken and dog), we constructed a phylogenetic tree which showed that equine VEGFα belonged to the same clade of the pig VEGFα. Analysis for synonymous (Ks) and non-synonymous substitution ratios (Ka) revealed that the horse VEGFα underwent positive selection. RNA was extracted from blood samples before and after exercise and different tissue samples of three horses. Expression analyses using reverse transcription-polymerase chain reaction (RT-PCR) and quantitative-polymerase chain reaction (qPCR) showed ubiquitous expression of VEGFα mRNA in skeletal muscle, kidney, thyroid, lung, appendix, colon, spinal cord, and heart tissues. Analysis of differential expression of VEGFα gene in blood leukocytes after exercise indicated a unimodal pattern. These results will be useful in developing biomarkers that can predict the recovery capacity of racing horses.