A high-density SNP panel reveals extensive diversity, frequent recombination and multiple recombination hotspots within the chicken major histocompatibility complex B region between BG2 and CD1A1

Mapping quantitative trait loci regions associated with Marek's disease on chicken autosomes by means of selective DNA pooling

Marek's Disease (MD), which can result in neurological damage and tumour formation, has large effects on the economy and animal welfare of the poultry industry worldwide. Previously, we mapped autosomal MD QTL regions (QTLRs) by individual genotyping of an F6 population from a full-sib advanced intercross line. We further mapped MD QTLRs on the chicken Z chromosome (GGZ) using the same F6 population, and by selective DNA pooling (SDP) of 8 elite egg production lines. Here we used SDP of the same pools used on GGZ to map autosomal MD QTLRs. Thirty-seven QTLRs were found. Seven of the QTLRs were tested by all sires from the same 8 lines, individually genotyped for QTLR markers. Five of the tested QTLRs were confirmed. Linkage disequilibrium (LD) was calculated for all QTLR markers on the same chromosome, and complex LD blocks were found. Distribution of P and LD values were used to assess the QTLR causative elements. Allele substitution effects were calculated based on both pooled SNP microarray genotypes, and individual genotypes of QTLRs markers. Substantial allele effect and contribution to the phenotypic and genotypic variation were obtained. The results explain part of the MD response, and provide targets for mitigating MD.

Different MHC class I cell surface expression levels in diverse chicken lines, associations with B blood group, and proposed relationship to antigen-binding repertoire

The Major Histocompatibility Complex (MHC) is a cluster of genes with primarily immune-related functions. The MHC class I genes are responsible for self- versus non-self-recognition and viral antigen presentation to T lymphocytes. The chicken MHC class I protein binds its cognate antigen(s) over a repertoire spectrum ranging from promiscuous (generalist) to fastidious (specialist). The MHC class I protein expression level at the cell surface is inversely related to the promiscuity of its peptide-binding repertoire. In our study, erythrocytes from 6 diverse and highly inbred lines of chickens, a closed broiler line, and a highly advanced intercross line were evaluated for MHC class I antigen expression level by flow cytometry using monoclonal antibodies to chicken MHC class I molecules. In chickens, the B blood group antigens include the MHC class I antigen expressed from the MHC. Thus, the B blood group has historically been used as a genetic marker for Marek's Disease virus response. Erythrocytes of the inbred lines were blood typed by serology. The B21 blood type is widely recognized as relatively resistant to Marek's disease and regarded as an MHC class I generalist with low MHC class I expression. The Spanish line, which types serologically as B21.1 (similar to B21), was the lowest MHC class I expressing line. The two sublines (B5.1 and B15.2) of the Fayoumi breed, which significantly differed in their MHC class I expression, also differ in response to multiple pathogens. These defined genetic lines of chickens, with distinct MHC class I expression levels, provide an excellent platform to further interrogate the hypothesis of high or low MHC class I expression (antigenic specialists or generalists, respectively) determining diverse responses to pathogens.

Association of Major Histocompatibility Complex Polymorphism With Acute Phase Response in Broiler Chicken

BACKGROUND

Stress associated with changes in host immunity occurs in response to altered environmental conditions, endogenous imbalances, infectious agents and harmful stimuli. The importance of genetic diversity in chickens has increased due to individual immune differences towards resistance and susceptibility to various stimuli.

OBJECTIVES

This study aimed to investigate the association of major histocompatibility complex (MHC) polymorphism with acute phase response (APR) in Ross 308 broiler chickens.

METHODS

The allelic diversity of the LEI0258 microsatellite marker was determined in 120 Ross broilers. In addition, acute phase proteins (APPs), including serum amyloid A (SAA) and alpha-1-acid glycoprotein (AGP), were analysed as markers of the APR. Furthermore, leukocyte count and the heterophil/lymphocyte ratio (H/L ratio) were examined. The antibody response to the Newcastle disease vaccine (NDV) was also measured to assess humoral mediated immunity. Lastly, the correlation between immune responses and MHC alleles was investigated to identify the most effective alleles in a stress-related situation.

RESULTS

A total of six alleles, ranging from 195 to 448 bp, were identified. Association study revealed a significant influence of MHC alleles on APPs in Ross population (p < 0.05). Notably, Allele 448 had a significant correlation with SAA concentration and the H/L ratio. Allele 207 displayed a positive association with AGP concentration, whereas Allele 195 showed a negative association. Furthermore, a significant association was observed between Allele 448 and basopenia, as well as between Allele 195 and monocytosis.

CONCLUSIONS

Results confirmed the significance of MHC as a candidate gene marker for immune responses, which supports its use for vaccine design, genetic improvement of disease-resistant traits and resource conservation in commercial broiler chickens.

Dual-Approach Co-expression Analysis Framework (D-CAF) Enables Identification of Novel Circadian Regulation From Multi-Omic Timeseries Data

The circadian clock is a central driver of many biological and behavioral processes, regulating the levels of many genes and proteins, termed clock controlled genes and proteins (CCGs/CCPs), to impart biological timing at the molecular level. While transcriptomic and proteomic data has been analyzed to find potential CCGs and CCPs, multi-omic modeling of circadian data, which has the potential to enhance the understanding of circadian control of biological timing, remains relatively rare due to several methodological hurdles. To address this gap, a Dual-approach Co-expression Analysis Framework (D-CAF) was created to perform perturbation-robust co-expression analysis on time-series measurements of both transcripts and proteins. Applying this D-CAF framework to previously gathered transcriptomic and proteomic data from mouse macrophages gathered over circadian time, we identified small, highly significant clusters of oscillating transcripts and proteins in the unweighted similarity matrices and larger, less significant clusters of of oscillating transcripts and proteins using the weighted similarity network. Functional enrichment analysis of these clusters identified novel immunological response pathways that appear to be under circadian control. Overall, our findings suggest that D-CAF is a tool that can be used by the circadian community to integrate multi-omic circadian data to improve our understanding of the mechanisms of circadian regulation of molecular processes.

A Dominance Hypothesis Argument for Historical Genetic Gains and the Fixation of Heterosis in Octoploid Strawberry

Heterosis was the catalyst for the domestication of cultivated strawberry (Fragaria × ananassa), an interspecific hybrid species that originated in the 1700s. The hybrid origin was discovered because the phenotypes of spontaneous hybrids transgressed those of their parent species. The transgressions included fruit yield increases and other genetic gains in the twentieth century that sparked the global expansion of strawberry production. The importance of heterosis to the agricultural success of the hybrid species, however, has remained a mystery. Here we show that heterosis has disappeared (become fixed) among improved hybrids within a population (the California population) that has been under long-term selection for increased fruit yield, weight, and firmness. We found that the highest yielding hybrids are among the most highly inbred (59-79%), which seems counterintuitive for a highly heterozygous, outbreeder carrying heavy genetic loads. Although faint remnants of heterosis were discovered, the between-parent allele frequency differences and dispersed favorable dominant alleles necessary for heterosis have decreased nearly genome-wide within the California population. Conversely, heterosis was prevalent and significant among wide hybrids, especially for fruit count, a significant driver of genetic gains for fruit yield. We attributed the disappearance (fixation) of heterosis within the California population to increased homozygosity of favorable dominant alleles and inbreeding associated with selection, random genetic drift, and selective sweeps. Despite historical inbreeding, the highest yielding hybrids reported to-date are estimated to be heterozygous for 20,370-44,280 of 97,000-108,000 genes in the octoploid genome, the equivalent of an entire diploid genome or more.

Identification of new major histocompatibility complex-B Haplotypes in Bangladesh native chickens

OBJECTIVE

The major histocompatibility complex in chicken demonstrates a great range of variations within varities, breeds, populations and that can eventually influence their immuneresponses. The preset study was conducted to understand the major histocompatibility complex-B (MHC-B) variability in five major populations of Bangladesh native chicken: Aseel, Hilly, Junglefowl, Non-descript Deshi, and Naked Neck.

METHODS

These five major populations of Bangladesh native chicken were analyzed with a subset of 89 single nucleotide polymorphisms (SNPs) in the high-density MHC-B SNP panel and Kompetitive Allele-Specific polymerase chain reaction genotyping was applied. To explore haplotype diversity within these populations, the results were analyzed both manually and computationally using PHASE 2.1 program. The phylogenetic investigations were also performed using MrBayes program.

RESULTS

A total of 136 unique haplotypes were identified within these five Bangladesh chicken populations, and only one was shared (between Hilly and Naked Neck). Phylogenetic analysis showed no distinct haplotype clustering among the five populations, although they were shared in distinct clades; notably, the first clade lacked Naked Neck haplotypes.

CONCLUSION

