Enhanced bovine genome annotation through integration of transcriptomics and epi-transcriptomics datasets facilitates genomic biology.

BACKGROUND: The accurate identification of the functional elements in the bovine genome is a fundamental requirement for high-quality analysis of data informing both genome biology and genomic selection. Functional annotation of the bovine genome was performed to identify a more complete catalog of transcript isoforms across bovine tissues. RESULTS: A total of 160,820 unique transcripts (50% protein coding) representing 34,882 unique genes (60% protein coding) were identified across tissues. Among them, 118,563 transcripts (73% of the total) were structurally validated by independent datasets (PacBio isoform sequencing data, Oxford Nanopore Technologies sequencing data, de novo assembled transcripts from RNA sequencing data) and comparison with Ensembl and NCBI gene sets. In addition, all transcripts were supported by extensive data from different technologies such as whole transcriptome termini site sequencing, RNA Annotation and Mapping of Promoters for the Analysis of Gene Expression, chromatin immunoprecipitation sequencing, and assay for transposase-accessible chromatin using sequencing. A large proportion of identified transcripts (69%) were unannotated, of which 86% were produced by annotated genes and 14% by unannotated genes. A median of two 5' untranslated regions were expressed per gene. Around 50% of protein-coding genes in each tissue were bifunctional and transcribed both coding and noncoding isoforms. Furthermore, we identified 3,744 genes that functioned as noncoding genes in fetal tissues but as protein-coding genes in adult tissues. Our new bovine genome annotation extended more than 11,000 annotated gene borders compared to Ensembl or NCBI annotations. The resulting bovine transcriptome was integrated with publicly available quantitative trait loci data to study tissue-tissue interconnection involved in different traits and construct the first bovine trait similarity network. CONCLUSIONS: These validated results show significant improvement over current bovine genome annotations.

A combinatorial approach implementing new database structures to facilitate practical data curation management of QTL, association, correlation and heritability data on trait variants.

A precise description of traits is essential in genetics and genomics studies to facilitate comparative genetics and meta-analyses. It is an ongoing challenge in research and production environments to unambiguously and consistently compare traits of interest from data collected under various conditions. Despite previous efforts to standardize trait nomenclature, it remains a challenge to fully and accurately capture trait nomenclature granularity in a way that ensures long-term data sustainability in terms of the data curation processes, data management logistics and the ability to make meaningful comparisons across studies. In the Animal Quantitative Trait Loci Database and the Animal Trait Correlation Database, we have recently introduced a new method to extend livestock trait ontologies by using trait modifiers and qualifiers to define traits that differ slightly in how they are measured, examined or combined with other traits or factors. Here, we describe the implementation of a system in which the extended trait data, with modifiers, are managed at the experiment level as 'trait variants'. This has helped us to streamline the management and curation of such trait information in our database environment. Database URL  https://www.animalgenome.org/PGNET/.

Transcriptome profile of skeletal muscle using different sources of dietary fatty acids in male pigs.

Pork is of great importance in world trade and represents the largest source of fatty acids in the human diet. Lipid sources such as soybean oil (SOY), canola (CO), and fish oil (FO) are used in pig diets and influence blood parameters and the ratio of deposited fatty acids. In this study, the main objective was to evaluate changes in gene expression in porcine skeletal muscle tissue resulting from the dietary oil sources and to identify metabolic pathways and biological process networks through RNA-Seq. The addition of FO in the diet of pigs led to intramuscular lipid with a higher FA profile composition of C20:5 n-3, C22:6 n-3, and SFA (C16:0 and C18:0). Blood parameters for the FO group showed lower cholesterol and HDL content compared with CO and SOY groups. Skeletal muscle transcriptome analyses revealed 65 differentially expressed genes (DEG, FDR 10%) between CO vs SOY, and 32 DEG for CO vs FO, and 531 DEG for SOY vs FO comparison. Several genes, including AZGP1, PDE3B, APOE, PLIN1, and LIPS, were found to be down-regulated in the diet of the SOY group compared to the FO group. The enrichment analysis revealed DEG involved in lipid metabolism, metabolic diseases, and inflammation between the oil groups, with specific gene functions in each group and altered blood parameters. The results provide mechanisms to help us understand the behavior of genes according to fatty acids.

