Search

Although the route to generate microRNAs (miRNAs) is often depicted as a linear series of sequential and constitutive cleavages, we now appreciate multiple alternative pathways as well as diverse strategies to modulate their processing and function. Here, we identify an unusually profound regulatory role of conserved loop sequences in vertebrate pre-mir-144, which are essential for its cleavage by the Dicer RNase III enzyme in human and zebrafish models. Our data indicate that pre-mir-144 dicing is positively regulated via its terminal loop, and involves the ILF3 complex (NF90 and its partner NF45/ILF2). We provide further evidence that this regulatory switch involves reshaping of the pre-mir-144 apical loop into a structure that is appropriate for Dicer cleavage. In light of our recent findings that mir-144 promotes the nuclear biogenesis of its neighbor mir-451, these data extend the complex hierarchy of nuclear and cytoplasmic regulatory events that can control the maturation of clustered miRNAs.

The lack of genetically diverse preclinical animal models in basic biology and efficacy testing has been cited as a potential cause of failure in clinical trials. We developed and characterized five diverse RAG1 null mouse strains as models that allow xenografts to grow. In these strains, we characterized the growth of breast cancer, leukemia and glioma cell lines. We found a wide range of growth characteristics that were far more dependent on strain than tumor type. For the breast cancer cell line, we characterized the spectrum of xenograft/tumor growth at structural, histological, cellular and molecular levels across each strain, and found that each strain captures unique structural components of the stroma. Furthermore, we showed that the increase in tumor-infiltrating myeloid CD45+ cells and the amount of circulating cytokine IL-6 and chemokine KC (also known as CXCL1) is associated with a higher tumor size in different strains. This resource is available to study established human xenografts, as well as difficult-to-xenograft tumors and growth of hematopoietic stems cells, and to decipher the role of myeloid cells in the development of spontaneous cancers.

Patient-derived xenografts (PDX) are tumor-in-mouse models for cancer. PDX collections, such as the NCI PDXNet, are powerful resources for preclinical therapeutic testing. However, variations in experimental and analysis procedures have limited interpretability. To determine the robustness of PDX studies, the PDXNet tested temozolomide drug response for three prevalidated PDX models (sensitive, resistant, and intermediate) across four blinded PDX Development and Trial Centers using independently selected standard operating procedures. Each PDTC was able to correctly identify the sensitive, resistant, and intermediate models, and statistical evaluations were concordant across all groups. We also developed and benchmarked optimized PDX informatics pipelines, and these yielded robust assessments across xenograft biological replicates. These studies show that PDX drug responses and sequence results are reproducible across diverse experimental protocols. In addition, we share the range of experimental procedures that maintained robustness, as well as standardized cloud-based workflows for PDX exome-sequencing and RNA-sequencing analyses and for evaluating growth.Significance:The PDXNet Consortium shows that PDX drug responses and sequencing results are reproducible across diverse experimental protocols, establishing the potential for multisite preclinical studies to translate into clinical trials.

To cure a patient's cancer is to eradicate invasive cells from the ecosystem of the body. However, the ecologic complexity of this challenge is not well understood. Here we show how results from eradications of invasive mammalian species from islands—one of the few contexts in which invasive species have been regularly cleared—inform new research directions for treating cancer. We first summarize the epidemiologic characteristics of island invader eradications and cancer treatments by analyzing recent datasets from the Database of Invasive Island Species Eradications and The Cancer Genome Atlas, detailing the superior successes of island eradication projects. Next, we compare how genetic and environmental factors impact success in each system. These comparisons illuminate a number of promising cancer research and treatment directions, such as heterogeneity engineering as motivated by gene drives and adaptive therapy; multiscale analyses of how population heterogeneity potentiates treatment resistance; and application of ecological data mining techniques to high-throughput cancer data. We anticipate that interdisciplinary comparisons between tumor progression and invasive species would inspire development of novel paradigms to cure cancer.

Significance Mutation of CLP1 causes pontocerebellar hypoplasia type 10 (PCH10), a neurodegenerative disorder associated with intellectual and motor disability. We made two mouse models of PCH10: one homozygous for the mutation found in patients and one heterozygous for this mutation and a null allele. Mutant mice had motor impairments and neurodegeneration in the spinal cord and cerebellum. Mutants also had altered tRNA metabolism; however, it is not clear whether these alterations contribute to pathogenesis. In addition, mutation of Clp1 resulted in altered poly(A) site selection and gene expression, suggesting that the role of CLP1 in mRNA 3′ end processing could be a promising avenue for future research into the pathogenesis of PCH10., Homozygous mutation of the RNA kinase CLP1 (cleavage factor polyribonucleotide kinase subunit 1) causes pontocerebellar hypoplasia type 10 (PCH10), a pediatric neurodegenerative disease. CLP1 is associated with the transfer RNA (tRNA) splicing endonuclease complex and the cleavage and polyadenylation machinery, but its function remains unclear. We generated two mouse models of PCH10: one homozygous for the disease-associated Clp1 mutation, R140H, and one heterozygous for this mutation and a null allele. Both models exhibit loss of lower motor neurons and neurons of the deep cerebellar nuclei. To explore whether Clp1 mutation impacts tRNA splicing, we profiled the products of intron-containing tRNA genes. While mature tRNAs were expressed at normal levels in mutant mice, numerous other products of intron-containing tRNA genes were dysregulated, with pre-tRNAs, introns, and certain tRNA fragments up-regulated, and other fragments down-regulated. However, the spatiotemporal patterns of dysregulation do not correlate with pathogenicity for most altered tRNA products. To elucidate the effect of Clp1 mutation on precursor messenger RNA (pre-mRNA) cleavage, we analyzed poly(A) site (PAS) usage and gene expression in Clp1R140H/− spinal cord. PAS usage was shifted from proximal to distal sites in the mutant mouse, particularly in short and closely spaced genes. Many such genes were also expressed at lower levels in the Clp1R140H/− mouse, possibly as a result of impaired transcript maturation. These findings are consistent with the hypothesis that select genes are particularly dependent upon CLP1 for proper pre-mRNA cleavage, suggesting that impaired mRNA 3′ processing may contribute to pathogenesis in PCH10.

