Large-Scale Whole Genome Sequencing Study Reveals Genetic Architecture and Key Variants for Breast Muscle Weight in Native Chickens

A large structural variant collection in Holstein cattle and associated database for variant discovery, characterization, and application

BACKGROUND

Structural variants (SVs) such as deletions, duplications, and insertions are known to contribute to phenotypic variation but remain challenging to identify and genotype. A more complete, accessible, and assessable collection of SVs will assist efforts to study SV function in cattle and to incorporate SV genotyping into animal evaluation.

RESULTS

In this work we produced a large and deeply characterized collection of SVs in Holstein cattle using two popular SV callers (Manta and Smoove) and publicly available Illumina whole-genome sequence (WGS) read sets from 310 samples (290 male, 20 female, mean 20X coverage). Manta and Smoove identified 31 K and 68 K SVs, respectively. In total the SVs cover 5% (Manta) and 6% (Smoove) of the reference genome, in contrast to the 1% impacted by SNPs and indels. SV genotypes from each caller were confirmed to accurately recapitulate animal relationships estimated using WGS SNP genotypes from the same dataset, with Manta genotypes outperforming Smoove, and deletions outperforming duplications. To support efforts to link the SVs to phenotypic variation, overlapping and tag SNPs were identified for each SV, using genotype sets extracted from the WGS results corresponding to two bovine SNP chips (BovineSNP50 and BovineHD). 9% (Manta) and 11% (Smoove) of the SVs were found to have overlapping BovineHD panel SNPs, while 21% (Manta) and 9% (Smoove) have BovineHD panel tag SNPs. A custom interactive database ( https://svdb-dc.pslab.ca ) containing the identified sequence variants with extensive annotations, gene feature information, and BAM file content for all SVs was created to enable the evaluation and prioritization of SVs for further study. Illustrative examples involving the genes POPDC3, ORM1, G2E3, FANCI, TFB1M, FOXC2, N4BP2, GSTA3, and COPA show how this resource can be used to find well-supported genic SVs, determine SV breakpoints, design genotyping approaches, and identify processed pseudogenes masquerading as deletions.

CONCLUSIONS

The resources developed through this study can be used to explore sequence variation in Holstein cattle and to develop strategies for studying SVs of interest. The lack of overlapping and tag SNPs from commonly used SNP chips for most of the SVs suggests that other genotyping approaches will be needed (for example direct genotyping) to understand their potential contributions to phenotype. The included SV genotype assessments point to challenges in characterizing SVs, especially duplications, using short-read data and support ongoing efforts to better characterize cattle genomes through long-read sequencing. Lastly, the identification of previously known functional SVs and additional CDS-overlapping SVs supports the phenotypic relevance of this dataset.

New insights into the genetic loci related to egg weight and age at first egg traits in broiler breeder

Egg weight (EW) and age at first egg (AFE) are economically important traits in breeder chicken production. The genetic basis of these traits, however, is far from understood, especially for broiler breeders. In this study, genetic parameter estimation, genome-wide association analysis, meta-analysis, and selective sweep analysis were carried out to identify genetic loci associated with EW and AFE in 6,842 broiler breeders. The study found that the heritability of EW ranged from 0.42 to 0.44, while the heritability of AFE was estimated at 0.33 in the maternal line. Meta-analysis and selective sweep analysis identified two colocalized regions on GGA4 that significantly influenced EW at 32 wk (EW32W) and at 43 wk (EW43W) with both paternal and maternal lines. The genes AR, YIPF6, and STARD8 were located within the significant region (GGA4: 366.86-575.50 kb), potentially affecting EW through the regulation of follicle development, cell proliferation, and lipid transfer etc. The promising genes LCORL and NCAPG were positioned within the significant region (GGA4:75.35-75.42 Mb), potentially influencing EW through pleiotropic effects on growth and development. Additionally, 3 significant regions were associated with AFE on chromosomes GGA7, GGA19, and GGA27. All of these factors affected the AFE by influencing ovarian development. In our study, the genomic information from both paternal and maternal lines was used to identify genetic regions associated with EW and AFE. Two genomic regions and eight genes were identified as the most likely candidates affecting EW and AFE. These findings contribute to a better understanding of the genetic basis of egg production traits in broiler breeders and provide new insights into future technology development.

Integrating miRNA and full-length transcriptome profiling to elucidate the mechanism of muscle growth in Muscovy ducks reveals key roles for miR-301a-3p/ANKRD1

BACKGROUND

The popularity of Muscovy ducks is attributed not only to their conformation traits but also to their slightly higher content of breast and leg meat, as well as their stronger-tasting meat compared to that of typical domestic ducks. However, there is a lack of comprehensive systematic research on the development of breast muscle in Muscovy ducks. In addition, since the number of skeletal muscle myofibers is established during the embryonic period, this study conducted a full-length transcriptome sequencing and microRNA sequencing of the breast muscle. Muscovy ducks at four developmental stages, namely Embryonic Day 21 (E21), Embryonic Day 27 (E27), Hatching Day (D0), and Post-hatching Day 7 (D7), were used to isolate total RNA for analysis.

RESULTS

A total of 68,161 genes and 472 mature microRNAs were identified. In order to uncover deeper insights into the regulation of mRNA by miRNAs, we conducted an integration of the differentially expressed miRNAs (known as DEMs) with the differentially expressed genes (referred to as DEGs) across various developmental stages. This integration allowed us to make predictions regarding the interactions between miRNAs and mRNA. Through this analysis, we identified a total of 274 DEGs that may serve as potential targets for the 68 DEMs. In the predicted miRNA‒mRNA interaction networks, let-7b, miR-133a-3p, miR-301a-3p, and miR-338-3p were the hub miRNAs. In addition, multiple DEMs also showed predicted target relationships with the DEGs associated with skeletal system development. These identified DEGs and DEMs as well as their predicted interaction networks involved in the regulation of energy homeostasis and muscle development were most likely to play critical roles in facilitating the embryo-to-hatchling transition. A candidate miRNA, miR-301a-3p, exhibited increased expression during the differentiation of satellite cells and was downregulated in the breast muscle tissues of Muscovy ducks at E21 compared to E27. A dual-luciferase reporter assay suggested that the ANKRD1 gene, which encodes a transcription factor, is a direct target of miR-301a-3p.

CONCLUSIONS

miR-301a-3p suppressed the posttranscriptional activity of ANKRD1, which is an activator of satellite cell proliferation, as determined with gain- and loss-of-function experiments. miR-301a-3p functions as an inducer of myogenesis by targeting the ANKRD1 gene in Muscovy ducks. These results provide novel insights into the early developmental process of black Muscovy breast muscles and will improve understanding of the underlying molecular mechanisms.

Whole-genome variants dataset of 209 local chickens from China

Compared to commercial chickens, local breeds exhibit better in meat quality and flavour, but the productivity (e.g., growth rate, body weight) of local chicken breeds is rather low. Genetic analysis based on whole-genome sequencing contributes to elucidating the genetic markers or putative candidate genes related to some economic traits, facilitating the improvement of production performance, the acceleration of breeding progress, and the conservation of genetic resources. Here, a total of 209 local chickens from 13 breeds were investigated, and the observation of approximately 91.4% high-quality sequences (Q30 > 90%) and a mapping rate over 99% for each individual indicated good results of this study, as confirmed by a genome coverage of 97.6%. Over 19 million single nucleotide polymorphisms (SNPs) and 1.98 million insertion-deletions (InDels) were identified using the reference genome (GRCg7b), further contributing to the public database. This dataset provides valuable resources for studying genetic diversity and adaptation and for the cultivation of new chicken breeds/lines.

Large-scale genomic and transcriptomic analyses elucidate the genetic basis of high meat yield in chickens

INTRODUCTION

Investigating the genetic markers and genomic signatures related to chicken meat production by combing multi-omics methods could provide new insights into modern chicken breeding technology systems.

OBJECT

Chicken is one of the most efficient and environmentally friendly livestock, especially the fast-growing white-feathered chicken (broiler), which is well known for high meat yield, but the underlying genetic basis is poorly understood.

METHOD

We generated whole-genome resequencing of three purebred broilers (n = 748) and six local breeds/lines (n = 114), and sequencing data of twelve chicken breeds (n = 199) were obtained from the NCBI database. Additionally, transcriptome sequencing of six tissues from two chicken breeds (n = 129) at two developmental stages was performed. A genome-wide association study combined with cis-eQTL mapping and the Mendelian randomization was applied.

RESULT

We identified > 17 million high-quality SNPs, of which 21.74% were newly identified, based on 21 chicken breeds/lines. A total of 163 protein-coding genes underwent positive selection in purebred broilers, and 83 genes were differentially expressed between purebred broilers and local chickens. Notably, muscle development was proven to be the major difference between purebred broilers and local chickens, or ancestors, based on genomic and transcriptomic evidence from multiple tissues and stages. The MYH1 gene family showed the top selection signatures and muscle-specific expression in purebred broilers. Furthermore, we found that the causal gene SOX6 influenced breast muscle yield and also related to myopathy occurrences. A refined haplotype was provided, which had a significant effect on SOX6 expression and phenotypic changes.

CONCLUSION

Our study provides a comprehensive atlas comprising the typical genomic variants and transcriptional characteristics for muscle development and suggests a new regulatory target (SOX6-MYH1s axis) for breast muscle yield and myopathy, which could aid in the development of genome-scale selective breeding aimed at high meat yield in broiler chickens.

Genetic Architecture of Abdominal Fat Deposition Revealed by a Genome-Wide Association Study in the Laying Chicken

Fat has a high energy density, and excessive fatness has been recognized as a problem for egg production and the welfare of chickens. The identification of a genetic polymorphism controlling fat deposition would be helpful to select against excessive fatness in the laying hen. This study aimed to estimate genomic heritability and identify the genetic architecture of abdominal fat deposition in a population of chickens from a Dongxiang blue-shelled local breed crossbred with the White Leghorn. A genome-wide association study was conducted on abdominal fat percentage, egg production and body weights using a sample of 1534 hens genotyped with a 600 K Chicken Genotyping Array. The analysis yielded a heritability estimate of 0.19 ± 0.04 for abdominal fat percentage; 0.56 ± 0.04 for body weight at 72 weeks; 0.11 ± 0.03 for egg production; and 0.24 ± 0.04 for body weight gain. The genetic correlation of abdominal fat percentage with egg production between 60 and 72 weeks of age was -0.35 ± 0.18. This implies a potential trade-off between these two traits related to the allocation of resources. Strong positive genetic correlations were found between fat deposition and weight traits. A promising locus close to COL12A1 on chromosome 3, associated with abdominal fat percent, was found in the present study. Another region located around HTR2A on chromosome 1, where allele substitution was predicted to be associated with body weight gain, accounted for 2.9% of phenotypic variance. Another region located on chromosome 1, but close to SOX5, was associated with egg production. These results may be used to influence the balanced genetic selection for laying hens.

Multi-Omics Analysis of Genes Encoding Proteins Involved in Alpha-Linolenic Acid Metabolism in Chicken

Alpha-linolenic acid (ALA, ω-3) is an antioxidant that reduces triglyceride (TG) levels in blood, a component of cell membranes and a precursor compound of eicosapentaenoic acid (EPA, ω-3) and eicosatrienoic acid (DHA, ω-3). Fatty acid content is a quantitative trait regulated by multiple genes, and the key genes regulating fatty acid metabolism have not been systematically identified. This study aims at investigating the protein-encoding genes regulating ω-3 polyunsaturated fatty acid (PUFA) content in chicken meat. We integrated genomics, transcriptomics and lipidomics data of Jingxing yellow chicken (JXY) to explore the interactions and associations among multiple genes involved in the regulation of fatty acid metabolism. Several key genes and pathways regulating ω-3 fatty acid metabolism in chickens were identified. The upregulation of GRB10 inhibited the mTOR signaling pathway, thereby improving the content of EPA and DHA. The downregulation of FGFR3 facilitated the conversion of ALA to EPA. Additionally, we analyzed the effects of ALA supplementation dose on glycerol esters (GLs), phospholipid (PL) and fatty acyl (FA) contents, as well as the regulatory mechanisms of nutritional responses in FFA metabolism. This study provides a basis for identifying genes and pathways that regulate the content of FFAs, and offers a reference for nutritional regulation systems in production.

RNA-Seq-based discovery of genetic variants and allele-specific expression of two layer lines and broiler chicken

Recent advances in the selective breeding of broilers and layers have made poultry production one of the fastest-growing industries. In this study, a transcriptome variant calling approach from RNA-seq data was used to determine population diversity between broilers and layers. In total, 200 individuals were analyzed from three different chicken populations (Lohmann Brown (LB), n = 90), Lohmann Selected Leghorn (LSL, n = 89), and Broiler (BR, n = 21). The raw RNA-sequencing reads were pre-processed, quality control checked, mapped to the reference genome, and made compatible with Genome Analysis ToolKit for variant detection. Subsequently, pairwise fixation index (F ST) analysis was performed between broilers and layers. Numerous candidate genes were identified, that were associated with growth, development, metabolism, immunity, and other economically significant traits. Finally, allele-specific expression (ASE) analysis was performed in the gut mucosa of LB and LSL strains at 10, 16, 24, 30, and 60 weeks of age. At different ages, the two-layer strains showed significantly different allele-specific expressions in the gut mucosa, and changes in allelic imbalance were observed across the entire lifespan. Most ASE genes are involved in energy metabolism, including sirtuin signaling pathways, oxidative phosphorylation, and mitochondrial dysfunction. A high number of ASE genes were found during the peak of laying, which were particularly enriched in cholesterol biosynthesis. These findings indicate that genetic architecture as well as biological processes driving particular demands relate to metabolic and nutritional requirements during the laying period shape allelic heterogeneity. These processes are considerably affected by breeding and management, whereby elucidating allele-specific gene regulation is an essential step towards deciphering the genotype to phenotype map or functional diversity between the chicken populations. Additionally, we observed that several genes showing significant allelic imbalance also colocalized with the top 1% of genes identified by the FST approach, suggesting a fixation of genes in cis-regulatory elements.

Molecular Characterization, Expression Profile, and A 21-bp Indel within the ASB9 Gene and Its Associations with Chicken Production Traits

A growing number of studies have shown that members of the ankyrin repeat and suppressors of cytokine signaling (SOCS) box-containing protein (ASB) family are extensively involved in biological processes such as cell growth, tissue development, insulin signaling, ubiquitination, protein degradation, and skeletal muscle membrane protein formation, while the specific biological role of ankyrin-repeat and SOCS box protein 9 (ASB9) remains unclear. In this study, a 21 bp indel in the intron of ASB9 was identified for the first time in 2641 individuals from 11 different breeds and an F₂ resource population, and differences were observed among individuals with different genotypes (II, ID, and DD). An association study of a cross-designed F₂ resource population revealed that the 21-bp indel was significantly related to growth and carcass traits. The significantly associated growth traits were body weight (BW) at 4, 6, 8, 10, and 12 weeks of age; sternal length (SL) at 4, 8, and 12 weeks of age; body slope length (BSL) at 4, 8, and 12 weeks of age; shank girth (SG) at 4 and 12 weeks of age; tibia length (TL) at 12 weeks of age; and pelvic width (PW) at 4 weeks of age (p < 0.05). This indel was also significantly correlated with carcass traits including semievisceration weight (SEW), evisceration weight (EW), claw weight (CLW), breast muscle weight (BMW), leg weight (LeW), leg muscle weight (LMW), claw rate (CLR), and shedding weight (ShW) (p < 0.05). In commercial broilers, the II genotype was the dominant genotype and underwent extensive selection. Interestingly, the ASB9 gene was expressed at significantly higher levels in the leg muscles of Arbor Acres broilers than those of Lushi chickens, while the opposite was true for the breast muscles. In summary, the 21-bp indel in the ASB9 gene significantly influenced the expression of the ASB9 gene in muscle tissue and was associated with multiple growth and carcass traits in the F₂ resource population. These findings suggested that the 21-bp indel within the ASB9 gene could be used in marker-assisted selection breeding for traits related to chicken growth.

Identification of New Genes and Genetic Variant Loci Associated with Breast Muscle Development in the Mini-Cobb F2 Chicken Population Using a Genome-Wide Association Study

Native chicken has become a favorite choice for consumers in many Asian countries recently, not only for its potential nutritional value but also for its deep ties to local food culture. However, low growth performance and limited meat production restrict their economic potential. Conducting a genome-wide association study (GWAS) for chicken-breast muscle development will help identify loci or candidate genes for different traits and potentially provide new insight into this phenotype in chickens and other species. To improve native chicken growth performance, especially breast muscle development, we performed a GWAS to explore the potential genetic mechanisms of breast muscle development in an F2 population constructed by reciprocal crosses between a fast-growing broiler chicken (Cobb500) and a slow-growing native chicken (Daweishan mini chicken). The results showed that 11 SNPs, which exceeded the 10% genome significance level (p = 1.79 × 10-8) were considered associated with breast muscle development traits, where six SNPS, NC_006126.5: g.3138376T>G, NC_006126.5: g.3138452A>G, NC_006088.5: g.73837197A>G, NC_006088.5: g.159574275A>G, NC_006089.5: g.80832197A>G, and NC_006127.5: g.48759869G>T was first identified in this study. In total, 13 genes near the SNPs were chosen as candidate genes, and none of them had previously been studied for their role in breast muscle development. After grouping the F2 population according to partial SNPs, significant differences in breast muscle weight were found among different genotypes (p < 0.05), and the expression levels of ALOX5AP, USPL1, CHRNA9, and EFNA5 among candidate genes were also significantly different (p < 0.05). The results of this study will contribute to the future exploration of the potential genetic mechanisms of breast muscle development in domestic chickens and also support the expansion of the market for native chicken in the world.

Fatty acid metabolism-related genes are associated with flavor-presenting aldehydes in Chinese local chicken

Aldehydes are primary volatile organic compounds (VOCs) in local Chinese chicken meat and contribute green grass, fatty, citrus, and bitter almond aromas to chicken meat. To understand the genetic basis of these aldehyde VOC aromas, we used approximately 500 Chinese Jingxing Yellow (JXY) chickens to conduct genome-wide association studies (GWAS) on the flavor traits with the data of single nucleotide polymorphisms (SNPs) and insertions and deletions (INDELs). In total, 501 association variants (253 SNPs and 248 INDELs) were found to be suggestively (SNPs: p-value < 2.77e-06 and INDELs: p-value < 3.78e-05) associated with total aldehydes (the sum of nine aldehydes), hexanal, heptanal, benzaldehyde, (E,E)-2,4-nonadienal, octanal, (E)-2-decenal, nonanal, decanal, and octadecanal. Of them, six SNPs and 23 INDELs reached a genome-wide significance level (SNPs: p-value < 1.38e-07 and INDELs: p-value < 1.89e-06). Potential candidate aldehyde genes were functionally annotated for lipid metabolism, especially fatty acid-related pathways and phospholipid-related gene ontology (GO) terms. Moreover, the GWAS analysis of total aldehydes, hexanal, and nonanal generated the most significant signals, and phenotypic content differed between different genotypes at candidate gene-related loci. For total aldehydes and hexanal traits, candidate genes were annotated based on the significant and suggestive variants on chromosomes 3 and 8 with highly polymorphic linkage blocks. The following candidate genes were also identified: GALM, MAP4K3, GPCPD1, RPS6KA2, CRLS1, ASAP1, TRMT6, SDC1, PUM2, ALDH9A1, MGST3, GMEB1, MECR, LDLRAP1, GPAM and ACSL5. We also found that polyunsaturated fatty acids (PUFAs) (C18:2n6c linoleic acid and C18:3n3 linolenic acid) were significantly correlated with total aldehydes and hexanal contents. PUFAs are important aldehyde precursors, and consistently, our results suggested that candidate genes involved in fatty acid pathways and phospholipid GO terms were identified in association loci. This work provides an understanding of the genetic basis of aldehyde formation, which is a key flavor-forming compound.