Whole-exome sequencing reveals overlap between macrophage activation syndrome in systemic juvenile idiopathic arthritis and familial hemophagocytic lymphohistiocytosis.

OBJECTIVE: Macrophage activation syndrome (MAS), a life-threatening complication of systemic juvenile idiopathic arthritis (JIA), resembles familial hemophagocytic lymphohistiocytosis (HLH), a constellation of autosomal-recessive immune disorders resulting from deficiency in cytolytic pathway proteins. We undertook this study to test our hypothesis that MAS predisposition in systemic JIA could be attributed to rare gene sequence variants affecting the cytotolytic pathway. METHODS: Whole-exome sequencing was used in 14 patients with systemic JIA and MAS and in their parents to identify protein-altering single-nucleotide polymorphisms/indels in known HLH-associated genes. To discover new candidate genes, the entire whole-exome sequencing data were filtered to identify protein-altering, rare recessive homozygous, compound heterozygous, and de novo variants with the potential to affect the cytolytic pathway. RESULTS: Heterozygous protein-altering rare variants in the known genes (LYST,MUNC13-4, and STXBP2) were found in 5 of 14 patients with systemic JIA and MAS (35.7%). This was in contrast to only 4 variants in 4 of 29 patients with systemic JIA without MAS (13.8%). Homozygosity and compound heterozygosity analysis applied to the entire whole-exome sequencing data in systemic JIA/MAS revealed 3 recessive pairs in 3 genes and compound heterozygotes in 73 genes. We also identified 20 heterozygous rare protein-altering variants that occurred in at least 2 patients. Many of the identified genes encoded proteins with a role in actin and microtubule reorganization and vesicle-mediated transport. "Cellular assembly and organization" was the top cellular function category based on Ingenuity Pathways Analysis (P < 3.10 × 10(-5) ). CONCLUSION: Whole-exome sequencing performed in patients with systemic JIA and MAS identified rare protein-altering variants in known HLH-associated genes as well as in new candidate genes.

Comprehensive gene expression profiling of rat lung reveals distinct acute and chronic responses to cigarette smoke inhalation.

Chronic obstructive pulmonary disease (COPD) is a smoking-related disease that lacks effective therapies due partly to the poor understanding of disease pathogenesis. The aim of this study was to identify molecular pathways that could be responsible for the damaging consequences of smoking. To do this, we employed Gene Set Enrichment Analysis to analyze differences in global gene expression, which we then related to the pathological changes induced by cigarette smoke (CS). Sprague-Dawley rats were exposed to whole body CS for 1 day and for various periods up to 8 mo. Gene Set Enrichment Analysis of microarray data identified that metabolic processes were most significantly increased early in the response to CS. Gene sets involved in stress response and inflammation were also upregulated. CS exposure increased neutrophil chemokines, cytokines, and proteases (MMP-12) linked to the pathogenesis of COPD. After a transient acute response, the CS-exposed rats developed a distinct molecular signature after 2 wk, which was followed by the chronic phase of the response. During this phase, gene sets related to immunity and defense progressively increased and predominated at the later time points in smoke-exposed rats. Chronic CS inhalation recapitulated many of the phenotypic changes observed in COPD patients including oxidative damage to macrophages, a slowly resolving inflammation, epithelial damage, mucus hypersecretion, airway fibrosis, and emphysema. As such, it appears that metabolic pathways are central to dealing with the stress of CS exposure; however, over time, inflammation and stress response gene sets become the most significantly affected in the chronic response to CS.

Estrogen-related receptor alpha is essential for the expression of antioxidant protection genes and mitochondrial function.

