Extracellular volatilomic alterations induced by hypoxia in breast cancer cells.

INTRODUCTION: The metabolic shift induced by hypoxia in cancer cells has not been explored at volatilomic level so far. The volatile organic metabolites (VOMs) constitute an important part of the metabolome and their investigation could provide us crucial aspects of hypoxia driven metabolic reconfiguration in cancer cells. OBJECTIVE: To identify the altered volatilomic response induced by hypoxia in metastatic/aggressive breast cancer (BC) cells. METHODS: BC cells were cultured under normoxic and hypoxic conditions and VOMs were extracted using HS-SPME approach and profiled by standard GC-MS system. Univariate and multivariate statistical approaches (p < 0.05, Log2 FC ≥ 0.58/≤ - 0.58, PC1 > 0.13/< - 0.13) were applied to select the VOMs differentially altered after hypoxic treatment. Metabolic pathway analysis was also carried out in order to identify altered metabolic pathways induced by the hypoxia in the selected BC cells. RESULTS: Overall, 20 VOMs were found to be significantly altered (p < 0.05, PC1 > 0.13/< - 0.13) upon hypoxic exposure to BC cells. Further, cell line specific volatilomic alterations were extracted by comparative metabolic analysis of aggressive (MDA-MB-231) vs. non-aggressive (MCF-7) cells incubated under hypoxia and normoxia. In this case, 15 and 12 VOMs each were found to be significantly altered in aggressive cells when exposed to hypoxic and normoxic condition respectively. Out of these, 9 VOMs were found to be uniquely associated with hypoxia, 6 were specific to normoxia and 6 were found common to both the conditions. Formic acid was identified as the most prominent molecule with higher abundance levels in aggressive as compared to non-aggressive cells in both conditions. Furthermore, metabolic pathway analyses revealed that fatty acid biosynthesis and nicotinate and nicotinamide metabolism were significantly altered in aggressive as compared to non-aggressive cells in normoxia and hypoxia respectively. CONCLUSIONS: Higher formate overflow was observed in aggressive cells compared to non-aggressive cells incubated under both the conditions, reinforcing its correlation with aggressive and invasive cancer type. Moreover, under hypoxia, aggressive cells preferred to be bioenergetically more efficient whereas, under normoxia, fatty acid biosynthesis was favoured when compared to non-aggressive cells.

Volatilomic insight of head and neck cancer via the effects observed on saliva metabolites.

Head and neck cancer (HNC) is a heterogeneous malignant disease with distinct global distribution. Metabolic adaptations of HNC are significantly gaining clinical interests nowadays. Here, we investigated effects of HNC on differential expression of volatile metabolites in human saliva. We applied headspace solid phase microextraction coupled with gas chromatography-mass spectrometry analysis of saliva samples collected from 59 human subjects (HNC - 32, Control - 27). We identified and quantified 48 volatile organic metabolites (VOMs) and observed profound effects of HNC on these metabolites. These effects were VOM specific and significantly differed in the biologically comparable healthy controls. HNC induced changes in salivary VOM composition were well attributed to in vivo metabolic effects. A panel of 15 VOMs with variable importance in projection (VIP) score >1, false discovery rate (FDR) corrected p-value < 0.05 and log2 fold change (log2 FC) value of ≥0.58/≤-0.58 were regarded as discriminatory metabolites of pathophysiological importance. Afterwards, receiver operator characteristic curve (ROC) projected certain VOMs viz., 1,4-dichlorobenzene, 1,2-decanediol, 2,5-bis1,1-dimethylethylphenol and E-3-decen-2-ol with profound metabolic effects of HNC and highest class segregation potential. Moreover, metabolic pathways analysis portrayed several dysregulated pathways in HNC, which enhanced our basic understanding on salivary VOM changes. Our observations could redefine several known/already investigated systemic phenomenons (e.g. biochemical pathways). These findings will inspire further research in this direction and may open unconventional avenues for non-invasive monitoring of HNC and its therapy in the future.

A non-invasive approach to explore the discriminatory potential of the urinary volatilome of invasive ductal carcinoma of the breast.