Ribosome-associated quality control pathways respond to defects in translational elongation to recycle arrested ribosomes and degrade aberrant polypeptides and mRNAs. Loss of a tRNA gene leads to ribosomal pausing that is resolved by the translational GTPase GTPBP2, and in its absence causes neuron death. Here, we show that loss of the homologous protein GTPBP1 during tRNA deficiency in the mouse brain also leads to codon-specific ribosome pausing and neurodegeneration, suggesting that these non-redundant GTPases function in the same pathway to mitigate ribosome pausing. As observed in

Colonization of visceral organs with melanoma in humanized NSG-SGM3 mice is dependent upon human CD33+CD11b+CD117+ progenitor cells. A gene signature of KIT/CD117–expressing CD33+ subset correlates with decreased overall survival in TCGA melanoma samples and represents a novel candidate biomarker., Metastasis of melanoma significantly worsens prognosis; thus, therapeutic interventions that prevent metastasis could improve patient outcomes. Here, we show using humanized mice that colonization of distant visceral organs with melanoma is dependent upon a human CD33+CD11b+CD117+ progenitor cell subset comprising, Graphical Abstract

Translation-dependent quality control pathways such as no-go decay (NGD), non-stop decay (NSD), and nonsense-mediated decay (NMD) govern protein synthesis and proteostasis by resolving non-translating ribosomes and preventing the production of potentially toxic peptides derived from faulty and aberrant mRNAs. However, how translation is altered and the in vivo defects that arise in the absence of these pathways are poorly understood. Here, we show that the NGD/NSD factorsPeloandHbs1lare critical in mice for cerebellar neurogenesis but expendable for survival of these neurons after development. Analysis of mutant mouse embryonic fibroblasts revealed translational pauses, alteration of signaling pathways, and translational reprogramming. Similar effects on signaling pathways, including mTOR activation, the translatome and mouse cerebellar development were observed upon deletion of the NMD factorUpf2. Our data reveal that these quality control pathways that function to mitigate errors at distinct steps in translation can evoke similar cellular responses.

Ribosome-associated quality control pathways respond to defects in translational elongation to recycle arrested ribosomes and degrade aberrant polypeptides and mRNAs. Loss of a tRNA gene leads to ribosomal pausing that is resolved by the translational GTPase GTPBP2, and in its absence causes neuron death. Here, we show that loss of the homologous protein GTPBP1 during tRNA deficiency in the mouse brain also leads to codon-specific ribosome pausing and neurodegeneration, suggesting that these non-redundant GTPases function in the same pathway to mitigate ribosome pausing. As observed in Gtpbp2-/- mice (Ishimura et al., 2016), GCN2-mediated activation of the integrated stress response (ISR) was apparent in the Gtpbp1-/- brain. We observed decreased mTORC1 signaling which increased neuronal death, whereas ISR activation was neuroprotective. Our data demonstrate that GTPBP1 functions as an important quality control mechanism during translation elongation and suggest that translational signaling pathways intricately interact to regulate neuronal homeostasis during defective elongation.

Development of candidate cancer treatments is a resource-intensive process, with the research community continuing to investigate options beyond static genomic characterization. Toward this goal, we have established the genomic landscapes of 536 patient-derived xenograft (PDX) models across 25 cancer types, together with mutation, copy number, fusion, transcriptomic profiles, and NCI-MATCH arms. Compared with human tumors, PDXs typically have higher purity and fit to investigate dynamic driver events and molecular properties via multiple time points from same case PDXs. Here, we report on dynamic genomic landscapes and pharmacogenomic associations, including associations between activating oncogenic events and drugs, correlations between whole-genome duplications and subclone events, and the potential PDX models for NCI-MATCH trials. Lastly, we provide a web portal having comprehensive pan-cancer PDX genomic profiles and source code to facilitate identification of more druggable events and further insights into PDXs' recapitulation of human tumors.