The present study discovered a set of unique MHC-B haplotypes in Bangladesh chickens that could possibly cause varied immune reponses. However, further investigations are required to evaluate their relationships with global chicken populations.

Genome-wide investigation to assess copy number variants in the Italian local chicken population

BACKGROUND

Copy number variants (CNV) hold significant functional and evolutionary importance. Numerous ongoing CNV studies aim to elucidate the etiology of human diseases and gain insights into the population structure of livestock. High-density chips have enabled the detection of CNV with increased resolution, leading to the identification of even small CNV. This study aimed to identify CNV in local Italian chicken breeds and investigate their distribution across the genome.

RESULTS

Copy number variants were mainly distributed across the first six chromosomes and primarily associated with loss type CNV. The majority of CNV in the investigated breeds were of types 0 and 1, and the minimum length of CNV was significantly larger than that reported in previous studies. Interestingly, a high proportion of the length of chromosome 16 was covered by copy number variation regions (CNVR), with the major histocompatibility complex being the likely cause. Among the genes identified within CNVR, only those present in at least five animals across breeds (n = 95) were discussed to reduce the focus on redundant CNV. Some of these genes have been associated to functional traits in chickens. Notably, several CNVR on different chromosomes harbor genes related to muscle development, tissue-specific biological processes, heat stress resistance, and immune response. Quantitative trait loci (QTL) were also analyzed to investigate potential overlapping with the identified CNVR: 54 out of the 95 gene-containing regions overlapped with 428 QTL associated to body weight and size, carcass characteristics, egg production, egg components, fat deposition, and feed intake.

CONCLUSIONS

The genomic phenomena reported in this study that can cause changes in the distribution of CNV within the genome over time and the comparison of these differences in CNVR of the local chicken breeds could help in preserving these genetic resources.

Replacement of Vitamin E by an Extract from an Olive Oil by-Product, Rich in Hydroxytyrosol, in Broiler Diets: Effects on Liver Traits, Oxidation, Lipid Profile, and Transcriptome

The study examines the effect of replacing vitamin E (VE) with a liquid obtained from alpeorujo, an olive oil by-product rich in hydroxytyrosol (HT), as an antioxidant in broiler chicken feeds on the gene expression, lipid profile, and oxidation in the liver. There were five diets that differed only in the substitution of supplemental VE (0 to 40 mg/kg with differences of 10 mg/kg) by HT (30 to 0 mg/kg with differences of 7.5 mg/kg). A linear decrease (p < 0.05) in α-tocopherol concentration in the liver was observed with the replacement of VE by HT. There were no significant changes in triglyceride, cholesterol, or TBARS concentrations. The hepatic transcriptome showed 378 differentially expressed genes between broilers fed HT15 (20 mg/kg VE and 15 mg/kg HT) and HT0 (40 mg/kg VE) diets (p < 0.05 and fold change less or higher than 1.3). Significant changes in cell cycle, cell nucleus activity, neuroactivity, and necroptosis pathways and functions were observed. It is concluded that the olive oil by-product, rich in HT, could be used to spare VE as an antioxidant in broiler diets without affecting liver lipid and TBARS concentrations. The differential gene expression analysis showed a potential role of olive polyphenols in enhancing the chicken immune response.

Discovery of novel MHC-B haplotypes in Chantecler chickens

The major histocompatibility complex (MHC) is a highly polymorphic cluster of genes which contribute to immune response. Located on chromosome 16, the chicken MHC has great influence over disease resistance and susceptibility. Through the use of a high-density SNP panel which encompasses the MHC-B region, haplotypes can be easily identified. This study aims to use an MHC-B SNP panel to evaluate the MHC-B variability in the Chantecler breed. This breed is native to Quebec, Canada, and is a dual-purpose breed known for its strong resistance to extreme cold temperatures. The Chantecler breed faced a near extinction event in the 1970s, which most likely resulted in a genetic bottleneck and loss of diversity. Despite this, SNP haplotype diversity was observed among 4 Chantecler populations. A total of 8 haplotypes were observed. Of these haplotypes, 6 were previously defined in other breeds, and the other 2 were unique to the Chantecler. Within the populations, the number of haplotypes ranged from 4 to 7, with 3 haplotypes, including the novel BSNP-Chant01, being present in all the groups. This study shows existence of reasonable diversity in the MHC-B region of the Chantecler breed and our results further contribute to understanding the variability of this region in chickens.

Genetic architecture of acute hyperthermia resistance in juvenile rainbow trout (Oncorhynchus mykiss) and genetic correlations with production traits

BACKGROUND

Selective breeding is a promising solution to reduce the vulnerability of fish farms to heat waves, which are predicted to increase in intensity and frequency. However, limited information about the genetic architecture of acute hyperthermia resistance in fish is available. Two batches of sibs from a rainbow trout commercial line were produced: the first (N = 1382) was phenotyped for acute hyperthermia resistance at nine months of age and the second (N = 1506) was phenotyped for main production traits (growth, body length, muscle fat content and carcass yield) at 20 months of age. Fish were genotyped on a 57 K single nucleotide polymorphism (SNP) array and their genotypes were imputed to high-density based on the parent's genotypes from a 665 K SNP array.

RESULTS

The heritability estimate of resistance to acute hyperthermia was 0.29 ± 0.05, confirming the potential of selective breeding for this trait. Since genetic correlations of acute hyperthermia resistance with the main production traits near harvest age were all close to zero, selecting for acute hyperthermia resistance should not impact the main production traits, and vice-versa. A genome-wide association study revealed that resistance to acute hyperthermia is a highly polygenic trait, with six quantitative trait loci (QTL) detected, but explaining less than 5% of the genetic variance. Two of these QTL, including the most significant one, may explain differences in acute hyperthermia resistance across INRAE isogenic lines of rainbow trout. Differences in mean acute hyperthermia resistance phenotypes between homozygotes at the most significant SNP was 69% of the phenotypic standard deviation, showing promising potential for marker-assisted selection. We identified 89 candidate genes within the QTL regions, among which the most convincing functional candidates are dnajc7, hsp70b, nkiras2, cdk12, phb, fkbp10, ddx5, cygb1, enpp7, pdhx and acly.

CONCLUSIONS

This study provides valuable insight into the genetic architecture of acute hyperthermia resistance in juvenile rainbow trout. We show that the selection potential for this trait is substantial and selection for this trait should not be too detrimental to improvement of other traits of interest. Identified functional candidate genes provide new knowledge on the physiological mechanisms involved in acute hyperthermia resistance, such as protein chaperoning, oxidative stress response, homeostasis maintenance and cell survival.

Structure and evolution of the squamate major histocompatibility complex as revealed by two Anolis lizard genomes

The major histocompatibility complex (MHC) is an important genomic region for adaptive immunity and has long been studied in ecological and evolutionary contexts, such as disease resistance and mate and kin selection. The MHC has been investigated extensively in mammals and birds but far less so in squamate reptiles, the third major radiation of amniotes. We localized the core MHC genomic region in two squamate species, the green anole (Anolis carolinensis) and brown anole (A. sagrei), and provide the first detailed characterization of the squamate MHC, including the presence and ordering of known MHC genes in these species and comparative assessments of genomic structure and composition in MHC regions. We find that the Anolis MHC, located on chromosome 2 in both species, contains homologs of many previously-identified mammalian MHC genes in a single core MHC region. The repetitive element composition in anole MHC regions was similar to those observed in mammals but had important distinctions, such as higher proportions of DNA transposons. Moreover, longer introns and intergenic regions result in a much larger squamate MHC region (11.7 Mb and 24.6 Mb in the green and brown anole, respectively). Evolutionary analyses of MHC homologs of anoles and other representative amniotes uncovered generally monophyletic relationships between species-specific homologs and a loss of the peptide-binding domain exon 2 in one of two mhc2β gene homologs of each anole species. Signals of diversifying selection in each anole species was evident across codons of mhc1, many of which appear functionally relevant given known structures of this protein from the green anole, chicken, and human. Altogether, our investigation fills a major gap in understanding of amniote MHC diversity and evolution and provides an important foundation for future squamate-specific or vertebrate-wide investigations of the MHC.

Application of next-generation sequencing for the high-resolution typing of MHC-B in Korean native chicken

The major histocompatibility complex-B (MHC-B) region of chicken is crucially important in their immunogenesis and highly diverse among different breeds, lines, and even populations. Because it determines the resistance/susceptibility to numerous infectious diseases, it is important to analyze this genomic region, particularly classical class I and II genes, to determine the variation and diversity that ultimately affect antigen presentation. This study investigated five lines of indigenous Korean native chicken (KNC) and the Ogye breed using next-generation sequencing (NGS) data with Geneious Prime-based assembly and variant calling with the Genome Analysis Toolkit (GATK) best practices pipeline. The consensus sequences of MHC-B (BG1-BF2) were obtained for each chicken line/breed and their variants were analyzed. All of the Korean native chicken lines possessed an excessive number of variants, including an ample amount of high-impact variants that provided useful information regarding modified major histocompatibility complex molecules. The study confirmed that next-generation sequencing techniques can effectively be used to detect MHC variabilities and the KNC lines are highly diverse for the MHC-B region, suggesting a substantial divergence from red junglefowl.

