High-throughput characterization of functional variants highlights heterogeneity and polygenicity underlying lung cancer susceptibility.

Genome-wide association studies (GWASs) have identified numerous lung cancer risk-associated loci. However, decoding molecular mechanisms of these associations is challenging since most of these genetic variants are non-protein-coding with unknown function. Here, we implemented massively parallel reporter assays (MPRAs) to simultaneously measure the allelic transcriptional activity of risk-associated variants. We tested 2,245 variants at 42 loci from 3 recent GWASs in East Asian and European populations in the context of two major lung cancer histological types and exposure to benzo(a)pyrene. This MPRA approach identified one or more variants (median 11 variants) with significant effects on transcriptional activity at 88% of GWAS loci. Multimodal integration of lung-specific epigenomic data demonstrated that 63% of the loci harbored multiple potentially functional variants in linkage disequilibrium. While 22% of the significant variants showed allelic effects in both A549 (adenocarcinoma) and H520 (squamous cell carcinoma) cell lines, a subset of the functional variants displayed a significant cell-type interaction. Transcription factor analyses nominated potential regulators of the functional variants, including those with cell-type-specific expression and those predicted to bind multiple potentially functional variants across the GWAS loci. Linking functional variants to target genes based on four complementary approaches identified candidate susceptibility genes, including those affecting lung cancer cell growth. CRISPR interference of the top functional variant at 20q13.33 validated variant-to-gene connections, including RTEL1, SOX18, and ARFRP1. Our data provide a comprehensive functional analysis of lung cancer GWAS loci and help elucidate the molecular basis of heterogeneity and polygenicity underlying lung cancer susceptibility.

Predicting functional UTR variants by integrating region-specific features.

The untranslated region (UTR) of messenger ribonucleic acid (mRNA), including the 5'UTR and 3'UTR, plays a critical role in regulating gene expression and translation. Variants within the UTR can lead to changes associated with human traits and diseases; however, computational prediction of UTR variant effect is challenging. Current noncoding variant prediction mainly focuses on the promoters and enhancers, neglecting the unique sequence of the UTR and thereby limiting their predictive accuracy. In this study, using consolidated datasets of UTR variants from disease databases and large-scale experimental data, we systematically analyzed more than 50 region-specific features of UTR, including functional elements, secondary structure, sequence composition and site conservation. Our analysis reveals that certain features, such as C/G-related sequence composition in 5'UTR and A/T-related sequence composition in 3'UTR, effectively differentiate between nonfunctional and functional variant sets, unveiling potential sequence determinants of functional UTR variants. Leveraging these insights, we developed two classification models to predict functional UTR variants using machine learning, achieving an area under the curve (AUC) value of 0.94 for 5'UTR and 0.85 for 3'UTR, outperforming all existing methods. Our models will be valuable for enhancing clinical interpretation of genetic variants, facilitating the prediction and management of disease risk.

AStruct: detection of allele-specific RNA secondary structure in structuromic probing data.

BACKGROUND: Uncovering functional genetic variants from an allele-specific perspective is of paramount importance in advancing our understanding of gene regulation and genetic diseases. Recently, various allele-specific events, such as allele-specific gene expression, allele-specific methylation, and allele-specific binding, have been explored on a genome-wide scale due to the development of high-throughput sequencing methods. RNA secondary structure, which plays a crucial role in multiple RNA-associated processes like RNA modification, translation and splicing, has emerged as an essential focus of relevant research. However, tools to identify genetic variants associated with allele-specific RNA secondary structures are still lacking. RESULTS: Here, we develop a computational tool called 'AStruct' that enables us to detect allele-specific RNA secondary structure (ASRS) from RT-stop based structuromic probing data. AStruct shows robust performance in both simulated datasets and public icSHAPE datasets. We reveal that single nucleotide polymorphisms (SNPs) with higher AStruct scores are enriched in coding regions and tend to be functional. These SNPs are highly conservative, have the potential to disrupt sites involved in m6A modification or protein binding, and are frequently associated with disease. CONCLUSIONS: AStruct is a tool dedicated to invoke allele-specific RNA secondary structure events at heterozygous SNPs in RT-stop based structuromic probing data. It utilizes allelic variants, base pairing and RT-stop information under different cell conditions to detect dynamic and functional ASRS. Compared to sequence-based tools, AStruct considers dynamic cell conditions and outperforms in detecting functional variants. AStruct is implemented in JAVA and is freely accessible at: https://github.com/canceromics/AStruct .

