Large-scale proximity extension assay reveals CSF midkine and DOPA decarboxylase as supportive diagnostic biomarkers for Parkinson's disease

BACKGROUND

There is a need for biomarkers to support an accurate diagnosis of Parkinson's disease (PD). Cerebrospinal fluid (CSF) has been a successful biofluid for finding neurodegenerative biomarkers, and modern highly sensitive multiplexing methods offer the possibility to perform discovery studies. Using a large-scale multiplex proximity extension assay (PEA) approach, we aimed to discover novel diagnostic protein biomarkers allowing accurate discrimination of PD from both controls and atypical Parkinsonian disorders (APD).

METHODS

CSF from patients with PD, corticobasal syndrome (CBS), progressive supranuclear palsy (PSP), multiple system atrophy and controls, were analysed with Olink PEA panels. Three cohorts were used in this study, comprising 192, 88 and 36 cases, respectively. All samples were run on the Cardiovascular II, Oncology II and Metabolism PEA panels.

RESULTS

Our analysis revealed that 26 and 39 proteins were differentially expressed in the CSF of test and validation PD cohorts, respectively, compared to controls. Among them, 6 proteins were changed in both cohorts. Midkine (MK) was increased in PD with the strongest effect size and results were validated with ELISA. Another most increased protein in PD, DOPA decarboxylase (DDC), which catalyses the decarboxylation of DOPA (L-3,4-dihydroxyphenylalanine) to dopamine, was strongly correlated with dopaminergic treatment. Moreover, Kallikrein 10 was specifically changed in APD compared with both PD and controls, but unchanged between PD and controls. Wnt inhibitory factor 1 was consistently downregulated in CBS and PSP patients in two independent cohorts.

CONCLUSIONS

Using the large-scale PEA approach, we have identified potential novel PD diagnostic biomarkers, most notably MK and DDC, in the CSF of PD patients.

Differential associations of APOE-ε2 and APOE-ε4 alleles with PET-measured amyloid-β and tau deposition in older individuals without dementia

PURPOSE

To examine associations between the APOE-ε2 and APOE-ε4 alleles and core Alzheimer's disease (AD) pathological hallmarks as measured by amyloid-β (Aβ) and tau PET in older individuals without dementia.

METHODS

We analyzed data from 462 ADNI participants without dementia who underwent Aβ ([18F]florbetapir or [18F]florbetaben) and tau ([18F]flortaucipir) PET, structural MRI, and cognitive testing. Employing APOE-ε3 homozygotes as the reference group, associations between APOE-ε2 and APOE-ε4 carriership with global Aβ PET and regional tau PET measures (entorhinal cortex (ERC), inferior temporal cortex, and Braak-V/VI neocortical composite regions) were investigated using linear regression models. In a subset of 156 participants, we also investigated associations between APOE genotype and regional tau accumulation over time using linear mixed models. Finally, we assessed whether Aβ mediated the cross-sectional and longitudinal associations between APOE genotype and tau.

RESULTS

Compared to APOE-ε3 homozygotes, APOE-ε2 carriers had lower global Aβ burden (βstd [95% confidence interval (CI)]: - 0.31 [- 0.45, - 0.16], p = 0.034) but did not differ on regional tau burden or tau accumulation over time. APOE-ε4 participants showed higher Aβ (βstd [95%CI]: 0.64 [0.42, 0.82], p < 0.001) and tau burden (βstd range: 0.27-0.51, all p < 0.006). In mediation analyses, APOE-ε4 only retained an Aβ-independent effect on tau in the ERC. APOE-ε4 showed a trend towards increased tau accumulation over time in Braak-V/VI compared to APOE-ε3 homozygotes (βstd [95%CI]: 0.10 [- 0.02, 0.18], p = 0.11), and this association was fully mediated by baseline Aβ.

CONCLUSION

Our data suggest that the established protective effect of the APOE-ε2 allele against developing clinical AD is primarily linked to resistance against Aβ deposition rather than tau pathology.

False lumen pressure estimation in type B aortic dissection using 4D flow cardiovascular magnetic resonance: comparisons with aortic growth

BACKGROUND

Chronic type B aortic dissection (TBAD) is associated with poor long-term outcome, and accurate risk stratification tools remain lacking. Pressurization of the false lumen (FL) has been recognized as central in promoting aortic growth. Several surrogate imaging-based metrics have been proposed to assess FL hemodynamics; however, their relationship to enlarging aortic dimensions remains unclear. We investigated the association between aortic growth and three cardiovascular magnetic resonance (CMR)-derived metrics of FL pressurization: false lumen ejection fraction (FLEF), maximum systolic deceleration rate (MSDR), and FL relative pressure (FL ΔPmax).

METHODS

CMR/CMR angiography was performed in 12 patients with chronic dissection of the descending thoracoabdominal aorta, including contrast-enhanced CMR angiography and time-resolved three-dimensional phase-contrast CMR (4D Flow CMR). Aortic growth rate was calculated as the change in maximal aortic diameter between baseline and follow-up imaging studies over the time interval, with patients categorized as having either 'stable' (< 3 mm/year) or 'enlarging' (≥ 3 mm/year) growth. Three metrics relating to FL pressurization were defined as: (1) FLEF: the ratio between retrograde and antegrade flow at the TBAD entry tear, (2) MSDR: the absolute difference between maximum and minimum systolic acceleration in the proximal FL, and (3) FL ΔPmax: the difference in absolute pressure between aortic root and distal FL.

