Germline variants in tumor suppressor genes (TSGs) can result in RNA mis-splicing and predisposition to cancer. However, identification of variants that impact splicing remains a challenge, ...contributing to a substantial proportion of patients with suspected hereditary cancer syndromes remaining without a molecular diagnosis. To address this, we used capture RNA-sequencing (RNA-seq) to generate a splicing profile of 18 TSGs (
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
, and
) in 345 whole-blood samples from healthy donors. We subsequently demonstrated that this approach can detect mis-splicing by comparing splicing profiles from the control dataset to profiles generated from whole blood of individuals previously identified with pathogenic germline splicing variants in these genes. To assess the utility of our TSG splicing profile to prospectively identify pathogenic splicing variants, we performed concurrent capture DNA and RNA-seq in a cohort of 1000 patients with suspected hereditary cancer syndromes. This approach improved the diagnostic yield in this cohort, resulting in a 9.1% relative increase in the detection of pathogenic variants, demonstrating the utility of performing simultaneous DNA and RNA genetic testing in a clinical context.
Background
Transforming growth factor-βs (TGF-βs) are a group of multifunctional proteins that have neuroprotective roles in various experimental models. We previously reported that intrathecal ...(i.t.) injections of TGF-β1 significantly inhibit neuropathy-induced thermal hyperalgesia, spinal microglia and astrocyte activation, as well as upregulation of tumor necrosis factor-α. However, additional cellular mechanisms for the antinociceptive effects of TGF-β1, such as the mitogen-activated protein kinase (MAPK) pathway, have not been elucidated. During persistent pain, activation of MAPKs, especially p38 and extracellular signal-regulated kinase (ERK), have crucial roles in the induction and maintenance of pain hypersensitivity, via both nontranscriptional and transcriptional regulation. In the present study, we used a chronic constriction injury (CCI) rat model to explore the role of spinal p38 and ERK in the analgesic effects of TGF-β1.
Methods
We investigated the cellular mechanisms of the antinociceptive effects of i.t. injections of TGF-β1 in CCI induced neuropathic rats by spinal immunohistofluorescence analyses.
Results
The results demonstrated that the antinociceptive effects of TGF-β1 (5 ng) were maintained at greater than 50 % of the maximum possible effect in rats with CCI for at least 6 h after a single i.t. administration. Thus, we further examined these alterations in spinal p38 and ERK from 0.5 to 6 h after i.t. administration of TGF-β1. TGF-β1 significantly attenuated CCI-induced upregulation of phosphorylated p38 (phospho-p38) and phosphorylated ERK (phospho-ERK) expression in the dorsal horn of the lumbar spinal cord. Double immunofluorescence staining illustrated that upregulation of spinal phospho-p38 was localized to neurons, activated microglial cells, and activated astrocytes in rats with CCI. Additionally, increased phospho-ERK occurred in activated microglial cells and activated astrocytes. Furthermore, i.t. administration of TGF-β1 markedly inhibited phospho-p38 upregulation in neurons, microglial cells, and astrocytes. However, i.t. injection of TGF-β1 also reduced phospho-ERK upregulation in microglial cells and astrocytes.
Conclusions
The present results demonstrate that suppressing p38 and ERK activity affects TGF-β1-induced analgesia during neuropathy.
Characterizing the conformations of protein in the transition state ensemble (TSE) is important for studying protein folding. A promising approach pioneered by Vendruscolo et al. Nature (London) 409, ...641 (2001) to study TSE is to generate conformations that satisfy all constraints imposed by the experimentally measured φ values that provide information about the native likeness of the transition states. Faísca et al. J. Chem. Phys. 129, 095108 (2008) generated conformations of TSE based on the criterion that, starting from a TS conformation, the probabilities of folding and unfolding are about equal through Markov Chain Monte Carlo (MCMC) simulations. In this study, we use the technique of constrained sequential Monte Carlo method Lin et al., J. Chem. Phys. 129, 094101 (2008); Zhang et al. Proteins 66, 61 (2007) to generate TSE conformations of acylphosphatase of 98 residues that satisfy the φ-value constraints, as well as the criterion that each conformation has a folding probability of 0.5 by Monte Carlo simulations. We adopt a two stage process and first generate 5000 contact maps satisfying the φ-value constraints. Each contact map is then used to generate 1000 properly weighted conformations. After clustering similar conformations, we obtain a set of properly weighted samples of 4185 candidate clusters. Representative conformation of each of these cluster is then selected and 50 runs of Markov chain Monte Carlo (MCMC) simulation are carried using a regrowth move set. We then select a subset of 1501 conformations that have equal probabilities to fold and to unfold as the set of TSE. These 1501 samples characterize well the distribution of transition state ensemble conformations of acylphosphatase. Compared with previous studies, our approach can access much wider conformational space and can objectively generate conformations that satisfy the φ-value constraints and the criterion of 0.5 folding probability without bias. In contrast to previous studies, our results show that transition state conformations are very diverse and are far from nativelike when measured in cartesian root-mean-square deviation (cRMSD): the average cRMSD between TSE conformations and the native structure is 9.4 Å for this short protein, instead of 6 Å reported in previous studies. In addition, we found that the average fraction of native contacts in the TSE is 0.37, with enrichment in native-like β-sheets and a shortage of long range contacts, suggesting such contacts form at a later stage of folding. We further calculate the first passage time of folding of TSE conformations through calculation of physical time associated with the regrowth moves in MCMC simulation through mapping such moves to a Markovian state model, whose transition time was obtained by Langevin dynamics simulations. Our results indicate that despite the large structural diversity of the TSE, they are characterized by similar folding time. Our approach is general and can be used to study TSE in other macromolecules.
The liver CYP1A2 enzyme may metabolize antidepressant escitalopram (S-CIT) to S-desmethylcitalopram (S-DCIT) and S-didesmethylcitalopram (S-DDCIT). This study tested whether genetic polymorphisms in ...the CYP1A2 gene are associated with the treatment responses to S-CIT.
Ten SNPs in CYP1A2 were selected and genotyped in 158 patients under S-CIT treatment. The serum levels of S-CIT and its metabolites were measured by HPLC.
CYP1A2 SNPs rs2069521, rs2069526, rs4646425 and rs4646427 are significantly associated with the metabolic ratios of S-DDCIT/S-DCIT (p = 0.002, 0.018, 0.008 and 0.004, respectively) at week 2 of treatment. Carriers of the allele types associated with higher S-DDCIT/S-DCIT ratios had more severe side effects.
These results suggest that genetic variants in CYP1A2 may be indicators for S-CIT metabolism and that the fast metabolizers may experience more severe adverse reactions in the early stages of S-CIT treatment. Original submitted 27 December 2012; Revision submitted 15 May 2013.
Inferring three-dimensional structural information of biomacromolecules such as proteins from limited experimental data is an important and challenging task. Nuclear Overhauser effect measurements ...based on nucleic magnetic resonance, disulfide linking, and electron paramagnetic resonance labeling studies can all provide useful partial distance constraint characteristic of the conformations of proteins. In this study, we describe a general approach for reconstructing conformations of biomolecules that are consistent with given distance constraints. Such constraints can be in the form of upper bounds and lower bounds of distances between residue pairs, contact maps based on specific contact distance cutoff values, or indirect distance constraints such as experimental phi-value measurement. Our approach is based on the framework of sequential Monte Carlo method, a chain growth-based method. We have developed a novel growth potential function to guide the generation of conformations that satisfy given distance constraints. This potential function incorporates not only the distance information of current residue during growth but also the distance information of future residues by introducing global distance upper bounds between residue pairs and the placement of reference points. To obtain protein conformations from indirect distance constraints in the form of experimental phi-values, we first generate properly weighted contact maps satisfying phi-value constraints, we then generate conformations from these contact maps. We show that our approach can faithfully generate conformations that satisfy the given constraints, which approach the native structures when distance constraints for all residue pairs are given.
Paroxetine is a drug of choice in the treatment of major depressive disorder (MDD). Its metabolism has recently been reported to be mediated through the CYP enzymes 1A2 and 2D6. In our current study, ...we tested whether genetic polymorphisms in CYP1A2 are associated with the treatment efficacy and side effects of paroxetine.
A total of 241 MDD patients who had taken paroxetine continually for 8 weeks were recruited, and their steady state paroxetine concentrations were measured at weeks 2, 4 and 8. The genotypes of these patients were then assessed for the presence of nine SNPs, which were selected from either the HapMap Chinese ethnic group, the literature report or through their functional role in the CYP1A2 gene.
The allele types for SNPs rs4646425 (permutation p = 0.03), rs2472304 (permutation p = 0.01) and rs2470890 (permutation p = 0.004) demonstrated significant associations with paroxetine treatment remission at week 8. Response rates in the Hamilton Rating Scale for Depression (HAM-D) and for The Hamilton Rating Scale for Anxiety (HAM-A) were significantly associated with the SNPs rs4646425 (p = 0.0126 and 0.0088 for HAM-D and HAM-A, respectively) and rs4646427 (p = 0.0067 and 0.0196 for HAM-D and HAM-A, respectively). The inducible SNP rs762551 had a significant association with paroxetine dose at week 4 (permutation p = 0.012). We did not find an association between these SNPs and the side effects or serum concentrations of paroxetine.
Genetic variants in the CYP1A2 region may be indicators of treatment response in MDD patients to paroxetine.
We study folding dynamics of proteinlike sequences on a square lattice using a physical move set that exhausts all possible conformational changes. By analytically solving the master equation, we ...follow the time-dependent probabilities of occupancy of all 802 075 conformations of 16 mers over 7 orders of time span. We find that (i) folding rates of these proteinlike sequences of the same length can differ by 4 orders of magnitude, (ii) folding rates of sequences of the same conformation can differ by a factor of 190, and (iii) parameters of the native structures, designability, and thermodynamic properties are weak predictors of the folding rates; rather, a basin analysis of the kinematic energy landscape defined by the moves can provide an excellent account of the observed folding rates.
Disorders of cobalamin deficiency are a heterogeneous group of disorders with at least 19 autosomal recessive-associated genes. Familial samples of an infant who died due to presumed cobalamin ...deficiency were referred for clinical exome sequencing. The patient died before obtaining a blood sample or skin biopsy, autopsy was declined, and DNA yielded from the newborn screening blood spot was insufficient for diagnostic testing. Whole-exome sequencing of the mother, father, and unaffected sister and tailored bioinformatics analysis was applied to search for mutations in underlying disorders with recessive inheritance. This approach identified alterations within two genes, each of which was carried by one parent. The mother carried a missense alteration in the MTR gene (c.3518C>T; p.P1173L) which was absent in the father and the sister. The father carried a translational frameshift alteration in the LMBRD1 gene (c.1056delG; p.L352Lfs*18) which was absent in the mother and present in the heterozygous state in the sister. These mutations in the MTR (MIM# 156570) and LMBRD1 (MIM# 612625) genes have been described in patients with disorders of cobalamin metabolism complementation groups cblG and cblF, respectively. The child’s clinical presentation and biochemical results demonstrated overlap with both cblG and cblF. Sanger sequencing using DNA from the infant’s blood spot confirmed the inheritance of the two alterations in compound heterozygous form. We present the first example of exome sequencing leading to a diagnosis in the absence of the affected patient. Furthermore, the data support the possibility for potential digenic inheritance associated with cobalamin deficiency.