Generation of biological association networks: A novel strategy to detect new targets in cancer therapy

The aim of this work was to design a novel strategy to detect new targets for anticancer treatments. The rationale was to build Biological Association Networks from differentially expressed genes in drug-resistant cells to identify important nodes within the Networks. These nodes may represent putative targets to attack in cancer therapy, as a way to destabilize the gene network developed by the resistant cells to escape from the drug pressure. As a model we used cells resistant to methotrexate (MTX), an inhibitor of DHFR. Selected node-genes were analyzed at the transcriptional level and from a genotypic point of view. In colon cancer cells, DHFR, the AKR1 family, PKC¿, S100A4, DKK1, and CAV1 were overexpressed while E-cadherin was lost. In breast cancer cells, the UGT1A family was overexpressed, whereas EEF1A1 was overexpressed in pancreatic cells. Interference RNAs directed against these targets sensitized cells towards MTX.

Cellular source and mechanisms of high transcriptome complexity in the Mammalian testis.

Understanding the extent of genomic transcription and its functional relevance is a central goal in genomics research. However, detailed genome-wide investigations of transcriptome complexity in major mammalian organs have been scarce. Here, using extensive RNA-seq data, we show that transcription of the genome is substantially more widespread in the testis than in other organs across representative mammals. Furthermore, we reveal that meiotic spermatocytes and especially postmeiotic round spermatids have remarkably diverse transcriptomes, which explains the high transcriptome complexity of the testis as a whole. The widespread transcriptional activity in spermatocytes and spermatids encompasses protein-coding and long noncoding RNA genes but also poorly conserves intergenic sequences, suggesting that it may not be of immediate functional relevance. Rather, our analyses of genome-wide epigenetic data suggest that this prevalent transcription, which most likely promoted the birth of new genes during evolution, is facilitated by an overall permissive chromatin in these germ cells that results from extensive chromatin remodeling.

Circulating miRNAs: a new generation of anti-doping biomarkers

MicroRNAs (miRNAs) are small non-coding RNAs that regulate a variety of biological processes. Cell-free miRNAs detected in blood plasma are used as specific and sensitive markers of physiological processes and some diseases. Circulating miRNAs are highly stable in body fluids, for example plasma. Therefore, profiles of circulating miRNAs have been investigated for potential use as novel, non-invasive anti-doping biomarkers. This review describes the biological mechanisms underlying the variation of circulating miRNAs, revealing that they have great potential as a new class of biomarker for detection of doping substances. The latest developments in extraction and profiling technology, and the technical design of experiments useful for anti-doping, are also discussed. Longitudinal measurements of circulating miRNAs in the context of the athlete biological passport are proposed as an efficient strategy for the use of these new markers. The review also emphasizes potential challenges for the translation of circulating miRNAs from research into practical anti-doping applications.

Reduction of choroidal neovascularization in mice by adeno-associated virus-delivered anti-vascular endothelial growth factor short hairpin RNA.

BACKGROUND: Strategies leading to the long-term suppression of inappropriate ocular angiogenesis are required to avoid the need for repetitive monthly injections for treatment of diseases of the eye, such as age-related macular degeneration (AMD). The present study aimed to develop a strategy for the sustained repression of vascular endothelial growth factor (VEGF), which is identified as the key player in exudative AMD. METHODS: We have employed short hairpin (sh)RNAs combined with adeno-associated virus (AAV) delivery to obtain the targeted expression of potent gene-regulatory molecules. Anti-VEGF shRNAs were analyzed in human retinal pigment epithelial (RPE) cells using Renilla luciferase screening. For in vivo delivery of the most potent shRNA, self-complementary AAV vectors were packaged in serotype 8 capsids (scAAV2/8-hU6-sh9). In vivo efficacy was evaluated either by injection of scAAV2/8-hU6-sh9 into murine hind limb muscles or in a laser-induced murine model of choroidal neovascularization (CNV) following scAAV2/8-hU6-sh9 subretinal delivery. RESULTS: Plasmids encoding anti-VEGF shRNAs showed efficient knockdown of human VEGF in RPEs. Intramuscular administration led to localized expression and 91% knockdown of endogenous murine (m)VEGF. Subsequently, the ability of AAV2/8-encoded shRNAs to impair vessel formation was evaluated in the murine model of CNV. In this model, the sizes of the CNV were significantly reduced (up to 48%) following scAAV2/8-hU6-sh9 subretinal delivery. CONCLUSIONS: Using anti-VEGF vectors, we have demonstrated efficient silencing of endogenous mVEGF and showed that subretinal administration of scAAV2/8-hU6-sh9 has the ability to impair vessel formation in an AMD animal model. Thus, AAV-encoded shRNA can be used for the inhibition of neovascularization, leading to the development of sustained anti-VEGF therapy. Copyright © 2012 John Wiley & Sons, Ltd.

