A pharmacogenetic study implicates NINJ2 in the response to Interferon-β in multiple sclerosis

BACKGROUND

Multiple sclerosis (MS) is a disease in which biomarker identification is fundamental to predict response to treatments and to deliver the optimal drug to patients. We previously found an association between rs7298096, a polymorphism upstream to the NINJ2 gene, and the 4-year response to interferon-β (IFNβ) treatment in MS patients.

OBJECTIVES

To analyse the association between rs7298096 and time to first relapse (TTFR) during IFNβ therapy in MS patients and to better investigate its functional role.

METHODS

Survival analysis was applied in three MS cohorts from different countries (n = 1004). We also studied the role of the polymorphism on gene expression using GTEx portal and a luciferase assay. We interrogated GEO datasets to explore the relationship between NINJ2 expression, IFNβ and TTFR.

RESULTS

Rs7298096_AA patients show a shorter TTFR than rs7298096_G-carriers (P_{meta-analysis} = 3 × 10^-4, hazard ratio = 1.41). Moreover, rs7298096_AA is associated with a higher NINJ2 expression in blood (p = 7.0 × 10^-6), which was confirmed in vitro (p = 0.009). Finally, NINJ2 expression is downregulated by IFNβ treatment and related to TTFR.

CONCLUSIONS

Rs7298096 could influence MS disease activity during IFNβ treatment by modulating NINJ2 expression in blood. The gene encodes for an adhesion molecule involved in inflammation and endothelial cells activation, supporting its role in MS.

Biological Processes Modulating Longevity across Primates: A Phylogenetic Genome-Phenome Analysis

Aging is a complex process affecting different species and individuals in different ways. Comparing genetic variation across species with their aging phenotypes will help understanding the molecular basis of aging and longevity. Although most studies on aging have so far focused on short-lived model organisms, recent comparisons of genomic, transcriptomic, and metabolomic data across lineages with different lifespans are unveiling molecular signatures associated with longevity. Here, we examine the relationship between genomic variation and maximum lifespan across primate species. We used two different approaches. First, we searched for parallel amino-acid mutations that co-occur with increases in longevity across the primate linage. Twenty-five such amino-acid variants were identified, several of which have been previously reported by studies with different experimental setups and in different model organisms. The genes harboring these mutations are mainly enriched in functional categories such as wound healing, blood coagulation, and cardiovascular disorders. We demonstrate that these pathways are highly enriched for pleiotropic effects, as predicted by the antagonistic pleiotropy theory of aging. A second approach was focused on changes in rates of protein evolution across the primate phylogeny. Using the phylogenetic generalized least squares, we show that some genes exhibit strong correlations between their evolutionary rates and longevity-associated traits. These include genes in the Sphingosine 1-phosphate pathway, PI3K signaling, and the Thrombin/protease-activated receptor pathway, among other cardiovascular processes. Together, these results shed light into human senescence patterns and underscore the power of comparative genomics to identify pathways related to aging and longevity.

Effect of the pulsation type (alternate or simultaneous) on milk yield and health status of the mammary gland of Murciano-Granadina goats

Recent research into how the pulsation type affects goat milking concluded that alternate pulsation (AP) and simultaneous pulsation (SP) similarly affect the milking efficiency (milking duration and milking fractioning) and teat end condition after milking when both pulsations are employed in short-term experiments. However, the use of SP in high pipeline milking machines led to higher vacuum fluctuations, which, in the long term, may entail a potential risk factor for the mammary gland health status. The aim of this work was to study how AP and SP affect the milk yield, mammary gland health status, and milk composition of goats in the long term, during a complete lactation period, in a high pipeline milking machine. With this objective, 100 Murciano-Granadina goats with similar parturition dates (4 ± 1 postpartum weeks) were divided into 2 similar groups (50 goats per group) according their parity, milk yield, mammary gland health status, and milking duration (after a 15-d pre-experimental period). Throughout the experimental period (7 mo), one group was milked using SP and the other group using AP. Every month, samplings were carried out of the variables related to sanitary status of the mammary gland, milk fractioning, milking duration, average and maximum milk flows, milk composition, teat end status, vacuum fluctuations at the short milk tube during milking, and pulsation tests. Results showed that the use of SP in a high pipeline milking machine showed higher vacuum fluctuations without involving any difference in the other variables studied. We concluded that the use of AP in the milking of Murciano-Granadina goats offers little improvement of the milking performance and increases the investment required for installation and maintenance of the milking machines.