Adaptive Control of the Meiotic Recombination Landscape by DNA Site-dependent Hotspots With Implications for Evolution

Meiosis is an essential component of the sexual life cycle in eukaryotes. The independent assortment of chromosomes in meiosis increases genetic diversity at the level of whole chromosomes and meiotic recombination increases genetic diversity within chromosomes. The resulting variability fuels evolution. Interestingly, global mapping of recombination in diverse taxa revealed dramatic changes in its frequency distribution between closely related species, subspecies, and even isolated populations of the same species. New insight into mechanisms for these evolutionarily rapid changes has come from analyses of environmentally induced plasticity of recombination in fission yeast. Many different DNA sites, and where identified their binding/activator proteins, control the positioning of recombination at hotspots. Each different class of hotspots functions as an independently controlled rheostat that modulates rates of recombination over a broad dynamic range in response to changing conditions. Together, this independent modulation can rapidly and dramatically alter the global frequency distribution of recombination. This process likely contributes substantially to (i.e., can largely explain) evolutionarily rapid, Prdm9-independent changes in the recombination landscape. Moreover, the precise control mechanisms allow cells to dynamically favor or disfavor newly arising combinations of linked alleles in response to changing extracellular and intracellular conditions, which has striking implications for the impacts of meiotic recombination on evolution.

Genome-wide SNP analysis reveals the selection signatures of two indigenous buffalo breeds in Sichuan

Sichuan Province spawned abundant of indigenous buffalo varieties, which probably harbor valuable gene resources beneficial to the genetic improvement of buffalo. However, limited genetic information was publicly available. To better understand their selection signatures between different populations, we performed a restriction site-associated DNA sequencing (RADseq) to explore genome-wide SNPs among two indigenous breeds of Sichuan buffaloes. As a result, a total of 2,110,077 high-quality SNPs were finally obtained. Population genetic analysis indicated a obviously genetic differentiation between two breeds. The detection of selective genes showed that 995 and 910 protein-coding genes underwent positive selection in Yibin buffalo (GYBS) and Dechang buffalo (XCS). Further functional analysis revealed distinctly discrepant selection in two breeds. Candidate genes that positively selected from Yibin buffaloes have mainly occurred in functions closely related to meat quality, complex living environment adaption capability, and disease resistance. While they were significantly enriched in cell proliferation and cell components in Dechang buffalo, indicating the selection pressure primarily derived from the requirement of organism growth and development speed during breed formation. Our dataset constitutes a promising reservoir of genome-wide SNP markers of Sichuan buffaloes and provides potentially traits selected in different local populations. Such comprehensive genetic resources offer an unprecedented opportunity for genetic association analysis of economically important traits and precision breeding programs in buffaloes.

Chasing genetic correlation breakers to stimulate population resilience to climate change

Global climate change introduces new combinations of environmental conditions, which is expected to increase stress on plants. This could affect many traits in multiple ways that are as yet unknown but will likely require the modification of existing genetic relationships among functional traits potentially involved in local adaptation. Theoretical evolutionary studies have determined that it is an advantage to have an excess of recombination events under heterogeneous environmental conditions. Our study, conducted on a population of radiata pine (Pinus radiata D. Don), was able to identify individuals that show high genetic recombination at genomic regions, which potentially include pleiotropic or collocating QTLs responsible for the studied traits, reaching a prediction accuracy of 0.80 in random cross-validation and 0.72 when whole family was removed from the training population and predicted. To identify these highly recombined individuals, a training population was constructed from correlation breakers, created through tandem selection of parents in the previous generation and their consequent mating. Although the correlation breakers showed lower observed heterogeneity possibly due to direct selection in both studied traits, the genomic regions with statistically significant differences in the linkage disequilibrium pattern showed higher level of heretozygosity, which has the effect of decomposing unfavourable genetic correlation. We propose undertaking selection of correlation breakers under current environmental conditions and using genomic predictions to increase the frequency of these ’recombined’ individuals in future plantations, ensuring the resilience of planted forests to changing climates. The increased frequency of such individuals will decrease the strength of the population-level genetic correlations among traits, increasing the opportunity for new trait combinations to be developed in the future.

A Genome-Wide Association Study for Tolerance to Paratuberculosis Identifies Candidate Genes Involved in DNA Packaging, DNA Damage Repair, Innate Immunity, and Pathogen Persistence

Although the genetic susceptibility to diseases has been extensively studied, the genetic loci and the primary molecular and cellular mechanisms that control disease tolerance are still largely unknown. Bovine paratuberculosis (PTB) is an enteritis caused by Mycobacterium avium subsp. paratuberculosis (MAP). PTB affects cattle worldwide and represents a major issue on animal health. In this study, the associations between host genetic and PTB tolerance were investigated using the genotypes from 277 Spanish Holstein cows with two distinct phenotypes: cases) infected animals with positive PCR and bacteriological culture results but without lesions in gut tissues (N= 24), and controls) animals with negative PCR and culture results but with PTB-associated lesions (N= 253). DNA from peripheral blood of the study population was genotyped with the Bovine EuroG MD Bead Chip, and the corresponding genotypes were imputed to whole-genome sequencing (WGS) data. A genome-wide association study was performed using the WGS data and the defined phenotypes in a case-control approach. A total of 142 single nucleotide polymorphisms (SNPs) were associated (false discovery rate ≤ 0.05, P values between 1.5 × 10^-7 and 5.7 × 10^-7) with tolerance (heritability= 0.55). The 40 SNPs with P-values < 5 × 10^-7 defined 9 QTLs and 98 candidate genes located on BTA4, BTA9, BTA16, BTA25, and BTA26. Some of the QTLs identified in this study overlap with QTLs previously associated with PTB, bovine tuberculosis, mastitis, somatic cell score, bovine diarrhea virus persistent infection, tick resistance, and length of productive life. Two candidate genes with important roles in DNA damage response (ERCC4 and RMI2) were identified on BTA25. Functional analysis using the 98 candidate genes revealed a significant enrichment of the DNA packaging process (TNP2/PRMI1/PRM2/PRM3). In addition, the TNF-signaling (bta04668; TRAF5/CREB5/CASP7/CHUK) and the toxoplasmosis (bta05145; TGFβ2/CHUK/CIITA/SOCS1) pathways were significantly enriched. Interestingly, the nuclear Factor NF-κβ Inhibitor Kinase Alpha (CHUK), a key molecule in the regulation of the NF-κB pathway, was enriched in both pathways. Taken together, our results define a distinct immunogenetic profile in the PTB-tolerant animals designed to control bacterial growth, modulate inflammation, limit tissue damage and increase repair, thus reducing the severity of the disease.

Recombination, selection and the evolution of tandem gene arrays

Multigene families—immunity genes or sensory receptors, for instance—are often subject to diversifying selection. Allelic diversity may be favored not only through balancing or frequency-dependent selection at individual loci but also by associating different alleles in multicopy gene families. Using a combination of analytical calculations and simulations, we explored a population genetic model of epistatic selection and unequal recombination, where a trade-off exists between the benefit of allelic diversity and the cost of copy abundance. Starting from the neutral case, where we showed that gene copy number is Gamma distributed at equilibrium, we derived also the mean and shape of the limiting distribution under selection. Considering a more general model, which includes variable population size and population substructure, we explored by simulations mean fitness and some summary statistics of the copy number distribution. We determined the relative effects of selection, recombination, and demographic parameters in maintaining allelic diversity and shaping the mean fitness of a population. One way to control the variance of copy number is by lowering the rate of unequal recombination. Indeed, when encoding recombination by a rate modifier locus, we observe exactly this prediction. Finally, we analyzed the empirical copy number distribution of 3 genes in human and estimated recombination and selection parameters of our model.

Genetic Modeling and Genomic Analyses of Yearling Temperament in American Angus Cattle and Its Relationship With Productive Efficiency and Resilience Traits

