Search

Purpose: HLA-A29 is the main susceptibility factor for birdshot chorioretinitis (BSCR). Our study assessed the impact of the second HLA-A allele alongside HLA-A29 on BSCR severity and susceptibility, focusing on HLA-A29 homozygous patients and those with alleles from the HLA-Aw19 group., Methods: We included 120 additional cases to our previous analysis of 286 patients with BSCR, all HLA-A29 positive. Patients were categorized based on the second allele being also HLA-A29 (A29/nonA29 vs. A29/A29) or belonging to the HLA-Aw19 family, including A29, A30, A31, and A33 (A29/nonAw19 vs. A29/Aw19). HLA-A32 was analyzed separately (A29/nonA32 vs. A29/A32). The prevalence of these groups among patients with BSCR was compared with their frequencies in a sample of 151,997 French subjects. Disease severity was approximated by assessing disease onset and visual function at the last visit and was compared between patient groups., Results: When comparing the HLA-A29-positive patients with BSCR to HLA-A29-positive French subjects, we found an overrepresentation of HLA-A29 for the second HLA-A allele (χ² = 4.34; P = 0.037; odds ratio, 1.57; confidence interval, 1.01-2.44). Within the HLA-Aw19 broad antigen family, HLA-A32 was found to be under-represented (χ² = 6.15; P = 0.013; odds ratio, 0.40; confidence interval, 0.19-0.85). The nature of the second HLA-A allele did not impact disease severity., Conclusions: Homozygosity for HLA-A29 increased the risk of developing BSCR without affecting disease severity. The under-representation of HLA-A32 in patients with BSCR suggests a potential protective role.

Microglia are the resident immune cells of the central nervous system and are involved in brain development, homeostasis, and disease. New imaging and genomics technologies are revealing microglial complexity across developmental and functional states, brain regions, and diseases. We curated a set of publicly available gene expression datasets from human microglia spanning disease and health to identify sets of genes reflecting physiological and pathological microglial states. We also integrated multiple human microglial single-cell RNA-seq datasets in Alzheimer's disease (AD), multiple sclerosis (MS), and Parkinson's disease, and identified a distinct microglial transcriptional signature shared across diseases. Analysis of germ-line DNA identified genes with variants associated with AD and MS that are overrepresented in microglial gene sets, including the disease-associated transcriptional signature. This work points to genes that are dysregulated in disease states and provides a resource for the analysis of diseases in which microglia are implicated by genetic evidence., (© 2024. The Author(s).)

The genetic factors of stroke in South Asians are largely unexplored. Exome-wide sequencing and association analysis (ExWAS) in 75 K Pakistanis identified NM_000435.3(NOTCH3):c.3691 C > T, encoding the missense amino acid substitution p.Arg1231Cys, enriched in South Asians (alternate allele frequency = 0.58% compared to 0.019% in Western Europeans), and associated with subcortical hemorrhagic stroke [odds ratio (OR) = 3.39, 95% confidence interval (CI) = [2.26, 5.10], p = 3.87 × 10 -9 ), and all strokes (OR [CI] = 2.30 [1.77, 3.01], p = 7.79 × 10 -10 ). NOTCH3 p.Arg231Cys was strongly associated with white matter hyperintensity on MRI in United Kingdom Biobank (UKB) participants (effect [95% CI] in SD units = 1.1 [0.61, 1.5], p = 3.0 × 10 -6 ). The variant is attributable for approximately 2.0% of hemorrhagic strokes and 1.1% of all strokes in South Asians. These findings highlight the value of diversity in genetic studies and have major implications for genomic medicine and therapeutic development in South Asian populations., (© 2024. The Author(s).)