A regulatory element associated to NAFLD in the promoter of DIO1 controls LDL-C, HDL-C and triglycerides in hepatic cells.

BACKGROUND: Genome-wide association studies (GWAS) have identified genetic variants linked to fat metabolism and related traits, but rarely pinpoint causative variants. This limitation arises from GWAS not considering functional implications of noncoding variants that can affect transcription factor binding and potentially regulate gene expression. The aim of this study is to investigate a candidate noncoding functional variant within a genetic locus flagged by a GWAS SNP associated with non-alcoholic fatty liver disease (NAFLD), a condition characterized by liver fat accumulation in non-alcohol consumers. METHODS: CRISPR-Cas9 gene editing in HepG2 cells was used to modify the regulatory element containing the candidate functional variant linked to NAFLD. Global gene expression in mutant cells was assessed through RT-qPCR and targeted transcriptomics. A phenotypic assay measured lipid droplet accumulation in the CRISPR-Cas9 mutants. RESULTS: The candidate functional variant, rs2294510, closely linked to the NAFLD-associated GWAS SNP rs11206226, resided in a regulatory element within the DIO1 gene's promoter region. Altering this element resulted in changes in transcription factor binding sites and differential expression of candidate target genes like DIO1, TMEM59, DHCR24, and LDLRAD1, potentially influencing the NAFLD phenotype. Mutant HepG2 cells exhibited increased lipid accumulation, a hallmark of NAFLD, along with reduced LDL-C, HDL-C and elevated triglycerides. CONCLUSIONS: This comprehensive approach, that combines genome editing, transcriptomics, and phenotypic assays identified the DIO1 promoter region as a potential enhancer. Its activity could regulate multiple genes involved in the NAFLD phenotype or contribute to defining a polygenic risk score for enhanced risk assessment in NAFLD patients.

Archaic humans have contributed to large-scale variation in modern human T cell receptor genes.

Human T cell receptors (TCRs) are critical for mediating immune responses to pathogens and tumors and regulating self-antigen recognition. Yet, variations in the genes encoding TCRs remain insufficiently defined. Detailed analysis of expressed TCR alpha, beta, gamma, and delta genes in 45 donors from four human populations-African, East Asian, South Asian, and European-revealed 175 additional TCR variable and junctional alleles. Most of these contained coding changes and were present at widely differing frequencies in the populations, a finding confirmed using DNA samples from the 1000 Genomes Project. Importantly, we identified three Neanderthal-derived, introgressed TCR regions including a highly divergent TRGV4 variant, which mediated altered butyrophilin-like molecule 3 (BTNL3) ligand reactivity and was frequent in all modern Eurasian population groups. Our results demonstrate remarkable variation in TCR genes in both individuals and populations, providing a strong incentive for including allelic variation in studies of TCR function in human biology.

Detection of potential functional variants based on systems-biology: the case of feed efficiency in beef cattle.

BACKGROUND: Potential functional variants (PFVs) can be defined as genetic variants responsible for a given phenotype. Ultimately, these are the best DNA markers for animal breeding and selection, especially for polygenic and complex phenotypes. Herein, we described the identification of PFVs for complex phenotypes (in this case, Feed Efficiency in beef cattle) using a systems-biology driven approach based on RNA-seq data from physiologically relevant organs. RESULTS: The systems-biology coupled with deep molecular phenotyping by RNA-seq of liver, muscle, hypothalamus, pituitary, and adrenal glands of animals with high and low feed efficiency (FE) measured by residual feed intake (RFI) identified 2,000,936 uniquely variants. Among them, 9986 variants were significantly associated with FE and only 78 had a high impact on protein expression and were considered as PFVs. A set of 169 significant uniquely variants were expressed in all five organs, however, only 27 variants had a moderate impact and none of them a had high impact on protein expression. These results provide evidence of tissue-specific effects of high-impact PFVs. The PFVs were enriched (FDR < 0.05) for processing and presentation of MHC Class I and II mediated antigens, which are an important part of the adaptive immune response. The experimental validation of these PFVs was demonstrated by the increased prediction accuracy for RFI using the weighted G matrix (ssGBLUP+wG; Acc = 0.10 and b = 0.48) obtained in the ssGWAS in comparison to the unweighted G matrix (ssGBLUP; Acc = 0.29 and b = 1.10). CONCLUSION: Here we identified PFVs for FE in beef cattle using a strategy based on systems-biology and deep molecular phenotyping. This approach has great potential to be used in genetic prediction programs, especially for polygenic phenotypes.

