Genetic association studies of alterations in protein function expose recessive effects on cancer predisposition

Genome wide association study and meta-analysis identified multiple new risk loci for freckles in 4813 Chinese individuals

Freckle is a prevalent pigmentary dermatosis with an obvious hereditary component. Dozens of freckles risk loci have been discovered through research on multiple traits or other diseases, rather than as an independent trait. To discover novel variants associated with freckles, we performed GWAS and meta-analysis in 4813 Chinese individuals. We conducted GWAS and meta-analysis of two cohorts: 197 patients and 1603 controls (Cohort I), and 336 patients and 2677 controls (Cohort II), both from China. Then we performed linkage disequilibrium (LD) analysis, eQTL study, and enrichment analysis with association results for functional implications. Finally, we discovered 59 new SNPs and 13 novel susceptibility genes associated with freckles (Pmeta <5 × 10-8), which has enriched the genetic research on freckles.

An Inherited Allele Confers Prostate Cancer Progression and Drug Resistance via RFX6/HOXA10-Orchestrated TGFβ Signaling

Genetic and epigenetic alterations are cancer hallmark characteristics. However, the role of inherited cancer predisposition alleles in co-opting lineage factor epigenetic reprogramming and tumor progression remains elusive. Here the FinnGen cohort phenome-wide analysis, along with multiple genome-wide association studies, has consistently identified the rs339331-RFX6/6q22 locus associated with prostate cancer (PCa) risk across diverse populations. It is uncovered that rs339331 resides in a reprogrammed androgen receptor (AR) binding site in PCa tumors, with the T risk allele enhancing AR chromatin occupancy. RFX6, an AR-regulated gene linked to rs339331, exhibits synergistic prognostic value for PCa recurrence and metastasis. This comprehensive in vitro and in vivo studies demonstrate the oncogenic functions of RFX6 in promoting PCa cell proliferation and metastasis. Mechanistically, RFX6 upregulates HOXA10 that profoundly correlates with adverse PCa outcomes and is pivotal in RFX6-mediated PCa progression, facilitating the epithelial-mesenchymal transition (EMT) and modulating the TGFβ/SMAD signaling axis. Clinically, HOXA10 elevation is associated with increased EMT scores, tumor advancement and PCa recurrence. Remarkably, reducing RFX6 expression restores enzalutamide sensitivity in resistant PCa cells and tumors. This findings reveal a complex interplay of genetic and epigenetic mechanisms in PCa pathogenesis and drug resistance, centered around disrupted prostate lineage AR signaling and abnormal RFX6 expression.

Revealing the genetic complexity of hypothyroidism: integrating complementary association methods

Hypothyroidism is a common endocrine disorder whose prevalence increases with age. The disease manifests itself when the thyroid gland fails to produce sufficient thyroid hormones. The disorder includes cases of congenital hypothyroidism (CH), but most cases exhibit hormonal feedback dysregulation and destruction of the thyroid gland by autoantibodies. In this study, we sought to identify causal genes for hypothyroidism in large populations. The study used the UK-Biobank (UKB) database, reporting on 13,687 cases of European ancestry. We used GWAS compilation from Open Targets (OT) and tuned protocols focusing on genes and coding regions, along with complementary association methods of PWAS (proteome-based) and TWAS (transcriptome-based). Comparing summary statistics from numerous GWAS revealed a limited number of variants associated with thyroid development. The proteome-wide association study method identified 77 statistically significant genes, half of which are located within the Chr6-MHC locus and are enriched with autoimmunity-related genes. While coding GWAS and PWAS highlighted the centrality of immune-related genes, OT and transcriptome-wide association study mostly identified genes involved in thyroid developmental programs. We used independent populations from Finland (FinnGen) and the Taiwan cohort to validate the PWAS results. The higher prevalence in females relative to males is substantiated as the polygenic risk score prediction of hypothyroidism relied mostly from the female group genetics. Comparing results from OT, TWAS, and PWAS revealed the complementary facets of hypothyroidism's etiology. This study underscores the significance of synthesizing gene-phenotype association methods for this common, intricate disease. We propose that the integration of established association methods enhances interpretability and clinical utility.

