The evolution of computational research in a data-centric world.

Computational data-centric research techniques play a prevalent and multi-disciplinary role in life science research. In the past, scientists in wet labs generated the data, and computational researchers focused on creating tools for the analysis of those data. Computational researchers are now becoming more independent and taking leadership roles within biomedical projects, leveraging the increased availability of public data. We are now able to generate vast amounts of data, and the challenge has shifted from data generation to data analysis. Here we discuss the pitfalls, challenges, and opportunities facing the field of data-centric research in biology. We discuss the evolving perception of computational data-driven research and its rise as an independent domain in biomedical research while also addressing the significant collaborative opportunities that arise from integrating computational research with experimental and translational biology. Additionally, we discuss the future of data-centric research and its applications across various areas of the biomedical field.

The Translational Landscape of the Human Heart.

Gene expression in human tissue has primarily been studied on the transcriptional level, largely neglecting translational regulation. Here, we analyze the translatomes of 80 human hearts to identify new translation events and quantify the effect of translational regulation. We show extensive translational control of cardiac gene expression, which is orchestrated in a process-specific manner. Translation downstream of predicted disease-causing protein-truncating variants appears to be frequent, suggesting inefficient translation termination. We identify hundreds of previously undetected microproteins, expressed from lncRNAs and circRNAs, for which we validate the protein products in vivo. The translation of microproteins is not restricted to the heart and prominent in the translatomes of human kidney and liver. We associate these microproteins with diverse cellular processes and compartments and find that many locate to the mitochondria. Importantly, dozens of microproteins are translated from lncRNAs with well-characterized noncoding functions, indicating previously unrecognized biology.

Complement genes contribute sex-biased vulnerability in diverse disorders.

Many common illnesses, for reasons that have not been identified, differentially affect men and women. For instance, the autoimmune diseases systemic lupus erythematosus (SLE) and Sjögren's syndrome affect nine times more women than men1, whereas schizophrenia affects men with greater frequency and severity relative to women2. All three illnesses have their strongest common genetic associations in the major histocompatibility complex (MHC) locus, an association that in SLE and Sjögren's syndrome has long been thought to arise from alleles of the human leukocyte antigen (HLA) genes at that locus3-6. Here we show that variation of the complement component 4 (C4) genes C4A and C4B, which are also at the MHC locus and have been linked to increased risk for schizophrenia7, generates 7-fold variation in risk for SLE and 16-fold variation in risk for Sjögren's syndrome among individuals with common C4 genotypes, with C4A protecting more strongly than C4B in both illnesses. The same alleles that increase risk for schizophrenia greatly reduce risk for SLE and Sjögren's syndrome. In all three illnesses, C4 alleles act more strongly in men than in women: common combinations of C4A and C4B generated 14-fold variation in risk for SLE, 31-fold variation in risk for Sjögren's syndrome, and 1.7-fold variation in schizophrenia risk among men (versus 6-fold, 15-fold and 1.26-fold variation in risk among women, respectively). At a protein level, both C4 and its effector C3 were present at higher levels in cerebrospinal fluid and plasma8,9 in men than in women among adults aged between 20 and 50 years, corresponding to the ages of differential disease vulnerability. Sex differences in complement protein levels may help to explain the more potent effects of C4 alleles in men, women's greater risk of SLE and Sjögren's syndrome and men's greater vulnerability to schizophrenia. These results implicate the complement system as a source of sexual dimorphism in vulnerability to diverse illnesses.

Fibroblasts as an in vitro model of circadian genetic and genomic studies.

Bipolar disorder (BD) is a heritable disorder characterized by shifts in mood that manifest in manic or depressive episodes. Clinical studies have identified abnormalities of the circadian system in BD patients as a hallmark of underlying pathophysiology. Fibroblasts are a well-established in vitro model for measuring circadian patterns. We set out to examine the underlying genetic architecture of circadian rhythm in fibroblasts, with the goal to assess its contribution to the polygenic nature of BD disease risk. We collected, from primary cell lines of 6 healthy individuals, temporal genomic features over a 48 h period from transcriptomic data (RNA-seq) and open chromatin data (ATAC-seq). The RNA-seq data showed that only a limited number of genes, primarily the known core clock genes such as ARNTL, CRY1, PER3, NR1D2 and TEF display circadian patterns of expression consistently across cell cultures. The ATAC-seq data identified that distinct transcription factor families, like those with the basic helix-loop-helix motif, were associated with regions that were increasing in accessibility over time. Whereas known glucocorticoid receptor target motifs were identified in those regions that were decreasing in accessibility. Further evaluation of these regions using stratified linkage disequilibrium score regression analysis failed to identify a significant presence of them in the known genetic architecture of BD, and other psychiatric disorders or neurobehavioral traits in which the circadian rhythm is affected. In this study, we characterize the biological pathways that are activated in this in vitro circadian model, evaluating the relevance of these processes in the context of the genetic architecture of BD and other disorders, highlighting its limitations and future applications for circadian genomic studies.