The Upregulation of Cathepsin G Is Associated with Resistance to Bovine Paratuberculosis.

An in silico genomic-transcriptomic combined approach allowed the identification of a polymorphism (cis-eQTL-rs41976219) in the Bos taurus genome associated with the CTSG mRNA expression in bovine blood samples, which suggests that individual genetic variation might modulate the CTSG transcriptional response. In the current study, a sandwich ELISA is used to measure the CTSG protein levels in supernatants of monocyte-derived macrophages (MDMs) isolated from cows with the AA (n = 5) and AC (n = 11) genotypes for the rs41976219 and infected ex vivo with MAP. Cows with the AC genotype have significantly higher CTSG protein levels (1.85 ng/mL) in the supernatants of enriched CD14+-MDMs after 2 h of infection when compared with infected CD14+-MDMs from cows with the AA genotype (1.68 ng/mL). Statistically significant differences in the intracellular MAP load at 7 d p.i. are observed between animals with the AA (2.16 log CFUs) and AC (1.44 log CFUs) genotypes. Finally, the association between the rs41976219 allelic variants and resistance to PTB is tested in a larger cattle population (n = 943) classified according to the presence (n = 442) or absence (n = 501) of PTB-associated lesions. The presence of the two minor alleles in the rs41976219 (CC) is more frequent among healthy cows than in cows with PTB-associated lesions in gut tissues (2.2% vs. 1.4%, OR = 0.61). In agreement with this, the CTSG levels in plasma samples of cows without lesions in gut tissues and with the CC (n = 8) genotype are significantly higher than in the plasmas of cows with the AA + AC (n = 36) genotypes.

Searching for Genetic Biomarkers for Hereditary Angioedema Due to C1-Inhibitor Deficiency (C1-INH-HAE).

Existing evidence indicates that modifier genes could change the phenotypic outcome of the causal SERPING1 variant and thus explain the expression variability of hereditary angioedema due to C1-inhibitor deficiency (C1-INH-HAE). To further examine this hypothesis, we investigated the presence or absence of 18 functional variants of genes encoding proteins involved in the metabolism and function of bradykinin, the main mediator of C1-INH-HAE attacks, in relation to three distinct phenotypic traits of patients with C1-INH-HAE, i.e., the age at disease onset, the need for long-term prophylaxis (LTP), and the severity of the disease. Genetic analyses were performed by a validated next-generation sequencing platform. In total, 233 patients with C1-INH-HAE from 144 unrelated families from five European countries were enrolled in the study. Already described correlations between five common functional variants [F12-rs1801020, KLKB1-rs3733402, CPN1-rs61751507, and two in SERPING1 (rs4926 and rs28362944)] and C1-INH-HAE severity were confirmed. Furthermore, significant correlations were found between either the age at disease onset, the LTP, or the severity score of the disease and a series of other functional variants (F13B-rs6003, PLAU-rs2227564, SERPINA1-rs28929474, SERPINA1-rs17580, KLK1-rs5515, SERPINE1-rs6092, and F2-rs1799963). Interestingly, correlations uncovered in the entire cohort of patients were different from those discovered in the cohort of patients carrying missense causal SERPING1 variants. Our findings indicate that variants other than the SERPING1 causal variants act as independent modifiers of C1-INH-HAE severity and could be tested as possible prognostic biomarkers.

Breeding Sustainable Beef Cows: Reducing Weight and Increasing Productivity.