RESULTS

FLEF was higher in enlarging TBAD (49.0 ± 17.9% vs. 10.0 ± 11.9%, p = 0.002), whereas FL ΔPmax was lower (32.2 ± 10.8 vs. 57.2 ± 12.5 mmHg/m, p = 0.017). MSDR and conventional anatomic variables did not differ significantly between groups. FLEF showed positive (r = 0.78, p = 0.003) correlation with aortic growth rate whereas FL ΔPmax showed negative correlation (r = - 0.64, p = 0.026). FLEF and FL ΔPmax remained as independent predictors of aortic growth rate after adjusting for baseline aortic diameter.

CONCLUSION

Comparative analysis of three 4D flow CMR metrics of TBAD FL pressurization demonstrated that those that focusing on retrograde flow (FLEF) and relative pressure (FL ΔPmax) independently correlated with growth and differentiated patients with enlarging and stable descending aortic dissections. These results emphasize the highly variable nature of aortic hemodynamics in TBAD patients, and suggest that 4D Flow CMR derived metrics of FL pressurization may be useful to separate patients at highest and lowest risk for progressive aortic growth and complications.

Transcriptional bursts explain autosomal random monoallelic expression and affect allelic imbalance

Transcriptional bursts render substantial biological noise in cellular transcriptomes. Here, we investigated the theoretical extent of allelic expression resulting from transcriptional bursting and how it compared to the amount biallelic, monoallelic and allele-biased expression observed in single-cell RNA-sequencing (scRNA-seq) data. We found that transcriptional bursting can explain the allelic expression patterns observed in single cells, including the frequent observations of autosomal monoallelic gene expression. Importantly, we identified that the burst frequency largely determined the fraction of cells with monoallelic expression, whereas the burst size had little effect on monoallelic observations. The high consistency between the bursting model predictions and scRNA-seq observations made it possible to assess the heterogeneity of a group of cells as their deviation in allelic observations from the expected. Finally, both burst frequency and size contributed to allelic imbalance observations and reinforced that studies of allelic imbalance can be confounded from the inherent noise in transcriptional bursting. Altogether, we demonstrate that allele-level transcriptional bursting renders widespread, although predictable, amounts of monoallelic and biallelic expression in single cells and cell populations.

Genes are transcribed into RNA and further translated into proteins. The maternal and paternal copy of each gene are typically transcribed independently, and transcription itself occur in discrete stochastic bursts (transcriptional bursts). Pioneering single-cell analysis of RNA across cells revealed abundant fluctuations in the amounts of maternal and paternal RNA in cells, with frequent observations of RNA from only the maternal or paternal gene copy (monoallelic expression). In this study, we investigated to which extent the observed monoallelic expression across single cells can be explained by transcriptional bursting. We demonstrate that the process of transcriptional bursting is sufficient to explain the amount of monoallelic expression, and we further demonstrate that the frequency of bursts mainly determines the frequency of monoallelic observations. Furthermore, we show that transcriptional bursts may lead to false positive observations of monoallelic expression across cell populations. Therefore, stochastic transcription renders large fluctuations in allelic origin of RNA in cells over time, including frequent monoallelic observations when profiling single cells.

Sources of variation in cell-type RNA-Seq profiles

Cell-type specific gene expression profiles are needed for many computational methods operating on bulk RNA-Seq samples, such as deconvolution of cell-type fractions and digital cytometry. However, the gene expression profile of a cell type can vary substantially due to both technical factors and biological differences in cell state and surroundings, reducing the efficacy of such methods. Here, we investigated which factors contribute most to this variation. We evaluated different normalization methods, quantified the variance explained by different factors, evaluated the effect on deconvolution of cell type fractions, and examined the differences between UMI-based single-cell RNA-Seq and bulk RNA-Seq. We investigated a collection of publicly available bulk and single-cell RNA-Seq datasets containing B and T cells, and found that the technical variation across laboratories is substantial, even for genes specifically selected for deconvolution, and this variation has a confounding effect on deconvolution. Tissue of origin is also a substantial factor, highlighting the challenge of using cell type profiles derived from blood with mixtures from other tissues. We also show that much of the differences between UMI-based single-cell and bulk RNA-Seq methods can be explained by the number of read duplicates per mRNA molecule in the single-cell sample. Our work shows the importance of either matching or correcting for technical factors when creating cell-type specific gene expression profiles that are to be used together with bulk samples.

Integration of social cues and individual experiences during instrumental avoidance learning

