Tissue-Specific Evolution of Protein Coding Genes in Human and Mouse.

Protein-coding genes evolve at different rates, and the influence of different parameters, from gene size to expression level, has been extensively studied. While in yeast gene expression level is the major causal factor of gene evolutionary rate, the situation is more complex in animals. Here we investigate these relations further, especially taking in account gene expression in different organs as well as indirect correlations between parameters. We used RNA-seq data from two large datasets, covering 22 mouse tissues and 27 human tissues. Over all tissues, evolutionary rate only correlates weakly with levels and breadth of expression. The strongest explanatory factors of purifying selection are GC content, expression in many developmental stages, and expression in brain tissues. While the main component of evolutionary rate is purifying selection, we also find tissue-specific patterns for sites under neutral evolution and for positive selection. We observe fast evolution of genes expressed in testis, but also in other tissues, notably liver, which are explained by weak purifying selection rather than by positive selection.

Assessing Recent Selection and Functionality at Long Noncoding RNA Loci in the Mouse Genome.

Long noncoding RNAs (lncRNAs) are one of the most intensively studied groups of noncoding elements. Debate continues over what proportion of lncRNAs are functional or merely represent transcriptional noise. Although characterization of individual lncRNAs has identified approximately 200 functional loci across the Eukarya, general surveys have found only modest or no evidence of long-term evolutionary conservation. Although this lack of conservation suggests that most lncRNAs are nonfunctional, the possibility remains that some represent recent evolutionary innovations. We examine recent selection pressures acting on lncRNAs in mouse populations. We compare patterns of within-species nucleotide variation at approximately 10,000 lncRNA loci in a cohort of the wild house mouse, Mus musculus castaneus, with between-species nucleotide divergence from the rat (Rattus norvegicus). Loci under selective constraint are expected to show reduced nucleotide diversity and divergence. We find limited evidence of sequence conservation compared with putatively neutrally evolving ancestral repeats (ARs). Comparisons of sequence diversity and divergence between ARs, protein-coding (PC) exons and lncRNAs, and the associated flanking regions, show weak, but significantly lower levels of sequence diversity and divergence at lncRNAs compared with ARs. lncRNAs conserved deep in the vertebrate phylogeny show lower within-species sequence diversity than lncRNAs in general. A set of 74 functionally characterized lncRNAs show levels of diversity and divergence comparable to PC exons, suggesting that these lncRNAs are under substantial selective constraints. Our results suggest that, in mouse populations, most lncRNA loci evolve at rates similar to ARs, whereas older lncRNAs tend to show signals of selection similar to PC genes.

Embryonic gene expression of Coregonus palaea (whitefish) under pathogen stress as analyzed by high-throughput RNA-sequencing.

Most fishes produce free-living embryos that are exposed to environmental stressors immediately following fertilization, including pathogenic microorganisms. Initial immune protection of embryos involves the chorion, as a protective barrier, and maternally-allocated antimicrobial compounds. At later developmental stages, host-genetic effects influence susceptibility and tolerance, suggesting a direct interaction between embryo genes and pathogens. So far, only a few host genes could be identified that correlate with embryonic survival under pathogen stress in salmonids. Here, we utilized high-throughput RNA-sequencing in order to describe the transcriptional response of a non-model fish, the Alpine whitefish Coregonus palaea, to infection, both in terms of host genes that are likely manipulated by the pathogen, and those involved in an early putative immune response. Embryos were produced in vitro, raised individually, and exposed at the late-eyed stage to a virulent strain of the opportunistic fish pathogen Pseudomonas fluorescens. The pseudomonad increased embryonic mortality and affected gene expression substantially. For example, essential, upregulated metabolic pathways in embryos under pathogen stress included ion binding pathways, aminoacyl-tRNA-biosynthesis, and the production of arginine and proline, most probably mediated by the pathogen for its proliferation. Most prominently downregulated transcripts comprised the biosynthesis of unsaturated fatty acids, the citrate cycle, and various isoforms of b-cell transcription factors. These factors have been shown to play a significant role in host blood cell differentiation and renewal. With regard to specific immune functions, differentially expressed transcripts mapped to the complement cascade, MHC class I and II, TNF-alpha, and T-cell differentiation proteins. The results of this study reveal insights into how P. fluorescens impairs the development of whitefish embryos and set a foundation for future studies investigating host pathogen interactions in fish embryos.

