Canadian COVID-19 host genetics cohort replicates known severity associations

The HostSeq initiative recruited 10,059 Canadians infected with SARS-CoV-2 between March 2020 and March 2023, obtained clinical information on their disease experience and whole genome sequenced (WGS) their DNA. We analyzed the WGS data for genetic contributors to severe COVID-19 (considering 3,499 hospitalized cases and 4,975 non-hospitalized after quality control). We investigated the evidence for replication of loci reported by the International Host Genetics Initiative (HGI); analyzed the X chromosome; conducted rare variant gene-based analysis and polygenic risk score testing. Population stratification was adjusted for using meta-analysis across ancestry groups. We replicated two loci identified by the HGI for COVID-19 severity: the LZTFL1/SLC6A20 locus on chromosome 3 and the FOXP4 locus on chromosome 6 (the latter with a variant significant at P < 5E-8). We found novel significant associations with MRAS and WDR89 in gene-based analyses, and constructed a polygenic risk score that explained 1.01% of the variance in severe COVID-19. This study provides independent evidence confirming the robustness of previously identified COVID-19 severity loci by the HGI and identifies novel genes for further investigation.

Dose-dependent Association of Alcohol Consumption With Obesity and Type 2 Diabetes: Mendelian Randomization Analyses

CONTEXT

Effects of modest alcohol consumption remain controversial. Mendelian randomization (MR) can help to mitigate biases due to confounding and reverse causation in observational studies, and evaluate the potential causal role of alcohol consumption.

OBJECTIVE

This work aimed to evaluate dose-dependent effect of alcohol consumption on obesity and type 2 diabetes.

METHODS

Assessing 408 540 participants of European ancestry in the UK Biobank, we first tested the association between self-reported alcohol intake frequency and 10 anthropometric measurements, obesity, and type 2 diabetes. We then conducted MR analyses both in the overall population and in subpopulations stratified by alcohol intake frequency.

RESULTS

Among individuals having more than 14 drinks per week, a 1-drink-per-week increase in genetically predicted alcohol intake frequency was associated with a 0.36-kg increase in fat mass (SD = 0.03 kg), a 1.08-fold increased odds of obesity (95% CI, 1.06-1.10), and a 1.10-fold increased odds of type 2 diabetes (95% CI, 1.06-1.13). These associations were stronger in women than in men. Furthermore, no evidence was found supporting the association between genetically increased alcohol intake frequency and improved health outcomes among individuals having 7 or fewer drinks per week, as MR estimates largely overlapped with the null. These results withstood multiple sensitivity analyses assessing the validity of MR assumptions.

CONCLUSION

As opposed to observational associations, MR results suggest there may not be protective effects of modest alcohol consumption on obesity traits and type 2 diabetes. Heavy alcohol consumption could lead to increased measures of obesity as well as increased risk of type 2 diabetes.

A Review of the Epigenetic Clock: Emerging Biomarkers for Asthma and Allergic Disease

DNA methylation (DNAm) is a dynamic, age-dependent epigenetic modification that can be used to study interactions between genetic and environmental factors. Environmental exposures during critical periods of growth and development may alter DNAm patterns, leading to increased susceptibility to diseases such as asthma and allergies. One method to study the role of DNAm is the epigenetic clock-an algorithm that uses DNAm levels at select age-informative Cytosine-phosphate-Guanine (CpG) dinucleotides to predict epigenetic age (EA). The difference between EA and calendar age (CA) is termed epigenetic age acceleration (EAA) and reveals information about the biological capacity of an individual. Associations between EAA and disease susceptibility have been demonstrated for a variety of age-related conditions and, more recently, phenotypes such as asthma and allergic diseases, which often begin in childhood and progress throughout the lifespan. In this review, we explore different epigenetic clocks and how they have been applied, particularly as related to childhood asthma. We delve into how in utero and early life exposures (e.g., smoking, air pollution, maternal BMI) result in methylation changes. Furthermore, we explore the potential for EAA to be used as a biomarker for asthma and allergic diseases and identify areas for further study.

Converging evidence from exome sequencing and common variants implicates target genes for osteoporosis

In this study, we leveraged the combined evidence of rare coding variants and common alleles to identify therapeutic targets for osteoporosis. We undertook a large-scale multiancestry exome-wide association study for estimated bone mineral density, which showed that the burden of rare coding alleles in 19 genes was associated with estimated bone mineral density (P < 3.6 × 10-7). These genes were highly enriched for a set of known causal genes for osteoporosis (65-fold; P = 2.5 × 10-5). Exome-wide significant genes had 96-fold increased odds of being the top ranked effector gene at a given GWAS locus (P = 1.8 × 10-10). By integrating proteomics Mendelian randomization evidence, we prioritized CD109 (cluster of differentiation 109) as a gene for which heterozygous loss of function is associated with higher bone density. CRISPR-Cas9 editing of CD109 in SaOS-2 osteoblast-like cell lines showed that partial CD109 knockdown led to increased mineralization. This study demonstrates that the convergence of common and rare variants, proteomics and CRISPR can highlight new bone biology to guide therapeutic development.

Development of risk prediction models for depression combining genetic and early life risk factors

Background

Both genetic and early life risk factors play important roles in the pathogenesis and progression of adult depression. However, the interplay between these risk factors and their added value to risk prediction models have not been fully elucidated.

Methods

Leveraging a meta-analysis of major depressive disorder genome-wide association studies (N = 45,591 cases and 97,674 controls), we developed and optimized a polygenic risk score for depression using LDpred in a model selection dataset from the UK Biobank (N = 130,092 European ancestry individuals). In a UK Biobank test dataset (N = 278,730 European ancestry individuals), we tested whether the polygenic risk score and early life risk factors were associated with each other and compared their associations with depression phenotypes. Finally, we conducted joint predictive modeling to combine this polygenic risk score with early life risk factors by stepwise regression, and assessed the model performance in identifying individuals at high risk of depression.

Results

In the UK Biobank test dataset, the polygenic risk score for depression was moderately associated with multiple early life risk factors. For instance, a one standard deviation increase in the polygenic risk score was associated with 1.16-fold increased odds of frequent domestic violence (95% CI: 1.14-1.19) and 1.09-fold increased odds of not having access to medical care as a child (95% CI: 1.05-1.14). However, the polygenic risk score was more strongly associated with depression phenotypes than most early life risk factors. A joint predictive model integrating the polygenic risk score, early life risk factors, age and sex achieved an AUROC of 0.6766 for predicting strictly defined major depressive disorder, while a model without the polygenic risk score and a model without any early life risk factors had an AUROC of 0.6593 and 0.6318, respectively.

Conclusion

We have developed a polygenic risk score to partly capture the genetic liability to depression. Although genetic and early life risk factors can be correlated, joint predictive models improved risk stratification despite limited improvement in magnitude, and may be explored as tools to better identify individuals at high risk of depression.

Novel insights into systemic sclerosis using a sensitive computational method to analyze whole-genome bisulfite sequencing data

BACKGROUND

Abnormal DNA methylation is thought to contribute to the onset and progression of systemic sclerosis. Currently, the most comprehensive assay for profiling DNA methylation is whole-genome bisulfite sequencing (WGBS), but its precision depends on read depth and it may be subject to sequencing errors. SOMNiBUS, a method for regional analysis, attempts to overcome some of these limitations. Using SOMNiBUS, we re-analyzed WGBS data previously analyzed using bumphunter, an approach that initially fits single CpG associations, to contrast DNA methylation estimates by both methods.

METHODS

Purified CD4+ T lymphocytes of 9 SSc and 4 control females were sequenced using WGBS. We separated the resulting sequencing data into regions with dense CpG data, and differentially methylated regions (DMRs) were inferred with the SOMNiBUS region-level test, adjusted for age. Pathway enrichment analysis was performed with ingenuity pathway analysis (IPA). We compared the results obtained by SOMNiBUS and bumphunter.

RESULTS

Of 8268 CpG regions of ≥ 60 CpGs eligible for analysis with SOMNiBUS, we identified 131 DMRs and 125 differentially methylated genes (DMGs; p-values less than Bonferroni-corrected threshold of 6.05-06 controlling family-wise error rate at 0.05; 1.6% of the regions). In comparison, bumphunter identified 821,929 CpG regions, 599 DMRs (of which none had ≥ 60 CpGs) and 340 DMGs (q-value of 0.05; 0.04% of all regions). The top ranked gene identified by SOMNiBUS was FLT4, a lymphangiogenic orchestrator, and the top ranked gene on chromosome X was CHST7, known to catalyze the sulfation of glycosaminoglycans in the extracellular matrix. The top networks identified by IPA included connective tissue disorders.

CONCLUSIONS

SOMNiBUS is a complementary method of analyzing WGBS data that enhances biological insights into SSc and provides novel avenues of investigation into its pathogenesis.

Circulating proteins to predict COVID-19 severity