Learning to avoid harmful consequences can be a costly trial-and-error process. In such situations, social information can be leveraged to improve individual learning outcomes. Here, we investigated how participants used their own experiences and others' social cues to avoid harm. Participants made repeated choices between harmful and safe options, each with different probabilities of generating shocks, while also seeing the image of a social partner. Some partners made predictive gaze cues towards the harmful choice option while others cued an option at random, and did so using neutral or fearful facial expressions. We tested how learned social information about partner reliability transferred across contexts by letting participants encounter the same partner in multiple trial blocks while facing novel choice options. Participants' decisions were best explained by a reinforcement learning model that independently learned the probabilities of options being safe and of partners being reliable and combined these combined these estimates to generate choices. Advice from partners making a fearful facial expression influenced participants' decisions more than advice from partners with neutral expressions. Our results showed that participants made better decisions when facing predictive partners and that they cached and transferred partner reliability estimates into new blocks. Using simulations we show that participants' transfer of social information into novel contexts is better adapted to variable social environments where social partners may change their cuing strategy or become untrustworthy. Finally, we found no relation between autism questionnaire scores and performance in our task, but do find autism trait related differences in learning rate parameters.

Metabolic regulation of neurodifferentiation in the adult brain

Understanding the mechanisms behind neurodifferentiation in adults will be an important milestone in our quest to identify treatment strategies for cognitive disorders observed during our natural ageing or disease. It is now clear that the maturation of neural stem cells to neurones, fully integrated into neuronal circuits requires a complete remodelling of cellular metabolism, including switching the cellular energy source. Mitochondria are central for this transition and are increasingly seen as the regulatory hub in defining neural stem cell fate and neurodevelopment. This review explores our current knowledge of metabolism during adult neurodifferentiation.

Segregation within school classes: Detecting social clustering in choice data

We suggest a new method for detecting patterns of social clustering based on choice data. The method compares similar subjects within and between cohorts and thereby allows us to isolate the effect of peer influence from that of exogenous factors. Using this method on Norwegian register data, we address the question of whether students tend to cluster socially based on similar background. We find that common background correlates with making the same choices of curricular tracks, and that both exogenous preferences and peer influence matter. This applies to immigrant students from the same country, and, to some extent, to descendants of immigrants, but not to students from culturally similar countries. There are also small effects related to parents' education and income.

Defects of mitochondrial RNA turnover lead to the accumulation of double-stranded RNA in vivo

The RNA helicase SUV3 and the polynucleotide phosphorylase PNPase are involved in the degradation of mitochondrial mRNAs but their roles in vivo are not fully understood. Additionally, upstream processes, such as transcript maturation, have been linked to some of these factors, suggesting either dual roles or tightly interconnected mechanisms of mitochondrial RNA metabolism. To get a better understanding of the turn-over of mitochondrial RNAs in vivo, we manipulated the mitochondrial mRNA degrading complex in Drosophila melanogaster models and studied the molecular consequences. Additionally, we investigated if and how these factors interact with the mitochondrial poly(A) polymerase, MTPAP, as well as with the mitochondrial mRNA stabilising factor, LRPPRC. Our results demonstrate a tight interdependency of mitochondrial mRNA stability, polyadenylation and the removal of antisense RNA. Furthermore, disruption of degradation, as well as polyadenylation, leads to the accumulation of double-stranded RNAs, and their escape out into the cytoplasm is associated with an altered immune-response in flies. Together our results suggest a highly organised and inter-dependable regulation of mitochondrial RNA metabolism with far reaching consequences on cellular physiology.

Although a number of factors have been implemented in the turnover of mitochondrial (mt) DNA-derived transcripts, their exact functions and interplay with one another is not entirely clear. Several of these factors have been proposed to co-ordinately regulate both transcript maturation, as well as degradation, but the order of events during mitochondrial RNA turnover is less well understood. Using a range of different genetically modified Drosophila melanogaster models, we studied the involvement of the RNA helicase SUV3, the polynucleotide phosphorylase PNPase, the leucine-rich pentatricopeptide repeat motif-containing protein LRPPRC, and the mitochondrial RNA poly(A) polymerase MTPAP, in stabilisation, polyadenylation, and degradation of mitochondrial transcripts. Our results show a tight collaborative activity of these factors in vivo and reveal a clear hierarchical order of events leading to mitochondrial mRNA maturation. Furthermore, we demonstrate that the loss of SUV3, PNPase, or MTPAP leads to the accumulation of mitochondrial-derived antisense RNA in the cytoplasm of cells, which is associated with an altered immune-response in flies.

How adaptive plasticity evolves when selected against

Adaptive plasticity allows organisms to cope with environmental change, thereby increasing the population's long-term fitness. However, individual selection can only compare the fitness of individuals within each generation: if the environment changes more slowly than the generation time (i.e., a coarse-grained environment) a population will not experience selection for plasticity even if it is adaptive in the long-term. How does adaptive plasticity then evolve? One explanation is that, if competing alleles conferring different degrees of plasticity persist across multiple environments, natural selection between genetic lineages could select for adaptive plasticity (lineage selection). We show that adaptive plasticity can evolve even in the absence of such lineage selection. Instead, we propose that adaptive plasticity in coarse-grained environments evolves as a by-product of inefficient short-term natural selection: populations that rapidly evolve their phenotypes in response to selective pressures follow short-term optima, with the result that they have reduced long-term fitness across environments. Conversely, populations that accumulate limited genetic change within each environment evolve long-term adaptive plasticity even when plasticity incurs short-term costs. These results remain qualitatively similar regardless of whether we decrease the efficiency of natural selection by increasing the rate of environmental change or decreasing mutation rate, demonstrating that both factors act via the same mechanism. We demonstrate how this mechanism can be understood through the concept of learning rate. Our work shows how plastic responses that are costly in the short term, yet adaptive in the long term, can evolve as a by-product of inefficient short-term selection, without selection for plasticity at either the individual or lineage level.

