Very-long-chain acyl-CoA dehydrogenase (VLCAD) catalyzes the initial rate-limiting step in mitochondrial fatty acid beta-oxidation. VLCAD deficiency is clinically heterogenous, with three major phenotypes: a severe childhood form, with early onset, high mortality, and high incidence of cardiomyopathy; a milder childhood form, with later onset, usually with hypoketotic hypoglycemia as the main presenting feature, low mortality, and rare cardiomyopathy; and an adult form, with isolated skeletal muscle involvement, rhabdomyolysis, and myoglobinuria, usually triggered by exercise or fasting. To examine whether these different phenotypes are due to differences in the VLCAD genotype, we investigated 58 different mutations in 55 unrelated patients representing all known clinical phenotypes and correlated the mutation type with the clinical phenotype. Our results show a clear relationship between the nature of the mutation and the severity of disease. Patients with the severe childhood phenotype have mutations that result in no residual enzyme activity, whereas patients with the milder childhood and adult phenotypes have mutations that may result in residual enzyme activity. This clear genotype-phenotype relationship is in sharp contrast to what has been observed in medium-chain acyl-CoA dehydrogenase deficiency, in which no correlation between genotype and phenotype can be established.

Methotrexate is used to treat autoimmune diseases and malignancies, including acute lymphoblastic leukemia (ALL). Inter-individual variation in clearance of methotrexate results in heterogeneous systemic exposure, clinical efficacy, and toxicity. In a genome-wide association study of children with ALL, we identified SLCO1B1 as harboring multiple common polymorphisms associated with methotrexate clearance. The extent of influence of rare versus common variants on pharmacogenomic phenotypes remains largely unexplored. We tested the hypothesis that rare variants in SLCO1B1 could affect methotrexate clearance and compared the influence of common versus rare variants in addition to clinical covariates on clearance. From deep resequencing of SLCO1B1 exons in 699 children, we identified 93 SNPs, 15 of which were nonsynonymous (NS). Three of these NS SNPs were common, with a minor allele frequency (MAF) >5%, one had low frequency (MAF 1%-5%), and 11 were rare (MAF <1%). NS SNPs (common or rare) predicted to be functionally damaging were more likely to be found among patients with the lowest methotrexate clearance than patients with high clearance. We verified lower function in vitro of four SLCO1B1 haplotypes that were associated with reduced methotrexate clearance. In a multivariate stepwise regression analysis adjusting for other genetic and non-genetic covariates, SLCO1B1 variants accounted for 10.7% of the population variability in clearance. Of that variability, common NS variants accounted for the majority, but rare damaging NS variants constituted 17.8% of SLCO1B1's effects (1.9% of total variation) and had larger effect sizes than common NS variants. Our results show that rare variants are likely to have an important effect on pharmacogenetic phenotypes.

No abstract

The idea that point mutations in exons may affect splicing is intriguing and adds an additional layer of complexity when evaluating their possible effects. Even in the best-studied examples, the molecular mechanisms are not fully understood. Here, we use patient cells, model minigenes, and in vitro assays to show that a missense mutation in exon 5 of the medium-chain acyl-CoA dehydrogenase (MCAD) gene primarily causes exon skipping by inactivating a crucial exonic splicing enhancer (ESE), thus leading to loss of a functional protein and to MCAD deficiency. This ESE functions by antagonizing a juxtaposed exonic splicing silencer (ESS) and is necessary to define a suboptimal 3' splice site. Remarkably, a synonymous polymorphic variation in MCAD exon 5 inactivates the ESS, and, although this has no effect on splicing by itself, it makes splicing immune to deleterious mutations in the ESE. Furthermore, the region of MCAD exon 5 that harbors these elements is nearly identical to the exon 7 region of the survival of motor neuron (SMN) genes that contains the deleterious silent mutation in SMN2, indicating a very similar and finely tuned interplay between regulatory elements in these two genes. Our findings illustrate a mechanism for dramatic context-dependent effects of single-nucleotide polymorphisms on gene-expression regulation and show that it is essential that potential deleterious effects of mutations on splicing be evaluated in the context of the relevant haplotype.