Predicting COVID-19 severity is difficult, and the biological pathways involved are not fully understood. To approach this problem, we measured 4701 circulating human protein abundances in two independent cohorts totaling 986 individuals. We then trained prediction models including protein abundances and clinical risk factors to predict COVID-19 severity in 417 subjects and tested these models in a separate cohort of 569 individuals. For severe COVID-19, a baseline model including age and sex provided an area under the receiver operator curve (AUC) of 65% in the test cohort. Selecting 92 proteins from the 4701 unique protein abundances improved the AUC to 88% in the training cohort, which remained relatively stable in the testing cohort at 86%, suggesting good generalizability. Proteins selected from different COVID-19 severity were enriched for cytokine and cytokine receptors, but more than half of the enriched pathways were not immune-related. Taken together, these findings suggest that circulating proteins measured at early stages of disease progression are reasonably accurate predictors of COVID-19 severity. Further research is needed to understand how to incorporate protein measurement into clinical care.

Genetic predisposition to depression and inflammation impacts symptom burden and survival in patients with head and neck cancer: A longitudinal study

OBJECTIVE

The primary purpose of this study was to investigate the contribution of genetic predispositions to depression and inflammation, as measured through polygenic risk scores, on symptom burden (physical and psychological) in patients with head and neck cancer in the immediate post-treatment period (i.e., at three months post-diagnosis), as well as on 3-, 6-, 12-, 24- and 36-month survival.

METHODS

Prospective longitudinal study of 223 adults (72 % participation) newly diagnosed with a first occurrence of primary head and neck cancer, paired with genetic data (Illumina PsychArray), validated psychometric measures, Structured Clinical Interviews for DSM Disorders (SCID-I), and medical chart reviews.

RESULTS

Symptom burden at 3 months was predicted by (R² adj. = 0.38, p < 0.001): a baseline SCID-I Anxiety Disorder (b = 1.69, B = 0.23, 95%CI = 0.43-2.94; p = 0.009), baseline levels of HADS anxiety (b = 0.20, B = 0.29, 95%CI = 0.07-0.34; p = 0.003), the polygenic risk score (PRS) for depression (b = 0.66, B = 0.18, 95%CI = 0.003-1.32; p = 0.049), and cumulated dose of radiotherapy (b = 0.002, B = 0.46, 95%CI = 0.001-0.003; p < 0.001). When controlling for factors known to be associated with cancer survival, patients with a higher PRS associated with depression and inflammation, respectively, presented higher risk of death within 36 months (b = 1.75, Exp(B) = 5.75, 95%CI = 1.55-21.27, p = 0.009 and b = 0.14, Exp(B) = 1.15, 95%CI = 1.01-1.30, p = 0.03).

CONCLUSIONS

Our results outline three potential pathways of symptom burden in patients with head and neck cancer: a genetic predisposition towards depression; an initial anxiety disorder upon being diagnosed with cancer or high levels of anxiety upon diagnosis; and a dose-related response to radiotherapy. One may want to investigate early interventions in these areas to alleviate symptom burden in patients faced with a life-threatening disease, as well as consider targeting genetic predisposition towards depression and inflammation implicated in survival. The high prevalence of distress in patients with head and neck cancer is an opportunity to study genetic predispositions, which could potentially be broadly generalized to other cancers and diseases.

Molecular Genetic Characteristics of FANCI, a Proposed New Ovarian Cancer Predisposing Gene

FANCI was recently identified as a new candidate ovarian cancer (OC)-predisposing gene from the genetic analysis of carriers of FANCI c.1813C>T; p.L605F in OC families. Here, we aimed to investigate the molecular genetic characteristics of FANCI, as they have not been described in the context of cancer. We first investigated the germline genetic landscape of two sisters with OC from the discovery FANCI c.1813C>T; p.L605F family (F1528) to re-affirm the plausibility of this candidate. As we did not find other conclusive candidates, we then performed a candidate gene approach to identify other candidate variants in genes involved in the FANCI protein interactome in OC families negative for pathogenic variants in BRCA1, BRCA2, BRIP1, RAD51C, RAD51D, and FANCI, which identified four candidate variants. We then investigated FANCI in high-grade serous ovarian carcinoma (HGSC) from FANCI c.1813C>T carriers and found evidence of loss of the wild-type allele in tumour DNA from some of these cases. The somatic genetic landscape of OC tumours from FANCI c.1813C>T carriers was investigated for mutations in selected genes, copy number alterations, and mutational signatures, which determined that the profiles of tumours from carriers were characteristic of features exhibited by HGSC cases. As other OC-predisposing genes such as BRCA1 and BRCA2 are known to increase the risk of other cancers including breast cancer, we investigated the carrier frequency of germline FANCI c.1813C>T in various cancer types and found overall more carriers among cancer cases compared to cancer-free controls (p = 0.007). In these different tumour types, we also identified a spectrum of somatic variants in FANCI that were not restricted to any specific region within the gene. Collectively, these findings expand on the characteristics described for OC cases carrying FANCI c.1813C>T; p.L605F and suggest the possible involvement of FANCI in other cancer types at the germline and/or somatic level.

Circulating Proteins Influencing Psychiatric Disease: A Mendelian Randomization Study

BACKGROUND

There is a pressing need for novel drug targets for psychiatric disorders. Circulating proteins are potential candidates because they are relatively easy to measure and modulate and play important roles in signaling.

METHODS

We performed two-sample Mendelian randomization analyses to estimate the associations between circulating protein abundances and risk of 10 psychiatric disorders. Genetic variants associated with 1611 circulating protein abundances identified in 6 large-scale proteomic studies were used as genetic instruments. Effects of the circulating proteins on psychiatric disorders were estimated by Wald ratio or inverse variance-weighted ratio tests. Horizontal pleiotropy, colocalization, and protein-altering effects were examined to validate the assumptions of Mendelian randomization.

RESULTS

Nine circulating protein-to-disease associations withstood multiple sensitivity analyses. Among them, 2 circulating proteins had associations replicated in 3 proteomic studies. A 1 standard deviation increase in the genetically predicted circulating TIMP4 level was associated with a reduced risk of anorexia nervosa (minimum odds ratio [OR] = 0.83; 95% CI, 0.76-0.91) and bipolar disorder (minimum OR = 0.88; 95% CI, 0.82-0.94). A 1 standard deviation increase in the genetically predicted circulating ESAM level was associated with an increased risk of schizophrenia (maximum OR = 1.32; 95% CI, 1.22-1.43). In addition, 58 suggestive protein-to-disease associations warrant validation with observational or experimental evidence. For instance, a 1 standard deviation increase in the ERLEC1-201-to-ERLEC1-202 splice variant ratio was associated with a reduced risk of schizophrenia (OR = 0.94; 95% CI, 0.90-0.97).

CONCLUSIONS

Prioritized circulating proteins appear to influence the risk of psychiatric disease and may be explored as intervention targets.

Genetic analyses of DNA repair pathway associated genes implicate new candidate cancer predisposing genes in ancestrally defined ovarian cancer cases

Not all familial ovarian cancer (OC) cases are explained by pathogenic germline variants in known risk genes. A candidate gene approach involving DNA repair pathway genes was applied to identify rare recurring pathogenic variants in familial OC cases not associated with known OC risk genes from a population exhibiting genetic drift. Whole exome sequencing (WES) data of 15 OC cases from 13 families tested negative for pathogenic variants in known OC risk genes were investigated for candidate variants in 468 DNA repair pathway genes. Filtering and prioritization criteria were applied to WES data to select top candidates for further analyses. Candidates were genotyped in ancestry defined study groups of 214 familial and 998 sporadic OC or breast cancer (BC) cases and 1025 population-matched controls and screened for additional carriers in 605 population-matched OC cases. The candidate genes were also analyzed in WES data from 937 familial or sporadic OC cases of diverse ancestries. Top candidate variants in ERCC5, EXO1, FANCC, NEIL1 and NTHL1 were identified in 5/13 (39%) OC families. Collectively, candidate variants were identified in 7/435 (1.6%) sporadic OC cases and 1/566 (0.2%) sporadic BC cases versus 1/1025 (0.1%) controls. Additional carriers were identified in 6/605 (0.9%) OC cases. Tumour DNA from ERCC5, NEIL1 and NTHL1 variant carriers exhibited loss of the wild-type allele. Carriers of various candidate variants in these genes were identified in 31/937 (3.3%) OC cases of diverse ancestries versus 0-0.004% in cancer-free controls. The strategy of applying a candidate gene approach in a population exhibiting genetic drift identified new candidate OC predisposition variants in DNA repair pathway genes.

Genomic atlas of the plasma metabolome prioritizes metabolites implicated in human diseases

Metabolic processes can influence disease risk and provide therapeutic targets. By conducting genome-wide association studies of 1,091 blood metabolites and 309 metabolite ratios, we identified associations with 690 metabolites at 248 loci and associations with 143 metabolite ratios at 69 loci. Integrating metabolite-gene and gene expression information identified 94 effector genes for 109 metabolites and 48 metabolite ratios. Using Mendelian randomization (MR), we identified 22 metabolites and 20 metabolite ratios having estimated causal effect on 12 traits and diseases, including orotate for estimated bone mineral density, α-hydroxyisovalerate for body mass index and ergothioneine for inflammatory bowel disease and asthma. We further measured the orotate level in a separate cohort and demonstrated that, consistent with MR, orotate levels were positively associated with incident hip fractures. This study provides a valuable resource describing the genetic architecture of metabolites and delivers insights into their roles in common diseases, thereby offering opportunities for therapeutic targets.

