Identification and prioritization of myeloid malignancy germline variants in a large cohort of adult patients with AML.

Inherited predisposition to myeloid malignancies is more common than previously appreciated. We analyzed the whole-exome sequencing data of paired leukemia and skin biopsy samples from 391 adult patients from the Beat AML 1.0 consortium. Using the 2015 American College of Medical Genetics and Genomics (ACMG) guidelines for variant interpretation, we curated 1547 unique variants from 228 genes. The pathogenic/likely pathogenic (P/LP) germline variants were identified in 53 acute myeloid leukemia (AML) patients (13.6%) in 34 genes, including 6.39% (25/391) of patients harboring P/LP variants in genes considered clinically actionable (tier 1). 41.5% of the 53 patients with P/LP variants were in genes associated with the DNA damage response. The most frequently mutated genes were CHEK2 (8 patients) and DDX41 (7 patients). Pathogenic germline variants were also found in new candidate genes (DNAH5, DNAH9, DNMT3A, and SUZ12). No strong correlation was found between the germline mutational rate and age of AML onset. Among 49 patients who have a reported history of at least one family member affected with hematological malignancies, 6 patients harbored known P/LP germline variants and the remaining patients had at least one variant of uncertain significance, suggesting a need for further functional validation studies. Using CHEK2 as an example, we show that three-dimensional protein modeling can be one of the effective methodologies to prioritize variants of unknown significance for functional studies. Further, we evaluated an in silico approach that applies ACMG curation in an automated manner using the tool for assessment and (TAPES) prioritization in exome studies, which can minimize manual curation time for variants. Overall, our findings suggest a need to comprehensively understand the predisposition potential of many germline variants in order to enable closer monitoring for disease management and treatment interventions for affected patients and families.

Functional genomic landscape of acute myeloid leukaemia.

The implementation of targeted therapies for acute myeloid leukaemia (AML) has been challenging because of the complex mutational patterns within and across patients as well as a dearth of pharmacologic agents for most mutational events. Here we report initial findings from the Beat AML programme on a cohort of 672 tumour specimens collected from 562 patients. We assessed these specimens using whole-exome sequencing, RNA sequencing and analyses of ex vivo drug sensitivity. Our data reveal mutational events that have not previously been detected in AML. We show that the response to drugs is associated with mutational status, including instances of drug sensitivity that are specific to combinatorial mutational events. Integration with RNA sequencing also revealed gene expression signatures, which predict a role for specific gene networks in the drug response. Collectively, we have generated a dataset-accessible through the Beat AML data viewer (Vizome)-that can be leveraged to address clinical, genomic, transcriptomic and functional analyses of the biology of AML.

Secondary fusion proteins as a mechanism of BCR::ABL1 kinase-independent resistance in chronic myeloid leukaemia.

Drug resistance in chronic myeloid leukaemia (CML) may occur via mutations in the causative BCR::ABL1 fusion or BCR::ABL1-independent mechanisms. We analysed 48 patients with BCR::ABL1-independent resistance for the presence of secondary fusion genes by RNA sequencing. We identified 10 of the most frequently detected secondary fusions in 21 patients. Validation studies, cell line models, gene expression analysis and drug screening revealed differences with respect to proliferation rate, differentiation and drug sensitivity. Notably, expression of RUNX1::MECOM led to resistance to ABL1 tyrosine kinase inhibitors in vitro. These results suggest secondary fusions contribute to BCR::ABL1-independent resistance and may be amenable to combined therapies.

Genomic landscape of neutrophilic leukemias of ambiguous diagnosis.

Chronic neutrophilic leukemia (CNL), atypical chronic myeloid leukemia (aCML), and myelodysplastic/myeloproliferative neoplasms, unclassifiable (MDS/MPN-U) are a group of rare and heterogeneous myeloid disorders. There is strong morphologic resemblance among these distinct diagnostic entities as well as a lack of specific molecular markers and limited understanding of disease pathogenesis, which has made diagnosis challenging in certain cases. The treatment has remained empirical, resulting in dismal outcomes. We, therefore, performed whole-exome and RNA sequencing of these rare hematologic malignancies and present the most complete survey of the genomic landscape of these diseases to date. We observed a diversity of combinatorial mutational patterns that generally do not cluster within any one diagnosis. Gene expression analysis reveals enrichment, but not cosegregation, of clinical and genetic disease features with transcriptional clusters. In conclusion, these groups of diseases represent a continuum of related diseases rather than discrete diagnostic entities.