Worldwide, breast invasive ductal carcinoma (IDC) accounts for the majority of the reported cases of this form of cancer. IDC effective management, as for any form of cancer, would greatly benefit from early diagnosis. This, however, due to various socio-economic reasons, is very far for the reality in developing countries like India, where cancer diagnosis is often carried out at late stages when disease management is troublesome. With the present work, we aim to evaluate a simple analytical methodology to identify a set of volatile organic compounds (VOCs) in urine samples, as a biosignature for IDC. Using solid-phase microextraction followed by gas chromatography/mass spectrometry, a panel of 14 urinary VOCs was found to discriminate IDC (n = 65) from a healthy control (HC) group (n = 70) through multivariate statistical treatments. Furthermore, metabolic pathway analysis revealed various dysregulated pathways involved in IDC patients hinting that their detailed investigations could lead to novel mechanistic insights into the disease pathophysiology. In addition, we validated the expression pattern of five of these VOCs namely 2-ethyl-1-hexanol, isolongifolenone, furan, dodecanoic acid, 2-methoxy-phenol in another external cohort of 59 urinary samples (IDC = 32 and HC = 27) and found their expression pattern to be consistent with the primary sample set. To our knowledge, this is the first study exploring breast IDC volatome alterations in Indian patients.

Interaction and localization diversities of global and local hubs in human protein-protein interaction networks.

Hubs, the highly connected nodes in protein-protein interaction networks (PPINs), are associated with several characteristic properties and are known to perform vital roles in cells. We defined two classes of hubs, global (housekeeping) and local (tissue-specific) hubs. These two categories of hubs are distinct from each other with respect to their abundance, structure and function. However, how distinct are the spatial expression pattern and other characteristics of their interacting partners is still not known. Our investigations revealed that the partners of the local hubs compared with those of global hubs are conserved across the tissues in which they are expressed. Partners of local hubs show diverse subcellular localizations as compared with the partners of global hubs. We examined the nature of interacting domains in both categories of hubs and found that they are promiscuous in global hubs but not so in local hubs. Deletion of some of the local and global hubs has an impact on the characteristic path length of the network indicating that those hubs are inter-modular in nature. Our present study has, therefore, shed further light on the characteristic features of the local and global hubs in human PPIN. This knowledge of different topological aspects of hubs with regard to their types and subtypes is essential as it helps in better understanding of roles of hub proteins in various cellular processes under various conditions including those caused by host-pathogen interactions and therefore useful in prioritizing targets for drug design and repositioning.

Spectrum of mutations in Glutaryl-CoA dehydrogenase gene in glutaric aciduria type I--Study from South India.

BACKGROUND: Glutaric aciduria type I is an autosomal recessive organic acid disorder. The primary defect is the deficiency of Glutaryl-CoA dehydrogenase (EC number 1.3.99.7) enzyme that is involved in the catabolic pathways of the amino acids l-lysine, l-hydroxylysine, and l-tryptophan. It is a treatable neuro-metabolic disorder. Early diagnosis and treatment helps in preventing brain damage. METHODS: The Glutaryl-CoA dehydrogenase gene (GCDH) gene was sequenced to identify disease causing mutations by direct sequencing of all the exons in twelve patients who were biochemically confirmed with GA I. RESULTS: We identified eleven mutations of which nine are homozygous mutations, one heterozygous and two synonymous mutations. Among the eleven mutations, four mutations p.Q162R, p.P286S, p.W225X in two families and p.V410M are novel. A milder clinical presentation is observed in those families who are either heterozygous or with a benign synonymous SNP. Multiple sequence alignment (MSA) of GCDH with its homologues revealed that the observed novel mutations are not tolerated by protein structure and function. CONCLUSIONS: The present study indicates genetic heterogeneity in GCDH gene mutations among South Indian population. Genetic analysis is useful in prenatal diagnosis and prevention. Mutation analysis is a useful tool in the absence of non-availability of enzyme assay in GA I.

Novel mutations of the arylsulphatase B (ARSB) gene in Indian patients with mucopolysaccharidosis type VI.

