Participation bias in the UK Biobank distorts genetic associations and downstream analyses

While volunteer-based studies such as the UK Biobank have become the cornerstone of genetic epidemiology, the participating individuals are rarely representative of their target population. To evaluate the impact of selective participation, here we derived UK Biobank participation probabilities on the basis of 14 variables harmonized across the UK Biobank and a representative sample. We then conducted weighted genome-wide association analyses on 19 traits. Comparing the output from weighted genome-wide association analyses (neffective = 94,643 to 102,215) with that from standard genome-wide association analyses (n = 263,464 to 283,749), we found that increasing representativeness led to changes in SNP effect sizes and identified novel SNP associations for 12 traits. While heritability estimates were less impacted by weighting (maximum change in h2, 5%), we found substantial discrepancies for genetic correlations (maximum change in rg, 0.31) and Mendelian randomization estimates (maximum change in βSTD, 0.15) for socio-behavioural traits. We urge the field to increase representativeness in biobank samples, especially when studying genetic correlates of behaviour, lifestyles and social outcomes

Mendelian randomization integrating GWAS and eQTL data reveals genetic determinants of complex and clinical traits

Genome-wide association studies (GWAS) have identified thousands of variants associated with complex traits, but their biological interpretation often remains unclear. Most of these variants overlap with expression QTLs, indicating their potential involvement in regulation of gene expression. Here, we propose a transcriptome-wide summary statistics-based Mendelian Randomization approach (TWMR) that uses multiple SNPs as instruments and multiple gene expression traits as exposures, simultaneously. Applied to 43 human phenotypes, it uncovers 3,913 putatively causal gene-trait associations, 36% of which have no genome-wide significant SNP nearby in previous GWAS. Using independent association summary statistics, we find that the majority of these loci were missed by GWAS due to power issues. Noteworthy among these links is educational attainment-associated BSCL2, known to carry mutations leading to a Mendelian form of encephalopathy. We also find pleiotropic causal effects suggestive of mechanistic connections. TWMR better accounts for pleiotropy and has the potential to identify biological mechanisms underlying complex traits

Mendelian randomization integrating GWAS and eQTL data reveals genetic determinants of complex and clinical traits

Genome-wide association studies (GWAS) have identified thousands of variants associated with complex traits, but their biological interpretation often remains unclear. Most of these variants overlap with expression QTLs, indicating their potential involvement in regulation of gene expression. Here, we propose a transcriptome-wide summary statistics-based Mendelian Randomization approach (TWMR) that uses multiple SNPs as instruments and multiple gene expression traits as exposures, simultaneously. Applied to 43 human phenotypes, it uncovers 3,913 putatively causal gene-trait associations, 36% of which have no genome-wide significant SNP nearby in previous GWAS. Using independent association summary statistics, we find that the majority of these loci were missed by GWAS due to power issues. Noteworthy among these links is educational attainment-associated BSCL2, known to carry mutations leading to a Mendelian form of encephalopathy. We also find pleiotropic causal effects suggestive of mechanistic connections. TWMR better accounts for pleiotropy and has the potential to identify biological mechanisms underlying complex traits.Peer reviewe

Limited evidence for blood eQTLs in human sexual dimorphism

The genetic underpinning of sexual dimorphism is very poorly understood. The prevalence of many diseases differs between men and women, which could be in part caused by sex-specific genetic effects. Nevertheless, only a few published genome-wide association studies (GWAS) were performed separately in each sex. The reported enrichment of expression quantitative trait loci (eQTLs) among GWAS-associated SNPs suggests a potential role of sex-specific eQTLs in the sex-specific genetic mechanism underlying complex traits. To explore this scenario, we combined sex-specific whole blood RNA-seq eQTL data from 3447 European individuals included in BIOS Consortium and GWAS data from UK Biobank. Next, to test the presence of sex-biased causal effect of gene expression on complex traits, we performed sex-specific transcriptome-wide Mendelian randomization (TWMR) analyses on the two most sexually dimorphic traits, waist-to-hip ratio (WHR) and testosterone levels. Finally, we performed power analysis to calculate the GWAS sample size needed to observe sex-specific trait associations driven by sex-biased eQTLs. Among 9 million SNP-gene pairs showing sex-combined associations, we found 18 genes with significant sex-biased cis-eQTLs (FDR 5%). Our phenome-wide association study of the 18 top sex-biased eQTLs on >700 traits unraveled that these eQTLs do not systematically translate into detectable sex-biased trait-associations. In addition, we observed that sex-specific causal effects of gene expression on complex traits are not driven by sex-specific eQTLs. Power analyses using real eQTL- and causal-effect sizes showed that millions of samples would be necessary to observe sex-biased trait associations that are fully driven by sex-biased cis-eQTLs. Compensatory effects may further hamper their detection. Our results suggest that sex-specific eQTLs in whole blood do not translate to detectable sex-specific trait associations of complex diseases, and vice versa that the observed sex-specific trait associations cannot be explained by sex-specific eQTLs

Successful Pregnancies, Births, and Children Development Following Oocyte Cryostorage in Female Cancer Patients During 25 Years of Fertility Preservation