Common and novel haplotype structures between different types of cancer

BACKGROUND

Background: Genome-wide association studies (GWAS) have identified hundreds of genetic variants associated with cancer risk. GWAS data are important for cancer prevention and understanding the underlying mechanisms of cancer.

AIMS

This study aimed to investigate the genetic association between different types of cancer using GWAS data and a bioinformatics approach.

METHODS AND RESULTS

The significant GWAS variants associated with more than one cancer type were identified. Common linkage disequilibrium (LD) variants between different types of cancer were identified by 1000 genomes phase 3 LD data. Haplotype blocks were identified by analyzing 1000 Genomes phase 3 genotyping data in the GWAS populations. Subsequent analyses included functional SNP analyses and TCGA gene expression. The results associated with significant GWAS variants (P<5E-8) showed the following haplotype associations in European population: GT rs4808075-rs8170 haplotype on BABAM1 with breast and ovarian cancers, GC rs16857609-rs11693806 haplotype on DIRC3 with breast and thyroid cancers, GCG rs380286-rs401681-rs31487 haplotype on CLPTM1L with skin and lung cancers, GGG rs4430796-rs11651052-rs11263763 haplotype on HNF1B with prostate and endometrial cancers, and GT rs10505477-rs6983267 haplotype on CASC8 associated with colorectal and prostate cancers. All these genes had significantly different expressions in tumor tissues (P<1E-3). In addition, the rs11693806 variant is located in the hsa-miR-873-5p binding site and has an enhancing effect on the hsa-miR-873-5p:DIRC3 interaction.

CONCLUSION

These novel haplotype structures and miRNA:lncRNA interactions are important for understanding the common genetic link between cancers. These results can potentially be used in genetic panels.

Novel genetic association between obesity, colorectal cancer, and inflammatory bowel disease

Purpose

Obesity/overweight is an important risk factor for CRC and IBD. The aim of this study was to investigate the role of common genetic factors and haplotypes associated with obesity, CRC and IBD.

Methods

Significant GWAS variants associated with CRC, IBD or obesity were extracted from the GWAS catalog. The common variants between CRC-IBD, CRC-obesity or IBD-obesity were identified. Finally, the haplotypic structure between these diseases was identified, and SNP function analysis, gene-gene expression, protein-protein interactions, gene survival analysis and pathway analysis were performed with the results.

Results

