Identification of candidate genes and mutations in QTL regions for immune responses in chicken

Genetic estimates and genome-wide association studies of antibody response in Tanzanian dairy cattle

Identifying the genetic determinants of host defence against infectious pathogens is central to enhancing disease resilience and therapeutic efficacy in livestock. Here, we investigated immune response heritability to important infectious diseases affecting smallholder dairy cattle using variance component analysis. We also conducted genome-wide association studies (GWAS) to identify genetic variants that may help understand the underlying biology of these health traits. By assessing 668,911 single-nucleotide polymorphisms (SNPs) genotyped in 2,045 crossbred cattle sampled from six regions of Tanzania, we identified high levels of interregional admixture and European introgression, which may increase infectious disease susceptibility relative to indigenous breeds. Heritability estimates were low to moderate, ranging from 0.03 (SE ± 0.06) to 0.44 (SE ± 0.07), depending on the health trait. GWAS results revealed several loci associated with seropositivity to the viral diseases Rift Valley fever and bovine viral diarrhoea, the protozoan parasites Neospora caninum and Toxoplasma gondii, and the bacterial pathogens Brucella sp, Leptospira hardjo, and Coxiella burnetii. The identified quantitative trait loci mapped to genes involved in immune defence, tumour suppression, neurological processes, and cell exocytosis. We propose that our results provide a basis for future understanding of the cellular pathways contributing to general and taxon-specific infection responses, and for advancing selective breeding and therapeutic target design.

Antigen specificity affects analysis of natural antibodies

Natural antibodies are used to compare immune systems across taxa, to study wildlife disease ecology, and as selection markers in livestock breeding. These immunoglobulins are present prior to immune stimulation. They are described as having low antigen specificity or polyreactive binding and are measured by binding to self-antigens or novel exogenous proteins. Most studies use only one or two antigens to measure natural antibodies and ignore potential effects of antigen specificity in analyses. It remains unclear how different antigen-specific natural antibodies are related or how diversity among natural antibodies may affect analyses of these immunoglobulins. Using genetically distinct lines of chickens as a model system, we tested the hypotheses that (1) antigen-specific natural antibodies are independent of each other and (2) antigen specificity affects the comparison of natural antibodies among animals. We used blood cell agglutination and enzyme-linked immunosorbent assays to measure levels of natural antibodies binding to four antigens: (i) rabbit erythrocytes, (ii) keyhole limpet hemocyanin, (iii) phytohemagglutinin, or (iv) ovalbumin. We observed that levels of antigen specific natural antibodies were not correlated. There were significant differences in levels of natural antibodies among lines of chickens, indicating genetic variation for natural antibody production. However, line distinctions were not consistent among antigen specific natural antibodies. These data show that natural antibodies are a pool of relatively distinct immunoglobulins, and that antigen specificity may affect interpretation of natural antibody function and comparative immunology.

Genetic Variation in Natural and Induced Antibody Responses in Layer Chickens

Selection of livestock for disease resistance is challenging due to the difficulty in obtaining reliable phenotypes. Antibodies are immunological molecules that provide direct and indirect defenses against infection and link the activities of both the innate and adaptive compartments of the immune system. As a result, antibodies have been used as a trait in selection for immune defense. The goal of this study was to identify genomic regions associated with natural and induced antibodies in chickens using low-pass sequencing. Enzyme-linked immunosorbent assays were used to quantify innate (natural) antibodies binding KLH, OVA, and PHA and induced (adaptive) antibodies binding IBD, IBV, NDV, and REO. We collected plasma from four White Leghorn (WL), two White Plymouth Rock (WPR), and two Rhode Island Red (RIR) lines. Samples numbers ranged between 198 and 785 per breed. GWAS was performed within breed on data pre-adjusted for Line-Hatch-Sex effects using GCTA. A threshold of p = 10-6 was used to select genes for downstream annotation and enrichment analysis with SNPEff and Panther. Significant enrichment was found for the defense/immunity protein, immunoglobulin receptor superfamily, and the antimicrobial response protein in RIR; and the immunoglobulin receptor superfamily, defense/immunity protein, and protein modifying enzyme in WL. However, none were present in WPR, but some of the selected SNP were annotated in immune pathways. This study provides new insights regarding the genetics of the antibody response in layer chickens.

