Divergent selection for natural antibodies in poultry in the presence of a major gene.

BACKGROUND: Natural antibodies (NAb) are antibodies that are present in a healthy individual without requiring previous exposure to an exogenous antigen. Selection for high NAb levels might contribute to improved general disease resistance. Our aim was to analyse the genetic background of NAb based on a divergent selection experiment in poultry, and in particular the effect of a polymorphism in the TLR1A gene. METHODS: The study population consisted of a base population from a commercial pure-bred elite white leghorn layer line and seven generations of birds from a High and Low selection line. Birds were selected for total KLH-binding NAb titer (IgTotal). An enzyme-linked immunosorbent assay was performed to determine NAb titers in blood plasma for IgTotal and the antibody isotypes IgM and IgG. NAb titers were available for 10,878 birds. Genotypes for a polymorphism in TLR1A were determined for chickens in generations 5, 6 and 7. The data were analysed using mixed linear animal models. RESULTS: The heritability estimate for IgM was 0.30 and higher than that for IgG and IgTotal (0.12). Maternal environmental effects explained 2 to 3% of the phenotypic variation in NAb. Selection for IgTotal resulted in a genetic difference between the High and Low line of 2.4 titer points (5.1 genetic standard deviation) in generation 7. For IgM, the selection response was asymmetrical and higher in the Low than the High line. The frequency of the TLR1A C allele was 0.45 in the base population and 0.66 and 0.04 in generation 7 of the High and Low line, respectively. The TLR1A polymorphism had large and significant effects on IgTotal and IgM. Estimated genotypic effects suggest full dominance of the TLR1A C allele. Significant TLR1A by generation interactions were detected for IgM and IgTotal. CONCLUSIONS: The effect of a polymorphism in the TLR1A gene on IgTotal and IgM NAb was confirmed. Furthermore, we provide experimental verification of changes in allele frequencies at a major gene with dominant gene action on a quantitative trait that is subjected to mass selection. TLR1A by generation interactions indicate sensitivity to environmental factors.

Genome-wide association study identifies loci influencing natural antibody titers in milk of Dutch Holstein-Friesian cattle.

Natural antibodies (NAb) are produced without any antigenic stimulation as a part of the innate immune system and provide a first line of defense against pathogens. Hence, they may be a useful trait when estimating an animal's potential immune competence and in selection for disease resistance. The aim of this study was to identify genomic regions associated with different NAb traits in milk and potentially describe candidate genes. Milk samples from 1,695 first-lactation Holstein Friesian cows with titer measurements for keyhole limpet hemocyanin, lipopolysaccharide, lipoteichoic acid, and peptidoglycan-binding total NAb and isotypes IgG1, IgM, and IgA were used. Genome-wide association study analyses were performed using imputed 777K SNP genotypes, accounting for relationships using pedigree information. Functional enrichment analysis was performed on the significantly associated genomic regions to look for candidate genes. For IgM NAb, significant associations (false discovery rate <0.05) were found on Bos taurus autosome (BTA) 17, 18, and 21 with candidate genes related to immunoglobulin structure and early B cell development. For IgG1, associations were found on BTA3, and we confirmed a quantitative trait loci on BTA21 previously reported for IgG NAb in serum. Our results provide new insights into the regulation of milk NAb that will help unravel the complex relationship between milk immunoglobulins and disease resistance in dairy cattle.

Genomic Region Containing Toll-Like Receptor Genes Has a Major Impact on Total IgM Antibodies Including KLH-Binding IgM Natural Antibodies in Chickens.

Natural antibodies (NAb) are antigen binding antibodies present in individuals without a previous exposure to this antigen. Keyhole limpet hemocyanin (KLH)-binding NAb levels were previously associated with survival in chickens. This suggests that selective breeding for KLH-binding NAb may increase survival by means of improved general disease resistance. Genome-wide association studies (GWAS) were performed to identify genes underlying genetic variation in NAb levels. The studied population consisted of 1,628 adolescent layer chickens with observations for titers of KLH-binding NAb of the isotypes IgM, IgA, IgG, the total KLH-binding (IgT) NAb titers, total antibody concentrations of the isotypes IgM, IgA, IgG, and the total antibodies concentration in plasma. GWAS were performed using 57,636 single-nucleotide polymorphisms (SNP). One chromosomal region on chromosome 4 was associated with KLH-binding IgT NAb, and total IgM concentration, and especially with KLH-binding IgM NAb. The region of interest was fine mapped by imputing the region of the study population to whole genome sequence, and subsequently performing an association study using the imputed sequence variants. 16 candidate genes were identified, of which FAM114A1, Toll-like receptor 1 family member B (TLR1B), TLR1A, Krüppel-like factor 3 (KLF3) showed the strongest associations. SNP located in coding regions of the candidate genes were checked for predicted changes in protein functioning. One SNP (at 69,965,939 base pairs) received the maximum impact score from two independent prediction tools, which makes this SNP the most likely causal variant. This SNP is located in TLR1A, which suggests a fundamental role of TLR1A on regulation of IgM levels (i.e., KLH-binding IgM NAb, and total IgM concentration), or B cells biology, or both. This study contributes to increased understanding of (genetic) regulation of KLH-binding NAb levels, and total antibody concentrations.

