Statistical power and significance testing in large-scale genetic studies

Hepatic WDR23 proteostasis mediates insulin homeostasis by regulating insulin-degrading enzyme capacity

Maintaining insulin homeostasis is critical for cellular and organismal metabolism. In the liver, insulin is degraded by the activity of the insulin-degrading enzyme (IDE). Here, we establish a hepatic regulatory axis for IDE through WDR23-proteostasis. Wdr23KO mice have increased IDE expression, reduced circulating insulin, and defective insulin responses. Genetically engineered human cell models lacking WDR23 also increase IDE expression and display dysregulated phosphorylation of insulin signaling cascade proteins, IRS-1, AKT2, MAPK, FoxO, and mTOR, similar to cells treated with insulin, which can be mitigated by chemical inhibition of IDE. Mechanistically, the cytoprotective transcription factor NRF2, a direct target of WDR23-Cul4 proteostasis, mediates the enhanced transcriptional expression of IDE when WDR23 is ablated. Moreover, an analysis of human genetic variation in WDR23 across a large naturally aging human cohort in the US Health and Retirement Study reveals a significant association of WDR23 with altered hemoglobin A1C (HbA1c) levels in older adults, supporting the use of WDR23 as a new molecular determinant of metabolic health in humans.

Sexually selected differences in warbler plumage are related to a putative inversion on the Z chromosome

Large structural variants in the genome, such as inversions, may play an important role in producing population structure and local adaptation to the environment through suppression of recombination. However, relatively few studies have linked inversions to phenotypic traits that are sexually selected and may play a role in reproductive isolation. Here, we found that geographic differences in the sexually selected plumage of a warbler, the common yellowthroat (Geothlypis trichas), are largely due to differences in the Z (sex) chromosome (males are ZZ), which contains at least one putative inversion spanning 40% (31/77 Mb) of its length. The inversions on the Z chromosome vary dramatically east and west of the Appalachian Mountains, which provides evidence of cryptic population structure within the range of the most widespread eastern subspecies (G. t. trichas). In an eastern (New York) and western (Wisconsin) population of this subspecies, female prefer different male ornaments; larger black facial masks are preferred in Wisconsin and larger yellow breasts are preferred in New York. The putative inversion also contains genes related to vision, which could influence mating preferences. Thus, structural variants on the Z chromosome are associated with geographic differences in male ornaments and female choice, which may provide a mechanism for maintaining different patterns of sexual selection in spite of gene flow between populations of the same subspecies.

Influence of Clinical and Genetic Factors on Propofol Dose Requirements: A Genome-wide Association Study

BACKGROUND

Propofol is a widely used intravenous hypnotic. Dosing is based mostly on weight, with great interindividual variation in consumption. Suggested factors affecting propofol requirements include age, sex, ethnicity, anxiety, alcohol consumption, smoking, and concomitant valproate use. Genetic factors have not been widely explored.

METHODS

This study considered 1,000 women undergoing breast cancer surgery under propofol and remifentanil anesthesia. Depth of anesthesia was monitored with State Entropy (GE Healthcare, Finland). Propofol requirements during surgery were recorded. DNA from blood was genotyped with a genome-wide array. A multivariable linear regression model was used to assess the relevance of clinical variables and select those to be used as covariates in a genome-wide association study. Imputed genotype data were used to explore selected loci further. In silico functional annotation was used to explore possible consequences of the discovered genetic variants. Additionally, previously reported genetic associations from candidate gene studies were tested.

RESULTS

Body mass index, smoking status, alcohol use, remifentanil dose (ln[mg · kg-1 · min-1]), and average State Entropy during surgery remained statistically significant in the multivariable model. Two loci reached genome-wide significance (P < 5 × 10-8). The most significant associations were for single-nucleotide polymorphisms rs997989 (30 kb from ROBO3), likely affecting expression of another nearby gene, FEZ1, and rs9518419, close to NALCN (sodium leak channel); rs10512538 near KCNJ2 encoding the Kir2.1 potassium channel showed suggestive association (P = 4.7 × 10-7). None of these single-nucleotide polymorphisms are coding variants but possibly affect the regulation of nearby genes. None of the single-nucleotide polymorphisms previously reported as affecting propofol pharmacokinetics or pharmacodynamics showed association in the data.

CONCLUSIONS

In this first genome-wide association study exploring propofol requirements, This study discovered novel genetic associations suggesting new biologically relevant pathways for propofol and general anesthesia. The roles of the gene products of ROBO3/FEZ1, NALCN, and KCNJ2 in propofol anesthesia warrant further studies.

EDITOR’S PERSPECTIVE

Enhancing portability of trans-ancestral polygenic risk scores through tissue-specific functional genomic data integration

Portability of trans-ancestral polygenic risk scores is often confounded by differences in linkage disequilibrium and genetic architecture between ancestries. Recent literature has shown that prioritizing GWAS SNPs with functional genomic evidence over strong association signals can improve model portability. We leveraged three RegulomeDB-derived functional regulatory annotations-SURF, TURF, and TLand-to construct polygenic risk models across a set of quantitative and binary traits highlighting functional mutations tagged by trait-associated tissue annotations. Tissue-specific prioritization by TURF and TLand provide a significant improvement in model accuracy over standard polygenic risk score (PRS) models across all traits. We developed the Trans-ancestral Iterative Tissue Refinement (TITR) algorithm to construct PRS models that prioritize functional mutations across multiple trait-implicated tissues. TITR-constructed PRS models show increased predictive accuracy over single tissue prioritization. This indicates our TITR approach captures a more comprehensive view of regulatory systems across implicated tissues that contribute to variance in trait expression.

Assessing the lack of diversity in genetics research across neurodegenerative diseases: A systematic review of the GWAS Catalog and literature

The under-representation of non-European cohorts in neurodegenerative disease genome-wide association studies (GWAS) hampers precision medicine efforts. Despite the inherent genetic and phenotypic diversity in these diseases, GWAS research consistently exhibits a disproportionate emphasis on participants of European ancestry. This study reviews GWAS up to 2022, focusing on non-European or multi-ancestry neurodegeneration studies. We conducted a systematic review of GWAS results and publications up to 2022, focusing on non-European or multi-ancestry neurodegeneration studies. Rigorous article inclusion and quality assessment methods were employed. Of 123 neurodegenerative disease (NDD) GWAS reviewed, 82% predominantly featured European ancestry participants. A single European study identified over 90 risk loci, compared to a total of 50 novel loci in identified in all non-European or multi-ancestry studies. Notably, only six of the loci have been replicated. The significant under-representation of non-European ancestries in NDD GWAS hinders comprehensive genetic understanding. Prioritizing genomic diversity in future research is crucial for advancing NDD therapies and understanding. HIGHLIGHTS: Eighty-two percent of neurodegenerative genome-wide association studies (GWAS) focus on Europeans. Only 6 of 50 novel neurodegenerative disease (NDD) genetic loci have been replicated. Lack of diversity significantly hampers understanding of NDDs. Increasing diversity in NDD genetic research is urgently required. New initiatives are aiming to enhance diversity in NDD research.

Using Alternative Definitions of Controls to Increase Statistical Power in GWAS

Genome-wide association studies (GWAS) are often underpowered due to small effect sizes of common single nucleotide polymorphisms (SNPs) on phenotypes and extreme multiple testing thresholds. The most common approach for increasing statistical power is to increase sample size. We propose an alternative strategy of redefining case-control outcomes into ordinal case-subthreshold-asymptomatic variables. While maintaining the clinical case threshold, we subdivide controls into two groups: individuals who are symptomatic but do not meet the clinical criteria for diagnosis (subthreshold) and individuals who are effectively asymptomatic. We conducted a simulation study to examine the impact of effect size, minor allele frequency, population prevalence, and the prevalence of the subthreshold group on statistical power to detect genetic associations in three scenarios: a standard case-control, an ordinal, and a case-asymptomatic control analysis. Our results suggest the ordinal model consistently provides the greatest statistical power while the case-control model the least. Power in the case-asymptomatic control model reflects the case-control or ordinal model depending on the population prevalence and size of the subthreshold category. We then analyzed a major depression phenotype from the UK Biobank to corroborate our simulation results. Overall, the ordinal model improves statistical power in GWAS consistent with increasing the sample size by approximately 10%.

Genetic parameter estimation and molecular foundation of chicken beak shape

