Why highly expressed proteins evolve slowly

Much recent work has explored molecular and population-genetic constraints on the rate of protein sequence evolution. The best predictor of evolutionary rate is expression level, for reasons which have remained unexplained. Here, we hypothesize that selection to reduce the burden of protein misfolding will favor protein sequences with increased robustness to translational missense errors. Pressure for translational robustness increases with expression level and constrains sequence evolution. Using several sequenced yeast genomes, global expression and protein abundance data, and sets of paralogs traceable to an ancient whole-genome duplication in yeast, we rule out several confounding effects and show that expression level explains roughly half the variation in Saccharomyces cerevisiae protein evolutionary rates. We examine causes for expression's dominant role and find that genome-wide tests favor the translational robustness explanation over existing hypotheses that invoke constraints on function or translational efficiency. Our results suggest that proteins evolve at rates largely unrelated to their functions, and can explain why highly expressed proteins evolve slowly across the tree of life.Comment: 40 pages, 3 figures, with supporting informatio

An organoid biobank for childhood kidney cancers that captures disease and tissue heterogeneity

Kidney tumours are among the most common solid tumours in children, comprising distinct subtypes differing in many aspects, including cell-of-origin, genetics, and pathology. Pre-clinical cell models capturing the disease heterogeneity are currently lacking. Here, we describe the first paediatric cancer organoid biobank. It contains tumour and matching normal kidney organoids from over 50 children with different subtypes of kidney cancer, including Wilms tumours, malignant rhabdoid tumours, renal cell carcinomas, and congenital mesoblastic nephromas. Paediatric kidney tumour organoids retain key properties of native tumours, useful for revealing patient-specific drug sensitivities. Using single cell RNA-sequencing and high resolution 3D imaging, we further demonstrate that organoid cultures derived from Wilms tumours consist of multiple different cell types, including epithelial, stromal and blastemal-like cells. Our organoid biobank captures the heterogeneity of paediatric kidney tumours, providing a representative collection of well-characterised models for basic cancer research, drug-screening and personalised medicine

Exome sequencing identifies novel AD-associated genes

The genetic component of Alzheimer’s disease (AD) has been mainly assessed using Genome Wide Association Studies (GWAS), which do not capture the risk contributed by rare variants. Here, we compared the gene-based burden of rare damaging variants in exome sequencing data from 32,558 individuals —16,036 AD cases and 16,522 controls— in a two-stage analysis. Next to known genes TREM2, SORL1 and ABCA7, we observed a significant association of rare, predicted damaging variants in ATP8B4 and ABCA1 with AD risk, and a suggestive signal in ADAM10. Next to these genes, the rare variant burden in RIN3, CLU, ZCWPW1 and ACE highlighted these genes as potential driver genes in AD-GWAS loci. Rare damaging variants in these genes, and in particular loss-of-function variants, have a large effect on AD-risk, and they are enriched in early onset AD cases. The newly identified AD-associated genes provide additional evidence for a major role for APP-processing, Aβ-aggregation, lipid metabolism and microglial function in AD

BRCA1-mutated and basal-like breast cancers have similar aCGH profiles and a high incidence of protein truncating TP53 mutations

<p>Abstract</p> <p>Background</p> <p>Basal-like breast cancers (BLBC) are aggressive breast cancers for which, so far, no targeted therapy is available because they typically lack expression of hormone receptors and HER2. Phenotypic features of BLBCs, such as clinical presentation and early age of onset, resemble those of breast tumors from <it>BRCA1</it>-mutation carriers. The genomic instability of <it>BRCA1</it>-mutated tumors can be effectively targeted with DNA-damaging agents and poly-(ADP-ribose) polymerase 1 (PARP1) inhibitors. Molecular similarities between BLBCs and <it>BRCA1</it>-mutated tumors may therefore provide predictive markers for therapeutic response of BLBCs.</p> <p>Methods</p> <p>There are several known molecular features characteristic for <it>BRCA1</it>-mutated breast tumors: 1) increased numbers of genomic aberrations, 2) a distinct pattern of genomic aberrations, 3) a high frequency of <it>TP53 </it>mutations and 4) a high incidence of complex, protein-truncating <it>TP53 </it>mutations. We compared the frequency of <it>TP53 </it>mutations and the pattern and amount of genomic aberrations between <it>BRCA1</it>-mutated breast tumors, BLBCs and luminal breast tumors by <it>TP53 </it>gene sequencing and array-based comparative genomics hybridization (aCGH) analysis.</p> <p>Results</p> <p>We found that the high incidence of protein truncating <it>TP53 </it>mutations and the pattern and amount of genomic aberrations specific for BRCA1-mutated breast tumors are also characteristic for BLBCs and different from luminal breast tumors.</p> <p>Conclusions</p> <p>Complex, protein truncating TP53 mutations in BRCA1-mutated tumors may be a direct consequence of genomic instability caused by BRCA1 loss, therefore, the presence of these types of TP53 mutations in sporadic BLBCs might be a hallmark of BRCAness and a potential biomarker for sensitivity to PARP inhibition. Also, our data suggest that a small subset of genomic regions may be used to identify BRCA1-like BLBCs. BLBCs share molecular features that were previously found to be specific for BRCA1-mutated breast tumors. These features might be useful for the identification of tumors with increased sensitivity to (high-dose or dose-dense) alkylating agents and PARP inhibitors.</p