Genome-wide association analyses of risk tolerance and risky behaviors in over 1 million individuals identify hundreds of loci and shared genetic influences

Humans vary substantially in their willingness to take risks. In a combined sample of over 1 million individuals, we conducted genome-wide association studies (GWAS) of general risk tolerance, adventurousness, and risky behaviors in the driving, drinking, smoking, and sexual domains. Across all GWAS, we identified hundreds of associated loci, including 99 loci associated with general risk tolerance. We report evidence of substantial shared genetic influences across risk tolerance and the risky behaviors: 46 of the 99 general risk tolerance loci contain a lead SNP for at least one of our other GWAS, and general risk tolerance is genetically correlated ([Formula: see text] ~ 0.25 to 0.50) with a range of risky behaviors. Bioinformatics analyses imply that genes near SNPs associated with general risk tolerance are highly expressed in brain tissues and point to a role for glutamatergic and GABAergic neurotransmission. We found no evidence of enrichment for genes previously hypothesized to relate to risk tolerance.

Study of arterial stiffness in patients with an acute coronary event and chronic obstructive pulmonary disease confirmed by spirometry

BACKGROUND AND OBJECTIVES

Chronic obstructive pulmonary disease (COPD) worsens the prognosis for patients with an acute coronary event (ACE) treated with percutaneous coronary intervention. Objective To assess the effect of COPD on arterial stiffness in patients with an ACE.

METHODS

The study included patients with an ACE treated with percutaneous coronary intervention. At 1 month, postbronchodilation spirometry was performed, and arterial stiffness and markers of myocardial damage (troponin T and ProBNP) were measured.

RESULTS

We included 68 patients, 33% of whom had COPD (59% undiagnosed). The patients with COPD presented higher arterial stiffness values after adjusting for age and blood pressure readings. The troponin T and ProBNP levels were higher in the patients with COPD.

CONCLUSIONS

Arterial stiffness is greater in patients with an ACE if they have concomitant COPD. These findings can help explain the poorer prognosis of patients with both conditions.

128 Next-generation RNA sequencing of horse adipose and endometrial mesenchymal stem cells from the same donors unveils striking differences in their transcriptomic pattern