The bird beak is mainly functioned as feeding and attacking, and its shape has extremely important significance for survival and reproduction. In chickens, since beak shape could lead to some disadvantages including pecking and waste of feed, it is important to understand the inheritance of chicken beak shape. In the present study, we firstly established 4 indicators to describe the chicken beak shapes, including upper beak length (UL), lower beak length (LL), distance between upper and lower beak tips (DB) and upper beak curvature (BC). And then, we measured the 4 beak shape indicators as well as some production traits including body weight (BW), shank length (SL), egg weight (EW), eggshell strength (ES) of a layer breed, Rhode Island Red (RIR), in order to estimate genetic parameters of chicken beak shape. The heritabilities of UL and LL were 0.41 and 0.37, and the heritabilities of DB and BC were 0.22 and 0.21, indicating that beak shape was a highly or mediumly heritable. There were significant positive genetic and phenotypic correlations among UL, LL, and DB. And UL was positively correlated with body weight (BW18) and shank length (SL18) at 18 weeks of age in genetics, and DB was positively correlated with BC in terms of genetics and phenotype. We also found that layers of chicken cages played a role on beak shape, which could be attributed to the difference of lightness in different cage layers. By a genome-wide association study (GWAS) for the chicken UL, we identified 9 significant candidate genes associated with UL in RIR. For the variants with low minor allele frequencies (MAF <0.01) and outside of high linkage disequilibrium (LD) regions, we also conducted rare variant association studies (RVA) and GWAS to find the association between genotype and phenotype. We also analyzed transcriptomic data from multiple tissues of chicken embryos and revealed that all of the 9 genes were highly expressed in beak of chicken embryos, indicating their potential function for beak development. Our results provided the genetic foundation of chicken beak shape, which could help chicken breeding on beak related traits.

Exploiting integrative metabolomics to study host-parasite interactions in Plasmodium infections

Despite years of research, malaria remains a significant global health burden, with poor diagnostic tests and increasing antimalarial drug resistance challenging diagnosis and treatment. While 'single-omics'-based approaches have been instrumental in gaining insight into the biology and pathogenicity of the Plasmodium parasite and its interaction with the human host, a more comprehensive understanding of malaria pathogenesis can be achieved through 'multi-omics' approaches. Integrative methods, which combine metabolomics, lipidomics, transcriptomics, and genomics datasets, offer a holistic systems biology approach to studying malaria. This review highlights recent advances, future directions, and challenges involved in using integrative metabolomics approaches to interrogate the interactions between Plasmodium and the human host, paving the way towards targeted antimalaria therapeutics and control intervention methods.

Causal relationships between Sjögren's syndrome and Parkinson's disease: A Mendelian randomization study

BACKGROUND

Epidemiological and observational studies have indicated an association between Sjögren's syndrome (SS) and Parkinson's disease (PD). However, consistent conclusions have not been reached due to various limitations. In order to determine whether SS and PD are causally related, we conducted a Mendelian randomization study (MR) with two samples.

METHODS

Data for SS derived from the FinnGen consortium's R9 release (2495 cases and 365 533 controls). Moreover, data for PD were acquired from the publicly available GWAS of European ancestry, which involved 33 674 cases and 449 056 controls. The inverse variance weighted, along with four other effective methodologies, were employed to comprehensively infer the causal relationships between SS and PD. To assess the estimation's robustness, a number of sensitivity studies were performed. To determine the probability of reverse causality, we performed a reverse MR analysis.

RESULTS

There was no evidence of a significant causal effect of SS on PD risks based on the MR [odds ratio (OR) = 1.03; 95% confidence interval (CI) = 0.95-1.11; p = .45]. Similarly, no evidence supported the causal effects of PD on SS (OR = 0.92; 95% CI = 0.81-1.04; p = .20). These findings held up under rigorous sensitivity analysis.

CONCLUSIONS

MR bidirectional analysis did not reveal any cause-and-effect relationship between SS and PD, or vice versa. Further study of the mechanisms that may underlie the probable causal association between SS and PD is needed.

Maternal metabolic profiling across body mass index groups: An exploratory longitudinal study

INTRODUCTION

Increased BMI has been identified as a risk factor for most pregnancy complications, but the underlying metabolic factors mediating the detrimental effects of BMI are largely unknown. We aimed to compare metabolic profiles in overweight/obese women (body mass index [BMI] ≥ 25 kg/m2 ) and normal weight/underweight women (BMI < 25 kg/m2 ) across gestation. We also explored how gestational weight gain (GWG) affected maternal metabolic profiles.

MATERIAL AND METHODS

Exploratory nested case-control study based on a prospective longitudinal cohort of women who were healthy prior to pregnancy and gave birth at Oslo University Hospital from 2002 to 2008. The sample consisted of 48 women who were overweight/obese and 59 normal-weight/underweight women. Plasma samples from four time points in pregnancy (weeks 14-16, 22-24, 30-32 and 36-38) were analyzed by nuclear magnetic resonance spectroscopy and 91 metabolites were measured. Linear regression models were fitted for each of the metabolites at each time point.

RESULTS

Overweight or obese women had higher levels of lipids in very-low-density lipoprotein (VLDL), total triglycerides, triglycerides in VLDL, total fatty acids, monounsaturated fatty acids, saturated fatty acids, leucine, valine, and total branched-chain amino acids in pregnancy weeks 14-16 compared to underweight and normal-weight women. Docosahexaenoic acid and degree of unsaturation were significantly lower in overweight/obese women in pregnancy weeks 36-38. In addition, overweight or obese women had higher particle concentration of XXL-VLDL and glycoprotein acetyls (GlycA) at weeks 14-16 and 30-32. GWG did not seem to affect the metabolic profile, regardless of BMI group when BMI was treated as a dichotomous variable, ≥25 kg/m2 (yes/no).

CONCLUSIONS

Overweight or obese women had smaller pregnancy-related metabolic alterations than normal-weight/underweight women. There was a trend toward higher triglyceride and VLDL particle concentration in overweight/obese women. As this was a hypothesis-generating study, the similarities with late-onset pre-eclampsia warrant further investigation. The unfavorable development of fatty acid composition in overweight/obese women, with possible implication for the offspring, should also be studied further in the future.

Genetic Factors and Long-term Treatment-Related Neurocognitive Deficits, Anxiety, and Depression in Childhood Leukemia Survivors: An Exome-Wide Association Study

BACKGROUND

An increased risk of neurocognitive deficits, anxiety, and depression has been reported in childhood cancer survivors.

METHODS

We analyzed associations of neurocognitive deficits, as well as anxiety and depression, with common and rare genetic variants derived from whole-exome sequencing data of acute lymphoblastic leukemia (ALL) survivors from the PETALE cohort. In addition, significant associations were assessed using stratified and multivariable analyses. Next, top-ranking common associations were analyzed in an independent SJLIFE replication cohort of ALL survivors.

RESULTS

Significant associations were identified in the entire discovery cohort (N = 229) between the AK8 gene and changes in neurocognitive function, whereas PTPRZ1, MUC16, TNRC6C-AS1 were associated with anxiety. Following stratification according to sex, the ZNF382 gene was linked to a neurocognitive deficit in males, whereas APOL2 and C6orf165 were associated with anxiety and EXO5 with depression. Following stratification according to prognostic risk groups, the modulatory effect of rare variants on depression was additionally found in the CYP2W1 and PCMTD1 genes. In the replication SJLIFE cohort (N = 688), the male-specific association in the ZNF382 gene was not significant; however, a P value<0.05 was observed when the entire SJLIFE cohort was analyzed. ZNF382 was significant in males in the combined cohorts as shown by meta-analyses as well as the depression-associated gene EXO5.

CONCLUSIONS

Further research is needed to confirm whether the current findings, along with other known risk factors, may be valuable in identifying patients at increased risk of these long-term complications.

IMPACT

Our results suggest that specific genes may be related to increased neuropsychological consequences.

Causal relationship between gut microbiota and polycystic ovary syndrome: a literature review and Mendelian randomization study

Introduction

Numerous studies have suggested an association between gut microbiota and polycystic ovarian syndrome (PCOS). However, the causal relationship between these two factors remains unclear.

Methods

A review of observational studies was conducted to compare changes in gut microbiota between PCOS patients and controls. The analysis focused on four levels of classification, namely, phylum, family, genus, and species/genus subgroups. To further investigate the causal relationship, Mendelian randomization (MR) was employed using genome-wide association study (GWAS) data on gut microbiota from the MiBioGen consortium, as well as GWAS data from a large meta-analysis of PCOS. Additionally, a reverse MR was performed, and the results were verified through sensitivity analyses.

Results

The present review included 18 observational studies that met the inclusion and exclusion criteria. The abundance of 64 gut microbiota taxa significantly differed between PCOS patients and controls. Using the MR method, eight bacteria were identified as causally associated with PCOS. The protective effects of the genus Sellimonas on PCOS remained significant after applying Bonferroni correction. No significant heterogeneity or horizontal pleiotropy was found in the instrumental variables (IVs). Reverse MR analyses did not reveal a significant causal effect of PCOS on gut microbiota.

