Accurate and rapid identification of perturbed pathways through the analysis of genome-wide expression profiles facilitates the generation of biological hypotheses. We propose a statistical framework for determining whether a specified group of genes for a pathway has a coordinated association with a phenotype of interest. Several issues on proper hypothesis-testing procedures are clarified. In particular, it is shown that the differences in the correlation structure of each set of genes can lead to a biased comparison among gene sets unless a normalization procedure is applied. We propose statistical tests for two important but different aspects of association for each group of genes. This approach has more statistical power than currently available methods and can result in the discovery of statistically significant pathways that are not detected by other methods. This method is applied to data sets involving diabetes, inflammatory myopathies, and Alzheimer's disease, using gene sets we compiled from various public databases. In the case of inflammatory myopathies, we have correctly identified the known cytotoxic T lymphocyte-mediated autoimmunity in inclusion body myositis. Furthermore, we predicted the presence of dendritic cells in inclusion body myositis and of an IFN-␣͞␤ response in dermatomyositis, neither of which was previously described. These predictions have been subsequently corroborated by immunohistochemistry.microarrays ͉ gene ontology ͉ normalization ͉ correlated data ͉ inflammatory myopathies

SummaryAlthough autism is a highly heritable neurodevelopmental disorder, attempts to identify specific susceptibility genes have thus far met with limited success 1. Genome-wide association studies (GWAS) using half a million or more markers, particularly those with very large sample sizes achieved through meta-analysis, have shown great success in mapping genes for other complex genetic traits (http://www.genome.gov/26525384). Consequently, we initiated a linkage and association mapping study using half a million genome-wide SNPs in a common set of 1,031 multiplex autism families (1,553 affected offspring). We identified regions of suggestive and significant linkage on chromosomes 6q27 and 20p13, respectively. Initial analysis did not yield genome-wide significant associations; however, genotyping of top hits in additional families revealed a SNP on chromosome 5p15 (between SEMA5A and TAS2R1) that was significantly associated with autism (P = 2 × 10−7). We also demonstrated that expression of SEMA5A is reduced in brains from autistic patients, further implicating SEMA5A as an autism susceptibility gene. The linkage regions reported here provide targets for rare variation screening while the discovery of a single novel association demonstrates the action of common variants.

The molecular chaperone HSP90 aids the maturation of a diverse but select set of metastable protein clients, many of which are key to a variety of signal transduction pathways. HSP90 function has been best investigated in animal and fungal systems, where inhibition of the chaperone has exceptionally diverse effects, ranging from reversing oncogenic transformation to preventing the acquisition of drug resistance. Inhibition of HSP90 in the model plant Arabidopsis thaliana uncovers novel morphologies dependent on normally cryptic genetic variation and increases stochastic variation inherent to developmental processes. The biochemical activity of HSP90 is strictly conserved between animals and plants. However, the substrates and pathways dependent on HSP90 in plants are poorly understood. Progress has been impeded by the necessity of reliance on light-sensitive HSP90 inhibitors due to redundancy in the A. thaliana HSP90 gene family. Here we present phenotypic and genome-wide expression analyses of A. thaliana with constitutively reduced HSP90 levels achieved by RNAi targeting. HSP90 reduction affects a variety of quantitative life-history traits, including flowering time and total seed set, increases morphological diversity, and decreases the developmental stability of repeated characters. Several morphologies are synergistically affected by HSP90 and growth temperature. Genome-wide expression analyses also suggest a central role for HSP90 in the genesis and maintenance of plastic responses. The expression results are substantiated by examination of the response of HSP90-reduced plants to attack by caterpillars of the generalist herbivore Trichoplusia ni. HSP90 reduction potentiates a more robust herbivore defense response. In sum, we propose that HSP90 exerts global effects on the environmental responsiveness of plants to many different stimuli. The comprehensive set of HSP90-reduced lines described here is a vital instrument to further examine the role of HSP90 as a central interface between organism, development, and environment.

National Human Genome Research Institute, Doris Duke Charitable Foundation, National Health Service Blood and Transplant, National Institute for Health Research, and Wellcome Trust.

The anterior heart field (AHF) comprises a population of mesodermal progenitor cells that are added to the nascent linear heart to give rise to the majority of the right ventricle, interventricular septum, and outflow tract in mammals and birds. The zinc finger transcription factor GATA4 functions as an integral member of the cardiac transcription factor network in the derivatives of the AHF. In addition to its role in cardiac differentiation, GATA4 is also required for cardiomyocyte replication, although the transcriptional targets of GATA4 required for proliferation have not been previously identified. In the present study, we disrupted Gata4 function exclusively in the AHF and its derivatives. Gata4 AHF knockout mice die by embryonic day 13.5 and exhibit hypoplasia of the right ventricular myocardium and interventricular septum and display profound ventricular septal defects. Loss of Gata4 function in the AHF results in decreased myocyte proliferation in the right ventricle, and we identified numerous cell cycle genes that are dependent on Gata4 by microarray analysis. We show that GATA4 is required for cyclin D2, cyclin A2, and Cdk4 expression in the right ventricle and that the Cyclin D2 and Cdk4 promoters are bound and activated by GATA4 via multiple consensus GATA binding sites in each gene's proximal promoter. These findings establish Cyclin D2 and Cdk4 as direct transcriptional targets of GATA4 and support a model in which GATA4 controls cardiomyocyte proliferation by coordinately regulating numerous cell cycle genes.The cardiac lineage in mammals is initially specified from the anterior lateral mesoderm at embryonic day 7.5 (E7.5) in the mouse. The nascent cardiac mesoderm migrates anteriolaterally and fuses ventrally in the embryo to form a linear tube. The linear tube elongates through the addition of cells from the second heart field to the arterial and venous poles (1,12,28). A more restricted, anterior subset of these cells are added only to the arterial pole from the pharyngeal and splanchnic mesoderm. These cells, referred as the anterior heart field (AHF), give rise to the outflow tract, right ventricle, and ventricular septum (1,9,11,27,81). As cells from the AHF are added, the heart bends toward the ventral side, undergoes rightward looping, expands dramatically, and is eventually remodeled into the mature, four-chambered organ (13, 66).Embryonic cardiomyocytes differentiate as they continue to proliferate (48, 52). At early stages in development, cardiomyocytes have a high proliferation rate, which decreases progressively in late gestation (67). The high rate of cell cycle activity during the early stages of cardiomyocyte differentiation contributes to the growth of the future chambers within the linear tube during looping morphogenesis (42). The trabecular myocardium has a high rate of proliferation at this stage. As ventricular volumes increase, the trabeculations become compressed within the ventricular wall, resulting in a significant increase in the thickness of the compact myocardium (...

Inclusion body myositis is a late onset treatment-refractory autoimmune disease of skeletal muscle associated with a blood autoantibody (anti-cN1A), an HLA autoimmune haplotype, and muscle pathology characterized by cytotoxic CD8+ T cell destruction of myofibres. Here, we report on translational studies of inclusion body myositis patient muscle compared with a diverse set of other muscle disease samples. Using available microarray data on 411 muscle samples from patients with inclusion body myositis (n = 40), other muscle diseases (n = 265), and without neuromuscular disease (normal, n = 106), we identified a signature of T-cell cytotoxicity in inclusion body myositis muscle coupled with a signature of highly differentiated CD8 T-cell effector memory and terminally differentiated effector cells. Further, we examined killer cell lectin-like receptor G1 (KLRG1) as a marker of this population of cells, demonstrated the correlation of KLRG1 gene expression with lymphocyte cytotoxicity across 28 870 human tissue samples, and identified the presence of KLRG1 on pathogenic inclusion body myositis muscle invading T cells and an increase in KLRG1 expressing T cells in inclusion body myositis blood. We examined inclusion body myositis muscle T-cell proliferation by Ki67 immunohistochemistry demonstrating that diseased muscle-invading T cells are minimally or non-proliferative, in accordance with known properties of highly differentiated or terminally differentiated T cells. We found low expression of KLRG1 on infection-protective human lymphoid tissue central memory T cells and autoimmune-protective human blood regulatory T cells. Targeting highly differentiated cytotoxic T cells could be a favourable approach to treatment of inclusion body myositis.

The developmental and epileptic encephalopathies (DEEs) are heterogeneous disorders with a strong genetic contribution, but the underlying genetic etiology remains unknown in a significant proportion of individuals. To explore whether statistical support for genetic etiologies can be generated on the basis of phenotypic features, we analyzed whole-exome sequencing data and phenotypic similarities by using Human Phenotype Ontology (HPO) in 314 individuals with DEEs. We identified a de novo c.508C>T (p.Arg170Trp) variant in AP2M1 in two individuals with a phenotypic similarity that was higher than expected by chance (p ¼ 0.003) and a phenotype related to epilepsy with myoclonic-atonic seizures. We subsequently found the same de novo variant in two individuals with neurodevelopmental disorders and generalized epilepsy in a cohort of 2,310 individuals who underwent diagnostic whole-exome sequencing. AP2M1 encodes the m-subunit of the adaptor protein complex 2 (AP-2), which is involved in clathrin-mediated endocytosis (CME) and synaptic vesicle recycling. Modeling of protein dynamics indicated that the p.Arg170Trp variant impairs the conformational activation and thermodynamic entropy of the AP-2 complex. Functional complementation of both the m-subunit carrying the p.Arg170Trp variant in human cells and astrocytes derived from AP-2m conditional knockout mice revealed a significant impairment of CME of transferrin. In contrast, stability, expression levels, membrane recruitment, and localization were not impaired, suggesting a functional alteration of the AP-2 complex as the underlying disease mechanism. We establish a recurrent pathogenic variant in AP2M1 as a cause of DEEs with distinct phenotypic features, and we implicate dysfunction of the early steps of endocytosis as a disease mechanism in epilepsy.

No abstract