Cross-species Conservation Research Articles

Cholinergic urethral brush cells are widespread throughout placental mammals

We previously identified a population of cholinergic epithelial cells in murine, human and rat urethrae that exhibits a structural marker of brush cells (villin) and expresses components of the canonical taste transduction signaling cascade (α-gustducin, phospholipase Cβ2 (PLCβ2), transient receptor potential cation channel melanostatin 5 (TRPM5)). These cells serve as sentinels, monitoring the chemical composition of the luminal content for potentially hazardous compounds such as bacteria, and initiate protective reflexes counteracting further ingression. In order to elucidate cross-species conservation of the urethral chemosensory pathway we investigated the occurrence and molecular make-up of urethral brush cells in placental mammals. We screened 11 additional species, at least one in each of the five mammalian taxonomic units primates, carnivora, perissodactyla, artiodactyla and rodentia, for immunohistochemical labeling of the acetylcholine synthesizing enzyme, choline acetyltransferase (ChAT), villin, and taste cascade components (α-gustducin, PLCβ2, TRPM5). Corresponding to findings in previously investigated species, urethral epithelial cells with brush cell shape were immunolabeled in all 11 mammals. In 8 species, immunoreactivities against all marker proteins and ChAT were observed, and double-labeling immunofluorescence confirmed the cholinergic nature of villin-positive and chemosensory (TRPM5-positive) cells. In cat and horse, these cells were not labeled by the ChAT antiserum used in this study, and unspecific reactions of the secondary antiserum precluded conclusions about ChAT-expression in the bovine epithelium. These data indicate that urethral brush cells are widespread throughout the mammalian kingdom and evolved not later than about 64.5millionyears ago.

New insights on the transcriptional regulation of CD69 gene through a potent enhancer located in the conserved non-coding sequence 2

The CD69 type II C-type lectin is one of the earliest indicators of leukocyte activation acting in lymphocyte migration and cytokine secretion. CD69 expression in hematopoietic lineage undergoes rapid changes depending on the cell-lineage, the activation state or the localization of the cell where it is expressed, suggesting a complex and tightly controlled regulation. Here we provide new insights on the transcriptional regulation of CD69 gene in mammal species. Through in silico studies, we analyzed several regulatory features of the 4 upstream conserved non-coding sequences (CNS 1–4) previously described, confirming a major function of CNS2 in the transcriptional regulation of CD69. In addition, multiple transcription binding sites are identified in the CNS2 region by DNA cross-species conservation analysis. By functional approaches we defined a core region of 226bp located within CNS2 as the main enhancer element of CD69 transcription in the hematopoietic cells analyzed. By chromatin immunoprecipitation, binding of RUNX1 to the core-CNS2 was shown in a T cell line. In addition, we found an activating but not essential role of RUNX1 in CD69 gene transcription by site-directed mutagenesis and RNA silencing, probably through the interaction with this potent enhancer specifically in the hematopoietic lineage. In summary, in this study we contribute with new evidences to the landscape of the transcriptional regulation of the CD69 gene.

The distribution pattern of genetic variation in the transcript isoforms of the alternatively spliced protein-coding genes in the human genome.

By enabling the transcription of multiple isoforms from the same gene locus, alternative-splicing mechanisms greatly expand the diversity of the human transcriptome and proteome. Currently, the alternatively spliced transcripts from each protein-coding gene locus in the human genome can be classified as either principal or non-principal isoforms, providing that they differ with respect to cross-species conservation or biological features. By mapping the variants from the 1000 Genomes Project onto the coding region of each isoform, an interesting pattern of the genetic variation distributions of the coding regions for these two types of transcript isoforms was revealed on a whole-genome scale: compared with the principal isoform-specific coding regions, the non-principal isoform-specific coding regions are significantly enriched in amino acid-changing variants, particularly those that have a strong impact on protein function and have higher derived allele frequencies, suggesting that non-principal isoform-specific substitutions are less likely to be related to phenotype changes or disease. The results herein can help us better understand the potential consequences of alternatively spliced products from a population perspective.