Smaller total brain volume but not subcortical structure volume related to common genetic risk for ADHD.

BACKGROUND: Mechanistic endophenotypes can inform process models of psychopathology and aid interpretation of genetic risk factors. Smaller total brain and subcortical volumes are associated with attention-deficit hyperactivity disorder (ADHD) and provide clues to its development. This study evaluates whether common genetic risk for ADHD is associated with total brain volume (TBV) and hypothesized subcortical structures in children. METHODS: Children 7-15 years old were recruited for a case-control study (N = 312, N = 199 ADHD). Children were assessed with a multi-informant, best-estimate diagnostic procedure and motion-corrected MRI measured brain volumes. Polygenic scores were computed based on discovery data from the Psychiatric Genomics Consortium (N = 19 099 ADHD, N = 34 194 controls) and the ENIGMA + CHARGE consortium (N = 26 577). RESULTS: ADHD was associated with smaller TBV, and altered volumes of caudate, cerebellum, putamen, and thalamus after adjustment for TBV; however, effects were larger and statistically reliable only in boys. TBV was associated with an ADHD polygenic score [ß = -0.147 (-0.27 to -0.03)], and mediated a small proportion of the effect of polygenic risk on ADHD diagnosis (average ACME = 0.0087, p = 0.012). This finding was stronger in boys (average ACME = 0.019, p = 0.008). In addition, we confirm genetic variation associated with whole brain volume, via an intracranial volume polygenic score. CONCLUSION: Common genetic risk for ADHD is not expressed primarily as developmental alterations in subcortical brain volumes, but appears to alter brain development in other ways, as evidenced by TBV differences. This is among the first demonstrations of this effect using molecular genetic data. Potential sex differences in these effects warrant further examination.

Whole-genome sequencing identifies EN1 as a determinant of bone density and fracture.

The extent to which low-frequency (minor allele frequency (MAF) between 1-5%) and rare (MAF ≤ 1%) variants contribute to complex traits and disease in the general population is mainly unknown. Bone mineral density (BMD) is highly heritable, a major predictor of osteoporotic fractures, and has been previously associated with common genetic variants, as well as rare, population-specific, coding variants. Here we identify novel non-coding genetic variants with large effects on BMD (ntotal = 53,236) and fracture (ntotal = 508,253) in individuals of European ancestry from the general population. Associations for BMD were derived from whole-genome sequencing (n = 2,882 from UK10K (ref. 10); a population-based genome sequencing consortium), whole-exome sequencing (n = 3,549), deep imputation of genotyped samples using a combined UK10K/1000 Genomes reference panel (n = 26,534), and de novo replication genotyping (n = 20,271). We identified a low-frequency non-coding variant near a novel locus, EN1, with an effect size fourfold larger than the mean of previously reported common variants for lumbar spine BMD (rs11692564(T), MAF = 1.6%, replication effect size = +0.20 s.d., Pmeta = 2 × 10(-14)), which was also associated with a decreased risk of fracture (odds ratio = 0.85; P = 2 × 10(-11); ncases = 98,742 and ncontrols = 409,511). Using an En1(cre/flox) mouse model, we observed that conditional loss of En1 results in low bone mass, probably as a consequence of high bone turnover. We also identified a novel low-frequency non-coding variant with large effects on BMD near WNT16 (rs148771817(T), MAF = 1.2%, replication effect size = +0.41 s.d., Pmeta = 1 × 10(-11)). In general, there was an excess of association signals arising from deleterious coding and conserved non-coding variants. These findings provide evidence that low-frequency non-coding variants have large effects on BMD and fracture, thereby providing rationale for whole-genome sequencing and improved imputation reference panels to study the genetic architecture of complex traits and disease in the general population.

Evaluating chronic emotional dysregulation and irritability in relation to ADHD and depression genetic risk in children with ADHD.

