Large-scale organoid study suggests effects of trisomy 21 on early fetal neurodevelopment are more subtle than variability between isogenic lines and experiments

This study examines cortical organoids generated from a panel of isogenic trisomic and disomic iPSC lines (subclones) as a model of early fetal brain development in Down syndrome (DS). An initial experiment comparing organoids from one trisomic and one disomic line showed many genome-wide transcriptomic differences and modest differences in cell-type proportions, suggesting there may be a neurodevelopmental phenotype that is due to trisomy of chr21. To better control for multiple sources of variation, we undertook a highly robust study of ∼1,200 organoids using an expanded panel of six all-isogenic lines, three disomic, and three trisomic. The power of this experimental design was indicated by strong detection of the ∼1.5-fold difference in chr21 genes. However, the numerous expression differences in non-chr21 genes seen in the smaller experiment fell away, and the differences in cell-type representation between lines did not correlate with trisomy 21. Results suggest that the initial smaller experiment picked up differences between small organoid samples and individual isogenic lines, which "averaged out" in the larger panel of isogenic lines. Our results indicate that even when organoid and batch variability are better controlled for, variation between isogenic cell lines (even subclones) may obscure, or be conflated with, subtle neurodevelopmental phenotypes that may be present in ∼2nd trimester DS brain development. Interestingly, despite this variability between organoid batches and lines, and the "fetal stage" of these organoids, an increase in secreted Aβ40 peptide levels-an Alzheimer-related cellular phenotype-was more strongly associated with trisomy 21 status than were neurodevelopmental shifts in cell-type composition.

A Targeted Approach for Evaluating DUX4-Regulated Proteins as Potential Serum Biomarkers for Facioscapulohumeral Muscular Dystrophy Using Immunoassay Proteomics

BACKGROUND

Facioscapulohumeral muscular dystrophy (FSHD) is a progressive myopathy caused by misexpression of the double homeobox 4 (DUX4) embryonic transcription factor in skeletal muscle. Identifying quantitative and minimally invasive FSHD biomarkers to report on DUX4 activity will significantly accelerate therapeutic development.

OBJECTIVE

The goal of this study was to analyze secreted proteins known to be induced by DUX4 using the commercially available Olink Proteomics platform in order to identify potential blood-based molecular FSHD biomarkers.

METHODS

We used high-throughput, multiplex immunoassays from Olink Proteomics to measure the levels of several known DUX4-induced genes in a cellular myoblast model of FSHD, in FSHD patient-derived myotube cell cultures, and in serum from individuals with FSHD. Levels of other proteins on the Olink Proteomics panels containing these DUX4 targets were also examined in secondary exploratory analysis.

RESULTS

Placental alkaline phosphatase (ALPP) levels correlated with DUX4 expression in both cell-based FSHD systems but did not distinguish FSHD patient serum from unaffected controls.

CONCLUSIONS

ALPP, as measured with the Olink Proteomics platform, is not a promising FSHD serum biomarker candidate but could be utilized to evaluate DUX4 activity in discovery research efforts.

Poly(GR) and poly(GA) in cerebrospinal fluid as potential biomarkers for C9ORF72-ALS/FTD

GGGGCC repeat expansion in C9ORF72, which can be translated in both sense and antisense directions into five dipeptide repeat (DPR) proteins, including poly(GP), poly(GR), and poly(GA), is the most common genetic cause of amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD). Here we developed sensitive assays that can detect poly(GA) and poly(GR) in the cerebrospinal fluid (CSF) of patients with C9ORF72 mutations. CSF poly(GA) and poly(GR) levels did not correlate with age at disease onset, disease duration, or rate of decline of ALS Functional Rating Scale, and the average levels of these DPR proteins were similar in symptomatic and pre-symptomatic patients with C9ORF72 mutations. However, in a patient with C9ORF72-ALS who was treated with antisense oligonucleotide (ASO) targeting the aberrant C9ORF72 transcript, CSF poly(GA) and poly(GR) levels decreased approximately 50% within 6 weeks, indicating they may serve as sensitive fluid-based biomarkers in studies directed against the production of GGGGCC repeat RNAs or DPR proteins.

iMyoblasts for ex vivo and in vivo investigations of human myogenesis and disease modeling

Skeletal muscle myoblasts (iMyoblasts) were generated from human induced pluripotent stem cells (iPSCs) using an efficient and reliable transgene-free induction and stem cell selection protocol. Immunofluorescence, flow cytometry, qPCR, digital RNA expression profiling, and scRNA-Seq studies identify iMyoblasts as a PAX3+/MYOD1+ skeletal myogenic lineage with a fetal-like transcriptome signature, distinct from adult muscle biopsy myoblasts (bMyoblasts) and iPSC-induced muscle progenitors. iMyoblasts can be stably propagated for >12 passages or 30 population doublings while retaining their dual commitment for myotube differentiation and regeneration of reserve cells. iMyoblasts also efficiently xenoengrafted into irradiated and injured mouse muscle where they undergo differentiation and fetal-adult MYH isoform switching, demonstrating their regulatory plasticity for adult muscle maturation in response to signals in the host muscle. Xenograft muscle retains PAX3+ muscle progenitors and can regenerate human muscle in response to secondary injury. As models of disease, iMyoblasts from individuals with Facioscapulohumeral Muscular Dystrophy revealed a previously unknown epigenetic regulatory mechanism controlling developmental expression of the pathological DUX4 gene. iMyoblasts from Limb-Girdle Muscular Dystrophy R7 and R9 and Walker Warburg Syndrome patients modeled their molecular disease pathologies and were responsive to small molecule and gene editing therapeutics. These findings establish the utility of iMyoblasts for ex vivo and in vivo investigations of human myogenesis and disease pathogenesis and for the development of muscle stem cell therapeutics.

Applying genome-wide CRISPR-Cas9 screens for therapeutic discovery in facioscapulohumeral muscular dystrophy

The emergence of CRISPR-Cas9 gene-editing technologies and genome-wide CRISPR-Cas9 libraries enables efficient unbiased genetic screening that can accelerate the process of therapeutic discovery for genetic disorders. Here, we demonstrate the utility of a genome-wide CRISPR-Cas9 loss-of-function library to identify therapeutic targets for facioscapulohumeral muscular dystrophy (FSHD), a genetically complex type of muscular dystrophy for which there is currently no treatment. In FSHD, both genetic and epigenetic changes lead to misexpression of DUX4, the FSHD causal gene that encodes the highly cytotoxic DUX4 protein. We performed a genome-wide CRISPR-Cas9 screen to identify genes whose loss-of-function conferred survival when DUX4 was expressed in muscle cells. Genes emerging from our screen illuminated a pathogenic link to the cellular hypoxia response, which was revealed to be the main driver of DUX4-induced cell death. Application of hypoxia signaling inhibitors resulted in increased DUX4 protein turnover and subsequent reduction of the cellular hypoxia response and cell death. In addition, these compounds proved successful in reducing FSHD disease biomarkers in patient myogenic lines, as well as improving structural and functional properties in two zebrafish models of FSHD. Our genome-wide perturbation of pathways affecting DUX4 expression has provided insight into key drivers of DUX4-induced pathogenesis and has identified existing compounds with potential therapeutic benefit for FSHD. Our experimental approach presents an accelerated paradigm toward mechanistic understanding and therapeutic discovery of a complex genetic disease, which may be translatable to other diseases with well-established phenotypic selection assays.

Salivary metabolite levels in perinatally HIV-infected youth with periodontal disease

INTRODUCTION

Salivary metabolite profiles are altered in adults with HIV compared to their uninfected counterparts. Less is known about youth with HIV and how oral disorders that commonly accompany HIV infection impact salivary metabolite levels.