Regulation of antibiotic production in root-colonizing Peudomonas spp. and relevance for biological control of plant disease.

Certain strains of fluorescent pseudomonads are important biological components of agricultural soils that are suppressive to diseases caused by pathogenic fungi on crop plants. The biocontrol abilities of such strains depend essentially on aggressive root colonization, induction of systemic resistance in the plant, and the production of diffusible or volatile antifungal antibiotics. Evidence that these compounds are produced in situ is based on their chemical extraction from the rhizosphere and on the expression of antibiotic biosynthetic genes in the producer strains colonizing plant roots. Well-characterized antibiotics with biocontrol properties include phenazines, 2,4-diacetylphloroglucinol, pyoluteorin, pyrrolnitrin, lipopeptides, and hydrogen cyanide. In vitro, optimal production of these compounds occurs at high cell densities and during conditions of restricted growth, involving (i) a number of transcriptional regulators, which are mostly pathway-specific, and (ii) the GacS/GacA two-component system, which globally exerts a positive effect on the production of extracellular metabolites at a posttranscriptional level. Small untranslated RNAs have important roles in the GacS/GacA signal transduction pathway. One challenge in future biocontrol research involves development of new strategies to overcome the broad toxicity and lack of antifungal specificity displayed by most biocontrol antibiotics studied so far.

Cognate microRNA-Gene Regulatory Networks Mediating Glioblastoma Oncogenic Properties