BACKGROUND: A central nosological problem concerns the etiological relationship of emotional dysregulation with ADHD. Molecular genetic risk scores provide a novel method for informing this question. METHODS: Participants were 514 community-recruited children of Northern European descent age 7-11 defined as ADHD or non-ADHD by detailed research evaluation. Parents-rated ADHD on standardized ratings and child temperament on the Temperament in Middle Childhood Questionnaire (TMCQ) and reported on ADHD and comorbid disorders by semi-structured clinical interview. Categorical and dimensional variables were created for ADHD, emotional dysregulation (implicating disruption of regulation of both anger-irritability and of positive valence surgency-sensation seeking), and irritability alone (anger dysregulation). Genome-wide polygenic risk scores (PRS) were computed for ADHD and depression genetic liability. Structural equation models and computationally derived emotion profiles guided analysis. RESULTS: The ADHD PRS was associated in variable-centered analyses with irritability (ß = .179, 95% CI = 0.087-0.280; ΔR2  = .034, p < .0002), but also with surgency/sensation seeking (B = .146, 95%CI = 0.052-0.240, ΔR2 =.022, p = .002). In person-centered analysis, the ADHD PRS was elevated in the emotion dysregulation ADHD group versus other ADHD children (OR = 1.44, 95% CI = 1.03-2.20, Nagelkerke ΔR2  = .013, p = .033) but did not differentiate irritable from surgent ADHD profiles. All effects were independent of variation in ADHD severity across traits or groups. The depression PRS was related to oppositional defiant disorder but not to ADHD emotion dysregulation. CONCLUSIONS: Irritability-anger and surgency-sensation seeking, as forms of negative and positively valenced dysregulated affect in ADHD populations, both relate principally to ADHD genetic risk and not mood-related genetic risk.

Polymorphisms in oxidative stress pathway genes and prostate cancer risk.

PURPOSE: Age-related factors including oxidative stress play an important role in prostate carcinogenesis. We hypothesize that germline single-nucleotide polymorphisms (SNPs) in oxidative stress pathway are associated with prostate cancer (PCa) risk. In this study, we aim to examine which of these SNPs is associated with PCa. METHODS: Participants included in this analyses came from the "Genetic Susceptibility, Environment and Prostate Cancer Risk Study" conducted at the Veterans Affairs Portland Health Care System. After applying exclusion criteria, 231 PCa cases and 382 prostate biopsy-negative controls who had genotyping data on twenty-two single-nucleotide polymorphisms (SNPs) in six genes (MAPK14, NRF2, CAT, GPX1, GSTP1, SOD2, and XDH) associated with oxidative stress pathway were included in the analyses. The genotyping of SNPs was conducted by the Illumina BeadXpress VeraCode platform. We investigated these SNPs in relation to overall and aggressive PCa risk using logistic regression models controlling for relevant covariates. RESULTS: One SNP in the MAPK14 (rs851023) was significantly associated with incident PCa risk. Compared to men carrying two copies of allele A, the presence of one or two copies of the G allele was associated with decreased risk of PCa [OR (95% CI) 0.19 (0.06-0.51)]. There was no statistically significant association between other SNPs in the NRF2, CAT, GPX1, GSTP1, SOD2, and XDH genes and PCa risk. CONCLUSIONS: The MAPK14 gene SNP rs851023 was associated with PCa and aggressive PCa risk after multiple comparison adjustment. Further studies in other populations or functional studies are needed to validate the finding.

Gender-Specific Effects of Selection for Drinking in the Dark on the Network Roles of Coding and Noncoding RNAs.

BACKGROUND: Transcriptional differences between heterogeneous stock mice and high drinking-in-the-dark selected mouse lines have previously been described based on microarray technology coupled with network-based analysis. The network changes were reproducible in 2 independent selections and largely confined to 2 distinct network modules; in contrast, differential expression appeared more specific to each selected line. This study extends these results by utilizing RNA-Seq technology, allowing evaluation of the relationship between genetic risk and transcription of noncoding RNA (ncRNA); we additionally evaluate sex-specific transcriptional effects of selection. METHODS: Naïve mice (N = 24/group and sex) were utilized for gene expression analysis in the ventral striatum; the transcriptome was sequenced with the Illumina HiSeq platform. Differential gene expression and the weighted gene co-expression network analysis were implemented largely as described elsewhere, resulting in the identification of genes that change expression level or (co)variance structure. RESULTS: Across both sexes, we detect selection effects on the extracellular matrix and synaptic signaling, although the identity of individual genes varies. A majority of nc RNAs cluster in a single module of relatively low density in both the male and female network. The most strongly differentially expressed transcript in both sexes was Gm22513, a small nuclear RNA with unknown function. Associated with selection, we also found a number of network hubs that change edge strength and connectivity. At the individual gene level, there are many sex-specific effects; however, at the annotation level, results are more concordant. CONCLUSIONS: In addition to demonstrating sex-specific effects of selection on the transcriptome, the data point to the involvement of extracellular matrix genes as being associated with the binge drinking phenotype.

NFAT5 and SLC4A10 Loci Associate with Plasma Osmolality.

