Correlation among traits is a fundamental feature of biological systems that remains difficult to study. To address this problem, we developed a flexible approach that allows us to identify factors associated with inter-individual variation in correlation. We use data from three human cohorts to study the effects of genetic and environmental variation on correlations among mRNA transcripts and among NMR metabolites. We first show that environmental exposures (infection and disease) lead to a systematic loss of correlation, which we define as 'decoherence'. Using longitudinal data, we show that decoherent metabolites are better predictors of whether someone will develop metabolic syndrome than metabolites commonly used as biomarkers of this disease. Finally, we demonstrate that correlation itself is under genetic control by mapping hundreds of 'correlation quantitative trait loci (QTLs)'. Together, this work furthers our understanding of how and why coordinated biological processes break down, and points to a potential role for decoherence in disease.

Deep brain stimulation (DBS) electrode implant trajectories are stereotactically defined using preoperative neuroimaging. To validate the correct trajectory, microelectrode recordings (MERs) or local field potential recordings can be used to extend neuroanatomical information (defined by MRI) with neurophysiological activity patterns recorded from micro- and macroelectrodes probing the surgical target site. Currently, these two sources of information (imaging vs. electrophysiology) are analyzed separately, while means to fuse both data streams have not been introduced.

Imprinted gene expression occurs during seed development in plants and is associated with differential DNA methylation of parental alleles, particularly at proximal transposable elements (TEs). Imprinting variability could contribute to observed parent-of-origin effects on seed development. We investigated intraspecific variation in imprinting, coupled with analysis of DNA methylation and small RNAs, among three Arabidopsis strains with diverse seed phenotypes. The majority of imprinted genes were parentally biased in the same manner among all strains. However, we identified several examples of allele-specific imprinting correlated with intraspecific epigenetic variation at a TE. We successfully predicted imprinting in additional strains based on methylation variability. We conclude that there is standing variation in imprinting even in recently diverged genotypes due to intraspecific epiallelic variation. Our data demonstrate that epiallelic variation and genomic imprinting intersect to produce novel gene expression patterns in seeds.

Humans usually assess things not according to their absolute value, but relative to reference points - a main tenant of Prospect Theory. For example, people rate a new salary relative to previous salaries and salaries of their peers, rather than absolute income. We demonstrate a similar effect in an insect: ants expecting to find low-quality food showed higher acceptance of medium-quality food than ants expecting medium quality, and vice versa for high expectations. Further experiments demonstrate that these contrast effects arise from cognitive rather than mere sensory or pre-cognitive perceptual causes. Social information gained inside the nest can also serve as a reference point: the quality of food received from other ants affected the perceived value of food found later. Value judgement is a key element in decision making, and thus relative value perception strongly influences which option is chosen and ultimately how all animals make decisions.

Impaired spermatogenesis and male infertility are common manifestations associated with mitochondrial diseases, yet the underlying mechanisms linking these conditions remain elusive. In this study, we demonstrate that mice deficient for the mitochondrial intra-membrane rhomboid protease PARL, a recently reported model of the mitochondrial encephalopathy Leigh syndrome, develop early testicular atrophy caused by a complete arrest of spermatogenesis during meiotic prophase I, followed by degeneration and death of arrested spermatocytes. This process is independent of neurodegeneration. Interestingly, genetic modifications of PINK1, PGAM5, and TTC19 - three major substrates of PARL with important roles in mitochondrial homeostasis - fail to reproduce or modify this severe phenotype, indicating that the spermatogenic arrest arises from distinct molecular pathways. We further observed severe abnormalities in mitochondrial ultrastructure in PARL-deficient spermatocytes, along with prominent electron transfer chain defects, disrupted coenzyme Q (CoQ) biosynthesis, and metabolic rewiring. These mitochondrial defects are associated with a germ cell-specific decrease in GPX4 expression leading arrested spermatocytes to ferroptosis - a regulated cell death modality characterized by uncontrolled lipid peroxidation. Our results suggest that mitochondrial defects induced by PARL depletion act as an initiating trigger for ferroptosis in primary spermatocytes through simultaneous effects on GPX4 and CoQ - two major inhibitors of ferroptosis. These findings shed new light on the potential role of ferroptosis in the pathogenesis of mitochondrial diseases and male infertility warranting further investigation.