Cattle temperament has been considered by farmers as a key breeding goal due to its relevance for cattlemen's safety, animal welfare, resilience, and longevity and its association with many economically important traits (e.g., production and meat quality). The definition of proper statistical models, accurate variance component estimates, and knowledge on the genetic background of the indicator trait evaluated are of great importance for accurately predicting the genetic merit of breeding animals. Therefore, 266,029 American Angus cattle with yearling temperament records (1-6 score) were used to evaluate statistical models and estimate variance components; investigate the association of sex and farm management with temperament; assess the weighted correlation of estimated breeding values for temperament and productive, reproductive efficiency and resilience traits; and perform a weighted single-step genome-wide association analysis using 69,559 animals genotyped for 54,609 single-nucleotide polymorphisms. Sex and extrinsic factors were significantly associated with temperament, including conception type, age of dam, birth season, and additional animal-human interactions. Similar results were observed among models including only the direct additive genetic effect and when adding other maternal effects. Estimated heritability of temperament was equal to 0.39 on the liability scale. Favorable genetic correlations were observed between temperament and other relevant traits, including growth, feed efficiency, meat quality, and reproductive traits. The highest approximated genetic correlations were observed between temperament and growth traits (weaning weight, 0.28; yearling weight, 0.28). Altogether, we identified 11 genomic regions, located across nine chromosomes including BTAX, explaining 3.33% of the total additive genetic variance. The candidate genes identified were enriched in pathways related to vision, which could be associated with reception of stimulus and/or cognitive abilities. This study encompasses large and diverse phenotypic, genomic, and pedigree datasets of US Angus cattle. Yearling temperament is a highly heritable and polygenic trait that can be improved through genetic selection. Direct selection for temperament is not expected to result in unfavorable responses on other relevant traits due to the favorable or low genetic correlations observed. In summary, this study contributes to a better understanding of the impact of maternal effects, extrinsic factors, and various genomic regions associated with yearling temperament in North American Angus cattle.

Variation in volatile organic compounds in Atlantic salmon mucus is associated with resistance to salmon lice infection

Salmon lice are ectoparasites that threaten wild and farmed salmonids. Artificial selection of salmon for resistance to the infectious copepodid lice stage currently relies on in vivo challenge trials on thousands of salmon a year. We challenged 5750 salmon with salmon lice (Lepeophtheirus salmonis) from two distinct farmed strains of salmon in two separate trials. We found that volatile organic compounds (VOC), 1-penten-3-ol, 1-octen-3-ol and 6-methyl-5-hepten-2-one in the mucus of the salmon host after salmon lice infection, were significantly associated with lice infection numbers across a range of water temperatures (5 °C, 10 °C, 17 °C). Some VOCs (benzene, 1-octen-3-ol and 3,5,5-trimethyl-2-hexene) were significantly different between lines divergently selected for salmon lice resistance. In a combined population assessment, selected VOCs varied between families in the range of 47- 59% indicating a genetic component and were positively correlated to the salmon hosts estimated breeding values 0.59–0.74. Mucosal VOC phenotypes could supplement current breeding practices and have the potential to be a more direct and ethical proxy for salmon lice resistance provided they can be measured prior to lice infestation.

Genome-wide association study of disease resilience traits from a natural polymicrobial disease challenge model in pigs identifies the importance of the major histocompatibility complex region

Infectious diseases cause tremendous financial losses in the pork industry, emphasizing the importance of disease resilience, which is the ability of an animal to maintain performance under disease. Previously, a natural polymicrobial disease challenge model was established, in which pigs were challenged in the late nursery phase by multiple pathogens to maximize expression of genetic differences in disease resilience. Genetic analysis found that performance traits in this model, including growth rate, feed and water intake, and carcass traits, as well as clinical disease phenotypes, were heritable and could be selected for to increase disease resilience of pigs. The objectives of the current study were to identify genomic regions that are associated with disease resilience in this model, using genome-wide association studies and fine-mapping methods, and to use gene set enrichment analyses to determine whether genomic regions associated with disease resilience are enriched for previously published quantitative trait loci, functional pathways, and differentially expressed genes subject to physiological states. Multiple quantitative trait loci were detected for all recorded performance and clinical disease traits. The major histocompatibility complex region was found to explain substantial genetic variance for multiple traits, including for growth rate in the late nursery (12.8%) and finisher (2.7%), for several clinical disease traits (up to 2.7%), and for several feeding and drinking traits (up to 4%). Further fine mapping identified 4 quantitative trait loci in the major histocompatibility complex region for growth rate in the late nursery that spanned the subregions for class I, II, and III, with 1 single-nucleotide polymorphism in the major histocompatibility complex class I subregion capturing the largest effects, explaining 0.8–27.1% of genetic variance for growth rate and for multiple clinical disease traits. This single-nucleotide polymorphism was located in the enhancer of TRIM39 gene, which is involved in innate immune response. The major histocompatibility complex region was pleiotropic for growth rate in the late nursery and finisher, and for treatment and mortality rates. Growth rate in the late nursery showed strong negative genetic correlations in the major histocompatibility complex region with treatment or mortality rates (−0.62 to −0.85) and a strong positive genetic correlation with growth rate in the finisher (0.79). Gene set enrichment analyses found genomic regions associated with resilience phenotypes to be enriched for previously identified disease susceptibility and immune capacity quantitative trait loci, for genes that were differentially expressed following bacterial or virus infection and immune response, and for gene ontology terms related to immune and inflammatory response. In conclusion, the major histocompatibility complex and other quantitative trait loci that harbor immune-related genes were identified to be associated with disease resilience traits in a large-scale natural polymicrobial disease challenge. The major histocompatibility complex region was pleiotropic for growth rate under challenge and for clinical disease traits. Four quantitative trait loci were identified across the class I, II, and III subregions of the major histocompatibility complex for nursery growth rate under challenge, with 1 single-nucleotide polymorphism in the major histocompatibility complex class I subregion capturing the largest effects. The major histocompatibility complex and other quantitative trait loci identified play an important role in host response to infectious diseases and can be incorporated in selection to improve disease resilience, in particular the identified single-nucleotide polymorphism in the major histocompatibility complex class I subregion.

The Chicken Embryo Model: A Novel and Relevant Model for Immune-Based Studies

Dysregulation of the immune system is associated with many pathologies, including cardiovascular diseases, diabetes, and cancer. To date, the most commonly used models in biomedical research are rodents, and despite the various advantages they offer, their use also raises numerous drawbacks. Recently, another in vivo model, the chicken embryo and its chorioallantoic membrane, has re-emerged for various applications. This model has many benefits compared to other classical models, as it is cost-effective, time-efficient, and easier to use. In this review, we explain how the chicken embryo can be used as a model for immune-based studies, as it gradually develops an embryonic immune system, yet which is functionally similar to humans'. We mainly aim to describe the avian immune system, highlighting the differences and similarities with the human immune system, including the repertoire of lymphoid tissues, immune cells, and other key features. We also describe the general in ovo immune ontogeny. In conclusion, we expect that this review will help future studies better tailor their use of the chicken embryo model for testing specific experimental hypotheses or performing preclinical testing.

Innate immune genes of the chicken MHC and related regions

Compared to the major histocompatibility complex (MHC) of typical mammals, the chicken BF/BL region is small and simple, with most of the genes playing central roles in the adaptive immune response. However, some genes of the chicken MHC are almost certainly involved in innate immunity, such as the complement component C4 and the lectin-like receptor/ligand gene pair BNK and Blec. The poorly expressed classical class I molecule BF1 is known to be recognised by natural killer (NK) cells and, analogous to mammalian immune responses, the classical class I molecules BF1 and BF2, the CD1 homologs and the butyrophilin homologs called BG may be recognised by adaptive immune lymphocytes with semi-invariant receptors in a so-called adaptate manner. Moreover, the TRIM and BG regions next to the chicken MHC, along with the genetically unlinked Y and olfactory/scavenger receptor regions on the same chromosome, have multigene families almost certainly involved in innate and adaptate responses. On this chicken microchromosome, the simplicity of the adaptive immune gene systems contrasts with the complexity of the gene systems potentially involved in innate immunity.

Comparison of major histocompatibility complex-B variability in Sri Lankan indigenous chickens with five global chicken populations using MHC-B SNP panel

In the present study, we investigated the major histocompatibility complex (MHC)-B haplotypes diversity of Sri Lankan indigenous chickens from three different geographical sites consisting of highly mixed populations using 90 SNPs in the MHC-B region. A total of 48 haplotypes were identified. Those included 37 novel haplotypes and 11 previously identified 'standard' haplotypes. The MHC-linked marker, LEI0258, had 23 alleles showing less diversity than defined by MHC-B SNP haplotypes. Among those identified haplotypes, five standard haplotypes-BSNP-O02, BSNP-M01, BSNP-A04, BSNP-K03, BSNP-T04-were most commonly observed, suggesting past introgression of imported breeds. Comparison of the MHC-B haplotypes of Sri Lankan and four other global populations with previously defined haplotypes indicated the sharing of 23 standard haplotypes with common origins. Novel haplotypes are population-specific and not shared among the geographical boundaries. Backyard indigenous chickens are unselected, highly crossbred, and generally thrive under dynamic environmental conditions. Hence free-range production systems may be responsible for maintaining high diversity in the MHC-B region with novel haplotypes.