Malignant gliomas, including the most common and fatal form glioblastoma (GBM, WHO grade IV astrocytoma), remain a challenge to treat. In the United States and Europe, more than 30,000 patients per year are newly diagnosed with GBM. Despite ongoing trials, the best currently available multimodal treatment approaches include surgical resection followed by concomitant and adjuvant radiation (RT) and temozolomide (TMZ) therapy, resulting in a low median overall survival (OS) rate ranging from 12.2 - 15.9 months. The important role of genetic and epigenetic changes in DNA, RNA, and protein alteration as well as epigenetic changes secondary to the tumor microenvironment and outside selection pressure (therapeutic interventions), are increasingly being recognized. In GBM treatment, the focus is shifting toward a more patient-centered (personalized) therapy. In this regard, in particular, microRNAs are being increasingly studied. MicroRNAs are non¬protein coding small RNAs that serve as negative gene regulators by binding to a specific sequence in the promoter region of a target gene, thus regulating gene expression. A single microRNA potentially targets hundreds of genes; thus, microRNAs and their cognate target genes have important roles as tumor suppressors and oncogenes as well as regulators of various cancer- specific cellular features, such as proliferation, apoptosis, invasion, and metastasis. The identification of distinct microRNA-gene regulatory networks in GBM patients can be expected to provide novel therapeutic insights by identifying candidate patients for targeted therapies. To this end, in this work we identified and validated clinically relevant and meaningful novel gene- microRNA regulatory networks that correlated with MR tumor phenotypes, histopathology, and patient survival and response rates to therapy. - Le traitement des gliomes malins, y compris sous leur forme la plus commune et meurtrière, le glioblastome (GBM, ou astrocytome de grade IV selon l'OMS), demeure à ce jour un défi. Aux États-Unis et en Europe, un nouveau diagnostic de GBM est prononcé dans plus de 30Ό00 cas par an. En dépit de tests en cours, les meilleures approches thérapeutiques combinées actuellement disponibles comprennent la résection chirurgicale de la tumeur, suivie d'une radiothérapie adjuvante ainsi que d'un traitement au temozolomide (RT/TMZ), thérapies dont résulte une médiane de survie globale basse (overall survival, OS), comprise entre 12.2 et 15.9 mois. On reconnaît de plus en plus le rôle majeur de l'ADN, de l'ARN et de l'altération des protéines ainsi que des modifications épigénétiques, secondaires par rapport au microenvironnement de la tumeur et à la pression de sélection extérieure (les interventions thérapeutiques). Dans le traitement du GBM, le centre d'intérêt se déplace vers une thérapie centrée sur le cas individuel du patient. Dans ce but, en particulier les microARN sont de plus en plus analysés. Les microARN sont de petits ARN non-codants (les protéines) qui servent de régulateurs négatifs de gènes en s'attachant à une séquence spécifique dans la région promotrice d'un gène-cible, régulant ainsi l'expression du gène. Un seul microARN cible potentiellement des centaines de gènes; on a ainsi découvert que les microARN et leurs gènes-cibles apparentés ont une fonction importante en tant que suppresseurs de tumeurs et d'oncogènes, ainsi que comme régulateurs de diverses caractéristiques cellulaires spécifiques du cancer, comme la prolifération, l'apoptose, l'invasion et la métastase. On peut s'attendre à ce que l'identification de réseaux microARN régulateurs de gènes, distincts selon les patients de GBM, fournisse une approche thérapeutique inédite par la détermination des patients susceptibles de réagir favorablement à des thérapies ciblées.

Networking of differentially expressed genes in human cancer cells resistant to methotrexate

BACKGROUND: The need for an integrated view of data obtained from high-throughput technologies gave rise to network analyses. These are especially useful to rationalize how external perturbations propagate through the expression of genes. To address this issue in the case of drug resistance, we constructed biological association networks of genes differentially expressed in cell lines resistant to methotrexate (MTX). METHODS: Seven cell lines representative of different types of cancer, including colon cancer (HT29 and Caco2), breast cancer (MCF-7 and MDA-MB-468), pancreatic cancer (MIA PaCa-2), erythroblastic leukemia (K562) and osteosarcoma (Saos-2), were used. The differential expression pattern between sensitive and MTX-resistant cells was determined by whole human genome microarrays and analyzed with the GeneSpring GX software package. Genes deregulated in common between the different cancer cell lines served to generate biological association networks using the Pathway Architect software. RESULTS: Dikkopf homolog-1 (DKK1) is a highly interconnected node in the network generated with genes in common between the two colon cancer cell lines, and functional validations of this target using small interfering RNAs (siRNAs) showed a chemosensitization toward MTX. Members of the UDP-glucuronosyltransferase 1A (UGT1A) family formed a network of genes differentially expressed in the two breast cancer cell lines. siRNA treatment against UGT1A also showed an increase in MTX sensitivity. Eukaryotic translation elongation factor 1 alpha 1 (EEF1A1) was overexpressed among the pancreatic cancer, leukemia and osteosarcoma cell lines, and siRNA treatment against EEF1A1 produced a chemosensitization toward MTX. CONCLUSIONS: Biological association networks identified DKK1, UGT1As and EEF1A1 as important gene nodes in MTX-resistance. Treatments using siRNA technology against these three genes showed chemosensitization toward MTX.

Allele-specific silencing of mutant huntingtin in rodent brain and human stem cells.

Huntington's disease (HD) is an autosomal dominant neurodegenerative disorder resulting from polyglutamine expansion in the huntingtin (HTT) protein and for which there is no cure. Although suppression of both wild type and mutant HTT expression by RNA interference is a promising therapeutic strategy, a selective silencing of mutant HTT represents the safest approach preserving WT HTT expression and functions. We developed small hairpin RNAs (shRNAs) targeting single nucleotide polymorphisms (SNP) present in the HTT gene to selectively target the disease HTT isoform. Most of these shRNAs silenced, efficiently and selectively, mutant HTT in vitro. Lentiviral-mediated infection with the shRNAs led to selective degradation of mutant HTT mRNA and prevented the apparition of neuropathology in HD rat's striatum expressing mutant HTT containing the various SNPs. In transgenic BACHD mice, the mutant HTT allele was also silenced by this approach, further demonstrating the potential for allele-specific silencing. Finally, the allele-specific silencing of mutant HTT in human embryonic stem cells was accompanied by functional recovery of the vesicular transport of BDNF along microtubules. These findings provide evidence of the therapeutic potential of allele-specific RNA interference for HD.

Combining RT-PCR-seq and RNA-seq to catalog all genic elements encoded in the human genome.

Within the ENCODE Consortium, GENCODE aimed to accurately annotate all protein-coding genes, pseudogenes, and noncoding transcribed loci in the human genome through manual curation and computational methods. Annotated transcript structures were assessed, and less well-supported loci were systematically, experimentally validated. Predicted exon-exon junctions were evaluated by RT-PCR amplification followed by highly multiplexed sequencing readout, a method we called RT-PCR-seq. Seventy-nine percent of all assessed junctions are confirmed by this evaluation procedure, demonstrating the high quality of the GENCODE gene set. RT-PCR-seq was also efficient to screen gene models predicted using the Human Body Map (HBM) RNA-seq data. We validated 73% of these predictions, thus confirming 1168 novel genes, mostly noncoding, which will further complement the GENCODE annotation. Our novel experimental validation pipeline is extremely sensitive, far more than unbiased transcriptome profiling through RNA sequencing, which is becoming the norm. For example, exon-exon junctions unique to GENCODE annotated transcripts are five times more likely to be corroborated with our targeted approach than with extensive large human transcriptome profiling. Data sets such as the HBM and ENCODE RNA-seq data fail sampling of low-expressed transcripts. Our RT-PCR-seq targeted approach also has the advantage of identifying novel exons of known genes, as we discovered unannotated exons in ~11% of assessed introns. We thus estimate that at least 18% of known loci have yet-unannotated exons. Our work demonstrates that the cataloging of all of the genic elements encoded in the human genome will necessitate a coordinated effort between unbiased and targeted approaches, like RNA-seq and RT-PCR-seq.

Mechanism of hcnA mRNA recognition in the Gac/Rsm signal transduction pathway of Pseudomonas fluorescens.

In the plant-beneficial bacterium Pseudomonas fluorescens CHA0, the expression of antifungal exoproducts is controlled by the GacS/GacA two-component system. Two RNA binding proteins (RsmA, RsmE) ensure effective translational repression of exoproduct mRNAs. At high cell population densities, GacA induces three small RNAs (RsmX, RsmY, RsmZ) which sequester both RsmA and RsmE, thereby relieving translational repression. Here we systematically analyse the features that allow the RNA binding proteins to interact strongly with the 5' untranslated leader mRNA of the P. fluorescens hcnA gene (encoding hydrogen cyanide synthase subunit A). We obtained evidence for three major RsmA/RsmE recognition elements in the hcnA leader, based on directed mutagenesis, RsmE footprints and toeprints, and in vivo expression data. Two recognition elements were found in two stem-loop structures whose existence in the 5' leader region was confirmed by lead(II) cleavage analysis. The third recognition element, which overlapped the hcnA Shine-Dalgarno sequence, was postulated to adopt either an open conformation, which would favour ribosome binding, or a stem-loop structure, which may form upon interaction with RsmA/RsmE and would inhibit access of ribosomes. Effective control of hcnA expression by the Gac/Rsm system appears to result from the combination of the three appropriately spaced recognition elements.

MiR-SNPs as markers of toxicity and clinical outcome in Hodgkin Lymphoma patients

Background: In recent years, microRNA (miRNA) pathways have emerged as a crucial system for the regulation of tumorogenesis. miR-SNPs are a novel class of single nucleotide polymorphisms that can affect miRNA pathways. Design and Methods: We analyzed eight miR-SNPs by allelic discrimination in 141 patients with Hodgkin lymphoma and correlated the results with treatment-related toxicity, response, disease-free survival (DFS) and overall survival (OS). Results: The KRT81 (rs3660) GG genotype was associated with an increased risk of neurological toxicity (P=0.016), while patients with XPO5 (rs11077) AA or CC genotypes had a higher rate of bleomycin-associated pulmonary toxicity (P=0.048). Both miR-SNPs emerged as independent factors in the multivariate analysis. The XPO5 AA and CC genotypes were also associated with a lower response rate (P=0.036). XPO5 (P=0.039) and TRBP (rs784567) (P=0.022) genotypes emerged as prognostic markers for DFS, and XPO5 was also associated with OS (P=0.033). In the multivariate analysis, only XPO5 emerged as an independent prognostic factor for DFS (HR: 2.622; 95%CI 1.039-6.620; P=0.041). Given the influence of XPO5 and TRBP as individual markers, we then investigated the combined effect of these miR-SNPs. Patients with both the XPO5 AA/CC and TRBP TT/TC genotypes had the shortest DFS (P=0.008) and OS (P=0.008). Conclusion: miR-SNPs can add useful prognostic information on treatment-related toxicity and clinical outcome in Hodgkin lymphoma and can be used to identify patients likely to be chemoresistant or to relapse.

Architectural features of bacterial genomes and their applications

Résumé -Caractéristiques architecturales des génomes bactériens et leurs applications Les bactéries possèdent généralement un seul chromosome circulaire. A chaque génération, ce chromosome est répliqué bidirectionnellement, par deux complexes enzymatiques de réplication se déplaçant en sens opposé depuis l'origine de réplication jusqu'au terminus, situé à l'opposé. Ce mode de réplication régit l'architecture du chromosome -l'orientation des gènes par rapport à la réplication, notamment - et est en grande partie à l'origine des pressions qui provoquent la variation de la composition en nucléotides du génome, hors des contraintes liées à la structure et à la fonction des protéines codées sur le chromosome. Le but de cette thèse est de contribuer à quantifier les effets de la réplication sur l'architecture chromosomique, en s'intéressant notamment aux gènes des ARN ribosomiques, cruciaux pour la bactérie. D'un autre côté, cette architecture est spécifique à l'espèce et donne ainsi une «identité génomique » aux gènes. Il est démontré ici qu'il est possible d'utiliser des marqueurs «naïfs » de cette identité pour détecter, notamment dans le génome du staphylocoque doré, des îlots de pathogénicité, qui concentrent un grand nombre de facteurs de virulence de la bactérie. Ces îlots de pathogénicité sont mobiles, et peuvent passer d'une bactérie à une autre, mais conservent durant un certain temps l'identité génomique de leur hôte précédent, ce qui permet de les reconnaître dans leur nouvel hôte. Ces méthodes simples, rapides et fiables seront de la plus haute importance lorsque le séquençage des génomes entiers sera rapide et disponible à très faible coût. Il sera alors possible d'analyser instantanément les déterminants pathogéniques et de résistance aux antibiotiques des agents pathogènes. Summary The bacterial genome is a highly organized structure, which may be referred to as the genome architecture, and is mainly directed by DNA replication. This thesis provides significant insights in the comprehension of the forces that shape bacterial chromosomes, different in each genome and contributing to confer them an identity. First, it shows the importance of the replication in directing the orientation of prokaryotic ribosomal RNAs, and how it shapes their nucleotide composition in a tax on-specific manner. Second, it highlights the pressure acting on the orientation of the genes in general, a majority of which are transcribed in the same direction as replication. Consequently, apparent infra-arm genome rearrangements, involving an exchange of the leading/lagging strands and shown to reduce growth rate, are very likely artifacts due to an incorrect contig assembly. Third, it shows that this genomic identity can be used to detect foreign parts in genomes, by establishing this identity for a given host and identifying the regions that deviate from it. This property is notably illustrated with Staphylococcus aureus: known pathogenicity islands and phages, and putative ancient pathogenicity islands concentrating many known pathogenicity-related genes are highlighted; the analysis also detects, incidentally, proteins responsible for the adhesion of S. aureus to the hosts' cells. In conclusion, the study of nucleotide composition of bacterial genomes provides the opportunity to better understand the genome-level pressures that shape DNA sequences, and to identify genes and regions potentially related to pathogenicity with fast, simple and reliable methods. This will be of crucial importance when whole-genome sequencing will be a rapid, inexpensive and routine tool.

Circulating microRNAs as novel biomarkers for diabetes mellitus.

Diabetes mellitus is characterized by insulin secretion from pancreatic β cells that is insufficient to maintain blood glucose homeostasis. Autoimmune destruction of β cells results in type 1 diabetes mellitus, whereas conditions that reduce insulin sensitivity and negatively affect β-cell activities result in type 2 diabetes mellitus. Without proper management, patients with diabetes mellitus develop serious complications that reduce their quality of life and life expectancy. Biomarkers for early detection of the disease and identification of individuals at risk of developing complications would greatly improve the care of these patients. Small non-coding RNAs called microRNAs (miRNAs) control gene expression and participate in many physiopathological processes. Hundreds of miRNAs are actively or passively released in the circulation and can be used to evaluate health status and disease progression. Both type 1 diabetes mellitus and type 2 diabetes mellitus are associated with distinct modifications in the profile of miRNAs in the blood, which are sometimes detectable several years before the disease manifests. Moreover, circulating levels of certain miRNAs seem to be predictive of long-term complications. Technical and scientific obstacles still exist that need to be overcome, but circulating miRNAs might soon become part of the diagnostic arsenal to identify individuals at risk of developing diabetes mellitus and its devastating complications.

Results of An International Study of Quantitative BCR/ABL Assays Using Defined RNA Samples

Molecular monitoring of BCR/ABL transcripts by real time quantitative reverse transcription PCR (qRT-PCR) is an essential technique for clinical management of patients with BCR/ABL-positive CML and ALL. Though quantitative BCR/ABL assays are performed in hundreds of laboratories worldwide, results among these laboratories cannot be reliably compared due to heterogeneity in test methods, data analysis, reporting, and lack of quantitative standards. Recent efforts towards standardization have been limited in scope. Aliquots of RNA were sent to clinical test centers worldwide in order to evaluate methods and reporting for e1a2, b2a2, and b3a2 transcript levels using their own qRT-PCR assays. Total RNA was isolated from tissue culture cells that expressed each of the different BCR/ABL transcripts. Serial log dilutions were prepared, ranging from 100 to 10-5, in RNA isolated from HL60 cells. Laboratories performed 5 independent qRT-PCR reactions for each sample type at each dilution. In addition, 15 qRT-PCR reactions of the 10-3 b3a2 RNA dilution were run to assess reproducibility within and between laboratories. Participants were asked to run the samples following their standard protocols and to report cycle threshold (Ct), quantitative values for BCR/ABL and housekeeping genes, and ratios of BCR/ABL to housekeeping genes for each sample RNA. Thirty-seven (n=37) participants have submitted qRT-PCR results for analysis (36, 37, and 34 labs generated data for b2a2, b3a2, and e1a2, respectively). The limit of detection for this study was defined as the lowest dilution that a Ct value could be detected for all 5 replicates. For b2a2, 15, 16, 4, and 1 lab(s) showed a limit of detection at the 10-5, 10-4, 10-3, and 10-2 dilutions, respectively. For b3a2, 20, 13, and 4 labs showed a limit of detection at the 10-5, 10-4, and 10-3 dilutions, respectively. For e1a2, 10, 21, 2, and 1 lab(s) showed a limit of detection at the 10-5, 10-4, 10-3, and 10-2 dilutions, respectively. Log %BCR/ABL ratio values provided a method for comparing results between the different laboratories for each BCR/ABL dilution series. Linear regression analysis revealed concordance among the majority of participant data over the 10-1 to 10-4 dilutions. The overall slope values showed comparable results among the majority of b2a2 (mean=0.939; median=0.9627; range (0.399 - 1.1872)), b3a2 (mean=0.925; median=0.922; range (0.625 - 1.140)), and e1a2 (mean=0.897; median=0.909; range (0.5174 - 1.138)) laboratory results (Fig. 1-3)). Thirty-four (n=34) out of the 37 laboratories reported Ct values for all 15 replicates and only those with a complete data set were included in the inter-lab calculations. Eleven laboratories either did not report their copy number data or used other reporting units such as nanograms or cell numbers; therefore, only 26 laboratories were included in the overall analysis of copy numbers. The median copy number was 348.4, with a range from 15.6 to 547,000 copies (approximately a 4.5 log difference); the median intra-lab %CV was 19.2% with a range from 4.2% to 82.6%. While our international performance evaluation using serially diluted RNA samples has reinforced the fact that heterogeneity exists among clinical laboratories, it has also demonstrated that performance within a laboratory is overall very consistent. Accordingly, the availability of defined BCR/ABL RNAs may facilitate the validation of all phases of quantitative BCR/ABL analysis and may be extremely useful as a tool for monitoring assay performance. Ongoing analyses of these materials, along with the development of additional control materials, may solidify consensus around their application in routine laboratory testing and possible integration in worldwide efforts to standardize quantitative BCR/ABL testing.

