Impact of a nationwide school policy on body mass index in Danish school children: An interrupted time series analysis.

Objectives: A new school policy mandating 45 min physical activity daily during school was introduced in Denmark in 2014. We aimed to evaluate the effect of this policy on BMI in school-aged children. It was hypothesized that the school policy would decrease BMI, especially in the obese fraction of the population (90th percentile BMI). Study design: This register-based study was conducted as a natural experiment. Methods: Analyses were based on data from The National Child Health Register that contains nationwide data on height and weight from mandatory preventive health examinations completed by school nurses or medical doctors during pre-preparatory classes (0th-3rd grade) and lower secondary education (7th-9th grade). A total of 401,517 children were included in the analyses with annual repeated cross-sectional data covering the period from 2012 to 2018. The effect of the school policy was evaluated using an interrupted time series approach comparing pre- and post-policy slopes in BMI, stratified by sex and age-group. Results: In boys, no significant differences were observed in mean BMI slopes from pre-to post-policy in either age-group. In girls, post-policy slopes were significantly higher compared to pre-policy in both age-groups (0th-3rd grade: ß:0·034 kg/m2, 95%-CI: (0·024; 0·043), p-value: <0·001; 7th-9th grade: ß:0·066 kg/m2, 95%-CI: (0·028; 0·103), p-value: 0·001). No significant differences in slopes were observed in BMI at the 90th percentile from pre-to post-policy for both sexes and across both age-groups. Adjustment for leisure-time physical activity as a potential time-varying confounder did not alter the findings. Conclusions: In conclusion, we did not detect a significant decrease in BMI levels among school-aged children following the introduction of a nationwide school policy specifying daily physical activity in school. If anything, a small positive change in BMI was observed in girls. More research is needed to understand whether structural changes similar to this requirement are able to prevent overweight and obesity in children and adolescents.

Incidence and assessment of delirium following open cardiac surgery; A systematic review and meta-analysis.

AIM: This systematic review and meta-analysis sought i) to provide an overview of the incidence of delirium following open cardiac surgery and ii) to investigate how incidences of delirium are associated with different assessment tools. METHODS AND RESULTS: A systematic search of studies investigating delirium following open cardiac surgery were conducted in Medline (Ovid), EMBASE, PsycINFO, CiNAHL and the Cochrane Database. Only studies with patients diagnosed or screened with a validated tool were included. Studies published from 2005 to 2021 were included in the meta-analysis.Of 7,126 individual studies retrieved, 106 met the inclusion criteria for the meta-analysis, hereof 31% of high quality. The weighted pooled incidence of delirium following open cardiac surgery across all studies was 23% (95% CI 20-26%), however we found a considerable heterogeneity (I2 = 99%), which could not be explained by subgroups or further sensitivity analyses. The most commonly applied screening tool for delirium is CAM/CAM-ICU. The lowest estimates of delirium were found by applying the Delirium Observation Scale (incidence 14%, 95% CI 8-20%), and the highest estimates in studies using "other" screening tools (Organic Brain Symptom Scale, Delirium Symptom Interview) pooled incidence of 43%, (95% CI 19 - 66%), however, only two studies applied these. CONCLUSION: Delirium following open cardiac surgery remains a complication with a high incidence of overall 23%, when applying a validated tool for screening or diagnosis. Nevertheless, this systematic review and meta-analyses highlight the significant inconsistency in current evidence regarding assessment tools and regimens. REGISTRATION: Prospero CRD42020215519.

A robust and fast two-sample test of equal correlations with an application to differential co-expression.

A robust and fast two-sample test for equal Pearson correlation coefficients (PCCs) is important in solving many biological problems, including, for example, analysis of differential co-expression. However, few existing methods for this test can achieve robustness against deviation from normal distributions, accuracy under small sample sizes, and computational efficiency simultaneously. Here, we propose a new method for testing differential correlation using a saddlepoint approximation of the residual bootstrap (DICOSAR). To achieve robustness, accuracy, and efficiency, DICOSAR combines the ideas underlying the pooled residual bootstrap, the signed root of a likelihood ratio statistic, and a multivariate saddlepoint approximation. Through a comprehensive simulation study and a real data analysis of gene co-expression, we demonstrate that DICOSAR is accurate and robust in controlling the type I error rate for detecting differential correlation and provides a faster alternative to the bootstrap and permutation methods. We further show that DICOSAR can also be used for testing differential correlation matrices. These results suggest that DICOSAR provides an analytical approach to facilitate rapid testing for the equality of PCCs in large-scale analysis.