Earlier we successfully isolated and characterised endometrial (eMSC) and adipose (aMSC) mesenchymal stem cells from the same donors. Mesenchymal stem cells share biological traits but display different surface marker phenotype and migration ability. Here we extended our research to their mRNA signature using next-generation sequencing. The RNA from cells (3 biological replicates from each cell type and 3 technical replicates) at 90% confluence was extracted using a total RNA extraction kit and sent for mRNA-Seq (Norgen, Ontario, Canada; Illumina Sequencing Platform NextSEqn 500). Raw 76-bp single-end reads were aligned against the EquCab3 genome using RNA-STAR aligner. Counts were filtrated at a minimum of 5. Pairwise comparisons between the cell types were the input for gene ontology enrichment analysis. Only genes differentially expressed (DE) with 5 folds change (FC; P<0.05) were analysed. For DE analysis, eMSC were set as control and compared with aMSC. Unsupervised hierarchical clustering of the global gene expression signatures was done to compare the samples from each line using principal component analysis (PCA) and EdgeR: v3.20.9. Gene expression was normalized using FPKM. The heat map was built using R studio with G-plot package. A total of 14,896 transcripts with at least 5 reads were found; of these, 1598 were DE: 627 up-regulated (FC range: 2 to 236×) and 971 down-regulated (FC range: 2 to 464×) in eMSC. There was a marked dispersion in the FC of up- and down-regulated genes (>50×: 8 and 13; >20×, <50×: 9 and 17; >10×, <20×: 29 and 63; >5×, <10×: 91 and 130 and >2×, <5×: 490 and 748, respectively). Only genes DE with FC at least 5× were used for gene ontology and PCA analysis. Though 14,058 genes were common to both cell types, specific set of genes were found only in eMSC (n=162) or aMSC (n=676). Among the top 50 genes overexpressed in eMSC, several genes key for stem cell growth, immune response, migration and angiogenesis were found: TRIL, CXCL8, PDGF-D, SEMA5A, PTGS2, FGD, LAMA2, IL36G. In the top 50 down-regulated genes, some pivotal for osteoblast, adipogenic and neural differentiation were dramatically down-regulated (GPM6B, SCARA5 and NOTCH3 or NEFM, respectively), but no genes involved in immune rejection or stem cell proliferation were found. In gene ontology, the categories represented the most were cellular, developmental, metabolic, and immune system processes, as well as biological regulation, response to stimuli, organellar biogenesis, locomotion, localization and biological adhesion. Heat map and PCA analysis showed that one individual cell line from each type diverged markedly from the shared pattern. Individual variability of the donors may impinge upon the results; nevertheless, striking differences in the mRNA portfolio of eMSC and aMSC were detected. The importance and potential biological role of several of the genes and processes named above will be discussed in detail elsewhere. This work was supported by grant FONDECYT REGULAR 1150757 and the Government of Chile.

Phenotypic evaluation of ESBL- and carbapenemaseproducing Escherichia coli and Klebsiella pneumoniae from a teaching hospital in the Philippines

Antimicrobial resistance is a worldwide public health concern. Rise in the number of antimicrobial resistant organisms, such as extended spectrum β-lactamase- (ESBL) and carbapenemase-producing Escherichia coli and Klebsiella pneumoniae, continue to burden millions of people worldwide. E. coli and K. pneumoniae were isolated and collected for four months from a teaching hospital in the Philippines. All isolates were subjected to ESBL and carbapenemase testing using the double disk synergy test and modified Hodge test, respectively. Their pattern of resistance among different classes of antimicrobial agents was also investigated using the Kirby-Bauer disk diffusion test. Among the 32 clinical isolates tested, 28.1% were positive for ESBL production and 6.3% were positive for carbapenemase production. Species-specific classification showed that E. coli (44.4%) has the highest rate of ESBL production whereas both E. coli (5.6%) and K. pneumoniae (7.1%) showed almost similar rates of carbapenemase production. Antimicrobial resistance pattern of drug resistant isolates showed that all organisms were resistant to ampicillin, and majority showed resistance towards ciprofloxacin, cefotaxime, ceftriaxone, and sulfamethoxazole/trimethoprim. ESBL production is seen highest among E. coli isolates while similar rates of carbapenemase production was observed to both E. coli and K. pneumoniae isolates. Overall, antimicrobial resistance continues to rise and poses a huge threat in public health worldwide. Efforts should be made in developing rapid tests for antimicrobial resistance and to search for effective treatment from infections caused by multidrug resistant organisms.

Effect of Collapsed Duplications on Diversity Estimates: What to Expect