BACKGROUND & OBJECTIVES: Mucopolysaccharidosis type VI (MPS VI) is a rare, autosomal recessive lysosomal storage disorder caused by deficient enzymatic activity of N-acetyl galactosamine-4-sulphatase resulting from mutations in the arylsulphatase B (ARSB) gene. The ARSB gene is located on chromosome 5q11-q13 and is composed of eight exons. More than hundred ARSB mutations have been reported so far, but the mutation spectrum of MPS VI in India is still unknown. Hence, the aim of the present study was to identify the mutational spectrum in patients with MPS VI in India and to study the genotype-phenotype association and functional outcomes of these mutations. METHODS: Molecular characterization of the ARSB gene by Sanger sequencing was done for 15 patients (aged 15 months to 11 yr) who were enzymatically confirmed to have MPS VI. Age of onset, clinical progression and enzyme activity levels in each patient were studied to look for genotype-phenotype association. Haplotype analysis performed for unrelated patients with the recurring mutation W450C, was suggestive of a founder effect. Sequence and structural analyses of the ARSB protein using standard software were carried out to determine the impact of detected mutations on the function of the ARSB protein. RESULTS: A total of 12 mutations were identified, of which nine were novel mutations namely, p.D53N, p.L98R, p.Y103SfsX9, p.W353X, p.H393R, p.F166fsX18, p.I220fsX5, p.W450L, and p.W450C, and three were known mutations (p.D54N, p.A237D and p.S320R). The nine novel sequence variants were confirmed not to be polymorphic variants by performing sequencing in 50 unaffected individuals from the same ethnic population. INTERPRETATION & CONCLUSIONS: Nine novel mutations were identified in MPS VI cases from India in the present study. The study also provides some insights into the genotype-phenotype association in MPS VI.

Recurrent and novel GLB1 mutations in India.

GM1 gangliosidosis is a lysosomal storage disorder caused by mutations in the GLB1 gene, leading to the deficiency of the enzyme ß-d-galactosidase. In this study, we report molecular findings in 50 Asian Indian families with GM1 gangliosidosis. We sequenced all the exons and flanking intronic sequences of GLB1 gene. We identified 33 different mutations (20 novel and 13 previously reported). The novel mutations include 12 missense (p.M1?, p.E129Q, p.G134R, p.L236P, p.G262E, p.L297F, p.Y331C, p.G414V, p.K493N, p.L514P, p.P597L, p.T600I), four splicing (c.246-2A>G, c.397-2A>G, c.552+1G>T, c.956-2A>G), three indels (p.R22Qfs*8, p.L24Cfs*47, p.I489Qfs*4) and one nonsense mutation (p.Q452*). Most common mutations identified in this study were c.75+2InsT (14%) and p.L337P (10%). Known mutations accounted for 67% of allele frequency in our cohort of patients, suggesting that these mutations in GLB1 are recurrent across different populations. Twenty three mutations were localized in the TIM barrel domain, ß-domain 1 and ß-domain 2. In silico sequence and structure analysis of GLB1 reveal that all the novel mutations affect the function and structure of the protein. We hereby report on the largest series of patients with GM1 gangliosidosis and the first from India.

GALNS mutations in Indian patients with mucopolysaccharidosis IVA.

Mucopolysaccharidosis IV A (Morquio syndrome A, MPS IVA) is a lysosomal storage disease caused by the deficiency of N-acetylgalactosamine-6-sulfatase (GALNS). The mutation spectrum in this condition is yet to be determined in Indians. We aimed to analyze the mutations in the GALNS gene in Asian Indians with MPS IVA. All the exons and the adjacent intronic regions of the gene were amplified and sequenced in sixty-eight unrelated Indian families. We identified 136 mutant alleles comprising of 40 different mutations. We report twenty-two novel mutations that comprise of seventeen missense (p.Asn32Thr, p.Leu36Arg, p.Pro52Leu, p.Pro77Ser, p.Cys79Arg, p.His142Pro, p.Tyr191Asp, p.Asn204Thr, p.Gly188Ser, p.Phe216Ser, p.Trp230Cys, p.Ala291Ser, p.Gly317Arg, p.His329Pro, p.Arg386Ser, p.Glu450Gly, p.Cys501Ser), three splice-site variants (c.120+1G>C, c.1003-3C>G, c.1139+1G>A), one nonsense mutation (p.Gln414*) and one frameshift mutation (p.Pro420Leufs*440). Eighteen mutations have been reported earlier. Among these p.Ser287Leu (8.82%), p.Phe216Ser (7.35%), p.Asn32Thr (6.61%) and p.Ala291Ser (5.88%) were the most frequent mutations in Indian patients but were rare in the mutational profiles reported in other populations. These results indicate that the Indian patients may have a distinct mutation spectrum compared to those of other populations. Mutant alleles in exon 1, 7 and 8 accounted for 44.8% of the mutations, and sequencing of these exons initially may be a cost-effective approach in Asian Indian patients. This is the largest study on molecular analysis of patients with MPS IVA reported in the literature, and the first report from India.