SLIMEr: probing flexibility of lipid metabolism in yeast with an improved constraint-based modeling framework

BACKGROUND

A recurrent problem in genome-scale metabolic models (GEMs) is to correctly represent lipids as biomass requirements, due to the numerous of possible combinations of individual lipid species and the corresponding lack of fully detailed data. In this study we present SLIMEr, a formalism for correctly representing lipid requirements in GEMs using commonly available experimental data.

RESULTS

SLIMEr enhances a GEM with mathematical constructs where we Split Lipids Into Measurable Entities (SLIME reactions), in addition to constraints on both the lipid classes and the acyl chain distribution. By implementing SLIMEr on the consensus GEM of Saccharomyces cerevisiae, we can represent accurate amounts of lipid species, analyze the flexibility of the resulting distribution, and compute the energy costs of moving from one metabolic state to another.

CONCLUSIONS

The approach shows potential for better understanding lipid metabolism in yeast under different conditions. SLIMEr is freely available at https://github.com/SysBioChalmers/SLIMEr .

An efficient proteome-wide strategy for discovery and characterization of cellular nucleotide-protein interactions

Metabolite-protein interactions define the output of metabolic pathways and regulate many cellular processes. Although diseases are often characterized by distortions in metabolic processes, efficient means to discover and study such interactions directly in cells have been lacking. A stringent implementation of proteome-wide Cellular Thermal Shift Assay (CETSA) was developed and applied to key cellular nucleotides, where previously experimentally confirmed protein-nucleotide interactions were well recaptured. Many predicted, but never experimentally confirmed, as well as novel protein-nucleotide interactions were discovered. Interactions included a range of different protein families where nucleotides serve as substrates, products, co-factors or regulators. In cells exposed to thymidine, a limiting precursor for DNA synthesis, both dose- and time-dependence of the intracellular binding events for sequentially generated thymidine metabolites were revealed. Interactions included known cancer targets in deoxyribonucleotide metabolism as well as novel interacting proteins. This stringent CETSA based strategy will be applicable for a wide range of metabolites and will therefore greatly facilitate the discovery and studies of interactions and specificities of the many metabolites in human cells that remain uncharacterized.

Genome-wide meta-analysis of 158,000 individuals of European ancestry identifies three loci associated with chronic back pain

Back pain is the #1 cause of years lived with disability worldwide, yet surprisingly little is known regarding the biology underlying this symptom. We conducted a genome-wide association study (GWAS) meta-analysis of chronic back pain (CBP). Adults of European ancestry were included from 15 cohorts in the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) consortium, and from the UK Biobank interim data release. CBP cases were defined as those reporting back pain present for ≥3-6 months; non-cases were included as comparisons ("controls"). Each cohort conducted genotyping using commercially available arrays followed by imputation. GWAS used logistic regression models with additive genetic effects, adjusting for age, sex, study-specific covariates, and population substructure. The threshold for genome-wide significance in the fixed-effect inverse-variance weighted meta-analysis was p<5×10(-8). Suggestive (p<5×10(-7)) and genome-wide significant (p<5×10(-8)) variants were carried forward for replication or further investigation in the remaining UK Biobank participants not included in the discovery sample. The discovery sample comprised 158,025 individuals, including 29,531 CBP cases. A genome-wide significant association was found for the intronic variant rs12310519 in SOX5 (OR 1.08, p = 7.2×10(-10)). This was subsequently replicated in 283,752 UK Biobank participants not included in the discovery sample, including 50,915 cases (OR 1.06, p = 5.3×10(-11)), and exceeded genome-wide significance in joint meta-analysis (OR 1.07, p = 4.5×10(-19)). We found suggestive associations at three other loci in the discovery sample, two of which exceeded genome-wide significance in joint meta-analysis: an intergenic variant, rs7833174, located between CCDC26 and GSDMC (OR 1.05, p = 4.4×10(-13)), and an intronic variant, rs4384683, in DCC (OR 0.97, p = 2.4×10(-10)). In this first reported meta-analysis of GWAS for CBP, we identified and replicated a genetic locus associated with CBP (SOX5). We also identified 2 other loci that reached genome-wide significance in a 2-stage joint meta-analysis (CCDC26/GSDMC and DCC).

Simulations to benchmark time-varying connectivity methods for fMRI

There is a current interest in quantifying time-varying connectivity (TVC) based on neuroimaging data such as fMRI. Many methods have been proposed, and are being applied, revealing new insight into the brain's dynamics. However, given that the ground truth for TVC in the brain is unknown, many concerns remain regarding the accuracy of proposed estimates. Since there exist many TVC methods it is difficult to assess differences in time-varying connectivity between studies. In this paper, we present tvc_benchmarker, which is a Python package containing four simulations to test TVC methods. Here, we evaluate five different methods that together represent a wide spectrum of current approaches to estimating TVC (sliding window, tapered sliding window, multiplication of temporal derivatives, spatial distance and jackknife correlation). These simulations were designed to test each method's ability to track changes in covariance over time, which is a key property in TVC analysis. We found that all tested methods correlated positively with each other, but there were large differences in the strength of the correlations between methods. To facilitate comparisons with future TVC methods, we propose that the described simulations can act as benchmark tests for evaluation of methods. Using tvc_benchmarker researchers can easily add, compare and submit their own TVC methods to evaluate its performance.