Rare coding variants that substantially affect function provide insights into the biology of a gene 1-3 . However, ascertaining the frequency of such variants requires large sample sizes 4-8 . Here we present a catalogue of human protein-coding variation, derived from exome sequencing of 983,578 individuals across diverse populations. In total, 23% of the Regeneron Genetics Center Million Exome (RGC-ME) data come from individuals of African, East Asian, Indigenous American, Middle Eastern and South Asian ancestry. The catalogue includes more than 10.4 million missense and 1.1 million predicted loss-of-function (pLOF) variants. We identify individuals with rare biallelic pLOF variants in 4,848 genes, 1,751 of which have not been previously reported. From precise quantitative estimates of selection against heterozygous loss of function (LOF), we identify 3,988 LOF-intolerant genes, including 86 that were previously assessed as tolerant and 1,153 that lack established disease annotation. We also define regions of missense depletion at high resolution. Notably, 1,482 genes have regions that are depleted of missense variants despite being tolerant of pLOF variants. Finally, we estimate that 3% of individuals have a clinically actionable genetic variant, and that 11,773 variants reported in ClinVar with unknown significance are likely to be deleterious cryptic splice sites. To facilitate variant interpretation and genetics-informed precision medicine, we make this resource of coding variation from the RGC-ME dataset publicly accessible through a variant allele frequency browser., (© 2024. The Author(s).)

Anterior Uveitis (AU) is the inflammation of the anterior part of the eye, the iris and ciliary body and is strongly associated with HLA-B*27. We report AU exome sequencing results from eight independent cohorts consisting of 3,850 cases and 916,549 controls. We identify common genome-wide significant loci in HLA-B (OR = 3.37, p = 1.03e-196) and ERAP1 (OR = 0.86, p = 1.1e-08), and find IPMK (OR = 9.4, p = 4.42e-09) and IDO2 (OR = 3.61, p = 6.16e-08) as genome-wide significant genes based on the burden of rare coding variants. Dividing the cohort into HLA-B*27 positive and negative individuals, we find ERAP1 haplotype is strongly protective only for B*27-positive AU (OR = 0.73, p = 5.2e-10). Investigation of B*27-negative AU identifies a common signal near HLA-DPB1 (rs3117230, OR = 1.26, p = 2.7e-08), risk genes IPMK and IDO2, and several additional candidate risk genes, including ADGFR5, STXBP2, and ACHE. Taken together, we decipher the genetics underlying B*27-positive and -negative AU and identify rare and common genetic signals for both subtypes of disease., (© 2023. The Author(s).)

Human genetic studies of smoking behavior have been thus far largely limited to common variants. Studying rare coding variants has the potential to identify drug targets. We performed an exome-wide association study of smoking phenotypes in up to 749,459 individuals and discovered a protective association in CHRNB2, encoding the β2 subunit of the α4β2 nicotine acetylcholine receptor. Rare predicted loss-of-function and likely deleterious missense variants in CHRNB2 in aggregate were associated with a 35% decreased odds for smoking heavily (odds ratio (OR) = 0.65, confidence interval (CI) = 0.56-0.76, P = 1.9 × 10 -8 ). An independent common variant association in the protective direction ( rs2072659 ; OR = 0.96; CI = 0.94-0.98; P = 5.3 × 10 -6 ) was also evident, suggesting an allelic series. Our findings in humans align with decades-old experimental observations in mice that β2 loss abolishes nicotine-mediated neuronal responses and attenuates nicotine self-administration. Our genetic discovery will inspire future drug designs targeting CHRNB2 in the brain for the treatment of nicotine addiction., (© 2023. The Author(s).)

De novo mutations in GNB1 , encoding the G β 1 subunit of G proteins, cause a neurodevelopmental disorder with global developmental delay and epilepsy, GNB1 encephalopathy. Here, we show that mice carrying a pathogenic mutation, K78R, recapitulate aspects of the disorder, including developmental delay and generalized seizures. Cultured mutant cortical neurons also display aberrant bursting activity on multi-electrode arrays. Strikingly, the antiepileptic drug ethosuximide (ETX) restores normal neuronal network behavior in vitro and suppresses spike-and-wave discharges (SWD) in vivo . ETX is a known blocker of T-type voltage-gated Ca 2+ channels and G protein-coupled potassium (GIRK) channels. Accordingly, we present evidence that K78R results in a gain-of-function (GoF) effect by increasing the activation of GIRK channels in cultured neurons and a heterologous model ( Xenopus oocytes)-an effect we show can be potently inhibited by ETX. This work implicates a GoF mechanism for GIRK channels in epilepsy, identifies a new mechanism of action for ETX in preventing seizures, and establishes this mouse model as a pre-clinical tool for translational research with predicative value for GNB1 encephalopathy., Competing Interests: DG is a founder of Actio Biosciences, a founder of and holds equity in Praxis, serves as a consultant to AstraZeneca, and has received research support from Janssen, Gilead, Biogen, AstraZeneca and UCB. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Colombo, Reddy, Petri, Williams, Shalomov, Dhindsa, Gelfman, Krizay, Bera, Yang, Peng, Makinson, Boland, Frankel, Goldstein and Dascal.)