SNP prioritization in targeted sequencing data associated with humoral immune responses in chicken

Our study aimed to identify single nucleotide polymorphisms (SNPs) with a significant impact on the innate immunity represented by antibody response against lipopolysaccharide (LPS) and lipoteichoid acid (LTA) and the adaptive immune response represented toward keyhole limpet hemocyanin (KLH) using the SNP prioritization method. Data set consisted of 288 F2 experimental individuals, created by crossing Green-legged Partridgelike and White Leghorn. The analyzed SNPs were located within 24 short genomic regions of GGA1, GGA2, GGA3, GGA4, GGA9, GGA10, GGA14, GGA18, and GGZ, pre-targeted based on literature references and database information. For the specific antibody response toward KLH at d 0 the most highly prioritized SNP for additive and dominance effects were located on GGA2 in the 3’UTR of MYD88. For the response at d 7, the most highly prioritized SNP pointed at the 3’UTR of MYD88, but potential causal additive variants were located within ADIPOQ and one in PROCR. The highest priority for additive and dominance effects in the antibody response toward lipoteichoic acid at d 0 was attributed to the same SNP, located on GGA2 in the 3’UTR region of MYD88. Two SNPs among the top-10 for additive effect were located in the exon of NOCT. SNPs selected for their additive effect on antibody response toward lipopolysaccharide at d 0 marked 3 genes – NOCT, MYD88, and SNX8, while SNPs selected for their dominance effect marked – NOCT, ADIPOQ, and MYD88. The top-10 variants identified in our study were located in different functional parts of the genome. In the context of causality three groups can be distinguished: variants located in exons of protein coding genes (ADIPOQ, NOCT, PROCR, SNX8), variants within exons of non-coding transcripts, and variants located in genes’ UTR regions. Variants from the first group influence protein structure and variants from both latter groups’ exhibit regulatory roles on DNA (UTR) or RNA (lncRNA).

Sensitivity analysis based on the random forest machine learning algorithm identifies candidate genes for regulation of innate and adaptive immune response of chicken

Two categories of immune responses—innate and adaptive immunity—have both polygenic backgrounds and a significant environmental component. The goal of the reported study was to define candidate genes and mutations for the immune traits of interest in chickens using machine learning–based sensitivity analysis for single-nucleotide polymorphisms (SNPs) located in candidate genes defined in quantitative trait loci regions. Here the adaptive immunity is represented by the specific antibody response toward keyhole limpet hemocyanin (KLH), whereas the innate immunity was represented by natural antibodies toward lipopolysaccharide (LPS) and lipoteichoic acid (LTA). The analysis consisted of 3 basic steps: an identification of candidate SNPs via feature selection, an optimisation of the feature set using recursive feature elimination, and finally a gene-level sensitivity analysis for final selection of models. The predictive model based on 5 genes (MAPK8IP3 CRLF3, UNC13D, ILR9, and PRCKB) explains 14.9% of variance for KLH adaptive response. The models obtained for LTA and LPS use more genes and have lower predictive power, explaining respectively 7.8 and 4.5% of total variance. In comparison, the linear models built on genes identified by a standard statistical analysis explain 1.5, 0.5, and 0.3% of variance for KLH, LTA, and LPS response, respectively. The present study shows that machine learning methods applied to systems with a complex interaction network can discover phenotype-genotype associations with much higher sensitivity than traditional statistical models. It adds contribution to evidence suggesting a role of MAPK8IP3 in the adaptive immune response. It also indicates that CRLF3 is involved in this process as well. Both findings need additional verification.