Brain fatty acid and transcriptome profiles of pig fed diets with different levels of soybean oil.

BACKGROUND: The high similarity in anatomical and neurophysiological processes between pigs and humans make pigs an excellent model for metabolic diseases and neurological disorders. Lipids are essential for brain structure and function, and the polyunsaturated fatty acids (PUFA) have anti-inflammatory and positive effects against cognitive dysfunction in neurodegenerative diseases. Nutrigenomics studies involving pigs and fatty acids (FA) may help us in better understanding important biological processes. In this study, the main goal was to evaluate the effect of different levels of dietary soybean oil on the lipid profile and transcriptome in pigs' brain tissue. RESULTS: Thirty-six male Large White pigs were used in a 98-day study using two experimental diets corn-soybean meal diet containing 1.5% soybean oil (SOY1.5) and corn-soybean meal diet containing 3.0% soybean oil (SOY3.0). No differences were found for the brain total lipid content and FA profile between the different levels of soybean oil. For differential expression analysis, using the DESeq2 statistical package, a total of 34 differentially expressed genes (DEG, FDR-corrected p-value < 0.05) were identified. Of these 34 DEG, 25 are known-genes, of which 11 were up-regulated (log2 fold change ranging from + 0.25 to + 2.93) and 14 were down-regulated (log2 fold change ranging from - 3.43 to -0.36) for the SOY1.5 group compared to SOY3.0. For the functional enrichment analysis performed using MetaCore with the 34 DEG, four pathway maps were identified (p-value < 0.05), related to the ALOX15B (log2 fold change - 1.489), CALB1 (log2 fold change - 3.431) and CAST (log2 fold change + 0.421) genes. A "calcium transport" network (p-value = 2.303e-2), related to the CAST and CALB1 genes, was also identified. CONCLUSION: The results found in this study contribute to understanding the pathways and networks associated with processes involved in intracellular calcium, lipid metabolism, and oxidative processes in the brain tissue. Moreover, these results may help a better comprehension of the modulating effects of soybean oil and its FA composition on processes and diseases affecting the brain tissue.

RNA-seq transcriptome profiling of pigs' liver in response to diet with different sources of fatty acids.

Pigs (Sus scrofa) are an animal model for metabolic diseases in humans. Pork is an important source of fatty acids (FAs) in the human diet, as it is one of the most consumed meats worldwide. The effects of dietary inclusion of oils such as canola, fish, and soybean oils on pig gene expression are mostly unknown. Our objective was to evaluate FA composition, identify changes in gene expression in the liver of male pigs fed diets enriched with different FA profiles, and identify impacted metabolic pathways and gene networks to enlighten the biological mechanisms' variation. Large White male pigs were randomly allocated to one of three diets with 18 pigs in each; all diets comprised a base of corn and soybean meal to which either 3% of soybean oil (SOY), 3% canola oil (CO), or 3% fish oil (FO) was added for a 98-day trial during the growing and finishing phases. RNA sequencing was performed on the liver samples of each animal by Illumina technology for differential gene expression analyses, using the R package DESeq2. The diets modified the FA profile, mainly in relation to polyunsaturated and saturated FAs. Comparing SOY vs. FO, 143 differentially expressed genes (DEGs) were identified as being associated with metabolism, metabolic and neurodegenerative disease pathways, inflammatory processes, and immune response networks. Comparing CO vs. SOY, 148 DEGs were identified, with pathways related to FA oxidation, regulation of lipid metabolism, and metabolic and neurodegenerative diseases. Our results help explain the behavior of genes with differential expression in metabolic pathways resulting from feeding different types of oils in pig diets.