Simple Summary The study goal is to demonstrate that oocyte cryopreservation is a feasible and efficient option for fertility preservation in cancer patients through the comparison of in vitro fertilization treatments in nononcological patients. The preservation of fertility in cancer patients is a crucial aspect of modern reproductive medicine. Amenorrhea and infertility often occur after cancer therapy, worsening the quality of life. Cryopreservation of oocytes in young cancer patients is a therapeutic option for preserving fertility. A prospective study was conducted on 508 cancer patients who underwent oocyte cryopreservation to preserve fertility between 1996 and 2021 including the COVID-19 pandemic period. Patients underwent ovarian stimulation, followed by egg retrieval, and oocytes were cryopreserved by slow freezing or vitrification. Sixty-four thawing/warming cycles were performed. Survival, fertilization, pregnancy, and birth rate over the thawing/warming cycles were obtained. The data were compared with those from a group of 1042 nononcological patients who cryopreserved supernumerary oocytes. An average of 8.8 +/- 6.9 oocytes were retrieved per cycle, and 6.1 +/- 4.2 oocytes were cryopreserved. With their own stored oocytes, 44 patients returned to attempt pregnancy. From a total of 194 thawed/warmed oocytes, 157 survived (80%). In total, 100 embryos were transferred in 57 transfer/cycles, and 18 pregnancies were achieved. The pregnancy rate per transfer and pregnancy rate per patient were 31% and 41%, respectively. No statistically significant differences were observed between oncological patients and nononcological patients. A total of 15 babies were born from oncological patients. Children born showed normal growth and development. One minor malformation was detected

CXCL1-CXCR1/2 signaling is induced in human temporal lobe epilepsy and contributes to seizures in a murine model of acquired epilepsy

Abstract CXCL1, a functional murine orthologue of the human chemokine CXCL8 (IL-8), and its CXCR1 and CXCR2 receptors were investigated in a murine model of acquired epilepsy developing following status epilepticus (SE) induced by intra-amygdala kainate. CXCL8 and its receptors were also studied in human temporal lobe epilepsy (TLE). The functional involvement of the chemokine in seizure generation and neuronal cell loss was assessed in mice using reparixin (formerly referred to as repertaxin), a non-competitive allosteric inhibitor of CXCR1/2 receptors. We found a significant increase in hippocampal CXCL1 level within 24 h of SE onset that lasted for at least 1 week. No changes were measured in blood. In analogy with human TLE, immunohistochemistry in epileptic mice showed that CXCL1 and its two receptors were increased in hippocampal neuronal cells. Additional expression of these molecules was found in glia in human TLE. Mice were treated with reparixin or vehicle during SE and for additional 6 days thereafter, using subcutaneous osmotic minipumps. Drug-treated mice showed a faster SE decay, a reduced incidence of acute symptomatic seizures during 48 h post-SE, and a delayed time to spontaneous seizures onset compared to vehicle controls. Upon reparixin discontinuation, mice developed spontaneous seizures similar to vehicle mice, as shown by EEG monitoring at 14 days and 2.5 months post-SE. In the same epileptic mice, reparixin reduced neuronal cell loss in the hippocampus vs vehicle-injected mice, as assessed by Nissl staining at completion of EEG monitoring. Reparixin administration for 2 weeks in mice with established chronic seizures, reduced by 2-fold on average seizure number vs pre-treatment baseline, and this effect was reversible upon drug discontinuation. No significant changes in seizure number were measured in vehicle-injected epileptic mice that were EEG monitored in parallel. Data show that CXCL1-IL-8 signaling is activated in experimental and human epilepsy and contributes to acute and chronic seizures in mice, therefore representing a potential new target to attain anti-ictogenic effects

Cognitive Syndromes and C9orf72 Mutation Are Not Related to Cerebellar Degeneration in Amyotrophic Lateral Sclerosis

ObjectiveThe notion that cerebellar pathology may contribute to cognitive impairment in ALS, especially in patients with C9orf72 repeated expansion, has been inconsistently reported. This study aimed exploring the relationship between cerebellar involvement, cognitive impairment and C9orf72 repeated expansion of patients with ALS.MethodsQuantitative in vivo assessment of cerebellar lobules has been investigated in 66 non-demented patients with ALS and 28 healthy controls (HCs). Pathologic C9orf72 repeated expansion was found in 13 patients. Mild cognitive and/or behavioral impairment was diagnosed in 22 C9orf72 negative ALS patients. Measures of cortical volume (CV) and cortical thickness (CT) of cerebellar lobules of all participants were used for Principal Component Analysis (PCA) to identify clusters of lobular measures highly correlated with each other. PCA outcomes were used for between group comparisons and correlation analyses with neuropsychological and clinical features.ResultsDisease severity measured with ALS functional rating scale and index of disease progression rate significantly correlated with CV reduction of the second PCA cluster loading CV measures of anterior lobules. In all patients, cognitive impairment, measured with verbal fluency, was related to CV reduction of the third cluster comprising posterior lobules. No specific cortical thinning or volume reduction of cerebellar clustering patterns could be detected in ALS subgroups.ConclusionOur data show that specific patterns of subregional cerebellar involvement are associated with physical disability or cognitive impairment in ALS, in line with the topographic organization of the cerebellum. However, there was no specific correlation between cerebellar degeneration and cognitive syndromes or C9orf72 mutations