The study of segmental duplications (SDs) and copy-number variants (CNVs) is of great importance in the fields of genomics and evolution. However, SDs and CNVs are usually excluded from genome-wide scans for natural selection. Because of high identity between copies, SDs and CNVs that are not included in reference genomes are prone to be collapsed-that is, mistakenly aligned to the same region-when aligning sequence data from single individuals to the reference. Such collapsed duplications are additionally challenging because concerted evolution between duplications alters their site frequency spectrum and linkage disequilibrium patterns. To investigate the potential effect of collapsed duplications upon natural selection scans we obtained expectations for four summary statistics from simulations of duplications evolving under a range of interlocus gene conversion and crossover rates. We confirm that summary statistics traditionally used to detect the action of natural selection on DNA sequences cannot be applied to SDs and CNVs since in some cases values for known duplications mimic selective signatures. As a proof of concept of the pervasiveness of collapsed duplications, we analyzed data from the 1,000 Genomes Project. We find that, within regions identified as variable in copy number, diversity between individuals with the duplication is consistently higher than between individuals without the duplication. Furthermore, the frequency of single nucleotide variants (SNVs) deviating from Hardy-Weinberg Equilibrium is higher in individuals with the duplication, which strongly suggests that higher diversity is a consequence of collapsed duplications and incorrect evaluation of SNVs within these CNV regions.

NLRP3 polymorphisms and response to interferon-beta in multiple sclerosis patients

We aimed to investigate whether NLR family, pyrin domain containing 3 (NLRP3) polymorphisms are associated with the response to interferon-beta (IFNβ) in multiple sclerosis (MS) patients. A total of 14 NLRP3 polymorphisms were genotyped in a cohort of 665 relapsing-remitting MS patients recruited across 5 centers and classified into responders and non-responders according to clinical-radiological criteria after 1 year of IFNβ treatment. A meta-analysis failed to demonstrate significant associations between the response to IFNβ and NLRP3 polymorphisms. These findings do not support a role of polymorphisms located in the NLRP3 gene and the response to IFNβ in MS patients.

FaST-LMM for Two-Way Epistasis Tests on High-Performance Clusters

We introduce a version of the epistasis test in FaST-LMM for clusters of multithreaded processors. This new software maintains the sensitivity of the original FaST-LMM while delivering acceleration that is close to linear on 12-16 nodes of two recent platforms, with respect to improved implementation of FaST-LMM presented in an earlier work. This efficiency is attained through several enhancements on the original single-node version of FaST-LMM, together with the development of a message passing interface (MPI)-based version that ensures a balanced distribution of the workload as well as a multigraphics processing unit (GPU) module that can exploit the presence of multiple GPUs per node.

A 3-way hybrid approach to generate a new high-quality chimpanzee reference genome (Pan_tro_3.0)

The chimpanzee is arguably the most important species for the study of human origins. A key resource for these studies is a high-quality reference genome assembly; however, as with most mammalian genomes, the current iteration of the chimpanzee reference genome assembly is highly fragmented. In the current iteration of the chimpanzee reference genome assembly (Pan_tro_2.1.4), the sequence is scattered across more then 183 000 contigs, incorporating more than 159 000 gaps, with a genome-wide contig N50 of 51 Kbp. In this work, we produce an extensive and diverse array of sequencing datasets to rapidly assemble a new chimpanzee reference that surpasses previous iterations in bases represented and organized in large scaffolds. To this end, we show substantial improvements over the current release of the chimpanzee genome (Pan_tro_2.1.4) by several metrics, such as increased contiguity by >750% and 300% on contigs and scaffolds, respectively, and closure of 77% of gaps in the Pan_tro_2.1.4 assembly gaps spanning >850 Kbp of the novel coding sequence based on RNASeq data. We further report more than 2700 genes that had putatively erroneous frame-shift predictions to human in Pan_tro_2.1.4 and show a substantial increase in the annotation of repetitive elements. We apply a simple 3-way hybrid approach to considerably improve the reference genome assembly for the chimpanzee, providing a valuable resource for the study of human origins. Furthermore, we produce extensive sequencing datasets that are all derived from the same cell line, generating a broad non-human benchmark dataset.

Replicability and Prediction: Lessons and Challenges from GWAS