Whole-genome SNP markers reveal runs of homozygosity in indigenous cattle breeds of Pakistan.

The runs of homozygosity (ROH) were identified in 14 Pakistani cattle breeds (n = 105) by genotyping with the Illumina 50 K SNP BeadChip. These breeds were categorized into Dairy, Dual, and Draft breeds based on their utility and production performance. We identified a total of 10,936 ROHs which mainly consisted of a high number of shorter segments (1-4 Mb). Dairy group exhibited the highest level of inbreeding (FROH: 0.078 ± 0.028) while the lowest (FROH: 0.002 ± 0.008) was observed in Dual group. In 48 genomic regions identified with a high frequency of ROH, 207 genes were detected in the three breed groups. A substantially higher number of ROH islands detected in dairy breeds indicated the impact of the positive selection pressure over the years. Important candidate genes and QTL were detected in the ROH islands associated with economic traits like milk production, reproduction, meat, carcass, and health traits in dairy cattle.

Differential Gene Expression Associated with Soybean Oil Level in the Diet of Pigs.

The aim of this study was to identify the differentially expressed genes (DEG) from the skeletal muscle and liver samples of animal models for metabolic diseases in humans. To perform the study, the fatty acid (FA) profile and RNA sequencing (RNA-Seq) data of 35 samples of liver tissue (SOY1.5, n = 17 and SOY3.0, n = 18) and 36 samples of skeletal muscle (SOY1.5, n = 18 and SOY3.0, n = 18) of Large White pigs were analyzed. The FA profile of the tissues was modified by the diet, mainly those related to monounsaturated (MUFA) and polyunsaturated (PUFA) FA. The skeletal muscle transcriptome analysis revealed 45 DEG (FDR 10%), and the functional enrichment analysis identified network maps related to inflammation, immune processes, and pathways associated with oxidative stress, type 2 diabetes, and metabolic dysfunction. For the liver tissue, the transcriptome profile analysis revealed 281 DEG, which participate in network maps related to neurodegenerative diseases. With this nutrigenomics study, we verified that different levels of soybean oil in the pig diet, an animal model for metabolic diseases in humans, affected the transcriptome profile of skeletal muscle and liver tissue. These findings may help to better understand the biological mechanisms that can be modulated by the diet.

Prune homolog 2 with BCH domain (PRUNE2) gene expression is associated with feed efficiency-related traits in Nelore steers.

Animal feeding is a critical factor in increasing producer profitability. Improving feed efficiency can help reduce feeding costs and reduce the environmental impact of beef production. Candidate genes previously identified for this trait in differential gene expression studies (e.g., case-control studies) have not examined continuous gene-phenotype variation, which is a limitation. The aim of this study was to investigate the association between the expression of five candidate genes in the liver, measured by quantitative real-time PCR and feed-related traits. We adopted a linear mixed model to associate liver gene expression from 52 Nelore steers with the following production traits: average daily gain (ADG), body weight (BW), dry matter intake (DMI), feed conversion ratio (FCR), feed efficiency (FE), Kleiber index (KI), metabolic body weight (MBW), residual feed intake (RFI), and relative growth ratio (RGR). The total expression of the prune homolog 2 (PRUNE2) gene was significantly associated with DMI, FCR, FE, and RFI (P < 0.05). Furthermore, we have identified a new transcript of PRUNE2 (TCONS_00027692, GenBank MZ041267) that was inversely correlated with FCR and FE (P < 0.05), in contrast to the originally identified PRUNE2 transcript. The cytochrome P450 subfamily 2B (CYP2B6), early growth response protein 1 (EGR1), collagen type I alpha 1 chain (COL1A1), and connective tissue growth factor (CTGF) genes were not associated with any feed efficiency-related traits (P > 0.05). The findings reported herein suggest that PRUNE2 expression levels affects feed efficiency-related traits variation in Nelore steers.