GD2 and GD3 gangliosides as diagnostic biomarkers for all stages and subtypes of epithelial ovarian cancer

Background

Ovarian cancer (OC) is the deadliest gynecological cancer, often diagnosed at advanced stages. A fast and accurate diagnostic method for early-stage OC is needed. The tumor marker gangliosides, GD2 and GD3, exhibit properties that make them ideal potential diagnostic biomarkers, but they have never before been quantified in OC. We investigated the diagnostic utility of GD2 and GD3 for diagnosis of all subtypes and stages of OC.

Methods

This retrospective study evaluated GD2 and GD3 expression in biobanked tissue and serum samples from patients with invasive epithelial OC, healthy donors, non-malignant gynecological conditions, and other cancers. GD2 and GD3 levels were evaluated in tissue samples by immunohistochemistry (n=299) and in two cohorts of serum samples by quantitative ELISA. A discovery cohort (n=379) showed feasibility of GD2 and GD3 quantitative ELISA for diagnosing OC, and a subsequent model cohort (n=200) was used to train and cross-validate a diagnostic model.

Results

GD2 and GD3 were expressed in tissues of all OC subtypes and FIGO stages but not in surrounding healthy tissue or other controls. In serum, GD2 and GD3 were elevated in patients with OC. A diagnostic model that included serum levels of GD2+GD3+age was superior to the standard of care (CA125, p<0.001) in diagnosing OC and early-stage (I/II) OC.

Conclusion

GD2 and GD3 expression was associated with high rates of selectivity and specificity for OC. A diagnostic model combining GD2 and GD3 quantification in serum had diagnostic power for all subtypes and all stages of OC, including early stage. Further research exploring the utility of GD2 and GD3 for diagnosis of OC is warranted.

Genetic determinants of polygenic prediction accuracy within a population

Genomic risk prediction is on the emerging path toward personalized medicine. However, the accuracy of polygenic prediction varies strongly in different individuals. Based on up to 352,277 European ancestry participants in the UK Biobank, we constructed polygenic risk scores for 15 physiological and biochemical quantitative traits. We identified a total of 185 polygenic prediction variability quantitative trait loci for 11 traits by Levene's test among 254,376 unrelated individuals. We validated the effects of prediction variability quantitative trait loci using an independent test set of 58,927 individuals. For instance, a score aggregating 51 prediction variability quantitative trait locus variants for triglycerides had the strongest Spearman correlation of 0.185 (P-value <1.0 × 10-300) with the squared prediction errors. We found a strong enrichment of complex genetic effects conferred by prediction variability quantitative trait loci compared to risk loci identified in genome-wide association studies, including 89 prediction variability quantitative trait loci exhibiting dominance effects. Incorporation of dominance effects into polygenic risk scores significantly improved polygenic prediction for triglycerides, low-density lipoprotein cholesterol, vitamin D, and platelet. In conclusion, we have discovered and profiled genetic determinants of polygenic prediction variability for 11 quantitative biomarkers. These findings may assist interpretation of genomic risk prediction in various contexts and encourage novel approaches for constructing polygenic risk scores with complex genetic effects.

Dementia with Lewy bodies post-mortem brains reveal differentially methylated CpG sites with biomarker potential

Dementia with Lewy bodies (DLB) is a common form of dementia with known genetic and environmental interactions. However, the underlying epigenetic mechanisms which reflect these gene-environment interactions are poorly studied. Herein, we measure genome-wide DNA methylation profiles of post-mortem brain tissue (Broadmann area 7) from 15 pathologically confirmed DLB brains and compare them with 16 cognitively normal controls using Illumina MethylationEPIC arrays. We identify 17 significantly differentially methylated CpGs (DMCs) and 17 differentially methylated regions (DMRs) between the groups. The DMCs are mainly located at the CpG islands, promoter and first exon regions. Genes associated with the DMCs are linked to "Parkinson's disease" and "metabolic pathway", as well as the diseases of "severe intellectual disability" and "mood disorders". Overall, our study highlights previously unreported DMCs offering insights into DLB pathogenesis with the possibility that some of these could be used as biomarkers of DLB in the future.

Considering strategies for SNP selection in genetic and polygenic risk scores

Genetic risk scores (GRS) and polygenic risk scores (PRS) are weighted sums of, respectively, several or many genetic variant indicator variables. Although they are being increasingly proposed for clinical use, the best ways to construct them are still actively debated. In this commentary, we present several case studies illustrating practical challenges associated with building or attempting to improve score performance when there is expected to be heterogeneity of disease risk between cohorts or between subgroups of individuals. Specifically, we contrast performance associated with several ways of selecting single nucleotide polymorphisms (SNPs) for inclusion in these scores. By considering GRS and PRS as predictors that are measured with error, insights into their strengths and weaknesses may be obtained, and SNP selection approaches play an important role in defining such errors.

Germline Missense Variants in CDC20 Result in Aberrant Mitotic Progression and Familial Cancer

CDC20 is a coactivator of the anaphase promoting complex/cyclosome (APC/C) and is essential for mitotic progression. APC/CCDC20 is inhibited by the spindle assembly checkpoint (SAC), which prevents premature separation of sister chromatids and aneuploidy in daughter cells. Although overexpression of CDC20 is common in many cancers, oncogenic mutations have never been identified in humans. Using whole-exome sequencing, we identified heterozygous missense CDC20 variants (L151R and N331K) that segregate with ovarian germ cell tumors in two families. Functional characterization showed these mutants retain APC/C activation activity but have impaired binding to BUBR1, a component of the SAC. Expression of L151R and N331K variants promoted mitotic slippage in HeLa cells and primary skin fibroblasts derived from carriers. Generation of mice carrying the N331K variant using CRISPR-Cas9 showed that, although homozygous N331K mice were nonviable, heterozygotes displayed accelerated oncogenicity of Myc-driven cancers. These findings highlight an unappreciated role for CDC20 variants as tumor-promoting genes.

SIGNIFICANCE

Two germline CDC20 missense variants that segregate with cancer in two families compromise the spindle assembly checkpoint and lead to aberrant mitotic progression, which could predispose cells to transformation. See related commentary by Villarroya-Beltri and Malumbres, p. 3432.

The dynamic changes and sex differences of 147 immune-related proteins during acute COVID-19 in 580 individuals

INTRODUCTION

Severe COVID-19 leads to important changes in circulating immune-related proteins. To date it has been difficult to understand their temporal relationship and identify cytokines that are drivers of severe COVID-19 outcomes and underlie differences in outcomes between sexes. Here, we measured 147 immune-related proteins during acute COVID-19 to investigate these questions.

METHODS

We measured circulating protein abundances using the SOMAscan nucleic acid aptamer panel in two large independent hospital-based COVID-19 cohorts in Canada and the United States. We fit generalized additive models with cubic splines from the start of symptom onset to identify protein levels over the first 14 days of infection which were different between severe cases and controls, adjusting for age and sex. Severe cases were defined as individuals with COVID-19 requiring invasive or non-invasive mechanical respiratory support.

RESULTS

580 individuals were included in the analysis. Mean subject age was 64.3 (sd 18.1), and 47% were male. Of the 147 proteins, 69 showed a significant difference between cases and controls (p < 3.4 × 10^-4). Three clusters were formed by 108 highly correlated proteins that replicated in both cohorts, making it difficult to determine which proteins have a true causal effect on severe COVID-19. Six proteins showed sex differences in levels over time, of which 3 were also associated with severe COVID-19: CCL26, IL1RL2, and IL3RA, providing insights to better understand the marked differences in outcomes by sex.

CONCLUSIONS

Severe COVID-19 is associated with large changes in 69 immune-related proteins. Further, five proteins were associated with sex differences in outcomes. These results provide direct insights into immune-related proteins that are strongly influenced by severe COVID-19 infection.

Single base-pair resolution analysis of DNA binding motif with MoMotif reveals an oncogenic function of CTCF zinc-finger 1 mutation

Defining the impact of missense mutations on the recognition of DNA motifs is highly dependent on bioinformatic tools that define DNA binding elements. However, classical motif analysis tools remain limited in their capacity to identify subtle changes in complex binding motifs between distinct conditions. To overcome this limitation, we developed a new tool, MoMotif, that facilitates a sensitive identification, at the single base-pair resolution, of complex, or subtle, alterations to core binding motifs, discerned from ChIP-seq data. We employed MoMotif to define the previously uncharacterized recognition motif of CTCF zinc-finger 1 (ZF1), and to further define the impact of CTCF ZF1 mutation on its association with chromatin. Mutations of CTCF ZF1 are exclusive to breast cancer and are associated with metastasis and therapeutic resistance, but the underlying mechanisms are unclear. Using MoMotif, we identified an extension of the CTCF core binding motif, necessitating a functional ZF1 to bind appropriately. Using a combination of ChIP-Seq and RNA-Seq, we discover that the inability to bind this extended motif drives an altered transcriptional program associated with the oncogenic phenotypes observed clinically. Our study demonstrates that MoMotif is a powerful new tool for comparative ChIP-seq analysis and characterising DNA-protein contacts.

