Endoplasmic reticulum calcium homeostasis is critical for cellular functions and is disrupted in diverse pathologies including neurodegeneration and cardiovascular disease. Owing to the high concentration of calcium within the ER, studying this subcellular compartment requires tools that are optimized for these conditions. To develop a single-fluorophore genetically encoded calcium indicator for this organelle, we targeted a low affinity variant of GCaMP3 to the ER lumen (GCaMPer (10.19)). A set of viral vectors was constructed to express GCaMPer in human neuroblastoma cells, rat primary cortical neurons, and human induced pluripotent stem cell-derived cardiomyocytes. We observed dynamic changes in GCaMPer (10.19) fluorescence in response to pharmacologic manipulations of the ER calcium store. Additionally, periodic calcium efflux from the ER was observed during spontaneous beating of cardiomyocytes. GCaMPer (10.19) has utility in imaging ER calcium in living cells and providing insight into luminal calcium dynamics under physiologic and pathologic states.

Fluorescent protein-based sensors for detecting neuronal activity have been developed largely based on non-neuronal screening systems. However, the dynamics of neuronal state variables (e.g., voltage, calcium, etc.) are typically very rapid compared to those of non-excitable cells. We developed an electrical stimulation and fluorescence imaging platform based on dissociated rat primary neuronal cultures. We describe its use in testing genetically-encoded calcium indicators (GECIs). Efficient neuronal GECI expression was achieved using lentiviruses containing a neuronal-selective gene promoter. Action potentials (APs) and thus neuronal calcium levels were quantitatively controlled by electrical field stimulation, and fluorescence images were recorded. Images were segmented to extract fluorescence signals corresponding to individual GECI-expressing neurons, which improved sensitivity over full-field measurements. We demonstrate the superiority of screening GECIs in neurons compared with solution measurements. Neuronal screening was useful for efficient identification of variants with both improved response kinetics and high signal amplitudes. This platform can be used to screen many types of sensors with cellular resolution under realistic conditions where neuronal state variables are in relevant ranges with respect to timing and amplitude.

Although most proteins conform to the classical one-structure/one-function paradigm, an increasing number of proteins with dual structures and functions are emerging. These fold-switching proteins remodel their secondary structures in response to cellular stimuli, fostering multi-functionality and tight cellular control. Accurate predictions of fold-switching proteins could both suggest underlying mechanisms for uncharacterized biological processes and reveal potential drug targets. Previously, we developed a prediction method for fold-switching proteins based on secondary structure predictions and structure-based thermodynamic calculations. Given the large number of genomic sequences without homologous experimentally characterized structures, however, we sought to predict fold-switching proteins from their sequences alone. To do this, we leveraged state-of-the-art secondary structure predictions, which require only amino acid sequences but are not currently designed to identify structural duality in proteins. Thus, we hypothesized that incorrect and inconsistent secondary structure predictions could be good initial predictors of fold-switching proteins. We found that secondary structure predictions of fold-switching proteins with solved structures are indeed less accurate than secondary structure predictions of non-fold-switching proteins with solved structures. These inaccuracies result largely from the conformations of fold-switching proteins that are underrepresented in the Protein Data Bank (PDB), and, consequently, the training sets of secondary structure predictors. Given that secondary structure predictions are homology-based, we hypothesized that decontextualizing the inaccurately-predicted regions of fold-switching proteins could weaken the homology relationships between these regions and their overpopulated structural representatives. Thus, we reran secondary structure predictions on these regions in isolation and found that they were significantly more inconsistent than in regions of non-fold-switching proteins. Thus, inconsistent secondary structure predictions can serve as a preliminary marker of fold switching. These findings have implications for genomics and the future development of secondary structure predictors.

Extant fold-switching proteins remodel their secondary structures and change their functions in response to cellular stimuli, regulating biological processes and affecting human health. In spite of their biological importance, these proteins remain understudied. Few representative examples of fold switchers are available in the Protein Data Bank, and they are difficult to predict. In fact, all 96 experimentally validated examples of extant fold switchers were stumbled upon by chance. Thus, predictive methods are needed to expedite the process of discovering and characterizing more of these shapeshifting proteins. Previous approaches require a solved structure or all-atom simulations, greatly constraining their use. Here, we propose a high-throughput sequence-based method for predicting extant fold switchers that transition from α-helix in one conformation to β-strand in the other. This method leverages two previous observations: (1) α-helix <-> β-strand prediction discrepancies from JPred4 are a robust predictor of fold switching, and (2) the fold-switching regions (FSRs) of some extant fold switchers have different secondary structure propensities when expressed in isolation (isolated FSRs) than when expressed within the context of their parent protein (contextualized FSRs). Combining these two observations, we ran JPred4 on the sequences of isolated and contextualized FSRs from 14 known extant fold switchers and found α-helix <->β-strand prediction discrepancies in every case. To test the overall robustness of this finding, we randomly selected regions of proteins not expected to switch folds (single-fold proteins) and found significantly fewer α-helix <-> β-strand prediction discrepancies (p < 4.2*10−20, Kolmogorov-Smirnov test). Combining these discrepancies with the overall percentage of predicted secondary structure, we developed a classifier that often robustly identifies extant fold switchers (Matthews Correlation Coefficient of 0.70). Although this classifier had a high false negative rate (6/14), its false positive rate was very low (1/211), suggesting that it can be used to predict a subset of extant fold switchers from billions of available genomic sequences.