Amyotrophic lateral sclerosis (ALS) is a progressive neurodegenerative disease involving cytotoxic conformations of Cu, Zn superoxide dismutase (SOD1). A major challenge in understanding ALS disease pathology has been the identification and atomic-level characterization of these conformers. Here, we use a combination of NMR methods to detect four distinct sparsely populated and transiently formed thermally accessible conformers in equilibrium with the native state of immature SOD1 (apoSOD1(2SH)). Structural models of two of these establish that they possess features present in the mature dimeric protein. In contrast, the other two are non-native oligomers in which the native dimer interface and the electrostatic loop mediate the formation of aberrant intermolecular interactions. Our results show that apoSOD1(2SH) has a rugged free energy landscape that codes for distinct kinetic pathways leading to either maturation or non-native association and provide a starting point for a detailed atomic-level understanding of the mechanisms of SOD1 oligomerization.

Translation-dependent quality control pathways such as no-go decay (NGD), non-stop decay (NSD), and nonsense-mediated decay (NMD) govern protein synthesis and proteostasis by resolving non-translating ribosomes and preventing the production of potentially toxic peptides derived from faulty and aberrant mRNAs. However, how translation is altered and the in vivo defects that arise in the absence of these pathways are poorly understood. Here, we show that the NGD/NSD factors Pelo and Hbs1l are critical in mice for cerebellar neurogenesis but expendable for survival of these neurons after development. Analysis of mutant mouse embryonic fibroblasts revealed translational pauses, alteration of signaling pathways, and translational reprogramming. Similar effects on signaling pathways, including mTOR activation, the translatome and mouse cerebellar development were observed upon deletion of the NMD factor Upf2. Our data reveal that these quality control pathways that function to mitigate errors at distinct steps in translation can evoke similar cellular responses.

A20 haploinsufficiency (HA20) is an autoinflammatory disease caused by heterozygous loss-of-function variations in TNFAIP3, the gene encoding the A20 protein. Diagnosis of HA20 is challenging due to its heterogeneous clinical presentation and the lack of pathognomonic symptoms. While the pathogenic effect of TNFAIP3 truncating variations is clearly established, that of missense variations is difficult to determine. Herein, we identified a novel TNFAIP3 variation, p.(Leu236Pro), located in the A20 ovarian tumor (OTU) domain and demonstrated its pathogenicity. In the patients' primary cells, we observed reduced A20 levels. Protein destabilization was predicted in silico for A20_Leu236Pro and enhanced proteasomal degradation was confirmed in vitro through a flow cytometry-based functional assay. By applying this approach to the study of another missense variant, A20_Leu275Pro, for which no functional characterization has been performed to date, we showed that this variant also undergoes enhanced proteasomal degradation. Moreover, we showed a disrupted ability of A20_Leu236Pro to inhibit the NF-κB pathway and to deubiquitinate its substrate TRAF6. Structural modeling revealed that two residues involved in OTU pathogenic missense variations (i.e. Glu192Lys and Cys243Tyr) establish common interactions with Leu236. Interpretation of newly identified missense variations is challenging, requiring, as illustrated here, functional demonstration of their pathogenicity. Together with functional studies, in silico structure analysis is a valuable approach that allowed us (i) to provide a mechanistic explanation for the haploinsufficiency resulting from missense variations and (ii) to unveil a region within the OTU domain critical for A20 function.

Shprintzen-Goldberg syndrome (SGS) is a multisystemic connective tissue disorder, with considerable clinical overlap with Marfan and Loeys-Dietz syndromes. These syndromes have commonly been associated with enhanced TGF-β signaling. In SGS patients, heterozygous point mutations have been mapped to the transcriptional co-repressor SKI, which is a negative regulator of TGF-β signaling that is rapidly degraded upon ligand stimulation. The molecular consequences of these mutations, however, are not understood. Here we use a combination of structural biology, genome editing, and biochemistry to show that SGS mutations in SKI abolish its binding to phosphorylated SMAD2 and SMAD3. This results in stabilization of SKI and consequently attenuation of TGF-β responses, both in knockin cells expressing an SGS mutation and in fibroblasts from SGS patients. Thus, we reveal that SGS is associated with an attenuation of TGF-β-induced transcriptional responses, and not enhancement, which has important implications for other Marfan-related syndromes.