Background Patient-derived xenograft (PDX) models are in vivo models of human cancer that have been used for translational cancer research and therapy selection for individual patients. The Jackson Laboratory (JAX) PDX resource comprises 455 models originating from 34 different primary sites (as of 05/08/2019). The models undergo rigorous quality control and are genomically characterized to identify somatic mutations, copy number alterations, and transcriptional profiles. Bioinformatics workflows for analyzing genomic data obtained from human tumors engrafted in a mouse host (i.e., Patient-Derived Xenografts; PDXs) must address challenges such as discriminating between mouse and human sequence reads and accurately identifying somatic mutations and copy number alterations when paired non-tumor DNA from the patient is not available for comparison. Results We report here data analysis workflows and guidelines that address these challenges and achieve reliable identification of somatic mutations, copy number alterations, and transcriptomic profiles of tumors from PDX models that lack genomic data from paired non-tumor tissue for comparison. Our workflows incorporate commonly used software and public databases but are tailored to address the specific challenges of PDX genomics data analysis through parameter tuning and customized data filters and result in improved accuracy for the detection of somatic alterations in PDX models. We also report a gene expression-based classifier that can identify EBV-transformed tumors. We validated our analytical approaches using data simulations and demonstrated the overall concordance of the genomic properties of xenograft tumors with data from primary human tumors in The Cancer Genome Atlas (TCGA). Conclusions The analysis workflows that we have developed to accurately predict somatic profiles of tumors from PDX models that lack normal tissue for comparison enable the identification of the key oncogenic genomic and expression signatures to support model selection and/or biomarker development in therapeutic studies. A reference implementation of our analysis recommendations is available at https://github.com/TheJacksonLaboratory/PDX-Analysis-Workflows. Electronic supplementary material The online version of this article (10.1186/s12920-019-0551-2) contains supplementary material, which is available to authorized users.

Despite advances in therapy, glioblastoma remains an incurable disease with a dismal prognosis. Recent studies have implicated cancer stem cells within glioblastoma (glioma stem cells, GSCs) as mediators of therapeutic resistance and tumor progression. In this study, we investigated the role of the transforming growth factor-β (TGF-β) superfamily, which has been found to play an integral role in the maintenance of stem cell homeostasis within multiple stem cell systems, as a mediator of stem-like cells in glioblastoma. We find that BMP and TGF-β signaling define divergent molecular and functional identities in glioblastoma, and mark relatively quiescent and proliferative GSCs, respectively. Treatment of GSCs with BMP inhibits cell proliferation, but does not abrogate their stem-ness, as measured by self-renewal and tumorigencity. Further, BMP pathway activation confers relative resistance to radiation and temozolomide chemotherapy. Our findings define a quiescent cancer stem cell population in glioblastoma that may be a cellular reservoir for tumor recurrence following cytotoxic therapy.

Purpose: Genomic and transcriptomic alterations during metastasis are considered to affect clinical outcome of colorectal cancers, but detailed clinical implications of metastatic alterations are not fully uncovered. We aimed to investigate the effect of metastatic evolution on in vivo treatment outcome, and identify genomic and transcriptomic alterations associated with drug responsiveness. Experimental Design: We developed and analyzed patient-derived xenograft (PDX) models from 35 patients with colorectal cancer including 5 patients with multiple organ metastases (MOMs). We performed whole-exome, DNA methylation, and RNA sequencing for patient and PDX tumors. With samples from patients with MOMs, we conducted phylogenetic and subclonal analysis and in vivo drug efficacy test on the corresponding PDX models. Results: Phylogenetic analysis using mutation, expression, and DNA methylation data in patients with MOMs showed that mutational alterations were closely connected with transcriptomic and epigenomic changes during the tumor evolution. Subclonal analysis revealed that initial primary tumors with larger number of subclones exhibited more dynamic changes in subclonal architecture according to metastasis, and loco-regional and distant metastases occurred in a parallel or independent fashion. The PDX models from MOMs demonstrated therapeutic heterogeneity for targeted treatment, due to subclonal acquisition of additional mutations or transcriptomic activation of bypass signaling pathway during tumor evolution. Conclusions: This study demonstrated in vivo therapeutic heterogeneity of colorectal cancers using PDX models, and suggests that acquired subclonal alterations in mutations or gene expression profiles during tumor metastatic processes can be associated with the development of drug resistance and therapeutic heterogeneity of colorectal cancers.

Patient-derived xenografts (PDXs) are resected human tumors engrafted into mice for preclinical studies and therapeutic testing. It has been proposed that the mouse host affects tumor evolution during PDX engraftment and propagation, affecting the accuracy of PDX modeling of human cancer. Here, we exhaustively analyze copy number alterations (CNAs) in 1,451 PDX and matched patient tumor (PT) samples from 509 PDX models. CNA inferences based on DNA sequencing and microarray data displayed substantially higher resolution and dynamic range than gene expression-based inferences, and they also showed strong CNA conservation from PTs through late-passage PDXs. CNA recurrence analysis of 130 colorectal and breast PT/PDX-early/PDX-late trios confirmed high-resolution CNA retention. We observed no significant enrichment of cancer-related genes in PDX-specific CNAs across models. Moreover, CNA differences between patient and PDX tumors were comparable to variations in multiregion samples within patients. Our study demonstrates the lack of systematic copy number evolution driven by the PDX mouse host.