Identifying genetic variants regulating MGMT gene expression - A study in monozygotic Danish twins.

Identifying genetic factors affecting the regulation of the O-6-Methylguanine-DNA Methyltransferase (MGMT) gene and estimating the genetic contribution of the MGMT gene through within-pair correlation in monozygotic twin pairs is of particular importance in various types of cancer such as glioblastoma. We used gene expression data in whole blood from 448 monozygotic twins from the Middle Age Danish Twins (MADT) study to investigate genetic regulation of the MGMT gene by performing a genome-wide association study (GWAS) of the variation in MGMT expression. Additionally, we estimated within-pair dependence measures of the expression values looking for the genetic influence of significant identified genes. We identified 243 single nucleotide polymorphisms (SNPs) significantly (p < 5e-8) associated with expression of MGMT, all located on chromosome 10 near the MGMT gene. Of the 243 SNPs, 7 are novel cis-eQTLs. By further looking into the suggestively significant SNPs (increasing cutoff to p = 1e-6), we identified 11 suggestive trans-eQTLs located on chromosome 17. These variants were in or proximal to a total of seven genes, which may regulate MGMT expression. The within-pair correlation of the expression of MGMT, TRIM37, and SEPT4 provided the upper bound genetic influence of these genes. Overall, identifying cis- or trans-acting genetic variations regulating the MGMT gene can pave the way for a better understanding of the MGMT gene function and ultimately in understanding the patient's sensitivity to therapeutic alkylating agents.

Circulating microRNAs related to lipid metabolism and solid tissue maintenance and morphology associate with mortality in elderly twins.

The lifespan of humans varies greatly between individuals. Here, we aimed to explore what biological roles miRNAs may have on old age mortality-variation. Circulating miRNAs were measured in plasma from 43 monozygotic twin pairs (73-95 years of age) and mortality analyses were applied using Cox regression survival analyses and linear regression analyses of lifespan. In general, nominally significant miRNAs were mainly upregulated with shorter lifespan, both in Cox analysis (72 % upregulated) and in linear regression analysis (81 % upregulated). A total of 29 miRNAs were associated to mortality at a nominal significance level (p < 0.05) in the survival analysis, but no miRNAs passed the FDR adjusted level of significance. Seven of the 29 miRNAs; hsa-miR-140-3p, hsa-miR-16-5p, hsa-miR-487b-3p, hsa-miR-19a-3p, hsa-let-7d-5p, hsa-miR-320a, hsa-miR-375, were nominally significant across two linear twin-paired analyses and the cox analysis. Pathway analyses of the 29 nominally significant miRNAs from the individual level analyses resulted in two nominally significant associated Reactome pathways (unadjusted p < 0.05); 'Negative regulation of FGFR signaling' and 'Neurotransmitter receptor binding and downstream transmission in the postsynaptic cell', and two significantly associated KEGG pathways; 'Linoleic acid metabolism' and 'Toxoplasmosis'. Additional pathway analyses and results of previous studies support that miRNAs linked to mortality at age 70 years or older play a role in lipid metabolism, tissues maintenance and morphology.

The heritability of BMI varies across the range of BMI-a heritability curve analysis in a twin cohort.

BACKGROUND: The heritability of traits such as body mass index (BMI), a measure of obesity, is generally estimated using family and twin studies, and increasingly by molecular genetic approaches. These studies generally assume that genetic effects are uniform across all trait values, yet there is emerging evidence that this may not always be the case. METHOD/SUBJECTS: This paper analyzes twin data using a recently developed measure of heritability called the heritability curve. Under the assumption that trait values in twin pairs are governed by a flexible Gaussian mixture distribution, heritability curves may vary across trait values. The data consist of repeated measures of BMI on 1506 monozygotic (MZ) and 2843 like-sexed dizygotic (DZ) adult twin pairs, gathered from multiple surveys in older Finnish Twin Cohorts. RESULTS: The heritability curve and BMI value-specific MZ and DZ pairwise correlations were estimated, and these varied across the range of BMI. MZ correlations were highest at BMI values from 21 to 24, with a stronger decrease for women than for men at higher values. Models with additive and dominance effects fit best at low and high BMI values, while models with additive genetic and common environmental effects fit best in the normal range of BMI. CONCLUSIONS: We demonstrate that twin and molecular genetic studies need to consider how genetic effects vary across trait values. Such variation may reconcile findings of traits with high heritability and major differences in mean values between countries or over time.

