Principles of metabolome conservation in animals.

Metabolite levels shape cellular physiology and disease susceptibility, yet the general principles governing metabolome evolution are largely unknown. Here, we introduce a measure of conservation of individual metabolite levels among related species. By analyzing multispecies tissue metabolome datasets in phylogenetically diverse mammals and fruit flies, we show that conservation varies extensively across metabolites. Three major functional properties, metabolite abundance, essentiality, and association with human diseases predict conservation, highlighting a striking parallel between the evolutionary forces driving metabolome and protein sequence conservation. Metabolic network simulations recapitulated these general patterns and revealed that abundant metabolites are highly conserved due to their strong coupling to key metabolic fluxes in the network. Finally, we show that biomarkers of metabolic diseases can be distinguished from other metabolites simply based on evolutionary conservation, without requiring any prior clinical knowledge. Overall, this study uncovers simple rules that govern metabolic evolution in animals and implies that most tissue metabolome differences between species are permitted, rather than favored by natural selection. More broadly, our work paves the way toward using evolutionary information to identify biomarkers, as well as to detect pathogenic metabolome alterations in individual patients.

Fully accessible fitness landscape of oncogene-negative lung adenocarcinoma.

Cancer genomes are almost invariably complex with genomic alterations cooperating during each step of carcinogenesis. In cancers that lack a single dominant oncogene mutation, cooperation between the inactivation of multiple tumor suppressor genes can drive tumor initiation and growth. Here, we shed light on how the sequential acquisition of genomic alterations generates oncogene-negative lung tumors. We couple tumor barcoding with combinatorial and multiplexed somatic genome editing to characterize the fitness landscapes of three tumor suppressor genes NF1, RASA1, and PTEN, the inactivation of which jointly drives oncogene-negative lung adenocarcinoma initiation and growth. The fitness landscape was surprisingly accessible, with each additional mutation leading to growth advantage. Furthermore, the fitness landscapes remained fully accessible across backgrounds with the inactivation of additional tumor suppressor genes. These results suggest that while predicting cancer evolution will be challenging, acquiring the multiple alterations that drive the growth of oncogene-negative tumors can be facilitated by the lack of constraints on mutational order.

Pathogen diversity drives the evolution of generalist MHC-II alleles in human populations.

Central players of the adaptive immune system are the groups of proteins encoded in the major histocompatibility complex (MHC), which shape the immune response against pathogens and tolerance to self-peptides. The corresponding genomic region is of particular interest, as it harbors more disease associations than any other region in the human genome, including associations with infectious diseases, autoimmune disorders, cancers, and neuropsychiatric diseases. Certain MHC molecules can bind to a much wider range of epitopes than others, but the functional implication of such an elevated epitope-binding repertoire has remained largely unclear. It has been suggested that by recognizing more peptide segments, such promiscuous MHC molecules promote immune response against a broader range of pathogens. If so, the geographical distribution of MHC promiscuity level should be shaped by pathogen diversity. Three lines of evidence support the hypothesis. First, we found that in pathogen-rich geographical regions, humans are more likely to carry highly promiscuous MHC class II DRB1 alleles. Second, the switch between specialist and generalist antigen presentation has occurred repeatedly and in a rapid manner during human evolution. Third, molecular positions that define promiscuity level of MHC class II molecules are especially diverse and are under positive selection in human populations. Taken together, our work indicates that pathogen load maintains generalist adaptive immune recognition, with implications for medical genetics and epidemiology.

Reading and spelling skills of prematurely born children in light of the underlying cognitive factors.