The processes by which tumors evolve are essential to the efficacy of treatment, but quantitative understanding of intratumoral dynamics has been limited. Although intratumoral heterogeneity is common, quantification of evolution is difficult from clinical samples because treatment replicates cannot be performed and because matched serial samples are infrequently available. To circumvent these problems we derived and assayed large sets of human triple-negative breast cancer xenografts and cell cultures from two patients, including 86 xenografts from cyclophosphamide, doxorubicin, cisplatin, docetaxel, or vehicle treatment cohorts as well as 45 related cell cultures. We assayed these samples via exome-seq and/or high-resolution droplet digital PCR, allowing us to distinguish complex therapy-induced selection and drift processes among endogenous cancer subclones with cellularity uncertainty

Epstein-Barr virus (EBV)-positive diffuse large B-cell lymphomas (EBV(+)-DLBLs) tend to occur in immunocompromised patients, such as the elderly or those undergoing solid organ transplantation. The pathogenesis and genomic characteristics of EBV(+)-DLBLs are largely unknown because of the limited availability of human samples and lack of experimental animal models. We observed the development of 25 human EBV(+)-DLBLs during the engraftment of gastric adenocarcinomas into immunodeficient mice. An integrated genomic analysis of the human-derived EBV(+)-DLBLs revealed enrichment of mutations in Rho pathway genes, including RHPN2, and Rho pathway transcriptomic activation. Targeting the Rho pathway using a Rho-associated protein kinase (ROCK) inhibitor, fasudil, markedly decreased tumor growth in EBV(+)-DLBL patient-derived xenograft (PDX) models. Thus, alterations in the Rho pathway appear to contribute to EBV-induced lymphomagenesis in immunosuppressed environments.

RNA editing refers to post-transcriptional processes that alter the base sequence of RNA. Recently, hundreds of new RNA editing targets have been reported. However, the mechanisms that determine the specificity and degree of editing are not well understood. We examined quantitative variation of site-specific editing in a genetically diverse multiparent population, Diversity Outbred mice, and mapped polymorphic loci that alter editing ratios globally for C-to-U editing and at specific sites for A-to-I editing. An allelic series in the C-to-U editing enzyme Apobec1 influences the editing efficiency of Apob and 58 additional C-to-U editing targets. We identified 49 A-to-I editing sites with polymorphisms in the edited transcript that alter editing efficiency. In contrast to the shared genetic control of C-to-U editing, most of the variable A-to-I editing sites were determined by local nucleotide polymorphisms in proximity to the editing site in the RNA secondary structure. Our results indicate that RNA editing is a quantitative trait subject to genetic variation and that evolutionary constraints have given rise to distinct genetic architectures in the two canonical types of RNA editing.

Ribosome stalling during translation has recently been shown to cause neurodegeneration, yet the signaling pathways triggered by stalled elongation complexes are unknown. To investigate these pathways we analyzed the brain of C57BL/6J-Gtpbp2nmf205-/- mice in which neuronal elongation complexes are stalled at AGA codons due to deficiencies in a tRNAArgUCU tRNA and GTPBP2, a mammalian ribosome rescue factor. Increased levels of phosphorylation of eIF2α (Ser51) were detected prior to neurodegeneration in these mice and transcriptome analysis demonstrated activation of ATF4, a key transcription factor in the integrated stress response (ISR) pathway. Genetic experiments showed that this pathway was activated by the eIF2α kinase, GCN2, in an apparent deacylated tRNA-independent fashion. Further we found that the ISR attenuates neurodegeneration in C57BL/6J-Gtpbp2nmf205-/- mice, underscoring the importance of cellular and stress context on the outcome of activation of this pathway. These results demonstrate the critical interplay between translation elongation and initiation in regulating neuron survival during cellular stress. DOI: http://dx.doi.org/10.7554/eLife.14295.001, eLife digest Information stored in DNA is used to make proteins in a two-step process. First, the DNA is copied to make molecules of messenger ribonucleic acid (or messenger RNA for short). Next, machines called ribosomes use the messenger RNAs as templates to assemble chains of amino acids – the building blocks of proteins – in a process called translation. Another type of RNA molecule called transfer RNA carries each amino acid to the ribosomes. If a specific transfer RNA is not available for translation at the right time, the ribosome might stall as it moves along the messenger RNA. At this point, the ribosome needs to be restarted or it will fall off the mRNA without finishing the protein. In 2014, a group of researchers reported that certain types of brain cells are very sensitive to ribosome stalling, and tend to die if translation does not continue. A protein called GTPBP2 was shown to play an important role in restarting stalled ribosomes in these cells. Here, Ishimura, Nagy et al. – including some of the researchers from the earlier work – investigated the molecular pathways that ribosome stalling triggers in brain cells using mutant mice that lacked the GTPBP2 protein. The experiments show that ribosome stalling activates an enzyme known as GCN2, which was already known to sense other types of malfunctions in cellular processes. Ishimura, Nagy et al. also show that GCN2 triggers stress responses in the cells by activating a communication system called the ATF4 pathway. This pathway protects the cells from damage, and its absence results in more rapid cell deterioration and death. The next challenges are to understand the exact mechanism by which GCN2 senses stalled ribosomes, and to find out how ribosome stalling causes the death of brain cells. DOI: http://dx.doi.org/10.7554/eLife.14295.002

