Genetical genomics of growth in a chicken model

Complex genetic architecture of the chicken Growth1 QTL region

The genetic complexity of polygenic traits represents a captivating and intricate facet of biological inheritance. Unlike Mendelian traits controlled by a single gene, polygenic traits are influenced by multiple genetic loci, each exerting a modest effect on the trait. This cumulative impact of numerous genes, interactions among them, environmental factors, and epigenetic modifications results in a multifaceted architecture of genetic contributions to complex traits. Given the well-characterized genome, diverse traits, and range of genetic resources, chicken (Gallus gallus) was employed as a model organism to dissect the intricate genetic makeup of a previously identified major Quantitative Trait Loci (QTL) for body weight on chromosome 1. A multigenerational advanced intercross line (AIL) of 3215 chickens whose genomes had been sequenced to an average of 0.4x was analyzed using genome-wide association study (GWAS) and variance-heterogeneity GWAS (vGWAS) to identify markers associated with 8-week body weight. Additionally, epistatic interactions were studied using the natural and orthogonal interaction (NOIA) model. Six genetic modules, two from GWAS and four from vGWAS, were strongly associated with the studied trait. We found evidence of both additive- and non-additive interactions between these modules and constructed a putative local epistasis network for the region. Our screens for functional alleles revealed a missense variant in the gene ribonuclease H2 subunit B (RNASEH2B), which has previously been associated with growth-related traits in chickens and Darwin's finches. In addition, one of the most strongly associated SNPs identified is located in a non-coding region upstream of the long non-coding RNA, ENSGALG00000053256, previously suggested as a candidate gene for regulating chicken body weight. By studying large numbers of individuals from a family material using approaches to capture both additive and non-additive effects, this study advances our understanding of genetic complexities in a highly polygenic trait and has practical implications for poultry breeding and agriculture.

Elucidation of the genetic determination of body weight and size in Chinese local chicken breeds by large-scale genomic analyses

BACKGROUND

Body weight and size are important economic traits in chickens. While many growth-related quantitative trait loci (QTLs) and candidate genes have been identified, further research is needed to confirm and characterize these findings. In this study, we investigate genetic and genomic markers associated with chicken body weight and size. This study provides new insights into potential markers for genomic selection and breeding strategies to improve meat production in chickens.

METHODS

We performed whole-genome resequencing of and Wenshang Barred (WB) chickens (n = 596) and three additional breeds with varying body sizes (Recessive White (RW), WB, and Luxi Mini (LM) chickens; (n = 50)). We then used selective sweeps of mutations coupled with genome-wide association study (GWAS) to identify genomic markers associated with body weight and size.

RESULTS

We identified over 9.4 million high-quality single nucleotide polymorphisms (SNPs) among three chicken breeds/lines. Among these breeds, 287 protein-coding genes exhibited positive selection in the RW and WB populations, while 241 protein-coding genes showed positive selection in the LM and WB populations. Genomic heritability estimates were calculated for 26 body weight and size traits, including body weight, chest breadth, chest depth, thoracic horn, body oblique length, keel length, pelvic width, shank length, and shank circumference in the WB breed. The estimates ranged from 0.04 to 0.67. Our analysis also identified a total of 2,522 genome-wide significant SNPs, with 2,474 SNPs clustered around two genomic regions. The first region, located on chromosome 4 (7.41-7.64 Mb), was linked to body weight after ten weeks and body size traits. LCORL, LDB2, and PPARGC1A were identified as candidate genes in this region. The other region, located on chromosome 1 (170.46-171.53 Mb), was associated with body weight from four to eighteen weeks and body size traits. This region contained CAB39L and WDFY2 as candidate genes. Notably, LCORL, LDB2, and PPARGC1A showed highly selective signatures among the three breeds of chicken with varying body sizes.

CONCLUSION

Overall this study provides a comprehensive map of genomic variants associated with body weight and size in chickens. We propose two genomic regions, one on chromosome 1 and the other on chromosome 4, that could helpful for developing genome selection breeding strategies to enhance meat yield in chickens.

Association analysis of production traits of Japanese quail (Coturnix japonica) using restriction-site associated DNA sequencing