A genome-wide scan to identify signatures of selection in two Iranian indigenous chicken ecotypes

Background

Various regions of the chicken genome have been under natural and artificial selection for thousands of years. The substantial diversity that exits among chickens from different geographic regions provides an excellent opportunity to investigate the genomic regions under selection which, in turn, will increase our knowledge about the mechanisms that underlie chicken diversity and adaptation. Several statistics have been developed to detect genomic regions that are under selection. In this study, we applied approaches based on differences in allele or haplotype frequencies (F_ST and hapFLK, respectively) between populations, differences in long stretches of consecutive homozygous sequences (ROH), and differences in allele frequencies within populations (composite likelihood ratio (CLR)) to identify inter- and intra-populations traces of selection in two Iranian indigenous chicken ecotypes, the Lari fighting chicken and the Khazak or creeper (short-leg) chicken.

Results

Using whole-genome resequencing data of 32 individuals from the two chicken ecotypes, approximately 11.9 million single nucleotide polymorphisms (SNPs) were detected and used in genomic analyses after quality processing. Examination of the distribution of ROH in the two populations indicated short to long ROH, ranging from 0.3 to 5.4 Mb. We found 90 genes that were detected by at least two of the four applied methods. Gene annotation of the detected putative regions under selection revealed candidate genes associated with growth (DCN, MEOX2 and CACNB1), reproduction (ESR1 and CALCR), disease resistance (S1PR1, ALPK1 and MHC-B), behavior pattern (AGMO, GNAO1 and PSEN1), and morphological traits (IHH and NHEJ1).

Conclusions

Our findings show that these two phenotypically different indigenous chicken populations have been under selection for reproduction, immune, behavioral, and morphology traits. The results illustrate that selection can play an important role in shaping signatures of differentiation across the genomic landscape of two chicken populations.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00664-9.

Identification of regeneration-involved growth factors in cartilage engineering procedure promotes its reconstruction

Aim: To fabricate mature cartilage for implantation, developmental biological processes and proteins should be understood and employed. Methods: A systems biology study of all protein-coding genes participating in cartilage regeneration resulted in a network graph with 11 nodes and 28 edges. Gene ontology and centrality analysis were performed based on the degree index. Results: The four most crucial biological processes along with the seven most interactive proteins involved in cartilage regeneration were identified. Some proteins, which are under serious discussion in cartilage developmental and disease processes, are included in regeneration. Conclusions: Findings positively correlate with the literature, supporting the use of the four most impressive proteins as growth factors applicable to cartilage tissue engineering, including COL2A1, SOX9, CTGF and TGFβ1.

Exonic SNP in MHC-DMB2 is associated with gene expression and humoral immunity in Japanese quails

The DMB2 gene is widely expressed at high levels in avian. This gene plays an important role in humoral immunity. The aim of this study was to investigate the effects of 361 G > C Single nucleotide polymorphism (SNP) on DMB2 protein structure and gene expression to determine how the 361 G > C SNP affects humoral immune response in Japanese quails. 0.2 mL of 5% sheep red blood cell (SRBC) was injected into breast muscle of 130 Japanese quails on 28 days. After DNA extraction, PCR was carried out to amplify a 333-base pair DNA fragment from the exon 2 of DMB2 gene. The pattern of all samples was determined through RFLP technique. PCR-RFLP results identified two alleles segregating (C, G) as three genotypes (CC, CG and GG) in Japanese Quails. The antibody response to SRBC with CC genotype was significantly higher than the CG and GG genotypes (P < 0.01). In silico analysis showed that the 361 G > C SNP has no effect on the physicochemical properties and 3D structure. The results of RT-qPCR indicated that the effect of genotype on gene expression is significant, so that the expression of CC genotype is more than CG and GG genotype. It can be inferred that the 361 G > C SNP in the exon 2 of MHC-DMB2 gene is not desirable. This mutation decreases humoral immune response by reducing DMB2 gene expression.

MHC-B variation in maternal and paternal synthetic lines of the Argentinian Campero INTA chicken

The Campero-INTA chicken of Argentina was developed to provide a robust bird that can survive under Argentinian pasture conditions with no significant additional nutrition, producing a source of animal protein for small producers or low-income families. In previous work, we described the AH paternal line of Campero and its Major Histocompatibility Complex B region (MHC-B) variation. In this work we analyzed the three remaining synthetic lines used to produce the Campero-INTA production bird: lines AS, A, and E. Because of the association between variation within the MHC of chickens and disease resistance, MHC variation within this breed is of particular interest. MHC variability within the lines used to produce the Campero-INTA chicken was examined using a 90 SNP panel encompassing the chicken MHC-B region plus the VNTR, LEI0258, located within the chicken MHC. Across all 4lines 12 haplotypes were found, with 7 of these being previously reported in North America/European breeds, reflecting the original breed sources for these birds. Three Campero unique haplotypes were found, 2 of which likely originated from MHC recombination events. MHC-B variation for all lines involved with production of the final Campero-INTA bird has now been determined.

Chickens as a simple system for scientific discovery: The example of the MHC

Chickens have played many roles in human societies over thousands of years, most recently as an important model species for scientific discovery, particularly for embryology, virology and immunology. In the last few decades, biomedical models like mice have become the most important model organism for understanding the mechanisms of disease, but for the study of outbred populations, they have many limitations. Research on humans directly addresses many questions about disease, but frank experiments into mechanisms are limited by practicality and ethics. For research into all levels of disease simultaneously, chickens combine many of the advantages of humans and of mice, and could provide an independent, integrated and overarching system to validate and/or challenge the dogmas that have arisen from current biomedical research. Moreover, some important systems are simpler in chickens than in typical mammals. An example is the major histocompatibility complex (MHC) that encodes the classical MHC molecules, which play crucial roles in the innate and adaptive immune systems. Compared to the large and complex MHCs of typical mammals, the chicken MHC is compact and simple, with single dominantly-expressed MHC molecules that can determine the response to infectious pathogens. As a result, some fundamental principles have been easier to discover in chickens, with the importance of generalist and specialist MHC alleles being the latest example.

hiphop: Improved paternity assignment among close relatives using a simple exclusion method for biallelic markers

Assignment of parentage with molecular markers is most difficult when the true parents have close relatives in the adult population. Here, we present an efficient solution to that problem by extending simple exclusion approaches to parentage analysis with single nucleotide polymorphic markers (SNPs). We augmented the previously published homozygote opposite test (hot), which counts mismatches due to the offspring and candidate parent having different homozygous genotypes, with an additional test. In this case, parents homozygous for the same SNP are incompatible with heterozygous offspring (i.e., "Homozygous Identical Parents, Heterozygous Offspring are Precluded": hiphop). We tested this approach in a cooperatively breeding bird, the superb fairy-wren, Malurus cyaneus, where rates of extra-pair paternity are exceptionally high, and where paternity assignment is challenging because breeding males typically have first-order adult relatives in their neighbourhood. Combining the tests and conditioning on the maternal genotype with a set of 1376 autosomal SNPs always allowed us to distinguish a single most likely sire from his relatives, and also to identify cases where the true sire must have been unsampled. In contrast, if just the hot test was used, we failed to identify a single most-likely sire in 2.5% of cases. Resampling enabled us to create guidelines for the number of SNPs required when first-order relatives coexist in the mating pool. Our method, implemented in the R package hiphop, therefore provides unambiguous parentage assignments even in systems with complex social organisation. We also identified a suite of Z- and W-linked SNPs that always identified sex correctly.

The dominantly expressed class II molecule from a resistant MHC haplotype presents only a few Marek's disease virus peptides by using an unprecedented binding motif

Viral diseases pose major threats to humans and other animals, including the billions of chickens that are an important food source as well as a public health concern due to zoonotic pathogens. Unlike humans and other typical mammals, the major histocompatibility complex (MHC) of chickens can confer decisive resistance or susceptibility to many viral diseases. An iconic example is Marek's disease, caused by an oncogenic herpesvirus with over 100 genes. Classical MHC class I and class II molecules present antigenic peptides to T lymphocytes, and it has been hard to understand how such MHC molecules could be involved in susceptibility to Marek's disease, given the potential number of peptides from over 100 genes. We used a new in vitro infection system and immunopeptidomics to determine peptide motifs for the 2 class II molecules expressed by the MHC haplotype B2, which is known to confer resistance to Marek's disease. Surprisingly, we found that the vast majority of viral peptide epitopes presented by chicken class II molecules arise from only 4 viral genes, nearly all having the peptide motif for BL2*02, the dominantly expressed class II molecule in chickens. We expressed BL2*02 linked to several Marek's disease virus (MDV) peptides and determined one X-ray crystal structure, showing how a single small amino acid in the binding site causes a crinkle in the peptide, leading to a core binding peptide of 10 amino acids, compared to the 9 amino acids in all other reported class II molecules. The limited number of potential T cell epitopes from such a complex virus can explain the differential MHC-determined resistance to MDV, but raises questions of mechanism and opportunities for vaccine targets in this important food species, as well as providing a basis for understanding class II molecules in other species including humans.