Disorders of water balance, an excess or deficit of total body water relative to body electrolyte content, are common and ascertained by plasma hypo- or hypernatremia, respectively. We performed a two-stage genome-wide association study meta-analysis on plasma sodium concentration in 45,889 individuals of European descent (stage 1 discovery) and 17,637 additional individuals of European descent (stage 2 replication), and a transethnic meta-analysis of replicated single-nucleotide polymorphisms in 79,506 individuals (63,526 individuals of European descent, 8765 individuals of Asian Indian descent, and 7215 individuals of African descent). In stage 1, we identified eight loci associated with plasma sodium concentration at P<5.0 × 10-6 Of these, rs9980 at NFAT5 replicated in stage 2 meta-analysis (P=3.1 × 10-5), with combined stages 1 and 2 genome-wide significance of P=5.6 × 10-10 Transethnic meta-analysis further supported the association at rs9980 (P=5.9 × 10-12). Additionally, rs16846053 at SLC4A10 showed nominally, but not genome-wide, significant association in combined stages 1 and 2 meta-analysis (P=6.7 × 10-8). NFAT5 encodes a ubiquitously expressed transcription factor that coordinates the intracellular response to hypertonic stress but was not previously implicated in the regulation of systemic water balance. SLC4A10 encodes a sodium bicarbonate transporter with a brain-restricted expression pattern, and variant rs16846053 affects a putative intronic NFAT5 DNA binding motif. The lead variants for NFAT5 and SLC4A10 are cis expression quantitative trait loci in tissues of the central nervous system and relevant to transcriptional regulation. Thus, genetic variation in NFAT5 and SLC4A10 expression and function in the central nervous system may affect the regulation of systemic water balance.

Functional and genomic context in pathway analysis of GWAS data.

Gene set analysis (GSA) is a promising tool for uncovering the polygenic effects associated with complex diseases. However, the available techniques reflect a wide variety of hypotheses about how genetic effects interact to contribute to disease susceptibility. The lack of consensus about the best way to perform GSA has led to confusion in the field and has made it difficult to compare results across methods. A clear understanding of the various choices made during GSA - such as how gene sets are defined, how single-nucleotide polymorphisms (SNPs) are assigned to genes, and how individual SNP-level effects are aggregated to produce gene- or pathway-level effects - will improve the interpretability and comparability of results across methods and studies. In this review we provide an overview of the various data sources used to construct gene sets and the statistical methods used to test for gene set association, as well as provide guidelines for ensuring the comparability of results.

HitWalker2: visual analytics for precision medicine and beyond.

UNLABELLED: The lack of visualization frameworks to guide interpretation and facilitate discovery is a potential bottleneck for precision medicine, systems genetics and other studies. To address this we have developed an interactive, reproducible, web-based prioritization approach that builds on our earlier work. HitWalker2 is highly flexible and can utilize many data types and prioritization methods based upon available data and desired questions, allowing it to be utilized in a diverse range of studies such as cancer, infectious disease and psychiatric disorders. AVAILABILITY AND IMPLEMENTATION: Source code is freely available at https://github.com/biodev/HitWalker2 and implemented using Python/Django, Neo4j and Javascript (D3.js and jQuery). We support major open source browsers (e.g. Firefox and Chromium/Chrome). CONTACT: wilmotb@ohsu.edu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online. Additional information/instructions are available at https://github.com/biodev/HitWalker2/wiki.

Sequencing of sporadic Attention-Deficit Hyperactivity Disorder (ADHD) identifies novel and potentially pathogenic de novo variants and excludes overlap with genes associated with autism spectrum disorder.

Attention-Deficit Hyperactivity Disorder (ADHD) has high heritability; however, studies of common variation account for <5% of ADHD variance. Using data from affected participants without a family history of ADHD, we sought to identify de novo variants that could account for sporadic ADHD. Considering a total of 128 families, two analyses were conducted in parallel: first, in 11 unaffected parent/affected proband trios (or quads with the addition of an unaffected sibling) we completed exome sequencing. Six de novo missense variants at highly conserved bases were identified and validated from four of the 11 families: the brain-expressed genes TBC1D9, DAGLA, QARS, CSMD2, TRPM2, and WDR83. Separately, in 117 unrelated probands with sporadic ADHD, we sequenced a panel of 26 genes implicated in intellectual disability (ID) and autism spectrum disorder (ASD) to evaluate whether variation in ASD/ID-associated genes were also present in participants with ADHD. Only one putative deleterious variant (Gln600STOP) in CHD1L was identified; this was found in a single proband. Notably, no other nonsense, splice, frameshift, or highly conserved missense variants in the 26 gene panel were identified and validated. These data suggest that de novo variant analysis in families with independently adjudicated sporadic ADHD diagnosis can identify novel genes implicated in ADHD pathogenesis. Moreover, that only one of the 128 cases (0.8%, 11 exome, and 117 MIP sequenced participants) had putative deleterious variants within our data in 26 genes related to ID and ASD suggests significant independence in the genetic pathogenesis of ADHD as compared to ASD and ID phenotypes. © 2017 Wiley Periodicals, Inc.