Identifying Causes of Fracture Beyond Bone Mineral Density: Evidence From Human Genetics

New therapies may help to prevent osteoporotic fractures other than through increasing bone mineral density (BMD). Because fracture risk has an important genetic component, we aim to identify loci increasing fracture risk that do not decrease BMD, using a recently-proposed structural equation model adapted to remove genetic influences of BMD on fracture risk. We used summary statistics of the largest genome-wide association studies (GWASs) for BMD and for fracture in these analyses. We next estimated the genetic correlation between the non-BMD or BMD-related genetic effects and other clinical risk factors for fracture. Last, based on white British participants in the UK Biobank, we conducted genetic risk score analyses to assess whether the aggregated genetic effects conferred increased major osteoporotic fracture risk. We found that only three loci affecting fracture risk exhibited genetic effects not mediated by BMD: SOST, CPED1-WNT16, and RSPO3, while these three loci simultaneously conferred BMD-related effects. No strong genetic associations between non-BMD or BMD-related effects and 16 clinical risk factors were observed. However, non-BMD effects might be genetic correlated with hip bone size. In the UK Biobank, a 1 standard deviation (1-SD) increase in the non-BMD genetic risk score conferred an odds ratio of 1.17 for incident major osteoporotic fracture, compared to 1.29 by a BMD-related genetic risk score. Our study suggests that the majority of common genetic predisposition toward fracture risk acts upon BMD. Although non-BMD genetic effects may exist, they are not strongly correlated with most traditional clinical risk factors. Risk loci harboring non-BMD genetic effects may influence other perspectives of bone quality, or confer effects that existing GWASs fail to capture, but they demonstrate weaker impact on fracture risk than BMD-related genetic effects. These findings suggest that most successful drug development programs for osteoporosis should focus on pathways identified through BMD-associated loci. © 2022 American Society for Bone and Mineral Research (ASBMR).

Polygenic risk score as a possible tool for identifying familial monogenic causes of complex diseases

PURPOSE

The study aimed to evaluate whether polygenic risk scores could be helpful in addition to family history for triaging individuals to undergo deep-depth diagnostic sequencing for identifying monogenic causes of complex diseases.

METHODS

Among 44,550 exome-sequenced European ancestry UK Biobank participants, we identified individuals with a clinically reported or computationally predicted monogenic pathogenic variant for breast cancer, bowel cancer, heart disease, diabetes, or Alzheimer disease. We derived polygenic risk scores for these diseases. We tested whether a polygenic risk score could identify rare pathogenic variant heterozygotes among individuals with a parental disease history.

RESULTS

Monogenic causes of complex diseases were more prevalent among individuals with a parental disease history than in the rest of the population. Polygenic risk scores showed moderate discriminative power to identify familial monogenic causes. For instance, we showed that prescreening the patients with a polygenic risk score for type 2 diabetes can prioritize individuals to undergo diagnostic sequencing for monogenic diabetes variants and reduce needs for such sequencing by up to 37%.

CONCLUSION

Among individuals with a family history of complex diseases, those with a low polygenic risk score are more likely to have monogenic causes of the disease and could be prioritized to undergo genetic testing.

A Bayesian hierarchical model for improving measurement of 5mC and 5hmC levels: Toward revealing associations between phenotypes and methylation states

5-hydroxymethylcytosine (5hmC) is a methylation state linked with gene regulation, commonly found in cells of the central nervous system. 5hmC is associated with demethylation of cytosines from 5-methylcytosine (5mC) to the unmethylated state. The presence of 5hmC can be inferred by a paired experiment involving bisulfite and oxidation-bisulfite treatments on the same sample, followed by a methylation assay using a platform such as the Illumina Infinium MethylationEPIC BeadChip (EPIC). Existing methods for analysis of the resulting EPIC data are not ideal. Most approaches ignore the correlation between the two experiments and any imprecision associated with DNA damage from the additional treatment. Estimates of 5mC/5hmC levels free from these limitations are desirable to reveal associations between methylation states and phenotypes. We propose a hierarchical Bayesian method called Constrained HYdroxy Methylation Estimation (CHYME) to simultaneously estimate 5mC/5hmC signals as well as any associations between these signals and covariates or phenotypes, while accounting for the potential impact of DNA damage and dependencies induced by the experimental design. Simulations show that CHYME has valid type 1 error and better power than a range of alternative methods, including the popular OxyBS method and linear models on transformed proportions. Other methods we examined suffer from hugely inflated type 1 error for inference on 5hmC proportions. We use CHYME to explore genome-wide associations between 5mC/5hmC levels and cause of death in postmortem prefrontal cortex brain tissue samples. These analyses indicate that CHYME is a useful tool to reveal phenotypic associations with 5mC/5hmC levels.

Capturing additional genetic risk from family history for improved polygenic risk prediction

Family history of complex traits may reflect transmitted rare pathogenic variants, intra-familial shared exposures to environmental and lifestyle factors, as well as a common genetic predisposition. We developed a latent factor model to quantify trait heritability in excess of that captured by a common variant-based polygenic risk score, but inferable from family history. For 941 children in the Avon Longitudinal Study of Parents and Children cohort, a joint predictor combining a polygenic risk score for height and mid-parental height was able to explain ^~55% of the total variance in sex-adjusted adult height z-scores, close to the estimated heritability. Marginal yet consistent risk prediction improvements were also achieved among ^~400,000 European ancestry participants for 11 complex diseases in the UK Biobank. Our work showcases a paradigm for risk calculation, and supports incorporation of family history into polygenic risk score-based genetic risk prediction models.

A functionally impaired missense variant identified in French Canadian families implicates FANCI as a candidate ovarian cancer-predisposing gene

Background

Familial ovarian cancer (OC) cases not harbouring pathogenic variants in either of the BRCA1 and BRCA2 OC-predisposing genes, which function in homologous recombination (HR) of DNA, could involve pathogenic variants in other DNA repair pathway genes.

Methods

Whole exome sequencing was used to identify rare variants in HR genes in a BRCA1 and BRCA2 pathogenic variant negative OC family of French Canadian (FC) ancestry, a population exhibiting genetic drift. OC cases and cancer-free individuals from FC and non-FC populations were investigated for carrier frequency of FANCI c.1813C>T; p.L605F, the top-ranking candidate. Gene and protein expression were investigated in cancer cell lines and tissue microarrays, respectively.

Results

In FC subjects, c.1813C>T was more common in familial (7.1%, 3/42) than sporadic (1.6%, 7/439) OC cases (P = 0.048). Carriers were detected in 2.5% (74/2950) of cancer-free females though female/male carriers were more likely to have a first-degree relative with OC (121/5249, 2.3%; Spearman correlation = 0.037; P = 0.011), suggesting a role in risk. Many of the cancer-free females had host factors known to reduce risk to OC which could influence cancer risk in this population. There was an increased carrier frequency of FANCI c.1813C>T in BRCA1 and BRCA2 pathogenic variant negative OC families, when including the discovery family, compared to cancer-free females (3/23, 13%; OR = 5.8; 95%CI = 1.7–19; P = 0.005). In non-FC subjects, 10 candidate FANCI variants were identified in 4.1% (21/516) of Australian OC cases negative for pathogenic variants in BRCA1 and BRCA2, including 10 carriers of FANCI c.1813C>T. Candidate variants were significantly more common in familial OC than in sporadic OC (P = 0.04). Localization of FANCD2, part of the FANCI-FANCD2 (ID2) binding complex in the Fanconi anaemia (FA) pathway, to sites of induced DNA damage was severely impeded in cells expressing the p.L605F isoform. This isoform was expressed at a reduced level, destabilized by DNA damaging agent treatment in both HeLa and OC cell lines, and exhibited sensitivity to cisplatin but not to a poly (ADP-ribose) polymerase inhibitor. By tissue microarray analyses, FANCI protein was consistently expressed in fallopian tube epithelial cells and only expressed at low-to-moderate levels in 88% (83/94) of OC samples.

Conclusions

This is the first study to describe candidate OC variants in FANCI, a member of the ID2 complex of the FA DNA repair pathway. Our data suggest that pathogenic FANCI variants may modify OC risk in cancer families.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13073-021-00998-5.

Block coordinate descent algorithm improves variable selection and estimation in error-in-variables regression

Medical research increasingly includes high‐dimensional regression modeling with a need for error‐in‐variables methods. The Convex Conditioned Lasso (CoCoLasso) utilizes a reformulated Lasso objective function and an error‐corrected cross‐validation to enable error‐in‐variables regression, but requires heavy computations. Here, we develop a Block coordinate Descent Convex Conditioned Lasso (BDCoCoLasso) algorithm for modeling high‐dimensional data that are only partially corrupted by measurement error. This algorithm separately optimizes the estimation of the uncorrupted and corrupted features in an iterative manner to reduce computational cost, with a specially calibrated formulation of cross‐validation error. Through simulations, we show that the BDCoCoLasso algorithm successfully copes with much larger feature sets than CoCoLasso, and as expected, outperforms the naïve Lasso with enhanced estimation accuracy and consistency, as the intensity and complexity of measurement errors increase. Also, a new smoothly clipped absolute deviation penalization option is added that may be appropriate for some data sets. We apply the BDCoCoLasso algorithm to data selected from the UK Biobank. We develop and showcase the utility of covariate‐adjusted genetic risk scores for body mass index, bone mineral density, and lifespan. We demonstrate that by leveraging more information than the naïve Lasso in partially corrupted data, the BDCoCoLasso may achieve higher prediction accuracy. These innovations, together with an R package, BDCoCoLasso, make error‐in‐variables adjustments more accessible for high‐dimensional data sets. We posit the BDCoCoLasso algorithm has the potential to be widely applied in various fields, including genomics‐facilitated personalized medicine research.