Temperature dependent sex determination (TSD) is the process by which the environmental temperature experienced during embryogenesis influences the sex of an organism, as in the red-eared slider turtle Trachemys scripta elegans. In accord with current paradigms of vertebrate sex determination, temperature is believed to exert its effects on sexual development in T. scripta entirely within the middle third of development, when the gonad is forming. However, whether temperature regulates the transcriptome in T. scripta early embryos in a manner that could influence secondary sex characteristics or establish a pro-male or pro-female environment has not been investigated. In addition, apart from a handful of candidate genes, very little is known about potential similarities between the expression cascade during TSD and the genetic cascade that drives mammalian sex determination. Here, we conducted an unbiased transcriptome-wide analysis of the effects of male- and female-promoting temperatures on the turtle embryo prior to gonad formation, and on the gonad during the temperature sensitive period. We found sexually dimorphic expression reflecting differences in steroidogenic enzymes and brain development prior to gonad formation. Within the gonad, we mapped a cascade of differential expression similar to the genetic cascade established in mammals. Using a Hidden Markov Model based clustering approach, we identified groups of genes that show heterochronic shifts between M. musculus and T. scripta. We propose a model in which multiple factors influenced by temperature accumulate during early gonadogenesis, and converge on the antagonistic regulation of aromatase to canalize sex determination near the end of the temperature sensitive window of development.

Down syndrome (DS) is a genetic disorder that causes cognitive impairment. The staggering effects associated with an extra copy of human chromosome 21 (HSA21) complicates mechanistic understanding of DS pathophysiology. We examined the neuron-astrocyte interplay in a fully recapitulated HSA21 trisomy cellular model differentiated from DS-patient-derived induced pluripotent stem cells (iPSCs). By combining calcium imaging with genetic approaches, we discovered the functional defects of DS astroglia and their effects on neuronal excitability. Compared with control isogenic astroglia, DS astroglia exhibited more-frequent spontaneous calcium fluctuations, which reduced the excitability of co-cultured neurons. Furthermore, suppressed neuronal activity could be rescued by abolishing astrocytic spontaneous calcium activity either chemically by blocking adenosine-mediated signaling or genetically by knockdown of inositol triphosphate (IP3) receptors or S100B, a calcium binding protein coded on HSA21. Our results suggest a mechanism by which DS alters the function of astrocytes, which subsequently disturbs neuronal excitability.

Cortical-feedback projections to primary sensory areas terminate most heavily in layer 1 (L1) of the neocortex, where they make synapses with tuft dendrites of pyramidal neurons. L1 input is thought to provide ‘contextual’ information, but the signals transmitted by L1 feedback remain uncharacterized. In the rodent somatosensory system, the spatially diffuse feedback projection from vibrissal motor cortex (vM1) to vibrissal somatosensory cortex (vS1, also known as the barrel cortex) may allow whisker touch to be interpreted in the context of whisker position to compute object location. When mice palpate objects with their whiskers to localize object features, whisker touch excites vS1 and later vM1 in a somatotopic manner. Here we use axonal calcium imaging to track activity in vM1–>vS1 afferents in L1 of the barrel cortex while mice performed whisker-dependent object localization. Spatially intermingled individual axons represent whisker movements, touch and other behavioural features. In a subpopulation of axons, activity depends on object location and persists for seconds after touch. Neurons in the barrel cortex thus have information to integrate movements and touches of multiple whiskers over time, key components of object identification and navigation by active touch.