Creatine deficiency syndomes : a new model of partial GAMT deficiency affecting brain cell development

Les syndromes de déficiences cérébrales en créatine (CCDS) sont dus à des mutations dans les gènes GATM et G AMT (codant pour les enzymes AGAT et G AMT de la voie de synthèse de créatine) ainsi que SLC6A8 (transporteur de créatine), et génèrent une absence ou une très forte baisse de créatine (Cr) dans le cerveau, mesurée par spectroscopic de résonance magnétique. Les patients CCDS développent des handicaps neurologiques sévères. Les patients AGAT et GAMT peuvent être traités avec des doses importantes de Cr, mais gardent dans la plupart des cas des séquelles neurologiques irréversibles. Aucun traitement efficace n'existe à ce jour pour la déficience en SLC6A8. Bien que de nombreux modèles aient été développés pour comprendre la Cr cérébrale en conditions physiologiques, les pathomécanismes des CCDS ne sont pas encore compris. Des souris transgéniques pour les gènes Gatm, Gamt et Slc6a8 ont été générées, mais elles ne miment que partiellement la pathologie humaine. Parmi les CCDS, la déficience en GAMT est la plus sévère, en raison de l'accumulation cérébrale de l'intermédiaire guanidinoacétate (GAA). Alors que la toxicité cérébrale du GAA a été étudiée par exposition directe au GAA d'animaux adultes sains, les mécanismes de la toxicité du GAA en condition de déficience en GAMT dans le cerveau en développement sont encore inconnus. Le but de ce projet était donc de développer un modèle de déficience en GAMT dans des cultures 3D primaires de cellules nerveuses de rat en agrégats par knock-down du gène GAMT, en utilisant un virus adéno-associé (AAV) induisant le mécanisme d'interférence à l'ARN (RNAi). Le virus scAAV2, à la multiplicité d'infection de 1000, s'est révélé le plus efficace pour transduire tous les types de cellules nerveuses des cultures (neurones, astrocytes, oligodendrocytes), et générer un knock-down maximal de la protéine GAMT de 85% (jour in vitro 18). Cette déficience partielle en GAMT s'est révélée insuffisante pour générer une déficience en Cr, mais a causé l'accumulation attendue de GAA, à des doses comparables aux niveaux observés dans le LCR des patients GAMT. Le GAA a induit une croissance axonale anarchique accompagnée d'une baisse de l'apoptose naturelle, suivis par une induction tardive de mort cellulaire non-apoptotique. Le co-traitement par la Cr a prévenu tous les effets toxiques du GAA. Ce travail montre que l'accumulation de GAA en absence de déficience en Cr est suffisante pour affecter le développement du tissu nerveux, et suggère que des formes de déficiences en GAMT supplémentaires, ne présentant pas de déficiences en Cr, pourraient être découvertes par mesure du GAA, en particulier à travers les programmes récemment proposés de dépistage néonatal de la déficience en GAMT. -- Cerebral creatine deficiency syndromes (CCDS) are caused by mutations in the genes GATM and GAMT (respectively coding for the two enzymes of the creatine synthetic pathway, AGAT and GAMT) as well as SLC6A8 (creatine transporter), and lead to the absence or very strong decrease of creatine (Cr) in the brain when measured by magnetic resonance spectroscopy. Affected patients show severe neurological impairments. While AGAT and GAMT deficient patients can be treated with high dosages of Cr, most remain with irreversible brain sequelae. No treatment has been successful so far for SLC6A8 deficiency. While many models have helped understanding the cerebral Cr pathways in physiological conditions, the pathomechanisms underlying CCDS are yet to be elucidated. Transgenic mice carrying mutations in the Gatm, Gamt and Slc6a8 genes have been developed, but only partially mimic the human pathology. Among CCDS, GAMT deficiency is the most severe, due to the CNS accumulation of the guanidinoacetate (GAA) intermediate. While brain toxicity of GAA has been explored through direct GAA exposure of adult healthy animals, the mechanisms underlying GAA toxicity in GAMT deficiency conditions on the developing CNS are yet unknown. The aim of this project was thus to develop and characterize a GAMT deficiency model in developing brain cells by gene knockdown, by adeno-associated virus (AAV)-driven RNA interference (RNAi) in rat 3D organotypic primary brain cell cultures in aggregates. scAAV2 with a multiplicity of infection of 1000 was shown as the most efficient serotype, was able to transduce all brain cell types (neurons, astrocytes, oligodendrocytes) and to induce a maximal GAMT protein knockdown of 85% (day in vitro 18). Metabolite analysis showed that partial GAMT knockdown was insufficient to induce Cr deficiency but generated the awaited GAA accumulation at concentrations comparable to the levels observed in cerebrospinal fluid of GAMT-deficient patients. Accumulated GAA induced axonal hypersprouting paralleled with inhibition of natural apoptosis, followed by a later induction in non-apoptotic cell death. Cr supplementation led to the prevention of all GAA-induced toxic effects. This work shows that GAA accumulation without Cr deficiency is sufficient to affect CNS development, and suggests that additional partial GAMT deficiencies, which may not show the classical brain Cr deficiency, may be discovered through GAA measurement including by recently proposed neonatal screening programs for GAMT deficiency.

