de novo variant calling identifies cancer mutation signatures in the 1000 Genomes Project.

Detection of de novo variants (DNVs) is critical for studies of disease-related variation and mutation rates. To accelerate DNV calling, we developed a graphics processing units-based workflow. We applied our workflow to whole-genome sequencing data from three parent-child sequenced cohorts including the Simons Simplex Collection (SSC), Simons Foundation Powering Autism Research (SPARK), and the 1000 Genomes Project (1000G) that were sequenced using DNA from blood, saliva, and lymphoblastoid cell lines (LCLs), respectively. The SSC and SPARK DNV callsets were within expectations for number of DNVs, percent at CpG sites, phasing to the paternal chromosome of origin, and average allele balance. However, the 1000G DNV callset was not within expectations and contained excessive DNVs that are likely cell line artifacts. Mutation signature analysis revealed 30% of 1000G DNV signatures matched B-cell lymphoma. Furthermore, we found variants in DNA repair genes and at Clinvar pathogenic or likely-pathogenic sites and significant excess of protein-coding DNVs in IGLL5; a gene known to be involved in B-cell lymphomas. Our study provides a new rapid DNV caller for the field and elucidates important implications of using sequencing data from LCLs for reference building and disease-related projects.

Identification of putative regulatory regions and transcription factors associated with intramuscular fat content traits.

BACKGROUND: 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. RESULTS: 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. CONCLUSION: 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.

Gene Co-expression Analysis Indicates Potential Pathways and Regulators of Beef Tenderness in Nellore Cattle.

Beef tenderness, a complex trait affected by many factors, is economically important to beef quality, industry, and consumer's palatability. In this study, RNA-Seq was used in network analysis to better understand the biological processes that lead to differences in beef tenderness. Skeletal muscle transcriptional profiles from 24 Nellore steers, selected by extreme estimated breeding values (EBVs) for shear force after 14 days of aging, were analyzed and 22 differentially expressed transcripts were identified. Among these were genes encoding ribosomal proteins, glutathione transporter ATP-binding cassette, sub-family C (CFTR/MRP), member 4 (ABCC4), and synaptotagmin IV (SYT4). Complementary co-expression analyses using Partial Correlation with Information Theory (PCIT), Phenotypic Impact Factor (PIF) and the Regulatory Impact Factor (RIF) methods identified candidate regulators and related pathways. The PCIT analysis identified ubiquitin specific peptidase 2 (USP2), growth factor receptor-bound protein 10 (GBR10), anoctamin 1 (ANO1), and transmembrane BAX inhibitor motif containing 4 (TMBIM4) as the most differentially hubbed (DH) transcripts. The transcripts that had a significant correlation with USP2, GBR10, ANO1, and TMBIM4 enriched for proteasome KEGG pathway. RIF analysis identified microRNAs as candidate regulators of variation in tenderness, including bta-mir-133a-2 and bta-mir-22. Both microRNAs have target genes present in the calcium signaling pathway and apoptosis. PIF analysis identified myoglobin (MB), enolase 3 (ENO3), and carbonic anhydrase 3 (CA3) as potentially having fundamental roles in tenderness. Pathways identified in our study impacted in beef tenderness included: calcium signaling, apoptosis, and proteolysis. These findings underscore some of the complex molecular mechanisms that control beef tenderness in Nellore cattle.

Genetic architecture of gene expression underlying variation in host response to porcine reproductive and respiratory syndrome virus infection.

It has been shown that inter-individual variation in host response to porcine reproductive and respiratory syndrome (PRRS) has a heritable component, yet little is known about the underlying genetic architecture of gene expression in response to PRRS virus (PRRSV) infection. Here, we integrated genome-wide genotype, gene expression, viremia level, and weight gain data to identify genetic polymorphisms that are associated with variation in inter-individual gene expression and response to PRRSV infection in pigs. RNA-seq analysis of peripheral blood samples collected just prior to experimental challenge (day 0) and at 4, 7, 11 and 14 days post infection from 44 pigs revealed 6,430 differentially expressed genes at one or more time points post infection compared to the day 0 baseline. We mapped genetic polymorphisms that were associated with inter-individual differences in expression at each day and found evidence of cis-acting expression quantitative trait loci (cis-eQTL) for 869 expressed genes (qval < 0.05). Associations between cis-eQTL markers and host response phenotypes using 383 pigs suggest that host genotype-dependent differences in expression of GBP5, GBP6, CCHCR1 and CMPK2 affect viremia levels or weight gain in response to PRRSV infection.