Schizophrenia has a multifactorial etiology, involving a polygenic architecture. The potential benefit of whole genome sequencing (WGS) in schizophrenia and other psychotic disorders is not well studied. We investigated the yield of clinical WGS analysis in 251 families with a proband diagnosed with schizophrenia (N = 190), schizoaffective disorder (N = 49), or other conditions involving psychosis (N = 48). Participants were recruited in Israel and USA, mainly of Jewish, Arab, and other European ancestries. Trio (parents and proband) WGS was performed for 228 families (90.8%); in the other families, WGS included parents and at least two affected siblings. In the secondary analyses, we evaluated the contribution of rare variant enrichment in particular gene sets, and calculated polygenic risk score (PRS) for schizophrenia. For the primary outcome, diagnostic rate was 6.4%; we found clinically significant, single nucleotide variants (SNVs) or small insertions or deletions (indels) in 14 probands (5.6%), and copy number variants (CNVs) in 2 (0.8%). Significant enrichment of rare loss-of-function variants was observed in a gene set of top schizophrenia candidate genes in affected individuals, compared with population controls (N = 6,840). The PRS for schizophrenia was significantly increased in the affected individuals group, compared to their unaffected relatives. Last, we were also able to provide pharmacogenomics information based on CYP2D6 genotype data for most participants, and determine their antipsychotic metabolizer status. In conclusion, our findings suggest that WGS may have a role in the setting of both research and genetic counseling for individuals with schizophrenia and other psychotic disorders and their families., (© 2021. The Author(s), under exclusive licence to Springer Nature Limited.)

Alopecia areata is a complex genetic disease that results in hair loss due to the autoimmune-mediated attack of the hair follicle. We previously defined a role for both rare and common variants in our earlier GWAS and linkage studies. Here, we identify rare variants contributing to Alopecia Areata using a whole exome sequencing and gene-level burden analyses approach on 849 Alopecia Areata patients compared to 15,640 controls. KRT82 is identified as an Alopecia Areata risk gene with rare damaging variants in 51 heterozygous Alopecia Areata individuals (6.01%), achieving genome-wide significance (p = 2.18E-07). KRT82 encodes a hair-specific type II keratin that is exclusively expressed in the hair shaft cuticle during anagen phase, and its expression is decreased in Alopecia Areata patient skin and hair follicles. Finally, we find that cases with an identified damaging KRT82 variant and reduced KRT82 expression have elevated perifollicular CD8 infiltrates. In this work, we utilize whole exome sequencing to successfully identify a significant Alopecia Areata disease-relevant gene, KRT82, and reveal a proposed mechanism for rare variant predisposition leading to disrupted hair shaft integrity., (© 2022. The Author(s).)