Comparison of Marek's disease virus challenge strains and bird types for vaccine licensing

Marek's disease virus (MDV) is an important poultry pathogen which is controlled through widespread vaccination with avirulent and attenuated strains, but continued evolution of field viruses to higher virulence has required ongoing improvement of available vaccine strains, and these vaccine strains also offer an attractive platform for designing recombinant vector vaccines with cross-protection against MDV and additional pathogens.  Recent reports of failures in vaccine licensing trials of positive controls to reach appropriately high levels of MD incidence prompted us to evaluate possible combinations of outbred specific pathogen-free (SPF) layer lines and alternative virulent challenge strains which could provide more consistent models for serotype-3 vectored vaccine development.  Choice of layer line and virulent MDV challenge strain each contributed to the ability of a challenge model to reach 80 percent virulence in unvaccinated positive control groups in the majority of trials without overwhelming serotype-3 vectored vaccine protection in vaccinated groups.  Conversely, reducing challenge virus dose by a factor of four, or vaccine dose by half, had no consistent effect across these models.  Although MDV strain 617A had the most potential as an alternative to strains that are currently approved for licensing trials, no combination of layer line and challenge virus consistently met the goals for a successful challenge model in all study replicates, indicating that high variability is an inherent difficulty in MDV challenge studies, at least when outbred birds are used.

Genome-wide fine-mapping identifies pleiotropic and functional variants that predict many traits across global cattle populations

The difficulty in finding causative mutations has hampered their use in genomic prediction. Here, we present a methodology to fine-map potentially causal variants genome-wide by integrating the functional, evolutionary and pleiotropic information of variants using GWAS, variant clustering and Bayesian mixture models. Our analysis of 17 million sequence variants in 44,000+ Australian dairy cattle for 34 traits suggests, on average, one pleiotropic QTL existing in each 50 kb chromosome-segment. We selected a set of 80k variants representing potentially causal variants within each chromosome segment to develop a bovine XT-50K genotyping array. The custom array contains many pleiotropic variants with biological functions, including splicing QTLs and variants at conserved sites across 100 vertebrate species. This biology-informed custom array outperformed the standard array in predicting genetic value of multiple traits across populations in independent datasets of 90,000+ dairy cattle from the USA, Australia and New Zealand.

From Chickens to Humans: The Importance of Peptide Repertoires for MHC Class I Alleles

In humans, killer immunoglobulin-like receptors (KIRs), expressed on natural killer (NK) and thymus-derived (T) cells, and their ligands, primarily the classical class I molecules of the major histocompatibility complex (MHC) expressed on nearly all cells, are both polymorphic. The variation of this receptor-ligand interaction, based on which alleles have been inherited, is known to play crucial roles in resistance to infectious disease, autoimmunity, and reproduction in humans. However, not all the variation in response is inherited, since KIR binding can be affected by a portion of the peptide bound to the class I molecules, with the particular peptide presented affecting the NK response. The extent to which the large multigene family of chicken immunoglobulin-like receptors (ChIRs) is involved in functions similar to KIRs is suspected but not proven. However, much is understood about the two MHC-I molecules encoded in the chicken MHC. The BF2 molecule is expressed at a high level and is thought to be the predominant ligand of cytotoxic T lymphocytes (CTLs), while the BF1 molecule is expressed at a much lower level if at all and is thought to be primarily a ligand for NK cells. Recently, a hierarchy of BF2 alleles with a suite of correlated properties has been defined, from those expressed at a high level on the cell surface but with a narrow range of bound peptides to those expressed at a lower level on the cell surface but with a very wide repertoire of bound peptides. Interestingly, there is a similar hierarchy for human class I alleles, although the hierarchy is not as wide. It is a question whether KIRs and ChIRs recognize class I molecules with bound peptide in a similar way, and whether fastidious to promiscuous hierarchy of class I molecules affect both T and NK cell function. Such effects might be different from those predicted by the similarities of peptide-binding based on peptide motifs, as enshrined in the idea of supertypes. Since the size of peptide repertoire can be very different for alleles with similar peptide motifs from the same supertype, the relative importance of these two properties may be testable.

The Chicken MHC: Insights into Genetic Resistance, Immunity, and Inflammation Following Infectious Bronchitis Virus Infections

The chicken immune system has provided an immense contribution to basic immunology knowledge by establishing major landmarks and discoveries that defined concepts widely used today. One of many special features on chickens is the presence of a compact and simple major histocompatibility complex (MHC). Despite its simplicity, the chicken MHC maintains the essential counterpart genes of the mammalian MHC, allowing for a strong association to be detected between the MHC and resistance or susceptibility to infectious diseases. This association has been widely studied for several poultry infectious diseases, including infectious bronchitis. In addition to the MHC and its linked genes, other non-MHC loci may play a role in the mechanisms underlying such resistance. It has been reported that innate immune responses, such as macrophage function and inflammation, might be some of the factors driving resistance or susceptibility, consequently influencing the disease outcome in an individual or a population. Information about innate immunity and genetic resistance can be helpful in developing effective preventative measures for diseases such as infectious bronchitis, to which a systemic antibody response is often not associated with disease protection. In this review, we summarize the importance of the chicken MHC in poultry disease resistance, particularly to infectious bronchitis virus (IBV) infections and the role played by innate immunity and inflammation on disease outcome. We highlight how future studies focusing on the MHC and non-MHC genes can potentially bring clarity to observed resistance in some chicken B haplotype lines.

Major histocompatibility complex B variability in Korean native chicken breeds

Adaptive genetic variations have direct influences on the fitness traits of the animal. The major histocompatibility complex B (MHC-B) region is responsible for adaptive and innate immune responses in chickens. In native Korean chicken breeds, no information on serologically defined B haplotypes is available. We investigated the MHC-B diversity in 5 restored lines of Korean native chicken and Ogye chicken breeds using a recently described MHC-B single-nucleotide polymorphism (SNP) panel and the MHC-linked LEI0258 variable number of tandem repeat marker. High SNP haplotype diversity was observed in Korean native chicken breeds with an average of 9.7 MHC-B SNP haplotypes per line. The total number of haplotypes ranged from 6 to 12 per line, and population-specific haplotypes ranged from 3 to 4. A total of 41 BSNP haplotypes, including 26 novel population-specific haplotypes and 15 common haplotypes, were reported over all populations. The 15 common haplotypes included 7 novel and 8 previously reported standard haplotypes. Selection and breeding evidence supports the observation of common haplotypes between the Korean native chicken and exotic breeds. Similarly, the LEI0258 marker showed allele variation, between 193 bp and 474 bp having 5 to 8 alleles per population. Some of these alleles (193, 249, 309, and 443 bp) were shared and more frequently observed. Comparison between SNP haplotypes and LEI0258 allele sizes for the same samples showed that some LEI0258 allele sizes correspond to more than one BSNP haplotype. The use of the MHC-B SNP panel greatly enhances the identification of MHC diversity compared with the sole use of the LEI0258 marker in native chicken populations.

Mapping QTL Associated with Resistance to Avian Oncogenic Marek's Disease Virus (MDV) Reveals Major Candidate Genes and Variants

Marek’s disease (MD) represents a significant global economic and animal welfare issue. Marek’s disease virus (MDV) is a highly contagious oncogenic and highly immune-suppressive α-herpes virus, which infects chickens, causing neurological effects and tumour formation. Though partially controlled by vaccination, MD continues to have a profound impact on animal health and on the poultry industry. Genetic selection provides an alternative and complementary method to vaccination. However, even after years of study, the genetic mechanisms underlying resistance to MDV remain poorly understood. The Major Histocompatability Complex (MHC) is known to play a role in disease resistance, along with a handful of other non-MHC genes. In this study, one of the largest to date, we used a multi-facetted approach to identify quantitative trait locus regions (QTLR) influencing resistance to MDV, including an F₆ population from a full-sib advanced intercross line (FSIL) between two elite commercial layer lines differing in resistance to MDV, RNA-seq information from virus challenged chicks, and genome wide association study (GWAS) from multiple commercial lines. Candidate genomic elements residing in the QTLR were further tested for association with offspring mortality in the face of MDV challenge in eight pure lines of elite egg-layer birds. Thirty-eight QTLR were found on 19 chicken chromosomes. Candidate genes, microRNAs, long non-coding RNAs and potentially functional mutations were identified in these regions. Association tests were carried out in 26 of the QTLR, using eight pure lines of elite egg-layer birds. Numerous candidate genomic elements were strongly associated with MD resistance. Genomic regions significantly associated with resistance to MDV were mapped and candidate genes identified. Various QTLR elements were shown to have a strong genetic association with resistance. These results provide a large number of significant targets for mitigating the effects of MDV infection on both poultry health and the economy, whether by means of selective breeding, improved vaccine design, or gene-editing technologies.