Switching from tenofovir disoproxil fumarate to tenofovir alafenamide in antiretroviral regimens for virologically suppressed adults with HIV-1 infection: a randomised, active-controlled, multicentre, open-label, phase 3, non-inferiority study.

BACKGROUND: Antiretroviral regimens containing tenofovir disoproxil fumarate have been associated with renal toxicity and reduced bone mineral density. Tenofovir alafenamide is a novel tenofovir prodrug that reduces tenofovir plasma concentrations by 90%, thereby decreasing off-target side-effects. We aimed to assess whether efficacy, safety, and tolerability were non-inferior in patients switched to a regimen containing tenofovir alafenamide versus in those remaining on one containing tenofovir disoproxil fumarate. METHODS: In this randomised, actively controlled, multicentre, open-label, non-inferiority trial, we recruited HIV-1-infected adults from Gilead clinical studies at 168 sites in 19 countries. Patients were virologically suppressed (HIV-1 RNA <50 copies per mL) with an estimated glomerular filtration rate of 50 mL per min or greater, and were taking one of four tenofovir disoproxil fumarate-containing regimens for at least 96 weeks before enrolment. With use of a third-party computer-generated sequence, patients were randomly assigned (2:1) to receive a once-a-day single-tablet containing elvitegravir 150 mg, cobicistat 150 mg, emtricitabine 200 mg, and tenofovir alafenamide 10 mg (tenofovir alafenamide group) or to carry on taking one of four previous tenofovir disoproxil fumarate-containing regimens (tenofovir disoproxil fumarate group) for 96 weeks. Randomisation was stratified by previous treatment regimen in blocks of six. Patients and treating physicians were not masked to the assigned study regimen; outcome assessors were masked until database lock. The primary endpoint was the proportion of patients who received at least one dose of study drug who had undetectable viral load (HIV-1 RNA <50 copies per mL) at week 48. The non-inferiority margin was 12%. This study was registered with ClinicalTrials.gov, number NCT01815736. FINDINGS: Between April 12, 2013 and April 3, 2014, we enrolled 1443 patients. 959 patients were randomly assigned to the tenofovir alafenamide group and 477 to the tenofovir disoproxil fumarate group. Viral suppression at week 48 was noted in 932 (97%) patients assigned to the tenofovir alafenamide group and in 444 (93%) assigned to the tenofovir disoproxil fumarate group (adjusted difference 4·1%, 95% CI 1·6-6·7), with virological failure noted in ten and six patients, respectively. The number of adverse events was similar between the two groups, but study drug-related adverse events were more common in the tenofovir alafenamide group (204 patients [21%] vs 76 [16%]). Hip and spine bone mineral density and glomerular filtration were each significantly improved in patients in the tenofovir alafenamide group compared with those in the tenofovir disoproxil fumarate group. INTERPRETATION: Switching to a tenofovir alafenamide-containing regimen from one containing tenofovir disoproxil fumarate was non-inferior for maintenance of viral suppression and led to improved bone mineral density and renal function. Longer term follow-up is needed to better understand the clinical impact of these changes. FUNDING: Gilead Sciences.

Further molecular profiling of tumors harboring therapeutic targets within non-small cell lung cancer