In vitro co-cultures of human gut bacterial species as predicted from co-occurrence network analysis

Network analysis of large metagenomic datasets generated by current sequencing technologies can reveal significant co-occurrence patterns between microbial species of a biological community. These patterns can be analyzed in terms of pairwise combinations between all species comprising a community. Here, we construct a co-occurrence network for abundant microbial species encompassing the three dominant phyla found in human gut. This was followed by an in vitro evaluation of the predicted microbe-microbe co-occurrences, where we chose species pairs Bifidobacterium adolescentis and Bacteroides thetaiotaomicron, as well as Faecalibacterium prausnitzii and Roseburia inulinivorans as model organisms for our study. We then delineate the outcome of the co-cultures when equal distributions of resources were provided. The growth behavior of the co-culture was found to be dependent on the types of microbial species present, their specific metabolic activities, and resulting changes in the culture environment. Through this reductionist approach and using novel in vitro combinations of microbial species under anaerobic conditions, the results of this work will aid in the understanding and design of synthetic community formulations.

Decentralized colonoscopic surveillance with high patient compliance prevents hereditary and familial colorectal cancer

Although colonoscopic surveillance is recommended both for individuals with known hereditary colorectal cancer (HCRC) syndromes and those with a more moderate familial colorectal cancer (FCRC) history, the evidence for the benefits of surveillance is limited and surveillance practices vary. This study evaluates the preventive effect for individuals with a family history of CRC of decentralized colonoscopic surveillance with the guidance of a cancer prevention clinic. We performed a population based prospective study of 261 patients with HCRC or FCRC, recorded in the colonoscopic surveillance registry at the Cancer genetics clinic, University Hospital of Umeå, Sweden. Colonoscopic surveillance was conducted every second (HCRC) or fifth (FCRC) year at local hospitals in Northern Sweden. Main outcome measures were findings of high-risk adenomas (HRA) or CRC, and patient compliance to surveillance. Estimations of the expected numbers of CRC without surveillance were made. During a total of 1256 person years of follow-up, one case of CRC was found. The expected numbers of cancers in the absence of surveillance was between 9.5 and 10.5, resulting in a standardized incidence ratio, observed versus expected cases of CRC, between 0.10 (CI 95 % 0.0012-0.5299) and 0.11 (CI 95 % 0.0014-0.5857). No CRC mortality was reported, but three patients needed surgical intervention. HRA were found in 5.9 % (14/237) of the initial and in 3.4 % (12/356) of the follow-up colonoscopies. Patient compliance to the surveillance program was 90 % as 597 of the planned 662 colonoscopies were performed. The study concludes that colonoscopic surveillance with high patient compliance to the program is effective in preventing CRC when using a decentralized method for colonoscopy surveillance with the guidance of a cancer prevention clinic.

On Edge Exchangeable Random Graphs

We study a recent model for edge exchangeable random graphs introduced by Crane and Dempsey; in particular we study asymptotic properties of the random simple graph obtained by merging multiple edges. We study a number of examples, and show that the model can produce dense, sparse and extremely sparse random graphs. One example yields a power-law degree distribution. We give some examples where the random graph is dense and converges a.s. in the sense of graph limit theory, but also an example where a.s. every graph limit is the limit of some subsequence. Another example is sparse and yields convergence to a non-integrable generalized graphon defined on ( 0 , ∞ ) .

Identification of a Sjögren's syndrome susceptibility locus at OAS1 that influences isoform switching, protein expression, and responsiveness to type I interferons

Sjögren's syndrome (SS) is a common, autoimmune exocrinopathy distinguished by keratoconjunctivitis sicca and xerostomia. Patients frequently develop serious complications including lymphoma, pulmonary dysfunction, neuropathy, vasculitis, and debilitating fatigue. Dysregulation of type I interferon (IFN) pathway is a prominent feature of SS and is correlated with increased autoantibody titers and disease severity. To identify genetic determinants of IFN pathway dysregulation in SS, we performed cis-expression quantitative trait locus (eQTL) analyses focusing on differentially expressed type I IFN-inducible transcripts identified through a transcriptome profiling study. Multiple cis-eQTLs were associated with transcript levels of 2'-5'-oligoadenylate synthetase 1 (OAS1) peaking at rs10774671 (PeQTL = 6.05 × 10-14). Association of rs10774671 with SS susceptibility was identified and confirmed through meta-analysis of two independent cohorts (Pmeta = 2.59 × 10-9; odds ratio = 0.75; 95% confidence interval = 0.66-0.86). The risk allele of rs10774671 shifts splicing of OAS1 from production of the p46 isoform to multiple alternative transcripts, including p42, p48, and p44. We found that the isoforms were differentially expressed within each genotype in controls and patients with and without autoantibodies. Furthermore, our results showed that the three alternatively spliced isoforms lacked translational response to type I IFN stimulation. The p48 and p44 isoforms also had impaired protein expression governed by the 3' end of the transcripts. The SS risk allele of rs10774671 has been shown by others to be associated with reduced OAS1 enzymatic activity and ability to clear viral infections, as well as reduced responsiveness to IFN treatment. Our results establish OAS1 as a risk locus for SS and support a potential role for defective viral clearance due to altered IFN response as a genetic pathophysiological basis of this complex autoimmune disease.

