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.     

Mitochondrial DNA haplogrouping of the brown bear,Ursus arctos (Carnivora: Ursidae) in Asia,based on a newly developed APLP analysis

Sequence analyses of the complete brown bear, Ursus arctos, mitochondrial DNA (mtDNA) genome have detected scattered single nucleotide polymorphisms (SNPs) that define distinct mtDNA haplogroups in phylogeographical studies. The degraded DNA in historical samples, such as stuffed or excavated specimens, however, is often not suitable for sequence analyses. To address this problem, we developed an amplified product length polymorphism (APLP) analysis for mtDNA‐haplogrouping U. arctos specimens by detecting haplogroup‐specific SNPs. We verified the validity and utility of this method by analysing up to 170‐year‐old skin samples from U. arctos specimens collected widely across continental Eurasia. We detected some of the same haplogroups as those occurring in eastern Hokkaido (Japan) and eastern Alaska in continental Eurasia (the Altai and the Caucasus). Our results show that U. arctos in eastern Hokkaido and eastern Alaska descended from a common ancestor in continental Eurasia, and suggest that U. arctos occupied several refugia in southern Asia during the Last Glacial Maximum. © 2014 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 111 , 627–635.     

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.     

Artificial Brassica napus flowering in Bangladesh

Summary Natural rapeseed (Brassica napus L.; AACC 2n=38), originated in the temperate climate of the Southwest European Mediterranean region, fails to complete its generative phase in the subtropics and is thus not cultivated in countries like Bangladesh. Adapted agroecotypes are available from the diploid representatives of its genome A (B. campestris/pekinensis, 2n=20) and C (B. oleracea/alboglabra, 2n=18). An artificial resynthesis based on carefully selected progenitor lines was expected to give a photoperiodically better adapted rapeseed. B. pekinensis x B. oleracea/alboglabra gave 2 hybrids and 87 matromorphous plants from 1,448 crossed flowers and the reciprocal combination gave no hybrid but 11 matromorphous plants from 2,228 pollinated flowers. The two true hybrids were vegetatively propagated and chromosome doubled. Part of the F₂ was grown in Sweden (all plants flowered and the most early ones were selected), part in Bangladesh (13 out of 706 plants flowered). The selected F₃ material flowered in Bangladesh and transgressions in earliness could be recorded, some lines were of definite agronomic potential. A correlation in earliness between reaction in Sweden (long day) and Bangladesh (short day) was observed.     

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.