Search

Venous thromboembolism (VTE) is a common disease with high heritability. However, only a small portion of the genetic variance of VTE can be explained by known genetic risk factors. Neutrophil extracellular traps (NETs) have been associated with prothrombotic activity. Therefore, the genetic basis of NETs could reveal novel risk factors for VTE. A recent genome-wide association study of plasma cell-free DNA (cfDNA) levels in the Genetic Analysis of Idiopathic Thrombophilia 2 (GAIT-2) Project showed a significant associated locus near ORM1. We aimed to further explore this candidate region by next-generation sequencing, copy number variation (CNV) quantification, and expression analysis using an extreme phenotype sampling design involving 80 individuals from the GAIT-2 Project. The RETROVE study with 400 VTE cases and 400 controls was used to replicate the results. A total of 105 genetic variants and a multiallelic CNV (mCNV) spanning ORM1 were identified in GAIT-2. Of these, 17 independent common variants, a region of 22 rare variants, and the mCNV were significantly associated with cfDNA levels. In addition, eight of these common variants and the mCNV influenced ORM1 expression. The association of the mCNV and cfDNA levels was replicated in RETROVE (p-value = 1.19 × 10−6). Additional associations between the mCNV and thrombin generation parameters were identified. Our results reveal that increased mCNV dosages in ORM1 decreased gene expression and upregulated cfDNA levels. Therefore, the mCNV in ORM1 appears to be a novel marker for cfDNA levels, which could contribute to VTE risk.

Schizophrenia is a polygenic psychiatric disorder with limited understanding about the mechanistic changes in gene expression regulation. To elucidate on this, we integrate interindividual variability of regulatory activity (ChIP-sequencing for H3K27ac histone mark) with gene expression and genotype data captured from the prefrontal cortex of 272 cases and controls. By measuring interindividual correlation among proximal chromatin peaks, we show that regulatory element activity is structured into 10,936 and 10,376 cis-regulatory domains in cases and controls, respectively. The schizophrenia-specific cis-regulatory domains are enriched for fetal-specific (p = 0.0014, OR = 1.52) and depleted of adult-specific regulatory activity (p = 3.04 × 10−50, OR = 0.57) and are enriched for SCZ heritability (p = 0.001). By studying the interplay among genetic variants, gene expression, and cis-regulatory domains, we ascertain that changes in coordinated regulatory activity tag alterations in gene expression levels (p = 3.43 × 10−5, OR = 1.65), unveil case-specific QTL effects, and identify regulatory machinery changes for genes affecting synaptic function and dendritic spine morphology in schizophrenia. Altogether, we show that accounting for coordinated regulatory activity provides a novel mechanistic approach to reduce the search space for unveiling genetically perturbed regulation of gene expression in schizophrenia.

Background: Human plasma contains a wide variety of circulating proteins. These proteins can be important clinical biomarkers in disease and also possible drug targets. Large scale genomics studies of circulating proteins can identify genetic variants that lead to relative protein abundance. Methods: We conducted a meta-analysis on genome-wide association studies of autosomal chromosomes in 22,997 individuals of primarily European ancestry across 12 cohorts to identify protein quantitative trait loci (pQTL) for 92 cardiometabolic associated plasma proteins. Results: We identified 503 (337 cis and 166 trans) conditionally independent pQTLs, including several novel variants not reported in the literature. We conducted a sex-stratified analysis and found that 118 (23.5%) of pQTLs demonstrated heterogeneity between sexes. The direction of effect was preserved but there were differences in effect size and significance. Additionally, we annotate trans-pQTLs with nearest genes and report plausible biological relationships. Using Mendelian randomization, we identified causal associations for 18 proteins across 19 phenotypes, of which 10 have additional genetic colocalization evidence. We highlight proteins associated with a constellation of cardiometabolic traits including angiopoietin-related protein 7 (ANGPTL7) and Semaphorin 3F (SEMA3F). Conclusion: Through large-scale analysis of protein quantitative trait loci, we provide a comprehensive overview of common variants associated with plasma proteins. We highlight possible biological relationships which may serve as a basis for further investigation into possible causal roles in cardiometabolic diseases.

This article is based on the address given by the author at the 2021 virtual meeting of the American Society of Human Genetics (ASHG). The video of the original address can be found at the ASHG website.

