Multiomic Metabolic Enrichment Network Analysis Reveals Metabolite–Protein Physical Interaction Subnetworks Altered in Cancer

Metabolism is recognized as an important driver of cancer progression and other complex diseases, but global metabolite profiling remains a challenge. Protein expression profiling is often a poor proxy since existing pathway enrichment models provide an incomplete mapping between the proteome and metabolism. To overcome these gaps, we introduce multiomic metabolic enrichment network analysis (MOMENTA), an integrative multiomic data analysis framework for more accurately deducing metabolic pathway changes from proteomics data alone in a gene set analysis context by leveraging protein interaction networks to extend annotated metabolic models. We apply MOMENTA to proteomic data from diverse cancer cell lines and human tumors to demonstrate its utility at revealing variation in metabolic pathway activity across cancer types, which we verify using independent metabolomics measurements. The novel metabolic networks we uncover in breast cancer and other tumors are linked to clinical outcomes, underscoring the pathophysiological relevance of the findings.     

A comparative Proteomics Analysis Identified Differentially Expressed Proteins in Pancreatic Cancer–Associated Stellate Cell Small Extracellular Vesicles

Human pancreatic stellate cells (HPSCs) are an essential stromal component and mediators of pancreatic ductal adenocarcinoma (PDAC) progression. Small extracellular vesicles (sEVs) are membrane-enclosed nanoparticles involved in cell-to-cell communications and are released from stromal cells within PDAC. A detailed comparison of sEVs from normal pancreatic stellate cells (HPaStec) and from PDAC-associated stellate cells (HPSCs) remains a gap in our current knowledge regarding stellate cells and PDAC. We hypothesized there would be differences in sEVs secretion and protein expression that might contribute to PDAC biology. To test this hypothesis, we isolated sEVs using ultracentrifugation followed by characterization by electron microscopy and Nanoparticle Tracking Analysis. We report here our initial observations. First, HPSC cells derived from PDAC tumors secrete a higher volume of sEVs when compared to normal pancreatic stellate cells (HPaStec). Although our data revealed that both normal and tumor-derived sEVs demonstrated no significant biological effect on cancer cells, we observed efficient uptake of sEVs by both normal and cancer epithelial cells. Additionally, intact membrane-associated proteins on sEVs were essential for efficient uptake. We then compared sEV proteins isolated from HPSCs and HPaStecs cells using liquid chromatography–tandem mass spectrometry. Most of the 1481 protein groups identified were shared with the exosome database, ExoCarta. Eighty-seven protein groups were differentially expressed (selected by 2-fold difference and adjusted p value ≤0.05) between HPSC and HPaStec sEVs. Of note, HPSC sEVs contained dramatically more CSE1L (chromosome segregation 1–like protein), a described marker of poor prognosis in patients with pancreatic cancer. Based on our results, we have demonstrated unique populations of sEVs originating from stromal cells with PDAC and suggest that these are significant to cancer biology. Further studies should be undertaken to gain a deeper understanding that could drive novel therapy.     

Microsecond Dynamics During the Binding-induced Folding of an Intrinsically Disordered Protein

Tau is an intrinsically disordered protein implicated in many neurodegenerative diseases. The repeat domain fragment of tau, tau-K18, is known to undergo a disorder to order transition in the presence of lipid micelles and vesicles, in which helices form in each of the repeat domains. Here, the mechanism of helical structure formation, induced by a phospholipid mimetic, sodium dodecyl sulfate (SDS) at sub-micellar concentrations, has been studied using multiple biophysical probes. A study of the conformational dynamics of the disordered state, using photoinduced electron transfer coupled to fluorescence correlation spectroscopy (PET-FCS) has indicated the presence of an intermediate state, I, in equilibrium with the unfolded state, U. The cooperative binding of the ligand (L), SDS, to I has been shown to induce the formation of a compact, helical intermediate (IL₅) within the dead time (∼37 µs) of a continuous flow mixer. Quantitative analysis of the PET-FCS data and the ensemble microsecond kinetic data, suggests that the mechanism of induction of helical structure can be described by a U ↔ I ↔ IL₅ ↔ FL₅ mechanism, in which the final helical state, FL₅, forms from IL₅ with a time constant of 50–200 µs. Finally, it has been shown that the helical conformation is an aggregation-competent state that can directly form amyloid fibrils.     

Lysine Succinylation of VBS Contributes to Sclerotia Development and Aflatoxin Biosynthesis in Aspergillus flavus

Aspergillus flavus is a common saprophytic and pathogenic fungus, and its secondary metabolic pathways are one of the most highly characterized owing to its aflatoxin (AF) metabolite affecting global economic crops and human health. Different natural environments can cause significant variations in AF synthesis. Succinylation was recently identified as one of the most critical regulatory post-translational modifications affecting metabolic pathways. It is primarily reported in human cells and bacteria with few studies on fungi. Proteomic quantification of lysine succinylation (Ksuc) exploring its potential involvement in secondary metabolism regulation (including AF production) has not been performed under natural conditions in A. flavus. In this study, a quantification method was performed based on tandem mass tag labeling and antibody-based affinity enrichment of succinylated peptides via high accuracy nano-liquid chromatography with tandem mass spectrometry to explore the succinylation mechanism affecting the pathogenicity of naturally isolated A. flavus strains with varying toxin production. Altogether, 1240 Ksuc sites in 768 proteins were identified with 1103 sites in 685 proteins quantified. Comparing succinylated protein levels between high and low AF-producing A. flavus strains, bioinformatics analysis indicated that most succinylated proteins located in the AF biosynthetic pathway were downregulated, which directly affected AF synthesis. Versicolorin B synthase is a key catalytic enzyme for heterochrome B synthesis during AF synthesis. Site-directed mutagenesis and biochemical studies revealed that versicolorin B synthase succinylation is an important regulatory mechanism affecting sclerotia development and AF biosynthesis in A. flavus. In summary, our quantitative study of the lysine succinylome in high/low AF-producing strains revealed the role of Ksuc in regulating AF biosynthesis. We revealed novel insights into the metabolism of AF biosynthesis using naturally isolated A. flavus strains and identified a rich source of metabolism-related enzymes regulated by succinylation.     