Leveraging existing data for prioritization of the ecological risks of human and veterinary pharmaceuticals to aquatic organisms.

Medicinal innovation has led to the discovery and use of thousands of human and veterinary drugs. With this comes the potential for unintended effects on non-target organisms exposed to pharmaceuticals inevitably entering the environment. The impracticality of generating whole-organism chronic toxicity data representative of all species in the environment has necessitated prioritization of drugs for focused empirical testing as well as field monitoring. Current prioritization strategies typically emphasize likelihood for exposure (i.e. predicted/measured environmental concentrations), while incorporating only rather limited consideration of potential effects of the drug to non-target organisms. However, substantial mammalian pharmacokinetic and mechanism/mode of action (MOA) data are produced during drug development to understand drug target specificity and efficacy for intended consumers. An integrated prioritization strategy for assessing risks of human and veterinary drugs would leverage available pharmacokinetic and toxicokinetic data for evaluation of the potential for adverse effects to non-target organisms. In this reiview, we demonstrate the utility of read-across approaches to leverage mammalian absorption, distribution, metabolism and elimination data; analyse cross-species molecular target conservation and translate therapeutic MOA to an adverse outcome pathway(s) relevant to aquatic organisms as a means to inform prioritization of drugs for focused toxicity testing and environmental monitoring.

Deep conservation of genes required for both Drosphila melanogaster and Caenorhabditis elegans sleep includes a role for dopaminergic signaling.

Cross-species conservation of sleep-like behaviors predicts the presence of conserved molecular mechanisms underlying sleep. However, limited experimental evidence of conservation exists. Here, this prediction is tested directly. During lethargus, Caenorhabditis elegans spontaneously sleep in short bouts that are interspersed with bouts of spontaneous locomotion. We identified 26 genes required for Drosophila melanogaster sleep. Twenty orthologous C. elegans genes were selected based on similarity. Their effect on C. elegans sleep and arousal during the last larval lethargus was assessed. The 20 most similar genes altered both the quantity of sleep and arousal thresholds. In 18 cases, the direction of change was concordant with Drosophila studies published previously. Additionally, we delineated a conserved genetic pathway by which dopamine regulates sleep and arousal. In C. elegans neurons, G-alpha S, adenylyl cyclase, and protein kinase A act downstream of D1 dopamine receptors to regulate these behaviors. Finally, a quantitative analysis of genes examined herein revealed that C. elegans arousal thresholds were directly correlated with amount of sleep during lethargus. However, bout duration varies little and was not correlated with arousal thresholds. The comprehensive analysis presented here suggests that conserved genes and pathways are required for sleep in invertebrates and, likely, across the entire animal kingdom. The genetic pathway delineated in this study implicates G-alpha S and previously known genes downstream of dopamine signaling in sleep. Quantitative analysis of various components of quiescence suggests that interdependent or identical cellular and molecular mechanisms are likely to regulate both arousal and sleep entry.

Comparative Relationships and Chromosome Evolution in Switchgrass (Panicum virgatum) and Its Genomic Model, Foxtail Millet (Setaria italica)