This study was designed to perform an association analysis and identify SNP markers associated with production traits of Japanese quail using restriction-site-associated DNA sequencing. Weekly body weight data from 805 quail were collected from hatching to 16 weeks of age. A total number of 3990 eggs obtained from 399 female quail were used to assess egg quality traits. Egg-related traits were measured at the beginning of egg production (first stage) and at 12 weeks of age (second stage). Five eggs were analyzed at each stage. Traits, such as egg weight, egg length and short axes, eggshell strength and weight, egg equator thickness, yolk weight, diameter, and colour, albumen weight, age of first egg, total number of laid eggs, and egg production rate, were assessed. A total of 383 SNPs and 1151 associations as well as 734 SNPs and 1442 associations were identified in relation to quail production traits using general linear model (GLM) and mixed linear model (MLM) approaches, respectively. The GLM-identified SNPs were located on chromosomes 1-13, 15, 17-20, 24, 26-28, and Z, underlying phenotypic traits, except for egg and albumen weight at the first stage and yolk yellowness at the second stage. The MLM-identified SNPs were positioned on defined chromosomes associated with phenotypic traits except for the egg long axis at the second stage of egg production. Finally, 35 speculated genes were identified as candidate genes for the targeted traits based on their nearest positions. Our findings provide a deeper understanding and allow a more precise genetic improvement of production traits of Galliformes, particularly in Japanese quail.

Sex proportion as a covariate increases the statistical test power in growth performance based experiments using as-hatched broilers

The availability of sexed day-old broiler chicks is becoming an issue as feather sexing is no longer possible. This has great implications for broiler researchers as the use of randomly distributed mixed-sex birds may result in a greater between-pen variation and thus less statistical power than the use of single-sex birds. The objective of this study was to evaluate the effect of including sex proportion as a covariate in an analysis of covariance (ANCOVA) on the statistical power compared to analysis of variance (ANOVA) where sex was not considered. The statistical parameters examined include mean square error (MSE), the F-statistic, model fit, model significance and observed power. A total of 4 separate experiments that used mixed-sex broilers with unequal numbers of male and female birds per pen were conducted during which performance of the birds was measured. The male % in each pen was recorded during each experiment and corrected for mortality. The performance results were analysed by ANOVA and the statistical parameters were then compared to ANCOVA where sex proportion was included as a covariate. The results showed that a set of assumptions first needed to be met to run ANCOVA. In addition, if the ANOVA results show a high level of model significance and power, then ANCOVA may not be necessary. In other circumstances where the assumptions are met and model significance and observed power are low, the inclusion of sex proportion as a covariate in the analysis will help to reduce MSE, increase the F-statistic value and improve the model significance, model fit and observed power. Therefore, it is suggested that sex proportion should be considered as a covariate in ANCOVA to improve statistical power in nutritional experiments when male and female broilers are unequally and randomly distributed in pens.

Growth performance, survivability and profitability of improved smallholder chicken genetics in Nigeria: A COVID-19 intervention study

The impact of COVID-19 pandemic on smallholder farming households (SFH) includes increased poverty, and loss of livelihoods. Provision of livestock to SFH is a helpful intervention to mitigate this impact. This study provided a total of 150 smallholder poultry farmers, randomly selected from three states (Kebbi, Nasarawa, and Imo) in Nigeria, with ten 5-week-old chickens (mixed sexes) each, of either FUNAAB Alpha or Noiler chicken genetics. The improved, dual-purpose chickens were evaluated for growth performance (GP), survivability and profitability. The birds were managed under semi-scavenging production system. Body weight, mortality, and cost of production (COP) were recorded every 4 weeks until 21 weeks of age. Profitability was a function of the COP, and the selling price for live-birds (cocks). Body weight of Noiler (1,927 g) birds was not significantly (p > .05) higher than FUNAAB Alpha (1,792 g) at 21 weeks. Agroecology and genetics had significant (p < .05) effects on GP and survivability. Survivability of FUNAAB Alpha was higher (p < .05) than Noiler, with Nasarawa (81%-96%), having the highest (p < .0001) survival rate compared to Imo (62%-81%), and Kebbi (58%-75%). At 21 weeks, the number of cocks and hens differed significantly (p < .05) within the states (Imo: 2.4 ± .2 and 5.4 ± .3; Kebbi: 2.6 ± .2 and 5.5 ± .3; and Nasarawa: 2.9 ± .2 and 5.8 ± .3). Nasarawa (NGN 7,808; USD 19) ranked best for profitability, followed by Kebbi (NGN 6,545; USD 16) and Imo (NGN 5,875; USD 14). Overall, this study demonstrates that provision of improved chickens to vulnerable SFH in Nigeria holds great potential for economic growth, and resilience during emergencies, such as the COVID-19 pandemic.