Prematurity is a serious risk factor for learning difficulties. Within the academic skills reading has the greatest impact on the prospects of the students; therefore, studying the reading skills in the risk populations is very important. The aim of our study was to investigate reading and spelling skills of prematurely born children. Our target group consisted of 8-11-year-old children (n = 23) who were born preterm with very low birthweights (VLBW). For comparison 57 full-term children (27 good readers and 30 dyslexics) were included in the study sample. To assess the reading and spelling abilities the Hungarian version of the 3DM (Dyslexia Differential Diagnosis) was used. Cognitive abilities were tested using the Hungarian adaptation of the WISC-IV and the Rey Complex Figure Test. The data were analyzed with a novel statistical approach using the R program. In the cognitive measures the mean performances of all three groups fell within the normal range. In the WISC-IV Full-scale IQ as well as in some other cognitive measures the good readers significantly outperformed both the dyslexics and the preterms. The findings of the study did not confirm our expectation that VLBW prematurity should lead to developmental disadvantages in the acquisition of reading and spelling skills since in the reading and spelling performances of the good readers and the preterms did not differ, while both the good readers and the preterms scored higher than the dyslexics. The results suggest that the cognitive assets of the preterm children contributing to their reading and spelling performances were their good spatial-visual memory, working memory, and processing speed. The identification of the cognitive mechanisms underlying reading and spelling abilities is of crucial importance for designing intervention for children with deficits in these academic skills.

Limited Evolutionary Conservation of the Phenotypic Effects of Antibiotic Resistance Mutations.

Multidrug-resistant clinical isolates are common in certain pathogens, but rare in others. This pattern may be due to the fact that mutations shaping resistance have species-specific effects. To investigate this issue, we transferred a range of resistance-conferring mutations and a full resistance gene into Escherichia coli and closely related bacteria. We found that resistance mutations in one bacterial species frequently provide no resistance, in fact even yielding drug hypersensitivity in close relatives. In depth analysis of a key gene involved in aminoglycoside resistance (trkH) indicated that preexisting mutations in other genes-intergenic epistasis-underlie such extreme differences in mutational effects between species. Finally, reconstruction of adaptive landscapes under multiple antibiotic stresses revealed that mutations frequently provide multidrug resistance or elevated drug susceptibility (i.e., collateral sensitivity) only with certain combinations of other resistance mutations. We conclude that resistance and collateral sensitivity are contingent upon the genetic makeup of the bacterial population, and such contingency could shape the long-term fate of resistant bacteria. These results underlie the importance of species-specific treatment strategies.

Total Chemical Synthesis of ISGylated-Ubiquitin Hybrid Chain Assisted by Acetamidomethyl Derivatives with Dual Functions.

Interferon-stimulated gene 15 (ISG15) is a member of the ubiquitin-like modifiers (ULM) family, which adopts a ß-grasp fold domain(s) similar to ubiquitin (Ub) with only minor sequence homology. ISG15 consists of two Ub-like domains and aids the immune system in neutralizing infections by numerous pathogens and plays an important role in defending cells against many viruses including influenza A. Recently, Ub was found to be a substrate for ISG15, which can be ISGylated on Lys29 and Lys48, while the former is more dominant. The discovery of such hybrid ISG15-Ub chains brought forward various fundamental questions regarding the nature and effect of this conjugation. To further investigate the role of hybrid ISG15-Ub chains, the pure homogeneous material of these chains is needed in workable quantities. By applying advanced chemical strategies for protein synthesis, we report the total chemical synthesis of a 231-residue ISG15-Lys29-Ub hybrid chain. During the synthesis we encountered insoluble peptide fragments, and therefore we developed a new reversible Acm based solubilizing tag to efficiently tackle this hurdle. This new Acm tag was compared with the known Arg based Acm solubilizing tag and was found to be more reliable in terms of incorporation and efficiency as demonstrated in the synthesis of the native ISG15-Ub hybrid chain.

Aspartic Acid Forming α-Ketoacid-Hydroxylamine (KAHA) Ligations with (S)-4,4-Difluoro-5-oxaproline.

The α-ketoacid-hydroxylamine (KAHA) ligation allows the coupling of unprotected peptide segments. Currently, the most applied hydroxylamine is the 5-membered cyclic hydroxylamine (S)-5-oxaproline, which forms a homoserine ester as the primary ligation product. In order to access native aspartic acid residues at the ligation site, we synthesized a 4,4-difluoro version of this monomer. Upon KAHA ligation, the resulting difluoro alcohol hydrolyzes to an aspartic acid residue with little or no formation of aspartamide. We applied this monomer for the synthesis of the hormone peptides glucagon and an insulin variant, and as well for segment ligation of the peptides UbcH5a and SUMO3.