A switchgrass (Panicum virgatum L.) consensus map was developed that combined data from two mapping populations and integrated recombination data from both parents of this largely obligate outcrossing species. The consensus map consisted of 1,321 loci and spanned 2,122 cM. An analysis of the distribution of genic and genomic markers across the length of a linkage group showed that genic markers were relatively overrepresented in distal regions, while genomic markers were relatively overrepresented in pericentromeric regions. Furthermore, genic markers located in distal and pericentromeric regions identified orthologs in the genome of the closely related species foxtail millet (Setaria italica (L.) P. Beauv.) at a similar rate, but fewer orthologs were found for genomic markers that mapped to pericentromeric regions compared to distal regions. Because cross-species sequence conservation is largely limited to genes, our data suggest that genomic markers that map to distal regions are more likely to be derived from genic regions than those that map to pericentromeric regions. A comparative analysis between loci on the switchgrass consensus map and their orthologs in the foxtail millet genome showed that the two species varied by at least nine inversions and one reciprocal translocation. Extending the comparative analysis to sorghum demonstrated that the majority of the rearrangements occurred in the foxtail millet genome in the past 13 million years. The reduction in chromosome number from 10 in the Andropogoneae to 9 in the Paniceae was achieved through a complex set of rearrangements involving three ancestral chromosomes orthologous to sorghum chromosomes 6, 8, and 9 and chromosomes III and VII in the lineage leading to switchgrass and foxtail millet. Upon insertion of ancestral chromosome 9 into ancestral chromosome 8, the centromere of the recipient chromosome was lost, and this was accompanied by a loss of repetitive DNA. The switchgrass–foxtail millet comparative map provides a guide of the rearrangements that need to be taken into account when using foxtail millet as a surrogate for switchgrass in genetic analyses.

Multiple evidence strands suggest that there may be as few as 19,000 human protein-coding genes.

Determining the full complement of protein-coding genes is a key goal of genome annotation. The most powerful approach for confirming protein-coding potential is the detection of cellular protein expression through peptide mass spectrometry (MS) experiments. Here, we mapped peptides detected in seven large-scale proteomics studies to almost 60% of the protein-coding genes in the GENCODE annotation of the human genome. We found a strong relationship between detection in proteomics experiments and both gene family age and cross-species conservation. Most of the genes for which we detected peptides were highly conserved. We found peptides for >96% of genes that evolved before bilateria. At the opposite end of the scale, we identified almost no peptides for genes that have appeared since primates, for genes that did not have any protein-like features or for genes with poor cross-species conservation. These results motivated us to describe a set of 2001 potential non-coding genes based on features such as weak conservation, a lack of protein features, or ambiguous annotations from major databases, all of which correlated with low peptide detection across the seven experiments. We identified peptides for just 3% of these genes. We show that many of these genes behave more like non-coding genes than protein-coding genes and suggest that most are unlikely to code for proteins under normal circumstances. We believe that their inclusion in the human protein-coding gene catalogue should be revised as part of the ongoing human genome annotation effort.

A surprising cross-species conservation in the genomic landscape of mouse and human oral cancer identifies a transcriptional signature predicting metastatic disease.

Improved understanding of the molecular basis underlying oral squamous cell carcinoma (OSCC) aggressive growth has significant clinical implications. Herein, cross-species genomic comparison of carcinogen-induced murine and human OSCCs with indolent or metastatic growth yielded results with surprising translational relevance. Murine OSCC cell lines were subjected to next-generation sequencing (NGS) to define their mutational landscape, to define novel candidate cancer genes, and to assess for parallels with known drivers in human OSCC. Expression arrays identified a mouse metastasis signature, and we assessed its representation in four independent human datasets comprising 324 patients using weighted voting and gene set enrichment analysis. Kaplan-Meier analysis and multivariate Cox proportional hazards modeling were used to stratify outcomes. A quantitative real-time PCR assay based on the mouse signature coupled to a machine-learning algorithm was developed and used to stratify an independent set of 31 patients with respect to metastatic lymphadenopathy. NGS revealed conservation of human driver pathway mutations in mouse OSCC, including in Trp53, mitogen-activated protein kinase, phosphoinositide 3-kinase, NOTCH, JAK/STAT, and Fat1-4. Moreover, comparative analysis between The Cancer Genome Atlas and mouse samples defined AKAP9, MED12L, and MYH6 as novel putative cancer genes. Expression analysis identified a transcriptional signature predicting aggressiveness and clinical outcomes, which were validated in four independent human OSCC datasets. Finally, we harnessed the translational potential of this signature by creating a clinically feasible assay that stratified patients with OSCC with a 93.5% accuracy. These data demonstrate surprising cross-species genomic conservation that has translational relevance for human oral squamous cell cancer. Clin Cancer Res; 20(11); 2873-84. ©2014 AACR.