Blood transcriptome comparison between sexes and their function in 4-week Rhode Island red chickens

Sex is a major biological factor in the development and physiology of a sexual reproductive organism, and its role in the growing process is needed to be investigated in various species. We compare blood transcriptome between 5 males and 5 females in 4-week-old Rhode Island Red chickens and perform functional annotation of differentially expressed genes (DEGs). The results are as follows. 141 and 109 DEGs were located in autosomes and sex chromosomes, respectively. The gene ontology (GO) terms are significantly (p < 0.05) enriched, which were limb development, inner ear development, positive regulation of dendrite development, the KEGG pathway the TGF-beta signaling pathway, and melanogenesis (p < 0.05). These pathways are related to morphological maintenance and growth of the tissues. In addition, the SMAD2W and the BMP5 were involved in the TGF-beta signaling pathway, and both play an important role in maintaining tissue development. The major DEGs related to the development of neurons and synapses include the up-regulated NRN1, GDF10, SLC1A1, BMP5, NBEA, and NRXN1. Also, 7 DEGs were validated using RT-qPCR with high correlation (r ² = 0.74). In conclusion, the differential expression of blood tissue in the early growing chicken was enriched in TGF-beta signaling and related to the development of neurons and synapses including SMAD2W and BMP5. These results suggest that blood in the early growing stage is differentially affected in tissue development, nervous system, and pigmentation by sex. For future research, experimental characterization of DEGs and a holistic investigation of various tissues and growth stages will be required.

High-density genotyping reveals candidate genomic regions for chicken body size in breeds of Asian origin

Body size is one of the main selection indices in chicken breeding. Although often investigated, knowledge of the underlying genetic mechanisms is incomplete. The aim of the current study was to identify genomic regions associated with body size differences between Asian Game and Asian Bantam type chickens. In this study, 94 and 107 chickens from 4 Asian Game and 5 Asian Bantam type breeds, respectively, were genotyped using the chicken 580K single nucleotide polymorphism (SNP) array. A genome-wide association study (GWAS) and principal component analyses (PCA) were performed to identify genomic regions associated with body size related-traits such as wing length, shank length, shank thickness, keel length, and body weight. Hierarchical clustering of genotype data showed a clear genetic difference between the investigated Asian Game and Asian Bantam chicken types. GWAS identified 16 genomic regions associated with wing length (2, FDR ≤ 0.018), shank thickness (6, FDR ≤ 0.008), keel length (5, FDR ≤ 0.023), and body weight (3, FDR ≤ 0.041). PCA showed that the first principal component (PC1) separated the 2 chicken types and significantly correlated with the measured body size related-traits (P ≤ 2.24e-40). SNPs contributing significantly to PC1 were subjected to a more detailed investigation. This analysis identified 11 regions potentially associated with differences in body size related-traits. A region on chromosome 4 (GGA4) (17.3-21.3 Mb) was detected in both analyses GWAS and PCA. This region harbors 60 genes. Among them are myotubularin 1 (MTM1) and secreted frizzled-related protein 2 (SFPR2) which can be considered as potential candidate genes for body size related-traits. Our results clearly show that the investigated Asian Game type chicken breeds are genetically different from the Asian Bantam breeds. A region on GGA4 between 17.3 and 21.3 Mb was identified which contributes to the phenotypic difference, though further validation of candidate genes is necessary.

Intermittent fasting induces chronic changes in the hepatic gene expression of Red Jungle Fowl (Gallus gallus)

Background