Studying the interplay between genetic variation, epigenetic changes and regulation of gene expression in immune cells is important to understand the modification of cellular states in various conditions, including immune diseases. Here, we built cis maps of regulatory regions with coordinated activity – Cis Regulatory Domains (CRDs) – in neutrophils, monocytes and T cells. For this, we leveraged (i) whole-genome sequencing (WGS), (ii) chromatin immunoprecipitation sequencing (ChIP-seq), (iii) DNA methylation (450k arrays), and (iv) transcriptional profiles (RNA-seq) from the BLUEPRINT consortium, for up to 200 individuals.Our study uncovers 9287, 7666 and 5480 histone CRDs (hCRDs) and 6053, 6112, 5701 methyl CRDs (mCRDs) in monocytes, neutrophils and T-cells, respectively. We discovered 15294 hCRD-gene and 6185 mCRD-gene associations (5% FDR). Only 33% of hCRD-gene associations and 37% of mCRD-gene associations were shared between cell-types, revealing the dynamic nature of regulatory interactions and how similarly located regulatory regions modulate the activity of different genes on different cell types. We mapped Quantitative Trait Loci associated with CRD activity (CRD-QTLs) and found that 89% and 70% of these hCRDs and mCRDs are under genetic control highlighting the importance of genetic variation to study the coordination of cellular regulatory programs. We found CRD-QTLs to be enriched in celltype-specific transcription factor binding sites, such as SPI1, STAT3, RFX1, SOX4, ATF3 for neutrophils and monocytes and TCF4 and BCL11A for T-cells, in line with the Human protein Atlas.We integrated PCHi-C data, which showed that most significant associations discovered within gene-CRD associations and co-expressed genes associated with the same CRD, involving large genomic distances, tend to happen between genomic regions in close spatial proximity. Finally, we mapped trans regulatory associations between CRDs, which enabled the discovery of 207 trans-eQTLs across cell types. Overlapping our hits with trans eQTLs from eQTLGen Consortium meta-analysis in whole blood revealed 52 trans-eQTLs shared between the two studies. Overall, we show that mapping functional regulatory units using population genomics data allows discovering important mechanisms in the regulation of gene expression in immune cells and gain a greater understanding of cell-type specific regulatory mechanisms of immunity.

Nutritional disorders cause secondary osteoporosis as well as musculoskeletal metabolism dysfunction. The Christian Orthodox Church’s fasting in Mediterranean countries such Greece and Cyprus, or M.C.O.C. diet consists of self-restraint from food and/or food categories for 180 ± 19 days of total, especially of animal protein and dairy products. This case-control study attempts to investigate the effect of this fasting pattern on musculoskeletal metabolism and bone density. One hundred fasters (or M.C.O.C. diet followers; 68 women and 32 men, mean 59 ± 6.5) with 32 years average fasting time and 100 non-fasters (66 females and 34 males, mean 58.1 ± 6.8; the control group of Mediterranean diet followers) over the age of 50, including menopausal women, were interviewed, as well as physically and laboratory examined along with DEXA measurements of the L2-4 vertebrae and hips. Nutrition data gathered through a 3-day food record during a non-fasting period, while energy intakes calculated on a daily food consumption basis. Given the overall low incidence of osteoporosis in the Mediterranean diet, it appears to be more than three times higher than that of M.C.O.C. diet despite the periodic restriction of food intake of animal origin into a slightly hypothermic pattern, which in turn is characterized by increased consumption of multicolored vegetable foods. Abstinence from dairy products and meat does not adversely affect musculoskeletal metabolism or bone density. M.C.O.C. diet seems to be “healthy” eating habit for the musculoskeletal system, as future studies expected to confirm.