Programs for sustainable beef production are established, but the specific role of beef cows in these systems is not well defined. This work characterized cows for two traits related to sustainability, cow weight (CW) and cumulative weight weaned (WtW). Cow weight indicates nutrient requirements and enteric methane emissions. Cumulative weight weaned reflects reproductive performance and avoidance of premature culling for characteristics related to animal health, welfare, and worker safety. Both traits were evaluated with random regression models with records from a crossbred population representing 18 breeds that conduct US national cattle evaluations. The genomic REML analyses included additive and dominance components, with relationships among 22,776 animals constructed from genotypes of 181,286 potentially functional variants imputed from a low-pass sequence. Projected to 8 years of age, the additive heritability estimate for CW was 0.57 and 0.11 for WtW. Dominance heritability was 0.02 for CW and 0.19 for WtW. Many variants with significant associations with CW were within previously described quantitative trait loci (QTL) for growth-related production, meat, and carcass traits. Significant additive WtW variants were covered by QTL for traits related to reproduction and structural soundness. All breeds contributed to groups of cows with high and low total genetic values (additive + dominance effects) for both traits. The high WtW cows and cows above the WtW mean but below the CW mean had larger heterosis values and fewer bases in runs of homozygosity. The high additive heritability of CW and dominance effects on WtW indicate that breeding to improve beef cow sustainability should involve selection to reduce CW and mate selection to maintain heterosis and reduce runs of homozygosity.

A systematic review and functional bioinformatics analysis of genes associated with Crohn's disease identify more than 120 related genes.

BACKGROUND: Crohn's disease is one of the two categories of inflammatory bowel diseases that affect the gastrointestinal tract. The heritability estimate has been reported to be 0.75. Several genes linked to Crohn's disease risk have been identified using a plethora of strategies such as linkage-based studies, candidate gene association studies, and lately through genome-wide association studies (GWAS). Nevertheless, to our knowledge, a compendium of all the genes that have been associated with CD is lacking. METHODS: We conducted functional analyses of a gene set generated from a systematic review where genes potentially related to CD found in the literature were analyzed and classified depending on the genetic evidence reported and putative biological function. For this, we retrieved and analyzed 2496 abstracts comprising 1067 human genes plus 22 publications regarding 133 genes from GWAS Catalog. Then, each gene was curated and categorized according to the type of evidence associated with Crohn's disease. RESULTS: We identified 126 genes associated with Crohn's disease risk by specific experiments. Additionally, 71 genes were recognized associated through GWAS alone, 18 to treatment response, 41 to disease complications, and 81 to related diseases. Bioinformatic analysis of the 126 genes supports their importance in Crohn's disease and highlights genes associated with specific aspects such as symptoms, drugs, and comorbidities. Importantly, most genes were not included in commercial genetic panels suggesting that Crohn's disease is genetically underdiagnosed. CONCLUSIONS: We identified a total of 126 genes from PubMed and 71 from GWAS that showed evidence of association to diagnosis, 18 to treatment response, and 41 to disease complications in Crohn's disease. This prioritized gene catalog can be explored at http://victortrevino.bioinformatics.mx/CrohnDisease .

Genetic epidemiological analysis of hypouricaemia from 4993 Japanese on non-functional variants of URAT1/SLC22A12 gene.

OBJECTIVES: Up to 0.3% of Japanese have hypouricaemia. Most cases appear to result from a hereditary disease, renal hypouricaemia (RHUC), which causes exercise-induced acute kidney injury and urolithiasis. However, to what extent RHUC accounts for hypouricaemia is not known. We therefore investigated its frequency and evaluated its risks by genotyping a general Japanese population. METHODS: A cohort of 4993 Japanese was examined by genotyping the non-functional variants R90H (rs121907896) and W258X (rs121907892) of URAT1/SLC22A12, the two most common causative variants of RHUC in Japanese. RESULTS: Participants' fractional excretion of uric acid and risk allele frequencies markedly increased at lower serum uric acid (SUA) levels. Ten participants (0.200%) had an SUA level ≤2.0 mg/dl and nine had R90H or W258X and were likely to have RHUC. Logistic regression analysis revealed these URAT1 variants to be significantly and independently associated with the risk of hypouricaemia and mild hypouricaemia (SUA ≤3.0 mg/dl) as well as sex, age and BMI, but these URAT1 variants were the only risks in the hypouricaemia population (SUA ≤2.0 mg/dl). W258X was only a risk in males with SUA ≤3.0 mg/dl. CONCLUSION: Our study accurately reveals the prevalence of RHUC and provides genetic evidence for its definition (SUA ≤2.0 mg/dl). We also show that individuals with SUA ≤3.0 mg/dl, especially males, are prone to RHUC. Our findings will help to promote a better epidemiological understanding of RHUC as well as more accurate diagnosis, especially in males with mild hypouricaemia.