Intermittent fasting (IF), the implementation of fasting periods of at least 12 consecutive hours on a daily to weekly basis, has received a lot of attention in recent years for imparting the life-prolonging and health-promoting effects of caloric restriction with no or only moderate actual restriction of caloric intake. IF is also widely practiced in the rearing of broiler breeders, the parent stock of meat-type chickens, who require strict feed restriction regimens to prevent the serious health problems associated with their intense appetites. Although intermittent fasting has been extensively used in this context to reduce feed competition and its resulting stress, the potential of IF in chickens as an alternative and complementary model to rodents has received less investigation. In both mammals and birds, the liver is a key component of the metabolic response to IF, responding to variations in energy balance. Here we use a microarray analysis to examine the liver transcriptomics of wild-type Red Jungle Fowl chickens fed either ad libitum, chronically restricted to around 70% of ad libitum daily or intermittently fasted (IF) on a 2:1 (2 days fed, 1 day fasted) schedule without actual caloric restriction. As red junglefowl are ancestral to domestic chicken breeds, these data serve as a baseline to which existing and future transcriptomic results from farmed birds such as broiler breeders can be compared.

Results

We find large effects of feeding regimen on liver transcriptomics, with most of the affected genes relating to energy metabolism. A cluster analysis shows that IF is associated with large and reciprocal changes in genes related to lipid and carbohydrate metabolism, but also chronic changes in genes related to amino acid metabolism (generally down-regulated) and cell cycle progression (generally up-regulated). The overall transcription pattern appears to be one of promoting high proliferative plasticity in response to fluctuations in available energy substrates. A small number of inflammation-related genes also show chronically changed expression profiles, as does one circadian rhythm gene.

Conclusions

The increase in proliferative potential suggested by the gene expression changes reported here indicates that birds and mammals respond similarly to intermittent fasting practices. Our findings therefore suggest that the health benefits of periodic caloric restriction are ubiquitous and not restricted to mammals alone. Whether a common fundamental mechanism, for example involving leptin, underpins these benefits remains to be elucidated.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08533-5.

Jejunal transcriptomic profiling of two layer strains throughout the entire production period

The jejunum plays crucial roles for the digestion and absorption of nutrients and minerals and for barrier functions that are essential for a healthy, productive life cycle of farm animals, including laying hens. Accordingly, knowledge of the molecular pathways that emerge in the intestine during development, and particularly at the beginning of laying activity, will help to derive strategies for improving nutrient efficiency in laying hens. In this study, jejunal samples were obtained from two high-yielding layer strains at five developmental stages (weeks 10, 16, 24, 30 and 60 of life) for RNA-sequencing, alongside the profiling of blood plasma parameters to approximate the dynamics of mineral homeostasis. The results reflected a marked distinction between the pre-laying and laying phase as inferred from levels of parathyroid hormone, triiodothyronine, estradiol, vitamin D, and calcium. Moreover, the expression patterns of the intestinal mucosa responded directly to the changing metabolic and nutritional profiles at the beginning of the laying phase in maturing high-yielding strains of laying hens. These comprise signaling events namely RANK/RANKL signaling and cellular senescence. Taken together, the timing of sexual maturity of laying hens demands closer examination to unravel metabolic requirements and associated endogenous mechanisms.

Shifts in growth, but not differentiation, foreshadow the formation of exaggerated forms under chicken domestication

Domestication provides an outstanding opportunity for biologists to explore the underpinnings of organismal diversification. In domesticated animals, selective breeding for exaggerated traits is expected to override genetic correlations that normally modulate phenotypic variation in nature. Whether this strong directional selection affects the sequence of tightly synchronized events by which organisms arise (ontogeny) is often overlooked. To address this concern, we compared the ontogeny of the red junglefowl (RJF) (Gallus gallus) to four conspecific lineages that underwent selection for traits of economic or ornamental value to humans. Trait differentiation sequences in embryos of these chicken breeds generally resembled the representative ancestral condition in the RJF, thus revealing that early ontogeny remains highly canalized even during evolution under domestication. This key finding substantiates that the genetic cost of domestication does not necessarily compromise early ontogenetic steps that ensure the production of viable offspring. Instead, disproportionate beak and limb growth (allometry) towards the end of ontogeny better explained phenotypes linked to intense selection for industrial-scale production over the last 100 years. Illuminating the spatial and temporal specificity of development is foundational to the enhancement of chicken breeds, as well as to ongoing research on the origins of phenotypic variation in wild avian species.

Characterization of Alternative Splicing (AS) Events during Chicken (Gallus gallus) Male Germ-Line Stem Cell Differentiation with Single-Cell RNA-seq

Simple Summary