OBJECTIVE

As part of the Adolescent Master Protocol multi-site cohort study of the Pediatric HIV/AIDS Cohort Study (PHACS) network we compared the salivary metabolome of youth with perinatally-acquired HIV (PHIV) and youth HIV-exposed, but uninfected (PHEU) and determined whether metabolites differ in PHIV versus PHEU.

METHODS

We used three complementary targeted and discovery-based liquid chromatography-tandem mass spectrometry (LC-MS/MS) workflows to characterize salivary metabolite levels in 20 PHIV and 20 PHEU youth with and without moderate periodontitis. We examined main effects associated with PHIV and periodontal disease, and the interaction between them.

RESULTS

We did not identify differences in salivary metabolite profiles that remained significant under stringent control for both multiple between-group comparisons and multiple metabolites. Levels of cadaverine, a known periodontitis-associated metabolite, were more abundant in individuals with periodontal disease with the difference being more pronounced in PHEU than PHIV. In the discovery-based dataset, we identified a total of 564 endogenous peptides in the metabolite extracts, showing that proteolytic processing and amino acid metabolism are important to consider in the context of HIV infection.

CONCLUSION

The salivary metabolite profiles of PHIV and PHEU youth were overall very similar. Individuals with periodontitis particularly among the PHEU youth had higher levels of cadaverine, suggesting that HIV infection, or its treatment, may influence the metabolism of oral bacteria.

Ly6a Differential Expression in Blood-Brain Barrier Is Responsible for Strain Specific Central Nervous System Transduction Profile of AAV-PHP.B

Adeno-associated virus (AAV) gene therapy for neurological diseases was revolutionized by the discovery that AAV9 crosses the blood-brain barrier (BBB) after systemic administration. Transformative results have been documented in various inherited diseases, but overall neuronal transduction efficiency is relatively low. The recent development of AAV-PHP.B with ∼60-fold higher efficiency than AAV9 in transducing the adult mouse brain was the major first step toward acquiring the ability to deliver genes to the majority of cells in the central nervous system (CNS). However, little is known about the mechanism utilized by AAV to cross the BBB, and how it may diverge across species. In this study, we show that AAV-PHP.B is ineffective for systemic CNS gene transfer in the inbred strains BALB/cJ, BALB/cByJ, A/J, NOD/ShiLtJ, NZO/HILtJ, C3H/HeJ, and CBA/J mice, but it is highly potent in C57BL/6J, FVB/NJ, DBA/2J, 129S1/SvImJ, and AKR/J mice and also the outbred strain CD-1. We used the power of classical genetics to uncover the molecular mechanisms AAV-PHP.B engages to transduce CNS at high efficiency, and by quantitative trait locus mapping we identify a 6 Mb region in chromosome 15 with an logarithm of the odds (LOD) score ∼20, including single nucleotide polymorphisms in the coding region of 9 different genes. Comparison of the publicly available data on the genome sequence of 16 different mouse strains, combined with RNA-seq data analysis of brain microcapillary endothelia, led us to conclude that the expression level of Ly6a is likely the determining factor for differential efficacy of AAV-PHP.B in transducing the CNS across different mouse strains.

Precise therapeutic gene correction by a simple nuclease-induced double-stranded break

Current programmable nuclease-based (e.g. CRISPR-Cas9) methods for precise correction of a disease-causing genetic mutation harness the Homology Directed Repair (HDR) pathway. However, this repair process requires co-delivery of an exogenous DNA donor to recode the sequence and can be inefficient in many cell types. Here, we show that disease-causing frameshift mutations resulting from microduplications can be efficiently reverted to the wild-type sequence simply by generating a double-strand break (DSB) near the center of the duplication. We demonstrate this in patient-derived cell lines for two diseases: Limb-Girdle Muscular Dystrophy 2G (LGMD2G)¹ and Hermansky-Pudlak Syndrome Type 1 (HPS1)². Clonal analysis of Streptococcus pyogenes Cas9 (SpyCas9) nuclease-treated LGMD2G iPSCs revealed that ~80% contained at least one wild-type allele and that this correction restored TCAP expression in LGMD2G iPSC-derived myotubes. Efficient genotypic correction was also observed upon SpyCas9 treatment of an HPS1 patient-derived B-lymphoblastoid cell line (B-LCL). Inhibition of PARP-1 (poly (ADP-ribose) polymerase) suppresses the nuclease-mediated collapse of the microduplication to the wild-type sequence, confirming that precise correction is mediated by the MMEJ (microhomology-mediated end joining) pathway. Analysis of editing by SpyCas9 and Lachnospiraceae bacterium ND2006 Cas12a (LbaCas12a) at non-pathogenic microduplications within the genome that range in length from 4 bp to 36 bp indicates that the correction strategy is broadly applicable to a wide range of microduplication lengths and can be initiated by a variety of nucleases. The simplicity, reliability and efficacy of this MMEJ-based therapeutic strategy should permit the development of nuclease-based gene correction therapies for a variety of diseases that are associated with microduplications.

Nuclear-Import Receptors Reverse Aberrant Phase Transitions of RNA-Binding Proteins with Prion-like Domains

RNA-binding proteins (RBPs) with prion-like domains (PrLDs) phase transition to functional liquids, which can mature into aberrant hydrogels composed of pathological fibrils that underpin fatal neurodegenerative disorders. Several nuclear RBPs with PrLDs, including TDP-43, FUS, hnRNPA1, and hnRNPA2, mislocalize to cytoplasmic inclusions in neurodegenerative disorders, and mutations in their PrLDs can accelerate fibrillization and cause disease. Here, we establish that nuclear-import receptors (NIRs) specifically chaperone and potently disaggregate wild-type and disease-linked RBPs bearing a NLS. Karyopherin-β2 (also called Transportin-1) engages PY-NLSs to inhibit and reverse FUS, TAF15, EWSR1, hnRNPA1, and hnRNPA2 fibrillization, whereas Importin-α plus Karyopherin-β1 prevent and reverse TDP-43 fibrillization. Remarkably, Karyopherin-β2 dissolves phase-separated liquids and aberrant fibrillar hydrogels formed by FUS and hnRNPA1. In vivo, Karyopherin-β2 prevents RBPs with PY-NLSs accumulating in stress granules, restores nuclear RBP localization and function, and rescues degeneration caused by disease-linked FUS and hnRNPA2. Thus, NIRs therapeutically restore RBP homeostasis and mitigate neurodegeneration.

A global analysis of IFT-A function reveals specialization for transport of membrane-associated proteins into cilia

Intraflagellar transport (IFT), which is essential for the formation and function of cilia in most organisms, is the trafficking of IFT trains (i.e. assemblies of IFT particles) that carry cargo within the cilium. Defects in IFT cause several human diseases. IFT trains contain the complexes IFT-A and IFT-B. To dissect the functions of these complexes, we studied a Chlamydomonas mutant that is null for the IFT-A protein IFT140. The mutation had no effect on IFT-B but destabilized IFT-A, preventing flagella assembly. Therefore, IFT-A assembly requires IFT140. Truncated IFT140, which lacks the N-terminal WD repeats of the protein, partially rescued IFT and supported formation of half-length flagella that contained normal levels of IFT-B but greatly reduced amounts of IFT-A. The axonemes of these flagella had normal ultrastructure and, as investigated by SDS-PAGE, normal composition. However, composition of the flagellar 'membrane+matrix' was abnormal. Analysis of the latter fraction by mass spectrometry revealed decreases in small GTPases, lipid-anchored proteins and cell signaling proteins. Thus, IFT-A is specialized for the import of membrane-associated proteins. Abnormal levels of the latter are likely to account for the multiple phenotypes of patients with defects in IFT140.This article has an associated First Person interview with the first author of the paper.