Accelerated brain aging as a biomarker for staging in bipolar disorder: an exploratory study.

BACKGROUND: Two established staging models outline the longitudinal progression in bipolar disorder (BD) based on episode recurrence or inter-episodic functioning. However, underlying neurobiological mechanisms and corresponding biomarkers remain unexplored. This study aimed to investigate if global and (sub)cortical brain structures, along with brain-predicted age difference (brain-PAD) reflect illness progression as conceptualized in these staging models, potentially identifying brain-PAD as a biomarker for BD staging. METHODS: In total, 199 subjects with bipolar-I-disorder and 226 control subjects from the Dutch Bipolar Cohort with a high-quality T1-weighted magnetic resonance imaging scan were analyzed. Global and (sub)cortical brain measures and brain-PAD (the difference between biological and chronological age) were estimated. Associations between individual brain measures and the stages of both staging models were explored. RESULTS: A higher brain-PAD (higher biological age than chronological age) correlated with an increased likelihood of being in a higher stage of the inter-episodic functioning model, but not in the model based on number of mood episodes. However, after correcting for the confounding factors lithium-use and comorbid anxiety, the association lost significance. Global and (sub)cortical brain measures showed no significant association with the stages. CONCLUSIONS: These results suggest that brain-PAD may be associated with illness progression as defined by impaired inter-episodic functioning. Nevertheless, the significance of this association changed after considering lithium-use and comorbid anxiety disorders. Further research is required to disentangle the intricate relationship between brain-PAD, illness stages, and lithium intake or anxiety disorders. This study provides a foundation for potentially using brain-PAD as a biomarker for illness progression.

Childhood abuse v. neglect and risk for major psychiatric disorders.

BACKGROUND: Childhood maltreatment (CM) is a strong risk factor for psychiatric disorders but serves in its current definitions as an umbrella for various fundamentally different childhood experiences. As first step toward a more refined analysis of the impact of CM, our objective is to revisit the relation of abuse and neglect, major subtypes of CM, with symptoms across disorders. METHODS: Three longitudinal studies of major depressive disorder (MDD, N = 1240), bipolar disorder (BD, N = 1339), and schizophrenia (SCZ, N = 577), each including controls (N = 881), were analyzed. Multivariate regression models were used to examine the relation between exposure to abuse, neglect, or their combination to the odds for MDD, BD, SCZ, and symptoms across disorders. Bidirectional Mendelian randomization (MR) was used to probe causality, using genetic instruments of abuse and neglect derived from UK Biobank data (N = 143 473). RESULTS: Abuse was the stronger risk factor for SCZ (OR 3.51, 95% CI 2.17-5.67) and neglect for BD (OR 2.69, 95% CI 2.09-3.46). Combined CM was related to increased risk exceeding additive effects of abuse and neglect for MDD (RERI = 1.4) and BD (RERI = 1.1). Across disorders, abuse was associated with hallucinations (OR 2.16, 95% CI 1.55-3.01) and suicide attempts (OR 2.16, 95% CI 1.55-3.01) whereas neglect was associated with agitation (OR 1.24, 95% CI 1.02-1.51) and reduced need for sleep (OR 1.64, 95% CI 1.08-2.48). MR analyses were consistent with a bidirectional causal effect of abuse with SCZ (IVWforward = 0.13, 95% CI 0.01-0.24). CONCLUSIONS: Childhood abuse and neglect are associated with different risks to psychiatric symptoms and disorders. Unraveling the origin of these differences may advance understanding of disease etiology and ultimately facilitate development of improved personalized treatment strategies.