In the early days of HIV treatment, drug resistance occurred rapidly and predictably in all patients, but under modern treatments, resistance arises slowly, if at all. The probability of resistance should be controlled by the rate of generation of resistance mutations. If many adaptive mutations arise simultaneously, then adaptation proceeds by soft selective sweeps in which multiple adaptive mutations spread concomitantly, but if adaptive mutations occur rarely in the population, then a single adaptive mutation should spread alone in a hard selective sweep. Here, we use 6717 HIV-1 consensus sequences from patients treated with first-line therapies between 1989 and 2013 to confirm that the transition from fast to slow evolution of drug resistance was indeed accompanied with the expected transition from soft to hard selective sweeps. This suggests more generally that evolution proceeds via hard sweeps if resistance is unlikely and via soft sweeps if it is likely.

Spinal Muscular Atrophy (SMA) is caused by diminished Survival of Motor Neuron (SMN) protein, leading to neuromuscular junction (NMJ) dysfunction and spinal motor neuron (MN) loss. Here, we report that reduced SMN function impacts the action of a pertinent microRNA and its mRNA target in MNs. Loss of the C. elegans SMN ortholog, SMN-1, causes NMJ defects. We found that increased levels of the C. elegans Gemin3 ortholog, MEL-46, ameliorates these defects. Increased MEL-46 levels also restored perturbed microRNA (miR-2) function in smn-1(lf) animals. We determined that miR-2 regulates expression of the C. elegans M2 muscarinic receptor (m2R) ortholog, GAR-2. GAR-2 loss ameliorated smn-1(lf) and mel-46(lf) synaptic defects. In an SMA mouse model, m2R levels were increased and pharmacological inhibition of m2R rescued MN process defects. Collectively, these results suggest decreased SMN leads to defective microRNA function via MEL-46 misregulation, followed by increased m2R expression, and neuronal dysfunction in SMA.

Mutations in pre-synaptic voltage-gated calcium channels can lead to familial hemiplegic migraine type 1 (FHM1). While mammalian studies indicate that the migraine brain is hyperexcitable due to enhanced excitation or reduced inhibition, the molecular and cellular mechanisms underlying this excitatory/inhibitory (E/I) imbalance are poorly understood. We identified a gain-of-function (gf) mutation in the Caenorhabditis elegans CaV2 channel α1 subunit, UNC-2, which leads to increased calcium currents. unc-2(zf35gf) mutants exhibit hyperactivity and seizure-like motor behaviors. Expression of the unc-2 gene with FHM1 substitutions R192Q and S218L leads to hyperactivity similar to that of unc-2(zf35gf) mutants. unc-2(zf35gf) mutants display increased cholinergic and decreased GABAergic transmission. Moreover, increased cholinergic transmission in unc-2(zf35gf) mutants leads to an increase of cholinergic synapses and a TAX-6/calcineurin-dependent reduction of GABA synapses. Our studies reveal mechanisms through which CaV2 gain-of-function mutations disrupt excitation-inhibition balance in the nervous system.

Dopamine and striatal dysfunctions play a key role in the pathophysiology of Parkinson's disease (PD) and Dystonia, but our understanding of the changes in the discharge rate and pattern of striatal projection neurons (SPNs) remains limited. Here, we recorded and examined multi-unit signals from the striatum of PD and dystonic patients undergoing deep brain stimulation surgeries. Contrary to earlier human findings, we found no drastic changes in the spontaneous discharge of the well-isolated and stationary SPNs of the PD patients compared to the dystonic patients or to the normal levels of striatal activity reported in healthy animals. Moreover, cluster analysis using SPN discharge properties did not characterize two well-separated SPN subpopulations, indicating no SPN subpopulation-specific (D1 or D2 SPNs) discharge alterations in the pathological state. Our results imply that small to moderate changes in spontaneous SPN discharge related to PD and Dystonia are likely amplified by basal ganglia downstream structures.