Estrogen-related receptor alpha (ERRalpha) is an important mediator of mitochondrial biogenesis and function. To investigate the transcriptional network controlling these phenomena, we investigated mitochondrial gene expression in embryonic fibroblasts isolated from ERRalpha null mice. Peroxisome proliferator-activated receptor gamma coactivator-1alpha (PGC-1alpha) stimulated mitochondrial gene expression program in control cells, but not in the ERRalpha null cells. Interestingly, the induction of levels of mitochondrial oxidative stress protection genes in response to increased PGC-1alpha levels was dependent on ERRalpha. Furthermore, we found that the PGC-1alpha-mediated induction of estrogen-related receptor gamma and nuclear respiratory factor 2 (NRF-2), was dependent on the presence of ERRalpha. Basal levels of NRF-2 were decreased in the absence of ERRalpha. The absence of ERRalpha resulted in a decrease in citrate synthase enzyme activity in response to PGC-1alpha overexpression. Our results indicate an essential role for ERRalpha as a key regulator of oxidative metabolism.

Extending the pathway analysis framework with a test for transcriptional variance implicates novel pathway modulation during myogenic differentiation.

MOTIVATION: We describe an extension of the pathway-based enrichment approach for analyzing microarray data via a robust test for transcriptional variance. The use of a variance test is intended to identify additional patterns of transcriptional regulation in which many genes in a pathway are up- and down-regulated. Such patterns may be indicative of the reciprocal regulation of pathway activators and inhibitors or of the differential regulation of separate biological sub-processes and should extend the number of detectable patterns of transcriptional modulation. RESULTS: We validated this new statistical approach on a microarray experiment that captures the temporal transcriptional profile of muscle differentiation in mouse C2C12 cells. Comparisons of the transcriptional state of myoblasts and differentiated myotubes via a robust variance test implicated several novel pathways in muscle cell differentiation previously overlooked by a standard enrichment analysis. Specifically, pathways involved in cell structure, calcium-mediated signaling and muscle-specific signaling were identified as differentially modulated based on their increased transcriptional variance. These biologically relevant results validate this approach and demonstrate the flexible nature of pathway-based methods of data analysis. AVAILABILITY: The software is available as Supplementary Material.

MAP1D, a novel methionine aminopeptidase family member is overexpressed in colon cancer.

N-terminal methionine removal is an important cellular process required for proper biological activity, subcellular localization, and eventual degradation of many proteins. The enzymes that catalyze this reaction are called Methionine Aminopeptidases (MAPs). To date, only two MAP family members, MAP1A and MAP2, have been well characterized and studied in mammals. In our studies, we have cloned a full length MAP1D gene. Expression and purification of full length recombinant protein shows that the sequence encodes an enzyme with MAP activity. MAP1D is overexpressed in colon cancer cell lines and in colon tumors as compared to matched normal tissue samples. Downregulation of MAP1D expression by shRNA in HCT-116 colon carcinoma cells reduces anchorage-independant growth in soft agar. These data suggest that MAP1D is a potentially oncogenic, novel member of the MAP gene family that may play an important role in colon tumorigenesis.

Identification of novel pathway regulation during myogenic differentiation.

Stem cell differentiation is governed by extracellular signals that activate intracellular networks (or pathways) to drive phenotypic specification. Using a novel gene clustering strategy we determined pathway relationships from a genome-wide transcriptional dataset of skeletal myoblast differentiation. Established myogenic pathways, including cell contractility and cell-cycle arrest, were predicted with extreme statistical significance (p approximately 0). In addition, gene sets associated with angiogenesis, neuronal activity, and mRNA splicing were regulated, exposing developmental and therapeutic implications. Acquisition of transcriptional data spanning the entire differentiation time course provided context for a dynamic landscape of functional pathway regulation. This novel perspective on myogenic cell differentiation revealed previously unrecognized patterns of regulation. We predict that similar analyses will facilitate ongoing efforts to define molecular mechanisms in other stem cell and developmental paradigms. Finally, by combining an iterative process of analysis with supplementation of novel pathways, this application may evolve into a powerful discovery tool.

Clustering protein sequences with a novel metric transformed from sequence similarity scores and sequence alignments with neural networks.