Conserved noncoding elements (CNEs) in vertebrate genomes often act as developmental enhancers, but a critical issue is how well orthologous CNE sequences retain the same activity in their respective species, a characteristic important for generalization of model organism studies. To quantify how well CNE enhancer activity has been preserved, we compared the anatomy-specific activities of 41 zebra fish CNEs in zebra fish embryos with the activities of orthologous human CNEs in mouse embryos. We found that 13/41 (∼30%) of the orthologous CNE pairs exhibit conserved positive activity in zebra fish and mouse. Conserved positive activity is only weakly associated with either sequence conservation or the absence of bases undergoing accelerated evolution. A stronger effect is that disparate activity is associated with transcription factor binding site divergence. To distinguish the contributions of cis- versus trans-regulatory changes, we analyzed 13 CNEs in a three-way experimental comparison: human CNE tested in zebra fish, human CNE tested in mouse, and orthologous zebra fish CNE tested in zebra fish. Both cis- and trans-changes affect a significant fraction of CNEs, although human and zebra fish sequences exhibit disparate activity in zebra fish (indicating cis regulatory changes) twice as often as human sequences show disparate activity when tested in mouse and zebra fish (indicating trans regulatory changes). In all four cases where the zebra fish and human CNE display a similar expression pattern in zebra fish, the human CNE also displays a similar expression pattern in mouse. This suggests that the endogenous enhancer activity of ∼30% of human CNEs can be determined from experiments in zebra fish alone, and to identify these CNEs, both the zebra fish and the human sequences should be tested.

Alterations in mitochondrial function have long been considered a hallmark of cancer. We compared the lipidome and electron transport chain activities of non-synaptic brain mitochondria in two inbred mouse strains, the C57BL/6J (B6) and the VM/Dk (VM). The VM strain is unique in expressing a high incidence of spontaneous brain tumors (1.5%) that are mostly gliomas. The incidence of gliomas is about 210-fold greater in VM mice than in B6 mice. Using shotgun lipidomics, we found that the mitochondrial content of ethanolamine glycerophospholipid, phosphatidylserine, and ceramide was higher, whereas the content of total choline glycerophospholipid was lower in the VM mice than in B6 mice. Total cardiolipin content was similar in the VM and the B6 mice, but the distribution of cardiolipin molecular species differed markedly between the strains. B6 non-synaptic mitochondria contained 95 molecular species of cardiolipin that were symmetrically distributed over 7 major groups based on mass charge. In contrast, VM non-synaptic mitochondria contained only 42 molecular species that were distributed asymmetrically. The activities of Complex I, I/III, and II/III enzymes were lower, whereas the activity of complex IV was higher in the mitochondria of VM mice than in B6 mice. The high glioma incidence and alterations in electron transport chain activities in VM mice compared to B6 mice could be related to the unusual composition of mitochondrial lipids in the VM mouse brain.

In higher eukaryotes, transfer RNAs (tRNAs) with the same anticodon are encoded by multiple nuclear genes, and little is known about how mutations in these genes affect translation and cellular homeostasis. Similarly, the surveillance systems that respond to such defects in higher eukaryotes are not clear. Here, we discover that loss of GTPBP2, a novel binding partner of the ribosome recycling protein Pelota, in mice with a mutation in a tRNA gene that is specifically expressed in the central nervous system causes ribosome stalling and widespread neurodegeneration. Our results not only define GTPBP2 as a ribosome rescue factor but also unmask the disease potential of mutations in nuclear-encoded tRNA genes.

Many conserved noncoding sequences function as transcriptional enhancers that regulate gene expression. Here, we report that protein-coding DNA also frequently contains enhancers functioning at the transcriptional level. We tested the enhancer activity of 31 protein-coding exons, which we chose based on strong sequence conservation between zebrafish and human, and occurrence in developmental genes, using a Tol2 transposable GFP reporter assay in zebrafish. For each exon we measured GFP expression in hundreds of embryos in 10 anatomies via a novel system that implements the voice-recognition capabilities of a cellular phone. We find that 24/31 (77%) exons drive GFP expression compared to a minimal promoter control, and 14/24 are anatomy-specific (expression in four anatomies or less). GFP expression driven by these coding enhancers frequently overlaps the anatomies where the host gene is expressed (60%), suggesting self-regulation. Highly conserved coding sequences and highly conserved noncoding sequences do not significantly differ in enhancer activity (coding: 24/31 vs. noncoding: 105/147) or tissue-specificity (coding: 14/24 vs. noncoding: 50/105). Furthermore, coding and noncoding enhancers display similar levels of the enhancer-related histone modification H3K4me1 (coding: 9/24 vs noncoding: 34/81). Meanwhile, coding enhancers are over three times as likely to contain an H3K4me1 mark as other exons of the host gene. Our work suggests that developmental transcriptional enhancers do not discriminate between coding and noncoding DNA and reveals widespread dual functions in protein-coding DNA.