Mega-analysis of association between obesity and cortical morphology in bipolar disorders: ENIGMA study in 2832 participants.

BACKGROUND: Obesity is highly prevalent and disabling, especially in individuals with severe mental illness including bipolar disorders (BD). The brain is a target organ for both obesity and BD. Yet, we do not understand how cortical brain alterations in BD and obesity interact. METHODS: We obtained body mass index (BMI) and MRI-derived regional cortical thickness, surface area from 1231 BD and 1601 control individuals from 13 countries within the ENIGMA-BD Working Group. We jointly modeled the statistical effects of BD and BMI on brain structure using mixed effects and tested for interaction and mediation. We also investigated the impact of medications on the BMI-related associations. RESULTS: BMI and BD additively impacted the structure of many of the same brain regions. Both BMI and BD were negatively associated with cortical thickness, but not surface area. In most regions the number of jointly used psychiatric medication classes remained associated with lower cortical thickness when controlling for BMI. In a single region, fusiform gyrus, about a third of the negative association between number of jointly used psychiatric medications and cortical thickness was mediated by association between the number of medications and higher BMI. CONCLUSIONS: We confirmed consistent associations between higher BMI and lower cortical thickness, but not surface area, across the cerebral mantle, in regions which were also associated with BD. Higher BMI in people with BD indicated more pronounced brain alterations. BMI is important for understanding the neuroanatomical changes in BD and the effects of psychiatric medications on the brain.

On definitions of signal duration, evaluated on close-range airgun signals.

In impact assessments for underwater noise, the duration of a transient signal is often expressed by the 90%-energy signal duration τ90 %. Consequently, the rms sound pressure is computed over this duration. Using a large set of measurements on marine-seismic airgun signals, it is shown that τ90 % is often very close to the interval between the primary and secondary pulse (the bubble period) or a small integer multiple thereof. In this situation τ90 % is a measure of the duration of the relative silence between primary and secondary peaks, which is not the intended measure. Rarely, τ90 % quantifies the duration of the main peak, leading to a much lower value of τ90 %. Since the number of peaks included in τ90 % is sensitive to the nature of the signal, relatively small differences in the signal lead to large differences in τ90 %, causing instability in any metric based on τ90 %, e.g., the rms sound pressure. Alternative metrics are proposed that do not exhibit these weaknesses. The consequences for the interpretation of sound pressure level of a transient signal, and the benefits of using a more stable metric than τ90 % are demonstrated.

Association between body mass index and subcortical brain volumes in bipolar disorders-ENIGMA study in 2735 individuals.

Individuals with bipolar disorders (BD) frequently suffer from obesity, which is often associated with neurostructural alterations. Yet, the effects of obesity on brain structure in BD are under-researched. We obtained MRI-derived brain subcortical volumes and body mass index (BMI) from 1134 BD and 1601 control individuals from 17 independent research sites within the ENIGMA-BD Working Group. We jointly modeled the effects of BD and BMI on subcortical volumes using mixed-effects modeling and tested for mediation of group differences by obesity using nonparametric bootstrapping. All models controlled for age, sex, hemisphere, total intracranial volume, and data collection site. Relative to controls, individuals with BD had significantly higher BMI, larger lateral ventricular volume, and smaller volumes of amygdala, hippocampus, pallidum, caudate, and thalamus. BMI was positively associated with ventricular and amygdala and negatively with pallidal volumes. When analyzed jointly, both BD and BMI remained associated with volumes of lateral ventricles  and amygdala. Adjusting for BMI decreased the BD vs control differences in ventricular volume. Specifically, 18.41% of the association between BD and ventricular volume was mediated by BMI (Z = 2.73, p = 0.006). BMI was associated with similar regional brain volumes as BD, including lateral ventricles, amygdala, and pallidum. Higher BMI may in part account for larger ventricles, one of the most replicated findings in BD. Comorbidity with obesity could explain why neurostructural alterations are more pronounced in some individuals with BD. Future prospective brain imaging studies should investigate whether obesity could be a modifiable risk factor for neuroprogression.

Diagnosis of bipolar disorders and body mass index predict clustering based on similarities in cortical thickness-ENIGMA study in 2436 individuals.