Oncogenic CSF3R mutations in chronic neutrophilic leukemia and atypical CML.

BACKGROUND: The molecular causes of many hematologic cancers remain unclear. Among these cancers are chronic neutrophilic leukemia (CNL) and atypical (BCR-ABL1-negative) chronic myeloid leukemia (CML), both of which are diagnosed on the basis of neoplastic expansion of granulocytic cells and exclusion of genetic drivers that are known to occur in other myeloproliferative neoplasms and myeloproliferative-myelodysplastic overlap neoplasms. METHODS: To identify potential genetic drivers in these disorders, we used an integrated approach of deep sequencing coupled with the screening of primary leukemia cells obtained from patients with CNL or atypical CML against panels of tyrosine kinase-specific small interfering RNAs or small-molecule kinase inhibitors. We validated candidate oncogenes using in vitro transformation assays, and drug sensitivities were validated with the use of assays of primary-cell colonies. RESULTS: We identified activating mutations in the gene encoding the receptor for colony-stimulating factor 3 (CSF3R) in 16 of 27 patients (59%) with CNL or atypical CML. These mutations segregate within two distinct regions of CSF3R and lead to preferential downstream kinase signaling through SRC family-TNK2 or JAK kinases and differential sensitivity to kinase inhibitors. A patient with CNL carrying a JAK-activating CSF3R mutation had marked clinical improvement after the administration of the JAK1/2 inhibitor ruxolitinib. CONCLUSIONS: Mutations in CSF3R are common in patients with CNL or atypical CML and represent a potentially useful criterion for diagnosing these neoplasms. (Funded by the Leukemia and Lymphoma Society and others.).

Methylomic analysis of salivary DNA in childhood ADHD identifies altered DNA methylation in VIPR2.

BACKGROUND: Peripheral epigenetic marks hold promise for understanding psychiatric illness and may represent fingerprints of gene-environment interactions. We conducted an initial examination of CpG methylation variation in children with or without attention-deficit/hyperactivity disorder (ADHD). METHODS: Children age 7-12 were recruited, screened, evaluated and assigned to ADHD or non-ADHD groups by defined research criteria. Two independent age-matched samples were examined, a discovery set (n = 92, all boys, half control, half ADHD) and a confirmation set (n = 20, half ADHD, all boys). 5-methylcytosine levels were quantified in salivary DNA using the Illumina 450 K HumanMethylation array. Genes for which multiple probes were nominally significant and had a beta difference of at least 2% were evaluated for biological relevance and prioritized for confirmation and sequence validation. Gene pathways were explored and described. RESULTS: Two genes met the criteria for confirmation testing, VIPR2 and MYT1L; both had multiple probes meeting cutoffs and strong biological relevance. Probes on VIPR2 passed FDR correction in the confirmation set and were confirmed through bisulfite sequencing. Enrichment analysis suggested involvement of gene sets or pathways related to inflammatory processes and modulation of monoamine and cholinergic neurotransmission. CONCLUSIONS: Although it is unknown to what extent CpG methylation seen in peripheral tissue reflect transcriptomic changes in the brain, these initial results indicate that peripheral DNA methylation markers in ADHD may be promising and suggest targeted hypotheses for future study in larger samples.

Pathway analysis in attention deficit hyperactivity disorder: An ensemble approach.

