Integration of high throughput DNA genotyping and RNA-sequencing data allows for the identification of genomic regions that control gene expression, known as expression quantitative trait loci (eQTL), on a whole genome scale. Intramuscular fat (IMF) content and carcass composition play important roles in metabolic and physiological processes in mammals because they influence insulin sensitivity and consequently prevalence of metabolic diseases such as obesity and type 2 diabetes. However, limited information is available on the genetic variants and mechanisms associated with IMF deposition in mammals. Thus, our hypothesis was that eQTL analyses could identify putative regulatory regions and transcription factors (TFs) associated with intramuscular fat (IMF) content traits. We performed an integrative eQTL study in skeletal muscle to identify putative regulatory regions and factors associated with intramuscular fat content traits. Data obtained from skeletal muscle samples of 192 animals was used for association analysis between 461,466 SNPs and the transcription level of 11,808 genes. This yielded 1268 cis- and 10,334 trans-eQTLs, among which we identified nine hotspot regions that each affected the expression of > 119 genes. These putative regulatory regions overlapped with previously identified QTLs for IMF content. Three of the hotspots respectively harbored the transcription factors USF1, EGR4 and RUNX1T1, which are known to play important roles in lipid metabolism. From co-expression network analysis, we further identified modules significantly correlated with IMF content and associated with relevant processes such as fatty acid metabolism, carbohydrate metabolism and lipid metabolism. This study provides novel insights into the link between genotype and IMF content as evident from the expression level. It thereby identifies genomic regions of particular importance and associated regulatory factors. These new findings provide new knowledge about the biological processes associated with genetic variants and mechanisms associated with IMF deposition in mammals.

Clinical mastitis (CM) is an inflammatory disease occurring in the mammary glands of lactating cows. CM is under genetic control, and a prominent CM resistance QTL located on chromosome 6 was reported in various dairy cattle breeds. Nevertheless, the biological mechanism underpinning this QTL has been lacking. Herein, we mapped, fine-mapped, and discovered the putative causal variant underlying this CM resistance QTL in the Dutch dairy cattle population. We identified a ~12 kb multi-allelic copy number variant (CNV), that is in perfect linkage disequilibrium with a lead SNP, as a promising candidate variant. By implementing a fine-mapping and through expression QTL mapping, we showed that the group-specific component gene (GC), a gene encoding a vitamin D binding protein, is an excellent candidate causal gene for the QTL. The multiplicated alleles are associated with increasedGCexpression and low CM resistance. Ample evidence from functional genomics data supports the presence of an enhancer within this CNV, which would exertcis-regulatory effect onGC. We observed that strong positive selection swept the region near the CNV, and haplotypes associated with the multiplicated allele were strongly selected for. Moreover, the multiplicated allele showed pleiotropic effects for increased milk yield and reduced fertility, hinting that a shared underlying biology for these effects may revolve around the vitamin D pathway. These findings together suggest a putative causal variant of a CM resistance QTL, where acis-regulatory element located within a CNV can alter gene expression and affect multiple economically important traits.

Chicken feed efficiency (FE) traits are the most important economic traits in broiler production. Several studies evaluating genetic factors affecting food consumption in chickens are available. However, most of these studies identified genomic regions containing putative quantitative trait loci for each trait separately. It is still a challenge to find common gene networks related to these traits. Therefore, here, a genome-wide association study (GWAS) was conducted to explore candidate genomic regions responsible for Feed Intake (FI), Body Weight Gain (BWG) and Feed Conversion Ratio (FCR) traits and their gene networks. A total of 1430 broilers from an experimental population was genotyped with the high density Affymetrix 600K SNP array. A total of 119 associated SNPs located in 20 chromosomes were identified, where some of them were common in more than one FE trait. In addition, novel genomic regions were prospected considering the SNPs dominance effects and sex interaction, identifying putative candidate genes only when these effects were fit in the model. Relevant candidate genes such as ATRNL1, PIK3C2A, PTPRN2, SORCS3 and gga-mir-1759 were highlighted in this study helping to elucidate the genomic architecture of feed efficiency traits. These results provide new insights on the mechanisms underlying the consumption and utilization of food in chickens.