Genome-wide association meta-analysis identifies 48 risk variants and highlights the role of the stria vascularis in hearing loss.

Hearing loss is one of the top contributors to years lived with disability and is a risk factor for dementia. Molecular evidence on the cellular origins of hearing loss in humans is growing. Here, we performed a genome-wide association meta-analysis of clinically diagnosed and self-reported hearing impairment on 723,266 individuals and identified 48 significant loci, 10 of which are novel. A large proportion of associations comprised missense variants, half of which lie within known familial hearing loss loci. We used single-cell RNA-sequencing data from mouse cochlea and brain and mapped common-variant genomic results to spindle, root, and basal cells from the stria vascularis, a structure in the cochlea necessary for normal hearing. Our findings indicate the importance of the stria vascularis in the mechanism of hearing impairment, providing future paths for developing targets for therapeutic intervention in hearing loss.

Causal relation between heart irradiation and survival of lung cancer patients after radiotherapy.

INTRODUCTION: In a recent study, setup uncertainties in the direction of the heart were shown to impact the overall survival of non-small cell lung cancer (NSCLC) patients after radiotherapy, indicating the causal effect between heart irradiation and survival. The current study aims to externally evaluate this observation within a patient cohort treated using daily IGRT. METHOD: NSCLC patients with locally-advanced disease and daily CBCT were included. For all treatment fractions, the distance between the isocenter and the heart was evaluated based on the clinical setup registrations. The variation in heart position between planning and treatment (DeltaDistance) was estimated from these registrations. The possible impact of DeltaDistance on survival was analysed by a multivariable Cox model of overall survival, allowing for a time-dependent impact of DeltaDistance to allow for toxicity latency. RESULTS: Daily CBCT information was available for 489 patients at Odense University Hospital. The primary Cox model contained GTV volume, patient age, performance status, and DeltaDistance. DeltaDistance significantly impacted overall survival approximately 50 months after radiotherapy. Subanalyses indicated that the observed effect is mainly present among the patients with the least clinical risk factors. CONCLUSION: Our results confirm the impact of setup variations in the direction of the heart on the survival of NSCLC patients, even within a cohort using daily CBCT setup guidance. This result indicates a causal effect between heart irradiation and survival. It will be challenging to reduce the setup uncertainty even further; thus, increased focus on dose constraints on the heart seems warranted.

Multivariate Generalized Linear Models for Twin and Family Data.

Multivariate twin and family studies are one of the most important tools to assess diseases inheritance as well as to study their genetic and environment interrelationship. The multivariate analysis of twin and family data is in general based on structural equation modelling or linear mixed models that essentially decomposes sources of covariation as originally suggested by  Fisher. In this paper, we propose a flexible and unified statistical modelling framework for analysing multivariate Gaussian and non-Gaussian twin and family data. The non-normality is taken into account by actually modelling the mean and variance relationship, while the covariance structure is modelled by means of a linear covariance model including the option to model the dispersion components as functions of known covariates in a regression model fashion. The marginal specification of our models allows us to extend classic models and biometric indices such as the bivariate heritability, genetic, environmental and phenotypic correlations to non-Gaussian data. We illustrate the proposed models through simulation studies and six data analyses and provide computational implementation in R through the package mglm4twin.

Familial risk and heritability of ischemic heart disease and stroke in Danish twins.