GATA1 pathogenic variants disrupt MYH10 silencing during megakaryopoiesis.

BACKGROUND: GATA1 is an essential transcription factor for both polyploidization and megakaryocyte (MK) differentiation. The polyploidization defect observed in GATA1 variant carriers is not well understood. OBJECTIVE: To extensively phenotype two pedigrees displaying different variants in the GATA1 gene and determine if GATA1 controls MYH10 expression levels, a key modulator of MK polyploidization. METHOD: A total of 146 unrelated propositi with constitutional thrombocytopenia were screened on a multigene panel. We described the genotype-phenotype correlation in GATA1 variant carriers and investigated the effect of these novel variants on MYH10 transcription using luciferase constructs. RESULTS: The clinical profile associated with the p.L268M variant localized in the C terminal zinc finger was unusual in that the patient displayed bleeding and severe platelet aggregation defects without early-onset thrombocytopenia. p.N206I localized in the N terminal zinc finger was associated, on the other hand, with severe thrombocytopenia (15G/L) in early life. High MYH10 levels were evidenced in platelets of GATA1 variant carriers. Analysis of MKs anti-GATA1 chromatin immunoprecipitation-sequencing data revealed two GATA1 binding sites, located in the 3' untranslated region and in intron 8 of the MYH10 gene. Luciferase reporter assays showed their respective role in the regulation of MYH10 gene expression. Both GATA1 variants significantly alter intron 8 driven MYH10 transcription. CONCLUSION: The discovery of an association between MYH10 and GATA1 is a novel one. Overall, this study suggests that impaired MYH10 silencing via an intronic regulatory element is the most likely cause of GATA1-related polyploidization defect.

Prediction of disease-associated functional variants in noncoding regions through a comprehensive analysis by integrating datasets and features.

One of the greatest challenges in human genetics is deciphering the link between functional variants in noncoding sequences and the pathophysiology of complex diseases. To address this issue, many methods have been developed to sort functional single-nucleotide variants (SNVs) for neutral SNVs in noncoding regions. In this study, we integrated well-established features and commonly used datasets and merged them into large-scale datasets based on a random forest model, which yielded promising performance and outperformed some cutting-edge approaches. Our analyses of feature importance and data coverage also provide certain clues for future research in enhancing the prediction of functional noncoding SNVs.

Polymorphisms of KLF3 gene coding region and identification of their functionality for abdominal fat in chickens.

KLF3 is a member of the Kruppel-like factor (KLF) family of transcription factors, and plays an important role in several biological processes, including adipogenesis, erythropoiesis and B-cell development. The purposes of this study are to search for polymorphisms of KLF3 coding region and to provide functional evidence for abdominal fat in chickens. A total of 168 SNPs in KLF3 coding region were detected in a unique chicken population, the Northeast Agricultural University broiler lines divergently selected for abdominal fat content (NEAUHLF). Of which three single nucleotide polymorphisms (g.3452T > C, g.8663A > G and g.10751G > A) were significantly correlated with abdominal fat weight (AFW) and abdominal fat percentage (AFP) of 329 birds from the 19th generation of NEAUHLF (FDR < 0.05). The reporter gene assay was performed to verify functionality of these three SNPs in both ICP-1 and DF1 cells. Results showed that the luciferase activity of G allele was significantly higher than that of A allele in g.10751G > A (p < 0.05). However, there were no significant differences between different alleles of others two SNPs in luciferase activity. Overall, KLF3 is an important candidate gene that affects chicken abdominal fat content, and the g.10751G > A is a functional variant that potential would be applied to marker-assisted selection (MAS) for selective breeding programme.

A potential role of fatty acid binding protein 4 in the pathophysiology of autism spectrum disorder.