Studies have shown that alternative splicing (AS) has been utilized in a wide variety of life processes. However, there are very few studies on AS during germ cell development. In this study, we preliminarily investigated the variation of variable shear events during the formation of chicken germ cells through the RNA-seq data analysis of embryonic stem cells (ESCs), gonad PGCs (gPGCs), and spermatogonia stem cells (SSCs), and the critical AS mode for several crucial stage-specific genes, which were identified during germ cell development. The results of this study lay a theoretical foundation for further analysis of the regulation mechanism of key genes involved in germ cell formation.

Abstract

Alternative splicing (AS) is a ubiquitous, co-transcriptional, and post-transcriptional regulation mechanism during certain developmental processes, such as germ cell differentiation. A thorough understanding of germ cell differentiation will help us to open new avenues for avian reproduction, stem cell biology, and advances in medicines for human consumption. Here, based on single-cell RNA-seq, we characterized genome-wide AS events in manifold chicken male germ cells: embryonic stem cells (ESCs), gonad primordial germ cells (gPGCs), and spermatogonia stem cells (SSCs). A total of 38,494 AS events from 15,338 genes were detected in ESCs, with a total of 48,955 events from 14,783 genes and 49,900 events from 15,089 genes observed in gPGCs and SSCs, respectively. Moreover, this distribution of AS events suggests the diverse splicing feature of ESCs, gPGCs, and SSCs. Finally, several crucial stage-specific genes, such as NANOG, POU5F3, LIN28B, BMP4, STRA8, and LHX9, were identified in AS events that were transmitted in ESCs, gPGCs, and SSCs. The gene expression results of the RNA-seq data were validated by qRT-PCR. In summary, we provided a comprehensive atlas of the genome-wide scale of the AS event landscape in male chicken germ-line cells and presented its distribution for the first time. This research may someday improve treatment options for men suffering from male infertility.

The methylation landscape and its role in domestication and gene regulation in the chicken

Domestication is one of the strongest examples of artificial selection and has produced some of the most extreme within-species phenotypic variation known. In the case of the chicken, it has been hypothesized that DNA methylation may play a mechanistic role in the domestication response. By inter-crossing wild-derived red junglefowl with domestic chickens, we mapped quantitative trait loci for hypothalamic methylation (methQTL), gene expression (eQTL) and behaviour. We find large, stable methylation differences, with 6,179 cis and 2,973 trans methQTL identified. Over 46% of the trans effects were genotypically controlled by five loci, mainly associated with increased methylation in the junglefowl genotype. In a third of eQTL, we find that there is a correlation between gene expression and methylation, while statistical causality analysis reveals multiple instances where methylation is driving gene expression, as well as the reverse. We also show that methylation is correlated with some aspects of behavioural variation in the inter-cross. In conclusion, our data suggest a role for methylation in the regulation of gene expression underlying the domesticated phenotype of the chicken.

Transcriptome Sequencing and Comparative Analysis of Amphoteric ESCs and PGCs in Chicken (Gallus gallus)

Chicken (Gallus gallus) pluripotent embryonic stem cells (ESCs) and primordial germ cells (PGCs) can be broadly applied in the research of developmental and embryonic biology, but the difference between amphoteric ESCs and PGCs is still elusive. This study determined the sex of collected samples by identifying specific sex markers via polymerase chain reaction (PCR) and fluorescence activated cell sorter (FACS). RNA-seq was utilized to investigate the transcriptomic profile of amphoteric ESCs and PGCs in chicken. The results showed no significant differentially expressed genes (DEGs) in amphoteric ESCs and 227 DEGs exhibited in amphoteric PGCs. Moreover, those 227 DEGs were mainly enriched in 17 gene ontology (GO) terms and 27 pathways according to Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. Furthermore, qRT-PCR was performed to verify RNA-seq results, and the results demonstrated that Notch1 was highly expressed in male PGCs. In summary, our results provided a knowledge base of chicken amphoteric ESCs and PGCs, which is helpful for future research in relevant biological processes.

Combined approaches to reveal genes associated with litter size in Yunshang black goats

