Comprehensive genotype–phenotype mapping of the green fluorescent protein shows that the local fitness peak is narrow, shaped by a high prevalence of epistatic interactions, providing for the loss of fluorescence when the joint effect of mutations exceeds a threshold. Fyodor Kondrashov and colleagues report comprehensive genotype–phenotype mapping across an entire protein, based on analysis of the fitness landscape of green fluorescent protein (GFP) using a molecular barcoding and sequencing approach. They find that the fitness landscape is characterized by locally narrow regions, combined with a high prevalence of epistatic interactions, providing for the loss of fluorescence when the joint effect of mutations exceeds a threshold. Fitness landscapes1,2 depict how genotypes manifest at the phenotypic level and form the basis of our understanding of many areas of biology2,3,4,5,6,7, yet their properties remain elusive. Previous studies have analysed specific genes, often using their function as a proxy for fitness2,4, experimentally assessing the effect on function of single mutations and their combinations in a specific sequence2,5,8,9,10,11,12,13,14,15 or in different sequences2,3,5,16,17,18. However, systematic high-throughput studies of the local fitness landscape of an entire protein have not yet been reported. Here we visualize an extensive region of the local fitness landscape of the green fluorescent protein from Aequorea victoria (avGFP) by measuring the native function (fluorescence) of tens of thousands of derivative genotypes of avGFP. We show that the fitness landscape of avGFP is narrow, with 3/4 of the derivatives with a single mutation showing reduced fluorescence and half of the derivatives with four mutations being completely non-fluorescent. The narrowness is enhanced by epistasis, which was detected in up to 30% of genotypes with multiple mutations and mostly occurred through the cumulative effect of slightly deleterious mutations causing a threshold-like decrease in protein stability and a concomitant loss of fluorescence. A model of orthologous sequence divergence spanning hundreds of millions of years predicted the extent of epistasis in our data, indicating congruence between the fitness landscape properties at the local and global scales. The characterization of the local fitness landscape of avGFP has important implications for several fields including molecular evolution, population genetics and protein design.

Introduction: This study was aimed at deciphering the secretome of adipose-derived mesenchymal stromal cells (ADSCs) cultured in standard and hypoxic conditions to reveal proteins, which may be responsible for regenerative action of these cells. Methods: Human ADSCs were isolated from 10 healthy donors and cultured for 3–4 passages. Cells were serum deprived and cell purity was assessed using multiple cell surface markers. Conditioned media was collected and analyzed using LC-MS with a focus on characterizing secreted proteins. Results: Purity of the ADSC assessed as CD90+/CD73+/CD105+/CD45-/CD31- cells was greater than 99 % and viability was greater than 97 %. More than 600 secreted proteins were detected in conditioned media of ADSCs. Of these 100 proteins were common to all cultures and included key molecules involved in tissue regeneration such as collagens and collagen maturation enzymes, matrix metalloproteases, matricellular proteins, macrophage-colony stimulating factor and pigment epithelium derived factor. Common set of proteins also included molecules, which contribute to regenerative processes but were not previously associated with ADSCs. These included olfactomedin-like 3, follistatin-like 1 and prosaposin. In addition, ADSCs from the different subjects secreted proteins, which were variable between different cultures. These included proteins with neurotrophic activities, which were not previously associated with ADSCs, such as mesencephalic astrocyte-derived neurotrophic factor, meteorin and neuron derived neurotrophic factor. Hypoxia resulted in secretion of 6 proteins, the most prominent included EGF-like repeats and discoidin I-like domains 3, adrenomedullin and ribonuclease 4 of RNase A family. It also caused the disappearance of 8 proteins, including regulator of osteogenic differentiation cartilage-associated protein. Conclusions: Human ADSCs with CD90+/CD73+/CD105+/CD45-/CD31-/PDGFRβ+/NG2+/CD146+(−) immunophenotype secrete a large array of proteins, the most represented group is comprised of extracellular matrix components. Number of secreted proteins is largely unaffected by prolonged hypoxia. Variability in the secretion of several proteins from cultured ADSCs of individual subjects suggests that these cells exist as a heterogeneous population containing functionally distinct subtypes, which differ in numbers between donors.