Genetic diversity of MHC-B in 12 chicken populations in Korea revealed by single-nucleotide polymorphisms

This study used a single-nucleotide polymorphism (SNP) panel to characterise the diversity in the major histocompatibility complex B region (MHC-B) in 12 chicken populations in Korea. Samples were genotyped for 96 MHC-B SNPs using an Illumina GoldenGate genotyping assay. The MHC-B SNP haplotypes were predicted using 58 informative SNPs and a coalescence-based Bayesian algorithm implemented by the PHASE program and a manual curation process. In total, 117 haplotypes, including 24 shared and 93 unique haplotypes, were identified. The unique haplotype numbers ranged from 0 in Rhode Island Red to 32 in the Korean native commercial chicken population 2 ("Hanhyup-3ho"). Population and haplotype principal component analysis (PCA) indicated no clear population structure based on the MHC haplotypes. Three haplotype clusters (A, B, C) segregated in these populations highlighted the relationship between the haplotypes in each cluster. The sequences from two clusters (B and C) overlapped, whereas the sequences from the third cluster (A) were very different. Overall, native breeds had high genetic diversity in the MHC-B region compared with the commercial breeds. This highlights their immune capabilities and genetic potential for resistance to many different pathogens.

Performances of Adaptive MultiBLUP, Bayesian regressions, and weighted-GBLUP approaches for genomic predictions in Belgian Blue beef cattle

BACKGROUND

Genomic selection has been successfully implemented in many livestock and crop species. The genomic best linear unbiased predictor (GBLUP) approach, assigning equal variance to all SNP effects, is one of the reference methods. When large-effect variants contribute to complex traits, it has been shown that genomic prediction methods that assign a higher variance to subsets of SNP effects can achieve higher prediction accuracy. We herein compared the efficiency of several such approaches, including the Adaptive MultiBLUP (AM-BLUP) that uses local genomic relationship matrices (GRM) to automatically identify and weight genomic regions with large effects, to predict genetic merit in Belgian Blue beef cattle.

RESULTS

We used a population of approximately 10,000 genotyped cows and their phenotypes for 14 traits, mostly related to muscular development and body dimensions. According to the trait, we found that 4 to 25% of the genetic variance could be associated with 2 to 12 genomic regions harbouring large-effect variants. Noteworthy, three previously identified recessive deleterious variants presented heterozygote advantage and were among the most significant SNPs for several traits. The AM-BLUP resulted in increased reliability of genomic predictions compared to GBLUP (+ 2%), but Bayesian methods proved more efficient (+ 3%). Overall, the reliability gains remained thus limited although higher gains were observed for skin thickness, a trait affected by two genomic regions having particularly large effects. Higher accuracies than those from the original AM-BLUP were achieved when applying the Bayesian Sparse Linear Mixed Model to pre-select groups of SNPs with large effects and subsequently use their estimated variance to build a weighted GRM. Finally, the single-step GBLUP performed best and could be further improved (+ 3% prediction accuracy) by using these weighted GRM.

CONCLUSIONS

The AM-BLUP is an attractive method to automatically identify and weight genomic regions with large effects on complex traits. However, the method was less accurate than Bayesian methods. Overall, weighted methods achieved modest accuracy gains compared to GBLUP. Nevertheless, the computational efficiency of the AM-BLUP might be valuable at higher marker density, including with whole-genome sequencing data. Furthermore, weighted GRM are particularly useful to account for large variance loci in the single-step GBLUP.

Assessing MHC-B diversity in Silkie chickens

The major histocompatibility complex (MHC) is a highly polymorphic region on chromosome 16, which contains numerous immune response genes, and is known to influence disease susceptibility and resistance in chickens. Variability of MHC-B haplotypes in various well-known and commercially utilized breeds has previously been identified. This study aims to understand MHC-B diversity in the Silkie breed using a high-density SNP panel that encompasses the chicken MHC-B region. DNA was obtained from 74 females and 27 males from a commercial Silkie breeder colony that is maintained through minimal genetic selection practices. A previously described panel of 90 SNPs, all located within the MHC-B region, was used to evaluate MHC-B variability in the commercial Silkie breeder colony. MHC-B haplotypes identified from the individual SNP information in the Silkie colony were compared to published haplotypes from the same region. Of the 27 haplotypes identified in the Silkie population, 8 have been previously described. Nineteen haplotypes are unique to the Silkie population and include one novel recombinant and 2 additional possible novel recombinants. Six haplotypes were found at a frequency greater than 5% of the population, of which 4 are novel. Finally, Hardy Weinberg Equilibrium (HWE) was calculated for the observed haplotypes, which were found to be in HWE. This study shows considerable MHC-B diversity in the Silkie breed and adds further information on variability of the MHC-B region in the chicken.

Advances in methodologies for detecting MHC-B variability in chickens

The chicken major histocompatibility B complex (MHC-B) region is of great interest owing to its very strong association with resistance to many diseases. Variation in the MHC-B was initially identified by hemagglutination of red blood cells with specific alloantisera. New technologies, developed to identify variation in biological materials, have been applied to the chicken MHC. Protein variation encoded by the MHC genes was examined by immunoprecipitation and 2-dimensional gel electrophoresis. Increased availability of DNA probes, PCR, and sequencing resulted in the application of DNA-based methods for MHC detection. The chicken reference genome, completed in 2004, allowed further refinements in DNA methods that enabled more rapid examination of MHC variation and extended such analyses to include very diverse chicken populations. This review progresses from the inception of MHC-B identification to the present, describing multiple methods, plus their advantages and disadvantages.

Mhc-B haplotypes in "Campero-Inta" chicken synthetic line

The major histocompatibility complex-B (MHC-B) in chickens is a cluster of genes located on chromosome 16. The chicken MHC-B is known to be highly associated with resistance to numerous diseases caused by viruses, bacteria, and parasitic pathogens. Since the level of resistance varies with MHC-B haplotypes, identification and classification of different haplotypes within lines is important for sustaining lines. The "Campero-INTA" chicken breed is a meat-type free-range poultry breed that was developed specifically for small producers in Argentina. Campero-INTA was started by selection in populations produced by crosses between a variety of established lines. MHC-B variation was examined in 65 samples obtained in 2002 using the VNTR marker LEI0258, a marker for MHC-B region. These samples plus and an additional 55 samples from 2018 were examined for variation using the MHC-B specific SNP panel that encompasses ∼230,000 bp of the MHC-B region. Eleven MHC-B SNP haplotypes with 6 LEI0258 alleles were identified in the 120 samples representing the Campero-INTA AH (male) line. Seven haplotypes originate from the breeds originally used in the development of Campero-INTA AH line. Two appear to be recombinant haplotypes. The origin of the remaining 2 is not known, but may be associated with genes introduced from crosses with the Fayoumi breed conducted more recently to sustain the line.

Development and optimization of a hybridization technique to type the classical class I and class II B genes of the chicken MHC

The classical class I and class II molecules of the major histocompatibility complex (MHC) play crucial roles in immune responses to infectious pathogens and vaccines as well as being important for autoimmunity, allergy, cancer and reproduction. These classical MHC genes are the most polymorphic known, with roughly 10,000 alleles in humans. In chickens, the MHC (also known as the BF-BL region) determines decisive resistance and susceptibility to infectious pathogens, but relatively few MHC alleles and haplotypes have been described in any detail. We describe a typing protocol for classical chicken class I (BF) and class II B (BLB) genes based on a hybridization method called reference strand-mediated conformational analysis (RSCA). We optimize the various steps, validate the analysis using well-characterized chicken MHC haplotypes, apply the system to type some experimental lines and discover a new chicken class I allele. This work establishes a basis for typing the MHC genes of chickens worldwide and provides an opportunity to correlate with microsatellite and with single nucleotide polymorphism (SNP) typing for approaches involving imputation.

A genomic inference of the White Plymouth Rock genealogy