Intensive artificial selection has been imposed in Yunshang black goats, the first black specialist mutton goat breed in China, with a breeding object of improving reproductive performance, which has contributed to reshaping of the genome including the characterization of SNP, ROH and haplotype. However, variation in reproductive ability exists in the present population. A WGS was implemented in two subpopulations (polytocous group, PG, and monotocous group, MG) with evident differences of litter size. Following the mapping to reference genome, and SNP calling and pruning, three approaches - GWAS, ROH analysis and detection of signatures of selection - were employed to unveil candidate genes responsible for litter size. Consequently, 12 candidate genes containing OSBPL8 with the minimum P-value were uncovered by GWAS. Differences were observed in the pattern of ROH between two subpopulations that shared similar low inbreeding coefficients. Two ROH hotspots and 12 corresponding genes emerged from ROH pool association analysis. Based on the nS_L statistic, 15 and 61 promising genes were disclosed under selection for MG and PG respectively. Of them, some promising genes participate in ovarian function (PPP2R5C, CDC25A, ESR1, RPS26 and SERPINBs), seasonal reproduction (DIO3, BTG1 and CRYM) and metabolism (OSBPL8, SLC39A5 and SERPINBs). Our study pinpointed some novel promising genes influencing litter size, provided a comprehensive insight into genetic makeup of litter size and might facilitate selective breeding in goats.

The genetic regulation of size variation in the transcriptome of the cerebrum in the chicken and its role in domestication and brain size evolution

Background

Large difference in cerebrum size exist between avian species and populations of the same species and is believed to reflect differences in processing power, i.e. in the speed and efficiency of processing information in this brain region. During domestication chickens developed a larger cerebrum compared to their wild progenitor, the Red jungle fowl. The underlying mechanisms that control cerebrum size and the extent to which genetic regulation is similar across brain regions is not well understood. In this study, we combine measurement of cerebrum size with genome-wide genetical genomics analysis to identify the genetic architecture of the cerebrum, as well as compare the regulation of gene expression in this brain region with gene expression in other regions of the brain (the hypothalamus) and somatic tissue (liver).

Results

We identify one candidate gene that putatively regulates cerebrum size (MTF2) as well as a large number of eQTL that regulate the transcriptome in cerebrum tissue, with the majority of these eQTL being trans-acting. The overall regulation of gene expression variation in the cerebrum was markedly different to the hypothalamus, with relatively few eQTL in common. In comparison, the cerebrum tissue shared more eQTL with a distant tissue (liver) than with a neighboring tissue (hypothalamus).

Conclusion

The candidate gene for cerebrum size (MTF2) has previously been linked to brain development making it a good candidate for further investigation as a regulator of inter-population variation in cerebrum size. The lack of shared eQTL between the two brain regions implies that genetic regulation of gene expression appears to be relatively independent between the two brain regions and suggest that coevolution between these two brain regions might be more functionally driven than developmental. These findings have relevance for current brain size evolution theories.

Genome-Wide Association Study in Two Cohorts from a Multi-generational Mouse Advanced Intercross Line Highlights the Difficulty of Replication Due to Study-Specific Heterogeneity

There has been extensive discussion of the "Replication Crisis" in many fields, including genome-wide association studies (GWAS). We explored replication in a mouse model using an advanced intercross line (AIL), which is a multigenerational intercross between two inbred strains. We re-genotyped a previously published cohort of LG/J x SM/J AIL mice (F34; n = 428) using a denser marker set and genotyped a new cohort of AIL mice (F39-43; n = 600) for the first time. We identified 36 novel genome-wide significant loci in the F34 and 25 novel loci in the F39-43 cohort. The subset of traits that were measured in both cohorts (locomotor activity, body weight, and coat color) showed high genetic correlations, although the SNP heritabilities were slightly lower in the F39-43 cohort. For this subset of traits, we attempted to replicate loci identified in either F34 or F39-43 in the other cohort. Coat color was robustly replicated; locomotor activity and body weight were only partially replicated, which was inconsistent with our power simulations. We used a random effects model to show that the partial replications could not be explained by Winner's Curse but could be explained by study-specific heterogeneity. Despite this heterogeneity, we performed a mega-analysis by combining F34 and F39-43 cohorts (n = 1,028), which identified four novel loci associated with locomotor activity and body weight. These results illustrate that even with the high degree of genetic and environmental control possible in our experimental system, replication was hindered by study-specific heterogeneity, which has broad implications for ongoing concerns about reproducibility.

Analysis of the Progeny of Sibling Matings Reveals Regulatory Variation Impacting the Transcriptome of Immune Cells in Commercial Chickens