Aims/hypothesis It is well established that physical activity, abdominal ectopic fat and glycaemic regulation are related but the underlying structure of these relationships is unclear. The previously proposed twin-cycle hypothesis (TC) provides a mechanistic basis for impairment in glycaemic control through the interactions of substrate availability, substrate metabolism and abdominal ectopic fat accumulation. Here, we hypothesise that the effect of physical activity in glucose regulation is mediated by the twin-cycle. We aimed to examine this notion in the Innovative Medicines Initiative Diabetes Research on Patient Stratification (IMI DIRECT) Consortium cohorts comprised of participants with normal or impaired glucose regulation (cohort 1: N ≤ 920) or with recently diagnosed type 2 diabetes (cohort 2: N ≤ 435). Methods We defined a structural equation model that describes the TC and fitted this within the IMI DIRECT dataset. A second model, twin-cycle plus physical activity (TC-PA), to assess the extent to which the effects of physical activity in glycaemic regulation are mediated by components in the twin-cycle, was also fitted. Beta cell function, insulin sensitivity and glycaemic control were modelled from frequently sampled 75 g OGTTs (fsOGTTs) and mixed-meal tolerance tests (MMTTs) in participants without and with diabetes, respectively. Abdominal fat distribution was assessed using MRI, and physical activity through wrist-worn triaxial accelerometry. Results are presented as standardised beta coefficients, SE and p values, respectively. Results The TC and TC-PA models showed better fit than null models (TC: χ2 = 242, p = 0.004 and χ2 = 63, p = 0.001 in cohort 1 and 2, respectively; TC-PA: χ2 = 180, p = 0.041 and χ2 = 60, p = 0.008 in cohort 1 and 2, respectively). The association of physical activity with glycaemic control was primarily mediated by variables in the liver fat cycle. Conclusions/interpretation These analyses partially support the mechanisms proposed in the twin-cycle model and highlight mechanistic pathways through which insulin sensitivity and liver fat mediate the association between physical activity and glycaemic control.

While the recurrent 22q11.2 deletion is one of the strongest genetic risk factors for schizophrenia (SCZ), variability of its associated neuropsychiatric endophenotypes reflects its incomplete penetrance for psychosis development. To assess whether this phenotypic variability is linked to common variants associated with SCZ, we studied the association between SCZ polygenic risk score (PRS) and longitudinally acquired phenotypic information of the Swiss 22q11.2DS cohort (n = 97, 50% females, mean age 17.7 yr, mean visit interval 3.8 yr). The SCZ PRS with the best predictive performance was ascertained in the Estonian Biobank (n = 201,146) with LDpred. The infinitesimal SCZ PRS model showed the strongest capacity in discriminating SCZ cases from controls with one SD difference in SCZ PRS corresponding to an odds ratio (OR) of 1.73 (95% CI 1.57–1.90, P = 1.47 × 10−29). In 22q11.2 patients, random-effects ordinal regression modelling using longitudinal data showed SCZ PRS to have the strongest effect on social anhedonia (OR = 2.09, P = 0.0002), and occupational functioning (OR = 1.82, P = 0.0003) within the negative symptoms course, and dysphoric mood (OR = 2.00, P = 0.002) and stress intolerance (OR = 1.76, P = 0.0002) within the general symptoms course. Genetic liability for SCZ was additionally associated with full scale cognitive decline (β = –0.25, P = 0.02) and with longitudinal volumetric reduction of the right and left hippocampi (β = –0.28, P = 0.005; β = –0.23, P = 0.02, respectively). Our results indicate that the polygenic contribution to SCZ acts upon the threshold-lowering first hit (i.e., the deletion). It modifies the endophenotypes of 22q11.2DS and augments the derailment of developmental trajectories of negative and general symptoms, cognition, and hippocampal volume.

Plasma cell-free DNA (cfDNA) is a surrogate marker of neutrophil extracellular traps (NETs) that contribute to immunothrombosis. There is growing interest about the mechanisms underlying NET formation and elevated cfDNA, but little is known about the factors involved. We aimed to identify genes involved in the regulation of cfDNA levels using data from the Genetic Analysis of Idiopathic Thrombophilia (GAIT-2) Project.Imputed genotypes, whole blood RNA-Seq data, and plasma cfDNA quantification were available for 935 of the GAIT-2 participants from 35 families with idiopathic thrombophilia. We performed heritability and GWAS analysis for cfDNA. The heritability of cfDNA was 0.26 (p = 3.7 × 10−6), while the GWAS identified a significant association (rs1687391, p = 3.55 × 10−10) near the ORM1 gene, on chromosome 9. An eQTL (expression quantitative trait loci) analysis revealed a significant association between the lead GWAS variant and the expression of ORM1 in whole blood (p = 6.14 × 10−9). Additionally, ORM1 expression correlated with levels of cfDNA (p = 4.38 × 10−4). Finally, genetic correlation analysis between cfDNA and thrombosis identified a suggestive association (ρ g = 0.43, p = 0.089).All in all, we show evidence of the role of ORM1 in regulating cfDNA levels in plasma, which might contribute to the susceptibility to thrombosis through mechanisms of immunothrombosis.