Conclusion

The differences in gut microbiota between PCOS patients and controls vary across observational studies. However, MR analyses identified specific gut microbiota taxa that are causally related to PCOS. Future studies should investigate the gut microbiota that showed significant results in the MR analyses, as well as the underlying mechanisms of this causal relationship and its potential clinical significance.

Using Alternative Definitions of Controls to Increase Statistical Power in GWAS

Genome-wide association studies (GWAS) are underpowered due to small effect sizes of single nucleotide polymorphisms (SNPs) on phenotypes and extreme multiple testing thresholds. The most common approach for increasing statistical power is to increase sample size. We propose an alternative strategy of redefining case-control outcomes into ordinal case-subthreshold-asymptomatic variables. While maintaining the clinical case threshold, we subdivide controls into two groups: individuals who are symptomatic but do not meet the clinical criteria for diagnosis (subthreshold) and individuals who are effectively asymptomatic. We conducted a simulation study to examine the impact of effect size, minor allele frequency, population prevalence, and the prevalence of the subthreshold group on statistical power to detect genetic associations in three scenarios: a standard case-control, an ordinal, and a case-asymptomatic control analysis. Our results suggest the ordinal model consistently provides the most statistical power while the case-control model the least. Power in the case-asymptomatic control model reflects the case-control or ordinal model depending on the population prevalence and size of the subthreshold category. We then analyzed a major depression phenotype from the UK Biobank to corroborate our simulation results. Overall, the ordinal model improves statistical power in GWAS consistent with increasing the sample size by approximately 10%.

Assessing the lack of diversity in genetics research across neurodegenerative diseases: a systematic review of the GWAS Catalog and literature

Importance

The under-representation of participants with non-European ancestry in genome-wide association studies (GWAS) is a critical issue that has significant implications, including hindering the progress of precision medicine initiatives. This issue is particularly significant in the context of neurodegenerative diseases (NDDs), where current therapeutic approaches have shown limited success. Addressing this under-representation is crucial to harnessing the full potential of genomic medicine in underserved communities and improving outcomes for NDD patients.

Objective

Our primary objective was to assess the representation of non-European ancestry participants in genetic discovery efforts related to NDDs. We aimed to quantify the extent of inclusion of diverse ancestry groups in NDD studies and determine the number of associated loci identified in more inclusive studies. Specifically, we sought to highlight the disparities in research efforts and outcomes between studies predominantly involving European ancestry participants and those deliberately targeting non-European or multi-ancestry populations across NDDs.

Evidence Review

We conducted a systematic review utilizing existing GWAS results and publications to assess the inclusion of diverse ancestry groups in neurodegeneration and neurogenetics studies. Our search encompassed studies published up to the end of 2022, with a focus on identifying research that deliberately included non-European or multi-ancestry cohorts. We employed rigorous methods for the inclusion of identified articles and quality assessment.

Findings

Our review identified a total of 123 NDD GWAS. Strikingly, 82% of these studies predominantly featured participants of European ancestry. Endeavors specifically targeting non-European or multi-ancestry populations across NDDs identified only 52 risk loci. This contrasts with predominantly European studies, which reported over 90 risk loci for a single disease. Encouragingly, over 65% of these discoveries occurred in 2020 or later, indicating a recent increase in studies deliberately including non-European cohorts.

Conclusions and relevance

Our findings underscore the pressing need for increased diversity in neurodegenerative research. The significant under-representation of non-European ancestry participants in NDD GWAS limits our understanding of the genetic underpinnings of these diseases. To advance the field of neurodegenerative research and develop more effective therapies, it is imperative that future investigations prioritize and harness the genomic diversity present within and across global populations.

Over-representation analysis of angiogenic factors in immunosuppressive mechanisms in neoplasms and neurological conditions during COVID-19

BACKGROUND

Recent studies emphasized the necessity to identify key (human) biological processes and pathways targeted by the Coronaviridae family of viruses, especially Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Coronavirus Disease (COVID-19) caused up to 33-55 % death rates in COVID-19 patients with malignant neoplasms and Alzheimer's disease. Given this scenario, we identified biological processes and pathways involved in various diseases which are most likely affected by COVID-19.

METHODS

The COVID-19 DisGeNET data set (v4.0) contains the associations between various diseases and human genes known to interact with viruses from Coronaviridae family and were obtained from the IntAct Coronavirus data set annotated with DisGeNET data. We constructed the disease-gene network to identify genes that are involved in various comorbid diseased states. Communities from the disease-gene network were identified using Louvain method and functional enrichment through over-representation analysis methodology was used to discover significant biological processes and pathways shared between COVID-19 and other diseases.

RESULT

The COVID-19 DisGeNET data set (v4.0) comprised of 828 human genes and 10,473 diseases (including various phenotypes) that together constituted nodes in the disease-gene network. Each of the 70,210 edges connects a human gene with an associated disease. The top 10 genes linked to most number of diseases were VEGFA, BCL2, CTNNB1, ALB, COX2, AGT, HLA-A, HMOX1, FGF2 and COMT. The most vulnerable group of patients thus discovered had comorbid conditions such as carcinomas, malignant neoplasms and Alzheimer's disease. Finally, we identified 15 potentially useful biological processes and pathways for improved therapies. Vascular endothelial growth factor (VEGF) is the key mediator of angiogenesis in cancer. It is widely distributed in the brain and plays a crucial role in brain inflammation regulating the level of angiopoietins. With a degree of 1899, VEGFA was associated with maximum number of diseases in the disease-gene network. Previous studies have indicated that increased levels of VEGFA in the blood results in dyspnea, Pulmonary Edema (PE), Acute Lung Injury (ALI) and Acute Respiratory Distress Syndrome (ARDS). In case of COVID-19 patients with neoplasms and other neurological symptoms, our results indicate VEGFA as a therapeutic target for inflammation suppression. As VEGFs are known to disproportionately affect cancer patients, improving endothelial permeability and vasodilation with anti-VEGF therapy could lead to suppression of inflammation and also improve oxygenation. As an outcome of our study, we make case for clinical investigations towards anti-VEGF therapies for such comorbid conditions affected by COVID-19 for better therapeutic outcomes.

Genomic Markers Associated with Cytomegalovirus DNAemia in Kidney Transplant Recipients

Human cytomegalovirus (CMV) is a major pathogen after solid organ transplantation, leading to high morbidity and mortality. Transplantation from a CMV-seropositive donor to a CMV-seronegative recipient (D+/R-) is associated with high risk of CMV disease. However, that risk is not uniform, suggesting a role for host factors in immune control of CMV. To identify host genetic factors that control CMV DNAemia post transplantation, we performed a whole-exome association study in two cohorts of D+/R- kidney transplant recipients. Quantitative CMV DNA was measured for at least one year following transplantation. Several CMV-protective single-nucleotide polymorphisms (SNPs) were identified in the first cohort (72 patients) but were not reproducible in the second cohort (126 patients). A meta-analysis of both cohorts revealed several SNPs that were significantly associated with protection from CMV DNAemia. The copy number variation of several genes was significantly different between recipients with and without CMV DNAemia. Amongst patients with CMV DNAemia in the second cohort, several variants of interest (p < 5 × 10-5), the most common of which was NLRC5, were associated with peak viral load. We provide new predictive genetic markers for protection of CMV DNAemia. These markers should be validated in larger cohorts.

From viral load to survival: Unveiling new genetic links for resistance against PMCV in Atlantic Salmon (Salmo salar)

The infectious agent piscine myocarditis virus (PMCV) causes cardiomyopathy syndrome (CMS) and is responsible for substantial mortality and economic losses in the Atlantic salmon (Salmo salar) farming industry. Previous research has demonstrated that breeding for resistance against PMCV is an effective approach to mitigate the disease's impact. In this study, a new quantitative trait locus (QTL) is described on chromosome 23, together with previously described QTLs on chromosomes 12 and 27. The findings are based on two genome-wide association studies conducted on two different year-classes of Atlantic salmon of the Rauma strain. In this study, we utilized data from an experimental challenge trial with the viral load as the phenotype and a field outbreak of CMS with survival data as the phenotype. The estimated SNP-based heritability was 0.55 and 0.44 in the two studies, respectively. In the infection trial, the top associated SNP on chromosome 23 accounted for approximately 46% of the genetic and 25.53% of the phenotypic variations in the viral load. In the field outbreak, we identified a QTL on the same genomic region of chromosome 23. The most significantly associated marker on this chromosome explained 13.57% and 5.97% of the genetic and phenotypic variations. The QTL on chromosome 23 is in proximity to delta-5 fatty acyl desaturase and fatty acid desaturase 2 genes, both of which play a role in the production of polyunsaturated fatty acids. This proximity is particularly interesting as it offers valuable insights into enhancing our understanding of resistance against PMCV.