Cardiolipin is a class of mitochondrial specific phospholipid, which is intricately involved in mitochondrial functionality. Differences in cardiolipin species exist in a variety of tissues and diseases. It has been demonstrated that the cardiolipin profile is a key modulator of the functions of many mitochondrial proteins. However, the chemical mechanism(s) leading to normal and/or pathological distribution of cardiolipin species remain elusive. Herein, we describe a novel approach for investigating the molecular mechanism of cardiolipin remodeling through a dynamic simulation. This approach applied data from shotgun lipidomic analyses of the heart, liver, brain, and lung mitochondrial lipidomes to model cardiolipin remodeling, including relative content, regiospecificity, and isomeric composition of cardiolipin species. Generated cardiolipin profiles were nearly identical to those determined by shotgun lipidomics. Importantly, the simulated isomeric compositions of cardiolipin species were further substantiated through product ion analysis. Finally, unique enzymatic activities involved in cardiolipin remodeling were assessed from the parameters used in the dynamic simulation of cardiolipin profiles. Collectively, we described, verified, and demonstrated a novel approach by integrating both lipidomic analysis and dynamic simulation to study cardiolipin biology. We believe this study provides a foundation to investigate cardiolipin metabolism and bioenergetic homeostasis in normal and disease states.

Plasmodium falciparum, the most important causative agent of human malaria, undergoes antigenic variation as a means of prolonging infection and ensuring transmission between hosts. Clonal variation is observed in the surface adhesins expressed on infected erythrocytes: primarily in the PfEMP1 adhesin encoded by the large var gene family. The sirtuin PfSIR2A was the first protein discovered to have a major influence on antigenic variation in P. falciparum. In the absence of PfSIR2A, normal silencing of the variantly-expressed var gene family is partially deregulated. To thoroughly investigate the role of PfSIR2A in controlling antigenic variation, multiple independent clones of wildtype and PfSIR2A-knockout (DeltaSir2a) parasites were generated. var gene expression was then measured qualitatively, quantitatively and longitudinally over extended periods in culture. DeltaSir2a parasites were found to activate about 10 specific var genes in every independent clone analyzed. The activated genes were biased towards the upsA, upsBA and upsEvar gene subclasses. The total var transcript level was two to three-fold higher in DeltaSir2a parasites than in wildtype parasites and at least one transcript - encoding the pregnancy malaria adhesin VAR2CSA - was successfully translated and expressed on the infected cell surface. In the absence of PfSIR2A, antigenic switching over time was also diminished, although not abolished. This work expands our understanding of clonal antigenic variation in this important human pathogen and demonstrates a central role for PfSIR2A in regulating both the variant expression of specific var gene subsets and the overall quantity of var gene expression.

Summary: There are thousands of strongly conserved non-coding elements (CNEs) in vertebrate genomes, and their functions remain largely unknown. However, without biologically relevant criteria for prioritizing them, selecting a particular CNE sequences to study can be haphazard. To address this problem, we present cneViewer—a database and webtool that systematizes information on conserved non-coding DNA elements in zebrafish. A key feature here is the ability to search for CNEs that may be relevant to tissue-specific gene regulation, based on known developmental expression patterns of nearby genes. cneViewer provides this and other organizing features that significantly facilitate experimental design and CNE analysis. Availability: http://cneviewer.zebrafishcne.org Contact: chuangj@bc.edu

The objective of this study was to characterize the lipidome and electron transport chain activities in purified non-synaptic (NS) and synaptic (Syn) mitochondria from C57BL/6J mouse cerebral cortex. Contamination from subcellular membranes, especially myelin, has hindered past attempts to accurately characterize the lipid composition of brain mitochondria. An improved Ficoll and sucrose discontinuous gradient method was employed that yielded highly enriched mitochondrial populations free of myelin contamination. The activities of Complexes I, II, III, and II/III were lower in Syn than in NS mitochondria, while Complexes I/III and IV activities were similar in both populations. Shotgun lipidomics showed that levels of cardiolipin (Ptd(2)Gro) were lower, whereas levels of ceramide and phosphatidylserine were higher in Syn than in NS mitochondria. Coenzyme Q(9) and Q(10) was also lower in Syn than in NS mitochondria. Gangliosides, phosphatidic acid, sulfatides, and cerebrosides were undetectable in brain mitochondria. The distribution of Ptd(2)Gro molecular species was similar in both populations and formed a unique pattern, consisting of seven major molecular species groups, when arranged according to mass to charge ratios. Remodeling involving choline and ethanolamine phosphoglycerides could explain Ptd(2)Gro heterogeneity. NS and Syn mitochondrial lipidomic heterogeneity could influence energy metabolism, which may contribute to metabolic compartmentation of the brain.