Aim: Our aim was to explore whether familial factors influence the risk of ischemic heart disease, stroke, and their co-occurrence. Methods: In total, 23,498 monozygotic and 39,540 same-sex dizygotic twins from the Danish Twin Registry were followed from 1977 to 2011 in the Danish National Patient Registry for ischemic heart disease and stroke. Time-to-event analyses accounting for censoring and competing risk of death were used to estimate familial risk (casewise concordance relative to the cumulative incidence) and heritability of ischemic heart disease, stroke, and the co-occurrence by age. Results: During follow-up, we observed 5561 and 4186 twin individuals with ischemic heart disease and stroke respectively, with 936 twin pairs concordant for ischemic heart disease and stroke. Familial risks were significant for both, with higher cumulative risks in monozygotic than in dizygotic twins. Estimates for heritability were significant for ischemic heart disease as well as for stroke diagnosed after the age of 80. The casewise concordance of ischemic heart disease in twins whose co-twin was diagnosed with stroke did not differ for monozygotic and dizygotic twins; however, from age 55 it was 10% higher than the cumulative risk in the overall twin cohort and was 25% higher at age 90. A similar pattern was seen for stroke following the co-twin's ischemic heart disease. Conclusions: As in previous studies, we found a higher heritability of ischemic heart disease than of stroke. There was a significant familial risk but no heritability for the co-occurrence of ischemic heart disease and stroke. The co-occurrence is therefore likely due to other shared familial than genetic factors, highlighting that preventive initiatives should target families rather than individuals.

Familial Risk and Heritability of Hematologic Malignancies in the Nordic Twin Study of Cancer.

We aimed to explore the genetic and environmental contributions to variation in the risk of hematologic malignancies and characterize familial dependence within and across hematologic malignancies. The study base included 316,397 individual twins from the Nordic Twin Study of Cancer with a median of 41 years of follow-up: 88,618 (28%) of the twins were monozygotic, and 3459 hematologic malignancies were reported. We estimated the cumulative incidence by age, familial risk, and genetic and environmental variance components of hematologic malignancies accounting for competing risk of death. The lifetime risk of any hematologic malignancy was 2.5% (95% CI 2.4-2.6%), as in the background population. This risk was elevated to 4.5% (95% CI 3.1-6.5%) conditional on hematologic malignancy in a dizygotic co-twin and was even greater at 7.6% (95% CI 4.8-11.8%) if a monozygotic co-twin had a hematologic malignancy. Heritability of the liability to develop any hematologic malignancy was 24% (95% CI 14-33%). This estimate decreased across age, from approximately 55% at age 40 to about 20-25% after age 55, when it seems to stabilize. In this largest ever studied twin cohort with the longest follow-up, we found evidence for familial risk of hematologic malignancies. The discovery of decreasing familial predisposition with increasing age underscores the importance of cancer surveillance in families with hematological malignancies.

Novel DNA methylation marker discovery by assumption-free genome-wide association analysis of cognitive function in twins.

Privileged by rapid increase in available epigenomic data, epigenome-wide association studies (EWAS) are to make a profound contribution to understand the molecular mechanism of DNA methylation in cognitive aging. Current statistical methods used in EWAS are dominated by models based on multiple assumptions, for example, linear relationship between molecular profiles and phenotype, normal distribution for the methylation data and phenotype. In this study, we applied an assumption-free method, the generalized correlation coefficient (GCC), and compare it to linear models, namely the linear mixed model and kinship model. We use DNA methylation associated with a cognitive score in 400 and 206 twins as discovery and replication samples respectively. DNA methylation associated with cognitive function using GCC, linear mixed model, and kinship model, identified 65 CpGs (p < 1e-04) from discovery sample displaying both nonlinear and linear correlations. Replication analysis successfully replicated 9 of these top CpGs. When combining results of GCC and linear models to cover diverse patterns of relationships, we identified genes like KLHDC4, PAPSS2, and MRPS18B as well as pathways including focal adhesion, axon guidance, and some neurological signaling. Genomic region-based analysis found 15 methylated regions harboring 11 genes, with three verified in gene expression analysis, also the 11 genes were related to top functional clusters including neurohypophyseal hormone and maternal aggressive behaviors. The GCC approach detects valuable methylation sites missed by traditional linear models. A combination of methylation markers from GCC and linear models enriched biological pathways sensible in neurological function that could implicate cognitive performance and cognitive aging.

Generalized correlation coefficient for genome-wide association analysis of cognitive ability in twins.