Autism spectrum disorder is a neurodevelopmental disorder characterized by difficulties in social communication and interaction, as well as repetitive and characteristic patterns of behaviour. Although the pathogenesis of autism spectrum disorder is unknown, being overweight or obesity during infancy and low weight at birth are known as risks, suggesting a metabolic aspect. In this study, we investigated adipose tissue development as a pathophysiological factor of autism spectrum disorder by examining the serum levels of adipokines and other metabolic markers in autism spectrum disorder children (n = 123) and typically developing children (n = 92) at 4-12 years of age. Among multiple measures exhibiting age-dependent trajectories, the leptin levels displayed different trajectory patterns between autism spectrum disorder and typically developing children, supporting an adipose tissue-dependent mechanism of autism spectrum disorder. Of particular interest, the levels of fatty acid binding protein 4 (FABP4) were significantly lower in autism spectrum disorder children than in typically developing subjects, at preschool age (4-6 years old: n = 21 for autism spectrum disorder and n = 26 for typically developing). The receiver operating characteristic curve analysis discriminated autism spectrum disorder children from typically developing children with a sensitivity of 94.4% and a specificity of 75.0%. We re-sequenced the exons of the FABP4 gene in a Japanese cohort comprising 659 autism spectrum disorder and 1000 control samples, and identified two rare functional variants in the autism spectrum disorder group. The Trp98Stop, one of the two variants, was transmitted to the proband from his mother with a history of depression. The disruption of the Fabp4 gene in mice evoked autism spectrum disorder-like behavioural phenotypes and increased spine density on apical dendrites of pyramidal neurons, which has been observed in the postmortem brains of autism spectrum disorder subjects. The Fabp4 knockout mice had an altered fatty acid composition in the cortex. Collectively, these results suggest that an 'adipo-brain axis' may underlie the pathophysiology of autism spectrum disorder, with FABP4 as a potential molecule for use as a biomarker.

Genome-wide Analyses Identify a Novel Risk Locus for Nonsyndromic Cleft Palate.

The 3 major subphenotypes observed in patients with nonsyndromic orofacial clefts (NSOFCs) are nonsyndromic cleft lip only (NSCLO), nonsyndromic cleft lip with palate (NSCLP), and nonsyndromic cleft palate only (NSCPO). However, the genetic architecture underlying NSCPO is largely unknown. Here we performed a 2-stage genome-wide association study (GWAS) on NSCPO and replication analyses of selected variants in other NSOFCs from the Chinese Han population. We identified a novel locus (15q24.3) and a known locus (1q32.2) where variants in or near the gene reached genome-wide significance (2.80 × 10-13 < P < 1.72 × 10-08) in a test for association with NSCPO in a case-control design. Although a variant from 15q24.3 was found to be significantly associated with both NSCPO and NSCLP, the direction of estimated effects on risk were opposite. Our functional annotation of the risk alleles within 15q24.3 coupled with previously established roles of the candidate genes within identified risk loci in periderm development, embryonic patterning, and/or regulation of cellular processes supports their involvement in palate development and the pathogenesis of cleft palate. Our study advances the understanding of the genetic basis of NSOFCs and provides novel insights into the pathogenesis of NSCPO.

Statistical considerations for the analysis of massively parallel reporter assays data.

Noncoding DNA contains gene regulatory elements that alter gene expression, and the function of these elements can be modified by genetic variation. Massively parallel reporter assays (MPRA) enable high-throughput identification and characterization of functional genetic variants, but the statistical methods to identify allelic effects in MPRA data have not been fully developed. In this study, we demonstrate how the baseline allelic imbalance in MPRA libraries can produce biased results, and we propose a novel, nonparametric, adaptive testing method that is robust to this bias. We compare the performance of this method with other commonly used methods, and we demonstrate that our novel adaptive method controls Type I error in a wide range of scenarios while maintaining excellent power. We have implemented these tests along with routines for simulating MPRA data in the Analysis Toolset for MPRA (@MPRA), an R package for the design and analyses of MPRA experiments. It is publicly available at http://github.com/redaq/atMPRA.

5' UTR variants in the quantitative trait gene Hnrnph1 support reduced 5' UTR usage and hnRNP H protein as a molecular mechanism underlying reduced methamphetamine sensitivity.