BACKGROUND: The sequencing of the human genome has enabled us to access a comprehensive list of genes (both experimental and predicted) for further analysis. While a majority of the approximately 30,000 known and predicted human coding genes are characterized and have been assigned at least one function, there remains a fair number of genes (about 12,000) for which no annotation has been made. The recent sequencing of other genomes has provided us with a huge amount of auxiliary sequence data which could help in the characterization of the human genes. Clustering these sequences into families is one of the first steps to perform comparative studies across several genomes. RESULTS: Here we report a novel clustering algorithm (CLUGEN) that has been used to cluster sequences of experimentally verified and predicted proteins from all sequenced genomes using a novel distance metric which is a neural network score between a pair of protein sequences. This distance metric is based on the pairwise sequence similarity score and the similarity between their domain structures. The distance metric is the probability that a pair of protein sequences are of the same Interpro family/domain, which facilitates the modelling of transitive homology closure to detect remote homologues. The hierarchical average clustering method is applied with the new distance metric. CONCLUSION: Benchmarking studies of our algorithm versus those reported in the literature shows that our algorithm provides clustering results with lower false positive and false negative rates. The clustering algorithm is applied to cluster several eukaryotic genomes and several dozens of prokaryotic genomes.

Expression and characterization of a 5-oxo-6E,8Z,11Z,14Z-eicosatetraenoic acid receptor highly expressed on human eosinophils and neutrophils.

Using a bioinformatics approach, we have isolated a novel G-protein-coupled receptor (GPCR), R527, and have demonstrated that this receptor shows no significant homology to previously deorphanized GPCRs. Quantitative reverse transcription-polymerase chain reaction analysis of the expression of GPCR R527 indicated a very high level of mRNA expression in eosinophils, with high expression also detected in neutrophils and lung macrophages. Stable cell lines were generated expressing this receptor together with the G-protein alpha-subunit G alpha(16). These cells were used to screen an agonist collection in a calcium mobilization assay and 5-oxo-6E,8Z,11Z,14Z-eicosatetraenoic acid (5-oxo-ETE) was identified as a putative ligand. 5(S)-hydroxyperoxy-6E,8Z,11Z,14Z-eicosatetraenoic acid was also shown to activate the receptor, whereas the leukotrienes LTB(4), LTC(4), LTD(4), and LTE(4) failed to elicit a response. In cAMP assays, pertussis toxin reversed the inhibitory effects of 5-oxo-ETE on cAMP production, indicating that the receptor is G alpha(i)-coupled. The GPCR R527 shows pharmacological properties similar to those of the previously described 5-oxo-ETE receptor expressed on eosinophils, neutrophils, and monocytes. These cell types show chemotactic responses to 5-oxo-ETE, and this eicosanoid has been proposed to play a key role in the inflammatory response. The molecular identification of a receptor binding 5-oxo-ETE will expand our understanding of the physiological role of this mediator and may provide new therapeutic opportunities.

On-resin macrocyclization of peptides via intramolecular SnAr reactions.

On-resin macrocyclization via an SNAr reaction is employed in the synthesis of tocinoic acid analogs. Specifically, an N-terminal nitrofluorobenzene is attacked by a nucleophilic C-terminal sidechain. The remaining nitro group can be reduced and acylated. NMR is used to compare the conformation of the new macrocyclic peptides to tocinoic acid.

Heterologously expressed inner lipoyl domain of dihydrolipoyl acetyltransferase inhibits ATP-dependent inactivation of pyruvate dehydrogenase complex. Identification of important amino acid residues.

The activity of the pyruvate dehydrogenase multienzyme complex (PDC), which catalyses the oxidation of pyruvate to acetyl-CoA within the mitochondrion, is diminished in animal models of diabetes. Studies with purified PDC components have suggested that the kinases responsible for inactivating the decarboxylase catalytic subunits of the complex are most efficient in their regulatory role when they are bound to dihydrolipoyl acetyltransferase (E2) subunits, which form the structural core of the complex. We report that the addition of an exogenous E2 subdomain (inner lipoyl domain) to an intact PDC inhibits ATP-dependent inactivation of the complex. By combining molecular modelling, site-directed mutagenesis and biophysical characterizations, we have also identified two amino acid residues in this subdomain (Ile229 and Phe231) that largely determine the magnitude of this effect.