Since the publication of the Wellcome Trust Case Control Consortium (WTCCC) landmark study a decade ago, genome-wide association studies (GWAS) have led to the discovery of thousands of risk variants involved in disease etiology. This success story has two angles that are often overlooked. First, GWAS findings are highly replicable. This is an unprecedented phenomenon in complex trait genetics, and indeed in many areas of science, which in past decades have been plagued by false positives. At a time of increasing concerns about the lack of reproducibility, we examine the biological and methodological reasons that account for the replicability of GWAS and identify the challenges ahead. In contrast to the exemplary success of disease gene discovery, at present GWAS findings are not useful for predicting phenotypes. We close with an overview of the prospects for individualized prediction of disease risk and its foreseeable impact in clinical practice.

Dual MET and ERBB inhibition overcomes intratumor plasticity in osimertinib-resistant-advanced non-small-cell lung cancer (NSCLC)

Background

Third-generation epidermal growth factor receptor tyrosine kinase inhibitors (EGFR-TKIs) such as osimertinib are the last line of targeted treatment of metastatic non-small-cell lung cancer (NSCLC) EGFR-mutant harboring T790M. Different mechanisms of acquired resistance to third-generation EGFR-TKIs have been proposed. It is therefore crucial to identify new and effective strategies to overcome successive acquired mechanisms of resistance.

Methods

For Amplicon-seq analysis, samples from the index patient (primary and metastasis lesions at different timepoints) as well as the patient-derived orthotopic xenograft tumors corresponding to the different treatment arms were used. All samples were formalin-fixed paraffin-embedded, selected and evaluated by a pathologist. For droplet digital PCR, 20 patients diagnosed with NSCLC at baseline or progression to different lines of TKI therapies were selected. Formalin-fixed paraffin-embedded blocks corresponding to either primary tumor or metastasis specimens were used for analysis. For single-cell analysis, orthotopically grown metastases were dissected from the brain of an athymic nu/nu mouse and cryopreserved at -80°C.

Results

In a brain metastasis lesion from a NSCLC patient presenting an EGFR T790M mutation, we detected MET gene amplification after prolonged treatment with osimertinib. Importantly, the combination of capmatinib (c-MET inhibitor) and afatinib (ErbB-1/2/4 inhibitor) completely suppressed tumor growth in mice orthotopically injected with cells derived from this brain metastasis. In those mice treated with capmatinib or afatinib as monotherapy, we observed the emergence of KRAS G12C clones. Single-cell gene expression analyses also revealed intratumor heterogeneity, indicating the presence of a KRAS-driven subclone. We also detected low-frequent KRAS G12C alleles in patients treated with various EGFR-TKIs.

Conclusion

Acquired resistance to subsequent EGFR-TKI treatment lines in EGFR-mutant lung cancer patients may induce genetic plasticity. We assess the biological insights of tumor heterogeneity in an osimertinib-resistant tumor with acquired MET-amplification and propose new treatment strategies in this situation.

Clinical impact of the subclonal architecture and mutational complexity in chronic lymphocytic leukemia

Genome studies of chronic lymphocytic leukemia (CLL) have revealed the remarkable subclonal heterogeneity of the tumors, but the clinical implications of this phenomenon are not well known. We assessed the mutational status of 28 CLL driver genes by deep-targeted next-generation sequencing and copy number alterations (CNA) in 406 previously untreated patients and 48 sequential samples. We detected small subclonal mutations (0.6-25% of cells) in nearly all genes (26/28), and they were the sole alteration in 22% of the mutated cases. CNA tended to be acquired early in the evolution of the disease and remained stable, whereas the mutational heterogeneity increased in a subset of tumors. The prognostic impact of different genes was related to the size of the mutated clone. Combining mutations and CNA, we observed that the accumulation of driver alterations (mutational complexity) gradually shortened the time to first treatment independently of the clonal architecture, IGHV status and Binet stage. Conversely, the overall survival was associated with the increasing subclonal diversity of the tumors but it was related to the age of patients, IGHV and TP53 status of the tumors. In conclusion, our study reveals that both the mutational complexity and subclonal diversity influence the evolution of CLL.