Individual epigenetic status of the pathogenic D4Z4 macrosatellite correlates with disease in facioscapulohumeral muscular dystrophy

Background

Both forms of facioscapulohumeral muscular dystrophy (FSHD) are associated with aberrant epigenetic regulation of the chromosome 4q35 D4Z4 macrosatellite. Chromatin changes due to large deletions of heterochromatin (FSHD1) or mutations in chromatin regulatory proteins (FSHD2) lead to relaxation of epigenetic repression and increased expression of the deleterious double homeobox 4 (DUX4) gene encoded within the distal D4Z4 repeat. However, many individuals with the genetic requirements for FSHD remain asymptomatic throughout their lives. Here we investigated family cohorts of FSHD1 individuals who were either affected (manifesting) or without any discernible weakness (nonmanifesting/asymptomatic) and their unaffected family members to determine if individual epigenetic status and stability of repression at the contracted 4q35 D4Z4 array in myocytes correlates with FSHD disease.

Results

Family cohorts were analyzed for DNA methylation on the distal pathogenic 4q35 D4Z4 repeat on permissive A-type subtelomeres. We found DNA hypomethylation in FSHD1-affected subjects, hypermethylation in healthy controls, and distinctly intermediate levels of methylation in nonmanifesting subjects. We next tested if these differences in DNA methylation had functional relevance by assaying DUX4-fl expression and the stability of epigenetic repression of DUX4-fl in myogenic cells. Treatment with drugs that alter epigenetic status revealed that healthy cells were refractory to treatment, maintaining stable repression of DUX4, while FSHD1-affected cells were highly responsive to treatment and thus epigenetically poised to express DUX4. Myocytes from nonmanifesting subjects had significantly higher levels of DNA methylation and were more resistant to DUX4 activation in response to epigenetic drug treatment than cells from FSHD1-affected first-degree relatives containing the same contraction, indicating that the epigenetic status of the contracted D4Z4 array is reflective of disease.

Conclusions

The epigenetic status of the distal 4qA D4Z4 repeat correlates with FSHD disease; FSHD-affected subjects have hypomethylation, healthy unaffected subjects have hypermethylation, and nonmanifesting subjects have characteristically intermediate methylation. Thus, analysis of DNA methylation at the distal D4Z4 repeat could be used as a diagnostic indicator of developing clinical FSHD. In addition, the stability of epigenetic repression upstream of DUX4 expression is a key regulator of disease and a viable therapeutic target.

Electronic supplementary material

The online version of this article (doi:10.1186/s13148-015-0072-6) contains supplementary material, which is available to authorized users.

Identifying diagnostic DNA methylation profiles for facioscapulohumeral muscular dystrophy in blood and saliva using bisulfite sequencing

Background

Facioscapulohumeral muscular dystrophy (FSHD) is linked to chromatin relaxation due to epigenetic changes at the 4q35 D4Z4 macrosatellite array. Molecular diagnostic criteria for FSHD are complex and involve analysis of high molecular weight (HMW) genomic DNA isolated from lymphocytes, followed by multiple restriction digestions, pulse-field gel electrophoresis (PFGE), and Southern blotting. A subject is genetically diagnosed as FSHD1 if one of the 4q alleles shows a contraction in the D4Z4 array to below 11 repeats, while maintaining at least 1 repeat, and the contraction is in cis with a disease-permissive A-type subtelomere. FSHD2 is contraction-independent and cannot be diagnosed or excluded by this common genetic diagnostic procedure. However, FSHD1 and FSHD2 are linked by epigenetic deregulation, assayed as DNA hypomethylation, of the D4Z4 array on FSHD-permissive alleles. We have developed a PCR-based assay that identifies the epigenetic signature for both types of FSHD, distinguishing FSHD1 from FSHD2, and can be performed on genomic DNA isolated from blood, saliva, or cultured cells.

Results

Samples were obtained from healthy controls or patients clinically diagnosed with FSHD, and include both FSHD1 and FSHD2. The genomic DNAs were subjected to bisulfite sequencing analysis for the distal 4q D4Z4 repeat with an A-type subtelomere and the DUX4 5’ promoter region. We compared genomic DNA isolated from saliva and blood from the same individuals and found similar epigenetic signatures. DNA hypomethylation was restricted to the contracted 4qA chromosome in FSHD1 patients while healthy control subjects were hypermethylated. Candidates for FSHD2 showed extreme DNA hypomethylation on the 4qA DUX4 gene body as well as all analyzed DUX4 5’ sequences. Importantly, our assay does not amplify the D4Z4 arrays with non-permissive B-type subtelomeres and accurately excludes the arrays with non-permissive A-type subtelomeres.

Conclusions

We have developed an assay to identify changes in DNA methylation on the pathogenic distal 4q D4Z4 repeat. We show that the DNA methylation profile of saliva reflects FSHD status. This assay can distinguish FSHD from healthy controls, differentiate FSHD1 from FSHD2, does not require HMW genomic DNA or PFGE, and can be performed on either cultured cells, tissue, blood, or saliva samples.

Electronic supplementary material

The online version of this article (doi:10.1186/1868-7083-6-23) contains supplementary material, which is available to authorized users.

PLAAC: a web and command-line application to identify proteins with prion-like amino acid composition

Prions are self-templating protein aggregates that stably perpetuate distinct biological states and are of keen interest to researchers in both evolutionary and biomedical science. The best understood prions are from yeast and have a prion-forming domain with strongly biased amino acid composition, most notably enriched for Q or N. PLAAC is a web application that scans protein sequences for domains with P: rion- L: ike A: mino A: cid C: omposition. Users can upload sequence files, or paste sequences directly into a textbox. PLAAC ranks the input sequences by several summary scores and allows scores along sequences to be visualized. Text output files can be downloaded for further analyses, and visualizations saved in PDF and PNG formats.

AVAILABILITY AND IMPLEMENTATION

http://plaac.wi.mit.edu/. The Ruby-based web framework and the command-line software (implemented in Java, with visualization routines in R) are available at http://github.com/whitehead/plaac under the MIT license. All software can be run under OS X, Windows and Unix.

Human skeletal muscle xenograft as a new preclinical model for muscle disorders

Development of novel therapeutics requires good animal models of disease. Disorders for which good animal models do not exist have very few drugs in development or clinical trial. Even where there are accepted, albeit imperfect models, the leap from promising preclinical drug results to positive clinical trials commonly fails, including in disorders of skeletal muscle. The main alternative model for early drug development, tissue culture, lacks both the architecture and, usually, the metabolic fidelity of the normal tissue in vivo. Herein, we demonstrate the feasibility and validity of human to mouse xenografts as a preclinical model of myopathy. Human skeletal muscle biopsies transplanted into the anterior tibial compartment of the hindlimbs of NOD-Rag1(null) IL2rγ(null) immunodeficient host mice regenerate new vascularized and innervated myofibers from human myogenic precursor cells. The grafts exhibit contractile and calcium release behavior, characteristic of functional muscle tissue. The validity of the human graft as a model of facioscapulohumeral muscular dystrophy is demonstrated in disease biomarker studies, showing that gene expression profiles of xenografts mirror those of the fresh donor biopsies. These findings illustrate the value of a new experimental model of muscle disease, the human muscle xenograft in mice, as a feasible and valid preclinical tool to better investigate the pathogenesis of human genetic myopathies and to more accurately predict their response to novel therapeutics.