Stool and Ruminal Microbiome Components Associated With Methane Emission and Feed Efficiency in Nelore Beef Cattle.

Background: The impact of extreme changes in weather patterns on the economy and human welfare is one of the biggest challenges our civilization faces. From anthropogenic contributions to climate change, reducing the impact of farming activities is a priority since it is responsible for up to 18% of global greenhouse gas emissions. To this end, we tested whether ruminal and stool microbiome components could be used as biomarkers for methane emission and feed efficiency in bovine by studying 52 Brazilian Nelore bulls belonging to two feed intervention treatment groups, that is, conventional and by-product-based diets. Results: We identified a total of 5,693 amplicon sequence variants (ASVs) in the Nelore bulls' microbiomes. A Differential abundance analysis with the ANCOM approach identified 30 bacterial and 15 archaeal ASVs as differentially abundant (DA) among treatment groups. An association analysis using Maaslin2 software and a linear mixed model indicated that bacterial ASVs are linked to the host's residual methane emission (RCH4) and residual feed intake (RFI) phenotype variation, suggesting their potential as targets for interventions or biomarkers. Conclusion: The feed composition induced significant differences in both abundance and richness of ruminal and stool microbial populations in ruminants of the Nelore breed. The industrial by-product-based dietary treatment applied to our experimental groups influenced the microbiome diversity of bacteria and archaea but not of protozoa. ASVs were associated with RCH4 emission and RFI in ruminal and stool microbiomes. While ruminal ASVs were expected to influence CH4 emission and RFI, the relationship of stool taxa, such as Alistipes and Rikenellaceae (gut group RC9), with these traits was not reported before and might be associated with host health due to their link to anti-inflammatory compounds. Overall, the ASVs associated here have the potential to be used as biomarkers for these complex phenotypes.

Effect of dietary soybean oil inclusion on liver-related transcription factors in a pig model for metabolic diseases.

Dietary fatty acids (FA) are components of the lipids, which contribute to membrane structure, energy input, and biological functions related to cellular signaling and transcriptome regulation. However, the consumers still associate dietary FA with fat deposition and increased occurrence of metabolic diseases such as obesity and atherosclerosis. Previous studies already demonstrated that some fatty acids are linked with inflammatory response, preventing metabolic diseases. To better understand the role of dietary FA on metabolic diseases, for the first time, a study to identify key transcription factors (TF) involved in lipid metabolism and inflammatory response by transcriptome analysis from liver samples of animal models was performed. The key TF were identified by functional enrichment analysis from the list of differentially expressed genes identified in liver samples between 35 pigs fed with 1.5% or 3.0% soybean oil. The functional enrichment analysis detected TF linked to lipid homeostasis and inflammatory response, such as RXRA, EGFR, and SREBP2 precursor. These findings demonstrated that key TF related to lipid metabolism could be modulated by dietary inclusion of soybean oil. It could contribute to nutrigenomics research field that aims to elucidate dietary interventions in animal and human health, as well as to drive food technology and science.

Variations in antibody repertoires correlate with vaccine responses.

An important challenge in vaccine development is to figure out why a vaccine succeeds in some individuals and fails in others. Although antibody repertoires hold the key to answering this question, there have been very few personalized immunogenomics studies so far aimed at revealing how variations in immunoglobulin genes affect a vaccine response. We conducted an immunosequencing study of 204 calves vaccinated against bovine respiratory disease (BRD) with the goal to reveal variations in immunoglobulin genes and somatic hypermutations that impact the efficacy of vaccine response. Our study represents the largest longitudinal personalized immunogenomics study reported to date across all species, including humans. To analyze the generated data set, we developed an algorithm for identifying variations of the immunoglobulin genes (as well as frequent somatic hypermutations) that affect various features of the antibody repertoire and titers of neutralizing antibodies. In contrast to relatively short human antibodies, cattle have a large fraction of ultralong antibodies that have opened new therapeutic opportunities. Our study reveals that ultralong antibodies are a key component of the immune response against the costliest disease of beef cattle in North America. The detected variants of the cattle immunoglobulin genes, which are implicated in the success/failure of the BRD vaccine, have the potential to direct the selection of individual cattle for ongoing breeding programs.

