Gitools: Analysis and Visualisation of Genomic Data Using Interactive Heat-Maps

Intuitive visualization of data and results is very important in genomics, especially when many conditions are to be analyzed and compared. Heat-maps have proven very useful for the representation of biological data. Here we present Gitools (http://www.gitools.org), an open-source tool to perform analyses and visualize data and results as interactive heat-maps. Gitools contains data import systems from several sources (i.e. IntOGen, Biomart, KEGG, Gene Ontology), which facilitate the integration of novel data with previous knowledge

A framework to rank genomic alterations as targets for cancer precision medicine: the ESMO Scale for Clinical Actionability of molecular Targets (ESCAT)

Background: In order to facilitate implementation of precision medicine in clinical management of cancer, there is a need to harmonise and standardise the reporting and interpretation of clinically relevant genomics data. / Methods: The European Society for Medical Oncology (ESMO) Translational Research and Precision Medicine Working Group (TR and PM WG) launched a collaborative project to propose a classification system for molecular aberrations based on the evidence available supporting their value as clinical targets. A group of experts from several institutions was assembled to review available evidence, reach a consensus on grading criteria and present a classification system. This was then reviewed, amended and finally approved by the ESMO TR and PM WG and the ESMO leadership. / Results: This first version of the ESMO Scale of Clinical Actionability for molecular Targets (ESCAT) defines six levels of clinical evidence for molecular targets according to the implications for patient management: tier I, targets ready for implementation in routine clinical decisions; tier II, investigational targets that likely define a patient population that benefits from a targeted drug but additional data are needed; tier III, clinical benefit previously demonstrated in other tumour types or for similar molecular targets; tier IV, preclinical evidence of actionability; tier V, evidence supporting co-targeting approaches; and tier X, lack of evidence for actionability. / Conclusions: The ESCAT defines clinical evidence-based criteria to prioritise genomic alterations as markers to select patients for targeted therapies. This classification system aims to offer a common language for all the relevant stakeholders in cancer medicine and drug development

Multiple splice defects in ABCA1 cause low HDL-C in a family with Hypoalphalipoproteinemia and premature coronary disease

<p>Abstract</p> <p>Background</p> <p>Mutations at splice junctions causing exon skipping are uncommon compared to exonic mutations, and two intronic mutations causing an aberrant phenotype have rarely been reported. Despite the high number of functional <it>ABCA1 </it>mutations reported to date, splice variants have been reported infrequently. We screened DNA from a 41 year-old male with low HDL-C (12 mg/dL [0.31 mmol/L]) and a family history of premature coronary heart disease (CHD) using polymerase chain reaction single-strand conformation polymorphism (SSCP) analysis.</p> <p>Methods</p> <p>Family members with low levels of HDL-C (n = 6) were screened by SSCP for mutations in <it>ABCA1</it>. Samples with altered SSCP patterns were sequenced directly using either an ABI 3700 or ABI3730Xl DNA Analyzer. To screen for splicing defects, cDNA was isolated from the proband's RNA and was sequenced as above. A series of minigenes were constructed to determine the contribution of normal and defective alleles.</p> <p>Results</p> <p>Two novel splice variants in <it>ABCA1 </it>were identified. The first mutation was a single base pair change (T->C) in IVS 7, 6 bps downstream from the exon7/intron7 junction. Amplification of cDNA and allelic subcloning identified skipping of Exon 7 that results in the elimination of 59 amino acids from the first extracellular loop of the ABCA1 protein. The second mutation was a single base pair change (G->C) at IVS 31 -1, at the intron/exon junction of exon 32. This mutation causes skipping of exon 32, resulting in 8 novel amino acids followed by a stop codon and a predicted protein size of 1496 AA, compared to normal (2261 AA). Bioinformatic studies predicted an impact on splicing as confirmed by <it>in vitro </it>assays of constitutive splicing.</p> <p>Conclusion</p> <p>In addition to carnitine-acylcarnitine translocase (CACT) deficiency and Hermansky-Pudlak syndrome type 3, this represents only the third reported case in which 2 different splice mutations has resulted in an aberrant clinical phenotype.</p

Association of IGF1 and KDM5A polymorphisms with performance, fatness and carcass traits in chickens

Two functional and positional candidate genes were selected in a region of chicken chromosome 1 (GGA1), based on their biological roles, and also where several quantitative trait loci (QTL) have been mapped and associated with performance, fatness and carcass traits in chickens. The insulin-like growth factor 1 (IGF1) gene has been associated with several physiological functions related to growth. The lysine (K)-specific demethylase 5A (KDM5A) gene participates in the epigenetic regulation of genes involved with the cell cycle. Our objective was to find associations of selected single-nucleotide polymorphisms (SNPs) in these genes with performance, fatness and carcass traits in 165 F chickens from a resource population. In the IGF1 gene, 17 SNPs were detected, and in the KDM5A gene, nine SNPs were detected. IGF1 SNP c. 47673G > A was associated with body weight and haematocrit percentage, and also with feed intake and percentages of abdominal fat and gizzard genotype × sex interactions. KDM5A SNP c. 34208C > T genotype × sex interaction affected body weight, feed intake, percentages of abdominal fat (p = 0. 0001), carcass, gizzard and haematocrit. A strong association of the diplotype × sex interaction (p < 0. 0001) with abdominal fat was observed, and also associations with body weight, feed intake, percentages of carcass, drums and thighs, gizzard and haematocrit. Our findings suggest that the KDM5A gene might play an important role in the abdominal fat deposition in chickens. The IGF1 and KDM5A genes are strong candidates to explain the QTL mapped in this region of GGA1

An Atlas of the Speed of Copy Number Changes in Animal Gene Families and Its Implications

The notion that gene duplications generating new genes and functions is commonly accepted in evolutionary biology. However, this assumption is more speculative from theory rather than well proven in genome-wide studies. Here, we generated an atlas of the rate of copy number changes (CNCs) in all the gene families of ten animal genomes. We grouped the gene families with similar CNC dynamics into rate pattern groups (RPGs) and annotated their function using a novel bottom-up approach. By comparing CNC rate patterns, we showed that most of the species-specific CNC rates groups are formed by gene duplication rather than gene loss, and most of the changes in rates of CNCs may be the result of adaptive evolution. We also found that the functions of many RPGs match their biological significance well. Our work confirmed the role of gene duplication in generating novel phenotypes, and the results can serve as a guide for researchers to connect the phenotypic features to certain gene duplications

Accumulation of Endogenous LITAF in Aggresomes

LITAF is a 161 amino acid cellular protein which includes a proline rich N-terminus and a conserved C-terminal domain known as the simple-like domain. Mutations in LITAF have been identified in Charcot-Marie tooth disease, a disease characterized by protein aggregates. Cells transfected with cellular LITAF reveal that LITAF is localized to late endosomes/lysosomes. Here we investigated the intracellular localization of endogenous LITAF. We demonstrated that endogenous LITAF accumulates at a discrete cytoplasmic site in BGMK cells that we identify as the aggresome. To determine the domain within LITAF that is responsible for the localization of LITAF to aggresomes, we created a construct that contained the C-terminal simple-like domain of LITAF and found that this construct also localizes to aggresomes. These data suggest the simple-like domain is responsible for targeting endogenous LITAF to the aggresome

Evolution of a Bacterial Regulon Controlling Virulence and Mg2+ Homeostasis

Related organisms typically rely on orthologous regulatory proteins to respond to a given signal. However, the extent to which (or even if) the targets of shared regulatory proteins are maintained across species has remained largely unknown. This question is of particular significance in bacteria due to the widespread effects of horizontal gene transfer. Here, we address this question by investigating the regulons controlled by the DNA-binding PhoP protein, which governs virulence and Mg2+ homeostasis in several bacterial species. We establish that the ancestral PhoP protein directs largely different gene sets in ten analyzed species of the family Enterobacteriaceae, reflecting both regulation of species-specific targets and transcriptional rewiring of shared genes. The two targets directly activated by PhoP in all ten species (the most distant of which diverged >200 million years ago), and coding for the most conserved proteins are the phoPQ operon itself and the lipoprotein-encoding slyB gene, which decreases PhoP protein activity. The Mg2+-responsive PhoP protein dictates expression of Mg2+ transporters and of enzymes that modify Mg2+-binding sites in the cell envelope in most analyzed species. In contrast to the core PhoP regulon, which determines the amount of active PhoP and copes with the low Mg2+ stress, the variable members of the regulon contribute species-specific traits, a property shared with regulons controlled by dissimilar regulatory proteins and responding to different signals

Role of Duplicate Genes in Robustness against Deleterious Human Mutations

It is now widely recognized that robustness is an inherent property of biological systems [1],[2],[3]. The contribution of close sequence homologs to genetic robustness against null mutations has been previously demonstrated in simple organisms [4],[5]. In this paper we investigate in detail the contribution of gene duplicates to back-up against deleterious human mutations. Our analysis demonstrates that the functional compensation by close homologs may play an important role in human genetic disease. Genes with a 90% sequence identity homolog are about 3 times less likely to harbor known disease mutations compared to genes with remote homologs. Moreover, close duplicates affect the phenotypic consequences of deleterious mutations by making a decrease in life expectancy significantly less likely. We also demonstrate that similarity of expression profiles across tissues significantly increases the likelihood of functional compensation by homologs

Transcriptome characterization by RNA sequencing identifies a major molecular and clinical subdivision in chronic lymphocytic leukemia

Chronic lymphocytic leukemia (CLL) has heterogeneous clinical and biological behavior. Whole-genome and -exome sequencing has contributed to the characterization of the mutational spectrum of the disease, but the underlying transcriptional profile is still poorly understood. We have performed deep RNA sequencing in different subpopulations of normal B-lymphocytes and CLL cells from a cohort of 98 patients, and characterized the CLL transcriptional landscape with unprecedented resolution. We detected thousands of transcriptional elements differentially expressed between the CLL and normal B cells, including protein-coding genes, noncoding RNAs, and pseudogenes. Transposable elements are globally derepressed in CLL cells. In addition, two thousand genes—most of which are not differentially expressed—exhibit CLL-specific splicing patterns. Genes involved in metabolic pathways showed higher expression in CLL, while genes related to spliceosome, proteasome, and ribosome were among the most down-regulated in CLL. Clustering of the CLL samples according to RNA-seq derived gene expression levels unveiled two robust molecular subgroups, C1 and C2. C1/C2 subgroups and the mutational status of the immunoglobulin heavy variable (IGHV) region were the only independent variables in predicting time to treatment in a multivariate analysis with main clinico-biological features. This subdivision was validated in an independent cohort of patients monitored through DNA microarrays. Further analysis shows that B-cell receptor (BCR) activation in the microenvironment of the lymph node may be at the origin of the C1/C2 differences