Whole genome diversity of inherited chromosomally integrated HHV-6 derived from healthy individuals of diverse geographic origin

Human herpesviruses 6-A and -B (HHV-6A, HHV-6B) are ubiquitous in human populations worldwide. These viruses have been associated with several diseases such as multiple sclerosis, Hodgkin's lymphoma or encephalitis. Despite of the need to understand the genetic diversity and geographic stratification of these viruses, the availability of complete viral sequences from different populations is still limited. Here, we present nine new inherited chromosomally integrated HHV-6 sequences from diverse geographical origin which were generated through target DNA enrichment on lymphoblastoid cell lines derived from healthy individuals. Integration with available HHV-6 sequences allowed the assessment of HHV-6A and -6B phylogeny, patterns of recombination and signatures of natural selection. Analysis of the intra-species variability showed differences between A and B diversity levels and revealed that the HHV-6B reference (Z29) is an uncommon sequence, suggesting the need for an alternative reference sequence. Signs of geographical variation are present and more defined in HHV-6A, while they appear partly masked by recombination in HHV-6B. Finally, we conducted a scan for signatures of selection in protein coding genes that yielded at least 6 genes (4 and 2 respectively for the A and B species) showing significant evidence for accelerated evolution, and 1 gene showing evidence of positive selection in HHV-6A.

miR-328 mediates a metabolic shift in colon cancer cells by targeting SLC2A1/GLUT1

Purpose

Increasing evidence shows that altered metabolism is a critical hallmark in colon cancer. There is a strong need to explore the molecular mechanisms underlying cancer metabolism. Whether the aberrant expression of microRNAs contributes to cancer metabolism is not fully understood. miR-328 is a putative potential target of SLC2A1, but the regulating mechanism between them remains unknown. We have examined whether miR-328 directly regulates SLC2A1/GLUT1 expression in colon cancer cells.

Methods

We performed in silico bioinformatic analyses to identify miR-328-mediated molecular pathways and targets. We also performed luciferase assays and western blot analyses in LOVO and SW480 colon cancer cell lines. In addition, we assessed miR-328 expression in 47 paired tumor and normal tissue specimens from resected colon cancer patients.

Results

Luciferase reporter assays showed that miR-328 directly targeted SLC2A1 3′-untranslated region (UTR), with a significant decrease in luciferase activity in both LOVO and SW480 cell lines. These results were validated by western blot. miR-328 expression was significantly downregulated in tumor tissue compared with paired normal tissue.

Conclusions

Our results show that miR-328 targets SLC2A1/GLUT1. We suggest that miR-328 may be involved in the orchestration of the Warburg effect in colon cancer cells. Furthermore, miR‐328 expression is reduced in colon cancer patients and thus inversely correlates with the classically reported upregulated SLC2A1/GLUT1 expression in tumors.

Similar genomic proportions of copy number variation within gray wolves and modern dog breeds inferred from whole genome sequencing

BACKGROUND

Whole genome re-sequencing data from dogs and wolves are now commonly used to study how natural and artificial selection have shaped the patterns of genetic diversity. Single nucleotide polymorphisms, microsatellites and variants in mitochondrial DNA have been interrogated for links to specific phenotypes or signals of domestication. However, copy number variation (CNV), despite its increasingly recognized importance as a contributor to phenotypic diversity, has not been extensively explored in canids.

RESULTS