Meat and egg-type chickens have been selected for several generations for different traits. Artificial and natural selection for different phenotypes can change frequency of genetic variants, leaving particular genomic footprints throghtout the genome. Thus, the aims of this study were to sequence 28 chickens from two Brazilian lines (meat and white egg-type) and use this information to characterize genome-wide genetic variations, identify putative regions under selection using Fst method, and find putative pathways under selection. A total of 13.93 million SNPs and 1.36 million INDELs were identified, with more variants detected from the broiler (meat-type) line. Although most were located in non-coding regions, we identified 7255 intolerant non-synonymous SNPs, 512 stopgain/loss SNPs, 1381 frameshift and 1094 non-frameshift INDELs that may alter protein functions. Genes harboring intolerant non-synonymous SNPs affected metabolic pathways related mainly to reproduction and endocrine systems in the white-egg layer line, and lipid metabolism and metabolic diseases in the broiler line. Fst analysis in sliding windows, using SNPs and INDELs separately, identified over 300 putative regions of selection overlapping with more than 250 genes. For the first time in chicken, INDEL variants were considered for selection signature analysis, showing high level of correlation in results between SNP and INDEL data. The putative regions of selection signatures revealed interesting candidate genes and pathways related to important phenotypic traits in chicken, such as lipid metabolism, growth, reproduction, and cardiac development. In this study, Fst method was applied to identify high confidence putative regions under selection, providing novel insights into selection footprints that can help elucidate the functional mechanisms underlying different phenotypic traits relevant to meat and egg-type chicken lines. In addition, we generated a large catalog of line-specific and common genetic variants from a Brazilian broiler and a white egg layer line that can be used for genomic studies involving association analysis with phenotypes of economic interest to the poultry industry.

MicroRNAs (miRNAs) are small noncoding RNAs of approximately 22 nucleotides, highly conserved among species, which modulate gene expression by cleaving messenger RNA target or inhibiting translation. MiRNAs are involved in the regulation of many processes including cell proliferation, differentiation, neurogenesis, angiogenesis, and apoptosis. Beef tenderness is an organoleptic characteristic of great influence in the acceptance of meat by consumers. Previous studies have shown that collagen level, marbling, apoptosis and proteolysis are among the many factors that affect beef tenderness. Considering that miRNAs can modulate gene expression, this study was designed to identify differentially expressed miRNAs that could be modulating biological processes involved with beef tenderness. Deep sequence analysis of miRNA libraries from longissimus thoracis muscle allowed the identification of 42 novel and 308 known miRNAs. Among the known miRNAs, seven were specifically expressed in skeletal muscle. Differential expression analysis between animals with high (H) and low (L) estimated breeding values for shear force (EBVSF) revealed bta-mir-182 and bta-mir-183 are up-regulated (q value < 0.05) in animals with L EBVSF, and bta-mir-338 is up-regulated in animals with H EBVSF. The number of bovine predicted targets for bta-mir-182, bta-mir-183 and bta-mir-338 were 811, 281 and 222, respectively, which correspond to 1204 unique target genes. Among these, four of them, MEF2C, MAP3K2, MTDH and TNRC6B were common targets of the three differentially expressed miRNAs. The functional analysis identified important pathways related to tenderness such as apoptosis and the calpain–calpastatin system. The results obtained indicate the importance of miRNAs in the regulatory mechanisms that influence muscle proteolysis and meat tenderness and contribute to our better understanding of the role of miRNAs in biological processes associated with beef tenderness.

Natural and artificial selection leads to changes in certain regions of the genome resulting in selection signatures that can reveal genes associated with the selected traits. Selection signatures may be identified using different methodologies, of which some are based on detecting contiguous sequences of homozygous identical-by-descent haplotypes, called runs of homozygosity (ROH), or estimating fixation index (FST) of genomic windows that indicates genetic differentiation. This study aimed to identify selection signatures in a paternal broiler TT line at generations 7th and 16th of selection and to investigate the genes annotated in these regions as well as the biological pathways involved. For such purpose, ROH and FST-based analysis were performed using whole genome sequence of twenty-eight chickens from two different generations. ROH analysis identified homozygous regions of short and moderate size. Analysis of ROH patterns revealed regions commonly shared among animals and changes in ROH abundance and size between the two generations. Results also suggest that whole genome sequencing (WGS) outperforms SNPchip data avoiding overestimation of ROH size and underestimation of ROH number; however, sequencing costs can limited the number of animals analyzed. FST-based analysis revealed genetic differentiation in several genomic windows. Annotation of the consensus regions of ROH and FST windows revealed new and previously identified genes associated with traits of economic interest, such as APOB, IGF1, IGFBP2, POMC, PPARG, and ZNF423. Over-representation analysis of the genes resulted in biological terms of skeletal muscle, matrilin proteins, adipose tissue, hyperglycemia, diabetes, Salmonella infections and tyrosine. Identification of ROH and FST-based analyses revealed selection signatures in TT line and genes that have important role in traits of economic interest. Changes in the genome of the chickens were observed between the 7th and 16th generations showing that ancient and recent selection in TT line may have acted over genomic regions affecting diseases and performance traits.