Preservation of Ranking Order in the Expression of Human Housekeeping Genes

Housekeeping (HK) genes fulfill the basic needs for a cell to survive and function properly. Their ubiquitous expression, originally thought to be constant, can vary from tissue to tissue, but this variation remains largely uncharacterized and it could not be explained by previously identified properties of HK genes such as short gene length and high GC content. By analyzing microarray expression data for human genes, we uncovered a previously unnoted characteristic of HK gene expression, namely that the ranking order of their expression levels tends to be preserved from one tissue to another. Further analysis by tensor product decomposition and pathway stratification identified three main factors of the observed ranking preservation, namely that, compared to those of non-HK (NHK) genes, the expression levels of HK genes show a greater degree of dispersion (less overlap), stableness (a smaller variation in expression between tissues), and correlation of expression. Our results shed light on regulatory mechanisms of HK gene expression that are probably different for different HK genes or pathways, but are consistent and coordinated in different tissues

Human brain harbors single nucleotide somatic variations in functionally relevant genes possibly mediated by oxidative stress

Somatic variation in DNA can cause cells to deviate from the preordained genomic path in both disease and healthy conditions. Here, using exome sequencing of paired tissue samples, we show that the normal human brain harbors somatic single base variations measuring up to 0.48% of the total variations. Interestingly, about 64% of these somatic variations in the brain are expected to lead to non-synonymous changes, and as much as 87% of these represent G:C>T:A transversion events. Further, the transversion events in the brain were mostly found in the frontal cortex, whereas the corpus callosum from the same individuals harbors the reference genotype. We found a significantly higher amount of 8-OHdG (oxidative stress marker) in the frontal cortex compared to the corpus callosum of the same subjects (p<0.01), correlating with the higher G:C>T:A transversions in the cortex. We found significant enrichment for axon guidance and related pathways for genes harbouring somatic variations. This could represent either a directed selection of genetic variations in these pathways or increased susceptibility of some loci towards oxidative stress. This study highlights that oxidative stress possibly influence single nucleotide somatic variations in normal human brain

A Fine-Structure Map of Spontaneous Mitotic Crossovers in the Yeast Saccharomyces cerevisiae

Homologous recombination is an important mechanism for the repair of DNA damage in mitotically dividing cells. Mitotic crossovers between homologues with heterozygous alleles can produce two homozygous daughter cells (loss of heterozygosity), whereas crossovers between repeated genes on non-homologous chromosomes can result in translocations. Using a genetic system that allows selection of daughter cells that contain the reciprocal products of mitotic crossing over, we mapped crossovers and gene conversion events at a resolution of about 4 kb in a 120-kb region of chromosome V of Saccharomyces cerevisiae. The gene conversion tracts associated with mitotic crossovers are much longer (averaging about 12 kb) than the conversion tracts associated with meiotic recombination and are non-randomly distributed along the chromosome. In addition, about 40% of the conversion events have patterns of marker segregation that are most simply explained as reflecting the repair of a chromosome that was broken in G1 of the cell cycle

Patterns and Mechanisms of Ancestral Histone Protein Inheritance in Budding Yeast

Tracking of ancestral histone proteins over multiple generations of genome replication in yeast reveals that old histones move along genes from 3′ toward 5′ over time, and that maternal histones move up to around 400 bp during genomic replication

Systematical Detection of Significant Genes in Microarray Data by Incorporating Gene Interaction Relationship in Biological Systems

Many methods, including parametric, nonparametric, and Bayesian methods, have been used for detecting differentially expressed genes based on the assumption that biological systems are linear, which ignores the nonlinear characteristics of most biological systems. More importantly, those methods do not simultaneously consider means, variances, and high moments, resulting in relatively high false positive rate. To overcome the limitations, the SWang test is proposed to determine differentially expressed genes according to the equality of distributions between case and control. Our method not only latently incorporates functional relationships among genes to consider nonlinear biological system but also considers the mean, variance, skewness, and kurtosis of expression profiles simultaneously. To illustrate biological significance of high moments, we construct a nonlinear gene interaction model, demonstrating that skewness and kurtosis could contain useful information of function association among genes in microarrays. Simulations and real microarray results show that false positive rate of SWang is lower than currently popular methods (T-test, F-test, SAM, and Fold-change) with much higher statistical power. Additionally, SWang can uniquely detect significant genes in real microarray data with imperceptible differential expression but higher variety in kurtosis and skewness. Those identified genes were confirmed with previous published literature or RT-PCR experiments performed in our lab