Cross-ancestry analyses identify new genetic loci associated with 25-hydroxyvitamin D

Vitamin D status-a complex trait influenced by environmental and genetic factors-is tightly associated with skin colour and ancestry. Yet very few studies have investigated the genetic underpinnings of vitamin D levels across diverse ancestries, and the ones that have, relied on small sample sizes, resulting in inconclusive results. Here, we conduct genome-wide association studies (GWAS) of 25 hydroxyvitamin D (25OHD)-the main circulating form of vitamin D-in 442,435 individuals from four broad genetically-determined ancestry groups represented in the UK Biobank: European (N = 421,867), South Asian (N = 9,983), African (N = 8,306) and East Asian (N = 2,279). We identify a new genetic determinant of 25OHD (rs146759773) in individuals of African ancestry, which was not detected in previous analysis of much larger European cohorts due to low minor allele frequency. We show genome-wide significant evidence of dominance effects in 25OHD that protect against vitamin D deficiency. Given that key events in the synthesis of 25OHD occur in the skin and are affected by pigmentation levels, we conduct GWAS of 25OHD stratified by skin colour and identify new associations. Lastly, we test the interaction between skin colour and variants associated with variance in 25OHD levels and identify two loci (rs10832254 and rs1352846) whose association with 25OHD differs in individuals of distinct complexions. Collectively, our results provide new insights into the complex relationship between 25OHD and skin colour and highlight the importance of diversity in genomic studies. Despite the much larger rates of vitamin D deficiency that we and others report for ancestry groups with dark skin (e.g., South Asian), our study highlights the importance of considering ancestral background and/or skin colour when assessing the implications of low vitamin D.

Separating Measurement Error and Signal in Environmental Data: Use of Replicates to Address Uncertainty

Measurement uncertainty has long been a concern in the characterizing and interpreting environmental and toxicological measurements. We compared statistical analysis approaches when there are replicates: a Naı̈ve approach that omits replicates, a Hybrid approach that inappropriately treats replicates as independent samples, and a Measurement Error Model (MEM) approach in a random effects analysis of variance (ANOVA) model that appropriately incorporates replicates. A simulation study assessed the effects of sample size and levels of replication, signal variance, and measurement error on estimates from the three statistical approaches. MEM results were superior overall with confidence intervals for the observed mean narrower on average than those from the Naı̈ve approach, giving improved characterization. The MEM approach also featured an unparalleled advantage in estimating signal and measurement error variance separately, directly addressing measurement uncertainty. These MEM estimates were approximately unbiased on average with more replication and larger sample sizes. Case studies illustrated analyzing normally distributed arsenic and log-normally distributed chromium concentrations in tap water and calculating MEM confidence intervals for the true, latent signal mean and latent signal geometric mean (i.e., with measurement error removed). MEM estimates are valuable for study planning; we used simulation to compare various sample sizes and levels of replication.

Large-scale plasma proteomics comparisons through genetics and disease associations

High-throughput proteomics platforms measuring thousands of proteins in plasma combined with genomic and phenotypic information have the power to bridge the gap between the genome and diseases. Here we performed association studies of Olink Explore 3072 data generated by the UK Biobank Pharma Proteomics Project1 on plasma samples from more than 50,000 UK Biobank participants with phenotypic and genotypic data, stratifying on British or Irish, African and South Asian ancestries. We compared the results with those of a SomaScan v4 study on plasma from 36,000 Icelandic people2, for 1,514 of whom Olink data were also available. We found modest correlation between the two platforms. Although cis protein quantitative trait loci were detected for a similar absolute number of assays on the two platforms (2,101 on Olink versus 2,120 on SomaScan), the proportion of assays with such supporting evidence for assay performance was higher on the Olink platform (72% versus 43%). A considerable number of proteins had genomic associations that differed between the platforms. We provide examples where differences between platforms may influence conclusions drawn from the integration of protein levels with the study of diseases. We demonstrate how leveraging the diverse ancestries of participants in the UK Biobank helps to detect novel associations and refine genomic location. Our results show the value of the information provided by the two most commonly used high-throughput proteomics platforms and demonstrate the differences between them that at times provides useful complementarity.

Critical assessment of on-premise approaches to scalable genome analysis

BACKGROUND

Plummeting DNA sequencing cost in recent years has enabled genome sequencing projects to scale up by several orders of magnitude, which is transforming genomics into a highly data-intensive field of research. This development provides the much needed statistical power required for genotype-phenotype predictions in complex diseases.

METHODS

In order to efficiently leverage the wealth of information, we here assessed several genomic data science tools. The rationale to focus on on-premise installations is to cope with situations where data confidentiality and compliance regulations etc. rule out cloud based solutions. We established a comprehensive qualitative and quantitative comparison between BCFtools, SnpSift, Hail, GEMINI, and OpenCGA. The tools were compared in terms of data storage technology, query speed, scalability, annotation, data manipulation, visualization, data output representation, and availability.

RESULTS

Tools that leverage sophisticated data structures are noted as the most suitable for large-scale projects in varying degrees of scalability in comparison to flat-file manipulation (e.g., BCFtools, and SnpSift). Remarkably, for small to mid-size projects, even lightweight relational database.

CONCLUSION

The assessment criteria provide insights into the typical questions posed in scalable genomics and serve as guidance for the development of scalable computational infrastructure in genomics.

ADuLT: An efficient and robust time-to-event GWAS

Proportional hazards models have been proposed to analyse time-to-event phenotypes in genome-wide association studies (GWAS). However, little is known about the ability of proportional hazards models to identify genetic associations under different generative models and when ascertainment is present. Here we propose the age-dependent liability threshold (ADuLT) model as an alternative to a Cox regression based GWAS, here represented by SPACox. We compare ADuLT, SPACox, and standard case-control GWAS in simulations under two generative models and with varying degrees of ascertainment as well as in the iPSYCH cohort. We find Cox regression GWAS to be underpowered when cases are strongly ascertained (cases are oversampled by a factor 5), regardless of the generative model used. ADuLT is robust to ascertainment in all simulated scenarios. Then, we analyse four psychiatric disorders in iPSYCH, ADHD, Autism, Depression, and Schizophrenia, with a strong case-ascertainment. Across these psychiatric disorders, ADuLT identifies 20 independent genome-wide significant associations, case-control GWAS finds 17, and SPACox finds 8, which is consistent with simulation results. As more genetic data are being linked to electronic health records, robust GWAS methods that can make use of age-of-onset information will help increase power in analyses for common health outcomes.

Trumpet plots: visualizing the relationship between allele frequency and effect size in genetic association studies