Background The patterns of mutation vary both within and across genomes. It has been shown for a few mammals that mutation rates vary within the genome, while for unknown reasons, the sensu stricto yeasts have uniform rates instead. The generality of these observations has been unknown. Here we examine silent site substitutions in a more expansive set (20 mammals, 27 fungi, 4 insects) to determine why some genomes demonstrate this mosaic distribution and why others are uniform. Results We applied several intragene and intergene correlation tests to measure regional substitution patterns. Assuming that silent sites are a reasonable approximation to neutrally mutating sequence, our results show that all multicellular eukaryotes exhibit mutational heterogeneity. In striking contrast, all fungi are mutationally uniform – with the exception of three Candida species: C. albicans, C. dubliniensis, and C. tropicalis. We speculate that aspects of replication timing may be responsible for distinguishing these species. Our analysis also reveals classes of genes whose silent sites behave anomalously with respect to the mutational background in many species, indicating prevalent selective pressures. Genes associated with nucleotide binding or gene regulation have consistently low silent substitution rates in every mammalian species, as well as multiple fungi. On the other hand, receptor genes repeatedly exhibit high silent substitution rates, suggesting they have been influenced by diversifying selection. Conclusion Our findings provide a framework for understanding the regional mutational properties of eukaryotes, revealing a sharp difference between fungi and multicellular species. They also elucidate common selective pressures acting on eukaryotic silent sites, with frequent evidence for both purifying and diversifying selection.

Given that a gene has a high (or low) synonymous substitution rate in one mammalian species, will it also have a high (or low) synonymous substitution rate in another mammalian species? Such similarities in the rate of synonymous substitution can reveal both selective pressures and neutral processes acting on mammalian gene sequences; however, the existence of such an effect has been a matter of disagreement. We resolve whether such synonymous substitution rate similarities exist using 7462 ortholog triplets aligned across rat, mouse, and human, a dataset two orders of magnitude larger than previous studies. We find that a gene’s synonymous substitution rate in the rat-mouse branch of the phylogeny is correlated with its rate in the branch connecting human and the rat-mouse ancestor. We confirm this for several different measures of synonymous substitution rate, including corrections for base composition and CpG dinucleotides, and we verify the results in the larger mouse-human-rat-dog phylogeny. This similarity of rates is most apparent for genes in which synonymous sites are well conserved across species, suggesting that a significant component of the effect is due to purifying selection. We observe rate correlations at a resolution as fine as a few hundred kilobases, and the genes with the most similar synonymous substitution rates are enriched for regulatory functions. Genes with above-average substitution rates also exhibit significant, though somewhat weaker, rate correlations, suggesting that some neutral processes may have persisted in the phylogeny as well.

The regulation and maintenance of the cellular lipidome through biosynthetic, remodeling, and catabolic mechanisms are critical for biological homeostasis during development, health and disease. These complex mechanisms control the architectures of lipid molecular species, which have diverse yet highly regulated fatty acid chains at both the sn1 and sn2 positions. Phosphatidylcholine (PC) and phosphatidylethanolamine (PE) serve as the predominant biophysical scaffolds in membranes, acting as reservoirs for potent lipid signals and regulating numerous enzymatic processes. Here we report the first rigorous computational dissection of the mechanisms influencing PC and PE molecular architectures from high-throughput shotgun lipidomic data. Using novel statistical approaches, we have analyzed multidimensional mass spectrometry-based shotgun lipidomic data from developmental mouse heart and mature mouse heart, lung, brain, and liver tissues. We show that in PC and PE, sn1 and sn2 positions are largely independent, though for low abundance species regulatory processes may interact with both the sn1 and sn2 chain simultaneously, leading to cooperative effects. Chains with similar biochemical properties appear to be remodeled similarly. We also see that sn2 positions are more regulated than sn1, and that PC exhibits stronger cooperative effects than PE. A key aspect of our work is a novel statistically rigorous approach to determine cooperativity based on a modified Fisher's exact test using Markov Chain Monte Carlo sampling. This computational approach provides a novel tool for developing mechanistic insight into lipidomic regulation.

Technical advances in lipidomic analysis have generated tremendous amounts of quantitative lipid molecular species data, whose value has not been fully explored. We describe a novel computational method to infer mechanisms of de novo lipid synthesis and remodeling from lipidomic data. We focus on the mitochondrial-specific lipid cardiolipin (CL), a polyglycerol phospholipid with four acyl chains. The lengths and degree of unsaturation of these acyl chains vary across CL molecules, and regulation of these differences is important for mitochondrial energy metabolism. We developed a novel mathematical approach to determine mechanisms controlling the steady-state distribution of acyl chain combinations in CL . We analyzed mitochondrial lipids from 18 types of steady-state samples, each with at least 3 replicates, from mouse brain, heart, lung, liver, tumor cells, and tumors grown in vitro. Using a mathematical model for the CL remodeling mechanisms and a maximum likelihood approach to infer parameters, we found that for most samples the four chain positions have an independent and identical distribution, indicating they are remodeled by the same processes. Furthermore, for most brain samples and liver, the distribution of acyl chains is well-fit by a simple linear combination of the pools of acyl chains in phosphatidylcholine (PC), phosphatidylethanolamine (PE), and phosphatidylglycerol (PG). This suggests that headgroup chemistry is the key determinant of acyl donation into CL, with chain length/saturation less important. This canonical remodeling behavior appears damaged in some tumor samples, which display a consistent excess of CL molecules having particular masses. For heart and lung, the "proportional incorporation" assumption is not adequate to explain the CL distribution, suggesting additional acyl CoA-dependent remodeling that is chain-type specific. Our findings indicate that CL remodeling processes can be described by a small set of quantitative relationships, and that bioinformatic approaches can help determine these processes from high-throughput lipidomic data.