A new approach for the study of lung smooth muscle phenotypes and its application in a murine model of allergic airway inflammation

Phenotypes of lung smooth muscle cells in health and disease are poorly characterized. This is due, in part, to a lack of methodologies that allow for the independent and direct isolation of bronchial smooth muscle cells (BSMCs) and vascular smooth muscle cells (VSMCs) from the lung. In this paper, we describe the development of a bi-fluorescent mouse that permits purification of these two cell populations by cell sorting. By subjecting this mouse to an acute allergen based-model of airway inflammation that exhibits many features of asthma, we utilized this tool to characterize the phenotype of so-called asthmatic BSMCs. First, we examined the biophysical properties of single BSMCs from allergen sensitized mice and found increases in basal tone and cell size that were sustained ex vivo. We then generated for the first time, a comprehensive characterization of the global gene expression changes in BSMCs isolated from the bi-fluorescent mice with allergic airway inflammation. Using statistical methods and pathway analysis, we identified a number of differentially expressed mRNAs in BSMCs from allergen sensitized mice that code for key candidate proteins underlying changes in matrix formation, contractility, and immune responses. Ultimately, this tool will provide direction and guidance for the logical development of new markers and approaches for studying human lung smooth muscle.

Exome sequencing to identify de novo mutations in sporadic ALS trios

Amyotrophic lateral sclerosis (ALS) is a devastating neurodegenerative disease whose causes are still poorly understood. To identify additional genetic risk factors, we assessed the role of de novo mutations in ALS by sequencing the exomes of 47 ALS patients and both of their unaffected parents (n = 141 exomes). We found that amino acid-altering de novo mutations were enriched in genes encoding chromatin regulators, including the neuronal chromatin remodeling complex (nBAF) component SS18L1 (also known as CREST). CREST mutations inhibited activity-dependent neurite outgrowth in primary neurons, and CREST associated with the ALS protein FUS. These findings expand our understanding of the ALS genetic landscape and provide a resource for future studies into the pathogenic mechanisms contributing to sporadic ALS.

Stress granules as crucibles of ALS pathogenesis

Amyotrophic lateral sclerosis (ALS) is a fatal human neurodegenerative disease affecting primarily motor neurons. Two RNA-binding proteins, TDP-43 and FUS, aggregate in the degenerating motor neurons of ALS patients, and mutations in the genes encoding these proteins cause some forms of ALS. TDP-43 and FUS and several related RNA-binding proteins harbor aggregation-promoting prion-like domains that allow them to rapidly self-associate. This property is critical for the formation and dynamics of cellular ribonucleoprotein granules, the crucibles of RNA metabolism and homeostasis. Recent work connecting TDP-43 and FUS to stress granules has suggested how this cellular pathway, which involves protein aggregation as part of its normal function, might be coopted during disease pathogenesis.

Telomere position effect regulates DUX4 in human facioscapulohumeral muscular dystrophy

Telomeres may regulate human disease by at least two independent mechanisms. First, replicative senescence occurs once short telomeres generate DNA-damage signals that produce a barrier to tumor progression. Second, telomere position effects (TPE) could change gene expression at intermediate telomere lengths in cultured human cells. Here we report that telomere length may contribute to the pathogenesis of facioscapulohumeral muscular dystrophy (FSHD). FSHD is a late-onset disease genetically residing only 25-60 kilobases from the end of chromosome 4q. We used a floxable telomerase to generate isogenic clones with different telomere lengths from affected patients and their unaffected siblings. DUX4, the primary candidate for FSHD pathogenesis, is upregulated over ten-fold in FSHD myoblasts and myotubes with short telomeres, and its expression is inversely proportional to telomere length. FSHD may be the first known human disease in which TPE contributes to age-related phenotype.

Mutations in prion-like domains in hnRNPA2B1 and hnRNPA1 cause multisystem proteinopathy and ALS

Algorithms designed to identify canonical yeast prions predict that around 250 human proteins, including several RNA-binding proteins associated with neurodegenerative disease, harbour a distinctive prion-like domain (PrLD) enriched in uncharged polar amino acids and glycine. PrLDs in RNA-binding proteins are essential for the assembly of ribonucleoprotein granules. However, the interplay between human PrLD function and disease is not understood. Here we define pathogenic mutations in PrLDs of heterogeneous nuclear ribonucleoproteins (hnRNPs) A2B1 and A1 in families with inherited degeneration affecting muscle, brain, motor neuron and bone, and in one case of familial amyotrophic lateral sclerosis. Wild-type hnRNPA2 (the most abundant isoform of hnRNPA2B1) and hnRNPA1 show an intrinsic tendency to assemble into self-seeding fibrils, which is exacerbated by the disease mutations. Indeed, the pathogenic mutations strengthen a 'steric zipper' motif in the PrLD, which accelerates the formation of self-seeding fibrils that cross-seed polymerization of wild-type hnRNP. Notably, the disease mutations promote excess incorporation of hnRNPA2 and hnRNPA1 into stress granules and drive the formation of cytoplasmic inclusions in animal models that recapitulate the human pathology. Thus, dysregulated polymerization caused by a potent mutant steric zipper motif in a PrLD can initiate degenerative disease. Related proteins with PrLDs should therefore be considered candidates for initiating and perhaps propagating proteinopathies of muscle, brain, motor neuron and bone.

Facioscapulohumeral muscular dystrophy: Are telomeres the end of the story?

Facioscapulohumeral muscular dystrophy (FSHD) is a progressive myopathy with a relatively late age of onset (usually in the late teens) compared with Duchenne and many other muscular dystrophies. The current FSHD disease model postulates that contraction of the D4Z4 array at chromosome 4q35 leads to a more open chromatin conformation in that region and allows transcription of the DUX4 gene. DUX4 mRNA is stable only when transcribed from certain haplotypes that contain a polyadenylation signal. DUX4 protein is hypothesized to cause FSHD by mediating cytotoxicity and impairing skeletal muscle differentiation. We recently showed in a cell culture model that DUX4 expression is regulated by telomere length, suggesting that telomere shortening during aging may be partially responsible for the delayed onset and progressive nature of FSHD. We here put our data in the context of other recent findings arguing that progressive telomere shortening may play a critical role in FSHD but is not the whole story and that the current disease model needs additional refinement.

Transcriptional profiling in facioscapulohumeral muscular dystrophy to identify candidate biomarkers

Facioscapulohumeral muscular dystrophy (FSHD) is a progressive neuromuscular disorder caused by contractions of repetitive elements within the macrosatellite D4Z4 on chromosome 4q35. The pathophysiology of FSHD is unknown and, as a result, there is currently no effective treatment available for this disease. To better understand the pathophysiology of FSHD and develop mRNA-based biomarkers of affected muscles, we compared global analysis of gene expression in two distinct muscles obtained from a large number of FSHD subjects and their unaffected first-degree relatives. Gene expression in two muscle types was analyzed using GeneChip Gene 1.0 ST arrays: biceps, which typically shows an early and severe disease involvement; and deltoid, which is relatively uninvolved. For both muscle types, the expression differences were mild: using relaxed cutoffs for differential expression (fold change ≥1.2; nominal P value <0.01), we identified 191 and 110 genes differentially expressed between affected and control samples of biceps and deltoid muscle tissues, respectively, with 29 genes in common. Controlling for a false-discovery rate of <0.25 reduced the number of differentially expressed genes in biceps to 188 and in deltoid to 7. Expression levels of 15 genes altered in this study were used as a "molecular signature" in a validation study of an additional 26 subjects and predicted them as FSHD or control with 90% accuracy based on biceps and 80% accuracy based on deltoids.