Across-line SNP association study for (innate) immune and behavioral traits in laying hens.

BACKGROUND: An association study between single nucleotide polymorphism markers (SNP) and (innate and adaptive) immune parameters but also feather condition score on the back, rump and belly of laying hens was performed. The immune parameters measured in blood samples were natural and acquired antibody titers and complement activity. Feather condition score as a measure of feather damage was determined, this parameter is closely related to feather pecking behavior in hens housed in groups.The aim of the study was to detect associations between genetic markers and immune parameters and feather condition score across nine lines of laying hens, focusing on the feather peckers as well as on the victims of feather pecking. METHODS: A novel approach based on across-line analysis and testing of the SNP-by-line interaction was performed. RESULTS: In total 59 significant associations between SNP and immune traits were detected. Previously identified QTL were confirmed and new associations of genes regulating immune function identified. The IL17A gene (chromosome 3) influences natural and acquired antibody titers and activation of classical and alternative complement pathways. The major histocompatibility complex on chromosome 16 showed significant association with natural and acquired antibody titers and classical complement activity. The IL12B and IRF1 genes on chromosome 13 were associated with natural antibody titers.The direct effect of the genotype of an individual on its feather condition and the associative effect of the genotype of the cage mates on the individual's feather condition were analyzed. The direct genetic effect can be described as the susceptibility to be pecked at, and the associative genetic effect as the propensity to perform feather pecking. Eleven significant associations were detected for the direct effect, and 81 for the associative effect. The serotonin receptor 2C (HTR2C) on chromosome 4 was highlighted in both analyses. CONCLUSIONS: Our results confirmed previously identified QTL and identified new associations of genes regulating immune function. The results for feather condition score supports existing evidence of involvement of the serotonergic system in feather pecking in laying hens. Immune regulatory genes were found to be associated to feather condition score, revealing relationships between the immune system and behavior.

Genome wide SNP discovery, analysis and evaluation in mallard (Anas platyrhynchos).

BACKGROUND: Next generation sequencing technologies allow to obtain at low cost the genomic sequence information that currently lacks for most economically and ecologically important organisms. For the mallard duck genomic data is limited. The mallard is, besides a species of large agricultural and societal importance, also the focal species when it comes to long distance dispersal of Avian Influenza. For large scale identification of SNPs we performed Illumina sequencing of wild mallard DNA and compared our data with ongoing genome and EST sequencing of domesticated conspecifics. This is the first study of its kind for waterfowl. RESULTS: More than one billion base pairs of sequence information were generated resulting in a 16× coverage of a reduced representation library of the mallard genome. Sequence reads were aligned to a draft domesticated duck reference genome and allowed for the detection of over 122,000 SNPs within our mallard sequence dataset. In addition, almost 62,000 nucleotide positions on the domesticated duck reference showed a different nucleotide compared to wild mallard. Approximately 20,000 SNPs identified within our data were shared with SNPs identified in the sequenced domestic duck or in EST sequencing projects. The shared SNPs were considered to be highly reliable and were used to benchmark non-shared SNPs for quality. Genotyping of a representative sample of 364 SNPs resulted in a SNP conversion rate of 99.7%. The correlation of the minor allele count and observed minor allele frequency in the SNP discovery pool was 0.72. CONCLUSION: We identified almost 150,000 SNPs in wild mallards that will likely yield good results in genotyping. Of these, ~101,000 SNPs were detected within our wild mallard sequences and ~49,000 were detected between wild and domesticated duck data. In the ~101,000 SNPs we found a subset of ~20,000 SNPs shared between wild mallards and the sequenced domesticated duck suggesting a low genetic divergence. Comparison of quality metrics between the total SNP set (122,000 + 62,000 = 184,000 SNPs) and the validated subset shows similar characteristics for both sets. This indicates that we have detected a large amount (~150,000) of accurately inferred mallard SNPs, which will benefit bird evolutionary studies, ecological studies (e.g. disentangling migratory connectivity) and industrial breeding programs.