Distal arthrogryposis type 5D with a novel ECEL1 gene mutation.

Distal arthrogryposis syndromes (DAs) show wide clinical variability and overlapping clinical findings with the other DAs classified by Bamshad et al. [1996]. Most of the DAs are inherited as autosomal dominant disorders. DA type 5D is a subtype of DA type 5 inherited as autosomal recessive disorder, clinically characterized by congenital distal joint contractures, knee extension contractures, congenital hip dislocation, club foot, ptosis and other eye findings, furrowed tongue, and scoliosis. Here, we report on a family with clinical features of DA type 5D with novel mutations in the ECEL1 gene.

Viral proteins that bridge unconnected proteins and components in the human PPI network.

Viruses, despite having small genomes and few proteins, make an array of interactions with host proteins as they solely depend on host machinery for their replication and reproduction. Hence, analysis of the Human-Virus Protein-Protein Interaction Network (Hu-Vir PPI network) helps us to gain certain insights into the molecular mechanisms underlying the hijacking of host cell machinery by viruses for their perpetuation. Here we report an analysis of the Human-Virus Bridged PPI Networks that has led us to identify viral articulation points (VAPs) which connect unconnected components of the Human-PPI (Hu-PPI) network. VAPs cross-link peripheral nodes to the giant component of the Hu-PPI network. VAPs interact with a number of relatively lower topologically central human proteins and are conserved among related viruses. The linked nodes comprise of those that are mostly expressed during viral infection, as well as those that are found exclusively in some metabolic pathways, indicating that the novel viral mediation of certain human protein-protein interactions may form the basis for virus-specific tuning of the host machinery. The functional importance of VAPs and their interaction partners in virus replication make them potential drug targets against viral infection. Our investigations also led to the discovery of an example of a Human Endogenous Retrovirus (HERV) encoded protein, syncytin, as an Articulation Point (AP) in the Hu-PPI network, suggesting that VAPs may be retained in a genome if they result in any beneficial function in the host.

Molecular analyses of novel ASAH1 mutations causing Farber lipogranulomatosis: analyses of exonic splicing enhancer inactivating mutation.

Farber lipogranulomatosis is a rare autosomal recessive lysosomal storage disorder caused by mutations in the ASAH1 gene. In the largest ever study, we identified and characterized ASAH1 mutations from 11 independent Farber disease (FD) families. A total of 13 different mutations were identified including 1 splice, 1 polypyrimidine tract (PPT) deletion and 11 missense mutations. Eleven mutations were exclusive to the Indian population. The IVS6+4A>G splice and IVS5-16delTTTTC PPT deletion mutations resulted in skipping of exon 6 precluding thereby the region responsible for cleavage of enzyme precursor. A missense mutation (p.V198A) resulted in skipping of exon 8 due to inactivation of an exonic splicing enhancer (ESE) element. This is the first report of mutations affecting PPT and ESE in the ASAH1 gene resulting in FD.

Molecular genetic analysis of MSUD from India reveals mutations causing altered protein truncation affecting the C-termini of E1α and E1ß.

Maple Syrup Urine Disease is a rare metabolic disorder caused by reduced/absent activity of the branched chain α-Ketoacid dehydrogenase enzyme complex. Mutations in BCKDHA, BCKDHB, and DBT, that encode important subunits of the enzyme complex namely E1α, E1ß, and E2, are the primary cause for the disease. We have performed the first molecular genetic analysis of MSUD from India on nine patients exhibiting classical MSUD symptoms. BCKDHA and BCKDHB mutations were identified in four and five patients, respectively including seven novel mutations namely the BCKDHA c.1249delC, c.1312T>C, and c.1561T>A and the BCKDHB c.401T>A, c.548G>A, c.964A>G, and c.1065delT. The BCKDHB c.970C>T (p.R324X) mutation was shown to trigger nonsense mediated decay-based degradation of the transcript. Seven of the total 11 mutations resulted in perturbations in the E1α or E1ß C-termini either through altered termination or through an amino acid change; these are expected to result in disruption of E1 enzyme complex assembly. Our study has therefore revealed that BCKDHA and BCKDHB mutations might be primarily responsible for MSUD in the Indian population.