No Evidence That Protein Noise-Induced Epigenetic Epistasis Constrains Gene Expression Evolution.

Changes in gene expression can affect phenotypes and therefore both its level and stochastic variability are frequently under selection. It has recently been proposed that epistatic interactions influence gene expression evolution: gene pairs where simultaneous knockout is more deleterious than expected should evolve reduced expression noise to avoid concurrent low expression of both proteins. In apparent support, yeast genes with many epistatic partners have low expression variation both among isogenic individuals and between species. However, the specific predictions and basic assumptions of this verbal model remain untested. Using bioinformatics analysis, we first demonstrate that the model's predictions are unsupported by available large-scale data. Based on quantitative biochemical modeling, we then show that epistasis between expression reductions (epigenetic epistasis) is not expected to aggravate the fitness cost of stochastic expression, which is in sharp contrast to the verbal argument. This nonintuitive result can be readily explained by the typical diminishing return of fitness on gene activity and by the fact that expression noise not only decreases but also increases the abundance of proteins. Overall, we conclude that stochastic variation in epistatic partners is unlikely to drive noise minimization or constrain gene expression divergence on a genomic scale.

The genomic landscape of compensatory evolution.

Adaptive evolution is generally assumed to progress through the accumulation of beneficial mutations. However, as deleterious mutations are common in natural populations, they generate a strong selection pressure to mitigate their detrimental effects through compensatory genetic changes. This process can potentially influence directions of adaptive evolution by enabling evolutionary routes that are otherwise inaccessible. Therefore, the extent to which compensatory mutations shape genomic evolution is of central importance. Here, we studied the capacity of the baker's yeast genome to compensate the complete loss of genes during evolution, and explored the long-term consequences of this process. We initiated laboratory evolutionary experiments with over 180 haploid baker's yeast genotypes, all of which initially displayed slow growth owing to the deletion of a single gene. Compensatory evolution following gene loss was rapid and pervasive: 68% of the genotypes reached near wild-type fitness through accumulation of adaptive mutations elsewhere in the genome. As compensatory mutations have associated fitness costs, genotypes with especially low fitnesses were more likely to be subjects of compensatory evolution. Genomic analysis revealed that as compensatory mutations were generally specific to the functional defect incurred, convergent evolution at the molecular level was extremely rare. Moreover, the majority of the gene expression changes due to gene deletion remained unrestored. Accordingly, compensatory evolution promoted genomic divergence of parallel evolving populations. However, these different evolutionary outcomes are not phenotypically equivalent, as they generated diverse growth phenotypes across environments. Taken together, these results indicate that gene loss initiates adaptive genomic changes that rapidly restores fitness, but this process has substantial pleiotropic effects on cellular physiology and evolvability upon environmental change. Our work also implies that gene content variation across species could be partly due to the action of compensatory evolution rather than the passive loss of genes.

Combinatorial in vivo genome editing identifies widespread epistasis during lung tumorigenesis.

Lung adenocarcinoma, the most common subtype of lung cancer, is genomically complex, with tumors containing tens to hundreds of non-synonymous mutations. However, little is understood about how genes interact with each other to enable tumorigenesis in vivo , largely due to a lack of methods for investigating genetic interactions in a high-throughput and multiplexed manner. Here, we employed a novel platform to generate tumors with all pairwise inactivation of ten tumor suppressor genes within an autochthonous mouse model of oncogenic KRAS-driven lung cancer. By quantifying the fitness of tumors with every single and double mutant genotype, we show that most tumor suppressor genetic interactions exhibited negative epistasis, with diminishing returns on tumor fitness. In contrast, Apc inactivation showed positive epistasis with the inactivation of several other genes, including dramatically synergistic effects on tumor fitness in combination with Lkb1 or Nf1 inactivation. This approach has the potential to expand the scope of genetic interactions that may be functionally characterized in vivo , which could lead to a better understanding of how complex tumor genotypes impact each step of carcinogenesis.