Microarray data mining using Bioconductor packages.

BACKGROUND: This paper describes the results of a Gene Ontology (GO) term enrichment analysis of chicken microarray data using the Bioconductor packages. By checking the enriched GO terms in three contrasts, MM8-PM8, MM8-MA8, and MM8-MM24, of the provided microarray data during this workshop, this analysis aimed to investigate the host reactions in chickens occurring shortly after a secondary challenge with either a homologous or heterologous species of Eimeria. The results of GO enrichment analysis using GO terms annotated to chicken genes and GO terms annotated to chicken-human orthologous genes were also compared. Furthermore, a locally adaptive statistical procedure (LAP) was performed to test differentially expressed chromosomal regions, rather than individual genes, in the chicken genome after Eimeria challenge. RESULTS: GO enrichment analysis identified significant (raw p-value < 0.05) GO terms for all three contrasts included in the analysis. Some of the GO terms linked to, generally, primary immune responses or secondary immune responses indicating the GO enrichment analysis is a useful approach to analyze microarray data. The comparisons of GO enrichment results using chicken gene information and chicken-human orthologous gene information showed more refined GO terms related to immune responses when using chicken-human orthologous gene information, this suggests that using chicken-human orthologous gene information has higher power to detect significant GO terms with more refined functionality. Furthermore, three chromosome regions were identified to be significantly up-regulated in contrast MM8-PM8 (q-value < 0.01). CONCLUSION: Overall, this paper describes a practical approach to analyze microarray data in farm animals where the genome information is still incomplete. For farm animals, such as chicken, with currently limited gene annotation, borrowing gene annotation information from orthologous genes in well-annotated species, such as human, will help improve the pathway analysis results substantially. Furthermore, LAP analysis approach is a relatively new and very useful way to be applied in microarray analysis.

Detection of QTL for innate: non-specific antibody levels binding LPS and LTA in two independent populations of laying hens.

In the current study results are presented of an experiment dealing with the Natural antibodies which are measured by level of homotopes LPS and LTA. Two independent populations were examined (F2 population descendant from a cross between chickens divergently selected for either High or Low specific Ab responses to SRBC (HL) and F2 cross descendant from lines expressed different behavior concerning feather pecking (FP)). In total 12 QTL were detected to non-specific antibody titers directed to LTA and LPS and at two ages after applying two statistical models in an F2 HL population. Similarly in an FP cross overall seven QTL were detected. Based on obtained results it might be concluded that different QTL regions are associated with immune responses to homotopes LPS and LTA in poultry.

Confirmation of quantitative trait loci affecting fatness in chickens.

In this report we describe the analysis of an advanced intercross line (AIL) to confirm the quantitative trait locus (QTL) regions found for fatness traits in a previous study. QTL analysis was performed on chromosomes 1, 3, 4, 15, 18, and 27. The AIL was created by random intercrossing in each generation from generation 2 (G(2)) onwards until generation 9 (G(9)) was reached. QTL for abdominal fat weight (AFW) and/or percentage abdominal fat (AF%) on chromosomes 1, 3 and 27 were confirmed in the G(9) population. In addition, evidence for QTL for body weight at the age of 5 (BW5) and 7 (BW7) weeks and for the percentage of intramuscular fat (IF%) were found on chromosomes 1, 3, 15, and 27. Significant evidence for QTL was detected on chromosome 1 for BW5 and BW7. Suggestive evidence was found on chromosome 1 for AFW, AF% and IF%, on chromosome 15 for BW5, and on chromosome 27 for AF% and IF%. Furthermore, evidence on the chromosome-wise level was found on chromosome 3 for AFW, AF%, and BW7 and on chromosome 27 for BW5. For chromosomes 4 and 18, test statistics did not exceed the significance threshold.

Estimation of the extent of linkage disequilibrium in seven regions of the porcine genome.

Linkage disequilibrium (LD) refers to the correlation among neighboring alleles, reflecting non-random patterns of association between alleles at (nearby) loci. A better understanding of LD in the porcine genome is of direct relevance for identification of genes and mutations with a certain effect on the traits of interest. Here, 215 SNPs in seven genomic regions were genotyped in individuals of three breeds. Pairwise linkage disequilibrium was calculated for all marker pairs. To estimate the extent of LD, all pairwise LD values were plotted against the distance between the markers. Based on SNP markers in four genomic regions analyzed in three panels from populations of Large White, Dutch Landrace, and Meishan origin, useful LD is estimated to extend for approximately 40 to 60 kb in the porcine genome.