We previously identified a 210 kb region on chromosome 11 (50.37-50.58 Mb, mm10) containing two protein-coding genes (Hnrnph1, Rufy1) that was necessary for reduced methamphetamine-induced locomotor activity in C57BL/6J congenic mice harboring DBA/2J polymorphisms. Gene editing of a small deletion in the first coding exon supported Hnrnph1 as a quantitative trait gene. We have since shown that Hnrnph1 mutants also exhibit reduced methamphetamine-induced reward, reinforcement, and dopamine release. However, the quantitative trait variants (QTVs) that modulate Hnrnph1 function at the molecular level are not known. Nine single nucleotide polymorphisms and seven indels distinguish C57BL/6J from DBA/2J within Hnrnph1, including four variants within the 5' untranslated region (UTR). Here, we show that a 114 kb introgressed region containing Hnrnph1 and Rufy1 was sufficient to cause a decrease in MA-induced locomotor activity. Gene-level transcriptome analysis of striatal tissue from 114 kb congenics vs Hnrnph1 mutants identified a nearly perfect correlation of fold-change in expression for those differentially expressed genes that were common to both mouse lines, indicating functionally similar effects on the transcriptome and behavior. Exon-level analysis (including noncoding exons) revealed decreased 5' UTR usage of Hnrnph1 and immunoblot analysis identified a corresponding decrease in hnRNP H protein in 114 kb congenic mice. Molecular cloning of the Hnrnph1 5' UTR containing all four variants (but none of them individually) upstream of a reporter induced a decrease in reporter signal in both HEK293 and N2a cells, thus, identifying a set of QTVs underlying molecular regulation of Hnrnph1.

Pathway Analysis of Genes Identified through Post-GWAS to Underpin Prostate Cancer Aetiology.

Understanding the functional role of risk regions identified by genome-wide association studies (GWAS) has made considerable recent progress and is referred to as the post-GWAS era. Annotation of functional variants to the genes, including cis or trans and understanding their biological pathway/gene network enrichments, is expected to give rich dividends by elucidating the mechanisms underlying prostate cancer. To this aim, we compiled and analysed currently available post-GWAS data that is validated through further studies in prostate cancer, to investigate molecular biological pathways enriched for assigned functional genes. In total, about 100 canonical pathways were significantly, at false discovery rate (FDR)< 0.05), enriched in assigned genes using different algorithms. The results have highlighted some well-known cancer signalling pathways, antigen presentation processes and enrichment in cell growth and development gene networks, suggesting risk loci may exert their functional effect on prostate cancer by acting through multiple gene sets and pathways. Additional upstream analysis of the involved genes identified critical transcription factors such as HDAC1 and STAT5A. We also investigated the common genes between post-GWAS and three well-annotated gene expression datasets to endeavour to uncover the main genes involved in prostate cancer development/progression. Post-GWAS generated knowledge of gene networks and pathways, although continuously evolving, if analysed further and targeted appropriately, will have an important impact on clinical management of the disease.

Identifying Putative Susceptibility Genes and Evaluating Their Associations with Somatic Mutations in Human Cancers.

Genome-wide association studies (GWASs) have identified hundreds of genetic risk variants for human cancers. However, target genes for the majority of risk loci remain largely unexplored. It is also unclear whether GWAS risk-loci-associated genes contribute to mutational signatures and tumor mutational burden (TMB) in cancer tissues. We systematically conducted cis-expression quantitative trait loci (cis-eQTL) analyses for 294 GWAS-identified variants for six major types of cancer-colorectal, lung, ovary, prostate, pancreas, and melanoma-by using transcriptome data from the Genotype-Tissue Expression (GTEx) Project, the Cancer Genome Atlas (TCGA), and other public data sources. By using integrative analysis strategies, we identified 270 candidate target genes, including 99 with previously unreported associations, for six cancer types. By analyzing functional genomic data, our results indicate that 180 genes (66.7% of 270) had evidence of cis-regulation by putative functional variants via proximal promoter or distal enhancer-promoter interactions. Together with our previously reported associations for breast cancer risk, our results show that 24 genes are shared by at least two cancer types, including four genes for both breast and ovarian cancer. By integrating mutation data from TCGA, we found that expression levels of 33 and 66 putative susceptibility genes were associated with specific mutational signatures and TMB of cancer-driver genes, respectively, at a Bonferroni-corrected p < 0.05. Together, these findings provide further insight into our understanding of how genetic risk variants might contribute to carcinogenesis through the regulation of susceptibility genes that are related to the biogenesis of somatic mutations.