A study on mutational dynamics of simple sequence repeats in relation to mismatch repair system in prokaryotic genomes.

Mutational bias toward expansion or contraction of simple sequence repeats (SSRs) is referred to as directionality of SSR evolution. In this communication, we report the mutational bias exhibited by mononucleotide SSRs occurring in the non-coding regions of several prokaryotic genomes. Our investigations revealed that the strains or species lacking mismatch repair (MMR) system generally show higher number of polymorphic SSRs than those species/strains having MMR system. An exception to this observation was seen in the mycobacterial genomes that are MMR deficient where only a few SSR tracts were seen with mutations. This low incidence of SSR mutations even in the MMR-deficient background could be attributed to the high fidelity of the DNA polymerases as a consequence of high generation time of the mycobacteria. MMR system-deficient species generally did not show any bias toward mononucleotide SSR expansions or contractions indicating a neutral evolution of SSRs in these species. The MMR-proficient species in which the observed mutations correspond to secondary mutations showed bias toward contraction of polymononucleotide tracts, perhaps, indicating low efficiency of MMR system to repair SSR-induced slippage errors on template strands. This bias toward deletion in the mononucleotide SSR tracts might be a probable reason behind scarcity for long poly A|T and G|C tracts in prokaryotic systems which are mostly MMR proficient. In conclusion, our study clearly demonstrates mutational dynamics of SSRs in relation to the presence/absence of MMR system in the prokaryotic system.

A founder ectodysplasin A receptor (EDAR) mutation results in a high frequency of the autosomal recessive form of hypohidrotic ectodermal dysplasia in India.

BACKGROUND: Hypohidrotic/anhidrotic ectodermal dysplasia (HED) is a rare Mendelian disorder affecting ectodermal tissues. The disease is primarily caused by inactivation of any one of three genes, namely ectodysplasin A1 (EDA-A1), which encodes a ligand belonging to the tumour necrosis factor (TNF) superfamily; ectodysplasin A receptor (EDAR), encoding the EDA-A1 receptor and ectodysplasin A receptor-associated death domain (EDARADD), encoding an adaptor protein. X-linked recessive (EDA-A1), the predominant form of HED, as well as autosomal recessive and dominant (EDAR and EDARADD) inheritance patterns have been identified in affected families. OBJECTIVES: To determine the common genes causing HED in India. METHODS: We performed mutation analysis on 26 HED families from India (including 30 patients). In addition, we carried out sequence and structural analysis of missense/nonsense and insertion/deletion mutations. RESULTS: Among the 26 families analysed, disease-causing EDAR mutations were identified in 12 (46%) while EDA-A1 mutations were detected in 11 (42%). Four novel mutations in EDAR and five in EDA-A1 were identified. More importantly, a possible founder EDAR mutation, namely c.1144G>A, was identified in five independent families, thus accounting for about one-fifth of affected families in whom mutation was detected. A majority of EDA-A1 mutations localized to the TNF-like domain while the location of EDAR mutations was more widespread. CONCLUSIONS: This is the first report of a founder EDAR mutation and of a significantly high frequency of autosomal recessive HED.

G-IMEx: A comprehensive software tool for detection of microsatellites from genome sequences.

Microsatellites are ubiquitous short tandem repeats found in all known genomes and are known to play a very important role in various studies and fields including DNA fingerprinting, paternity studies, evolutionary studies, virulence and adaptation of certain bacteria and viruses etc. Due to the sequencing of several genomes and the availability of enormous amounts of sequence data during the past few years, computational studies of microsatellites are of interest for many researchers. In this context, we developed a software tool called Imperfect Microsatellite Extractor (IMEx), to extract perfect, imperfect and compound microsatellites from genome sequences along with their complete statistics. Recently we developed a user-friendly graphical-interface using JAVA for IMEx to be used as a stand-alone software named G-IMEx. G-IMEx takes a nucleotide sequence as an input and the results are produced in both html and text formats. The Linux version of G-IMEx can be downloaded for free from http://www.cdfd.org.in/imex.