Type 2 diabetes and obesity induce similar transcriptional reprogramming in human myocytes

BACKGROUND

Skeletal muscle is one of the primary tissues involved in the development of type 2 diabetes (T2D). The close association between obesity and T2D makes it difficult to isolate specific effects attributed to the disease alone. Therefore, here we set out to identify and characterize intrinsic properties of myocytes, associated independently with T2D or obesity.

METHODS

We generated and analyzed RNA-seq data from primary differentiated myotubes from 24 human subjects, using a factorial design (healthy/T2D and non-obese/obese), to determine the influence of each specific factor on genome-wide transcription. This setup enabled us to identify intrinsic properties, originating from muscle precursor cells and retained in the corresponding myocytes. Bioinformatic and statistical methods, including differential expression analysis, gene-set analysis, and metabolic network analysis, were used to characterize the different myocytes.

RESULTS

We found that the transcriptional program associated with obesity alone was strikingly similar to that induced specifically by T2D. We identified a candidate epigenetic mechanism, H3K27me3 histone methylation, mediating these transcriptional signatures. T2D and obesity were independently associated with dysregulated myogenesis, down-regulated muscle function, and up-regulation of inflammation and extracellular matrix components. Metabolic network analysis identified that in T2D but not obesity a specific metabolite subnetwork involved in sphingolipid metabolism was transcriptionally regulated.

CONCLUSIONS

Our findings identify inherent characteristics in myocytes, as a memory of the in vivo phenotype, without the influence from a diabetic or obese extracellular environment, highlighting their importance in the development of T2D.

Selective mitochondrial DNA degradation following double-strand breaks

Mitochondrial DNA (mtDNA) can undergo double-strand breaks (DSBs), caused by defective replication, or by various endogenous or exogenous sources, such as reactive oxygen species, chemotherapeutic agents or ionizing radiations. MtDNA encodes for proteins involved in ATP production, and maintenance of genome integrity following DSBs is thus of crucial importance. However, the mechanisms involved in mtDNA maintenance after DSBs remain unknown. In this study, we investigated the consequences of the production of mtDNA DSBs using a human inducible cell system expressing the restriction enzyme PstI targeted to mitochondria. Using this system, we could not find any support for DSB repair of mtDNA. Instead we observed a loss of the damaged mtDNA molecules and a severe decrease in mtDNA content. We demonstrate that none of the known mitochondrial nucleases are involved in the mtDNA degradation and that the DNA loss is not due to autophagy, mitophagy or apoptosis. Our study suggests that a still uncharacterized pathway for the targeted degradation of damaged mtDNA in a mitophagy/autophagy-independent manner is present in mitochondria, and might provide the main mechanism used by the cells to deal with DSBs.

Human ROBO1 regulates white matter structure in corpus callosum

The axon guidance receptor, Robo1, controls the pathfinding of callosal axons in mice. To determine whether the orthologous ROBO1 gene is involved in callosal development also in humans, we studied polymorphisms in the ROBO1 gene and variation in the white matter structure in the corpus callosum using both structural magnetic resonance imaging and diffusion tensor magnetic resonance imaging. We found that five polymorphisms in the regulatory region of ROBO1 were associated with white matter density in the posterior part of the corpus callosum pathways. One of the polymorphisms, rs7631357, was also significantly associated with the probability of connections to the parietal cortical regions. Our results demonstrate that human ROBO1 may be involved in the regulation of the structure and connectivity of posterior part of corpus callosum.

Longitudinal analysis of hepatic transcriptome and serum metabolome demonstrates altered lipid metabolism following the onset of hyperglycemia in spontaneously diabetic biobreeding rats

Type 1 diabetes is associated with abberations of fat metabolism before and after the clinical onset of disease. It has been hypothesized that the absence of the effect of insulin in the liver contributes to reduced hepatic fat synthesis. We measured hepatic gene expression and serum metabolites before and after the onset of hyperglycemia in a BioBreeding rat model of type 1 diabetes. Functional pathway annotation identified that lipid metabolism was differentially expressed in hyperglycemic rats and that these pathways significantly overlapped with genes regulated by insulin. 17 serum metabolites significantly changed in concentration. All but 2 of the identified metabolites had previously been reported in type 1 diabetes, and carbohydrates were overall the most upregulated class of metabolites. We conclude that lack of insulin in the liver contributes to the changes in fat metabolism observed in type 1 diabetes. Further studies are needed to understand the clinical consequences of a lack of insulin in the liver in patients with type 1 diabetes.

Organic materials database: An open-access online database for data mining