Leveraging Cross-Species Transcription Factor Binding Site Patterns: From Diabetes Risk Loci to Disease Mechanisms

Genome-wide association studies have revealed numerous risk loci associated with diverse diseases. However, identification of disease-causing variants within association loci remains a major challenge. Divergence in gene expression due to cis-regulatory variants in noncoding regions is central to disease susceptibility. We show that integrative computational analysis of phylogenetic conservation with a complexity assessment of co-occurring transcription factor binding sites (TFBS) can identify cis-regulatory variants and elucidate their mechanistic role in disease. Analysis of established type 2 diabetes risk loci revealed a striking clustering of distinct homeobox TFBS. We identified the PRRX1 homeobox factor as a repressor of PPARG2 expression in adipose cells and demonstrate its adverse effect on lipid metabolism and systemic insulin sensitivity, dependent on the rs4684847 risk allele that triggers PRRX1 binding. Thus, cross-species conservation analysis at the level of co-occurring TFBS provides a valuable contribution to the translation of genetic association signals to disease-related molecular mechanisms.

Achilles tendon injuries in elite athletes: lessons in pathophysiology from their equine counterparts.

Superficial digital flexor tendon (SDFT) injury in equine athletes is one of the most well-accepted, scientifically supported companion animal models of human disease (i.e., exercise-induced Achilles tendon [AT] injury). The SDFT and AT are functionally and clinically equivalent (and important) energy-storing structures for which no equally appropriate rodent, rabbit, or other analogues exist. Access to equine tissues has facilitated significant advances in knowledge of tendon maturation and aging, determination of specific exercise effects (including early life), and definition of some of the earliest stages of subclinical pathology. Access to human surgical biopsies has provided complementary information on more advanced phases of disease. Importantly, equine SDFT injuries are only a model for acute ruptures in athletes, not the entire spectrum of human tendonopathy (including chronic tendon pain). In both, pathology begins with a potentially prolonged phase of accumulation of (subclinical) microdamage. Recent work has revealed remarkably similar genetic risk factors, including further evidence that tenocyte dysfunction plays an active role. Mice are convenient but not necessarily accurate models for multiple diseases, particularly at the cellular level. Mechanistic studies, including tendon cell responses to combinations of exercise-associated stresses, require a more thorough investigation of cross-species conservation of key stress pathway auditors. Molecular evidence has provided some context for the poor performance of mouse models; equines may provide better systems at this level. The use of horses may be additionally justifiable based on comparable species longevity, lifestyle factors, and selection pressure by similar infectious agents (e.g., herpesviruses) on general cell stress pathway evolution.

Death of p53-defective cells triggered by forced mitotic entry in the presence of DNA damage is not uniquely dependent on Caspase-2 or the PIDDosome

Much effort has been put in the discovery of ways to selectively kill p53-deficient tumor cells and targeting cell cycle checkpoint pathways has revealed promising candidates. Studies in zebrafish and human cell lines suggested that the DNA damage response kinase, checkpoint kinase 1 (Chk1), not only regulates onset of mitosis but also cell death in response to DNA damage in the absence of p53. This effect reportedly relies on ataxia telangiectasia mutated (ATM)-dependent and PIDDosome-mediated activation of Caspase-2. However, we show that genetic ablation of PIDDosome components in mice does not affect cell death in response to γ-irradiation. Furthermore, Chk1 inhibition largely failed to sensitize normal and malignant cells from p53−/− mice toward DNA damaging agents, and p53 status did not affect the death-inducing activity of DNA damage after Chk1 inhibition in human cancer cells. These observations argue against cross-species conservation of a Chk1-controlled cell survival pathway demanding further investigation of the molecular machinery responsible for cell death elicited by forced mitotic entry in the presence of DNA damage in different cell types and model organisms.

Cross-Species Withdrawal of MCL1 Facilitates Postpartum Uterine Involution in Both the Mouse and Baboon