COMIT is presented; an algorithm for detecting functional non-coding motifs in coding regions, separating nucleotide and amino acid effects., Coding nucleotide sequences contain myriad functions independent of their encoded protein sequences. We present the COMIT algorithm to detect functional noncoding motifs in coding regions using sequence conservation, explicitly separating nucleotide from amino acid effects. COMIT concurs with diverse experimental datasets, including splicing enhancers, silencers, replication motifs, and microRNA targets, and predicts many novel functional motifs. Intriguingly, COMIT scores are well-correlated to scores uncalibrated for amino acids, suggesting that nucleotide motifs often override peptide-level constraints.

Glioblastomas (GBMs) are deadly tumors of the central nervous system. Most GBM exhibit homozygous deletions of the CDKN2A and CDKN2B tumor suppressors at 9p21.3, although loss of CDKN2A/B alone is insufficient to drive gliomagenesis. MIR31HG, which encodes microRNA-31 (miR-31), is a novel non-coding tumor suppressor positioned adjacent to CDKN2A/B at 9p21.3. We have determined that miR-31 expression is compromised in >72% of all GBM, and for patients, this predicts significantly shortened survival times independent of CDKN2A/B status. We show that miR-31 inhibits NF-κB signaling by targeting TRADD, its upstream activator. Moreover, upon reintroduction, miR-31 significantly reduces tumor burden and lengthens survival times in animal models. As such, our work identifies loss of miR-31 as a novel non-coding tumor-driving event in GBM.

Background The rate at which neutral (non-functional) bases undergo substitution is highly dependent on their location within a genome. However, it is not clear how fast these location-dependent rates change, or to what extent the substitution rate patterns are conserved between lineages. To address this question, which is critical not only for understanding the substitution process but also for evaluating phylogenetic footprinting algorithms, we examine ancestral repeats: a predominantly neutral dataset with a significantly higher genomic density than other datasets commonly used to study substitution rate variation. Using this repeat data, we measure the extent to which orthologous ancestral repeat sequences exhibit similar substitution patterns in separate mammalian lineages, allowing us to ascertain how well local substitution rates have been preserved across species. Results We calculated substitution rates for each ancestral repeat in each of three independent mammalian lineages (primate – from human/macaque alignments, rodent – from mouse/rat alignments, and laurasiatheria – from dog/cow alignments). We then measured the correlation of local substitution rates among these lineages. Overall we found the correlations between lineages to be statistically significant, but too weak to have much predictive power (r 2 <5%). These correlations were found to be primarily driven by regional effects at the scale of several hundred kb or larger. A few repeat classes (e.g. 7SK, Charlie8, and MER121) also exhibited stronger conservation of rate patterns, likely due to the effect of repeat-specific purifying selection. These classes should be excluded when estimating local neutral substitution rates. Conclusion Although local neutral substitution rates have some correlations among mammalian species, these correlations have little predictive power on the scale of individual repeats. This indicates that local substitution rates have changed significantly among the lineages we have studied, and are likely to have changed even more for more diverged lineages. The correlations that do persist are too weak to be responsible for many of the highly conserved elements found by phylogenetic footprinting algorithms, leading us to conclude that such elements must be conserved due to selective forces.

Background Little is known, either experimentally or computationally, about the genomic sequence features that regulate malaria genes. A sequence conservation analysis of the malaria species P. falciparum, P. berghei, P. yoelii, and P. chabaudi could significantly advance knowledge of malaria gene regulation. Results We computationally identify intergenic sequences conserved beyond neutral expectations, using a conservation algorithm that accounts for the strong compositional biases in malaria genomes. We first quantify the composition-specific divergence at silent positions in coding sequence. Using this as a background, we examine gene 5' regions, identifying 610 blocks conserved far beyond neutral expectations across the three mouse malariae, and 81 blocks conserved as strongly across all four species (p < 10-6). Detailed analysis of these blocks indicates that only a minor fraction are likely to be previously unknown coding sequences. Analogous noncoding conserved blocks have been shown to regulate adjacent genes in other phylogenies, making the predicted blocks excellent candidates for novel regulatory functions. We also find three potential transcription factor binding motifs which exhibit strong conservation and overrepresentation among the rodent malariae. Conclusion A broader finding of our analysis is that less malaria intergenic sequence has been conserved by selection than in yeast or vertebrate genomes. This supports the hypothesis that transcriptional regulation is simpler in malaria than other eukaryotic species. We have built a public database containing all sequence alignments and functional predictions, and we expect this to be a valuable resource to the malaria research community.