Fully accessible fitness landscape of oncogene-negative lung adenocarcinoma.

Cancer genomes are almost invariably complex with genomic alterations cooperating during each step of carcinogenesis. In cancers that lack a single dominant oncogene mutation, cooperation between the inactivation of multiple tumor suppressor genes can drive tumor initiation and growth. Here, we shed light on how the sequential acquisition of genomic alterations generates oncogene-negative lung tumors. We couple tumor barcoding with combinatorial and multiplexed somatic genome editing to characterize the fitness landscapes of three tumor suppressor genes NF1, RASA1, and PTEN, the inactivation of which jointly drives oncogene-negative lung adenocarcinoma initiation and growth. The fitness landscape was surprisingly accessible, with each additional mutation leading to growth advantage. Furthermore, the fitness landscapes remained fully accessible across backgrounds with additional tumor suppressor mutations. These results suggest that while predicting cancer evolution will be challenging, acquiring the multiple alterations required for the growth of oncogene-negative tumors can be facilitated by the lack of constraints on mutational order.

Combinatorial Inactivation of Tumor Suppressors Efficiently Initiates Lung Adenocarcinoma with Therapeutic Vulnerabilities.

Lung cancer is the leading cause of cancer death worldwide, with lung adenocarcinoma being the most common subtype. Many oncogenes and tumor suppressor genes are altered in this cancer type, and the discovery of oncogene mutations has led to the development of targeted therapies that have improved clinical outcomes. However, a large fraction of lung adenocarcinomas lacks mutations in known oncogenes, and the genesis and treatment of these oncogene-negative tumors remain enigmatic. Here, we perform iterative in vivo functional screens using quantitative autochthonous mouse model systems to uncover the genetic and biochemical changes that enable efficient lung tumor initiation in the absence of oncogene alterations. Generation of hundreds of diverse combinations of tumor suppressor alterations demonstrates that inactivation of suppressors of the RAS and PI3K pathways drives the development of oncogene-negative lung adenocarcinoma. Human genomic data and histology identified RAS/MAPK and PI3K pathway activation as a common feature of an event in oncogene-negative human lung adenocarcinomas. These Onc-negativeRAS/PI3K tumors and related cell lines are vulnerable to pharmacologic inhibition of these signaling axes. These results transform our understanding of this prevalent yet understudied subtype of lung adenocarcinoma. SIGNIFICANCE: To address the large fraction of lung adenocarcinomas lacking mutations in proto-oncogenes for which targeted therapies are unavailable, this work uncovers driver pathways of oncogene-negative lung adenocarcinomas and demonstrates their therapeutic vulnerabilities.

Human microRNAs co-silence in well-separated groups and have different predicted essentialities.

BACKGROUND: Short regulating RNAs guide many cellular processes. Compared with transcription factor proteins they appear to provide more specialized control and their deletions are less frequently lethal. RESULTS: We find large differences between computationally predicted lists of human microRNA (miRNA)-target pairs. Instead of integrating these lists we use the two most accurate of them. Next, we construct the co-regulation network of human miRNAs as nodes by computing the correlation (link weight) between the gene silencing scores of individual miRNAs. In this network, we locate groups of tightly co-regulating nodes (modules). Despite explicitly allowing overlaps the co-regulation modules of miRNAs are well separated. We use the modules and miRNA co-expression data to define and compute miRNA essentiality. Instead of focusing on particular biological functions we identify a miRNA as essential, if it has a low co-expression with the miRNAs in its module. This may be thought of as having many workers performing the same tasks together in one place (non-essential miRNAs) as opposed to a single worker performing those tasks alone (essential miRNA). CONCLUSIONS: On the system level, we quantitatively confirm previous findings about the specialized control provided by miRNAs. For knock-out tests we list the groups of our predicted most and least essential miRNAs. In addition, we provide possible explanations for (i) the low number of individually essential miRNAs in Caenorhabdtits elegans and (ii) the high number of ubiquitous miRNAs influencing cell and tissue-specific miRNA expression patterns in mouse and human.