Despite a strong genetic background in cognitive function only a limited number of single nucleotide polymorphisms (SNPs) have been found in genome-wide association studies (GWASs). We hypothesize that this is partially due to mis-specified modeling concerning phenotype distribution as well as the relationship between SNP dosage and the level of the phenotype. To overcome these issues, we introduced an assumption-free method based on generalized correlation coefficient (GCC) in a GWAS of cognitive function in Danish and Chinese twins to compare its performance with traditional linear models. The GCC-based GWAS identified two significant SNPs in Danish samples (rs71419535, p = 1.47e-08; rs905838, p = 1.69e-08) and two significant SNPs in Chinese samples (rs2292999, p = 9.27e-10; rs17019635, p = 2.50e-09). In contrast, linear models failed to detect any genome-wide significant SNPs. The number of top significant genes overlapping between the two samples in the GCC-based GWAS was higher than when applying linear models. The GCC model identified significant genetic variants missed by conventional linear models, with more replicated genes and biological pathways related to cognitive function. Moreover, the GCC-based GWAS was robust in handling correlated samples like twin pairs. GCC is a useful statistical method for GWAS that complements traditional linear models for capturing genetic effects beyond the additive assumption.

A Major Limitation of the Direction of Causation Model: Non-Shared Environmental Confounding.

Determining (1) the direction of causation and (2) the size of causal effects between two constructs is a central challenge of the scientific study of humans. In the early 1990s, researchers in behavioral genetics invented what was termed the direction of causation (DoC) model to address exactly these two concerns. The model claims that for any two traits whose mode of inheritance is sufficiently different, the direction of causation can be ascertained using a sufficiently large genetically informative sample. Using a series of simulation studies, we demonstrate a major challenge to the DoC model, namely that it is extremely sensitive to even tiny amounts of non-shared confounding. Even under ideal conditions for the DoC model (a large sample, N = 10,000), a large causal relationship (e.g., a causal correlation of .50) with very different modes of inheritance between the two traits (e.g., a pure AE model for one trait and a pure CE model for another trait) and a modest degree (correlation of .10) of non-shared confounding between the two traits results in the choice of the wrong causal models and estimating the wrong causal effects.

Evaluation of the proliferation marker Ki-67 in gliomas: Interobserver variability and digital quantification.

BACKGROUND: The Ki-67 Labelling Index (LI) is used as an ancillary tool in glioma diagnostics. Interobserver variability has been reported and no precise guidelines are available. Nor is it known whether novel digital approaches would be an advantage. Our aim was to evaluate the inter- and intraobserver variability of the Ki-67 LI between two pathologists and between pathologists and digital quantification both in whole tumour slides and in hot spots using narrow but diagnostically relevant intervals. METHODS: In samples of 235 low and high grade gliomas, two pathologists (A and B) estimated the Ki-67 LI (5-10% intervals) for whole tumour slides and for hot spots. In 20 of the cases intraobserver variability was evaluated. For digital quantification (C) slides were scanned with subsequent systematic random sampling of viable tumour areas. A software classifier trained to identify positive and negative nuclei calculated the Ki-67 LI. The interobserver agreements were evaluated using kappa (κ) statistics. RESULTS: The observed proportions of agreement and κ values for Ki-67 LI for whole tumour slides were: A/B: 46% (κ = 0.32); A/C: 37% (κ = 0.26); B/C: 37% (κ = 0.26). For hot spots equivalent values were: A/B: 14% (κ = 0.04); A/C: 18% (κ = 0.09); B/C: 31% (κ = 0.21). CONCLUSIONS: Interobserver variability was pronounced between pathologists and for pathologists versus digital quantification when attempting to estimate a precise value of the Ki-67 LI. Ki-67 LI should therefore be used with caution and should not be over interpreted in the grading of gliomas. Digital quantification of Ki-67 LI in gliomas was feasible, but intra- and interlaboratory robustness need to be determined.

Does the sex of one's co-twin affect height and BMI in adulthood? A study of dizygotic adult twins from 31 cohorts.