We present an organic materials database (OMDB) hosting thousands of Kohn-Sham electronic band structures, which is freely accessible online at http://omdb.diracmaterials.org. The OMDB focus lies on electronic structure, density of states and other properties for purely organic and organometallic compounds that are known to date. The electronic band structures are calculated using density functional theory for the crystal structures contained in the Crystallography Open Database. The OMDB web interface allows users to retrieve materials with specified target properties using non-trivial queries about their electronic structure. We illustrate the use of the OMDB and how it can become an organic part of search and prediction of novel functional materials via data mining techniques. As a specific example, we provide data mining results for metals and semiconductors, which are known to be rare in the class of organic materials.

Chemical delivery array with millisecond neurotransmitter release

Technologies that restore or augment dysfunctional neural signaling represent a promising route to deeper understanding and new therapies for neurological disorders. Because of the chemical specificity and subsecond signaling of the nervous system, these technologies should be able to release specific neurotransmitters at specific locations with millisecond resolution. We have previously demonstrated an organic electronic lateral electrophoresis technology capable of precise delivery of charged compounds, such as neurotransmitters. However, this technology, the organic electronic ion pump, has been limited to a single delivery point, or several simultaneously addressed outlets, with switch-on speeds of seconds. We report on a vertical neurotransmitter delivery device, configured as an array with individually controlled delivery points and a temporal resolution of 50 ms. This is achieved by supplementing lateral electrophoresis with a control electrode and an ion diode at each delivery point to allow addressing and limit leakage. By delivering local pulses of neurotransmitters with spatiotemporal dynamics approaching synaptic function, the high-speed delivery array promises unprecedented access to neural signaling and a path toward biochemically regulated neural prostheses.

The Impact of Endurance Training on Human Skeletal Muscle Memory, Global Isoform Expression and Novel Transcripts

Regularly performed endurance training has many beneficial effects on health and skeletal muscle function, and can be used to prevent and treat common diseases e.g. cardiovascular disease, type II diabetes and obesity. The molecular adaptation mechanisms regulating these effects are incompletely understood. To date, global transcriptome changes in skeletal muscles have been studied at the gene level only. Therefore, global isoform expression changes following exercise training in humans are unknown. Also, the effects of repeated interventions on transcriptional memory or training response have not been studied before. In this study, 23 individuals trained one leg for three months. Nine months later, 12 of the same subjects trained both legs in a second training period. Skeletal muscle biopsies were obtained from both legs before and after both training periods. RNA sequencing analysis of all 119 skeletal muscle biopsies showed that training altered the expression of 3,404 gene isoforms, mainly associated with oxidative ATP production. Fifty-four genes had isoforms that changed in opposite directions. Training altered expression of 34 novel transcripts, all with protein-coding potential. After nine months of detraining, no training-induced transcriptome differences were detected between the previously trained and untrained legs. Although there were several differences in the physiological and transcriptional responses to repeated training, no coherent evidence of an endurance training induced transcriptional skeletal muscle memory was found. This human lifestyle intervention induced differential expression of thousands of isoforms and several transcripts from unannotated regions of the genome. It is likely that the observed isoform expression changes reflect adaptational mechanisms and processes that provide the functional and health benefits of regular physical activity.

Skeletal muscle is the most abundant tissue of the healthy human body. It is also highly adaptable to different environmental stimuli, e.g. regular exercise. Exercise training improves overall health and muscle function, and can be used to prevent and treat several common diseases e.g. cardiovascular disease and type II diabetes. Therefore, it is of great importance to understand the molecular mechanisms behind adaptation processes in human skeletal muscle. In this study, we show that different expression variants from the same gene can be regulated in different directions with training, implicating alternative protein functions from one single gene. Such findings are emblematic of the complex mechanisms regulating the effects of training. We also find that training changes the activity of functionally unknown parts of the genome, with the potential for new proteins involved in the health-enhancing effects of exercise. Additionally, our results challenge the belief of a skeletal muscle memory, where previous training can affect the response to a subsequent training period. Overall, we provide understanding of the skeletal muscle biology and novel insights into the mechanisms behind the massive benefits of regular exercise on the human skeletal muscle transcriptome, inspiring further studies for deeper investigation.

Mitochondrial Polyadenylation Is a One-Step Process Required for mRNA Integrity and tRNA Maturation

Polyadenylation has well characterised roles in RNA turnover and translation in a variety of biological systems. While polyadenylation on mitochondrial transcripts has been suggested to be a two-step process required to complete translational stop codons, its involvement in mitochondrial RNA turnover is less well understood. We studied knockdown and knockout models of the mitochondrial poly(A) polymerase (MTPAP) in Drosophila melanogaster and demonstrate that polyadenylation of mitochondrial mRNAs is exclusively performed by MTPAP. Further, our results show that mitochondrial polyadenylation does not regulate mRNA stability but protects the 3' terminal integrity, and that despite a lack of functioning 3' ends, these trimmed transcripts are translated, suggesting that polyadenylation is not required for mitochondrial translation. Additionally, loss of MTPAP leads to reduced steady-state levels and disturbed maturation of tRNA^Cys, indicating that polyadenylation in mitochondria might be important for the stability and maturation of specific tRNAs.