While the results showed several common variants between CRC and IBD, IBD and obesity, and CRC and obesity identified in previous GWAS, rs3184504 was the only common variant for CRC-IBD-obesity (P ≤ 5E-8). The result also identified a haplotypic block AGCAGT (r2 ≥ 0.8 and D'≥0.08) associated with the common variants of CRC-IBD-obesity. These variants are located on the SH2B3 gene, whose expression level decreases in both colon and rectal cancers (P ≤ 1E-3) and which has protein-protein interaction with inflammation- and cancer-associated genes.

Conclusion

The rs3184504 variant and the novel haplotype AGCAGT co-occurred in CRC, IBD, obesity, and inflammation. This novel haplotype could potentially be used in genetic panels to identify CRC/IBD susceptibility in obese patients.

rs2736098, a synonymous polymorphism, is associated with carcinogenesis and cell count in multiple tissue types by regulating TERT expression

rs2736098 is a synonymous polymorphism in TERT (telomerase reverse transcriptase), an enzyme involved in tumor onset of multiple tissues, and should play no roles in carcinogenesis. However, a search in cancer somatic mutation database indicated that the mutation frequency at rs2736098 is much higher than the average one for TERT. Moreover, there are significant H3K4me1 and H3K27Ac signals, two universal histone modifications for active enhancers, surrounding rs2736098. Therefore, we hypothesized that rs2736098 might be within an enhancer region, regulate TERT expression and influence cancer risk. Through luciferase assay, it was verified that the enhancer activity of rs2736098C allele is significantly higher than that of T in multiple tissues. Transfection of plasmids containing TERT coding region with two different alleles indicated that rs2736098C allele can induce a significantly higher TERT expression than T. By chromatin immunoprecipitation, it was observed that the fragment spanning rs2736098 can interact with USF1 (upstream transcription factor 1). The two alleles of rs2736098 present evidently different binding affinity with nuclear proteins. Database and literature search indicated that rs2736098 is significantly associated with carcinogenesis in multiple tissues and count of multiple cell types. All these facts indicated that rs2736098 is also an oncogenic polymorphism and plays important role in cell proliferation.

Epigenome-wide association study identifies neonatal DNA methylation associated with two-year attention problems in children born very preterm

Prior research has identified epigenetic predictors of attention problems in school-aged children but has not yet investigated these in young children, or children at elevated risk of attention problems due to preterm birth. The current study evaluated epigenome-wide associations between neonatal DNA methylation and attention problems at age 2 years in children born very preterm. Participants included 441 children from the Neonatal Neurobehavior and Outcomes in Very Preterm Infants (NOVI) Study, a multi-site study of infants born < 30 weeks gestational age. DNA methylation was measured from buccal swabs collected at NICU discharge using the Illumina MethylationEPIC Bead Array. Attention problems were assessed at 2 years of adjusted age using the attention problems subscale of the Child Behavior Checklist (CBCL). After adjustment for multiple testing, DNA methylation at 33 CpG sites was associated with child attention problems. Differentially methylated CpG sites were located in genes previously linked to physical and mental health, including several genes associated with ADHD in prior epigenome-wide and genome-wide association studies. Several CpG sites were located in genes previously linked to exposure to prenatal risk factors in the NOVI sample. Neonatal epigenetics measured at NICU discharge could be useful in identifying preterm children at risk for long-term attention problems and related psychiatric disorders, who could benefit from early prevention and intervention efforts.

A mammalian-specific Alex3/Gαq protein complex regulates mitochondrial trafficking, dendritic complexity, and neuronal survival

Mitochondrial dynamics and trafficking are essential to provide the energy required for neurotransmission and neural activity. We investigated how G protein-coupled receptors (GPCRs) and G proteins control mitochondrial dynamics and trafficking. The activation of Gαq inhibited mitochondrial trafficking in neurons through a mechanism that was independent of the canonical downstream PLCβ pathway. Mitoproteome analysis revealed that Gαq interacted with the Eutherian-specific mitochondrial protein armadillo repeat-containing X-linked protein 3 (Alex3) and the Miro1/Trak2 complex, which acts as an adaptor for motor proteins involved in mitochondrial trafficking along dendrites and axons. By generating a CNS-specific Alex3 knockout mouse line, we demonstrated that Alex3 was required for the effects of Gαq on mitochondrial trafficking and dendritic growth in neurons. Alex3-deficient mice had altered amounts of ER stress response proteins, increased neuronal death, motor neuron loss, and severe motor deficits. These data revealed a mammalian-specific Alex3/Gαq mitochondrial complex, which enables control of mitochondrial trafficking and neuronal death by GPCRs.

Investigation of common genetic risk factors between thyroid traits and breast cancer

Breast cancer (BC) risk is suspected to be linked to thyroid disorders, however observational studies exploring the association between BC and thyroid disorders gave conflicting results. We proposed an alternative approach by investigating the shared genetic risk factors between BC and several thyroid traits. We report a positive genetic correlation between BC and thyroxine (FT4) levels (corr = 0.13, p-value = 2.0 × 10-4) and a negative genetic correlation between BC and thyroid-stimulating hormone (TSH) levels (corr = -0.09, p-value = 0.03). These associations are more striking when restricting the analysis to estrogen receptor-positive BC. Moreover, the polygenic risk scores (PRS) for FT4 and hyperthyroidism are positively associated to BC risk (OR = 1.07, 95%CI: 1.00-1.13, p-value = 2.8 × 10-2 and OR = 1.04, 95%CI: 1.00-1.08, p-value = 3.8 × 10-2, respectively), while the PRS for TSH is inversely associated to BC risk (OR = 0.93, 95%CI: 0.89-0.97, p-value = 2.0 × 10-3). Using the PLACO method, we detected 49 loci associated to both BC and thyroid traits (p-value < 5 × 10-8), in the vicinity of 130 genes. An additional colocalization and gene-set enrichment analyses showed a convincing causal role for a known pleiotropic locus at 2q35 and revealed an additional one at 8q22.1 associated to both BC and thyroid cancer. We also found two new pleiotropic loci at 14q32.33 and 17q21.31 that were associated to both TSH levels and BC risk. Enrichment analyses and evidence of regulatory signals also highlighted brain tissues and immune system as candidates for obtaining associations between BC and TSH levels. Overall, our study sheds light on the complex interplay between BC and thyroid traits and provides evidence of shared genetic risk between those conditions.

Genome-wide prediction of disease variant effects with a deep protein language model

Predicting the effects of coding variants is a major challenge. While recent deep-learning models have improved variant effect prediction accuracy, they cannot analyze all coding variants due to dependency on close homologs or software limitations. Here we developed a workflow using ESM1b, a 650-million-parameter protein language model, to predict all ~450 million possible missense variant effects in the human genome, and made all predictions available on a web portal. ESM1b outperformed existing methods in classifying ~150,000 ClinVar/HGMD missense variants as pathogenic or benign and predicting measurements across 28 deep mutational scan datasets. We further annotated ~2 million variants as damaging only in specific protein isoforms, demonstrating the importance of considering all isoforms when predicting variant effects. Our approach also generalizes to more complex coding variants such as in-frame indels and stop-gains. Together, these results establish protein language models as an effective, accurate and general approach to predicting variant effects.

Explainable multi-task learning improves the parallel estimation of polygenic risk scores for many diseases through shared genetic basis

Many complex diseases share common genetic determinants and are comorbid in a population. We hypothesized that the co-occurrences of diseases and their overlapping genetic etiology can be exploited to simultaneously improve multiple diseases' polygenic risk scores (PRS). This hypothesis was tested using a multi-task learning (MTL) approach based on an explainable neural network architecture. We found that parallel estimations of the PRS for 17 prevalent cancers in a pan-cancer MTL model were generally more accurate than independent estimations for individual cancers in comparable single-task learning (STL) models. Such performance improvement conferred by positive transfer learning was also observed consistently for 60 prevalent non-cancer diseases in a pan-disease MTL model. Interpretation of the MTL models revealed significant genetic correlations between the important sets of single nucleotide polymorphisms used by the neural network for PRS estimation. This suggested a well-connected network of diseases with shared genetic basis.

Gene-based association study reveals a distinct female genetic signal in primary hypertension

Hypertension is a polygenic disease that affects over 1.2 billion adults aged 30-79 worldwide. It is a major risk factor for renal, cerebrovascular, and cardiovascular diseases. The heritability of hypertension is estimated to be high; nevertheless, our understanding of its underlying mechanisms remains scarce and incomplete. This study covered the entries from European ancestry from the UK-Biobank (UKB), with 74,090 cases diagnosed with essential (primary) hypertension and 200,734 controls. We compared the findings from large-scale genome-wide association studies (GWAS) to the gene-based method of proteome-wide association studies (PWAS). We focused on 70 statistically significant associated genes, most of which failed to reach significance in variant-based GWAS. A total of 30% of the PWAS-associated genes were validated against independent cohorts, including the Finnish Biobank. Furthermore, gene-based analyses that were performed on both sexes revealed sex-dependent genetics with a stronger genetic component associated with females. Analysis of systolic and diastolic blood pressure measurements confirms a strong genetic effect associated with females. We demonstrated that gene-based approaches provide insight into the underlying biology of hypertension. Specifically, the expression profiles of the identified genes exposed the enrichment of endothelial cells from multiple organs. Furthermore, females' top-ranked significant genes are involved in cellular immunity. We conclude that studying hypertension and blood pressure via gene-based association methods improves interpretability and exposes sex-dependent genetic effects, which enhances clinical utility.

The performance of multi-gene panels for breast/ovarian cancer predisposition

BRCA1 and BRCA2 are the most mutated genes in breast cancer. We analyzed 48 breast cancer subjects using two methods that differ in terms of number of genes investigated and strategy used (primers: Panel A - 12 genes - vs probes: Panel B - 48 genes). Both the panels and procedures identified "pathogenic" or "likely pathogenic" variants in TP53, ATM, CHEK2 and BARD1 besides BRCA1 and BRCA2. Panel B identified two other putatively pathogenic variants in RNASEL and in RAD50. Identification of variants other than the BRCA genes can be useful in patient management. A total of 121 variants were distributed within the 12 genes and were correctly detected by both panels. However, the number of calls without divergence, namely ± 0.10 difference of allelic frequency, was 78.3%, while calls with a divergence below 0.10 was 16.7%, thus indicating that only 5% (n = 275) of 5,412 calls had a divergence above 0.10. Although these panels differ from each other, both are useful in different situations, particularly when patients should be tested for genes other than BRCA1/2 (as occurs in patients affected by a so called hereditary syndrome) or for therapeutic purposes.

Long-Term Simulation of Microgravity Induces Changes in Gene Expression in Breast Cancer Cells

Microgravity changes the gene expression pattern in various cell types. This study focuses on the breast cancer cell lines MCF-7 (less invasive) and MDA-MB-231 (triple-negative, highly invasive). The cells were cultured for 14 days under simulated microgravity (s-µg) conditions using a random positioning machine (RPM). We investigated cytoskeletal and extracellular matrix (ECM) factors as well as focal adhesion (FA) and the transmembrane proteins involved in different cellular signaling pathways (MAPK, PAM and VEGF). The mRNA expressions of 24 genes of interest (TUBB, ACTB, COL1A1, COL4A5, LAMA3, ITGB1, CD44, VEGF, FLK1, EGFR, SRC, FAK1, RAF1, AKT1, ERK1, MAPK14, MAP2K1, MTOR, RICTOR, VCL, PXN, CDKN1, CTNNA1 and CTNNB1) were determined by quantitative real-time PCR (qPCR) and studied using STRING interaction analysis. Histochemical staining was carried out to investigate the morphology of the adherent cells (ADs) and the multicellular spheroids (MCSs) after RPM exposure. To better understand this experimental model in the context of breast cancer patients, a weighted gene co-expression network analysis (WGCNA) was conducted to obtain the expression profiles of 35 breast cell lines from the HMS LINCS Database. The qPCR-verified genes were searched in the mammalian phenotype database and the human genome-wide association studies (GWAS) Catalog. The results demonstrated the positive association between the real metastatic microtumor environment and MCSs with respect to the extracellular matrix, cytoskeleton, morphology, different cellular signaling pathway key proteins and several other components. In summary, the microgravity-engineered three-dimensional MCS model can be utilized to study breast cancer cell behavior and to assess the therapeutic efficacies of drugs against breast cancer in the future.

Parent-of-Origin inference for biobanks

Identical genetic variations can have different phenotypic effects depending on their parent of origin. Yet, studies focusing on parent-of-origin effects have been limited in terms of sample size due to the lack of parental genomes or known genealogies. We propose a probabilistic approach to infer the parent-of-origin of individual alleles that does not require parental genomes nor prior knowledge of genealogy. Our model uses Identity-By-Descent sharing with second- and third-degree relatives to assign alleles to parental groups and leverages chromosome X data in males to distinguish maternal from paternal groups. We combine this with robust haplotype inference and haploid imputation to infer the parent-of-origin for 26,393 UK Biobank individuals. We screen 99 phenotypes for parent-of-origin effects and replicate the discoveries of 6 GWAS studies, confirming signals on body mass index, type 2 diabetes, standing height and multiple blood biomarkers, including the known maternal effect at the MEG3/DLK1 locus on platelet phenotypes. We also report a novel maternal effect at the TERT gene on telomere length, thereby providing new insights on the heritability of this phenotype. All our summary statistics are publicly available to help the community to better characterize the molecular mechanisms leading to parent-of-origin effects and their implications for human health.

Open problems in human trait genetics

Genetic studies of human traits have revolutionized our understanding of the variation between individuals, and yet, the genetics of most traits is still poorly understood. In this review, we highlight the major open problems that need to be solved, and by discussing these challenges provide a primer to the field. We cover general issues such as population structure, epistasis and gene-environment interactions, data-related issues such as ancestry diversity and rare genetic variants, and specific challenges related to heritability estimates, genetic association studies, and polygenic risk scores. We emphasize the interconnectedness of these problems and suggest promising avenues to address them.