Every decision that we make involves a conflict between exploiting our current knowledge of an action's value or exploring alternative courses of action that might lead to a better, or worse outcome. The sub-cortical nuclei that make up the basal ganglia have been proposed as a neural circuit that may contribute to resolving this explore-exploit 'dilemma'. To test this hypothesis, we examined the effects of neuromodulating the basal ganglia's output nucleus, the globus pallidus interna, in patients who had undergone deep brain stimulation (DBS) for isolated dystonia. Neuromodulation enhanced the number of exploratory choices to the lower value option in a two-armed bandit probabilistic reversal-learning task. Enhanced exploration was explained by a reduction in the rate of evidence accumulation (drift rate) in a reinforcement learning drift diffusion model. We estimated the functional connectivity profile between the stimulating DBS electrode and the rest of the brain using a normative functional connectome derived from heathy controls. Variation in the extent of neuromodulation induced exploration between patients was associated with functional connectivity from the stimulation electrode site to a distributed brain functional network. We conclude that the basal ganglia's output nucleus, the globus pallidus interna, can adaptively modify decision choice when faced with the dilemma to explore or exploit.

Retinitis pigmentosa (RP) and associated inherited retinal diseases (IRDs) are caused by rod photoreceptor degeneration, necessitating therapeutics promoting rod photoreceptor survival. To address this, we tested compounds for neuroprotective effects in multiple zebrafish and mouse RP models, reasoning drugs effective across species and/or independent of disease mutation may translate better clinically. We first performed a large-scale phenotypic drug screen for compounds promoting rod cell survival in a larval zebrafish model of inducible RP. We tested 2934 compounds, mostly human-approved drugs, across six concentrations, resulting in 113 compounds being identified as hits. Secondary tests of 42 high-priority hits confirmed eleven lead candidates. Leads were then evaluated in a series of mouse RP models in an effort to identify compounds effective across species and RP models, that is, potential pan-disease therapeutics. Nine of 11 leads exhibited neuroprotective effects in mouse primary photoreceptor cultures, and three promoted photoreceptor survival in mouse rd1 retinal explants. Both shared and complementary mechanisms of action were implicated across leads. Shared target tests implicated parp1-dependent cell death in our zebrafish RP model. Complementation tests revealed enhanced and additive/synergistic neuroprotective effects of paired drug combinations in mouse photoreceptor cultures and zebrafish, respectively. These results highlight the value of cross-species/multi-model phenotypic drug discovery and suggest combinatorial drug therapies may provide enhanced therapeutic benefits for RP patients.

Vascular networks surrounding individual organs are important for their development, maintenance, and function; however, how these networks are assembled remains poorly understood. Here we show that CNS progenitors, referred to as radial glia, modulate vascular patterning around the spinal cord by acting as negative regulators. We found that radial glia ablation in zebrafish embryos leads to excessive sprouting of the trunk vessels around the spinal cord, and exclusively those of venous identity. Mechanistically, we determined that radial glia control this process via the Vegf decoy receptor sFlt1: sflt1 mutants exhibit the venous over-sprouting observed in radial glia-ablated larvae, and sFlt1 overexpression rescues it. Genetic mosaic analyses show that sFlt1 function in trunk endothelial cells can limit their over-sprouting. Together, our findings identify CNS-resident progenitors as critical angiogenic regulators that determine the precise patterning of the vasculature around the spinal cord, providing novel insights into vascular network formation around developing organs.