[Analysis of predictive tools for further axillary involvement in patients with sentinel-lymph-node-positive, small (< or =15 mm) invasive breast cancer]. / A hónalji nyirokcsomók további érintettségére vonatkozó modellek elemzése kisméretu (< / =15 mm) ôrszemnyirokcsomó-áttétes emlôrákokban.

Small breast cancers often require different treatment than larger ones. The frequency and predictability of further nodal involvement was evaluated in patients with positive sentinel lymph nodes and breast cancers < or =15 mm by means of 8 different predictive tools. Of 506 patients with such small tumors 138 with positive sentinel nodes underwent axillary dissection and 39 of these had non-sentinel node involvement too. The Stanford nomogram and the micrometastatic nomogram were the predictive tools identifying a small group of patients with low probability of further axillary involvement that might not require completion axillary lymph node dissection. Our data also suggest that the Tenon score can separate subsets of patients with a low and a higher risk of non-sentinel node metastasis. Predictive tools based on multivariate models can help in omitting completion axillary dissection in patients with low risk of non-sentinel lymph node metastasis based on their small tumor size.

Rapid decline of bacterial drug-resistance in an antibiotic-free environment through phenotypic reversion.

Antibiotic resistance typically induces a fitness cost that shapes the fate of antibiotic-resistant bacterial populations. However, the cost of resistance can be mitigated by compensatory mutations elsewhere in the genome, and therefore the loss of resistance may proceed too slowly to be of practical importance. We present our study on the efficacy and phenotypic impact of compensatory evolution in Escherichia coli strains carrying multiple resistance mutations. We have demonstrated that drug-resistance frequently declines within 480 generations during exposure to an antibiotic-free environment. The extent of resistance loss was found to be generally antibiotic-specific, driven by mutations that reduce both resistance level and fitness costs of antibiotic-resistance mutations. We conclude that phenotypic reversion to the antibiotic-sensitive state can be mediated by the acquisition of additional mutations, while maintaining the original resistance mutations. Our study indicates that restricting antimicrobial usage could be a useful policy, but for certain antibiotics only.

The role of radiotherapy in the conservative treatment of ductal carcinoma in situ of the breast.

Breast-conserving surgery (BCS) followed by radiotherapy (RT) has become the standard of care for the treatment of early-stage (St. I-II) invasive breast carcinoma. However, controversy exists regarding the value of RT in the conservative treatment of ductal carcinoma in situ (DCIS). In this article we review the role of RT in the management of DCIS. Retrospective and prospective trials and meta-analyses published between 1975 and 2007 in the MEDLINE database, and recent issues of relevant journals/handbooks relating to DCIS, BCS and RT were searched for. In retrospective series (10,194 patients) the 10-year rate of local recurrence (LR) with and without RT was reported in the range of 9-28% and 22-54%, respectively. In four large randomised controlled trials (NSABP-B-17, EORTC-10853, UKCCCR, SweDCIS; 4,568 patients) 50 Gy whole-breast RT significantly decreased the 5-year LR rate from 16-22% (annual LR rate: 2.6-5.0%) to 7-10% (annual LR rate: 1.3-1.9%). In a recent meta-analysis of randomised trials the addition of RT to BCS resulted in a 60% risk reduction of both invasive and in situ recurrences. In a multicentre retrospective study, an additional dose of 10 Gy to the tumour bed yielded a further 55% risk reduction compared to RT without boost. To date, no subgroups have been reliably identified that do not benefit from RT after BCS. In the NSABP-B-24 trial, the addition of tamoxifen (TAM) to RT reduced ipsilateral (11.1% vs. 7.7%) and contralateral (4.9% vs. 2.3%) breast events significantly. In contrast, in the UKCCCR study, TAM produced no significant reduction in all breast events. Based on available evidence obtained from retrospective and prospective trials, all patients with DCIS have potential benefit from RT after BCS. Further prospective studies are warranted to identify subgroups of low-risk patients with DCIS for whom RT can be safely omitted. Until long-term results of ongoing studies on outcomes of patients treated with BCS alone (with or without TAM or aromatase inhibitors) are available, RT should be routinely recommended after BCS for all patients except those with contraindication.