Genome-wide association study of susceptibility loci for breast cancer in Sardinian population

Abstract Background Despite progress in identifying genes associated with breast cancer, many more risk loci exist. Genome-wide association analyses in genetically-homogeneous populations, such as that of Sardinia (Italy), could represent an additional approach to detect low penetrance alleles. Methods We performed a genome-wide association study comparing 1431 Sardinian patients with non-familial, BRCA1/2-mutation-negative breast cancer to 2171 healthy Sardinian blood donors. DNA was genotyped using GeneChip Human Mapping 500 K Arrays or Genome-Wide Human SNP Arrays 6.0. To increase genomic coverage, genotypes of additional SNPs were imputed using data from HapMap Phase II. After quality control filtering of genotype data, 1367 cases (9 men) and 1658 controls (1156 men) were analyzed on a total of 2,067,645 SNPs. Results Overall, 33 genomic regions (67 candidate SNPs) were associated with breast cancer risk at the p < 10−6 level. Twenty of these regions contained defined genes, including one already associated with breast cancer risk: TOX3. With a lower threshold for preliminary significance to p < 10−5, we identified 11 additional SNPs in FGFR2, a well-established breast cancer-associated gene. Ten candidate SNPs were selected, excluding those already associated with breast cancer, for technical validation as well as replication in 1668 samples from the same population. Only SNP rs345299, located in intron 1 of VAV3, remained suggestively associated (p-value, 1.16x10−5), but it did not associate with breast cancer risk in pooled data from two large, mixed-population cohorts. Conclusions This study indicated the role of TOX3 and FGFR2 as breast cancer susceptibility genes in BRCA1/2-wild-type breast cancer patients from Sardinian population

Genome-wide association study of susceptibility loci for breast cancer in Sardinian population.

BACKGROUND: Despite progress in identifying genes associated with breast cancer, many more risk loci exist. Genome-wide association analyses in genetically-homogeneous populations, such as that of Sardinia (Italy), could represent an additional approach to detect low penetrance alleles. METHODS: We performed a genome-wide association study comparing 1431 Sardinian patients with non-familial, BRCA1/2-mutation-negative breast cancer to 2171 healthy Sardinian blood donors. DNA was genotyped using GeneChip Human Mapping 500 K Arrays or Genome-Wide Human SNP Arrays 6.0. To increase genomic coverage, genotypes of additional SNPs were imputed using data from HapMap Phase II. After quality control filtering of genotype data, 1367 cases (9 men) and 1658 controls (1156 men) were analyzed on a total of 2,067,645 SNPs. RESULTS: Overall, 33 genomic regions (67 candidate SNPs) were associated with breast cancer risk at the p <  0(-6) level. Twenty of these regions contained defined genes, including one already associated with breast cancer risk: TOX3. With a lower threshold for preliminary significance to p < 10(-5), we identified 11 additional SNPs in FGFR2, a well-established breast cancer-associated gene. Ten candidate SNPs were selected, excluding those already associated with breast cancer, for technical validation as well as replication in 1668 samples from the same population. Only SNP rs345299, located in intron 1 of VAV3, remained suggestively associated (p-value, 1.16 x 10(-5)), but it did not associate with breast cancer risk in pooled data from two large, mixed-population cohorts. CONCLUSIONS: This study indicated the role of TOX3 and FGFR2 as breast cancer susceptibility genes in BRCA1/2-wild-type breast cancer patients from Sardinian population

Parent-of-origin-specific allelic associations among 106 genomic loci for age at menarche.

Age at menarche is a marker of timing of puberty in females. It varies widely between individuals, is a heritable trait and is associated with risks for obesity, type 2 diabetes, cardiovascular disease, breast cancer and all-cause mortality. Studies of rare human disorders of puberty and animal models point to a complex hypothalamic-pituitary-hormonal regulation, but the mechanisms that determine pubertal timing and underlie its links to disease risk remain unclear. Here, using genome-wide and custom-genotyping arrays in up to 182,416 women of European descent from 57 studies, we found robust evidence (P < 5 × 10(-8)) for 123 signals at 106 genomic loci associated with age at menarche. Many loci were associated with other pubertal traits in both sexes, and there was substantial overlap with genes implicated in body mass index and various diseases, including rare disorders of puberty. Menarche signals were enriched in imprinted regions, with three loci (DLK1-WDR25, MKRN3-MAGEL2 and KCNK9) demonstrating parent-of-origin-specific associations concordant with known parental expression patterns. Pathway analyses implicated nuclear hormone receptors, particularly retinoic acid and γ-aminobutyric acid-B2 receptor signalling, among novel mechanisms that regulate pubertal timing in humans. Our findings suggest a genetic architecture involving at least hundreds of common variants in the coordinated timing of the pubertal transition