Transfer nuclear Overhauser effect study of the conformation of oxytocin bound to bovine neurophysin I.

This study reports the structure of the peptide hormone oxytocin bound to its carrier protein, neurophysin I, obtained by nuclear magnetic resonance techniques. At the pH value of 2.1 in our experiments, the ligand is in fast exchange with its carrier protein, allowing the use of transfer-NOE methods. The number of distance constraints for the peptide being limited, considerable attention has been paid to an accurate distance determination. The resulting accurate distance limits were used as input for a distance geometry calculation followed by a restrained molecular dynamics run. Convergence to a well-defined family of structures for oxytocin in its bound state was reached. Both the backbone and the side-chain conformations differ between the bound form and the crystal structure of free oxytocin [Wood, S. P., et al. (1986) Science 232, 633]. These differences, as well as other structural features of the bound form, are discussed in terms of interactions made with the carrier protein. Transfer-NOE experiments at low peptide protein ratios provide direct experimental evidence for contacts between the oxytocin Tyr2 residue and an aromatic residue of neurophysin. The resonance assignments of the aromatic groups [Whittaker, B. A., et al. (1985) Biochemistry 24, 2782] together with the recently published X-ray structure of the neurophysin II protein complexed with a dipeptide [Chen et al. (1991) Proc. Natl. Acad. Sci. U.S.A. 88, 4240] allow us to assign the aromatic signal on the protein to the neurophysin Phe22 residue.

NMR studies of structure and dynamics of isotope enriched proteins.

Structural studies of globular proteins by nmr can be enhanced by the use of isotope enrichment. We have been working with proteins enriched with 15N, and with both 15N and 13C. Due to the isotope enrichment we could assign several large proteins with up to 186 residues and could address structural questions. Furthermore, we can accurately measure heteronuclear and homonuclear vicinal coupling constants. This involves in part multidimensional multiple resonance experiments. This is important for characterization of minor conformational changes caused by mutations. We have also made use of isotope enrichment to study the internal mobility of proteins. We also have developed novel methods for measuring accurately 15N relaxation parameters, in particular transverse relaxation rates. This has led us toward a method for directly mapping spectral density functions of the rotational motions of N-H bond vectors in proteins. The protein systems that are discussed include the unlabeled proteins kistrin and cytochrome c551, and the labeled proteins eglin c, a flavodoxin, and human dihydrofolate reductase.

Sequence-specific 1H and 15N resonance assignments for human dihydrofolate reductase in solution.

Dihydrofolate reductase is an intracellular target enzyme for folate antagonists, including the anticancer drug methotrexate. In order to design novel drugs with altered binding properties, a detailed description of protein-drug interactions in solution is desirable to understand the specificity of drug binding. As a first step in this process, heteronuclear three-dimensional NMR spectroscopy has been used to make sequential resonance assignments for more than 90% of the residues in human dihydrofolate reductase complexed with methotrexate. Uniform enrichment of the 21.5-kDa protein with 15N was required to obtain the resonance assignments via heteronuclear 3D NMR spectroscopy since homonuclear 2D spectra did not provide sufficient 1H resonance dispersion. Medium- and long-range NOE's have been used to characterize the secondary structure of the binary ligand-enzyme complex in solution.

Methotrexate binds in a non-productive orientation to human dihydrofolate reductase in solution, based on NMR spectroscopy.

Dihydrofolate reductase (DHFR) is an intracellular target enzyme for folate antagonist drugs, including methotrexate. In order to compare the binding of methotrexate to human DHFR in solution with that observed in the crystalline state, NMR spectroscopy has been used to determine the conformation of the drug bound to human DHFR in solution. In agreement with what has been observed in the crystalline state, NOE's identified protein and methotrexate protons indicate that methotrexate binds in a non-productive orientation. In contrast to what has been reported for E. coli DHFR in solution, only one bound conformation of methotrexate is observed.