Inhibiting GPI anchor biosynthesis in fungi stresses the endoplasmic reticulum and enhances immunogenicity

In fungi, the anchoring of proteins to the plasma membrane via their covalent attachment to glycosylphosphatidylinositol (GPI) is essential and thus provides a valuable point of attack for the development of antifungal therapeutics. Unfortunately, studying the underlying biology of GPI-anchor synthesis is difficult, especially in medically relevant fungal pathogens because they are not genetically tractable. Compounding difficulties, many of the genes in this pathway are essential in Saccharomyces cerevisiae. Here, we report the discovery of a new small molecule christened gepinacin (for GPI acylation inhibitor) which selectively inhibits Gwt1, a critical acyltransferase required for the biosynthesis of fungal GPI anchors. After delineating the target specificity of gepinacin using genetic and biochemical techniques, we used it to probe key, therapeutically relevant consequences of disrupting GPI anchor metabolism in fungi. We found that, unlike all three major classes of antifungals in current use, the direct antimicrobial activity of this compound results predominantly from its ability to induce overwhelming stress to the endoplasmic reticulum. Gepinacin did not affect the viability of mammalian cells nor did it inhibit their orthologous acyltransferase. This enabled its use in co-culture experiments to examine Gwt1's effects on host-pathogen interactions. In isolates of Candida albicans, the most common fungal pathogen in humans, exposure to gepinacin at sublethal concentrations impaired filamentation and unmasked cell wall β-glucan to stimulate a pro-inflammatory cytokine response in macrophages. Gwt1 is a promising antifungal drug target, and gepanacin is a useful probe for studying how disrupting GPI-anchor synthesis impairs viability and alters host-pathogen interactions in genetically intractable fungi.

Daily scheduled high fat meals moderately entrain behavioral anticipatory activity, body temperature, and hypothalamic c-Fos activation

When fed in restricted amounts, rodents show robust activity in the hours preceding expected meal delivery. This process, termed food anticipatory activity (FAA), is independent of the light-entrained clock, the suprachiasmatic nucleus, yet beyond this basic observation there is little agreement on the neuronal underpinnings of FAA. One complication in studying FAA using a calorie restriction model is that much of the brain is activated in response to this strong hunger signal. Thus, daily timed access to palatable meals in the presence of continuous access to standard chow has been employed as a model to study FAA in rats. In order to exploit the extensive genetic resources available in the murine system we extended this model to mice, which will anticipate rodent high fat diet but not chocolate or other sweet daily meals (Hsu, Patton, Mistlberger, and Steele; 2010, PLoS ONE e12903). In this study we test additional fatty meals, including peanut butter and cheese, both of which induced modest FAA. Measurement of core body temperature revealed a moderate preprandial increase in temperature in mice fed high fat diet but entrainment due to handling complicated interpretation of these results. Finally, we examined activation patterns of neurons by immunostaining for the immediate early gene c-Fos and observed a modest amount of entrainment of gene expression in the hypothalamus of mice fed a daily fatty palatable meal.

Facioscapulohumeral muscular dystrophy family studies of DUX4 expression: evidence for disease modifiers and a quantitative model of pathogenesis

Facioscapulohumeral muscular dystrophy (FSHD), the most prevalent myopathy afflicting both children and adults, is predominantly associated with contractions in the 4q35-localized macrosatellite D4Z4 repeat array. Recent studies have proposed that FSHD pathology is caused by the misexpression of the DUX4 (double homeobox 4) gene resulting in production of a pathogenic protein, DUX4-FL, which has been detected in FSHD, but not in unaffected control myogenic cells and muscle tissue. Here, we report the analysis of DUX4 mRNA and protein expression in a much larger collection of myogenic cells and muscle biopsies derived from biceps and deltoid muscles of FSHD affected subjects and their unaffected first-degree relatives. We confirmed that stable DUX4-fl mRNA and protein were expressed in myogenic cells and muscle tissues derived from FSHD affected subjects, including several genetically diagnosed adult FSHD subjects yet to show clinical manifestations of the disease in the assayed muscles. In addition, we report DUX4-fl mRNA and protein expression in muscle biopsies and myogenic cells from genetically unaffected relatives of the FSHD subjects, although at a significantly lower frequency. These results establish that DUX4-fl expression per se is not sufficient for FSHD muscle pathology and indicate that quantitative modifiers of DUX4-fl expression and/or function and family genetic background are determinants of FSHD muscle disease progression.

Food anticipatory activity behavior of mice across a wide range of circadian and non-circadian intervals

When rodents are fed in a limited amount during the daytime, they rapidly redistribute some of their nocturnal activity to the time preceding the delivery of food. In rats, anticipation of a daily meal has been interpreted as a circadian rhythm controlled by a food-entrained oscillator (FEO) with circadian limits to entrainment. Lesion experiments place this FEO outside of the light-entrainable circadian pacemaker in the suprachiasmatic nucleus. Mice also anticipate a fixed daily meal, but circadian limits to entrainment and anticipation of more than 2 daily meals, have not been assessed. We used a video-based behavior recognition system to quantify food anticipatory activity in mice receiving 2, 3, or 6 daily meals at intervals of 12, 8, or 4-hours (h). Individual mice were able to anticipate as many as 4 of 6 daily meals, and anticipation persisted during meal omission tests. On the 6 meal schedule, pre-prandial activity and body temperature were poorly correlated, suggesting independent regulation. Mice showed a limited ability to anticipate an 18 h feeding schedule. Finally, mice showed concurrent circadian and sub-hourly anticipation when provided with 6 small meals, at 30 minute intervals, at a fixed time of day. These results indicate that mice can anticipate feeding opportunities at a fixed time of day across a wide range of intervals not previously associated with anticipatory behavior in studies of rats. The methods described here can be exploited to determine the extent to which timing of different intervals in mice relies on common or distinct neural and molecular mechanisms.

Evaluating the role of the FUS/TLS-related gene EWSR1 in amyotrophic lateral sclerosis

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease affecting motor neurons. Mutations in related RNA-binding proteins TDP-43, FUS/TLS and TAF15 have been connected to ALS. These three proteins share several features, including the presence of a bioinformatics-predicted prion domain, aggregation-prone nature in vitro and in vivo and toxic effects when expressed in multiple model systems. Given these commonalities, we hypothesized that a related protein, EWSR1 (Ewing sarcoma breakpoint region 1), might also exhibit similar properties and therefore could contribute to disease. Here, we report an analysis of EWSR1 in multiple functional assays, including mutational screening in ALS patients and controls. We identified three missense variants in EWSR1 in ALS patients, which were absent in a large number of healthy control individuals. We show that disease-specific variants affect EWSR1 localization in motor neurons. We also provide multiple independent lines of in vitro and in vivo evidence that EWSR1 has similar properties as TDP-43, FUS and TAF15, including aggregation-prone behavior in vitro and ability to confer neurodegeneration in Drosophila. Postmortem analysis of sporadic ALS cases also revealed cytoplasmic mislocalization of EWSR1. Together, our studies highlight a potential role for EWSR1 in ALS, provide a collection of functional assays to be used to assess roles of additional RNA-binding proteins in disease and support an emerging concept that a class of aggregation-prone RNA-binding proteins might contribute broadly to ALS and related neurodegenerative diseases.

The tip of the iceberg: RNA-binding proteins with prion-like domains in neurodegenerative disease