Crossing of populations has been, and still is, a central component in domestication and breed and variety formation. It is a way for breeders to utilize heterosis and to introduce new genetic variation into existing plant and livestock populations. During the mid-19th century, several chicken breeds that had been introduced to America from Europe and Asia became the founders for those formed in the USA. Historical records about the genealogy of these populations are often unclear and inconsistent. Here, we used genomics in an attempt to describe the ancestry of the White Plymouth Rock (WPR) chicken. In total, 150 chickens from the WPR and 8 other stocks that historical records suggested contributed to its formation were whole-genome re-sequenced. The admixture analyses of the autosomal and sex chromosomes showed that the WPR was likely founded as a cross between a paternal lineage that was primarily Dominique, and a maternal lineage where Black Java and Cochin contributed in essentially equal proportions. These results were consistent and provided quantification with the historical records that they were the main contributors to the WPR. The genomic analyses also revealed genome-wide contributions (<10% each) by Brahma, Langshan, and Black Minorca. When viewed on an individual chromosomal basis, contributions varied considerably among stocks.

Signatures of selection in the genome of Swedish warmblood horses selected for sport performance

Background

A growing demand for improved physical skills and mental attitude in modern sport horses has led to strong selection for performance in many warmblood studbooks. The aim of this study was to detect genomic regions with low diversity, and therefore potentially under selection, in Swedish Warmblood horses (SWB) by analysing high-density SNP data. To investigate if such signatures could be the result of selection for equestrian sport performance, we compared our SWB SNP data with those from Exmoor ponies, a horse breed not selected for sport performance traits.

Results

The genomic scan for homozygous regions identified long runs of homozygosity (ROH) shared by more than 85% of the genotyped SWB individuals. Such ROH were located on ECA4, ECA6, ECA7, ECA10 and ECA17. Long ROH were instead distributed evenly across the genome of Exmoor ponies in 77% of the chromosomes. Two population differentiation tests (F_ST and XP-EHH) revealed signatures of selection on ECA1, ECA4, and ECA6 in SWB horses.

Conclusions

Genes related to behaviour, physical abilities and fertility, appear to be targets of selection in the SWB breed. This study provides a genome-wide map of selection signatures in SWB horses, and ground for further functional studies to unravel the biological mechanisms behind complex traits in horses.

Exploration of fine-scale recombination rate variation in the domestic horse

Total genetic map length and local recombination landscapes typically vary within and across populations. As a first step to understanding the recombination landscape in the domestic horse, we calculated population recombination rates and identified likely recombination hotspots using approximately 1.8 million SNP genotypes for 485 horses from 32 distinct breeds. The resulting breed-averaged recombination map spans 2.36 Gb and accounts for 2939.07 cM. Recombination hotspots occur once per 23.8 Mb on average and account for ∼9% of the physical map length. Regions with elevated recombination rates in the entire cohort were enriched for genes in pathways involving interaction with the environment: immune system processes (specifically, MHC class I and class II genes), responses to stimuli, and serotonin receptor pathways. We found significant correlations between differences in local recombination rates and population differentiation quantified by F _ST Analysis of breed-specific maps revealed thousands of hotspot regions unique to particular breeds, as well as unique "coldspots," regions where a particular breed showed below-average recombination, whereas all other breeds had evidence of a hotspot. Finally, we identified relative enrichment (P = 5.88 × 10^-27) for the in silico-predicted recognition motif for equine PR/SET domain 9 (PRDM9) in recombination hotspots. These results indicate that selective pressures and PRDM9 function contribute to variation in recombination rates across the domestic horse genome.

Genetic diversity of 10 indigenous chicken ecotypes from Southern Highlands of Tanzania based on Major Histocompatibility Complex-linked microsatellite LEI0258 marker typing

Unraveling the genetic diversity of livestock species is central to understanding their value and importance for conservation and improvement in diverse production environments. In developing countries, information on genetic attributes of many livestock species is unfortunately scanty to support well-informed decision-making upon relevant management strategies. This study aimed at investigating allelic variability, genetic diversity, and genetic relationships of 10 indigenous chicken ecotypes from Southern Highlands of Tanzania using the Major Histocompatibility Complex-linked LEI0258 marker. A total of 400 DNA samples, 40 per ecotype, were genotyped by capillary electrophoresis. Thirty different alleles with sizes ranging from 197 to 569 bp were determined. The number of alleles ranged from 17 (Itunduma) to 21 (Mbeya), with an average of 19.20 alleles per ecotype. Allelic polymorphism was further evaluated through genotyping by Sanger sequencing. Thirty-three DNA samples with different fragment sizes were re-amplified and their alleles sequenced to depict polymorphism based on a combination of two repeat regions at 12 and 13 bp, respectively, and flanking regions with SNP and indels. The repeat region at 13 bp appeared 1 to 28 times, whereas the region at 12 bp appeared 3 to 19 times in all sequenced fragments. The numbers of indels and SNP determined were 7 and 9, respectively. From capillary electrophoresis, the Chunya and Msimbazi ecotypes exhibited the highest genetic diversity (0.937), whereas the lowest value (0.910) was observed from the Mbarali ecotype, with an average of 0.925. The Namtumbo and Wanging'ombe ecotypes showed high inbreeding coefficients (FIS > 0.05), whereas a high excess heterozygote value (FIS = -0.098) was observed from the Njombe ecotype. Two percent of the genetic diversity was due to differences among ecotypes, and the rest was due to differences among individuals within the ecotypes. Despite the overall low genetic differentiation, both fragment and sequencing analyses depicted a high allelic and genetic variability across 10 chicken ecotypes. These results therefore, underscore the importance of establishing appropriate conservation and management strategies to capitalize on observed variability and maintain genetic flexibility across diverse production environments.

Association of Candidate Genes with Response to Heat and Newcastle Disease Virus

Newcastle disease is considered the number one disease constraint to poultry production in low and middle-income countries, however poultry that is raised in resource-poor areas often experience multiple environmental challenges. Heat stress has a negative impact on production, and immune response to pathogens can be negatively modulated by heat stress. Candidate genes and regions chosen for this study were based on previously reported associations with response to immune stimulants, pathogens, or heat, including: TLR3, TLR7, MX, MHC-B (major histocompatibility complex, gene complex), IFI27L2, SLC5A1, HSPB1, HSPA2, HSPA8, IFRD1, IL18R1, IL1R1, AP2A2, and TOLLIP. Chickens of a commercial egg-laying line were infected with a lentogenic strain of NDV (Newcastle disease virus); half the birds were maintained at thermoneutral temperature and the other half were exposed to high ambient temperature before the NDV challenge and throughout the remainder of the study. Phenotypic responses to heat, to NDV, or to heat + NDV were measured. Selected SNPs (single nucleotide polymorphisms) within 14 target genes or regions were genotyped; and genotype effects on phenotypic responses to NDV or heat + NDV were tested in each individual treatment group and the combined groups. Seventeen significant haplotype effects, among seven genes and seven phenotypes, were detected for response to NDV or heat or NDV + heat. These findings identify specific genetic variants that are associated with response to heat and/or NDV which may be useful in the genetic improvement of chickens to perform favorably when faced with pathogens and heat stress.

Detection of copy number variations in brown and white layers based on genotyping panels with different densities

Background

Copy number variations (CNV) are an important source of genetic variation that has gained increasing attention over the last couple of years. In this study, we performed CNV detection and functional analysis for 18,719 individuals from four pure lines and one commercial cross of layer chickens. Samples were genotyped on four single nucleotide polymorphism (SNP) genotyping platforms, i.e. the Illumina 42K, Affymetrix 600K, and two different customized Affymetrix 50K chips. CNV recovered from the Affymetrix chips were identified by using the Axiom^® CNV Summary Tools and PennCNV software and those from the Illumina chip were identified by using the cnvPartition in the Genome Studio software.

Results

The mean number of CNV per individual varied from 0.50 to 4.87 according to line or cross and size of the SNP genotyping set. The length of the detected CNV across all datasets ranged from 1.2 kb to 3.2 Mb. The number of duplications exceeded the number of deletions for most lines. Between the lines, there were considerable differences in the number of detected CNV and their distribution. Most of the detected CNV had a low frequency, but 19 CNV were identified with a frequency higher than 5% in birds that were genotyped on the 600K panel, with the most common CNV being detected in 734 birds from three lines.

Conclusions

Commonly used SNP genotyping platforms can be used to detect segregating CNV in chicken layer lines. The sample sizes for this study enabled a detailed characterization of the CNV landscape within commercially relevant lines. The size of the SNP panel used affected detection efficiency, with more CNV detected per individual on the higher density 600K panel. In spite of the high level of inter-individual diversity and a large number of CNV observed within individuals, we were able to detect 19 frequent CNV, of which, 57.9% overlapped with annotated genes and 89% overlapped with known quantitative trait loci.

Electronic supplementary material

The online version of this article (10.1186/s12711-018-0428-4) contains supplementary material, which is available to authorized users.