Developmental biochemistry of cottonseed embryogenesis and germination. XIX. Sequences and genomic organization of the α globulin (vicilin) genes of cottonseed

The globulin storage protein genes of cotton are found to exist as gene tandems that contain a gene from each of the 2 globulin subfamilies separated by a spacer region of about 2700 or 3400 base pairs. Three different tandems have been identified by restriction endonuclease mapping of genomic DNA. A cDNA that is different from the genes of the tandems in map sites and/or in nucleotide sequence indicates that a fourth tandem probably exists in the cotton genome. Since the species of cotton used here (Gossypium hirsutum) is an amphidiploid, it is likely that two of the tandems are contributed from each genome.Considerable divergence in nucleotide sequence (18%) and in derived amino acid sequence (28%) is found when the 2 genes of a sequenced tandem are compared. The sequence of the cDNA closely resembles one of the genes in the tandem showing only a 4% divergence in nucleotides and a 4.2% divergence in amino acids. Thus the 2 genes of each tandem represent a relatively ancient gene duplication that has given rise to the two globulin subfamilies of cotton. Only one subfamily has a glycosylation site and the glycosylation of its derived proteins gives rise to the 2 molecular weight sets of globulins seen on gel electrophoresis.Other basic features of these genes and their derived proteins are presented.     

A BamHI family of highly repeated DNA sequences of Nicotiana tabacum

Summary HRS60.1, a monomer unit (184 bp) of a highly repeated nuclear DNA sequence of Nicotiana tabacum, has been cloned and sequenced. Following BamHI digestion of tobacco DNA, Southern hybridization with HRS60.1 revealed a ladder of hybridization bands corresponding to multiples of the basic monomer unit. If the tobacco DNA was digested with restriction endonucleases which have no target site in HRS60.1, the larger part of DNA homologous to HRS60.1 remained as uncleaved relic DNA. These results suggest a tandem arrangement of this DNA repeat unit. Four other clones of tobacco nuclear DNA cross-hybridized with HRS60.1, thus forming a HRS60-family. Sequencing their inserts has shown their strong mutual homology. HRS60-family comprised about 2% of the nuclear genome of N. tabacum. Computer comparisons with other tandem plant-repeated DNA sequences could not detect any other homologous sequence.     

An immunodominant site of γ-zein1 is in the region of tandem hexapeptide repeats

The immunochemical data from studies with polyclonal antisera to -zein₁, the 27 kD component of the maize prolamin, indicated that the region containing 8 tandem repeats of the sequence PPPVHL is an immunodominant site. In one case, the entire antibody repertoire of an antiserum recognized epitope(s) within this region. Three 17-mer oligopeptides corresponding to the predicted antigenic epitopes of -zein₁ were synthesized and reacted with three different anti--zein₁ sera in order to map antigenic sites in the intact protein. These antisera yielded positive reactions with a 17-mer peptide (peptide 37), which was not in a hydrophilic maximum but derived from the repeat region. The same antisera gave little or no reaction with other peptides (peptides 38 and 39), both of which were in a hydrophilic maximum. In addition, an antiserum to peptide 37 reacted strongly with both the homologous antigen and the intact -zein₁. Peptide 37 also blocked the binding of antisera to -zein₁ in competition assays. Subsequently, the shorter 6-mer (peptide 82) and 12-mer (peptide 80) versions of peptide 37 were synthesized, and both reacted with anti-peptide 37 serum and also with each of the three anti--zein₁ sera. In these reactions and in competition assays, the reactivity and the blocking ability increased in proportion to the length of the peptide. Based on these data, it was concluded that the repeat region of -zein₁ is the site of one or more continuous immunodominant epitopes. The data also suggest that the repeat region is exposed on the surface of the folded protein and probably occur as a mobile, random coil.     

Detection of linkage between marker loci and loci affecting quantitative traits in crosses between segregating populations

Summary By making use of pedigree information and information on marker-genotypes of the parent and F-1 individuals crossed to form an F-2 population, it is possible to carry out a linkage analysis between marker loci and loci affecting quantitative traits in a cross between segregating parent populations that are at fixation for alternative alleles at the QTL, but share the same alleles at the marker loci. For two-allele systems, depending on marker allele frequencies in the parent populations, 2–4 times as many F-2 offspring will have to be raised and scored for markers and quantitative traits in order to provide power equivalent to that obtained in a cross between fully inbred lines. Major savings in number of F-2 offspring raised can be achieved by scoring each parent pair for a large number of markers in each chromosomal region and scoring F-1 and F-2 offspring only for those markers for which the parents were homozygous for alternative alleles. For multiple allele systems, particularly when dealing with hypervariable loci, only 10%–20% additional F-2 offspring will have to be raised and scored to provide power equivalent to that obtained in a cross between inbred lines. When a resource population contains novel favorable alleles at quantitative trait loci that are not present (or rare) in a commercial population, analyses of this sort will enable the loci of interest to be identified, mapped and manipulated effectively in breeding programs.Contribution no. 2124-E, 1987 series from The Agricultural Research Organization, The Volcani Center, Bet Dagan, Israel