A successful postpartum involution permits the postnatal uterus to rapidly regain its prepregnancy function and size to ultimately facilitate an ensuing blastocyst implantation. This study investigates the molecular mechanisms that govern the initiation of the involution process by examining the signaling events that occur as the uterus transitions from the pregnant to postnatal state. Using mouse and baboon uteri, we found a remarkable cross-species conservation at the signal transduction level as the pregnant uterus initiates and progresses through the involution process. This study originated with the observation of elevated levels of caspase-3 activation in both the laboring mouse and baboon uterus, which we found to be apoptotic in nature as evidenced by the concurrent appearance of cleaved poly(ADP-ribose) polymerase. We previously defined a nonapoptotic and potential tocolytic role for uterine caspase-3 during pregnancy regulated by increased antiapoptotic signaling mediated by myeloid cell leukemia sequence 1 and X-linked inhibitor of apoptosis. In contrast, this study determined that diminished antiapoptotic signaling in the postpartum uterus allowed for both endometrial apoptotic and myometrial autophagic episodes, which we speculate are responsible for the rapid reduction in size of the postpartum uterus. Using our human telomerase immortalized myometrial cell line and the Simian virus-40 immortalized endometrial cell line (12Z), we demonstrated that the withdrawal of antiapoptotic signaling was also an upstream event for both the autophagic and apoptotic processes in the human uterine myocyte and endometrial epithelial cell.

Proteome Mapping of Adult Zebrafish Marrow Neutrophils Reveals Partial Cross Species Conservation to Human Peripheral Neutrophils

Neutrophil granulocytes are pivotal cells within the first line of host defense of the innate immune system. In this study, we have used a gel-based LC-MS/MS approach to explore the proteome of primary marrow neutrophils from adult zebrafish. The identified proteins originated from all major cellular compartments. Gene ontology analysis revealed significant association of proteins with different immune-related network and pathway maps. 75% of proteins identified in neutrophils were identified in neutrophils only when compared to neutrophil-free brain tissue. Moreover, cross-species comparison with human peripheral blood neutrophils showed partial conservation of immune-related proteins between human and zebrafish. This study provides the first zebrafish neutrophil proteome and may serve as a valuable resource for an understanding of neutrophil biology and innate immunity.

Genome-Wide and Species-Wide In Silico Screening for Intragenic MicroRNAs in Human, Mouse and Chicken