ESTIMATION OF THE EXTENT OF LINKAGE DISEQUILIBRIUM IN SEVEN REGIONS OF THE PORCINE GENOME.

Linkage disequilibrium (LD) refers to the correlation among neighboring alleles, reflecting non-random patterns of association between alleles at (nearby) loci. A better understanding of LD in the porcine genome is of direct relevance for identification of genes and mutations with a certain effect on the traits of interest. Here, 215 SNPs in seven genomic regions were genotyped in individuals of three breeds. Pairwise linkage disequilibrium was calculated for all marker pairs. To estimate the extent of LD, all pairwise LD values were plotted against the distance between the markers. Based on SNP markers in four genomic regions analyzed in three panels from populations of Large White, Dutch Landrace, and Meishan origin, useful LD is estimated to extend for approximately 40 to 60 kb in the porcine genome.

A radiation hybrid map of chicken Chromosome 4.

The mapping resolution of the physical map for chicken Chromosome 4 (GGA4) was improved by a combination of radiation hybrid (RH) mapping and bacterial artificial chromosome (BAC) mapping. The ChickRH6 hybrid panel was used to construct an RH map of GGA4. Eleven microsatellites known to be located on GGA4 were included as anchors to the genetic linkage map for this chromosome. Based on the known conserved synteny between GGA4 and human Chromosomes 4 and X, sequences were identified for the orthologous chicken genes from these human chromosomes by BLAST analysis. These sequences were subsequently used for the development of STS markers to be typed on the RH panel. Using a logarithm of the odds (LOD) threshold of 5.0, nine linkage groups could be constructed which were aligned with the genetic linkage map of this chromosome. The resulting RH map consisted of the 11 microsatellite markers and 50 genes. To further increase the number of genes on the map and to provide additional anchor points for the physical BAC map of this chromosome, BAC clones were identified for 22 microsatellites and 99 genes. The combined RH and BAC mapping approach resulted in the mapping of 61 genes on GGA4 increasing the resolution of the chicken-human comparative map for this chromosome. This enhanced comparative mapping resolution enabled the identification of multiple rearrangements between GGA4 and human Chromosomes 4q and Xp.

Comparative map between chicken chromosome 15 and human chromosomal region 12q24 and 22q11-q12.

The physical and comparative map of GGA15 was improved by the construction of 9 BAC contigs around loci previously mapped on GGA15 by linkage analysis. In total, 240 BAC clones were isolated, covering 30-35% of GGA15, and 120 STS were developed (104 STS derived from BAC end sequences and 18 STS derived within genes). Seventeen chicken orthologues of human genes located on human Chr 22q11-q12 were directly mapped within BAC contigs of GGA15. Furthermore, the partial sequences of the chicken BAC clones were compared with sequences present in the EMBL/GenBank databases and revealed matches to 26 genes, ESTs, and genomic clones located on HSA22q11-q12 and HSA12q24. These results provide a better alignment of GGA15 with the corresponding regions in human and mouse, and improve our knowledge of the evolution and dynamics of the vertebrate genome.

Early in vivo cytokine genes expression in chickens after challenge with Salmonella typhimurium lipopolysaccharide and modulation by dietary n--3 polyunsaturated fatty acids.

We studied the effects of Salmonella typhimurium lipopolysaccharide (LPS) on in vivo cytokine mRNA levels in chickens, and investigated whether these levels could be altered by different nutrients. Two hundred and forty chicks were assigned in a 2 x 4 factorial design of treatments. Factors were intravenous injection with S. typhimurium LPS, or saline (control), and four dietary fat sources: corn oil (CO), linseed oil (LO), menhaden oil and beef tallow (BT). Two hours after injection birds were killed and their spleens removed for RNA extraction. Quantitative real-time RT-PCR assays for mRNA of chicken IFN-gamma, IL-6, IL-8, IL-15, IL-18 and 28S rRNA were used to obtain the in vivo splenic cytokine profiles. Expression levels of IL-6, IL-8, IL-18 and IFN-gamma mRNA increased, but IL-15 mRNA decreased 2h after challenge with LPS compared with saline controls. In saline-injected control chickens, the dietary oil source did not affect the splenic mRNA level of any cytokine. In LPS challenged chickens IFN-gamma mRNA was significantly higher in the chickens fed the fish oil enriched diet compared with the LO, CO and BT enriched diets. The present data imply that avian IL-15 has, at least partially, a different function compared to its mammalian counterpart, and in addition, chicken innate immune responses might be affected differently by n-3 PUFA compared to mammals.