Here, we develop a new accurate probabilistic framework to create fine-scale genomic maps of segmental duplications (SDs), compare patterns of CNV across groups and investigate their role in the evolution of the domestic dog by using information from 34 canine genomes. Our analyses show that duplicated regions are enriched in genes and hence likely possess functional importance. We identify 86 loci with large CNV differences between dogs and wolves, enriched in genes responsible for sensory perception, immune response, metabolic processes, etc. In striking contrast to the observed loss of nucleotide diversity in domestic dogs following the population bottlenecks that occurred during domestication and breed creation, we find a similar proportion of CNV loci in dogs and wolves, suggesting that other dynamics are acting to particularly select for CNVs with potentially functional impacts.

CONCLUSIONS

This work is the first comparison of genome wide CNV patterns in domestic and wild canids using whole-genome sequencing data and our findings contribute to study the impact of novel kinds of genetic changes on the evolution of the domestic dog.

The impact of endogenous content, replicates and pooling on genome capture from faecal samples

Target-capture approach has improved over the past years, proving to be very efficient tool for selectively sequencing genetic regions of interest. These methods have also allowed the use of noninvasive samples such as faeces (characterized by their low quantity and quality of endogenous DNA) to be used in conservation genomic, evolution and population genetic studies. Here we aim to test different protocols and strategies for exome capture using the Roche SeqCap EZ Developer kit (57.5 Mb). First, we captured a complex pool of DNA libraries. Second, we assessed the influence of using more than one faecal sample, extract and/or library from the same individual, to evaluate its effect on the molecular complexity of the experiment. We validated our experiments with 18 chimpanzee faecal samples collected from two field sites as a part of the Pan African Programme: The Cultured Chimpanzee. Those two field sites are in Kibale National Park, Uganda (N = 9) and Loango National Park, Gabon (N = 9). We demonstrate that at least 16 libraries can be pooled, target enriched through hybridization, and sequenced allowing for the genotyping of 951,949 exome markers for population genetic analyses. Further, we observe that molecule richness, and thus, data acquisition, increase when using multiple libraries from the same extract or multiple extracts from the same sample. Finally, repeated captures significantly decrease the proportion of off-target reads from 34.15% after one capture round to 7.83% after two capture rounds, supporting our conclusion that two rounds of target enrichment are advisable when using complex faecal samples.

Potential damaging mutation in LRP5 from genome sequencing of the first reported chimpanzee with the Chiari malformation

The genus Pan is the closest related to humans (Homo sapiens) and it includes two species: Pan troglodytes (chimpanzees) and Pan paniscus (bonobos). Different characteristics, some of biomedical aspect, separate them from us. For instance, some common human medical conditions are rare in chimpanzees (menopause, Alzheimer disease) although it is unclear to which extent longevity plays an active role in these differences. However, both humans and chimpanzees present similar pathologies, thus, understanding traits in chimpanzees can help unravel the molecular basis of human conditions. Here, we sequenced the genome of Nico, a central chimpanzee diagnosed with a particular biomedical condition, the Chiari malformation. We performed a variant calling analysis comparing his genome to 25 whole genomes from healthy individuals (bonobos and chimpanzees), and after predicting the effects of the genetic variants, we looked for genes within the OMIM database. We found a novel, private, predicted as damaging mutation in Nico in LRP5, a gene related to bone density alteration pathologies, and we suggest a link between this mutation and his Chiari malformation as previously shown in humans. Our results reinforce the idea that a comparison between humans and chimpanzees can be established in this genetic frame of common diseases.

Integrating genomic alterations in diffuse large B-cell lymphoma identifies new relevant pathways and potential therapeutic targets