MicroRNAs (miRNAs) are non-coding RNAs (ncRNAs) involved in regulation of gene expression. Intragenic miRNAs, especially those exhibiting a high degree of evolutionary conservation, have been shown to be coordinately regulated and/or expressed with their host genes, either with synergistic or antagonistic correlation patterns. However, the degree of cross-species conservation of miRNA/host gene co-location is not known and co-expression information is incomplete and fragmented among several studies. Using the genomic resources (miRBase and Ensembl) we performed a genome-wide in silico screening (GWISS) for miRNA/host gene pairs in three well-annotated vertebrate species: human, mouse, and chicken. Approximately half of currently annotated miRNA genes resided within host genes: 53.0% (849/1,600) in human, 48.8% (418/855) in mouse, and 42.0% (210/499) in chicken, which we present in a central publicly available Catalog of intragenic miRNAs (http://www.integratomics-time.com/miR-host/catalog). The miRNA genes resided within either protein-coding or ncRNA genes, which include long intergenic ncRNAs (lincRNAs) and small nucleolar RNAs (snoRNAs). Twenty-seven miRNA genes were found to be located within the same host genes in all three species and the data integration from literature and databases showed that most (26/27) have been found to be co-expressed. Particularly interesting are miRNA genes located within genes encoding for miRNA silencing machinery (DGCR8, DICER1, and SND1 in human and Cnot3, Gdcr8, Eif4e, Tnrc6b, and Xpo5 in mouse). We furthermore discuss a potential for phenotype misattribution of miRNA host gene polymorphism or gene modification studies due to possible collateral effects on miRNAs hosted within them. In conclusion, the catalog of intragenic miRNAs and identified 27 miRNA/host gene pairs with cross-species conserved co-location, co-expression, and potential co-regulation, provide excellent candidates for further functional annotation of intragenic miRNAs in health and disease.

GWAS3D: Detecting human regulatory variants by integrative analysis of genome-wide associations, chromosome interactions and histone modifications.

Interpreting the genetic variants located in the regulatory regions, such as enhancers and promoters, is an indispensable step to understand molecular mechanism of complex traits. Recent studies show that genetic variants detected by genome-wide association study (GWAS) are significantly enriched in the regulatory regions. Therefore, detecting, annotating and prioritizing of genetic variants affecting gene regulation are critical to our understanding of genotype–phenotype relationships. Here, we developed a web server GWAS3D to systematically analyze the genetic variants that could affect regulatory elements, by integrating annotations from cell type-specific chromatin states, epigenetic modifications, sequence motifs and cross-species conservation. The regulatory elements are inferred from the genome-wide chromosome interaction data, chromatin marks in 16 different cell types and 73 regulatory factors motifs from the Encyclopedia of DNA Element project. Furthermore, we used these function elements, as well as risk haplotype, binding affinity, conservation and P-values reported from the original GWAS to reprioritize the genetic variants. Using studies from low-density lipoprotein cholesterol, we demonstrated that our reprioritizing approach was effective and cell type specific. In conclusion, GWAS3D provides a comprehensive annotation and visualization tool to help users interpreting their results. The web server is freely available at http://jjwanglab.org/gwas3d.

Prediction of clustered RNA-binding protein motif sites in the mammalian genome.

Sequence-specific interactions of RNA-binding proteins (RBPs) with their target transcripts are essential for post-transcriptional gene expression regulation in mammals. However, accurate prediction of RBP motif sites has been difficult because many RBPs recognize short and degenerate sequences. Here we describe a hidden Markov model (HMM)-based algorithm mCarts to predict clustered functional RBP-binding sites by effectively integrating the number and spacing of individual motif sites, their accessibility in local RNA secondary structures and cross-species conservation. This algorithm learns and quantifies rules of these features, taking advantage of a large number of in vivo RBP-binding sites obtained from cross-linking and immunoprecipitation data. We applied this algorithm to study two representative RBP families, Nova and Mbnl, which regulate tissue-specific alternative splicing through interacting with clustered YCAY and YGCY elements, respectively, and predicted their binding sites in the mouse transcriptome. Despite the low information content in individual motif elements, our algorithm made specific predictions for successful experimental validation. Analysis of predicted sites also revealed cases of extensive and distal RBP-binding sites important for splicing regulation. This algorithm can be readily applied to other RBPs to infer their RNA-regulatory networks. The software is freely available at http://zhanglab.c2b2.columbia.edu/index.php/MCarts.

Progesterone receptor expression in the brain of the socially monogamous and paternal male prairie vole

Differences in the social organization and behavior of male mammals are attributable to species differences in neurochemistry, including differential expression of steroid hormone receptors. However, the distribution of progestin receptors (PR) in a socially monogamous and spontaneously parental male rodent has never been examined. Here we determined if PR exists and is regulated by testicular hormones in forebrain sites traditionally influencing socioreproductive behaviors in male prairie voles (Microtus ochrogaster). We hypothesized that PR expression in male prairie voles would differ from that described in other male rodents because PR activity inhibits parental behaviors and social memory in laboratory mice and rats. Adult male prairie voles received a sham surgery, were gonadectomized, or were gonadectomized and implanted with a testosterone-filled capsule. PR immunoreactivity (PRir) was measured four weeks later in areas of the hypothalamus and extended amygdala. A group of gonadally intact female prairie voles was included to reveal possible sex differences. We found considerable PRir in all sites examined. Castration reduced PRir in males' medial preoptic nucleus, anteroventral periventricular nucleus, ventromedial hypothalamus, and posterodorsal medial amygdala, and it was maintained in these sites by testosterone. This is the first study to examine PR expression in brain sites involved in socioreproductive behaviors in a socially monogamous and spontaneously paternal male rodent. Our results mostly reveal cross-species conservation in the distribution and hormone sensitivity of PR expression. Because PR interferes with aspects of sociality in other male rodents, PR may eventually be found to have different neurobiological actions in male prairie voles.

Use of the Protein Ontology for Multi-Faceted Analysis of Biological Processes: A Case Study of the Spindle Checkpoint

As a member of the Open Biomedical Ontologies (OBO) foundry, the Protein Ontology (PRO) provides an ontological representation of protein forms and complexes and their relationships. Annotations in PRO can be assigned to individual protein forms and complexes, each distinguishable down to the level of post-translational modification, thereby allowing for a more precise depiction of protein function than is possible with annotations to the gene as a whole. Moreover, PRO is fully interoperable with other OBO ontologies and integrates knowledge from other protein-centric resources such as UniProt and Reactome. Here we demonstrate the value of the PRO framework in the investigation of the spindle checkpoint, a highly conserved biological process that relies extensively on protein modification and protein complex formation. The spindle checkpoint maintains genomic integrity by monitoring the attachment of chromosomes to spindle microtubules and delaying cell cycle progression until the spindle is fully assembled. Using PRO in conjunction with other bioinformatics tools, we explored the cross-species conservation of spindle checkpoint proteins, including phosphorylated forms and complexes; studied the impact of phosphorylation on spindle checkpoint function; and examined the interactions of spindle checkpoint proteins with the kinetochore, the site of checkpoint activation. Our approach can be generalized to any biological process of interest.

APPRIS: annotation of principal and alternative splice isoforms

Here, we present APPRIS (http://appris.bioinfo.cnio.es), a database that houses annotations of human splice isoforms. APPRIS has been designed to provide value to manual annotations of the human genome by adding reliable protein structural and functional data and information from cross-species conservation. The visual representation of the annotations provided by APPRIS for each gene allows annotators and researchers alike to easily identify functional changes brought about by splicing events. In addition to collecting, integrating and analyzing reliable predictions of the effect of splicing events, APPRIS also selects a single reference sequence for each gene, here termed the principal isoform, based on the annotations of structure, function and conservation for each transcript. APPRIS identifies a principal isoform for 85% of the protein-coding genes in the GENCODE 7 release for ENSEMBL. Analysis of the APPRIS data shows that at least 70% of the alternative (non-principal) variants would lose important functional or structural information relative to the principal isoform.

CSF1R mutations identified in three families with autosomal dominantly inherited leukoencephalopathy

Genetic and phenotypic heterogeneities are considerably high in adult-onset leukoencephalopathy, in which comprehensive mutational analyses of the candidate genes by conventional methods are too laborious. We applied exome sequencing to conduct a comprehensive mutational analysis of genes for autosomal dominant leukoencephalopathies. Genomic DNA samples from four patients of three families with autosomal dominantly inherited adult-onset leukodystrophy were subjected to exome sequencing. On the basis of the results, 21 patients with adult-onset sporadic leukodystrophy and one patient with pathologically proven HDLS were additionally screened for CSF1R mutations. Exome sequencing identified heterozygous CSF1R mutations (p.I794T and p.R777W) in two families. I794T has recently been reported as a causative mutation for hereditary diffuse leukoencephalopathy with spheroids (HDLS), and R777W is a novel mutation. Although mutational analysis of CSF1R in 21 sporadic cases revealed no mutations, another novel CSF1R mutation, p.C653Y, was identified in one patient with autopsy-proven HDSL. These variants were located in the PTK domain where the causative mutations cluster. Functional prediction of the mutant CSF1R as well as cross-species conservation of the affected amino acids supports the notion that these variants are pathogenic for HDLS. Exome sequencing is useful for a comprehensive mutational analysis of causative genes for hereditary leukoencephalopathies, and CSF1R should be considered a candidate gene for patients with autosomal dominant leukoencephalopathies.