Transposable elements (TEs) are interspersed repeats that contribute to more than half of the human genome, and TE-embedded regulatory sequences are increasingly recognized as major components of the human regulome. Perturbations of this system can contribute to tumorigenesis, but the impact of TEs on gene expression in cancer cells remains to be fully assessed. Here, we analyzed 275 normal colon and 276 colorectal cancer (CRC) samples from the SYSCOL colorectal cancer cohort and discovered 10,111 and 5,152 TE expression quantitative trait loci (eQTLs) in normal and tumor tissues, respectively. Amongst the latter, 376 were exclusive to CRC, likely driven by changes in methylation patterns. We identified that transcription factors are more enriched in tumor-specific TE-eQTLs than shared TE-eQTLs, indicating that TEs are more specifically regulated in tumor than normal. Using Bayesian Networks to assess the causal relationship between eQTL variants, TEs and genes, we identified that 1,758 TEs are mediators of genetic effect, altering the expression of 1,626 nearby genes significantly more in tumor compared to normal, of which 51 are cancer driver genes. We show that tumor-specific TE-eQTLs trigger the driver capability of TEs subsequently impacting expression of nearby genes. Collectively, our results highlight a global profile of a new class of cancer drivers, thereby enhancing our understanding of tumorigenesis and providing potential new candidate mechanisms for therapeutic target development.

Differences between sexes contribute to variation in the levels of fasting glucose and insulin. Epidemiological studies established a higher prevalence of impaired fasting glucose in men and impaired glucose tolerance in women, however, the genetic component underlying this phenomenon is not established. We assess sex-dimorphic (73,089/50,404 women and 67,506/47,806 men) and sex-combined (151,188/105,056 individuals) fasting glucose/fasting insulin genetic effects via genome-wide association study meta-analyses in individuals of European descent without diabetes. Here we report sex dimorphism in allelic effects on fasting insulin at IRS1 and ZNF12 loci, the latter showing higher RNA expression in whole blood in women compared to men. We also observe sex-homogeneous effects on fasting glucose at seven novel loci. Fasting insulin in women shows stronger genetic correlations than in men with waist-to-hip ratio and anorexia nervosa. Furthermore, waist-to-hip ratio is causally related to insulin resistance in women, but not in men. These results position dissection of metabolic and glycemic health sex dimorphism as a steppingstone for understanding differences in genetic effects between women and men in related phenotypes., Sex differences in fasting glucose and insulin have been identified, but the genetic loci underlying these differences have not. Here, the authors perform a meta-analysis of genome-wide association studies to detect sex-specific and sex-dimorphic loci associated with fasting glucose and insulin.

Schizophrenia is a polygenic psychiatric disorder with limited understanding about the mechanistic changes in gene expression regulation. To elucidate on this, we integrate interindividual variability of regulatory activity with gene expression and genotype data captured from the prefrontal cortex of 272 cases and controls. We show that regulatory element activity is structured into 10,936 and 10,376 cis-regulatory domains in cases and controls, respectively, which display distinct regulatory element coordination structures in both states. By studying the interplay among genetic variants, gene expression and cis-regulatory domains, we ascertain that changes in coordinated regulatory activity tag alterations in gene expression levels (p=8.62e-06, OR=1.60), unveil case-specific QTL effects, and identify regulatory machinery changes for genes affecting synaptic function and dendritic spine morphology in schizophrenia. Altogether, we show that accounting for coordinated regulatory activity provides a novel mechanistic approach to reduce the search space for unveiling genetically perturbed regulation of gene expression in schizophrenia.