Genome studies of diffuse large B-cell lymphoma (DLBCL) have revealed a large number of somatic mutations and structural alterations. However, the clinical significance of these alterations is still not well defined. In this study, we have integrated the analysis of targeted next-generation sequencing of 106 genes and genomic copy number alterations (CNA) in 150 DLBCL. The clinically significant findings were validated in an independent cohort of 111 patients. Germinal center B-cell and activated B-cell DLBCL had a differential profile of mutations, altered pathogenic pathways and CNA. Mutations in genes of the NOTCH pathway and tumor suppressor genes (TP53/CDKN2A), but not individual genes, conferred an unfavorable prognosis, confirmed in the independent validation cohort. A gene expression profiling analysis showed that tumors with NOTCH pathway mutations had a significant modulation of downstream target genes, emphasizing the relevance of this pathway in DLBCL. An in silico drug discovery analysis recognized 69 (46%) cases carrying at least one genomic alteration considered a potential target of drug response according to early clinical trials or preclinical assays in DLBCL or other lymphomas. In conclusion, this study identifies relevant pathways and mutated genes in DLBCL and recognizes potential targets for new intervention strategies.

Genetic factors affecting EBV copy number in lymphoblastoid cell lines derived from the 1000 Genome Project samples

Epstein-Barr virus (EBV), human herpes virus 4, has been classically associated with infectious mononucleosis, multiple sclerosis and several types of cancers. Many of these diseases show marked geographical differences in prevalence, which points to underlying genetic and/or environmental factors. Those factors may include a different susceptibility to EBV infection and viral copy number among human populations. Since EBV is commonly used to transform B-cells into lymphoblastoid cell lines (LCLs) we hypothesize that differences in EBV copy number among individual LCLs may reflect differential susceptibility to EBV infection. To test this hypothesis, we retrieved whole-genome sequenced EBV-mapping reads from 1,753 LCL samples derived from 19 populations worldwide that were sequenced within the context of the 1000 Genomes Project. An in silico methodology was developed to estimate the number of EBV copy number in LCLs and validated these estimations by real-time PCR. After experimentally confirming that EBV relative copy number remains stable over cell passages, we performed a genome wide association analysis (GWAS) to try detecting genetic variants of the host that may be associated with EBV copy number. Our GWAS has yielded several genomic regions suggestively associated with the number of EBV genomes per cell in LCLs, unraveling promising candidate genes such as CAND1, a known inhibitor of EBV replication. While this GWAS does not unequivocally establish the degree to which genetic makeup of individuals determine viral levels within their derived LCLs, for which a larger sample size will be needed, it potentially highlighted human genes affecting EBV-related processes, which constitute interesting candidates to follow up in the context of EBV related pathologies.

Estimate of the total costs of allergic rhinitis in specialized care based on real-world data: the FERIN Study

BACKGROUND

Despite the socioeconomic importance of allergic rhinitis (AR), very few prospective studies have been performed under conditions of clinical practice and with a sufficiently long observation period outside the clinical trial scenario. We prospectively estimated the direct and indirect costs of AR in patients attending specialized clinics in Spain.

METHODS

Patients were recruited at random from allergy outpatient clinics in 101 health centers throughout Spain over 12 months. We performed a multicenter, observational, prospective study under conditions of clinical practice. We analyzed direct costs from a funder perspective (healthcare costs) and from a societal perspective (healthcare and non-healthcare costs). Indirect costs (absenteeism and presenteeism [productivity lost in the workplace]) were also calculated. The cost of treating conjunctivitis was evaluated alongside that of AR.

RESULTS

The total mean cost of AR per patient-year (n = 498) was €2326.70 (direct, €553.80; indirect, €1772.90). Direct costs were significantly higher in women (€600.34 vs €484.46, P = 0.02). Total costs for intermittent AR were significantly lower than for persistent AR (€1484.98 vs €2655.86, P < 0.001). Total indirect costs reached €1772.90 (presenteeism, €1682.71; absenteeism, €90.19). The direct costs of AR in patients with intermittent asthma (€507.35) were lower than in patients with mild-persistent asthma (€719.07) and moderate-persistent asthma (€798.71) (P = 0.006).

CONCLUSIONS

The total cost of AR for society is considerable. Greater frequency of symptoms and more severe AR are associated with higher costs. Indirect costs are almost threefold direct costs, especially in presenteeism. A reduction in presenteeism would generate considerable savings for society.