The origins, evolution, and functional potential of alternative splicing in vertebrates.

Alternative splicing (AS) has the potential to greatly expand the functional repertoire of mammalian transcriptomes. However, few variant transcripts have been characterized functionally, making it difficult to assess the contribution of AS to the generation of phenotypic complexity and to study the evolution of splicing patterns. We have compared the AS of 309 protein-coding genes in the human ENCODE pilot regions against their mouse orthologs in unprecedented detail, utilizing traditional transcriptomic and RNAseq data. The conservation status of every transcript has been investigated, and each functionally categorized as coding (separated into coding sequence [CDS] or nonsense-mediated decay [NMD] linked) or noncoding. In total, 36.7% of human and 19.3% of mouse coding transcripts are species specific, and we observe a 3.6 times excess of human NMD transcripts compared with mouse; in contrast to previous studies, the majority of species-specific AS is unlinked to transposable elements. We observe one conserved CDS variant and one conserved NMD variant per 2.3 and 11.4 genes, respectively. Subsequently, we identify and characterize equivalent AS patterns for 22.9% of these CDS or NMD-linked events in nonmammalian vertebrate genomes, and our data indicate that functional NMD-linked AS is more widespread and ancient than previously thought. Furthermore, although we observe an association between conserved AS and elevated sequence conservation, as previously reported, we emphasize that 30% of conserved AS exons display sequence conservation below the average score for constitutive exons. In conclusion, we demonstrate the value of detailed comparative annotation in generating a comprehensive set of AS transcripts, increasing our understanding of AS evolution in vertebrates. Our data supports a model whereby the acquisition of functional AS has occurred throughout vertebrate evolution and is considered alongside amino acid change as a key mechanism in gene evolution.