Funder: Universität des Saarlandes (1036), Background: Understanding the influence of genetic variants on DNA methylation is fundamental for the interpretation of epigenomic data in the context of disease. There is a need for systematic approaches not only for determining methylation quantitative trait loci (methQTL), but also for discriminating general from cell type-specific effects. Results: Here, we present a two-step computational framework MAGAR (https://bioconductor.org/packages/MAGAR), which fully supports the identification of methQTLs from matched genotyping and DNA methylation data, and additionally allows for illuminating cell type-specific methQTL effects. In a pilot analysis, we apply MAGAR on data in four tissues (ileum, rectum, T cells, B cells) from healthy individuals and demonstrate the discrimination of common from cell type-specific methQTLs. We experimentally validate both types of methQTLs in an independent data set comprising additional cell types and tissues. Finally, we validate selected methQTLs located in the PON1, ZNF155, and NRG2 genes by ultra-deep local sequencing. In line with previous reports, we find cell type-specific methQTLs to be preferentially located in enhancer elements. Conclusions: Our analysis demonstrates that a systematic analysis of methQTLs provides important new insights on the influences of genetic variants to cell type-specific epigenomic variation.

Type 2 diabetes (T2D) and non-alcoholic fatty liver disease (NAFLD) often co-occur. Defining causal pathways underlying this relationship may help optimize the prevention and treatment of both diseases. Thus, we assessed the strength and magnitude of the putative causal pathways linking dysglycemia and fatty liver, using a combination of causal inference methods.Measures of glycemia, insulin dynamics, magnetic resonance imaging (MRI)-derived abdominal and liver fat content, serological biomarkers, lifestyle, and anthropometry were obtained in participants from the IMI DIRECT cohorts (n=795 with new onset T2D and 2234 individuals free from diabetes). UK Biobank (n=3641) was used for modelling and replication purposes. Bayesian networks were employed to infer causal pathways, with causal validation using two-sample Mendelian randomization.Bayesian networks fitted to IMI DIRECT data identified higher basal insulin secretion rate (BasalISR) and MRI-derived excess visceral fat (VAT) accumulation as the features of dysmetabolism most likely to cause liver fat accumulation; the unconditional probability of fatty liver (>5%) increased significantly when conditioning on high levels of BasalISR and VAT (by 23%, 32% respectively; 40% for both). Analyses in UK Biobank yielded comparable results. MR confirmed most causal pathways predicted by the Bayesian networks.Here, BasalISR had the highest causal effect on fatty liver predisposition, providing mechanistic evidence underpinning the established association of NAFLD and T2D. BasalISR may represent a pragmatic biomarker for NAFLD prediction in clinical practice.

Trait-associated genetic variants affect complex phenotypes primarily via regulatory mechanisms on the transcriptome. To investigate the genetics of gene expression, we performed cis- and trans-expression quantitative trait locus (eQTL) analyses using blood-derived expression from 31,684 individuals through the eQTLGen Consortium. We detected cis-eQTL for 88% of genes, and these were replicable in numerous tissues. Distal trans-eQTL (detected for 37% of 10,317 trait-associated variants tested) showed lower replication rates, partially due to low replication power and confounding by cell type composition. However, replication analyses in single-cell RNA-seq data prioritized intracellular trans-eQTL. Trans-eQTL exerted their effects via several mechanisms, primarily through regulation by transcription factors. Expression of 13% of the genes correlated with polygenic scores for 1,263 phenotypes, pinpointing potential drivers for those traits. In summary, this work represents a large eQTL resource, and its results serve as a starting point for in-depth interpretation of complex phenotypes.Analyses of expression profiles from whole blood of 31,684 individuals identify cis-expression quantitative trait loci (eQTL) effects for 88% of genes and trans-eQTL effects for 37% of trait-associated variants.

Influenza makes millions of people ill every year, placing a large burden on the healthcare system and the economy. Combination therapies with two or more compounds hindering different mechanisms of viral replication represent a suitable approach to not only enhance the effectiveness of the individual drugs, but also to reduce the likelihood of resistance emergence. To develop a novel treatment against influenza, we combine virucidal sialylated cyclodextrins with interferon lambda and demonstrate, in ex vivo human airway epithelia, that the two compounds inhibit the replication of a clinical H1N1 strain more efficiently when administered together rather than alone. We investigate the mechanism of action of the combined treatment and its effects on the host cells in absence of infection, by performing single cell RNA sequencing analysis. We find that both the single and combined treatments impair viral replication to different extents across distinct epithelial cell types. We also show that each cell type comprises multiple sub-types, whose proportions are altered by H1N1 infection, and assess the ability of the treatments to restore them. To the best of our knowledge this is the first study investigating the effectiveness of an antiviral therapy by transcriptomic studies at the single cell level.

Vaccines are thought to be the best available solution for controlling the ongoing SARS-CoV-2 pandemic. However, the emergence of vaccine-resistant strains may come too rapidly for current vaccine developments to alleviate the health, economic and social consequences of the pandemic. To quantify and characterize the risk of such a scenario, we created a SIR-derived model with initial stochastic dynamics of the vaccine-resistant strain to study the probability of its emergence and establishment. Using parameters realistically resembling SARS-CoV-2 transmission, we model a wave-like pattern of the pandemic and consider the impact of the rate of vaccination and the strength of non-pharmaceutical intervention measures on the probability of emergence of a resistant strain. As expected, we found that a fast rate of vaccination decreases the probability of emergence of a resistant strain. Counterintuitively, when a relaxation of non-pharmaceutical interventions happened at a time when most individuals of the population have already been vaccinated the probability of emergence of a resistant strain was greatly increased. Consequently, we show that a period of transmission reduction close to the end of the vaccination campaign can substantially reduce the probability of resistant strain establishment. Our results suggest that policymakers and individuals should consider maintaining non-pharmaceutical interventions and transmission-reducing behaviours throughout the entire vaccination period.

BackgroundUnderstanding the influence of genetic variants on DNA methylation is fundamental for the interpretation of epigenomic data in the context of disease. There is a need for systematic approaches not only for determining methylation quantitative trait loci (methQTL) but also for discriminating general from cell-type-specific effects.ResultsHere, we present a two-step computational framework MAGAR, which fully supports identification of methQTLs from matched genotyping and DNA methylation data, and additionally the identification of quantitative cell-type-specific methQTL effects. In a pilot analysis, we apply MAGAR on data in four tissues (ileum, rectum, T-cells, B-cells) from healthy individuals and demonstrate the discrimination of common from cell-type-specific methQTLs. We experimentally validate both types of methQTLs in an independent dataset comprising additional cell types and tissues. Finally, we validate selected methQTLs (PON1, ZNF155, NRG2) by ultra-deep local sequencing. In line with previous reports, we find cell-type-specific methQTLs to be preferentially located in enhancer elements.ConclusionsOur analysis demonstrates that a systematic analysis of methQTLs provides important new insights on the influences of genetic variants to cell-type-specific epigenomic variation.

Leydig cells (LC) are the main testicular androgen-producing cells. In eutherian mammals, two types of LCs emerge successively during testicular development, fetal Leydig cells (FLCs) and adult Leydig cells (ALCs). Both display significant differences in androgen production and regulation. Using bulk RNA sequencing, we compared the transcriptomes of both LC populations to characterise their specific transcriptional and functional features. Despite similar transcriptomic profiles, a quarter of the genes show significant variations in expression between FLCs and ALCs. Non-transcriptional events, such as alternative splicing was also observed, including a high rate of intron retention in FLCs compared to ALCs. The use of single-cell RNA sequencing data also allowed the identification of nine FLC-specific genes and 50 ALC-specific genes. Expression of the corticotropin-releasing hormone 1 (Crhr1) receptor and the ACTH receptor melanocortin type 2 receptor (Mc2r) specifically in FLCs suggests a dual regulation of steroidogenesis. The androstenedione synthesis by FLCs is stimulated by luteinizing hormone (LH), CRH and ACTH whereas the testosterone synthesis by ALCs is dependent exclusively on LH. Overall, our study provides a useful database to explore LC development and function.

Despite the importance of germ cell (GC) differentiation for sexual reproduction, the gene networks underlying their fate remain unclear. Here, we comprehensively characterize the gene expression dynamics during sex determination based on single-cell RNA sequencing of 14 914 XX and XY mouse GCs between embryonic days (E) 9.0 and 16.5. We found that XX and XY GCs diverge transcriptionally as early as E11.5 with upregulation of genes downstream of the bone morphogenic protein (BMP) and nodal/Activin pathways in XY and XX GCs, respectively. We also identified a sex-specific upregulation of genes associated with negative regulation of mRNA processing and an increase in intron retention consistent with a reduction in mRNA splicing in XY testicular GCs by E13.5. Using computational gene regulation network inference analysis, we identified sex-specific, sequential waves of putative key regulator genes during GC differentiation and revealed that the meiotic genes are regulated by positive and negative master modules acting in an antagonistic fashion. Finally, we found that rare adrenal GCs enter meiosis similarly to ovarian GCs but display altered expression of master genes controlling the female and male genetic programs, indicating that the somatic environment is important for GC function. Our data are available on a web platform and provide a molecular roadmap of GC sex determination at single-cell resolution, which will serve as a valuable resource for future studies of gonad development, function, and disease.

Genetic variants identified by genome-wide association studies can contribute to disease risk by altering the production and abundance of mRNA, proteins and other molecules. However, the interplay between molecular intermediaries that define the pathway from genetic variation to disease is not well understood. Here, we evaluated the shared genetic regulation of mRNA molecules, proteins and metabolites derived from whole blood from 3,029 human donors. We find abundant allelic heterogeneity, where multiple variants regulate a particular molecular phenotype, and pleiotropy, where a single variant was associated with multiple molecular phenotypes over multiple genomic regions. We find varying proportions of shared genetic regulation across phenotypes, highest between expression and proteins (66.6%). We were able to recapitulate a substantial proportion of gene expression genetic regulation in a diverse set of 44 tissues, with a median of 88% shared associations for blood expression and 22.3% for plasma proteins. Finally, the genetic and molecular associations were represented in networks including 2,828 known GWAS variants. One sub-network shows the trans relationship between rs149007767 and RTEN, and identifies GRB10 and IKZF1 as candidates mediating genes. Our work provides a roadmap to understanding molecular networks and deriving the underlying mechanism of action of GWAS variants across different molecular phenotypes.

BackgroundIn December 2020, the United Kingdom (UK) reported a SARS-CoV-2 Variant of Concern (VoC) which is now named B.1.1.7. Based on initial data from the UK and later data from other countries, this variant was estimated to have a transmission fitness advantage of around 40-80% [1, 2, 3].AimThis study aims to estimate the transmission fitness advantage and the effective reproductive number of B.1.1.7 through time based on data from Switzerland.MethodsWe generated whole genome sequences from 11.8% of all confirmed SARS-CoV-2 cases in Switzerland between 14 December 2020 and 11 March 2021. Based on these data, we determine the daily frequency of the B.1.1.7 variant and quantify the variant’s transmission fitness advantage on a national and a regional scale.ResultsWe estimate B.1.1.7 had a transmission fitness advantage of 43-52% compared to the other variants circulating in Switzerland during the study period. Further, we estimate B.1.1.7 had a reproductive number above 1 from 01 January 2021 until the end of the study period, compared to below 1 for the other variants. Specifically, we estimate the reproductive number for B.1.1.7 was 1.24 [1.07-1.41] from 01 January until 17 January 2021 and 1.18 [1.06-1.30] from 18 January until 01 March 2021 based on the whole genome sequencing data. From 10 March to 16 March 2021, once B.1.1.7 was dominant, we estimate the reproductive number was 1.14 [1.00-1.26] based on all confirmed cases. For reference, Switzerland applied more non-pharmaceutical interventions to combat SARS-CoV-2 on 18 January 2021 and lifted some measures again on 01 March 2021.ConclusionThe observed increase in B.1.1.7 frequency in Switzerland during the study period is as expected based on observations in the UK. In absolute numbers, B.1.1.7 increased exponentially with an estimated doubling time of around 2-3.5 weeks. To monitor the ongoing spread of B.1.1.7, our plots are available online.

Lifestyle choices significantly influence mental health in later life. In this study we investigated the effects of the Christian Orthodox Church (COC) fasting diet, which includes long-term regular abstinence from animal-based products for half the calendar year, on cognitive function and emotional wellbeing of healthy adults. Two groups of fasting and non-fasting individuals were evaluated regarding their cognitive performance and the presence of anxiety and depression using the Mini Mental Examination Scale, the Hamilton Anxiety Scale, and the Geriatric Depression Scale (GDS), respectively. Data on physical activity, smoking, and vitamin levels were collected and correlated with mental health scoring. Negative binomial regression was performed to examine differences in the GDS scores between the two groups. Significantly lower levels of anxiety (7.48 ± 4.98 vs. 9.71 ± 5.25, p &lt, 0.001) and depression (2.24 ± 1.77 vs. 3.5 ± 2.52, 0.001), along with better cognitive function (29.15 ± 0.79 vs. 28.64 ± 1.27, 0.001), were noticed in fasting compared to non-fasting individuals. GDS score was 31% lower (Incidence Rate Ratio: 0.69, 95% Confidence Interval: 0.56–0.85) in the fasting group compared to the control, while vitamin and ferrum levels did not differ. The COC fasting diet was found to have an independent positive impact on cognition and mood in middle-aged and elderly individuals.

Vaccines are thought to be the best available solution for controlling the ongoing SARS-CoV-2 pandemic1,2. However, the emergence of vaccine-resistant strains3–6 may come too rapidly for current vaccine developments to alleviate the health, economic and social consequences of the pandemic7,8. To quantify and characterize the risk of such a scenario, we created a SIR-derived model9,10 with initial stochastic dynamics of the vaccine-resistant strain to study the probability of its emergence and establishment. Using parameters realistically resembling SARS-CoV-2 transmission, we model a wave-like pattern of the pandemic and consider the impact of the rate of vaccination and the strength of non-pharmaceutical intervention measures on the probability of emergence of a resistant strain. We found a counterintuitive result that the highest probability for the establishment of the resistant strain comes at a time of reduced non-pharmaceutical intervention measures when most individuals of the population have been vaccinated. Consequently, we show that a period of transmission reduction close to the end of the vaccination campaign can substantially reduce the probability of resistant strain establishment. Our results suggest that policymakers and individuals should consider maintaining non-pharmaceutical interventions7,11,12 throughout the entire vaccination period.

Individuals experiencing socioeconomic disadvantage in childhood have a higher rate of inflammation-related diseases decades later. Little is known about the mechanisms linking early life experiences to the functioning of the immune system in adulthood. To address this, we explore the relationship across social-to-biological layers of early life social exposures on levels of adulthood inflammation and the mediating role of gene regulatory mechanisms, epigenetic and transcriptomic profiling from blood, in 2,329 individuals from two European cohort studies. Consistently across both studies, we find transcriptional activity explains a substantive proportion (78% and 26%) of the estimated effect of early life disadvantaged social exposures on levels of adulthood inflammation. Furthermore, we show that mechanisms other than cis DNA methylation may regulate those transcriptional fingerprints. These results further our understanding of social-to-biological transitions by pinpointing the role of gene regulation that cannot fully be explained by differential cis DNA methylation.

OBJECTIVE We investigated the processes underlying glycemic deterioration in type 2 diabetes (T2D). RESEARCH DESIGN AND METHODS A total of 732 recently diagnosed patients with T2D from the Innovative Medicines Initiative Diabetes Research on Patient Stratification (IMI DIRECT) study were extensively phenotyped over 3 years, including measures of insulin sensitivity (OGIS), β-cell glucose sensitivity (GS), and insulin clearance (CLIm) from mixed meal tests, liver enzymes, lipid profiles, and baseline regional fat from MRI. The associations between the longitudinal metabolic patterns and HbA1c deterioration, adjusted for changes in BMI and in diabetes medications, were assessed via stepwise multivariable linear and logistic regression. RESULTS Faster HbA1c progression was independently associated with faster deterioration of OGIS and GS and increasing CLIm; visceral or liver fat, HDL-cholesterol, and triglycerides had further independent, though weaker, roles (R2 = 0.38). A subgroup of patients with a markedly higher progression rate (fast progressors) was clearly distinguishable considering these variables only (discrimination capacity from area under the receiver operating characteristic = 0.94). The proportion of fast progressors was reduced from 56% to 8–10% in subgroups in which only one trait among OGIS, GS, and CLIm was relatively stable (odds ratios 0.07–0.09). T2D polygenic risk score and baseline pancreatic fat, glucagon-like peptide 1, glucagon, diet, and physical activity did not show an independent role. CONCLUSIONS Deteriorating insulin sensitivity and β-cell function, increasing insulin clearance, high visceral or liver fat, and worsening of the lipid profile are the crucial factors mediating glycemic deterioration of patients with T2D in the initial phase of the disease. Stabilization of a single trait among insulin sensitivity, β-cell function, and insulin clearance may be relevant to prevent progression.

International audience; AIM: The worldwide increase in obesity and type 2 diabetes (T2D) represents a major health challenge. Chronically altered lipids induced by obesity further promote the development of T2D, and the accumulation of toxic lipid metabolites in serum and peripheral organs may contribute to the diabetic phenotype. METHODS: To better understand the complex metabolic pattern of lean and obese T2D and non-T2D individuals, we analysed the lipid profile of human serum, skeletal muscle and visceral adipose tissue of two cohorts by systematic mass spectrometry-based lipid analysis. RESULTS: Lipid homeostasis was strongly altered in a disease- and tissue-specific manner, allowing us to define T2D signatures associated with obesity from those that were obesity independent. Lipid changes encompassed lyso-, diacyl- and ether-phospholipids. Moreover, strong changes in sphingolipids included cytotoxic 1-deoxyceramide accumulation in a disease-specific manner in serum and visceral adipose tissue. The high amounts of non-canonical 1-deoxyceramide present in human adipose tissue most likely come from cell-autonomous synthesis because 1-deoxyceramide production increased upon differentiation to adipocytes in mouse cell culture experiments. CONCLUSION: Taken together, the observed lipidome changes in obesity and T2D will facilitate the identification of T2D patient subgroups and represent an important step towards personalized medicine in diabetes.