AIMS: Rates of obesity have reached epidemic proportions, especially among people with psychiatric disorders. While the effects of obesity on the brain are of major interest in medicine, they remain markedly under-researched in psychiatry. METHODS: We obtained body mass index (BMI) and magnetic resonance imaging-derived regional cortical thickness, surface area from 836 bipolar disorders (BD) and 1600 control individuals from 14 sites within the ENIGMA-BD Working Group. We identified regionally specific profiles of cortical thickness using K-means clustering and studied clinical characteristics associated with individual cortical profiles. RESULTS: We detected two clusters based on similarities among participants in cortical thickness. The lower thickness cluster (46.8% of the sample) showed thinner cortex, especially in the frontal and temporal lobes and was associated with diagnosis of BD, higher BMI, and older age. BD individuals in the low thickness cluster were more likely to have the diagnosis of bipolar disorder I and less likely to be treated with lithium. In contrast, clustering based on similarities in the cortical surface area was unrelated to BD or BMI and only tracked age and sex. CONCLUSIONS: We provide evidence that both BD and obesity are associated with similar alterations in cortical thickness, but not surface area. The fact that obesity increased the chance of having low cortical thickness could explain differences in cortical measures among people with BD. The thinner cortex in individuals with higher BMI, which was additive and similar to the BD-associated alterations, may suggest that treating obesity could lower the extent of cortical thinning in BD.

Long-term follow-up including extensive complement analysis of a pediatric C3 glomerulopathy cohort.

BACKGROUND: C3 glomerulopathy (C3G) is a rare kidney disorder characterized by predominant glomerular depositions of complement C3. C3G can be subdivided into dense deposit disease (DDD) and C3 glomerulonephritis (C3GN). This study describes the long-term follow-up with extensive complement analysis of 29 Dutch children with C3G. METHODS: Twenty-nine C3G patients (19 DDD, 10 C3GN) diagnosed between 1992 and 2014 were included. Clinical and laboratory findings were collected at presentation and during follow-up. Specialized assays were used to detect rare variants in complement genes and measure complement-directed autoantibodies and biomarkers in blood. RESULTS: DDD patients presented with lower estimated glomerular filtration rate (eGFR). C3 nephritic factors (C3NeFs) were detected in 20 patients and remained detectable over time despite immunosuppressive treatment. At presentation, low serum C3 levels were detected in 84% of all patients. During follow-up, in about 50% of patients, all of them C3NeF-positive, C3 levels remained low. Linear mixed model analysis showed that C3GN patients had higher soluble C5b-9 (sC5b-9) and lower properdin levels compared to DDD patients. With a median follow-up of 52 months, an overall benign outcome was observed with only six patients with eGFR below 90 ml/min/1.73 m2 at last follow-up. CONCLUSIONS: We extensively described clinical and laboratory findings including complement features of an exclusively pediatric C3G cohort. Outcome was relatively benign, persistent low C3 correlated with C3NeF presence, and C3GN was associated with higher sC5b-9 and lower properdin levels. Prospective studies are needed to further elucidate the pathogenic mechanisms underlying C3G and guide personalized medicine with complement therapeutics.

Genetic analysis of activity, brain and behavioral associations in extended families with heavy genetic loading for bipolar disorder.

BACKGROUND: Disturbed sleep and activity are prominent features of bipolar disorder type I (BP-I). However, the relationship of sleep and activity characteristics to brain structure and behavior in euthymic BP-I patients and their non-BP-I relatives is unknown. Additionally, underlying genetic relationships between these traits have not been investigated. METHODS: Relationships between sleep and activity phenotypes, assessed using actigraphy, with structural neuroimaging (brain) and cognitive and temperament (behavior) phenotypes were investigated in 558 euthymic individuals from multi-generational pedigrees including at least one member with BP-I. Genetic correlations between actigraphy-brain and actigraphy-behavior associations were assessed, and bivariate linkage analysis was conducted for trait pairs with evidence of shared genetic influences. RESULTS: More physical activity and longer awake time were significantly associated with increased brain volumes and cortical thickness, better performance on neurocognitive measures of long-term memory and executive function, and less extreme scores on measures of temperament (impulsivity, cyclothymia). These associations did not differ between BP-I patients and their non-BP-I relatives. For nine activity-brain or activity-behavior pairs there was evidence for shared genetic influence (genetic correlations); of these pairs, a suggestive bivariate quantitative trait locus on chromosome 7 for wake duration and verbal working memory was identified. CONCLUSIONS: Our findings indicate that increased physical activity and more adequate sleep are associated with increased brain size, better cognitive function and more stable temperament in BP-I patients and their non-BP-I relatives. Additionally, we found evidence for pleiotropy of several actigraphy-behavior and actigraphy-brain phenotypes, suggesting a shared genetic basis for these traits.