The polyadenylation of cellular RNAs is a well-studied signal for gene expression, with a defined function in either RNA turnover or translation, in the majority of systems. In mammalian mitochondria the role of polyadenylation is less clear, and can to date only be attributed to completing the translational stop signal on several mitochondrial transcripts. Previous work though demonstrated that mitochondrial polyadenylation requires a certain length and shortening of the poly(A) tail signal has detrimental effects on mitochondrial function. In this study we deleted the mitochondrial polymerase responsible for polyadenylation in the fly, Drosophila melanogaster, and demonstrate that the mitochondrial poly(A) tail is essential for preserving the 3’ ends of mitochondrial transcripts, with no other enzyme capable of completing stop signals. Our study also shows that polyadenylation does not regulate transcript stability nor is it required for translation, but might be involved in the maturation of certain mitochondrial tRNAs. We therefore conclude that besides completing translational stop signals, mitochondrial polyadenylation protects the 3’ termini from degradation.

Modulation of Hemostatic and Inflammatory Responses by Leptospira Spp

Leptospirosis is a worldwide spread zoonotic and neglected infectious disease of human and veterinary concern that is caused by pathogenic Leptospira species. In severe infections, hemostatic impairments such as coagulation/fibrinolysis dysfunction are frequently observed. These complications often occur when the host response is controlled and/or modulated by the bacterial pathogen. In the present investigation, we aimed to analyze the modulation of the hemostatic and inflammatory host responses by the bacterial pathogen Leptospira. The effects of leptospires and their secreted products on stimulation of human intrinsic and extrinsic pathways of coagulation were investigated by means of altered clotting times, assembly and activation of contact system and induction of tissue factor. We show that both extrinsic and intrinsic coagulation cascades are modulated in response to Leptospira or leptospiral secreted proteins. We further find that the pro-inflammatory mediator bradykinin is released following contact activation at the bacterial surface and that pro-coagulant microvesicles are shed from monocytes in response to infection. Also, we show that human leptospirosis patients present higher levels of circulating pro-coagulant microvesicles than healthy individuals. Here we show that both pathways of the coagulation system are modulated by leptospires, possibly leading to altered hemostatic and inflammatory responses during the disease. Our results contribute to the understanding of the leptospirosis pathophysiological mechanisms and may open new routes for the discovery of novel treatments for the severe manifestations of the disease.

Leptospirosis is one of the most relevant and spread zoonotic and neglected infectious diseases affecting humans and other mammals, and is caused by pathogenic bacteria of the genus Leptospira. During infectious diseases, when bacterial pathogens control and/or modulate the host response, impaired hemostasis and inflammation are frequently observed. Here we studied the effects of leptospires and their secreted products on stimulation of human intrinsic and extrinsic pathways of coagulation, showing that both coagulation cascades are modulated in response to Leptospira or leptospiral secreted proteins. We further find that activation of the coagulation cascades culminates in the release of the pro-inflammatory mediator bradykinin and noted an induction of pro-coagulant microvesicles. These findings contribute to a better understanding of the local and systemic hemostastic complications during leptospirosis. Collectively, our results show how leptospires can affect host responses, possibly leading to altered host responses during the disease and giving rise to the leptospirosis symptomatology.

WUSCHEL-RELATED HOMEOBOX 2 is important for protoderm and suspensor development in the gymnosperm Norway spruce

BACKGROUND

Distinct expression domains of WUSCHEL-RELATED HOMEOBOX (WOX) gene family members are involved in patterning and morphogenesis of the early embryo in Arabidopsis. However, the role of WOX genes in other taxa, including gymnosperms, remains elusive. Here, we use somatic embryos and reverse genetics for studying expression and function of PaWOX2, the corresponding homolog of AtWOX2 in the gymnosperm Picea abies (Pa; Norway spruce).

RESULTS

The mRNA level of PaWOX2 was transiently up-regulated during early and late embryogeny. PaWOX2 mRNA in early and early late embryos was detected both in the embryonal mass and in the upper part of the suspensor. Down-regulation of PaWOX2 during development of early embryos resulted in aberrant early embryos, which failed to form a proper protoderm. Cells on the surface layer of the embryonal mass became vacuolated, and new embryogenic tissue differentiated from the embryonal mass. In addition, the aberrant early embryos lacked a distinct border between the embryonal mass, and the suspensor and the length of the suspensor cells was reduced. Down-regulation of PaWOX2 in the beginning of embryo development, before late embryos were formed, caused a significant decrease in the yield of mature embryos. On the contrary, down-regulation of PaWOX2 after late embryos were formed had no effect on further embryo development and maturation.

CONCLUSIONS

Our data suggest an evolutionarily conserved function of WOX2 in protoderm formation early during embryo development among seed plants. In addition, PaWOX2 might exert a unique function in suspensor expansion in gymnosperms.

Electronic plants

The roots, stems, leaves, and vascular circuitry of higher plants are responsible for conveying the chemical signals that regulate growth and functions. From a certain perspective, these features are analogous to the contacts, interconnections, devices, and wires of discrete and integrated electronic circuits. Although many attempts have been made to augment plant function with electroactive materials, plants' "circuitry" has never been directly merged with electronics. We report analog and digital organic electronic circuits and devices manufactured in living plants. The four key components of a circuit have been achieved using the xylem, leaves, veins, and signals of the plant as the template and integral part of the circuit elements and functions. With integrated and distributed electronics in plants, one can envisage a range of applications including precision recording and regulation of physiology, energy harvesting from photosynthesis, and alternatives to genetic modification for plant optimization.