Early or Late Gestational Exposure to Maternal Immune Activation Alters Neurodevelopmental Trajectories in Mice: An Integrated Neuroimaging, Behavioral, and Transcriptional Study

BACKGROUND

Exposure to maternal immune activation (MIA) in utero is a risk factor for neurodevelopmental disorders later in life. The impact of the gestational timing of MIA exposure on downstream development remains unclear.

METHODS

We characterized neurodevelopmental trajectories of mice exposed to the viral mimetic poly I:C (polyinosinic:polycytidylic acid) either on gestational day 9 (early) or on day 17 (late) using longitudinal structural magnetic resonance imaging from weaning to adulthood. Using multivariate methods, we related neuroimaging and behavioral variables for the time of greatest alteration (adolescence/early adulthood) and identified regions for further investigation using RNA sequencing.

RESULTS

Early MIA exposure was associated with accelerated brain volume increases in adolescence/early adulthood that normalized in later adulthood in the striatum, hippocampus, and cingulate cortex. Similarly, alterations in anxiety-like, stereotypic, and sensorimotor gating behaviors observed in adolescence normalized in adulthood. MIA exposure in late gestation had less impact on anatomical and behavioral profiles. Multivariate maps associated anxiety-like, social, and sensorimotor gating deficits with volume of the dorsal and ventral hippocampus and anterior cingulate cortex, among others. The most transcriptional changes were observed in the dorsal hippocampus, with genes enriched for fibroblast growth factor regulation, autistic behaviors, inflammatory pathways, and microRNA regulation.

CONCLUSIONS

Leveraging an integrated hypothesis- and data-driven approach linking brain-behavior alterations to the transcriptome, we found that MIA timing differentially affects offspring development. Exposure in late gestation leads to subthreshold deficits, whereas exposure in early gestation perturbs brain development mechanisms implicated in neurodevelopmental disorders.

Thousands of CpGs Show DNA Methylation Differences in ACPA-Positive Individuals

High levels of anti-citrullinated protein antibodies (ACPA) are often observed prior to a diagnosis of rheumatoid arthritis (RA). We undertook a replication study to confirm CpG sites showing evidence of differential methylation in subjects positive vs. negative for ACPA, in a new subset of 112 individuals sampled from the population cohort and biobank CARTaGENE in Quebec, Canada. Targeted custom capture bisulfite sequencing was conducted at approximately 5.3 million CpGs located in regulatory or hypomethylated regions from whole blood; library and protocol improvements had been instituted between the original and this replication study, enabling better coverage and additional identification of differentially methylated regions (DMRs). Using binomial regression models, we identified 19,472 ACPA-associated differentially methylated cytosines (DMCs), of which 430 overlapped with the 1909 DMCs reported by the original study; 814 DMRs of relevance were clustered by grouping adjacent DMCs into regions. Furthermore, we performed an additional integrative analysis by looking at the DMRs that overlap with RA related loci published in the GWAS Catalog, and protein-coding genes associated with these DMRs were enriched in the biological process of cell adhesion and involved in immune-related pathways.

Abstract 2056: The genomic landscape of carriers of rare variants in FANCI, a new candidate ovarian cancer predisposing gene

The potentially pathogenic variant (PPV), FANCI c.1813C&gt;T; p.L605F, in a new candidate ovarian cancer (OC) predisposing gene was discovered by whole exome sequencing (WES) of familial OC cases from the founder French Canadian (FC) population for discovering new OC predisposing genes. Modeling this variant in cellulo suggested this variant encodes an unstable protein. FANCI is an essential member upstream of the homologous recombination DNA repair pathway and intersects BRCA1 and BRCA2 function, proteins encoded by genes involved in hereditary OC. Investigating the FC founder population facilitates the discovery of PPVs as they are more likely to harbor recurrent variants due to common ancestors, increasing the likelihood of identifying candidate genes in cancer families.

To further support the candidacy of FANCI as a new OC predisposing gene, we investigated the germline landscape of c.1813C&gt;T carrier FC OC cases for co-occurring PPVs in known or proposed new (emerging) OC predisposing genes or other genes involved in similar pathways (DNA repair).

Using WES of peripheral blood lymphocyte DNA, the genomic landscape of 10 FANCI carriers were investigated for heterozygous PPVs in 276 DNA repair pathways genes, which included known (BRCA1, BRCA2, MSH2, MLH1, MSH6, PMS2) and emerging (BRIP1, RAD51C, RAD51D) OC predisposing genes. Top ranking candidate variants (minor allele frequency &lt;1%) were identified using 11 different in silico tools that assessed amino acid conversation or potential pathogenicity. Pathogenicity of known and emerging OC predisposing genes was assessed using BRCAExchange (www.brcaexchange.org) and ClinVar (www.ncbi.nlm.nih.gov/clinvar/). A similar analysis was done with WES data from 13 FC OC cases harboring pathogenic BRCA1/BRCA2 variants.

We identified 31 variants in 27 genes in FANCI c.1813C&gt;T carriers. A previously known carrier of a pathogenic BRCA1 variant (c.2836_2837del; p.Ile946GlnfsTer5) was also identified in a familial case. No carriers of other pathogenic variants in known or emerging OC predisposing genes were found. However, FANCI carriers were found to carry at least one other PPV in a DNA repair pathway gene (range 2-9; average=4.1). There were no other carriers of variants in common among all OC cases, though at least 2 cases carried a variant in the same gene. FC OC cases harboring BRCA1/BRCA2 variants carried at least one other PPV in a DNA repair pathway gene (range 2-7; average 4).

It is possible that the identified variants influence or modify risk in conjunction with FANCI, though no PPV was identified in all carriers. As new cancer predisposing genes are identified it will become increasingly important to characterize the genetic context in which variants are identified. This will allow for further insight to clinical translatability once penetrance has been established.