[Combined surgery and radiotherapy in the treatment of ductal carcinoma in situ of the breast: preliminary results of the Hungarian multicenter prospective randomised study]. / In situ duktális emlokarcinóma kombinált sebészi- és sugárkezelése: a magyarországi multicentrikus prospektív.

The aim of this work is to report the preliminary results of the Hungarian multicentric randomised DCIS study. Between 2000 and 2007, 278 patients with ductal carcinoma in situ (DCIS) treated by breast-conserving surgery were randomised according to predetermined risk groups. Low/intermediate-risk patients (n=29) were randomised to 50 Gy whole-breast irradiation (WBI) or observation. High-risk cases (n=235) were allocated to receive 50 Gy WBI vs. 50 Gy WBI plus 16 Gy tumour bed boost. Very high-risk patients (patients with involved surgical margins; n=14) were randomised to 50 Gy WBI plus 16 Gy tumour bed boost or reoperation (reexcision plus radiotherapy or mastectomy alone). Immunohistochemistry (IHC) was performed to detect the expression of potential molecular prognostic markers (ER, PR, Her2, p53, Bcl-2 and Ki-67). At a median follow-up of 36 months no recurrence was observed in the low/intermediate- and very high-risk patient groups. In the high-risk group, 4 (1.7%) local recurrences and 1 (0.4%) distant metastasis occurred. No patient died of breast cancer. In the high-risk group of patients, the 3- and 5-year probability of local recurrence was 1.1% and 3.1%, respectively. The positive immunostaining for Her2 (38%), p53 (37%) and Ki-67 (44%) correlated with a high nuclear grade. Significant inverse correlation was found between the expression of ER (77%), PR (67%), Bcl-2 (64%) and grade. Preliminary results suggest that breast-conserving surgery followed by radiotherapy yields an annual local recurrence rate of less than 1% in patients with DCIS. IHC of molecular prognostic markers can assist to gain insight into the biologic heterogeneity of DCIS.

Facile folding of insulin variants bearing a prosthetic C-peptide prepared by α-ketoacid-hydroxylamine (KAHA) ligation.

The chemical synthesis of insulin is an enduring challenge due to the hydrophobic peptide chains and construction of the correct intermolecular disulfide pattern. We report a new approach to the chemical synthesis of insulin using a short, traceless, prosthetic C-peptide that facilitates the formation of the correct disulfide pattern during folding and its removal by basic treatment. The linear precursor is assembled by an ester forming α-ketoacid-hydroxylamine (KAHA) ligation that provides access to the linear insulin precursors in good yield from two readily prepared segments. This convergent and flexible route provides access to various human, mouse, and guinea pig insulins containing a single homoserine mutation that shows no detrimental effect on the biological activities.

Hsp70-associated chaperones have a critical role in buffering protein production costs.

Proteins are necessary for cellular growth. Concurrently, however, protein production has high energetic demands associated with transcription and translation. Here, we propose that activity of molecular chaperones shape protein burden, that is the fitness costs associated with expression of unneeded proteins. To test this hypothesis, we performed a genome-wide genetic interaction screen in baker's yeast. Impairment of transcription, translation, and protein folding rendered cells hypersensitive to protein burden. Specifically, deletion of specific regulators of the Hsp70-associated chaperone network increased protein burden. In agreement with expectation, temperature stress, increased mistranslation and a chemical misfolding agent all substantially enhanced protein burden. Finally, unneeded protein perturbed interactions between key components of the Hsp70-Hsp90 network involved in folding of native proteins. We conclude that specific chaperones contribute to protein burden. Our work indicates that by minimizing the damaging impact of gratuitous protein overproduction, chaperones enable tolerance to massive changes in genomic expression.