BACKGROUND: The comparison of traits in twins from opposite-sex (OS) and same-sex (SS) dizygotic twin pairs is considered a proxy measure of prenatal hormone exposure. To examine possible prenatal hormonal influences on anthropometric traits, we compared mean height, body mass index (BMI), and the prevalence of being overweight or obese between men and women from OS and SS dizygotic twin pairs. METHODS: The data were derived from the COllaborative project of Development of Anthropometrical measures in Twins (CODATwins) database, and included 68,494 SS and 53,808 OS dizygotic twin individuals above the age of 20 years from 31 twin cohorts representing 19 countries. Zygosity was determined by questionnaires or DNA genotyping depending on the study. Multiple regression and logistic regression models adjusted for cohort, age, and birth year with the twin type as a predictor were carried out to compare height and BMI in twins from OS pairs with those from SS pairs and to calculate the adjusted odds ratios and 95% confidence intervals for being overweight or obese. RESULTS: OS females were, on average, 0.31 cm (95% confidence interval (CI) 0.20, 0.41) taller than SS females. OS males were also, on average, taller than SS males, but this difference was only 0.14 cm (95% CI 0.02, 0.27). Mean BMI and the prevalence of overweight or obesity did not differ between males and females from SS and OS twin pairs. The statistically significant differences between OS and SS twins for height were small and appeared to reflect our large sample size rather than meaningful differences of public health relevance. CONCLUSIONS: We found no evidence to support the hypothesis that prenatal hormonal exposure or postnatal socialization (i.e., having grown up with a twin of the opposite sex) has a major impact on height and BMI in adulthood.

Telomeres and the natural lifespan limit in humans.

An ongoing debate in demography has focused on whether the human lifespan has a maximal natural limit. Taking a mechanistic perspective, and knowing that short telomeres are associated with diminished longevity, we examined whether telomere length dynamics during adult life could set a maximal natural lifespan limit. We define leukocyte telomere length of 5 kb as the 'telomeric brink', which denotes a high risk of imminent death. We show that a subset of adults may reach the telomeric brink within the current life expectancy and more so for a 100-year life expectancy. Thus, secular trends in life expectancy should confront a biological limit due to crossing the telomeric brink.

Genetic susceptibility to bilateral tinnitus in a Swedish twin cohort.

PURPOSE: Genetic contributions to tinnitus have been difficult to determine due to the heterogeneity of the condition and its broad etiology. Here, we evaluated the genetic and nongenetic influences on self-reported tinnitus from the Swedish Twin Registry (STR). METHODS: Cross-sectional data from the STR was obtained. Casewise concordance rates (the risk of one twin being affected given that his/her twin partner has tinnitus) were compared for monozygotic (MZ) and dizygotic (DZ) twin pairs (N = 10,464 concordant and discordant twin pairs) and heritability coefficients (the proportion of the total variance attributable to genetic factors) were calculated using biometrical model fitting procedures. RESULTS: Stratification of tinnitus cases into subtypes according to laterality (unilateral versus bilateral) revealed that heritability of bilateral tinnitus was 0.56; however, it was 0.27 for unilateral tinnitus. Heritability was greater in men (0.68) than in women (0.41). However, when female pairs younger than 40 years of age were selected, heritability of 0.62 was achieved with negligible effects of shared environment. CONCLUSION: Unlike unilateral tinnitus, bilateral tinnitus is influenced by genetic factors and might constitute a genetic subtype. Overall, our study provides the initial evidence for a tinnitus phenotype with a genetic influence.Genet Med advance online publication 23 March 2017.

Early-life mortality risks in opposite-sex and same-sex twins: a Danish cohort study of the twin testosterone transfer hypothesis.

PURPOSE: To investigate the twin testosterone transfer (TTT) hypothesis by comparing early-life mortality risks of opposite-sex (OS) and same-sex (SS) twins during the first 15 years of life. METHODS: We performed a population-based cohort study to compare mortality in OS and SS twins. We included 68,629 live-born Danish twins from 1973 to 2009 identified through the Danish Twin Registry and performed piecewise stratified Cox regression and log-binomial regression. RESULTS: Among 1933 deaths, we found significantly higher mortality for twin boys than for twin girls. For both sexes, OS twins had lower mortality than SS twins; the difference persisted for the first year of life for boys and for the first week of life for girls. CONCLUSIONS: Although the mortality risk for OS boys was in the expected direction according to the TTT hypothesis, the results for OS girls pointed in the opposite direction, providing no clear evidence for the TTT hypothesis.