There is increasing recognition that the underlying genetic variation contributing to complex traits influences transcriptional regulation and can be detected at a population level as expression quantitative trait loci. At the level of an individual, allelic variation in transcriptional regulation of individual genes can be detected by measuring allele-specific expression in RNAseq data. We reasoned that extreme variants in gene expression could be identified by analysis of inbred progeny with shared grandparents. Commercial chickens have been intensively selected for production traits. Selection is associated with large blocks of linkage disequilibrium with considerable potential for co-selection of closely linked "hitch-hiker alleles" affecting traits unrelated to the feature being selected, such as immune function, with potential impact on the productivity and welfare of the animals. To test this hypothesis that there is extreme allelic variation in immune-associated genes we sequenced a founder population of commercial broiler and layer birds. These birds clearly segregated genetically based upon breed type. Each genome contained numerous candidate null mutations, protein-coding variants predicted to be deleterious and extensive non-coding polymorphism. We mated selected broiler-layer pairs then generated cohorts of F2 birds by sibling mating of the F1 generation. Despite the predicted prevalence of deleterious coding variation in the genomic sequence of the founders, clear detrimental impacts of inbreeding on survival and post-hatch development were detected in only one F2 sibship of 15. There was no effect on circulating leukocyte populations in hatchlings. In selected F2 sibships we performed RNAseq analysis of the spleen and isolated bone marrow-derived macrophages (with and without lipopolysaccharide stimulation). The results confirm the predicted emergence of very large differences in expression of individual genes and sets of genes. Network analysis of the results identified clusters of co-expressed genes that vary between individuals and suggested the existence of trans-acting variation in the expression in macrophages of the interferon response factor family that distinguishes the parental broiler and layer birds and influences the global response to lipopolysaccharide. This study shows that the impact of inbreeding on immune cell gene expression can be substantial at the transcriptional level, and potentially opens a route to accelerate selection using specific alleles known to be associated with desirable expression levels.

Genetical Genomics of Tonic Immobility in the Chicken

Identifying the molecular mechanisms of animal behaviour is an enduring goal for researchers. Gaining insight into these mechanisms enables us to gain a greater understanding of behaviour and their genetic control. In this paper, we perform Quantitative Trait Loci (QTL) mapping of tonic immobility behaviour in an advanced intercross line between wild and domestic chickens. Genes located within the QTL interval were further investigated using global expression QTL (eQTL) mapping from hypothalamus tissue, as well as causality analysis. This identified five candidate genes, with the genes PRDX4 and ACOT9 emerging as the best supported candidates. In addition, we also investigated the connection between tonic immobility, meat pH and struggling behaviour, as the two candidate genes PRDX4 and ACOT9 have previously been implicated in controlling muscle pH at slaughter. We did not find any phenotypic correlations between tonic immobility, struggling behaviour and muscle pH in a smaller additional cohort, despite these behaviours being repeatable within-test.

Changes in pituitary gene expression may underlie multiple domesticated traits in chickens

Domesticated animals share a unique set of morphological and behavioral traits, jointly referred to as the domesticated phenotype. Striking similarities amongst a range of unrelated domesticated species suggest that similar regulatory mechanisms may underlie the domesticated phenotype. These include color pattern, growth, reproduction, development and stress response. Although previous studies have focused on the brain to find mechanisms underlying domestication, the potential role of the pituitary gland as a target of domestication is highly overlooked. Here, we study gene expression in the pituitary gland of the domesticated White Leghorn chicken and its wild ancestor, the Red Junglefowl. By overlapping differentially expressed genes with a previously published list of functionally important genes in the pituitary gland, we narrowed down to 34 genes. Amongst them, expression levels of genes with inhibitory function on pigmentation (ASIP), main stimulators of metabolism and sexual maturity (TSHB and DIO2), and a potential inhibitor of broodiness (PRLR), were higher in the domesticated breed. Additionally, expression of 2 key inhibitors of the stress response (NR3C1, CRHR2) was higher in the domesticated breed. We suggest that changes in the transcription of important modulatory genes in the pituitary gland can account not only for domestication of the stress response in domestic chickens, but also for changes in pigmentation, development, and reproduction. Given the pivotal role of the pituitary gland in the regulation of multiple shared domesticated traits, we suggest that similar changes in pituitary transcriptome may contribute to the domesticated phenotype in other species as well.