Citation Format: Caitlin Fierheller, Wejdan M Alenezi, Corinne Serruya, Timothée Revil, Javad Nadaf, Anne-Marie Mes-Masson, Diane Provencher, William D Foulkes, Zaki El Haffaf, Celia M T Greenwood, Jean-Yves Masson, Jiannis Ragoussis, Patricia N Tonin. The genomic landscape of carriers of rare variants in FANCI, a new candidate ovarian cancer predisposing gene [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 2056.

Detecting cord blood cell type-specific epigenetic associations with gestational diabetes mellitus and early childhood growth

BACKGROUND

Epigenome-wide association studies (EWAS) have provided opportunities to understand the role of epigenetic mechanisms in development and pathophysiology of many chronic diseases. However, an important limitation of conventional EWAS is that profiles of epigenetic variability are often obtained in samples of mixed cell types. Here, we aim to assess whether changes in cord blood DNA methylation (DNAm) associated with gestational diabetes mellitus (GDM) exposure and early childhood growth markers occur in a cell type-specific manner.

RESULTS

We analyzed 275 cord blood samples collected at delivery from a prospective pre-birth cohort with genome-wide DNAm profiled by the Illumina MethylationEPIC array. We estimated proportions of seven common cell types in each sample using a cord blood-specific DNAm reference panel. Leveraging a recently developed approach named CellDMC, we performed cell type-specific EWAS to identify CpG loci significantly associated with GDM, or 3-year-old body mass index (BMI) z-score. A total of 1410 CpG loci displayed significant cell type-specific differences in methylation level between 23 GDM cases and 252 controls with a false discovery rate < 0.05. Gene Ontology enrichment analysis indicated that LDL transportation emerged from CpG specifically identified from B-cells DNAm analyses and the mitogen-activated protein kinase pathway emerged from CpG specifically identified from natural killer cells DNAm analyses. In addition, we identified four and six loci associated with 3-year-old BMI z-score that were specific to CD8+ T-cells and monocytes, respectively. By performing genome-wide permutation tests, we validated that most of our detected signals had low false positive rates.

CONCLUSION

Compared to conventional EWAS adjusting for the effects of cell type heterogeneity, the proposed approach based on cell type-specific EWAS could provide additional biologically meaningful associations between CpG methylation, prenatal maternal GDM or 3-year-old BMI. With careful validation, these findings may provide new insights into the pathogenesis, programming, and consequences of related childhood metabolic dysregulation. Therefore, we propose that cell type-specific analyses are worth cautious explorations.

A Polygenic Risk Score to Predict Future Adult Short Stature Among Children

CONTEXT

Adult height is highly heritable, yet no genetic predictor has demonstrated clinical utility compared to mid-parental height.

OBJECTIVE

To develop a polygenic risk score for adult height and evaluate its clinical utility.

DESIGN

A polygenic risk score was constructed based on meta-analysis of genomewide association studies and evaluated on the Avon Longitudinal Study of Parents and Children (ALSPAC) cohort.

SUBJECTS

Participants included 442 599 genotyped White British individuals in the UK Biobank and 941 genotyped child-parent trios of European ancestry in the ALSPAC cohort.

INTERVENTIONS

None.

MAIN OUTCOME MEASURES

Standing height was measured using stadiometer; Standing height 2 SDs below the sex-specific population average was considered as short stature.

RESULTS

Combined with sex, a polygenic risk score captured 71.1% of the total variance in adult height in the UK Biobank. In the ALSPAC cohort, the polygenic risk score was able to identify children who developed adulthood short stature with an area under the receiver operating characteristic curve (AUROC) of 0.84, which is close to that of mid-parental height. Combining this polygenic risk score with mid-parental height or only one of the child's parent's height could improve the AUROC to at most 0.90. The polygenic risk score could also substitute mid-parental height in age-specific Khamis-Roche height predictors and achieve an equally strong discriminative power in identifying children with a short stature in adulthood.

CONCLUSIONS

A polygenic risk score could be considered as an alternative or adjunct to mid-parental height to improve screening for children at risk of developing short stature in adulthood in European ancestry populations.

Genome-wide analysis of gene dosage in 24,092 individuals estimates that 10,000 genes modulate cognitive ability

Genomic copy number variants (CNVs) are routinely identified and reported back to patients with neuropsychiatric disorders, but their quantitative effects on essential traits such as cognitive ability are poorly documented. We have recently shown that the effect size of deletions on cognitive ability can be statistically predicted using measures of intolerance to haploinsufficiency. However, the effect sizes of duplications remain unknown. It is also unknown if the effect of multigenic CNVs are driven by a few genes intolerant to haploinsufficiency or distributed across tolerant genes as well. Here, we identified all CNVs > 50 kilobases in 24,092 individuals from unselected and autism cohorts with assessments of general intelligence. Statistical models used measures of intolerance to haploinsufficiency of genes included in CNVs to predict their effect size on intelligence. Intolerant genes decrease general intelligence by 0.8 and 2.6 points of intelligence quotient when duplicated or deleted, respectively. Effect sizes showed no heterogeneity across cohorts. Validation analyses demonstrated that models could predict CNV effect sizes with 78% accuracy. Data on the inheritance of 27,766 CNVs showed that deletions and duplications with the same effect size on intelligence occur de novo at the same frequency. We estimated that around 10,000 intolerant and tolerant genes negatively affect intelligence when deleted, and less than 2% have large effect sizes. Genes encompassed in CNVs were not enriched in any GOterms but gene regulation and brain expression were GOterms overrepresented in the intolerant subgroup. Such pervasive effects on cognition may be related to emergent properties of the genome not restricted to a limited number of biological pathways.

A 10-color flow cytometry panel for diagnosis and minimal residual disease in chronic lymphocytic leukemia

Diagnosis and minimal residual disease (MRD) monitoring of chronic lymphocytic leukemia (CLL) by flow cytometry currently requires multiple antibody panels. We added CD23 and CD200 to the EuroFlow^TM lymphoid screening tube (LST) to create a 10-color modified LST (mLST) capable of diagnosing typical CLL in a single tube. We then explored if the mLST could be used for MRD by comparing its performance to the European Research Initiative on CLL (ERIC) panel using spiked cryopreserved and fresh patient samples. Over 1 year of use in our clinical laboratory, the mLST diagnosed CLL without further immunophenotyping in 56% of samples with an abnormal clone. There was good agreement in MRD results between the mLST and ERIC panels. Therefore, the mLST can streamline CLL diagnosis by reducing technician time and the number of panels required. It may have the potential to screen for MRD in laboratories without access to dedicated panels (ERIC).

High Systolic Blood Pressure at Hospital Admission Is an Important Risk Factor in Models Predicting Outcome of COVID-19 Patients

BACKGROUND

The risk that coronavirus disease 2019 (COVID-19) patients develop critical illness that can be fatal depends on their age and immune status and may also be affected by comorbidities like hypertension. The goal of this study was to develop models that predict outcome using parameters collected at admission to the hospital.

METHODS AND RESULTS

This is a retrospective single-center cohort study of COVID-19 patients at the Seventh Hospital of Wuhan City, China. Forty-three demographic, clinical, and laboratory parameters collected at admission plus discharge/death status, days from COVID-19 symptoms onset, and days of hospitalization were analyzed. From 157 patients, 120 were discharged and 37 died. Pearson correlations showed that hypertension and systolic blood pressure (SBP) were associated with death and respiratory distress parameters. A penalized logistic regression model efficiently predicts the probability of death with 13 of 43 variables. A regularized Cox regression model predicts the probability of survival with 7 of above 13 variables. SBP but not hypertension was a covariate in both mortality and survival prediction models. SBP was elevated in deceased compared with discharged COVID-19 patients.

CONCLUSIONS

Using an unbiased approach, we developed models predicting outcome of COVID-19 patients based on data available at hospital admission. This can contribute to evidence-based risk prediction and appropriate decision-making at hospital triage to provide the most appropriate care and ensure the best patient outcome. High SBP, a cause of end-organ damage and an important comorbid factor, was identified as a covariate in both mortality and survival prediction models.

Candidate Markers of Olaparib Response from Genomic Data Analyses of Human Cancer Cell Lines

The benefit of PARP inhibitor olaparib in relapsed and advanced high-grade serous ovarian carcinoma (HGSOC) is well established especially in BRCA1/2 mutation carriers. Identification of additional biomarkers can help expand the population of patients most likely to benefit from olaparib treatment. To identify candidate markers of olaparib response we analyzed genomic and in vitro olaparib response data from two independent groups of cancer cell lines. Using pan-cancer cell lines (n = 896) from the Genomics of Drug Sensitivity in Cancer database, we applied linear regression methods to identify statistically significant gene predictors of olaparib response based on mRNA expression. We then analyzed whole exome sequencing and mRNA gene expression data from our collection of 18 HGSOC cell lines previously classified as sensitive, intermediate, or resistant based on in vitro olaparib response for mutations, copy number variation and differential expression of candidate olaparib response genes. We identify genes previously associated with olaparib response (SLFN11, ABCB1), and discover novel candidate olaparib sensitivity genes with known functions including interaction with PARP1 (PUM3, EEF1A1) and involvement in homologous recombination DNA repair (ELP4). Further investigations at experimental and clinical levels are required to validate novel candidates, and ultimately determine their efficacy as potential biomarkers of olaparib sensitivity.

Improved prediction of fracture risk leveraging a genome-wide polygenic risk score

BACKGROUND

Accurately quantifying the risk of osteoporotic fracture is important for directing appropriate clinical interventions. While skeletal measures such as heel quantitative speed of sound (SOS) and dual-energy X-ray absorptiometry bone mineral density are able to predict the risk of osteoporotic fracture, the utility of such measurements is subject to the availability of equipment and human resources. Using data from 341,449 individuals of white British ancestry, we previously developed a genome-wide polygenic risk score (PRS), called gSOS, that captured 25.0% of the total variance in SOS. Here, we test whether gSOS can improve fracture risk prediction.

METHODS

We examined the predictive power of gSOS in five genome-wide genotyped cohorts, including 90,172 individuals of European ancestry and 25,034 individuals of Asian ancestry. We calculated gSOS for each individual and tested for the association between gSOS and incident major osteoporotic fracture and hip fracture. We tested whether adding gSOS to the risk prediction models had added value over models using other commonly used clinical risk factors.

RESULTS

A standard deviation decrease in gSOS was associated with an increased odds of incident major osteoporotic fracture in populations of European ancestry, with odds ratios ranging from 1.35 to 1.46 in four cohorts. It was also associated with a 1.26-fold (95% confidence interval (CI) 1.13-1.41) increased odds of incident major osteoporotic fracture in the Asian population. We demonstrated that gSOS was more predictive of incident major osteoporotic fracture (area under the receiver operating characteristic curve (AUROC) = 0.734; 95% CI 0.727-0.740) and incident hip fracture (AUROC = 0.798; 95% CI 0.791-0.805) than most traditional clinical risk factors, including prior fracture, use of corticosteroids, rheumatoid arthritis, and smoking. We also showed that adding gSOS to the Fracture Risk Assessment Tool (FRAX) could refine the risk prediction with a positive net reclassification index ranging from 0.024 to 0.072.

CONCLUSIONS

We generated and validated a PRS for SOS which was associated with the risk of fracture. This score was more strongly associated with the risk of fracture than many clinical risk factors and provided an improvement in risk prediction. gSOS should be explored as a tool to improve risk stratification to identify individuals at high risk of fracture.

Understanding the impact of preprocessing pipelines on neuroimaging cortical surface analyses

BACKGROUND

The choice of preprocessing pipeline introduces variability in neuroimaging analyses that affects the reproducibility of scientific findings. Features derived from structural and functional MRI data are sensitive to the algorithmic or parametric differences of preprocessing tasks, such as image normalization, registration, and segmentation to name a few. Therefore it is critical to understand and potentially mitigate the cumulative biases of pipelines in order to distinguish biological effects from methodological variance.

METHODS

Here we use an open structural MRI dataset (ABIDE), supplemented with the Human Connectome Project, to highlight the impact of pipeline selection on cortical thickness measures. Specifically, we investigate the effect of (i) software tool (e.g., ANTS, CIVET, FreeSurfer), (ii) cortical parcellation (Desikan-Killiany-Tourville, Destrieux, Glasser), and (iii) quality control procedure (manual, automatic). We divide our statistical analyses by (i) method type, i.e., task-free (unsupervised) versus task-driven (supervised); and (ii) inference objective, i.e., neurobiological group differences versus individual prediction.

RESULTS

Results show that software, parcellation, and quality control significantly affect task-driven neurobiological inference. Additionally, software selection strongly affects neurobiological (i.e. group) and individual task-free analyses, and quality control alters the performance for the individual-centric prediction tasks.

CONCLUSIONS

This comparative performance evaluation partially explains the source of inconsistencies in neuroimaging findings. Furthermore, it underscores the need for more rigorous scientific workflows and accessible informatics resources to replicate and compare preprocessing pipelines to address the compounding problem of reproducibility in the age of large-scale, data-driven computational neuroscience.

Neuropsychiatric symptoms are early indicators of an upcoming metabolic decline in Alzheimer's disease

Background

Neuropsychiatric symptoms (NPS) are increasingly recognized as early non-cognitive manifestations in the Alzheimer’s disease (AD) continuum. However, the role of NPS as an early marker of pathophysiological progression in AD remains unclear. Dominantly inherited AD (DIAD) mutation carriers are young individuals who are destined to develop AD in future due to the full penetrance of the genetic mutation. Hence, the study of DIAD mutation carriers enables the evaluation of the associations between pure AD pathophysiology and metabolic correlates of NPS without the confounding effects of co-existing pathologies. In this longitudinal study, we aimed to identify regional brain metabolic dysfunctions associated with NPS in cognitively intact DIAD mutation carriers.

Methods

We stratified 221 cognitively intact participants from the Dominantly Inherited Alzheimer’s Network according to their mutation carrier status. The interactions of NPS measured by the Neuropsychiatric Inventory-Questionnaire (NPI-Q), age, and estimated years to symptom onset (EYO) as a function of metabolism measured by [¹⁸F]flurodeoxyglucose ([¹⁸F]FDG) positron emission tomography, were evaluated by the mixed-effects regression model with family-level random effects in DIAD mutation carriers and non-carriers. Exploratory factor analysis was performed to identify the neuropsychiatric subsyndromes in DIAD mutation carriers using the NPI-Q sub-components. Then the effects of interactions between specific neuropsychiatric subsyndromes and EYO on metabolism were evaluated with the mixed-effects regression model.

Results

A total of 119 mutation carriers and 102 non-carriers were studied. The interaction of higher NPI-Q and shorter EYO was associated with more rapid declines of global and regional [¹⁸F]FDG uptake in the posterior cingulate and ventromedial prefrontal cortices, the bilateral parietal lobes and the right insula in DIAD mutation carriers. The neuropsychiatric subsyndromes of agitation, disinhibition, irritability and depression interacted with the EYO to drive the [¹⁸F]FDG uptake decline in the DIAD mutation carriers. The interaction of NPI and EYO was not associated with [¹⁸F]FDG uptake in DIAD mutation non-carriers.

Conclusions

The NPS in cognitively intact DIAD mutation carriers may be a clinical indicator of subsequent metabolic decline in brain networks vulnerable to AD, which supports the emerging conceptual framework that NPS represent early manifestations of neuronal injury in AD. Further studies using different methodological approaches to identify NPS in preclinical AD are needed to validate our findings.

Maternal Prenatal Mood, Pregnancy-Specific Worries, and Early Child Psychopathology: Findings From the DREAM BIG Consortium

OBJECTIVE

Few studies have attempted to identify how distinct dimensions of maternal prenatal affective symptoms relate to offspring psychopathology. We defined latent dimensions of women's prenatal affective symptoms and pregnancy-specific worries to examine their association with early offspring psychopathology in three prenatal cohorts.

METHOD

Data were used from three cohorts of the DREAM-BIG consortium: Avon Longitudinal Study of Parents and Children (ALSPAC [N = 12,515]), Generation R (N = 6,803), and the Canadian prenatal cohort Maternal Adversity, Vulnerability, and Neurodevelopment (MAVAN [N = 578]). Maternal prenatal affective symptoms and pregnancy-specific worries were assessed using different measures in each cohort. Through confirmatory factor analyses, we determined whether comparable latent dimensions of prenatal maternal affective symptoms existed across the cohorts. We used structural equation models to examine cohort-specific associations between these dimensions and offspring psychopathology at 4 to 8 years of age (general psychopathology, specific internalizing and externalizing previously derived using confirmatory factor analyses). Cohort-based estimates were meta-analyzed using inverse variance-weighing.

RESULTS

Four prenatal maternal factors were similar in all cohorts: a general affective symptoms factor and three specific factors-an anxiety/depression factor, a somatic factor, and a pregnancy-specific worries factor. In meta-analyses, both the general affective symptoms factor and pregnancy-specific worries factor were independently associated with offspring general psychopathology. The general affective symptoms factor was further associated with offspring specific internalizing problems. There were no associations with specific externalizing problems.

CONCLUSION

These replicated findings of independent and adverse effects for prenatal general affective symptoms and pregnancy-specific worries on child mental health support the need for specific interventions in pregnancy.

Statistical power in COVID-19 case-control host genomic study design

The identification of genetic variation that directly impacts infection susceptibility to SARS-CoV-2 and disease severity of COVID-19 is an important step towards risk stratification, personalized treatment plans, therapeutic, and vaccine development and deployment. Given the importance of study design in infectious disease genetic epidemiology, we use simulation and draw on current estimates of exposure, infectivity, and test accuracy of COVID-19 to demonstrate the feasibility of detecting host genetic factors associated with susceptibility and severity in published COVID-19 study designs. We demonstrate that limited phenotypic data and exposure/infection information in the early stages of the pandemic significantly impact the ability to detect most genetic variants with moderate effect sizes, especially when studying susceptibility to SARS-CoV-2 infection. Our insights can aid in the interpretation of genetic findings emerging in the literature and guide the design of future host genetic studies.

Individuals with common diseases but with a low polygenic risk score could be prioritized for rare variant screening

PURPOSE

Identifying rare genetic causes of common diseases can improve diagnostic and treatment strategies, but incurs high costs. We tested whether individuals with common disease and low polygenic risk score (PRS) for that disease generated from less expensive genome-wide genotyping data are more likely to carry rare pathogenic variants.

METHODS

We identified patients with one of five common complex diseases among 44,550 individuals who underwent exome sequencing in the UK Biobank. We derived PRS for these five diseases, and identified pathogenic rare variant heterozygotes. We tested whether individuals with disease and low PRS were more likely to carry rare pathogenic variants.

RESULTS

While rare pathogenic variants conferred, at most, 5.18-fold (95% confidence interval [CI]: 2.32-10.13) increased odds of disease, a standard deviation increase in PRS, at most, increased the odds of disease by 5.25-fold (95% CI: 5.06-5.45). Among diseased patients, a standard deviation decrease in the PRS was associated with, at most, 2.82-fold (95% CI: 1.14-7.46) increased odds of identifying rare variant heterozygotes.

CONCLUSION

Rare pathogenic variants were more prevalent among affected patients with a low PRS. Therefore, prioritizing individuals for sequencing who have disease but low PRS may increase the yield of sequencing studies to identify rare variant heterozygotes.

Development of a polygenic risk score to improve screening for fracture risk: A genetic risk prediction study

BACKGROUND

Since screening programs identify only a small proportion of the population as eligible for an intervention, genomic prediction of heritable risk factors could decrease the number needing to be screened by removing individuals at low genetic risk. We therefore tested whether a polygenic risk score for heel quantitative ultrasound speed of sound (SOS)-a heritable risk factor for osteoporotic fracture-can identify low-risk individuals who can safely be excluded from a fracture risk screening program.

METHODS AND FINDINGS

A polygenic risk score for SOS was trained and selected in 2 separate subsets of UK Biobank (comprising 341,449 and 5,335 individuals). The top-performing prediction model was termed "gSOS", and its utility in fracture risk screening was tested in 5 validation cohorts using the National Osteoporosis Guideline Group clinical guidelines (N = 10,522 eligible participants). All individuals were genome-wide genotyped and had measured fracture risk factors. Across the 5 cohorts, the average age ranged from 57 to 75 years, and 54% of studied individuals were women. The main outcomes were the sensitivity and specificity to correctly identify individuals requiring treatment with and without genetic prescreening. The reference standard was a bone mineral density (BMD)-based Fracture Risk Assessment Tool (FRAX) score. The secondary outcomes were the proportions of the screened population requiring clinical-risk-factor-based FRAX (CRF-FRAX) screening and BMD-based FRAX (BMD-FRAX) screening. gSOS was strongly correlated with measured SOS (r2 = 23.2%, 95% CI 22.7% to 23.7%). Without genetic prescreening, guideline recommendations achieved a sensitivity and specificity for correct treatment assignment of 99.6% and 97.1%, respectively, in the validation cohorts. However, 81% of the population required CRF-FRAX tests, and 37% required BMD-FRAX tests to achieve this accuracy. Using gSOS in prescreening and limiting further assessment to those with a low gSOS resulted in small changes to the sensitivity and specificity (93.4% and 98.5%, respectively), but the proportions of individuals requiring CRF-FRAX tests and BMD-FRAX tests were reduced by 37% and 41%, respectively. Study limitations include a reliance on cohorts of predominantly European ethnicity and use of a proxy of fracture risk.

CONCLUSIONS

Our results suggest that the use of a polygenic risk score in fracture risk screening could decrease the number of individuals requiring screening tests, including BMD measurement, while maintaining a high sensitivity and specificity to identify individuals who should be recommended an intervention.

Bayesian Hyper-LASSO Classification for Feature Selection with Application to Endometrial Cancer RNA-seq Data

Feature selection is demanded in many modern scientific research problems that use high-dimensional data. A typical example is to identify gene signatures that are related to a certain disease from high-dimensional gene expression data. The expression of genes may have grouping structures, for example, a group of co-regulated genes that have similar biological functions tend to have similar expressions. Thus it is preferable to take the grouping structure into consideration to select features. In this paper, we propose a Bayesian Robit regression method with Hyper-LASSO priors (shortened by BayesHL) for feature selection in high dimensional genomic data with grouping structure. The main features of BayesHL include that it discards more aggressively unrelated features than LASSO, and it makes feature selection within groups automatically without a pre-specified grouping structure. We apply BayesHL in gene expression analysis to identify subsets of genes that contribute to the 5-year survival outcome of endometrial cancer (EC) patients. Results show that BayesHL outperforms alternative methods (including LASSO, group LASSO, supervised group LASSO, penalized logistic regression, random forest, neural network, XGBoost and knockoff) in terms of predictive power, sparsity and the ability to uncover grouping structure, and provides insight into the mechanisms of multiple genetic pathways leading to differentiated EC survival outcome.

A novel statistical method for modeling covariate effects in bisulfite sequencing derived measures of DNA methylation

Identifying disease-associated changes in DNA methylation can help us gain a better understanding of disease etiology. Bisulfite sequencing allows the generation of high-throughput methylation profiles at single-base resolution of DNA. However, optimally modeling and analyzing these sparse and discrete sequencing data is still very challenging due to variable read depth, missing data patterns, long-range correlations, data errors, and confounding from cell type mixtures. We propose a regression-based hierarchical model that allows covariate effects to vary smoothly along genomic positions and we have built a specialized EM algorithm, which explicitly allows for experimental errors and cell type mixtures, to make inference about smooth covariate effects in the model. Simulations show that the proposed method provides accurate estimates of covariate effects and captures the major underlying methylation patterns with excellent power. We also apply our method to analyze data from rheumatoid arthritis patients and controls. The method has been implemented in R package SOMNiBUS.

Simultaneous SNP selection and adjustment for population structure in high dimensional prediction models

Complex traits are known to be influenced by a combination of environmental factors and rare and common genetic variants. However, detection of such multivariate associations can be compromised by low statistical power and confounding by population structure. Linear mixed effects models (LMM) can account for correlations due to relatedness but have not been applicable in high-dimensional (HD) settings where the number of fixed effect predictors greatly exceeds the number of samples. False positives or false negatives can result from two-stage approaches, where the residuals estimated from a null model adjusted for the subjects' relationship structure are subsequently used as the response in a standard penalized regression model. To overcome these challenges, we develop a general penalized LMM with a single random effect called ggmix for simultaneous SNP selection and adjustment for population structure in high dimensional prediction models. We develop a blockwise coordinate descent algorithm with automatic tuning parameter selection which is highly scalable, computationally efficient and has theoretical guarantees of convergence. Through simulations and three real data examples, we show that ggmix leads to more parsimonious models compared to the two-stage approach or principal component adjustment with better prediction accuracy. Our method performs well even in the presence of highly correlated markers, and when the causal SNPs are included in the kinship matrix. ggmix can be used to construct polygenic risk scores and select instrumental variables in Mendelian randomization studies. Our algorithms are available in an R package available on CRAN (https://cran.r-project.org/package=ggmix).

MDiNE: a model to estimate differential co-occurrence networks in microbiome studies

MOTIVATION

The human microbiota is the collection of microorganisms colonizing the human body, and plays an integral part in human health. A growing trend in microbiome analysis is to construct a network to estimate the co-occurrence patterns among taxa through precision matrices. Existing methods do not facilitate investigation into how these networks change with respect to covariates.

RESULTS

We propose a new model called Microbiome Differential Network Estimation (MDiNE) to estimate network changes with respect to a binary covariate. The counts of individual taxa in the samples are modeled through a multinomial distribution whose probabilities depend on a latent Gaussian random variable. A sparse precision matrix over all the latent terms determines the co-occurrence network among taxa. The model fit is obtained and evaluated using Hamiltonian Monte Carlo methods. The performance of our model is evaluated through an extensive simulation study and is shown to outperform existing methods in terms of estimation of network parameters. We also demonstrate an application of the model to estimate changes in the intestinal microbial network topology with respect to Crohn's disease.

AVAILABILITY AND IMPLEMENTATION

MDiNE is implemented in a freely available R package: https://github.com/kevinmcgregor/mdine.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Polygenic risk for coronary heart disease acts through atherosclerosis in type 2 diabetes

Background

Type 2 diabetes increases the risk of coronary heart disease (CHD), yet the mechanisms involved remain poorly described. Polygenic risk scores (PRS) provide an opportunity to understand risk factors since they reflect etiologic pathways from the entire genome. We therefore tested whether a PRS for CHD influenced risk of CHD in individuals with type 2 diabetes and which risk factors were associated with this PRS.

Methods

We tested the association of a CHD PRS with CHD and its traditional clinical risk factors amongst individuals with type 2 diabetes in UK Biobank (N = 21,102). We next tested the association of the CHD PRS with atherosclerotic burden in a cohort of 352 genome-wide genotyped participants with type 2 diabetes who had undergone coronary angiograms.

Results

In the UK Biobank we found that the CHD PRS was strongly associated with CHD amongst individuals with type 2 diabetes (OR per standard deviation increase = 1.50; p = 1.5 × 10^− 59). But this CHD PRS was, at best, only weakly associated with traditional clinical risk factors, such as hypertension, hyperlipidemia, glycemic control, obesity and smoking. Conversely, in the angiographic cohort, the CHD PRS was strongly associated with multivessel stenosis (OR = 1.65; p = 4.9 × 10^− 4) and increased number of major stenotic lesions (OR = 1.35; p = 9.4 × 10^− 3).

Conclusions

Polygenic predisposition to CHD is strongly associated with atherosclerotic burden in individuals with type 2 diabetes and this effect is largely independent of traditional clinical risk factors. This suggests that genetic risk for CHD acts through atherosclerosis with little effect on most traditional risk factors, providing the opportunity to explore new biological pathways.

Suicidal ideation is common in autosomal dominant Alzheimer's disease at-risk persons

OBJECTIVES

To study the frequency of suicidal ideation and its association with clinical and neurobiological correlates among cognitively intact autosomal dominant Alzheimer's disease (ADAD) at-risk individuals.

METHODS/DESIGN

In a cross-sectional study of 183 ADAD at-risk individuals (91 mutation carriers and 92 noncarriers), we compared the frequency of suicidal ideation among carriers and noncarriers. Linear mixed-effects models with family-level random effects evaluated the relationships between geriatric depression scale (GDS), neuropsychiatric inventory-questionnaire (NPI-Q), and suicidal ideation scores among all ADAD at-risk individuals. An interaction term was added to the regression models to evaluate the interactions of suicidal ideation and mutation status on neuropsychiatric symptoms.

RESULTS

Twenty-six (14.20%) ADAD at-risk individuals (13 [14.28%] carriers and 13 [14.13%] noncarriers) had suicidal ideation. The frequency of suicidal ideation did not differ between carriers and noncarriers. Suicidal ideation was associated with higher GDS among all ADAD at-risk individuals. When stratified into mutation carrier status, noncarriers with suicidal ideation had higher GDS than carriers. There was no statistically significant association between suicidal ideation and NPI-Q among ADAD at-risk individuals. Awareness of mutation status, neuropsychological performances, and cerebrospinal fluid AD biomarkers were not associated with suicidal ideation among carriers and noncarriers.

CONCLUSIONS

Suicidal ideation is common among cognitively intact ADAD at-risk individuals. While ADAD at-risk individuals with suicidal ideation have greater depressive symptoms, noncarriers with suicidal ideation have higher GDS scores than carriers. Interestingly, awareness of the mutation status was not associated with suicidal ideation in our study. Early identification of suicidal thoughts can facilitate timely interventions to prevent suicidal behaviours. Keywords autosomal dominant Alzheimer's diseasedominantly inherited Alzheimer's networkneuropsychiatric symptomssuicidal ideation.