Cyanobacteria are photosynthetic bacteria with a unique CO2 concentrating mechanism (CCM), enhancing carbon fixation. Understanding the CCM requires a systems level perspective of how molecular components work together to enhance CO2 fixation. We present a mathematical model of the cyanobacterial CCM, giving the parameter regime (expression levels, catalytic rates, permeability of carboxysome shell) for efficient carbon fixation. Efficiency requires saturating the RuBisCO reaction, staying below saturation for carbonic anhydrase, and avoiding wasteful oxygenation reactions. We find selectivity at the carboxysome shell is not necessary; there is an optimal non-specific carboxysome shell permeability. We compare the efficacy of facilitated CO2 uptake, CO2 scavenging, and HCO3- transport with varying external pH. At the optimal carboxysome permeability, contributions from CO2 scavenging at the cell membrane are small. We examine the cumulative benefits of CCM spatial organization strategies: enzyme co-localization and compartmentalization.

Distinct anatomical and spectral channels are thought to play specialized roles in the communication within cortical networks. While activity in the alpha and beta frequency range (7 - 40 Hz) is thought to predominantly originate from infragranular cortical layers conveying feedback-related information, activity in the gamma range (>40 Hz) dominates in supragranular layers communicating feedforward signals. We leveraged high precision MEG to test this proposal, directly and non-invasively, in human participants performing visually cued actions. We found that visual alpha mapped onto deep cortical laminae, whereas visual gamma predominantly occurred more superficially. This lamina-specificity was echoed in movement-related sensorimotor beta and gamma activity. These lamina-specific pre- and post- movement changes in sensorimotor beta and gamma activity suggest a more complex functional role than the proposed feedback and feedforward communication in sensory cortex. Distinct frequency channels thus operate in a lamina-specific manner across cortex, but may fulfill distinct functional roles in sensory and motor processes.

Ossification of the posterior longitudinal ligament of the spine (OPLL) is an intractable disease leading to severe neurological deficits. Its etiology and pathogenesis are primarily unknown. The relationship between OPLL and comorbidities, especially type 2 diabetes (T2D) and high body mass index (BMI), has been the focus of attention; however, no trait has been proven to have a causal relationship. We conducted a meta-analysis of genome-wide association studies (GWASs) using 22,016 Japanese individuals and identified 14 significant loci, 8 of which were previously unreported. We then conducted a gene-based association analysis and a transcriptome-wide Mendelian randomization approach and identified three candidate genes for each. Partitioning heritability enrichment analyses observed significant enrichment of the polygenic signals in the active enhancers of the connective/bone cell group, especially H3K27ac in chondrogenic differentiation cells, as well as the immune/hematopoietic cell group. Single-cell RNA sequencing of Achilles tendon cells from a mouse Achilles tendon ossification model confirmed the expression of genes in GWAS and post-GWAS analyses in mesenchymal and immune cells. Genetic correlations with 96 complex traits showed positive correlations with T2D and BMI and a negative correlation with cerebral aneurysm. Mendelian randomization analysis demonstrated a significant causal effect of increased BMI and high bone mineral density on OPLL. We evaluated the clinical images in detail and classified OPLL into cervical, thoracic, and the other types. GWAS subanalyses identified subtype-specific signals. A polygenic risk score for BMI demonstrated that the effect of BMI was particularly strong in thoracic OPLL. Our study provides genetic insight into the etiology and pathogenesis of OPLL and is expected to serve as a basis for future treatment development.

Biological age, distinct from an individual's chronological age, has been studied extensively through predictive aging clocks. However, these clocks have limited accuracy in short time-scales. Here we trained deep learning models on fundus images from the EyePACS dataset to predict individuals' chronological age. Our retinal aging clocking, 'eyeAge', predicted chronological age more accurately than other aging clocks (mean absolute error of 2.86 and 3.30 years on quality-filtered data from EyePACS and UK Biobank, respectively). Additionally, eyeAge was independent of blood marker-based measures of biological age, maintaining an all-cause mortality hazard ratio of 1.026 even when adjusted for phenotypic age. The individual-specific nature of eyeAge was reinforced via multiple GWAS hits in the UK Biobank cohort. The top GWAS locus was further validated via knockdown of the fly homolog, Alk, which slowed age-related decline in vision in flies. This study demonstrates the potential utility of a retinal aging clock for studying aging and age-related diseases and quantitatively measuring aging on very short time-scales, opening avenues for quick and actionable evaluation of gero-protective therapeutics.