Systemic lupus erythematosus (SLE) is an inflammatory autoimmune disease with a strong genetic component. African-Americans (AA) are at increased risk of SLE, but the genetic basis of this risk is largely unknown. To identify causal variants in SLE loci in AA, we performed admixture mapping followed by fine mapping in AA and European-Americans (EA). Through genome-wide admixture mapping in AA, we identified a strong SLE susceptibility locus at 2q22-24 (LOD=6.28), and the admixture signal is associated with the European ancestry (ancestry risk ratio  1.5). Large-scale genotypic analysis on 19,726 individuals of African and European ancestry revealed three independently associated variants in the IFIH1 gene: an intronic variant, rs13023380 [P(meta) = 5.20×10(-14); odds ratio, 95% confidence interval = 0.82 (0.78-0.87)], and two missense variants, rs1990760 (Ala946Thr) [P(meta) = 3.08×10(-7); 0.88 (0.84-0.93)] and rs10930046 (Arg460His) [P(dom) = 1.16×10(-8); 0.70 (0.62-0.79)]. Both missense variants produced dramatic phenotypic changes in apoptosis and inflammation-related gene expression. We experimentally validated function of the intronic SNP by DNA electrophoresis, protein identification, and in vitro protein binding assays. DNA carrying the intronic risk allele rs13023380 showed reduced binding efficiency to a cellular protein complex including nucleolin and lupus autoantigen Ku70/80, and showed reduced transcriptional activity in vivo. Thus, in SLE patients, genetic susceptibility could create a biochemical imbalance that dysregulates nucleolin, Ku70/80, or other nucleic acid regulatory proteins. This could promote antibody hypermutation and auto-antibody generation, further destabilizing the cellular network. Together with molecular modeling, our results establish a distinct role for IFIH1 in apoptosis, inflammation, and autoantibody production, and explain the molecular basis of these three risk alleles for SLE pathogenesis.

The calcium-modulated photoactivatable ratiometric integrator CaMPARI (Fosque et al., 2015) facilitates the study of neural circuits by permanently marking cells active during user-specified temporal windows. Permanent marking enables measurement of signals from large swathes of tissue and easy correlation of activity with other structural or functional labels. One potential application of CaMPARI is labeling neurons postsynaptic to specific populations targeted for optogenetic stimulation, giving rise to all-optical functional connectivity mapping. Here, we characterized the response of CaMPARI to several common types of neuronal calcium signals in mouse acute cortical brain slices. Our experiments show that CaMPARI is effectively converted by both action potentials and sub-threshold synaptic inputs, and that conversion level is correlated to synaptic strength. Importantly, we found that conversion rate can be tuned: it is linearly related to light intensity. At low photoconversion light levels CaMPARI offers a wide dynamic range due to slower conversion rate; at high light levels conversion is more rapid and more sensitive to activity. Finally, we employed CaMPARI and optogenetics for functional circuit mapping in ex vivo acute brain slices, which preserve in vivo-like connectivity of axon terminals. With a single light source, we stimulated channelrhodopsin-2-expressing long-range posteromedial (POm) thalamic axon terminals in cortex and induced CaMPARI conversion in recipient cortical neurons. We found that POm stimulation triggers robust photoconversion of layer 5 cortical neurons and weaker conversion of layer 2/3 neurons. Thus, CaMPARI enables network-wide, tunable, all-optical functional circuit mapping that captures supra- and sub-threshold depolarization. This article is protected by copyright. All rights reserved.

Recent reports have associated NCF2, encoding a core component of the multi-protein NADPH oxidase (NADPHO), with systemic lupus erythematosus (SLE) susceptibility in individuals of European ancestry. To identify ethnicity-specific and -robust variants within NCF2, we assessed 145 SNPs in and around the NCF2 gene in 5325 cases and 21 866 controls of European-American (EA), African-American (AA), Hispanic (HS) and Korean (KR) ancestry. Subsequent imputation, conditional, haplotype and bioinformatic analyses identified seven potentially functional SLE-predisposing variants. Association with non-synonymous rs17849502, previously reported in EA, was detected in EA, HS and AA (P(EA) = 1.01 × 10(-54), PHS = 3.68 × 10(-10), P(AA) = 0.03); synonymous rs17849501 was similarly significant. These SNPs were monomorphic in KR. Novel associations were detected with coding variants at rs35937854 in AA (PAA = 1.49 × 10(-9)), and rs13306575 in HS and KR (P(HS) = 7.04 × 10(-7), P(KR) = 3.30 × 10(-3)). In KR, a 3-SNP haplotype was significantly associated (P = 4.20 × 10(-7)), implying that SLE predisposing variants were tagged. Significant SNP-SNP interaction (P = 0.02) was detected between rs13306575 and rs17849502 in HS, and a dramatically increased risk (OR = 6.55) with a risk allele at each locus. Molecular modeling predicts that these non-synonymous mutations could disrupt NADPHO complex assembly. The risk allele of rs17849501, located in a conserved transcriptional regulatory region, increased reporter gene activity, suggesting in vivo enhancer function. Our results not only establish allelic heterogeneity within NCF2 associated with SLE, but also emphasize the utility of multi-ethnic cohorts to identify predisposing variants explaining additional phenotypic variance ('missing heritability') of complex diseases like SLE.