Nonribosomal peptide synthetases and their biotechnological potential in Penicillium rubens.

Nonribosomal peptide synthetases (NRPS) are large multimodular enzymes that synthesize a diverse variety of peptides. Many of these are currently used as pharmaceuticals, thanks to their activity as antimicrobials (penicillin, vancomycin, daptomycin, echinocandin), immunosuppressant (cyclosporin) and anticancer compounds (bleomycin). Because of their biotechnological potential, NRPSs have been extensively studied in the past decades. In this review, we provide an overview of the main structural and functional features of these enzymes, and we consider the challenges and prospects of engineering NRPSs for the synthesis of novel compounds. Furthermore, we discuss secondary metabolism and NRP synthesis in the filamentous fungus Penicillium rubens and examine its potential for the production of novel and modified ß-lactam antibiotics.

Lessons learned from the Amsterdam Cohort Studies among people who use drugs: a historical perspective.

The Netherlands is well known for its early adoption of harm reduction (HR) programs at the height of its heroin crisis in the 1970s/1980s, including the implementation of the first needle and syringe program worldwide. In this manuscript, we describe how the Amsterdam Cohort Studies (ACS) among people who use drugs (PWUD) was conceived within the context of the Dutch HR approach, including the challenges scientists faced while establishing this cohort. This required striking a balance between public health and individual benefit, solving research dilemmas in the face of uncertainty, developing controversial innovative and cutting-edge interventions, which changed the prevention landscape for PWUD, and using longitudinal cohort data to provide unique insights. Studies from the ACS covering follow-up between 1985 and 2016 revealed that participation in both opioid agonist therapy and needle and syringe programs led to a major decrease in the risk of HIV and hepatitis B and C infection acquisition. ACS data have shown that the observed decrease in incidence also likely included shifts in drug markets and drug culture over time, selective mortality among those with the highest levels of risk behaviour, demographic changes of the PWUD population, and progression of the HIV and HCV epidemics. Moreover, HR programs in the Netherlands provided services beyond care for drug use, such as social support and welfare services, likely contributing to its success in curbing the HIV and viral hepatitis epidemics, increasing access and retention to HIV and HCV care and ultimately decreases in overdose mortality over time. Given the low coverage of HR programs in certain regions, it is unsurprising that continued HIV and HCV outbreaks occur and that transmission is ongoing in many countries worldwide. If we aim to reach the World Health Organization viral hepatitis and HIV elimination targets in 2030, as well as to improve the life of PWUD beyond infection risk, comprehensive HR programs need to be integrated as a part of prevention services, as in the Netherlands. We should use the evidence generated by longstanding cohorts, including the ACS, as a basis for which implementation and improved coverage of integrated HR services can be achieved for PWUD worldwide.

Integrative analysis of Dupuytren's disease identifies novel risk locus and reveals a shared genetic etiology with BMI.

Dupuytren's disease is a common inherited tissue-specific fibrotic disorder, characterized by progressive and irreversible fibroblastic proliferation affecting the palmar fascia of the hand. Although genome-wide association study (GWAS) have identified 24 genomic regions associated with Dupuytrens risk, the biological mechanisms driving signal at these regions remain elusive. We identify potential biological mechanisms for Dupuytren's disease by integrating the most recent, largest GWAS (3,871 cases and 4,686 controls) with eQTLs (47 tissue panels from five consortia, total n = 3,975) to perform a transcriptome-wide association study. We identify 43 tissue-specific gene associations with Dupuytren's risk, including one in a novel risk region. We also estimate the genome-wide genetic correlation between Dupuytren's disease and 45 complex traits and find significant genetic correlations between Dupuytren's disease and body mass index (BMI), type II diabetes, triglycerides, and high-density lipoprotein (HDL), suggesting a shared genetic etiology between these traits. We further examine local genetic correlation to identify 8 and 3 novel regions significantly correlated with BMI and HDL respectively. Our results are consistent with previous epidemiological findings showing that lower BMI increases risk for Dupuytren's disease. These 12 novel risk regions provide new insight into the biological mechanisms of Dupuytren's disease and serve as a starting point for functional validation.