Background: Treatment of NSCLC has been revolutionized in recent years with the introduction of several targeted therapies for selected genetically altered subtypes of NSCLC. A better understanding of molecular characteristics of NSCLC, which features common drug targets, may identify new therapeutic options. Methods: Over 6,700 non-small cell lung cancer cases referred to Caris Life Sciences between 2009 and 2014. Diagnoses and history were collected from referring physicians. Specific testing was performed per physician request and included a combination of sequencing (Sanger, NGS or pyrosequencing), protein expression (IHC), gene amplification/rearrangement (CISH or FISH), and/or RNA fragment analysis. Results: Tumors profiles from patients with hormone receptor positive disease (HER2, ER, PR, or AR positive by IHC) (n=629), HER2 mutations (n=8) ALK rearrangements (n=55), ROS1 rearrangement (n=17), cMET amplification or mutation (n=126), and cKIT mutation (n=11) were included in this analysis and compared to the whole cohort. Tumors with ALK rearrangement overexpressed AR in 18% of cases, and 7% presented with concomitant KRAS mutation. Lower rates of PTEN loss, as assessed by IHC, were observed in ALK positive (20%), ROS1 positive (9%) and cKIT mutated tumors (25%) compared to the overall NSCLC population (58%). cMET was overexpressed in 66% of ROS1 translocated and 57% of HER2 mutated tumors. cKIT mutations were found co-existing with APC (20%) and EGFR (20%) mutations. Pathway analysis revealed that hormone receptor positive disease carried more mutations in the ERK pathway (32%) compared to 9% in the mTOR pathway. 25% of patients with HER2 mutations harbored a co-existing mutation in the mTOR pathway. Conclusions: Pathway profiling reveals that NSCLC tumors present more often than reported with several concomitant alterations affecting the ERK or AKT pathway. Additionally, they are also characterized by the expression of potential biological modifiers of the cell cycle like hormonal receptors, representing a rationale for dual inhibition strategies in selected patients. Further refining of the understanding of NSCLC biomarker profile will optimize research for new treatment strategies.

Molecular profiling of non-small cell lung cancer by histologic subtype

Background: A substantial proportion of NSCLC has been shown to harbour specific molecular alterations affecting tumour proliferation and resulting in sensitivity to inhibition of the corresponding activated oncogenic pathway by targeted therapies. Comprehensive tumor profiling can diagnose such alterations and may identify new alterations opening additional treatment options for all distinct NSCLC subtypes. Methods: Over 6,700 non-small cell lung cancer cases referred to Caris Life Sciences between 2009 and 2014 were evaluated; clinical diagnoses and detailed tumor pathology were collected from referring physicians. Specific profiling was performed per physician request and included a combination of sequencing (Sanger, NGS or pyrosequencing), protein expression (IHC), gene amplification/rearrangement (CISH or FISH), and/or RNA fragment analysis within potential cancer-related genes and pathways. Results: Patients were grouped into cohorts according to histological subtype - adenocarcinoma (AD) (n=4,286), squamous cell carcinoma (SCC) (n=1,280), large cell carcinoma (LCC) (n=153) and bronchioalveolar carcinoma (BAC) (n=94). Protein overexpression of cMET (>2+ in >50% cells) was higher in AD (35.9%) compared to other subgroups (12-20%) while RRM1 and TOP2A levels were lower in AD. ALK or ROS1 were rearranged in 5.3% of patients with AD compared to 3.7% of patients with LCC and 1.2% of patients with SCC. EGFR mutations were found at low prevalence in both the LCC (0%) and SCC cohorts (2.8%) compared to 21% in AD. Similar lower rates of BRAF mutations were observed in the LCC and SCC cohorts compared to AD (0%, 1.1% and 5.1%). Pathway analysis showed activating mutations in the ERK pathway in 40% of patients with AD. Only 10-12% of patients with LCC or SCC had activating mutations in the ERK pathway. Conclusions: Despite the limitations of this retrospective series, we report comprehensive profiling of the largest cohort of NSCLC. Tumor profiling reveals that ADs may be more addicted to the ERK pathway than other histological subtypes. Drugs which target cMET may also have most utility in AD. Full analysis by histological subtype and additional correlative data on protein expression, gene copy number and mutations will be presented.