Our ability to analyze adaptive immunity and engineer its activity has long been constrained by our limited ability to identify native pairs of heavy–light antibody chains and alpha–beta T‐cell receptor (TCR) chains — both of which comprise coupled “halves of a key”, collectively capable of recognizing specific antigens. Here, we report a cell‐based emulsion RT‐PCR approach that allows the selective fusion of the native pairs of amplified TCR alpha and beta chain genes for complex samples. A new type of PCR suppression technique was developed that makes it possible to amplify the fused library with minimal noise for subsequent analysis by high‐throughput paired‐end Illumina sequencing. With this technique, single analysis of a complex blood sample allows identification of multiple native TCR chain pairs. This approach may be extended to identify native antibody chain pairs and, more generally, pairs of mRNA molecules that are coexpressed in the same living cells.

The role of tumor-infiltrating B-cells (TIBs) and intratumorally-produced antibodies in cancer-immunity interactions essentially remains terra incognita. In particular, it remains unexplored how driver mutations could be associated with distinct TIBs signatures and their role in tumor microenvironment. Here we analyzed associations of immunoglobulin isotypes and clonality with survival in TCGA RNA-Seq data for lung adenocarcinoma (LUAD), stratifying patients into 12 driver mutation and phenotypic tumor subgroups. We revealed several unexpected associations between TIBs behavior and prognosis. Abundance and high proportion of IgG1 isotype, and low proportion of IgA among all intratumorally produced immunoglobulins were specifically associated with improved overall survival for KRASmut but not KRASwt LUAD, revealing the first link between a driver mutation and B-cell response. We found specific IgG1 signature associated with long survival, which suggests that particular specificities of IgG1+ TIBs could be beneficial in KRASmut LUAD. In contrast to our previous observations for melanoma, highly clonal IgG1 production by plasma cells had no meaningful effect on prognosis, suggesting that IgG1+ TIBs may exert a beneficial effect in KRASmut cases in an alternative way, such as efficient presentation of cognate antigens or direct B cell attack on tumor cells. Notably, a high proportion of the IgG1 isotype is positively correlated with the non-silent mutation burden both in the general LUAD cohort and in most patient subgroups, supporting a role for IgG1+ TIBs in antigen presentation. Complementing the recent finding that the presence of stromal IgG4-producing cells is associated with a favorable prognosis for patients with stage I squamous cell carcinoma, we show that the abundance of IgG4-producing TIBs likewise has a strong positive effect on overall survival in STK11mut and proximal proliferative subgroups of LUAD patients. We hypothesize that the positive role of IgG4 antibodies in some of the lung cancer subtypes could be associated with reported inability of IgG4 isotype to form immune complexes, thus preventing immunosuppression via activation of the myeloid-derived suppressor cell (MDSC) phenotype. We discover prominent and distinct associations between TIBs antibody isotypes and survival in lung adenocarcinoma carrying specific driver mutations. These findings indicate that particular types of tumor-immunity relations could be beneficial in particular driver mutation context, which should be taken into account in developing strategies of cancer immunotherapy and combination therapies. Specificity of protective B cell populations in specific cancer subgroups could become a clue to efficient targeted immunotherapies for appropriate cohorts of patients.

Cytotoxins from cobra venom are known to manifest cytotoxicity in various cell types. It is widely accepted that the plasma membrane is a target of cytotoxins, but the mechanism of their action remains obscure. Using the confocal spectral imaging technique, we show for the first time that cytotoxins from cobra venom penetrate readily into living cancer cells and accumulate markedly in lysosomes. Cytotoxins CT1 and CT2 from Naja oxiana, CT3 from Naja kaouthia and CT1 from Naja haje are demonstrated to possess this property with respect to human lung adenocarcinoma A549 and promyelocytic leukaemia HL60 cells. Immobilized plasma membrane binding accompanies the internalization of CT3 from Naja kaouthia in the HL60 cells, but it is very weak for other cytotoxins. Detectable membrane binding is not a property of any of the cytotoxins tested in A549 cells. The kinetics and concentration-dependence of cytotoxin accumulation in lysosomes correlate well with their cytotoxic effects. On the basis of the results obtained, we propose that lysosomes are a primary target of the lytic action of cytotoxins. Plasma membrane permeabilization seems to be a downstream event relative to lysosome rupture. Direct damage to the plasma membrane may be a complementary mechanism, but its relative contribution to the cytotoxic action depends on the cytotoxin structure and cell type.