Double hits in schizophrenia.

The co-occurrence of a copy number variant (CNV) and a functional variant on the other allele may be a relevant genetic mechanism in schizophrenia. We hypothesized that the cumulative burden of such double hits-in particular those composed of a deletion and a coding single-nucleotide variation (SNV)-is increased in patients with schizophrenia. We combined CNV data with coding variants data in 795 patients with schizophrenia and 474 controls. To limit false CNV-detection, only CNVs called by two algorithms were included. CNV-affected genes were subsequently examined for coding SNVs, which we termed "CNV-SNVs." Correcting for total queried sequence, we assessed the CNV-SNV-burden and the combined predicted deleterious effect. We estimated P-values by permutation of the phenotype. We detected 105 CNV-SNVs; 67 in duplicated and 38 in deleted genic sequence. Although the difference in CNV-SNVs rates was not significant, the combined deleteriousness inferred by CNV-SNVs in deleted sequence was almost 4-fold higher in cases compared with controls (nominal P = 0.009). This effect may be driven by a higher number of CNV-SNVs and/or by a higher degree of predicted deleteriousness of CNV-SNVs. No such effect was observed for duplications. We provide early evidence that deletions co-occurring with a functional variant may be relevant, albeit of modest impact, for the genetic etiology of schizophrenia. Large-scale consortium studies are required to validate our findings. Sequence-based analyses would provide the best resolution for detection of CNVs as well as coding variants genome-wide.

Functionally distinct ERAP1 and ERAP2 are a hallmark of HLA-A29-(Birdshot) Uveitis.

Birdshot Uveitis (Birdshot) is a rare eye condition that affects HLA-A29-positive individuals and could be considered a prototypic member of the recently proposed 'MHC-I (major histocompatibility complex class I)-opathy' family. Genetic studies have pinpointed the endoplasmic reticulum aminopeptidase (ERAP1) and (ERAP2) genes as shared associations across MHC-I-opathies, which suggests ERAP dysfunction may be a root cause for MHC-I-opathies. We mapped the ERAP1 and ERAP2 haplotypes in 84 Dutch cases and 890 controls. We identified association at variant rs10044354, which mediated a marked increase in ERAP2 expression. We also identified and cloned an independently associated ERAP1 haplotype (tagged by rs2287987) present in more than half of the cases; this ERAP1 haplotype is also the primary risk and protective haplotype for other MHC-I-opathies. We show that the risk ERAP1 haplotype conferred significantly altered expression of ERAP1 isoforms in transcriptomic data (n = 360), resulting in lowered protein expression and distinct enzymatic activity. Both the association for rs10044354 (meta-analysis: odds ratio (OR) [95% CI]=2.07[1.58-2.71], P = 1.24 × 10(-7)) and rs2287987 (OR[95% CI]: =2.01[1.51-2.67], P = 1.41 × 10(-6)) replicated and showed consistent direction of effect in an independent Spanish cohort of 46 cases and 2103 controls. In both cohorts, the combined rs2287987-rs10044354 haplotype associated with Birdshot more strongly than either variant alone [meta-analysis: P=3.9 × 10(-9)]. Finally, we observed that ERAP2 protein expression is dependent on the ERAP1 background across three European populations (n = 3353). In conclusion, a functionally distinct combination of ERAP1 and ERAP2 are a hallmark of Birdshot and provide rationale for strategies designed to correct ERAP function for treatment of Birdshot and MHC-I-opathies more broadly.

Impact of Classical Strain Improvement of Penicillium rubens on Amino Acid Metabolism during ß-Lactam Production.