Prions are self-templating protein conformers that are naturally transmitted between individuals and promote phenotypic change. In yeast, prion-encoded phenotypes can be beneficial, neutral or deleterious depending upon genetic background and environmental conditions. A distinctive and portable 'prion domain' enriched in asparagine, glutamine, tyrosine and glycine residues unifies the majority of yeast prion proteins. Deletion of this domain precludes prionogenesis and appending this domain to reporter proteins can confer prionogenicity. An algorithm designed to detect prion domains has successfully identified 19 domains that can confer prion behavior. Scouring the human genome with this algorithm enriches a select group of RNA-binding proteins harboring a canonical RNA recognition motif (RRM) and a putative prion domain. Indeed, of 210 human RRM-bearing proteins, 29 have a putative prion domain, and 12 of these are in the top 60 prion candidates in the entire genome. Startlingly, these RNA-binding prion candidates are inexorably emerging, one by one, in the pathology and genetics of devastating neurodegenerative disorders, including: amyotrophic lateral sclerosis (ALS), frontotemporal lobar degeneration with ubiquitin-positive inclusions (FTLD-U), Alzheimer's disease and Huntington's disease. For example, FUS and TDP-43, which rank 1st and 10th among RRM-bearing prion candidates, form cytoplasmic inclusions in the degenerating motor neurons of ALS patients and mutations in TDP-43 and FUS cause familial ALS. Recently, perturbed RNA-binding proteostasis of TAF15, which is the 2nd ranked RRM-bearing prion candidate, has been connected with ALS and FTLD-U. We strongly suspect that we have now merely reached the tip of the iceberg. We predict that additional RNA-binding prion candidates identified by our algorithm will soon surface as genetic modifiers or causes of diverse neurodegenerative conditions. Indeed, simple prion-like transfer mechanisms involving the prion domains of RNA-binding proteins could underlie the classical non-cell-autonomous emanation of neurodegenerative pathology from originating epicenters to neighboring portions of the nervous system. This article is part of a Special Issue entitled RNA-Binding Proteins.

Using the Mouse Grimace Scale to reevaluate the efficacy of postoperative analgesics in laboratory mice

Postoperative pain management in animals is complicated greatly by the inability to recognize pain. As a result, the choice of analgesics and their doses has been based on extrapolation from greatly differing pain models or the use of measures with unclear relevance to pain. We recently developed the Mouse Grimace Scale (MGS), a facial-expression-based pain coding system adapted directly from scales used in nonverbal human populations. The MGS has shown to be a reliable, highly accurate measure of spontaneous pain of moderate duration, and therefore is particularly useful in the quantification of postoperative pain. In the present study, we quantified the relative intensity and duration of postoperative pain after a sham ventral ovariectomy (laparotomy) in outbred mice. In addition, we compiled dose-response data for 4 commonly used analgesics: buprenorphine, carprofen, ketoprofen, and acetaminophen. We found that postoperative pain in mice, as defined by facial grimacing, lasts for 36 to 48 h, and appears to show relative exacerbation during the early dark (active) photophase. We find that buprenorphine was highly effective in inhibiting postoperative pain-induced facial grimacing in mice at doses equal to or lower than current recommendations, that carprofen and ketoprofen are effective only at doses markedly higher than those currently recommended, and that acetaminophen was ineffective at any dose used. We suggest the revision of practices for postoperative pain management in mice in light of these findings.

Opposing effects of glutamine and asparagine govern prion formation by intrinsically disordered proteins

Sequences rich in glutamine (Q) and asparagine (N) residues often fail to fold at the monomer level. This, coupled to their unusual hydrogen-bonding abilities, provides the driving force to switch between disordered monomers and amyloids. Such transitions govern processes as diverse as human protein-folding diseases, bacterial biofilm assembly, and the inheritance of yeast prions (protein-based genetic elements). A systematic survey of prion-forming domains suggested that Q and N residues have distinct effects on amyloid formation. Here, we use cell biological, biochemical, and computational techniques to compare Q/N-rich protein variants, replacing Ns with Qs and Qs with Ns. We find that the two residues have strong and opposing effects: N richness promotes assembly of benign self-templating amyloids; Q richness promotes formation of toxic nonamyloid conformers. Molecular simulations focusing on intrinsic folding differences between Qs and Ns suggest that their different behaviors are due to the enhanced turn-forming propensity of Ns over Qs.

The Rat Grimace Scale: a partially automated method for quantifying pain in the laboratory rat via facial expressions

We recently demonstrated the utility of quantifying spontaneous pain in mice via the blinded coding of facial expressions. As the majority of preclinical pain research is in fact performed in the laboratory rat, we attempted to modify the scale for use in this species. We present herein the Rat Grimace Scale, and show its reliability, accuracy, and ability to quantify the time course of spontaneous pain in the intraplantar complete Freund's adjuvant, intraarticular kaolin-carrageenan, and laparotomy (post-operative pain) assays. The scale's ability to demonstrate the dose-dependent analgesic efficacy of morphine is also shown. In addition, we have developed software, Rodent Face Finder^®, which successfully automates the most labor-intensive step in the process. Given the known mechanistic dissociations between spontaneous and evoked pain, and the primacy of the former as a clinical problem, we believe that widespread adoption of spontaneous pain measures such as the Rat Grimace Scale might lead to more successful translation of basic science findings into clinical application.

Gene expression profiling of skeletal muscles treated with a soluble activin type IIB receptor

Inhibition of the myostatin signaling pathway is emerging as a promising therapeutic means to treat muscle wasting and degenerative disorders. Activin type IIB receptor (ActRIIB) is the putative myostatin receptor, and a soluble activin receptor (ActRIIB-Fc) has been demonstrated to potently inhibit a subset of transforming growth factor (TGF)-β family members including myostatin. To determine reliable and valid biomarkers for ActRIIB-Fc treatment, we assessed gene expression profiles for quadriceps muscles from mice treated with ActRIIB-Fc compared with mice genetically lacking myostatin and control mice. Expression of 134 genes was significantly altered in mice treated with ActRIIB-Fc over a 2-wk period relative to control mice (fold change > 1.5, P < 0.001), whereas the number of significantly altered genes in mice treated for 2 days was 38, demonstrating a time-dependent response to ActRIIB-Fc in overall muscle gene expression. The number of significantly altered genes in Mstn(-/-) mice relative to control mice was substantially higher (360), but for most of these genes the expression levels in the 2-wk treated mice were closer to the levels in the Mstn(-/-) mice than in control mice (P < 10⁻³⁰). Expression levels of 30 selected genes were further validated with quantitative real-time polymerase chain reaction (qPCR), and a correlation of ≥ 0.89 was observed between the fold changes from the microarray analysis and the qPCR analysis. These data suggest that treatment with ActRIIB-Fc results in overlapping but distinct gene expression signatures compared with myostatin genetic mutation. Differentially expressed genes identified in this study can be used as potential biomarkers for ActRIIB-Fc treatment, which is currently in clinical trials as a therapeutic agent for muscle wasting and degenerative disorders.

Prion-like disorders: blurring the divide between transmissibility and infectivity

Prions are proteins that access self-templating amyloid forms, which confer phenotypic changes that can spread from individual to individual within or between species. These infectious phenotypes can be beneficial, as with yeast prions, or deleterious, as with mammalian prions that transmit spongiform encephalopathies. However, the ability to form self-templating amyloid is not unique to prion proteins. Diverse polypeptides that tend to populate intrinsically unfolded states also form self-templating amyloid conformers that are associated with devastating neurodegenerative disorders. Moreover, two RNA-binding proteins, FUS and TDP-43, which form cytoplasmic aggregates in amyotrophic lateral sclerosis, harbor a 'prion domain' similar to those found in several yeast prion proteins. Can these proteins and the neurodegenerative diseases to which they are linked become 'infectious' too? Here, we highlight advances that define the transmissibility of amyloid forms connected with Alzheimer's disease, Parkinson's disease and Huntington's disease. Collectively, these findings suggest that amyloid conformers can spread from cell to cell within the brains of afflicted individuals, thereby spreading the specific neurodegenerative phenotypes distinctive to the protein being converted to amyloid. Importantly, this transmissibility mandates a re-evaluation of emerging neuronal graft and stem-cell therapies. In this Commentary, we suggest how these treatments might be optimized to overcome the transmissible conformers that confer neurodegeneration.