Bringing the Animal QTLdb and CorrDB into the future: meeting new challenges and providing updated services.

The Animal QTLdb (https://www.animalgenome.org/QTLdb) and CorrDB (https://www.animalgenome.org/CorrDB) are unique resources for livestock animal genetics and genomics research which have been used extensively by the international livestock genome research community. This is largely due to the active development of the databases over the years to keep up with the rapid advancement of genome sciences. The ongoing development has ensured that these databases provide researchers not only with continually updated data but also with new web tools to disseminate the data. Through our continued efforts, the databases have evolved from the original Pig QTLdb for cross-experiment QTL data comparisons to an Animal QTLdb hosting 220 401 QTL, SNP association and eQTL data linking phenotype to genotype for 2210 traits. In addition, there are 23 552 correlations for 866 traits and 4273 heritability data on 1069 traits in CorrDB. All these data were curated from 3157 publications that cover seven livestock species. Along with the continued data curation, new species, additional genome builds, and new functions and features have been built into the databases as well. Standardized procedures to support data mapping on multiple species/genome builds and the ability to browse data based on linked ontology terms are highlights of the recent developments.

Genome-Wide Association Study for Fatty Acid Composition in American Angus Cattle.

Livestock is an important commodity playing a major role in the global economy. Red meat plays an important role in human life, as it is a good source of animal protein and energy. The fatty acid content of beef has been shown to impact the eating experience and nutritional value of beef. Therefore, this study aimed to identify genomic regions which can account for genetic variation in meat fatty acid content. Genotypes imputed to the Illumina BovineHD 770K BeadChip were used in this study. Thirty-six 1-Mb genomic regions with a posterior probability of inclusion (PPI) greater than 0.90 were identified to be associated with variation in the content of at least one fatty acid. The genomic regions (1Mb) which were associated with more than one fatty acid trait with high genetic variance and harbored good candidate genes were on Chromosome (Chr) 6 (fatty acid binding protein 2), Chr 19 (thyroid hormone receptor alpha, fatty acid synthase), Chr 26 (stearoyl-CoA desaturase), and Chr 29 (thyroid hormone responsive, fatty acid desaturase 2, and fatty acid desaturase 3). Further studies are required to identify the causal variants within the identified genomic regions. Findings from the present study will help to increase understanding of the variation in fatty acid content of beef and help to enhance selection for beef with improved fatty acid composition.

Comparative untargeted metabolome analysis of ruminal fluid and feces of Nelore steers (Bos indicus).

We conducted a study to identify the fecal metabolite profile and its proximity to the ruminal metabolism of Nelore steers based on an untargeted metabolomic approach. Twenty-six Nelore were feedlot with same diet during 105 d. Feces and rumen fluid were collected before and at slaughter, respectively. The metabolomics analysis indicated 49 common polar metabolites in the rumen and feces. Acetate, propionate, and butyrate were the most abundant polar metabolites in both bio-samples. The rumen presented significantly higher concentrations of the polar compounds when compared to feces (P < 0.05); even though, fecal metabolites presented an accentuated representability of the ruminal fluid metabolites. All fatty acids present in the ruminal fluid were also observed in the feces, except for C20:2n6 and C20:4n6. The identified metabolites offer information on the main metabolic pathways (higher impact factor and P < 0.05), as synthesis and degradation of ketone bodies; the alanine, aspartate and glutamate metabolisms, the glycine, serine; and threonine metabolism and the pyruvate metabolism. The findings reported herein on the close relationship between the ruminal fluid and feces metabolic profiles may offer new metabolic information, in addition to facilitating the sampling for metabolism investigation in animal production and health routines.

Gene expression in tonsils in swine following infection with porcine reproductive and respiratory syndrome virus.

BACKGROUND: Porcine reproductive and respiratory syndrome (PRRS) is a threat to pig production worldwide. Our objective was to understand mechanisms of persistence of PRRS virus (PRRSV) in tonsil. Transcriptome data from tonsil samples collected at 42 days post infection (dpi) were generated by RNA-seq and NanoString on 51 pigs that were selected to contrast the two PRRSV isolates used, NVSL and KS06, high and low tonsil viral level at 42 dpi, and the favorable and unfavorable genotypes at a genetic marker (WUR) for the putative PRRSV resistance gene GBP5. RESULTS: The number of differentially expressed genes (DEGs) differed markedly between models with and without accounting for cell-type enrichments (CE) in the samples that were predicted from the RNA-seq data. This indicates that differences in cell composition in tissues that consist of multiple cell types, such as tonsil, can have a large impact on observed differences in gene expression. Based on both the NanoString and the RNA-seq data, KS06-infected pigs showed greater activation, or less inhibition, of immune response in tonsils at 42 dpi than NVSL-infected pigs, with and without accounting for CE. This suggests that the NVSL virus may be better than the KS06 virus at evading host immune response and persists in tonsils by weakening, or preventing, host immune responses. Pigs with high viral levels showed larger CE of immune cells than low viral level pigs, potentially to trigger stronger immune responses. Presence of high tonsil virus was associated with a stronger immune response, especially innate immune response through interferon signaling, but these differences were not significant when accounting for CE. Genotype at WUR was associated with different effects on immune response in tonsils of pigs during the persistence stage, depending on viral isolate and tonsil viral level. CONCLUSIONS: Results of this study provide insights into the effects of PRRSV isolate, tonsil viral level, and WUR genotype on host immune response and into potential mechanisms of PRRSV persistence in tonsils that could be targeted to improve strategies to reduce viral rebreaks. Finally, to understand transcriptome responses in tissues that consist of multiple cell types, it is important to consider differences in cell composition.

Systems Biology-Derived Genetic Signatures of Mastitis in Dairy Cattle: A New Avenue for Drug Repurposing.

Mastitis, a disease with high incidence worldwide, is the most prevalent and costly disease in the dairy industry. Gram-negative bacteria such as Escherichia coli (E. coli) are assumed to be among the leading agents causing acute severe infection with clinical signs. E. Coli, environmental mastitis pathogens, are the primary etiological agents of bovine mastitis in well-managed dairy farms. Response to E. Coli infection has a complex pattern affected by genetic and environmental parameters. On the other hand, the efficacy of antibiotics and/or anti-inflammatory treatment in E. coli mastitis is still a topic of scientific debate, and studies on the treatment of clinical cases show conflicting results. Unraveling the bio-signature of mastitis in dairy cattle can open new avenues for drug repurposing. In the current research, a novel, semi-supervised heterogeneous label propagation algorithm named Heter-LP, which applies both local and global network features for data integration, was used to potentially identify novel therapeutic avenues for the treatment of E. coli mastitis. Online data repositories relevant to known diseases, drugs, and gene targets, along with other specialized biological information for E. coli mastitis, including critical genes with robust bio-signatures, drugs, and related disorders, were used as input data for analysis with the Heter-LP algorithm. Our research identified novel drugs such as Glibenclamide, Ipratropium, Salbutamol, and Carbidopa as possible therapeutics that could be used against E. coli mastitis. Predicted relationships can be used by pharmaceutical scientists or veterinarians to find commercially efficacious medicines or a combination of two or more active compounds to treat this infectious disease.

Whole Genome Sequence Data Provides Novel Insights Into the Genetic Architecture of Meat Quality Traits in Beef.

Tenderness is a major quality attribute for fresh beef steaks in the United States, and meat quality traits in general are suitable candidates for genomic research. The objectives of the present analysis were to (1) perform genome-wide association (GWA) analysis for marbling, Warner-Bratzler shear force (WBSF), tenderness, and connective tissue using whole-genome data in an Angus population, (2) identify enriched pathways in each GWA analysis; (3) construct a protein-protein interaction network using the associated genes and (4) perform a µ-calpain proteolysis assessment for associated structural proteins. An Angus-sired population of 2,285 individuals was assessed. Animals were transported to a commercial packing plant and harvested at an average age of 457 ± 46 days. After 48 h postmortem, marbling was recorded by graders' visual appraisal. Two 2.54-cm steaks were sampled from each muscle for recording of WBSF, and tenderness, and connective tissue by a sensory panel. The relevance of additive effects on marbling, WBSF, tenderness, and connective tissue was evaluated on a genome-wide scale using a two-step mixed model-based approach in single-trait analysis. A tissue-restricted gene enrichment was performed for each GWA where all polymorphisms with an association p-value lower than 1 × 10-3 were included. The genes identified as associated were included in a protein-protein interaction network and a candidate structural protein assessment of proteolysis analyses. A total of 1,867, 3,181, 3,926, and 3,678 polymorphisms were significantly associated with marbling, WBSF, tenderness, and connective tissue, respectively. The associate region on BTA29 (36,432,655-44,313,046 bp) harbors 13 highly significant markers for meat quality traits. Enrichment for the GO term GO:0005634 (Nucleus), which includes transcription factors, was evident. The final protein-protein network included 431 interations between 349 genes. The 42 most important genes based on significance that encode structural proteins were included in a proteolysis analysis, and 81% of these proteins were potential µ-Calpain substrates. Overall, this comprehensive study unraveled genetic variants, genes and mechanisms of action responsible for the variation in meat quality traits. Our findings can provide opportunities for improving meat quality in beef cattle via marker-assisted selection.

Haplotype-based genome-wide association studies for carcass and growth traits in chicken.

There have been several genome-wide association study (GWAS) reported for carcass, growth, and meat traits in chickens. Most of these studies have been based on single SNPs GWAS. In contrast, haplotype-based GWAS reports have been limited. In the present study, 2 Northeast Agricultural University broiler lines divergently selected for abdominal fat content (NEAUHLF) and genotyped with the chicken 60K SNP chip were used to perform a haplotype-based GWAS. The lean and fat chicken lines were selected for abdominal fat content for 11 yr. Abdominal fat weight was significantly different between the 2 lines; however, there was no difference for body weight between the lean and fat lines. A total of 132 haplotype windows were significantly associated with abdominal fat weight. These significantly associated haplotype windows were primarily located on chromosomes 2, 4, 8, 10, and 26. Seven candidate genes, including SHH, LMBR1, FGF7, IL16, PLIN1, IGF1R, and SLC16A1, were located within these associated regions. These genes may play important roles in the control of abdominal fat content. Two regions on chromosomes 3 and 10 were significantly associated with testis weight. These 2 regions were previously detected by the single SNP GWAS using this same resource population. TCF21 on chromosome 3 was identified as a potentially important candidate gene for testis growth and development based on gene expression analysis and the reported function of this gene. TCF12, which was previously detected in our SNP by SNP interaction analysis, was located in a region on chromosome 10 that was significantly associated with testis weight. Six candidate genes, including TNFRSF1B, PLOD1, NPPC, MTHFR, EPHB2, and SLC35A3, on chromosome 21 may play important roles in bone development based on the known function of these genes. In addition, several regions were significantly associated with other carcass and growth traits, but no candidate genes were identified. The results of the present study may be helpful in understanding the genetic mechanisms of carcass and growth traits in chickens.