Despite a wealth of evidence for the role of genetics in attention deficit hyperactivity disorder (ADHD), specific and definitive genetic mechanisms have not been identified. Pathway analyses, a subset of gene-set analyses, extend the knowledge gained from genome-wide association studies (GWAS) by providing functional context for genetic associations. However, there are numerous methods for association testing of gene sets and no real consensus regarding the best approach. The present study applied six pathway analysis methods to identify pathways associated with ADHD in two GWAS datasets from the Psychiatric Genomics Consortium. Methods that utilize genotypes to model pathway-level effects identified more replicable pathway associations than methods using summary statistics. In addition, pathways implicated by more than one method were significantly more likely to replicate. A number of brain-relevant pathways, such as RhoA signaling, glycosaminoglycan biosynthesis, fibroblast growth factor receptor activity, and pathways containing potassium channel genes, were nominally significant by multiple methods in both datasets. These results support previous hypotheses about the role of regulation of neurotransmitter release, neurite outgrowth and axon guidance in contributing to the ADHD phenotype and suggest the value of cross-method convergence in evaluating pathway analysis results. © 2016 Wiley Periodicals, Inc.

plethy: management of whole body plethysmography data in R.

BACKGROUND: Characterization of respiratory phenotypes can enhance complex trait and genomic studies involving allergic/autoimmune and infectious diseases. Many aspects of respiration can be measured using devices known as plethysmographs that can measure thoracic movement. One such approach (the Buxco platform) performs unrestrained whole body plethysmography on mice which infers thoracic movements from pressure differences from the act of inhalation and exhalation. While proprietary software is available to perform basic statistical analysis as part of machine's bundled software, it is desirable to be able to incorporate these analyses into high-throughput pipelines and integrate them with other data types, as well as leverage the wealth of analytic and visualization approaches provided by the R statistical computing environment. RESULTS: This manuscript describes the plethy package which is an R/Bioconductor framework for pre-processing and analysis of plethysmography data with emphasis on larger scale longitudinal experiments. The plethy package was designed to facilitate quality control and exploratory data analysis. We provide a demonstration of the features of plethy using a dataset assessing the respiratory effects over time of SARS and Influenza infection in mice. CONCLUSION: The plethy package provides functionality for users to import, perform quality assessment and exploratory data analysis in a manner that allows interoperability with existing modelling tools. Our package is implemented in R and is freely available as part of the Bioconductor project http://www.bioconductor.org/packages/release/bioc/html/plethy.html .

Splicing landscape of the eight collaborative cross founder strains.

BACKGROUND: The Collaborative Cross (CC) is a large panel of genetically diverse recombinant inbred mouse strains specifically designed to provide a systems genetics resource for the study of complex traits. In part, the utility of the CC stems from the extensive genome-wide annotations of founder strain sequence and structural variation. Still missing, however, are transcriptome-specific annotations of the CC founder strains that could further enhance the utility of this resource. RESULTS: We provide a comprehensive survey of the splicing landscape of the 8 CC founder strains by leveraging the high level of alternative splicing within the brain. Using deep transcriptome sequencing, we found that a majority of the splicing landscape is conserved among the 8 strains, with ~65% of junctions being shared by at least 2 strains. We, however, found a large number of potential strain-specific splicing events as well, with an average of ~3000 and ~500 with ≥3 and ≥10 sequence read coverage, respectively, within each strain. To better understand strain-specific splicing within the CC founder strains, we defined criteria for and identified high-confidence strain-specific splicing events. These splicing events were defined as exon-exon junctions 1) found within only one strain, 2) with a read coverage ≥10, and 3) defined by a canonical splice site. With these criteria, a total of 1509 high-confidence strain-specific splicing events were identified, with the majority found within two of the wild-derived strains, CAST and PWK. Strikingly, the overwhelming majority, 94%, of these strain-specific splicing events are not yet annotated. Strain-specific splicing was also located within genomic regions recently reported to be over- and under-represented within CC populations. CONCLUSIONS: Phenotypic characterization of CC populations is increasing; thus these results will not only aid in further elucidating the transcriptomic architecture of the individual CC founder strains, but they will also help in guiding the utilization of the CC populations in the study of complex traits. This report is also the first to establish guidelines in defining and identifying strain-specific splicing across different mouse strains.

PLA2G6, encoding a phospholipase A2, is mutated in neurodegenerative disorders with high brain iron.

Neurodegenerative disorders with high brain iron include Parkinson disease, Alzheimer disease and several childhood genetic disorders categorized as neuroaxonal dystrophies. We mapped a locus for infantile neuroaxonal dystrophy (INAD) and neurodegeneration with brain iron accumulation (NBIA) to chromosome 22q12-q13 and identified mutations in PLA2G6, encoding a calcium-independent group VI phospholipase A2, in NBIA, INAD and the related Karak syndrome. This discovery implicates phospholipases in the pathogenesis of neurodegenerative disorders with iron dyshomeostasis.