Transcription analysis of the response of the porcine adrenal cortex to a single subclinical dose of lipopolysaccharide from Salmonella Enteritidis

Lipopolysaccharide (LPS) is a bacterial endotoxin which can participate in the induction of inflammatory responses. LPS may also play a significant role in some neurodegenerative, oncological and metabolic disorders. The aim of the current study was to determine the effect of a subclinical low single dose of LPS from Salmonella Enteritidis administrated in vivo on the transcriptome of porcine adrenal cortex cells, especially gene expression levels, long non-coding RNA (lncRNA) profiles, alternative splicing events and RNA editing sites using RNA-seq technology. The subclinical dose of LPS changed the expression of 354 genes, 27 lncRNA loci and other unclassified RNAs. An analysis of alternative splicing events revealed 104 genes with differentially expressed splice junction sites, and the single nucleotide variant calling approach supported the identification of 376 canonical RNA editing candidates and 7249 allele-specific expression variants. The obtained results suggest that the RIG-I-like receptor signaling pathway, may play a more important role than the Toll-like signaling pathway after the administration of a subclinical dose of LPS. Single subclinical dose of LPS can affect the expression profiles of genes coding peptide hormones, steroidogenic enzymes and transcriptional factors, and modulate the endocrine functions of the gland.

Genome-Wide Association Study of H/L Traits in Chicken

Presently, the heterophil-to-lymphocyte (H/L) ratio is being studied extensively as a disease resistance trait. Through intricate mechanisms to identify and destroy pathogenic microorganisms, heterophils play a pivotal role in the immune defense systems of avian species. To reveal the genetic basis and molecular mechanisms affecting the H/L ratio, phenotypic and H/L data from 1650 white feather chicken broilers were used in performing a genome-wide association study. A self-developed, chicken-specific 55K chip was used for heterophils, lymphocytes, and H/L classification, according to individual genomic DNA profiles. We identified five significant single nucleotide polymorphisms (SNPs) when the genome-wide significance threshold was set to 5% (p < 2.42 × 10^-6). A total of 15 SNPs obtained seemingly significant levels (p < 4.84 × 10^-5). Gene annotation indicated that CARD11 (Caspase recruitment domain family member 11), BRIX1 (Biogenesis of ribosomes BRX1), and BANP (BTG3 associated nuclear protein) play a role in H/L-associated cell regulation and potentially constitute candidate gene regions for cellular functions dependent on H/L ratios. These results lay the foundation for revealing the genetic basis of disease resistance and future marker-assisted selection for disease resistance.

Identification of quantitative trait loci affecting production and biochemical traits in a unique Japanese quail resource population

The objective of the current study was to identify QTL associated with body weight, growth rate, egg quality traits, concentration of selected blood plasma, and yolk lipids as well as concentration of selected macro- and microelements, color, pH, basic chemical composition, and drip loss of breast muscle of Japanese quail (Coturnix japonica). Twenty-two meat-type males (line F33) were crossed with twenty-two laying-type females (line S22) to produce a generation of F1 hybrids. The F2 generation was created by mating 44 randomly chosen F1 hybrids, which were full siblings. The birds were individually weighed from the first to eighth week of age. At the age of 19 wk, 2 to 4 eggs were individually collected from each female and an analysis of the egg quality traits was performed. At slaughter, blood and breast muscles were collected from 324 individuals of the resource population. The basic chemical composition, concentration of chosen macro- and microelements, color, pH, and drip loss were determined in the muscle samples. The concentration of chosen lipids was determined in egg yolk and blood plasma. In total, 30 microsatellite markers located on chromosome 1 and 2 were genotyped. QTL mapping including additive and dominance genetic effects revealed 6 loci on chromosome 1 of the Japanese quail affecting the egg number, egg production rate, egg weight, specific gravity, egg shell weight, concentration of Na in breast muscle. In turn, there were 9 loci on chromosome 2 affecting the body weight in the first, fourth, and sixth week of age, growth rate in the second and seventh week of age, specific gravity, concentration of K and Cu in breast muscle, and the levels of triacylglycerols in blood plasma. In this study, QTL with a potential effect on the Na, K, and Cu content in breast muscles in poultry and on specific gravity in the Japanese quail were mapped for the first time.