Spontaneous generation of prion infectivity in fatal familial insomnia knockin mice

A crucial tenet of the prion hypothesis is that misfolding of the prion protein (PrP) induced by mutations associated with familial prion disease is, in an otherwise normal mammalian brain, sufficient to generate the infectious agent. Yet this has never been demonstrated. We engineered knockin mice to express a PrP mutation associated with a distinct human prion disease, fatal familial insomnia (FFI). An additional substitution created a strong transmission barrier against pre-existing prions. The mice spontaneously developed a disease distinct from that of other mouse prion models and highly reminiscent of FFI. Unique pathology was transmitted from FFI mice to mice expressing wild-type PrP sharing the same transmission barrier. FFI mice were highly resistant to infection by pre-existing prions, confirming infectivity did not arise from contaminating agents. Thus, a single amino acid change in PrP is sufficient to induce a distinct neurodegenerative disease and the spontaneous generation of prion infectivity.

Bridging high-throughput genetic and transcriptional data reveals cellular responses to alpha-synuclein toxicity

Cells respond to stimuli by changes in various processes, including signaling pathways and gene expression. Efforts to identify components of these responses increasingly depend on mRNA profiling and genetic library screens. By comparing the results of these two assays across various stimuli, we found that genetic screens tend to identify response regulators, whereas mRNA profiling frequently detects metabolic responses. We developed an integrative approach that bridges the gap between these data using known molecular interactions, thus highlighting major response pathways. We used this approach to reveal cellular pathways responding to the toxicity of alpha-synuclein, a protein implicated in several neurodegenerative disorders including Parkinson's disease. For this we screened an established yeast model to identify genes that when overexpressed alter alpha-synuclein toxicity. Bridging these data and data from mRNA profiling provided functional explanations for many of these genes and identified previously unknown relations between alpha-synuclein toxicity and basic cellular pathways.

Combining guilt-by-association and guilt-by-profiling to predict Saccharomyces cerevisiae gene function

Background:

Learning the function of genes is a major goal of computational genomics. Methods for inferring gene function have typically fallen into two categories: 'guilt-by-profiling', which exploits correlation between function and other gene characteristics; and 'guilt-by-association', which transfers function from one gene to another via biological relationships.

Results:

We have developed a strategy ('Funckenstein') that performs guilt-by-profiling and guilt-by-association and combines the results. Using a benchmark set of functional categories and input data for protein-coding genes in Saccharomyces cerevisiae, Funckenstein was compared with a previous combined strategy. Subsequently, we applied Funckenstein to 2,455 Gene Ontology terms. In the process, we developed 2,455 guilt-by-profiling classifiers based on 8,848 gene characteristics and 12 functional linkage graphs based on 23 biological relationships.

Conclusion:

Funckenstein outperforms a previous combined strategy using a common benchmark dataset. The combination of 'guilt-by-profiling' and 'guilt-by-association' gave significant improvement over the component classifiers, showing the greatest synergy for the most specific functions. Performance was evaluated by cross-validation and by literature examination of the top-scoring novel predictions. These quantitative predictions should help prioritize experimental study of yeast gene functions.

The evolution of bet-hedging adaptations to rare scenarios

When faced with a variable environment, organisms may switch between different strategies according to some probabilistic rule. In an infinite population, evolution is expected to favor the rule that maximizes geometric mean fitness. If some environments are encountered only rarely, selection may not be strong enough for optimal switching probabilities to evolve. Here we calculate the evolution of switching probabilities in a finite population by analyzing fixation probabilities of alleles specifying switching rules. We calculate the conditions required for the evolution of phenotypic switching as a form of bet-hedging as a function of the population size N, the rate theta at which a rare environment is encountered, and the selective advantage s associated with switching in the rare environment. We consider a simplified model in which environmental switching and phenotypic switching are one-way processes, and mutation is symmetric and rare with respect to the timescale of fixation events. In this case, the approximate requirements for bet-hedging to be favored by a ratio of at least R are that sN>log(R) and thetaN>square root R .

The power of automated high-resolution behavior analysis revealed by its application to mouse models of Huntington's and prion diseases

Automated analysis of mouse behavior will be vital for elucidating the genetic determinants of behavior, for comprehensive analysis of human disease models, and for assessing the efficacy of various therapeutic strategies and their unexpected side effects. We describe a video-based behavior-recognition technology to analyze home-cage behaviors and demonstrate its power by discovering previously unrecognized features of two already extensively characterized mouse models of neurodegenerative disease. The severe motor abnormalities in Huntington's disease mice manifested in our analysis by decreased hanging, jumping, stretching, and rearing. Surprisingly, behaviors such as resting and grooming were also affected. Unexpectedly, mice with infectious prion disease showed profound increases in activity at disease onset: rearing increased 2.5-fold, walking 10-fold and jumping 30-fold. Strikingly, distinct behaviors were altered specifically during day or night hours. We devised a systems approach for multiple-parameter phenotypic characterization and applied it to defining disease onset robustly and at early time points.

The loss of adaptive plasticity during long periods of environmental stasis

Adaptive plasticity allows populations to adjust rapidly to environmental change. If this is useful only rarely, plasticity may undergo mutational degradation and be lost from a population. We consider a population of constant size N undergoing loss of plasticity at functional mutation rate m and with selective advantage s associated with loss. Environmental change events occur at rate theta per generation, killing all individuals that lack plasticity. The expected time until loss of plasticity in a fluctuating environment is always at least tau, the expected time until loss of plasticity in a static environment. When mN > 1 and N theta >> 1, we find that plasticity will be maintained for an average of at least 10(8) generations in a single population, provided tau > 18/theta. In a metapopulation, plasticity is retained under the more lenient condition tau > 1.3/theta, irrespective of mN, for a modest number of demes. We calculate both exact and approximate solutions for tau and find that it is linearly dependent only on the logarithm of N, and so, surprisingly, both the population size and the number of demes in the metapopulation make little difference to the retention of plasticity. Instead, tau is dominated by the term 1/(m+s/2).

Motifs, themes and thematic maps of an integrated Saccharomyces cerevisiae interaction network

BACKGROUND

Large-scale studies have revealed networks of various biological interaction types, such as protein-protein interaction, genetic interaction, transcriptional regulation, sequence homology, and expression correlation. Recurring patterns of interconnection, or 'network motifs', have revealed biological insights for networks containing either one or two types of interaction.

RESULTS

To study more complex relationships involving multiple biological interaction types, we assembled an integrated Saccharomyces cerevisiae network in which nodes represent genes (or their protein products) and differently colored links represent the aforementioned five biological interaction types. We examined three- and four-node interconnection patterns containing multiple interaction types and found many enriched multi-color network motifs. Furthermore, we showed that most of the motifs form 'network themes' -- classes of higher-order recurring interconnection patterns that encompass multiple occurrences of network motifs. Network themes can be tied to specific biological phenomena and may represent more fundamental network design principles. Examples of network themes include a pair of protein complexes with many inter-complex genetic interactions -- the 'compensatory complexes' theme. Thematic maps -- networks rendered in terms of such themes -- can simplify an otherwise confusing tangle of biological relationships. We show this by mapping the S. cerevisiae network in terms of two specific network themes.