Gene set-based signal detection analyses are used to detect an association between a trait and a set of genes by accumulating signals across the genes in the gene set. Since signal detection is concerned with identifying whether any of the genes in the gene set are non-null, a goodness-of-fit (GOF) test can be used to compare whether the observed distribution of gene-level tests within the gene set agrees with the theoretical null distribution. Here, we present a flexible gene set-based signal detection framework based on tail-focused GOF statistics. We show that the power of the various statistics in this framework depends critically on two parameters: the proportion of genes within the gene set that are non-null and the degree of separation between the null and alternative distributions of the gene-level tests. We give guidance on which statistic to choose for a given situation and implement the methods in a fast and user-friendly R package, wHC (https://github.com/mqzhanglab/wHC). Finally, we apply these methods to a whole exome sequencing study of amyotrophic lateral sclerosis., (© The Author(s) 2021. Published by Oxford University Press. All rights reserved. For Permissions, please email: journals.permissions@oup.com.)

Purpose: Birdshot chorioretinopathy (BSCR) is strongly associated with HLA-A29. This study was designed to elucidate the genetic modifiers of BSCR in HLA-A29 carriers., Methods: We sequenced the largest BSCR cohort to date, including 286 cases and 108 HLA-A29-positive controls to determine genome-wide common and rare variant associations. We further typed the HLA alleles of cases and 45,386 HLA-A29 controls of European ancestry to identify HLA alleles that associate with BSCR risk., Results: Carrying a second allele that belongs to the HLA-Aw19 broad antigen family (including HLA-A29, -A30, -A31, and -A33) increases the risk for BSCR (odds ratio [OR] = 4.44; P = 2.2e-03). This result was validated by comparing allele frequencies to large HLA-A29-controlled cohorts (n = 45,386; OR > 2.5; P < 1.3e-06). We also confirm that ERAP1 and ERAP2 haplotypes modulate disease risk. A meta-analysis with an independent dataset confirmed that ERAP1 and ERAP2 haplotypes modulate the risk for disease at a genome-wide significant level: ERAP1-rs27432 (OR = 2.46; 95% confidence interval [CI], 1.85-3.26; P = 4.07e-10), an expression quantitative trait locus (eQTL) decreasing ERAP1 expression; and ERAP2-rs10044354 (OR = 1.95; 95% CI, 1.55-2.44; P = 6.2e-09), an eQTL increasing ERAP2 expression. Furthermore, ERAP2-rs2248374 that disrupts ERAP2 expression is protective (OR = 0.56; 95% CI, 0.45-0.70; P = 2.39e-07). BSCR risk is additively increased when combining ERAP1/ERAP2 risk genotypes with two copies of HLA-Aw19 alleles (OR = 13.53; 95% CI, 3.79-54.77; P = 1.17e-05)., Conclusions: The genetic factors increasing BSCR risk demonstrate a pattern of increased processing, as well as increased presentation of ERAP2-specific peptides. This suggests a mechanism in which exceeding a peptide presentation threshold activates the immune response in choroids of A29 carriers.

Gene-set analyses are used to assess whether there is any evidence of association with disease among a set of biologically related genes. Such an analysis typically treats all genes within the sets similarly, even though there is substantial, external, information concerning the likely importance of each gene within each set. For example, for traits that are under purifying selection, we would expect genes showing extensive genic constraint to be more likely to be trait associated than unconstrained genes. Here we improve gene-set analyses by incorporating such external information into a higher-criticism-based signal detection analysis. We show that when this external information is predictive of whether a gene is associated with disease, our approach can lead to a significant increase in power. Further, our approach is particularly powerful when the signal is sparse, that is when only a small number of genes within the set are associated with the trait. We illustrate our approach with a gene-set analysis of amyotrophic lateral sclerosis (ALS) and implicate a number of gene-sets containing SOD1 and NEK1 as well as showing enrichment of small p values for gene-sets containing known ALS genes. We implement our approach in the R package wHC., (© 2020 Wiley Periodicals, Inc.)

Background: Studies have identified many common genetic associations that influence renal function and all-cause CKD, but these explain only a small fraction of variance in these traits. The contribution of rare variants has not been systematically examined., Methods: We performed exome sequencing of 3150 individuals, who collectively encompassed diverse CKD subtypes, and 9563 controls. To detect causal genes and evaluate the contribution of rare variants we used collapsing analysis, in which we compared the proportion of cases and controls carrying rare variants per gene., Results: The analyses captured five established monogenic causes of CKD: variants in PKD1 , PKD2 , and COL4A5 achieved study-wide significance, and we observed suggestive case enrichment for COL4A4 and COL4A3 . Beyond known disease-associated genes, collapsing analyses incorporating regional variant intolerance identified suggestive dominant signals in CPT2 and several other candidate genes. Biallelic mutations in CPT2 cause carnitine palmitoyltransferase II deficiency, sometimes associated with rhabdomyolysis and acute renal injury. Genetic modifier analysis among cases with APOL1 risk genotypes identified a suggestive signal in AHDC1 , implicated in Xia-Gibbs syndrome, which involves intellectual disability and other features. On the basis of the observed distribution of rare variants, we estimate that a two- to three-fold larger cohort would provide 80% power to implicate new genes for all-cause CKD., Conclusions: This study demonstrates that rare-variant collapsing analyses can validate known genes and identify candidate genes and modifiers for kidney disease. In so doing, these findings provide a motivation for larger-scale investigation of rare-variant risk contributions across major clinical CKD categories., (Copyright © 2019 by the American Society of Nephrology.)

Large-scale sequencing efforts in amyotrophic lateral sclerosis (ALS) have implicated novel genes using gene-based collapsing methods. However, pathogenic mutations may be concentrated in specific genic regions. To address this, we developed two collapsing strategies: One focuses rare variation collapsing on homology-based protein domains as the unit for collapsing, and the other is a gene-level approach that, unlike standard methods, leverages existing evidence of purifying selection against missense variation on said domains. The application of these two collapsing methods to 3093 ALS cases and 8186 controls of European ancestry, and also 3239 cases and 11,808 controls of diversified populations, pinpoints risk regions of ALS genes, including SOD1 , NEK1 , TARDBP , and FUS While not clearly implicating novel ALS genes, the new analyses not only pinpoint risk regions in known genes but also highlight candidate genes as well., (© 2019 Gelfman et al.; Published by Cold Spring Harbor Laboratory Press.)

Here we present an open-source R package 'meaRtools' that provides a platform for analyzing neuronal networks recorded on Microelectrode Arrays (MEAs). Cultured neuronal networks monitored with MEAs are now being widely used to characterize in vitro models of neurological disorders and to evaluate pharmaceutical compounds. meaRtools provides core algorithms for MEA spike train analysis, feature extraction, statistical analysis and plotting of multiple MEA recordings with multiple genotypes and treatments. meaRtools functionality covers novel solutions for spike train analysis, including algorithms to assess electrode cross-correlation using the spike train tiling coefficient (STTC), mutual information, synchronized bursts and entropy within cultured wells. Also integrated is a solution to account for bursts variability originating from mixed-cell neuronal cultures. The package provides a statistical platform built specifically for MEA data that can combine multiple MEA recordings and compare extracted features between different genetic models or treatments. We demonstrate the utilization of meaRtools to successfully identify epilepsy-like phenotypes in neuronal networks from Celf4 knockout mice. The package is freely available under the GPL license (GPL> = 3) and is updated frequently on the CRAN web-server repository. The package, along with full documentation can be downloaded from: https://cran.r-project.org/web/packages/meaRtools/., Competing Interests: I have read the journal's policy and the authors of this manuscript have the following competing interests: DBG reports Equity from Pairnomix and Praxis, personal fees from Eli Lilly (travel honorarium) and AstraZeneca and research grants from Janssen, Gilead, Biogen and UCB. QW declares affiliation with and employment by Simcere Diagnostics Co, Ltd. MAB declares affiliation with and employment by Axion BioSystems. All other authors have declared that no competing interests exist.

Identifying the underlying causes of disease requires accurate interpretation of genetic variants. Current methods ineffectively capture pathogenic non-coding variants in genic regions, resulting in overlooking synonymous and intronic variants when searching for disease risk. Here we present the Transcript-inferred Pathogenicity (TraP) score, which uses sequence context alterations to reliably identify non-coding variation that causes disease. High TraP scores single out extremely rare variants with lower minor allele frequencies than missense variants. TraP accurately distinguishes known pathogenic and benign variants in synonymous (AUC = 0.88) and intronic (AUC = 0.83) public datasets, dismissing benign variants with exceptionally high specificity. TraP analysis of 843 exomes from epilepsy family trios identifies synonymous variants in known epilepsy genes, thus pinpointing risk factors of disease from non-coding sequence data. TraP outperforms leading methods in identifying non-coding variants that are pathogenic and is therefore a valuable tool for use in gene discovery and the interpretation of personal genomes.While non-coding synonymous and intronic variants are often not under strong selective constraint, they can be pathogenic through affecting splicing or transcription. Here, the authors develop a score that uses sequence context alterations to predict pathogenicity of synonymous and non-coding genetic variants, and provide a web server of pre-computed scores.

Cultured neuronal networks monitored with microelectrode arrays (MEAs) have been used widely to evaluate pharmaceutical compounds for potential neurotoxic effects. A newer application of MEAs has been in the development of in vitro models of neurological disease. Here, we directly evaluated the utility of MEAs to recapitulate in vivo phenotypes of mature microRNA-128 (miR-128) deficiency, which causes fatal seizures in mice. We show that inhibition of miR-128 results in significantly increased neuronal activity in cultured neuronal networks derived from primary mouse cortical neurons. These results support the utility of MEAs in developing in vitro models of neuroexcitability disorders, such as epilepsy, and further suggest that MEAs provide an effective tool for the rapid identification of microRNAs that promote seizures when dysregulated., (© 2016 McSweeney et al.; Published by Cold Spring Harbor Laboratory Press.)

During development, enhancers play pivotal roles in regulating gene expression programs; however, their involvement in cancer progression has not been fully characterized. We performed an integrative analysis of DNA methylation, RNA-seq, and small RNA-seq profiles from thousands of patients, including 25 diverse primary malignances and seven body sites of metastatic melanoma. We found that enhancers are consistently the most differentially methylated regions (DMR) as cancer progresses from normal to primary tumors and then to metastases, compared to other genomic features. Remarkably, identification of enhancer DMRs (eDMRs) enabled classification of primary tumors according to physiological organ systems, and in metastasis eDMRs are the most correlated with patient outcome. To further understand the eDMR role in cancer progression, we developed a model to predict genes and microRNAs that are regulated by enhancer and not promotor methylation, which shows high accuracy with chromatin architecture methods and was experimentally validated. Interestingly, among all metastatic melanoma eDMRs, the most correlated with patient survival were eDMRs that "switched" their methylation patterns back and forth between normal, primary, and metastases and target cancer drivers, e.g., KIT We further demonstrated that eDMR target genes were modulated in melanoma by the bone metastasis microenvironment, suggesting that eDMRs respond to microenvironmental cues in metastatic niches. Our findings that aberrant methylation in cancer cells mostly affects enhancers, which contribute to tumor progression and cancer cell plasticity, will facilitate development of epigenetic anticancer approaches., (© 2016 Bell et al.; Published by Cold Spring Harbor Laboratory Press.)

The global impact of DNA methylation on alternative splicing is largely unknown. Using a genome-wide approach in wild-type and methylation-deficient embryonic stem cells, we found that DNA methylation can either enhance or silence exon recognition and affects the splicing of more than 20% of alternative exons. These exons are characterized by distinct genetic and epigenetic signatures. Alternative splicing regulation of a subset of these exons can be explained by heterochromatin protein 1 (HP1), which silences or enhances exon recognition in a position-dependent manner. We constructed an experimental system using site-specific targeting of a methylated/unmethylated gene and demonstrate a direct causal relationship between DNA methylation and alternative splicing. HP1 regulates this gene's alternative splicing in a methylation-dependent manner by recruiting splicing factors to its methylated form. Our results demonstrate DNA methylation's significant global influence on mRNA splicing and identify a specific mechanism of splicing regulation mediated by HP1., (Copyright © 2015 The Authors. Published by Elsevier Inc. All rights reserved.)

DNA methylation is known to regulate transcription and was recently found to be involved in exon recognition via cotranscriptional splicing. We recently observed that exon-intron architectures can be grouped into two classes: one with higher GC content in exons compared to the flanking introns, and the other with similar GC content in exons and introns. The first group has higher nucleosome occupancy on exons than introns, whereas the second group exhibits weak nucleosome marking of exons, suggesting another type of epigenetic marker distinguishes exons from introns when GC content is similar. We find different and specific patterns of DNA methylation in each of the GC architectures; yet in both groups, DNA methylation clearly marks the exons. Exons of the leveled GC architecture exhibit a significantly stronger DNA methylation signal in relation to their flanking introns compared to exons of the differential GC architecture. This is accentuated by a reduction of the DNA methylation level in the intronic sequences in proximity to the splice sites and shows that different epigenetic modifications mark the location of exons already at the DNA level. Also, lower levels of methylated CpGs on alternative exons can successfully distinguish alternative exons from constitutive ones. Three positions at the splice sites show high CpG abundance and accompany elevated nucleosome occupancy in a leveled GC architecture. Overall, these results suggest that DNA methylation affects exon recognition and is influenced by the GC architecture of the exon and flanking introns.

Despite evidence that at the interspecific scale, exonic splicing silencers (ESSs) are under negative selection in constitutive exons, little is known about the effects of slightly deleterious polymorphisms on these splicing regulators. Through the application of a modified version of the McDonald-Kreitman test, we compared the normalized proportions of human polymorphisms and human/rhesus substitutions affecting exonic splicing regulators (ESRs) on sequences of constitutive and alternative exons. Our results show a depletion of substitutions and an enrichment of SNPs associated with ESS gain in constitutive exons. Moreover, we show that this evolutionary pattern is also present in a set of ESRs previously involved in the transition from constitutive to skipped exons in the mammalian lineage. The similarity between these two sets of ESRs suggests that the transition from constitutive to skipped exons in mammals is more frequently associated with the inhibition than with the promotion of splicing signals. This is in accordance with the hypothesis of a constitutive origin of exon skipping and corroborates previous findings about the antagonistic role of certain exonic splicing enhancers.

During evolution segments of homeothermic genomes underwent a GC content increase. Our analyses reveal that two exon-intron architectures have evolved from an ancestral state of low GC content exons flanked by short introns with a lower GC content. One group underwent a GC content elevation that abolished the differential exon-intron GC content, with introns remaining short. The other group retained the overall low GC content as well as the differential exon-intron GC content, and is associated with longer introns. We show that differential exon-intron GC content regulates exon inclusion level in this group, in which disease-associated mutations often lead to exon skipping. This group's exons also display higher nucleosome occupancy compared to flanking introns and exons of the other group, thus "marking" them for spliceosomal recognition. Collectively, our results reveal that differential exon-intron GC content is a previously unidentified determinant of exon selection and argue that the two GC content architectures reflect the two mechanisms by which splicing signals are recognized: exon definition and intron definition., (Copyright © 2012 The Authors. Published by Elsevier Inc. All rights reserved.)

Exon-intron architecture is one of the major features directing the splicing machinery to the short exons that are located within long flanking introns. However, the evolutionary dynamics of exon-intron architecture and its impact on splicing is largely unknown. Using a comparative genomic approach, we analyzed 17 vertebrate genomes and reconstructed the ancestral motifs of both 3' and 5' splice sites, as also the ancestral length of exons and introns. Our analyses suggest that vertebrate introns increased in length from the shortest ancestral introns to the longest primate introns. An evolutionary analysis of splice sites revealed that weak splice sites act as a restrictive force keeping introns short. In contrast, strong splice sites allow recognition of exons flanked by long introns. Reconstruction of the ancestral state suggests these phenomena were not prevalent in the vertebrate ancestor, but appeared during vertebrate evolution. By calculating evolutionary rate shifts in exons, we identified cis-acting regulatory sequences that became fixed during the transition from early vertebrates to mammals. Experimental validations performed on a selection of these hexamers confirmed their regulatory function. We additionally revealed many features of exons that can discriminate alternative from constitutive exons. These features were integrated into a machine-learning approach to predict whether an exon is alternative. Our algorithm obtains very high predictive power (AUC of 0.91), and using these predictions we have identified and successfully validated novel alternatively spliced exons. Overall, we provide novel insights regarding the evolutionary constraints acting upon exons and their recognition by the splicing machinery.