Novel putative causal mutations associated with fat traits in Nellore cattle uncovered by eQTLs located in open chromatin regions.

Intramuscular fat (IMF) and backfat thickness (BFT) are critical economic traits impacting meat quality. However, the genetic variants controlling these traits need to be better understood. To advance knowledge in this area, we integrated RNA-seq and single nucleotide polymorphisms (SNPs) identified in genomic and transcriptomic data to generate a linkage disequilibrium filtered panel of 553,581 variants. Expression quantitative trait loci (eQTL) analysis revealed 36,916 cis-eQTLs and 14,408 trans-eQTLs. Association analysis resulted in three eQTLs associated with BFT and 24 with IMF. Functional enrichment analysis of genes regulated by these 27 eQTLs revealed noteworthy pathways that can play a fundamental role in lipid metabolism and fat deposition, such as immune response, cytoskeleton remodeling, iron transport, and phospholipid metabolism. We next used ATAC-Seq assay to identify and overlap eQTL and open chromatin regions. Six eQTLs were in regulatory regions, four in predicted insulators and possible CCCTC-binding factor DNA binding sites, one in an active enhancer region, and the last in a low signal region. Our results provided novel insights into the transcriptional regulation of IMF and BFT, unraveling putative regulatory variants.

Allele-specific expression reveals functional SNPs affecting muscle-related genes in bovine.

Single nucleotide polymorphisms showing allele-specific expression (ASE SNPs) are useful for cis-regulatory variants discovery. Despite this potential, there are expensive costs involved in genome-level ASE analysis for large sample sizes. If different data resolutions are available, genotype imputation can be used to mitigate this limitation. Aiming to increase the power to detect regulatory variants, we used a large dataset (>4 million) of imputed SNP genotypes and RNA-Seq data from 190 Nelore steers. Differences between major and minor allele expressions in muscle were tested with a Binomial Test. We identified 38,177 ASE SNPs (FDR ≤ 0.05) within 7304 linkage disequilibrium blocks. After that, we searched for aseQTLs (i.e., neighboring SNPs potentially regulating the ASE SNPs' allelic expression) by comparing the ASE of heterozygous to homozygous sample groups under a Wilcoxon Rank Sum test. We identified 21,543 aseQTLs potentially regulating 430 ASE SNPs (FDR ≤ 0.05). A total of 3333 cis-eQTLs (being 2098 ASE SNPs and 1075 aseQTLs) were associated with the expression of 758 transcripts (FDR ≤ 0.05), demonstrating the cis-regulatory effect of these ASE SNPs and aseQTLs. Data integration showed reproducibility with previous studies in bovine ASE and genomic imprinting. Furthermore, we identified 36,756 novel ASE regions due to the imputation approach. Comparisons with epigenetics data from Functional Annotation of Animal Genomes (FAANG) suggest a regulatory potential of the ASE-related SNPs. The affected genes were enriched in metabolic pathways essential for muscle homeostasis. These findings reinforce the potential of using ASE for discovering cis-regulatory SNPs that may affect muscle-related traits.

Simulation of group testing scenarios can boost COVID-19 screening power.

The COVID-19 has severely affected economies and health systems around the world. Mass testing could work as a powerful alternative to restrain disease dissemination, but the shortage of reagents is a limiting factor. A solution to optimize test usage relies on 'grouping' or 'pooling' strategies, which combine a set of individuals in a single reaction. To compare different group testing configurations, we developed the poolingr package, which performs an innovative hybrid in silico/in vitro approach to search for optimal testing configurations. We used 6759 viral load values, observed in 2389 positive individuals, to simulate a wide range of scenarios. We found that larger groups (>100) framed into multi-stage setups (up to six stages) could largely boost the power to detect spreaders. Although the boost was dependent on the disease prevalence, our method could point to cheaper grouping schemes to better mitigate COVID-19 dissemination through identification and quarantine recommendation for positive individuals.

Genome-wide detection of CNVs and their association with performance traits in broilers.

BACKGROUND: Copy number variations (CNVs) are a major type of structural genomic variants that underlie genetic architecture and phenotypic variation of complex traits, not only in humans, but also in livestock animals. We identified CNVs along the chicken genome and analyzed their association with performance traits. Genome-wide CNVs were inferred from Affymetrix® high density SNP-chip data for a broiler population. CNVs were concatenated into segments and association analyses were performed with linear mixed models considering a genomic relationship matrix, for birth weight, body weight at 21, 35, 41 and 42 days, feed intake from 35 to 41 days, feed conversion ratio from 35 to 41 days and, body weight gain from 35 to 41 days of age. RESULTS: We identified 23,214 autosomal CNVs, merged into 5042 distinct CNV regions (CNVRs), covering 12.84% of the chicken autosomal genome. One significant CNV segment was associated with BWG on GGA3 (q-value = 0.00443); one significant CNV segment was associated with BW35 (q-value = 0.00571), BW41 (q-value = 0.00180) and BW42 (q-value = 0.00130) on GGA3, and one significant CNV segment was associated with BW on GGA5 (q-value = 0.00432). All significant CNV segments were verified by qPCR, and a validation rate of 92.59% was observed. These CNV segments are located nearby genes, such as KCNJ11, MyoD1 and SOX6, known to underlie growth and development. Moreover, gene-set analyses revealed terms linked with muscle physiology, cellular processes regulation and potassium channels. CONCLUSIONS: Overall, this CNV-based GWAS study unravels potential candidate genes that may regulate performance traits in chickens. Our findings provide a foundation for future functional studies on the role of specific genes in regulating performance in chickens.

Widespread modulation of gene expression by copy number variation in skeletal muscle.

Copy number variation (CNV) is a frequently observed deviation from the diploid state due to duplication or deletion of genomic regions. Although intensively analyzed for association with diseases and production traits, the specific mechanisms and extent by which such variations affect the phenotype are incompletely understood. We present an integrative study on CNV and genome-wide gene expression in Brazilian Bos indicus cattle. We analyzed CNVs inferred from SNP-chip data for effects on gene expression measured with RNA-seq in skeletal muscle samples of 183 steers. Local effects, where expression changes coincided with CNVs in the respective genes, were restricted to immune genes. Distal effects were attributable to several high-impact CNVs that modulated remote expression in an orchestrated and intertwined fashion. These CNVs were located in the vicinity of major skeletal muscle pathway regulators and associated genes were enriched for proteolysis, autophagy, and muscle structure development. From association analysis between CNVs and several meat quality and production traits, we found CNV-associated expression effects to also manifest at the phenotype level. Based on genome sequences of the population founders, we further demonstrate that CNVs with impact on expression and phenotype are passed on from one generation to another.