CONCLUSION

Significantly enriched motifs in an integrated S. cerevisiae interaction network are often signatures of network themes, higher-order network structures that correspond to biological phenomena. Representing networks in terms of network themes provides a useful simplification of complex biological relationships.

Combining biological networks to predict genetic interactions

Genetic interactions define overlapping functions and compensatory pathways. In particular, synthetic sick or lethal (SSL) genetic interactions are important for understanding how an organism tolerates random mutation, i.e., genetic robustness. Comprehensive identification of SSL relationships remains far from complete in any organism, because mapping these networks is highly labor intensive. The ability to predict SSL interactions, however, could efficiently guide further SSL discovery. Toward this end, we predicted pairs of SSL genes in Saccharomyces cerevisiae by using probabilistic decision trees to integrate multiple types of data, including localization, mRNA expression, physical interaction, protein function, and characteristics of network topology. Experimental evidence demonstrated the reliability of this strategy, which, when extended to human SSL interactions, may prove valuable in discovering drug targets for cancer therapy and in identifying genes responsible for multigenic diseases.

Comment on "Subgraphs in random networks"

We point out biases in the algorithms used by Itzkovitz et a. [Phys. Rev. E 68, 026127 (2003)] to assess their approximate formulas for the average number of occurrences of certain subgraphs in random graphs with prescribed degree sequences.

Predicting phenotype from patterns of annotation

MOTIVATION

Predicting the outcome of specific experiments (such as the growth of a particular mutant strain in a particular medium) has the potential to allow researchers to devote resources to experiments with higher expected numbers of 'hits'.

RESULTS

We use decision trees to predict phenotypes associated with Saccharomyces cerevisiae genes on the basis of Gene Ontology (GO) functional annotations from the Saccharomyces Genome Database (SGD) and other phenotypic annotations from the Yeast Phenotype Catalog at the Munich Information Center for Protein Sequences (MIPS). We assess the methodology in three ways: (1) we use cross-validation on the phenotypic annotations listed in MIPS, and show ROC curves indicating the tradeoff between true-positive rate and false-positive rate; (2) we do a literature-search for 100 of the predicted gene-phenotype associations that are not listed in MIPS, and find evidence for 43 of them; (3) we use deletion strains to experimentally assess 61 predicted gene-phenotype associations not listed in MIPS; significantly more of these deletion strains show abnormal growth than would be expected by chance.

Characterizing gene sets with FuncAssociate

SUMMARY

FuncAssociate is a web-based tool to help researchers use Gene Ontology attributes to characterize large sets of genes derived from experiment. Distinguishing features of FuncAssociate include the ability to handle ranked input lists, and a Monte Carlo simulation approach that is more appropriate to determine significance than other methods, such as Bonferroni or idák p-value correction. FuncAssociate currently supports 10 organisms (Vibrio cholerae, Shewanella oneidensis, Saccharomyces cerevisiae, Schizosaccharomyces pombe, Arabidopsis thaliana, Caenorhaebditis elegans, Drosophila melanogaster, Mus musculus, Rattus norvegicus and Homo sapiens).

AVAILABILITY

FuncAssociate is freely accessible at http://llama.med.harvard.edu/Software.html. Source code (in Perl and C) is freely available to academic users 'as is'.

Predicting protein complex membership using probabilistic network reliability

Evidence for specific protein-protein interactions is increasingly available from both small- and large-scale studies, and can be viewed as a network. It has previously been noted that errors are frequent among large-scale studies, and that error frequency depends on the large-scale method used. Despite knowledge of the error-prone nature of interaction evidence, edges (connections) in this network are typically viewed as either present or absent. However, use of a probabilistic network that considers quantity and quality of supporting evidence should improve inference derived from protein networks. Here we demonstrate inference of membership in a partially known protein complex by using a probabilistic network model and an algorithm previously used to evaluate reliability in communication networks.

Predicting co-complexed protein pairs using genomic and proteomic data integration

BACKGROUND

Identifying all protein-protein interactions in an organism is a major objective of proteomics. A related goal is to know which protein pairs are present in the same protein complex. High-throughput methods such as yeast two-hybrid (Y2H) and affinity purification coupled with mass spectrometry (APMS) have been used to detect interacting proteins on a genomic scale. However, both Y2H and APMS methods have substantial false-positive rates. Aside from high-throughput interaction screens, other gene- or protein-pair characteristics may also be informative of physical interaction. Therefore it is desirable to integrate multiple datasets and utilize their different predictive value for more accurate prediction of co-complexed relationship.

RESULTS

Using a supervised machine learning approach--probabilistic decision tree, we integrated high-throughput protein interaction datasets and other gene- and protein-pair characteristics to predict co-complexed pairs (CCP) of proteins. Our predictions proved more sensitive and specific than predictions based on Y2H or APMS methods alone or in combination. Among the top predictions not annotated as CCPs in our reference set (obtained from the MIPS complex catalogue), a significant fraction was found to physically interact according to a separate database (YPD, Yeast Proteome Database), and the remaining predictions may potentially represent unknown CCPs.

CONCLUSIONS

We demonstrated that the probabilistic decision tree approach can be successfully used to predict co-complexed protein (CCP) pairs from other characteristics. Our top-scoring CCP predictions provide testable hypotheses for experimental validation.

Intensity-based protein identification by machine learning from a library of tandem mass spectra

Tandem mass spectrometry (MS/MS) has emerged as a cornerstone of proteomics owing in part to robust spectral interpretation algorithms. Widely used algorithms do not fully exploit the intensity patterns present in mass spectra. Here, we demonstrate that intensity pattern modeling improves peptide and protein identification from MS/MS spectra. We modeled fragment ion intensities using a machine-learning approach that estimates the likelihood of observed intensities given peptide and fragment attributes. From 1,000,000 spectra, we chose 27,000 with high-quality, nonredundant matches as training data. Using the same 27,000 spectra, intensity was similarly modeled with mismatched peptides. We used these two probabilistic models to compute the relative likelihood of an observed spectrum given that a candidate peptide is matched or mismatched. We used a 'decoy' proteome approach to estimate incorrect match frequency, and demonstrated that an intensity-based method reduces peptide identification error by 50-96% without any loss in sensitivity.

A non-parametric model for transcription factor binding sites

We introduce a non-parametric representation of transcription factor binding sites which can model arbitrary dependencies between positions. As two parameters are varied, this representation smoothly interpolates between the empirical distribution of binding sites and the standard position-specific scoring matrix (PSSM). In a test of generalization to unseen binding sites using 10-fold cross-validation on known binding sites for 95 TRANSFAC transcription factors, this representation outperforms PSSMs on between 65 and 89 of the 95 transcription factors, depending on the choice of the two adjustable parameters. We also discuss how the non- parametric representation may be incorporated into frameworks for finding binding sites given only a collection of unaligned promoter regions.

Predicting gene function from patterns of annotation

The Gene Ontology (GO) Consortium has produced a controlled vocabulary for annotation of gene function that is used in many organism-specific gene annotation databases. This allows the prediction of gene function based on patterns of annotation. For example, if annotations for two attributes tend to occur together in a database, then a gene holding one attribute is likely to hold the other as well. We modeled the relationships among GO attributes with decision trees and Bayesian networks, using the annotations in the Saccharomyces Genome Database (SGD) and in FlyBase as training data. We tested the models using cross-validation, and we manually assessed 100 gene-attribute associations that were predicted by the models but that were not present in the SGD or FlyBase databases. Of the 100 manually assessed associations, 41 were judged to be true, and another 42 were judged to be plausible.