Profilin oligomerization and its effect on poly (L-proline) binding and phosphorylation.

Profilin is a cytoskeletal protein that interacts specifically with actin, phosphoinositides and poly (l-proline). Experimental results and in silico studies revealed that profilin exists as dimer and tetramer. Profilin oligomers possess weak affinity to poly (l-proline) due to unavailability of binding sites in dimers and tetramers. Phosphorylation studies indicate that profilin dimers are not phosphorylated while teramers are preferentially phosphorylated over monomers. In silico studies revealed that PKC phosphorylation site, S137 is buried in dimer while it is accessible in tetramer.

Support Vector Machine-based classification of protein folds using the structural properties of amino acid residues and amino acid residue pairs.

MOTIVATION: Fold recognition is a key step in the protein structure discovery process, especially when traditional sequence comparison methods fail to yield convincing structural homologies. Although many methods have been developed for protein fold recognition, their accuracies remain low. This can be attributed to insufficient exploitation of fold discriminatory features. RESULTS: We have developed a new method for protein fold recognition using structural information of amino acid residues and amino acid residue pairs. Since protein fold recognition can be treated as a protein fold classification problem, we have developed a Support Vector Machine (SVM) based classifier approach that uses secondary structural state and solvent accessibility state frequencies of amino acids and amino acid pairs as feature vectors. Among the individual properties examined secondary structural state frequencies of amino acids gave an overall accuracy of 65.2% for fold discrimination, which is better than the accuracy by any method reported so far in the literature. Combination of secondary structural state frequencies with solvent accessibility state frequencies of amino acids and amino acid pairs further improved the fold discrimination accuracy to more than 70%, which is approximately 8% higher than the best available method. In this study we have also tested, for the first time, an all-together multi-class method known as Crammer and Singer method for protein fold classification. Our studies reveal that the three multi-class classification methods, namely one versus all, one versus one and Crammer and Singer method, yield similar predictions. AVAILABILITY: Dataset and stand-alone program are available upon request.

Escherichia coli RNA polymerase mutations located near the upstream edge of an RNA:DNA hybrid and the beginning of the RNA-exit channel are defective for transcription antitermination by the N protein from lambdoid phage H-19B.

Transcription antitermination is an important mechanism that can control regulation of gene expression. The N protein of lambdoid phages modifies the transcription elongation complex (EC) and helps it to overcome downstream terminators. In this modified EC, the C-terminal domain of N makes specific interactions with RNA polymerase (RNAP). The interacting surface of RNAP for N is unknown. Here, we report five mutations in the beta (G1045D) and beta' (P251S, P254L, R270C and G336S) subunits of RNAP that are specifically defective for antitermination by N protein of the lambdoid phage, H-19B. A mutation in the C-terminal domain of N, L108F, suppresses the defect of beta'-P254L. Purified mutant holoenzymes exhibit less processive antitermination. The amino acid substitutions in the mutant RNAPs cluster very close to the RNA:DNA hybrid at the beginning of the RNA-exit channel of the EC. We suggest that the action of H-19B N is exerted through the region defined by these amino acids. Wild-type N stabilizes the EC at terminator sites and in this modified EC a part of the terminator hairpin may form but appears to be unstable. We propose that the action of N close to the active center alters the RNAP-nucleic acid interactions around the RNA:DNA hybrid, which impairs proper folding of the terminator hairpin or stabilizes the weak RNA:DNA hybrid, or both.

SilkSatDb: a microsatellite database of the silkworm, Bombyx mori.

The SilkSatDb (silkmoth microsatellite database) (http://www.cdfd.org.in/silksatdb) is a relational database of microsatellites extracted from the available expressed sequence tags and whole genome shotgun sequences of the silkmoth, Bombyx mori. The database has been rendered with a simple and robust web-based search facility, developed using PHP. The SilkSatDb also stores information on primers developed and validated in the laboratory. Users can retrieve information on the microsatellite and the protocols used, along with informative figures and polymorphism status of those microsatellites. In addition, the interface is coupled with Autoprimer, a primer-designing program, using which users can design primers for the loci of interest.