Recent advances in genome-wide association and sequencing studies have shown that the genetic architecture of complex traits and diseases involves a combination of rare and common genetic variants distributed throughout the genome. One way to better understand this architecture is to visualize genetic associations across a wide range of allele frequencies. However, there is currently no standardized or consistent graphical representation for effectively illustrating these results. Here we propose a standardized approach for visualizing the effect size of risk variants across the allele frequency spectrum. The proposed plots have a distinctive trumpet shape: with the majority of variants having high frequency and small effects, and a small number of variants having lower frequency and larger effects. To demonstrate the utility of trumpet plots in illustrating the relationship between the number of variants, their frequency, and the magnitude of their effects in shaping the genetic architecture of complex traits and diseases, we generated trumpet plots for more than one hundred traits in the UK Biobank. To facilitate their broader use, we developed an R package, 'TrumpetPlots' (available at the Comprehensive R Archive Network) and R Shiny application, 'Shiny Trumpets' (available at https://juditgg.shinyapps.io/shinytrumpets/) that allows users to explore these results and submit their own data.

EpiMix is an integrative tool for epigenomic subtyping using DNA methylation

DNA methylation (DNAme) is a major epigenetic factor influencing gene expression with alterations leading to cancer and immunological and cardiovascular diseases. Recent technological advances have enabled genome-wide profiling of DNAme in large human cohorts. There is a need for analytical methods that can more sensitively detect differential methylation profiles present in subsets of individuals from these heterogeneous, population-level datasets. We developed an end-to-end analytical framework named "EpiMix" for population-level analysis of DNAme and gene expression. Compared with existing methods, EpiMix showed higher sensitivity in detecting abnormal DNAme that was present in only small patient subsets. We extended the model-based analyses of EpiMix to cis-regulatory elements within protein-coding genes, distal enhancers, and genes encoding microRNAs and long non-coding RNAs (lncRNAs). Using cell-type-specific data from two separate studies, we discover epigenetic mechanisms underlying childhood food allergy and survival-associated, methylation-driven ncRNAs in non-small cell lung cancer.

E-value: a superior alternative to P-value and its adjustments in DNA methylation studies

DNA methylation plays a crucial role in transcriptional regulation. Reduced representation bisulfite sequencing (RRBS) is a technique of increasing use for analyzing genome-wide methylation profiles. Many computational tools such as Metilene, MethylKit, BiSeq and DMRfinder have been developed to use RRBS data for the detection of the differentially methylated regions (DMRs) potentially involved in epigenetic regulations of gene expression. For DMR detection tools, as for countless other medical applications, P-values and their adjustments are among the most standard reporting statistics used to assess the statistical significance of biological findings. However, P-values are coming under increasing criticism relating to their questionable accuracy and relatively high levels of false positive or negative indications. Here, we propose a method to calculate E-values, as likelihood ratios falling into the null hypothesis over the entire parameter space, for DMR detection in RRBS data. We also provide the R package 'metevalue' as a user-friendly interface to implement E-value calculations into various DMR detection tools. To evaluate the performance of E-values, we generated various RRBS benchmarking datasets using our simulator 'RRBSsim' with eight samples in each experimental group. Our comprehensive benchmarking analyses showed that using E-values not only significantly improved accuracy, area under ROC curve and power, over that of P-values or adjusted P-values, but also reduced false discovery rates and type I errors. In applications using real RRBS data of CRL rats and a clinical trial on low-salt diet, the use of E-values detected biologically more relevant DMRs and also improved the negative association between DNA methylation and gene expression.

Genome-wide Association Study for AKI

Key Points

Two genetic variants in the DISP1-TLR5 gene locus were associated with risk of AKI. DISP1 and TLR5 were differentially regulated in kidney biopsy tissue from patients with AKI compared with no AKI.

Background

Although common genetic risks for CKD are well established, genetic factors influencing risk for AKI in hospitalized patients are poorly understood.

Methods

We conducted a genome-wide association study in 1369 participants in the Assessment, Serial Evaluation, and Subsequent Sequelae of AKI Study; a multiethnic population of hospitalized participants with and without AKI matched on demographics, comorbidities, and kidney function before hospitalization. We then completed functional annotation of top-performing variants for AKI using single-cell RNA sequencing data from kidney biopsies in 12 patients with AKI and 18 healthy living donors from the Kidney Precision Medicine Project.

Results

No genome-wide significant associations with AKI risk were found in Assessment, Serial Evaluation, and Subsequent Sequelae of AKI (P < 5×10 −8 ). The top two variants with the strongest association with AKI mapped to the dispatched resistance-nodulation-division (RND) transporter family member 1 (DISP1) gene and toll-like receptor 5 (TLR5) gene locus, rs17538288 (odds ratio, 1.55; 95% confidence interval, 1.32 to 182; P = 9.47×10 −8 ) and rs7546189 (odds ratio, 1.53; 95% confidence interval, 1.30 to 1.81; P = 4.60×10 −7 ). In comparison with kidney tissue from healthy living donors, kidney biopsies in patients with AKI showed differential DISP1 expression in proximal tubular epithelial cells (adjusted P = 3.9× 10−2) and thick ascending limb of the loop of Henle (adjusted P = 8.7× 10−3) and differential TLR5 gene expression in thick ascending limb of the loop of Henle (adjusted P = 4.9× 10−30).

Conclusions

AKI is a heterogeneous clinical syndrome with various underlying risk factors, etiologies, and pathophysiology that may limit the identification of genetic variants. Although no variants reached genome-wide significance, we report two variants in the intergenic region between DISP1 and TLR5 , suggesting this region as a novel risk for AKI susceptibility.

Increase in power by obtaining 10 or more controls per case when type-1 error is small in large-scale association studies

BACKGROUND

The rule of thumb that there is little gain in statistical power by obtaining more than 4 controls per case, is based on type-1 error α = 0.05. However, association studies that evaluate thousands or millions of associations use smaller α and may have access to plentiful controls. We investigate power gains, and reductions in p-values, when increasing well beyond 4 controls per case, for small α.

METHODS

We calculate the power, the median expected p-value, and the minimum detectable odds-ratio (OR), as a function of the number of controls/case, as α decreases.

RESULTS

As α decreases, at each ratio of controls per case, the increase in power is larger than for α = 0.05. For α between 10-6 and 10-9 (typical for thousands or millions of associations), increasing from 4 controls per case to 10-50 controls per case increases power. For example, a study with power = 0.2 (α = 5 × 10-8) with 1 control/case has power = 0.65 with 4 controls/case, but with 10 controls/case has power = 0.78, and with 50 controls/case has power = 0.84. For situations where obtaining more than 4 controls per case provides small increases in power beyond 0.9 (at small α), the expected p-value can decrease by orders-of-magnitude below α. Increasing from 1 to 4 controls/case reduces the minimum detectable OR toward the null by 20.9%, and from 4 to 50 controls/case reduces by an additional 9.7%, a result which applies regardless of α and hence also applies to "regular" α = 0.05 epidemiology.

CONCLUSIONS

At small α, versus 4 controls/case, recruiting 10 or more controls/cases can increase power, reduce the expected p-value by 1-2 orders of magnitude, and meaningfully reduce the minimum detectable OR. These benefits of increasing the controls/case ratio increase as the number of cases increases, although the amount of benefit depends on exposure frequencies and true OR. Provided that controls are comparable to cases, our findings suggest greater sharing of comparable controls in large-scale association studies.

Invited review: Good practices in genome-wide association studies to identify candidate sequence variants in dairy cattle

Genotype data from dairy cattle selection programs have greatly facilitated GWAS to identify variants related to economic traits. Results can enhance the accuracy of genomic prediction, analyze more complex models that go beyond additive effects, elucidate the genetic architecture of a trait, and finally, decipher the underlying biology of traits. The entire process, comprising data generation, quality control, statistical analyses, interpretation of association results, and linking results to biology should be designed and executed to minimize the generation of false-positive and false-negative associations and misleading links to biological processes. This review aims to provide general guidelines for data analysis that address data quality control, association tests, adjustment for population stratification, and significance evaluation to improve the reliability of conclusions. We also provide guidance on post-GWAS strategy and the interpretation of results. These guidelines are tailored to dairy cattle, which are characterized by long-range linkage disequilibrium, large half-sib families, and routinely collected phenotypes, requiring different approaches than those applied in human GWAS. We discuss common limitations and challenges that have been overlooked in the analysis and interpretation of GWAS to identify candidate sequence variants in dairy cattle.

Next-Generation Sequencing Data-Based Association Testing of a Group of Genetic Markers for Complex Responses Using a Generalized Linear Model Framework

Association testing has been widely used to study the relationship between genetic variants and phenotypes. Most association testing methods are genotype-based, i.e. first estimate genotype and then regress phenotype on estimated genotype and other variables. Directly testing methods based on next generation sequencing (NGS) data without genotype calling have been proposed and shown advantage over genotype-based methods in the scenarios when genotype calling is not accurate. NGS data-based single-variant testing have been proposed including our previously proposed single-variant testing method, i.e. UNC combo method [1]. NGS data-based group testing methods for continuous phenotype have also been proposed by us using a linear model framework which can handle continuous responses [2]. In this paper, we extend our linear model-based framework to a generalized linear model-based framework so that the methods can handle other types of responses especially binary responses which is commonly-faced in association studies. We have conducted extensive simulation studies to evaluate the performance of different estimators and compare our estimators with their corresponding genotype-based methods. We found that all methods have Type I errors controlled, and our NGS data-based testing methods have better performance than their corresponding genotype-based methods in the literature for other types of responses including binary responses (logistic regression) and count responses (Poisson regression especially when sequencing depth is low. In conclusion, we have extended our previous linear model (LM) framework to a generalized linear model (GLM) framework and derived NGS data-based testing methods for a group of genetic variants. Compared with our previously proposed LM-based methods [2], the new GLM-based methods can handle more complex responses (for example, binary responses and count responses) in addition to continuous responses. Our methods have filled the literature gap and shown advantage over their corresponding genotype-based methods in the literature.

Identification and Genomic Localization of Autosomal sdY Locus in a Population of Atlantic Salmon (Salmo salar)

The determination of sex in salmonid fishes is controlled by genetic mechanisms, with males being the heterogametic sex. The master sex-determining gene, the sexually dimorphic gene on the Y chromosome (sdY), is a conserved gene across various salmonid species. Nevertheless, variations in the genomic location of sdY have been observed both within and between species. Furthermore, different studies have reported discordances in the association between the sdY and the phenotypic gender. While some males seem to lack this locus, there have been reports of females carrying sdY. Although the exact reasons behind this discordance remain under investigation, some recent studies have proposed the existence of an autosomal, non-functional copy of sdY as a potential cause. In this study, we confirmed the presence of this autosomal sdY in the SalmoBreed strain of Atlantic salmon using a genotyping platform through a novel approach that allows for high-throughput screening of a large number of individuals. We further characterized the segregation profile of this locus across families and found the ratio of genetically assigned female-to-male progeny to be in accordance with the expected profile of a single autosomal sdY locus. Additionally, our mapping efforts localized this locus to chromosome 3 and suggested a putative copy on chromosome 6.

Community detection in Epstein-Barr virus associated carcinomas and role of tyrosine kinase in etiological mechanisms for oncogenesis

BACKGROUND

Epstein-Barr virus (EBV) affects more than 90% of global population. The role of the virus in causing infectious mononucleosis (IM) affecting B-cells and epithelial cells and in the development of EBV associated cancers is well documented. Investigating the associated interactions can pave way for the discovery of novel therapeutic targets for EBV associated lymphoproliferative (Burkitt's Lymphoma and Hodgkin's Lymphoma) and non-lymphoproliferative diseases (Gastric cancer and Nasopharyngeal cancer).

METHODS

Based on the DisGeNET (v7.0) data set, we constructed a disease-gene network to identify genes that are involved in various carcinomas, viz. Gastric cancer (GC), Nasopharyngeal cancer (NPC), Hodgkin's lymphoma (HL) and Burkitt's lymphoma (BL). We identified communities in the disease-gene network and performed functional enrichment using over-representation analysis to detect significant biological processes/pathways and the interactions between them.

RESULT

We identified the modular communities to explore the relation of this common causative pathogen (EBV) with different carcinomas such as GC, NPC, HL and BL. Through network analysis we identified the top 10 genes linked with EBV associated carcinomas as CASP10, BRAF, NFKBIA, IFNA2, GSTP1, CSF3, GATA3, UBR5, AXIN2 and POLE. Further, the tyrosine-protein kinase (ABL1) gene was significantly over-represented in 3 out of 9 critical biological processes, viz. in regulatory pathways in cancer, the TP53 network and the Imatinib and chronic myeloid leukemia biological processes. Consequently, the EBV pathogen appears to target critical pathways involved in cellular growth arrest/apoptosis. We make our case for BCR-ABL1 tyrosine-kinase inhibitors (TKI) for further clinical investigations in the inhibition of BCR-mediated EBV activation in carcinomas for better prognostic and therapeutic outcomes.

Both fine-grained and coarse-grained spatial patterns of neural activity measured by functional MRI show preferential encoding of pain in the human brain

How pain emerges from human brain remains an unresolved question in pain neuroscience. Neuroimaging studies have suggested that all brain areas activated by painful stimuli were also activated by tactile stimuli, and vice versa. Nonetheless, pain-preferential spatial patterns of voxel-level activation in the brain have been observed when distinguishing painful and tactile brain activations using multivariate pattern analysis (MVPA). According to two hypotheses, the neural activity pattern preferentially encoding pain could exist at a global, coarse-grained, regional level, corresponding to the "pain connectome" hypothesis proposing that pain-preferential information may be encoded by the synchronized activity across multiple distant brain regions, and/or exist at a local, fine-grained, voxel level, corresponding to the "intermingled specialized/preferential neurons" hypothesis proposing that neurons responding specially or preferentially to pain could be present and intermingled with non-pain neurons within a voxel. Here, we systematically investigated the spatial scales of pain-distinguishing information in the human brain measured by fMRI using machine learning techniques, and found that pain-distinguishing information could be detected at both coarse-grained spatial scales across widely distributed brain regions and fine-grained spatial scales within many local areas. Importantly, the spatial distribution of pain-distinguishing information in the brain varies across individuals and such inter-individual variations may be related to a person's trait about pain perception, particularly the pain vigilance and awareness. These results provide new insights into the long-standing question of how pain is represented in the human brain and help the identification of characteristic neuroimaging measurements of pain.

Genetic variants associated with syncope implicate neural and autonomic processes

AIMS

Syncope is a common and clinically challenging condition. In this study, the genetics of syncope were investigated to seek knowledge about its pathophysiology and prognostic implications.

METHODS AND RESULTS

This genome-wide association meta-analysis included 56 071 syncope cases and 890 790 controls from deCODE genetics (Iceland), UK Biobank (United Kingdom), and Copenhagen Hospital Biobank Cardiovascular Study/Danish Blood Donor Study (Denmark), with a follow-up assessment of variants in 22 412 cases and 286 003 controls from Intermountain (Utah, USA) and FinnGen (Finland). The study yielded 18 independent syncope variants, 17 of which were novel. One of the variants, p.Ser140Thr in PTPRN2, affected syncope only when maternally inherited. Another variant associated with a vasovagal reaction during blood donation and five others with heart rate and/or blood pressure regulation, with variable directions of effects. None of the 18 associations could be attributed to cardiovascular or other disorders. Annotation with regard to regulatory elements indicated that the syncope variants were preferentially located in neural-specific regulatory regions. Mendelian randomization analysis supported a causal effect of coronary artery disease on syncope. A polygenic score (PGS) for syncope captured genetic correlation with cardiovascular disorders, diabetes, depression, and shortened lifespan. However, a score based solely on the 18 syncope variants performed similarly to the PGS in detecting syncope risk but did not associate with other disorders.

CONCLUSION

The results demonstrate that syncope has a distinct genetic architecture that implicates neural regulatory processes and a complex relationship with heart rate and blood pressure regulation. A shared genetic background with poor cardiovascular health was observed, supporting the importance of a thorough assessment of individuals presenting with syncope.

Incorporating knowledge of disease-defining hub genes and regulatory network into a machine learning-based model for predicting treatment response in lupus nephritis after the first renal flare

BACKGROUND

Identifying candidates responsive to treatment is important in lupus nephritis (LN) at the renal flare (RF) because an effective treatment can lower the risk of progression to end-stage kidney disease. However, machine learning (ML)-based models that address this issue are lacking.

METHODS

Transcriptomic profiles based on DNA microarray data were extracted from the GSE32591 and GSE112943 datasets. Comprehensive bioinformatics analyses were performed to identify disease-defining genes (DDGs). Peripheral blood samples (GSE81622, GSE99967, and GSE72326) were used to evaluate the effect of DDGs. Single-sample gene set enrichment analysis (ssGSEA) scores of the DDGs were calculated and correlated with specific immunology genes listed in the nCounter panel. GSE60681 and GSE69438 were used to examine the ability of the DDGs to discriminate LN from other renal diseases. K-means clustering was used to obtain the separate gene sets. The clustering results were extended to data derived using the nCounter technique. The least absolute shrinkage and selection operator (LASSO) algorithm was used to identify genes with high predictive value for treatment response after the first RF in each cluster. LASSO models with tenfold validation were built in GSE200306 and assessed by receiver operating characteristic (ROC) analysis with area under curve (AUC). The models were validated by using an independent dataset (GSE113342).

RESULTS

Forty-five hub genes specific to LN were identified. Eight optimal disease-defining clusters (DDCs) were identified in this study. Th1 and Th2 cell differentiation pathway was significantly enriched in DDC-6. LCK in DDC-6, whose expression positively correlated with various subsets of T cell infiltrations, was found to be differentially expressed between responders and non-responders and was ranked high in regulatory network analysis. Based on DDC-6, the prediction model had the best performance (AUC: 0.75; 95% confidence interval: 0.44-1 in the testing set) and high precision (0.83), recall (0.71), and F1 score (0.77) in the validation dataset.

CONCLUSIONS

Our study demonstrates that incorporating knowledge of biological phenotypes into the ML model is feasible for evaluating treatment response after the first RF in LN. This knowledge-based incorporation improves the model's transparency and performance. In addition, LCK may serve as a biomarker for T-cell infiltration and a therapeutic target in LN.

Artificial Neural Network Study on the Pyrolysis of Polypropylene with a Sensitivity Analysis

Among machine learning (ML) studies, artificial neural network (ANN) analysis is the most widely used technique in pyrolysis research. In this work, the pyrolysis of polypropylene (PP) polymers was established using a thermogravimetric analyzer (TGA) with five sets of heating rates (5-40 K min^-1). TGA data was used to exploit an ANN network by achieving a feed-forward backpropagation optimization technique in order to predict the weight-left percentage. Two important ANN model input variables were identified as the heating rate (K min^-1) and temperature (K). For the range of TGA values, a 2-10-10-1 network with two hidden layers (Logsig-Tansig) was concluded to be the best structure for predicting the weight-left percentage. The ANN demonstrated a good agreement between the experimental and calculated values, with a high correlation coefficient (R) of greater than 0.9999. The final network was then simulated with the new input data set for effective performance. In addition, a sensitivity analysis was performed to identify the uncertainties associated with the relationship between the output and input parameters. Temperature was found to be a more sensitive input parameter than the heating rate on the weight-left percentage calculation.

Sample Size Calculation in Genetic Association Studies: A Practical Approach

Genetic association studies, testing the relationship between genetic variants and disease status, are useful tools for identifying genes that grant susceptibility to complex disorders. In such studies, an inadequate sample size may provide unreliable results: a small sample is unable to accurately describe the population, whereas a large sample makes the study expensive and complex to run. However, in genetic association studies, the sample size calculation is often overlooked or inadequately assessed for the small number of parameters included. In light of this, herein we list and discuss the role of the statistical and genetic parameters to be considered in the sample size calculation, show examples reporting incorrect estimation and, by using a genetic software program, we provide a practical approach for the assessment of the adequate sample size in a hypothetical study aimed at analyzing a gene-disease association.

EpiMix: an integrative tool for epigenomic subtyping using DNA methylation

DNA methylation (DNAme) is a major epigenetic factor influencing gene expression with alterations leading to cancer, immunological, and cardiovascular diseases. Recent technological advances enable genome-wide quantification of DNAme in large human cohorts. So far, existing methods have not been evaluated to identify differential DNAme present in large and heterogeneous patient cohorts. We developed an end-to-end analytical framework named "EpiMix" for population-level analysis of DNAme and gene expression. Compared to existing methods, EpiMix showed higher sensitivity in detecting abnormal DNAme that was present in only small patient subsets. We extended the model-based analyses of EpiMix to cis-regulatory elements within protein-coding genes, distal enhancers, and genes encoding microRNAs and lncRNAs. Using cell-type specific data from two separate studies, we discovered novel epigenetic mechanisms underlying childhood food allergy and survival-associated, methylation-driven non-coding RNAs in non-small cell lung cancer.

A statistical genetic investigation of psychiatric resilience

Background: Although trauma exposure (TE) is a transdiagnostic risk factor for many psychiatric disorders, not everyone who experiences TE develops a psychiatric disorder. Resilience may explain this heterogeneity; thus, it is critical to understand the etiologic underpinnings of resilience.Objective: The present study sought to examine the genetic underpinnings of psychiatric resilience using genome-wide association studies (GWAS), genome-wide complex trait analysis (GCTA), and polygenic risk score (PRS) analyses.Method: Participants were 6,634 trauma exposed college students attending a diverse, public university in the Mid Atlantic. GWAS and GCTA analyses were conducted, and using GWAS summary statistics from large genetic consortia, PRS analyses examined the shared genetic risk between resilience and various phenotypes.Results: Results demonstrate that nine single-nucleotide polymorphisms (SNPs) met the suggestive of significance threshold, heritability estimates for resilience were non-significant, and that there is genetic overlap between resilience and AD, as well as resilience and PTSD.Conclusion: Mixed findings from the present study suggest additional research to elucidate the etiological underpinnings of resilience, ideally with larger samples less biased by variables such as heterogeneity (i.e. clinical vs. population based) and population stratification. Genetic investigations of resilience have the potential to elucidate the molecular bases of stress-related psychopathology, suggesting new avenues for prevention and intervention efforts.

Genetics and epigenetics of human aggression

There is substantial variation between humans in aggressive behavior, with its biological etiology and molecular genetic basis mostly unknown. This review chapter offers an overview of genomic and omics studies revealing the genetic contribution to aggression and first insights into associations with epigenetic and other omics (e.g., metabolomics) profiles. We allowed for a broad phenotype definition including studies on "aggression," "aggressive behavior," or "aggression-related traits," "antisocial behavior," "conduct disorder," and "oppositional defiant disorder." Heritability estimates based on family and twin studies in children and adults of this broadly defined phenotype of aggression are around 50%, with relatively small fluctuations around this estimate. Next, we review the genome-wide association studies (GWAS) which search for associations with alleles and also allow for gene-based tests and epigenome-wide association studies (EWAS) which seek to identify associations with differently methylated regions across the genome. Both GWAS and EWAS allow for construction of Polygenic and DNA methylation scores at an individual level. Currently, these predict a small percentage of variance in aggression. We expect that increases in sample size will lead to additional discoveries in GWAS and EWAS, and that multiomics approaches will lead to a more comprehensive understanding of the molecular underpinnings of aggression.

Association between PTPN1 polymorphisms and obesity-related phenotypes in European adolescents: influence of physical activity

BACKGROUND

To study the associations of Protein Tyrosine Phosphatase-N1 (PTPN1) polymorphisms with obesity-related phenotypes in European adolescents, and the influence of physical activity on these relationships.

METHODS

Five polymorphisms of PTPN1 were genotyped in 1057 European adolescents (12-18 years old). We measured several phenotypes related to obesity, such as adiposity markers, and biochemical and clinical parameters. Physical activity was objectively measured by accelerometry.

RESULTS

The T, A, T, T and G alleles of the rs6067472, rs10485614, rs2143511, rs6020608 and rs968701 polymorphisms, respectively, were associated with lower levels of obesity-related phenotypes (i.e., body mass index, body fat percentage, hip circumference, fat mass index, systolic blood pressure and leptin) in European adolescents. In addition, the TATTG haplotype was associated with lower body fat percentage and fat mass index compared to the AACCA haplotype. Finally, when physical activity levels were considered, alleles of the rs6067472, rs2143511, rs6020608 and rs968701 polymorphisms were only associated with lower adiposity in active adolescents.

CONCLUSIONS

PTPN1 polymorphisms were associated with adiposity in European adolescents. Specifically, alleles of these polymorphisms were associated with lower adiposity only in physically active adolescents. Therefore, meeting the recommendations of daily physical activity may reduce obesity risk by modulating the genetic predisposition to obesity.

IMPACT

Using gene-phenotype and gene*environment analyses, we detected associations between polymorphisms of the Protein Tyrosine Phosphatase-N1 (PTPN1) gene and obesity-related phenotypes, suggesting a mechanism that can be modulated by physical activity. This study shows that genetic variability of PTPN1 is associated with adiposity, while physical activity seems to modulate the genetic predisposition. This brings insights about the mechanisms by which physical activity positively influences obesity.

Opportunities and challenges for the use of common controls in sequencing studies

Genome-wide association studies using large-scale genome and exome sequencing data have become increasingly valuable in identifying associations between genetic variants and disease, transforming basic research and translational medicine. However, this progress has not been equally shared across all people and conditions, in part due to limited resources. Leveraging publicly available sequencing data as external common controls, rather than sequencing new controls for every study, can better allocate resources by augmenting control sample sizes or providing controls where none existed. However, common control studies must be carefully planned and executed as even small differences in sample ascertainment and processing can result in substantial bias. Here, we discuss challenges and opportunities for the robust use of common controls in high-throughput sequencing studies, including study design, quality control and statistical approaches. Thoughtful generation and use of large and valuable genetic sequencing data sets will enable investigation of a broader and more representative set of conditions, environments and genetic ancestries than otherwise possible.

Polygenic power calculator: Statistical power and polygenic prediction accuracy of genome-wide association studies of complex traits

Power calculation is a necessary step when planning genome-wide association studies (GWAS) to ensure meaningful findings. Statistical power of GWAS depends on the genetic architecture of phenotype, sample size, and study design. While several computer programs have been developed to perform power calculation for single SNP association testing, it might be more appropriate for GWAS power calculation to address the probability of detecting any number of associated SNPs. In this paper, we derive the statistical power distribution across causal SNPs under the assumption of a point-normal effect size distribution. We demonstrate how key outcome indices of GWAS are related to the genetic architecture (heritability and polygenicity) of the phenotype through the power distribution. We also provide a fast, flexible and interactive power calculation tool which generates predictions for key GWAS outcomes including the number of independent significant SNPs, the phenotypic variance explained by these SNPs, and the predictive accuracy of resulting polygenic scores. These results could also be used to explore the future behaviour of GWAS as sample sizes increase further. Moreover, we present results from simulation studies to validate our derivation and evaluate the agreement between our predictions and reported GWAS results.

Joint analysis of functionally related genes yields further candidates associated with Tetralogy of Fallot

Although several genes involved in the development of Tetralogy of Fallot have been identified, no genetic diagnosis is available for the majority of patients. Low statistical power may have prevented the identification of further causative genes in gene-by-gene survey analyses. Thus, bigger samples and/or novel analytic approaches may be necessary. We studied if a joint analysis of groups of functionally related genes might be a useful alternative approach. Our reanalysis of whole-exome sequencing data identified 12 groups of genes that exceedingly contribute to the burden of Tetralogy of Fallot. Further analysis of those groups showed that genes with high-impact variants tend to interact with each other. Thus, our results strongly suggest that additional candidate genes may be found by studying the protein interaction network of known causative genes. Moreover, our results show that the joint analysis of functionally related genes can be a useful complementary approach to classical single-gene analyses.

Genetic determinants of apixaban plasma levels and their relationship to bleeding and thromboembolic events

Apixaban is a direct oral anticoagulant, a factor Xa inhibitor, used for the prevention of ischemic stroke in patients with atrial fibrillation. Despite using recommended dosing a few patients might still experience bleeding or lack of efficacy that might be related to inappropriate drug exposure. We conducted a genome-wide association study using data from 1,325 participants in the pivotal phase three trial of apixaban with the aim to identify genetic factors affecting the pharmacokinetics of apixaban. A candidate gene analysis was also performed for pre-specified variants in ABCB1, ABCG2, CYP3A4, CYP3A5, and SULT1A1, with a subsequent analysis of all available polymorphisms within the candidate genes. Significant findings were further evaluated to assess a potential association with clinical outcome such as bleeding or thromboembolic events. No variant was consistently associated with an altered apixaban exposure on a genome-wide level. The candidate gene analyses showed a statistically significant association with a well-known variant in the drug transporter gene ABCG2 (c.421G &gt; T, rs2231142). Patients carrying this variant had a higher exposure to apixaban [area under the curve (AUC), beta = 151 (95% CI 59–243), p = 0.001]. On average, heterozygotes displayed a 5% increase of AUC and homozygotes a 17% increase of AUC, compared with homozygotes for the wild-type allele. Bleeding or thromboembolic events were not significantly associated with ABCG2 rs2231142. This large genome-wide study demonstrates that genetic variation in the drug transporter gene ABCG2 is associated with the pharmacokinetics of apixaban. However, the influence of this finding on drug exposure was small, and further studies are needed to better understand whether it is of relevance for ischemic and bleeding events.

Genome-wide association study of liver enzyme elevation in an extended cohort of rheumatoid arthritis patients starting low-dose methotrexate

Aim: A follow-up genome-wide association study (GWAS) in an extended cohort of rheumatoid arthritis (RA) patients starting low-dose methotrexate (MTX) treatment was performed to identify further genetic variants associated with alanine aminotransferase (ALT) elevation. Patients & methods: A GWAS was performed on 346 RA patients. Two outcomes within the first 6 months of MTX treatment were assessed: ALT >1.5-times the upper level of normal (ULN) and maximum level of ALT. Results: SPATA9 (rs72783407) was significantly associated with maximum level of ALT (p = 2.58 × 10^-8) and PLCG2 (rs60427389) was tentatively associated with ALT >1.5 × ULN. Conclusion: Associations with SNPs in genes related to male fertility (SPATA9) and inflammatory processes (PLCG2) were identified.

Genetic associations with lymphomas in Polish patients: A pooled-DNA genome-wide association analysis

Several single nucleotide polymorphisms (SNPs) associated with susceptibility to Hodgkin lymphoma (HL) and diffuse large B-cell lymphoma (DLBCL) have been identified. The aim of this study was to identify susceptibility loci for HL and DLBCL in Polish patients. Altogether, DLBCL (n = 218 and HL patients (n = 224) and healthy individuals (n = 1181) were recruited. Lymphoma diagnosis was based on standard criteria. Genome-wide association study (GWAS) was performed using pooled-DNA samples on llumina Infinium Omni2.5 Exome-8 v1.3, and selected loci were replicated by TaqMan SNP genotyping of individuals. GWAS detected thirteen and seven SNPs associated with DLBCL and HL, respectively. In the replication study, six and seven SNPs reached significance after correction for multiple testing in the DLBCL and HL cohorts, respectively. One and four SNPs associated with DLBCL and HL, respectively, were localized within, and two SNPs-near the major histocompatibility complex (MHC) region. In conclusion, the majority of loci associated with HL and DLBCL aetiology in previous studies have potential roles in immune function. Our pooled-DNA GWAS enabled the identification of several susceptibility loci for DLBCL and HL in the Polish population; some of them were mapped within or adjacent to the MHC, and other associated SNPs were located outside the MHC.

Kawasaki Disease Arab Initiative [Kawarabi]: Establishment and Results of a Multicenter Survey

Studies on Kawasaki disease (KD) in Arab countries are scarce, often providing incomplete data. This along with the benefits of multicenter research collaboratives led to the creation of the KD Arab Initiative [Kawarabi] consortium. An anonymous survey was completed among potential collaborative Arab medical institutions to assess burden of KD in those countries and resources available to physicians. An online 32-item survey was distributed to participating institutions after conducting face validity. One survey per institution was collected. Nineteen physicians from 12 countries completed the survey representing 19 out of 20 institutions (response rate of 95%). Fifteen (79%) institutions referred to the 2017 American Heart Association guidelines when managing a patient with KD. Intravenous immunoglobulin (IVIG) is not readily available at 2 institutions (11%) yet available in the country. In one center (5%), IVIG is imported on-demand. The knowledge and awareness among countries' general population was graded (0 to 10) at median/interquartiles (IQR) 3 (2-5) and at median/IQR 7 (6-8) in the medical community outside their institution. Practice variations in KD management and treatment across Arab countries require solid proactive collaboration. The low awareness and knowledge estimates about KD among the general population contrasted with a high level among the medical community. The Kawarabi collaborative will offer a platform to assess disease burden of KD, among Arab population, decrease practice variation and foster population-based knowledge.

The Impact of Donor and Recipient Genetic Variation on Outcomes After Solid Organ Transplantation: A Scoping Review and Future Perspectives

At the outset of solid organ transplantation, genetic variation between donors and recipients was recognized as a major player in mechanisms such as allograft tolerance and rejection. Genome-wide association studies have been very successful in identifying novel variant-trait associations, but have been difficult to perform in the field of solid organ transplantation due to complex covariates, era effects, and poor statistical power for detecting donor-recipient interactions. To overcome a lack of statistical power, consortia such as the International Genetics and Translational Research in Transplantation Network have been established. Studies have focused on the consequences of genetic dissimilarities between donors and recipients and have reported associations between polymorphisms in candidate genes or their regulatory regions with transplantation outcomes. However, knowledge on the exact influence of genetic variation is limited due to a lack of comprehensive characterization and harmonization of recipients' or donors' phenotypes and validation using an experimental approach. Causal research in genetics has evolved from agnostic discovery in genome-wide association studies to functional annotation and clarification of underlying molecular mechanisms in translational studies. In this overview, we summarize how the recent advances and progresses in the field of genetics and genomics have improved the understanding of outcomes after solid organ transplantation.

Single-Nucleotide Polymorphisms in the Beta-Tubulin Gene and Its Relationship with Treatment Response to Albendazole in Human Soil-Transmitted Helminths in Southern Mozambique

Soil-transmitted helminth (STH) cornerstone control strategy is mass drug administration (MDA) with benzimidazoles. However, MDA might contribute to selection pressure for anthelmintic resistance, as occurred in livestock. The aim of this study is to evaluate the treatment response to albendazole and the relationship with the presence of putative benzimidazole resistance single-nucleotide polymorphisms (SNPs) in the β-tubulin gene of STH in Southern Mozambique. After screening 819 participants, we conducted a cohort study with 184 participants infected with STH in Manhiça district, Southern Mozambique. A pretreatment and a posttreatment stool samples were collected and the STH infection was identified by duplicate Kato-Katz and quantitative polymerase chain reaction (qPCR). Cure rate and egg reduction rates were calculated. Putative benzimidazole resistance SNPs (F167Y, F200T, and E198A) in Trichuris trichiura and Necator americanus were assessed by pyrosequencing. Cure rates by duplicate Kato-Katz and by qPCR were 95.8% and 93.6% for Ascaris lumbricoides, 28% and 7.8% for T. trichiura, and 88.9% and 56.7% for N. americanus. Egg reduction rate by duplicate Kato-Katz was 85.4% for A. lumbricoides, 34.9% for T. trichiura, and 40.5% for N. americanus. Putative benzimidazole resistance SNPs in the β-tubulin gene were detected in T. trichiura (23%) and N. americanus (21%) infected participants at pretreatment. No statistical difference was observed between pretreatment and posttreatment frequencies for none of the SNPs. Although treatment response to albendazole was low, particularly in T. trichiura, the putative benzimidazole resistance SNPs were not higher after treatment in the population studied. New insights are needed for a better understanding and monitoring of human anthelmintic resistance.