Circular RNA expression profiles of peripheral blood mononuclear cells in rheumatoid arthritis patients, based on microarray chip technology

Rheumatoid arthritis (RA) is a systemic autoimmune disease characterized by chronic synovial inflammation and finally leads to variable degrees of bone and cartilage erosion. The diagnosis of RA is not an accurate indicator, but a series of scores and the mechanisms underlying it remain only partially understood. The present study explored whether circular RNAs (circRNAs) contribute to the RA pathophysiological mechanism. Total RNA from peripheral blood mononuclear cells of 10 RA patients and 10 healthy controls were extracted and circRNA expression profiling was followed by microarray analysis. In addition, circRNA interactions with microRNAs were performed and microRNA response elements were listed to identify differentially expressed binding site targets in RA. Reverse transcription-quantitative polymerase chain reaction amplification (RT-qPCR) was used to verify the differential expression of circRNAs. A total of 584 circRNAs were differentially expressed in RA patients vs. healthy controls, by circRNA microarray, including 255 circRNAs which were significantly upregulated and 329 downregulated among the RA samples. RT-qPCR validation demonstrated that the expression levels of hsa_circRNA_104194, hsa_circRNA_104593, hsa_circRNA_103334, hsa_circRNA_101407 and hsa_circRNA_102594 were consistent with the results from the microarray analysis. The current study presented differentially expressed circRNAs and their corresponding microRNA binding sites in RA. circRNAs may exhibit a role in the regulation of expression of symbol genes that influence the occurrence and development of RA.

Genome-wide association studies of immune, disease and production traits in indigenous chicken ecotypes

Background

The majority of chickens in sub-Saharan Africa are indigenous ecotypes, well adapted to the local environment and raised in scavenging production systems. Although they are generally resilient to disease challenge, routine vaccination and biosecurity measures are rarely applied and infectious diseases remain a major cause of mortality and reduced productivity. Management and genetic improvement programmes are hampered by lack of routine data recording. Selective breeding based on genomic technologies may provide the means to enhance sustainability. In this study, we investigated the genetic architecture of antibody response to four major infectious diseases [infectious bursal disease (IBDV), Marek’s disease (MDV), fowl typhoid (SG), fowl cholera (PM)] and resistance to Eimeria and cestode parasitism, along with two production traits [body weight and body condition score (BCS)] in two distinct indigenous Ethiopian chicken ecotypes. We conducted variance component analyses, genome-wide association studies, and pathway and selective sweep analyses.

Results

The large majority of birds was found to have antibody titres for all pathogens and were infected with both parasites, suggesting almost universal exposure. We derived significant moderate to high heritabilities for IBDV, MDV and PM antibody titres, cestodes infestation, body weight and BCS. We identified single nucleotide polymorphisms (SNPs) with genome-wide significance for each trait. Based on these associations, we identified for each trait, pathways, networks and functional gene clusters that include plausible candidate genes. Selective sweep analyses revealed a locus on chromosome 18 associated with viral antibody titres and resistance to Eimeria parasitism that is within a positive selection signal. We found no significant genetic correlations between production, immune and disease traits, implying that selection for altered antibody response and/or disease resistance will not affect production.

Conclusions

We confirmed the presence of genetic variability and identified SNPs significantly associated with immune, disease and production traits in indigenous village chickens. Results underpin the feasibility of concomitant genetic improvement for enhanced antibody response, resistance to parasitism and productivity within and across indigenous chicken ecotypes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-016-0252-7) contains supplementary material, which is available to authorized users.