To produce high levels of ß-lactams, the filamentous fungus Penicillium rubens (previously named Penicillium chrysogenum) has been subjected to an extensive classical strain improvement (CSI) program during the last few decades. This has led to the accumulation of many mutations that were spread over the genome. Detailed analysis reveals that several mutations targeted genes that encode enzymes involved in amino acid metabolism, in particular biosynthesis of l-cysteine, one of the amino acids used for ß-lactam production. To examine the impact of the mutations on enzyme function, the respective genes with and without the mutations were cloned and expressed in Escherichia coli, purified, and enzymatically analyzed. Mutations severely impaired the activities of a threonine and serine deaminase, and this inactivates metabolic pathways that compete for l-cysteine biosynthesis. Tryptophan synthase, which converts l-serine into l-tryptophan, was inactivated by a mutation, whereas a mutation in 5-aminolevulinate synthase, which utilizes glycine, was without an effect. Importantly, CSI caused increased expression levels of a set of genes directly involved in cysteine biosynthesis. These results suggest that CSI has resulted in improved cysteine biosynthesis by the inactivation of the enzymatic conversions that directly compete for resources with the cysteine biosynthetic pathway, consistent with the notion that cysteine is a key component during penicillin production.IMPORTANCEPenicillium rubens is an important industrial producer of ß-lactam antibiotics. High levels of penicillin production were enforced through extensive mutagenesis during a classical strain improvement (CSI) program over 70 years. Several mutations targeted amino acid metabolism and resulted in enhanced l-cysteine biosynthesis. This work provides a molecular explanation for the interrelation between secondary metabolite production and amino acid metabolism and how classical strain improvement has resulted in improved production strains.

Whole blood transcriptome analysis in bipolar disorder reveals strong lithium effect.

BACKGROUND: Bipolar disorder (BD) is a highly heritable mood disorder with complex genetic architecture and poorly understood etiology. Previous transcriptomic BD studies have had inconsistent findings due to issues such as small sample sizes and difficulty in adequately accounting for confounders like medication use. METHODS: We performed a differential expression analysis in a well-characterized BD case-control sample (Nsubjects = 480) by RNA sequencing of whole blood. We further performed co-expression network analysis, functional enrichment, and cell type decomposition, and integrated differentially expressed genes with genetic risk. RESULTS: While we observed widespread differential gene expression patterns between affected and unaffected individuals, these effects were largely linked to lithium treatment at the time of blood draw (FDR < 0.05, Ngenes = 976) rather than BD diagnosis itself (FDR < 0.05, Ngenes = 6). These lithium-associated genes were enriched for cell signaling and immune response functional annotations, among others, and were associated with neutrophil cell-type proportions, which were elevated in lithium users. Neither genes with altered expression in cases nor in lithium users were enriched for BD, schizophrenia, and depression genetic risk based on information from genome-wide association studies, nor was gene expression associated with polygenic risk scores for BD. CONCLUSIONS: These findings suggest that BD is associated with minimal changes in whole blood gene expression independent of medication use but emphasize the importance of accounting for medication use and cell type heterogeneity in psychiatric transcriptomic studies. The results of this study add to mounting evidence of lithium's cell signaling and immune-related mechanisms.

Exploring the clinical utility of two staging models for bipolar disorder.

OBJECTIVE: To assess the clinical utility of two staging models for bipolar disorder by examining distribution, correlation, and the relationship to external criteria. These are primarily defined by the recurrence of mood episodes (model A), or by intra-episodic functioning (model B). METHODS: In the Dutch Bipolar Cohort, stages according to models A and B were assigned to all patients with bipolar-I-disorder (BD-I; N = 1396). The dispersion of subjects over the stages was assessed and the association between the two models calculated. For both models, change in several clinical markers were concordant with the stage was investigated. RESULTS: Staging was possible in 87% of subjects for model A and 75% for model B. For model A, 1079 participants (93%) were assigned to stage 3c (recurrent episodes). Subdividing stage 3c with cut-offs at 5 and 10 episodes resulted in subgroups containing 242, 510, and 327 subjects. For model B, most participants were assigned to stage II (intra-episodic symptoms, N = 431 (41%)) or stage III (inability to work, N = 451 (43%)). A low association between models was found. For both models, the clinical markers "age at onset," "treatment resistance," and "episode acceleration" changed concordant with the stages. CONCLUSION: The majority of patients with BD-I clustered in recurrent stage 3 of Model A. Model B showed a larger dispersion. The stepwise change in several clinical markers supports the construct validity of both models. Combining the two staging models and sub-differentiating the recurrent stage into categories with cut-offs at 5 and 10 lifetime episodes improves the clinical utility of staging for individual patients.