Dramatic down-regulation of oxidoreductases in human hepatocellular carcinoma hepG2 cells: proteomics and gene ontology unveiling new frontiers in cancer enzymology

Background

Oxidoreductases are enzymes that catalyze many redox reactions in normal and neoplastic cells. Their actions include catalysis of the transformation of free, neutral oxygen gas into oxygen free radicals, superoxide, hydroperoxide, singlet oxygen and hydrogen peroxide. These activated forms of oxygen contribute to oxidative stress that modifies lipids, proteins, DNA and carbohydrates. On the other hand, oxidoreductases constitute one of the most important free radical scavenger systems typified by catalase, superoxide dismutase and glutathione peroxidase. In this work, proteomics, Gene Ontology mapping and Directed Acyclic Graphs (DAG) are employed to detect and quantify differential oxidoreductase enzyme expressions between HepG2 cells and normal human liver tissues.

Results

For the set of bioinformatics calculations whose BLAST searches are performed using the BLAST program BLASTP 2.2.13 [Nov-27-2005], DAG of the Gene Ontology's Molecular Function annotations show that oxidoreductase activity parent node of the liver proteome contains 331 annotated protein sequences, 7 child nodes and an annotation score of 188.9, whereas that of HepG2 cells has 188 annotated protein sequences, 3 child nodes and an annotation score of only 91.9. Overwhelming preponderance of oxidoreductases in the liver is additionally supported by the isomerase DAGs: nearly all the reactions described in the normal liver isomerase DAG are oxidoreductase isomerization reactions, whereas only one of the three child nodes in the HepG2 isomerase DAG is oxidoreductase. Upon normalization of the annotation scores to the parent Molecular Function nodes, oxidoreductases are down-regulated in HepG2 cells by 58%. Similarly, for the set of bioinformatics calculations whose BLAST searches are carried out using BLASTP 2.2.15 [Oct-15-2006], oxidoreductases are down-regulated in HepG2 cells by 56%.

Conclusion

Proteomics and Gene Ontology reveal, for the first time, differential enzyme activities between HepG2 cells and normal human liver tissues, which may be a promising new prognostic marker of Hepatocellular carcinoma. Two independent sets of bioinformatics calculations that employ two BLAST program versions, and searched different databases, arrived at essentially the same conclusion: oxidoreductases are down-regulated in HepG2 cells by approximately 57%, when compared to normal human liver tissues. Down-regulation of oxidoreductases in hepatoma is additionally supported by Gene Ontology analysis of isomerises.     

Sample prep for proteomics of breast cancer: proteomics and gene ontology reveal dramatic differences in protein solubilization preferences of radioimmunoprecipitation assay and urea lysis buffers

Background

An important step in the proteomics of solid tumors, including breast cancer, consists of efficiently extracting most of proteins in the tumor specimen. For this purpose, Radio-Immunoprecipitation Assay (RIPA) buffer is widely employed. RIPA buffer's rapid and highly efficient cell lysis and good solubilization of a wide range of proteins is further augmented by its compatibility with protease and phosphatase inhibitors, ability to minimize non-specific protein binding leading to a lower background in immunoprecipitation, and its suitability for protein quantitation.

Results

In this work, the insoluble matter left after RIPA buffer extraction of proteins from breast tumors are subjected to another extraction step, using a urea-based buffer. It is shown that RIPA and urea lysis buffers fractionate breast tissue proteins primarily on the basis of molecular weights. The average molecular weight of proteins that dissolve exclusively in urea buffer is up to 60% higher than in RIPA. Gene Ontology (GO) and Directed Acyclic Graphs (DAG) are used to map the collective biological and biophysical attributes of the RIPA and urea proteomes. The Cellular Component and Molecular Function annotations reveal protein solubilization preferences of the buffers, especially the compartmentalization and functional distributions. It is shown that nearly all extracellular matrix proteins (ECM) in the breast tumors and matched normal tissues are found, nearly exclusively, in the urea fraction, while they are mostly insoluble in RIPA buffer. Additionally, it is demonstrated that cytoskeletal and extracellular region proteins are more soluble in urea than in RIPA, whereas for nuclear, cytoplasmic and mitochondrial proteins, RIPA buffer is preferred. Extracellular matrix proteins are highly implicated in cancer, including their proteinase-mediated degradation and remodelling, tumor development, progression, adhesion and metastasis. Thus, if they are not efficiently extracted by RIPA buffer, important information may be missed in cancer research.

Conclusion

For proteomics of solid tumors, a two-step extraction process is recommended. First, proteins in the tumor specimen should be extracted with RIPA buffer. Second, the RIPA-insoluble material should be extracted with the urea-based buffer employed in this work.     

A comparison of alternative methods to compute conditional genotype probabilities for genetic evaluation with finite locus models

An increased availability of genotypes at marker loci has prompted the development of models that include the effect of individual genes. Selection based on these models is known as marker-assisted selection (MAS). MAS is known to be efficient especially for traits that have low heritability and non-additive gene action. BLUP methodology under non-additive gene action is not feasible for large inbred or crossbred pedigrees. It is easy to incorporate non-additive gene action in a finite locus model. Under such a model, the unobservable genotypic values can be predicted using the conditional mean of the genotypic values given the data. To compute this conditional mean, conditional genotype probabilities must be computed. In this study these probabilities were computed using iterative peeling, and three Markov chain Monte Carlo (MCMC) methods – scalar Gibbs, blocking Gibbs, and a sampler that combines the Elston Stewart algorithm with iterative peeling (ESIP). The performance of these four methods was assessed using simulated data. For pedigrees with loops, iterative peeling fails to provide accurate genotype probability estimates for some pedigree members. Also, computing time is exponentially related to the number of loci in the model. For MCMC methods, a linear relationship can be maintained by sampling genotypes one locus at a time. Out of the three MCMC methods considered, ESIP, performed the best while scalar Gibbs performed the worst.     

Haplotype inference in crossbred populations without pedigree information

Background

Current methods for haplotype inference without pedigree information assume random mating populations. In animal and plant breeding, however, mating is often not random. A particular form of nonrandom mating occurs when parental individuals of opposite sex originate from distinct populations. In animal breeding this is called crossbreeding and hybridization in plant breeding. In these situations, association between marker and putative gene alleles might differ between the founding populations and origin of alleles should be accounted for in studies which estimate breeding values with marker data. The sequence of alleles from one parent constitutes one haplotype of an individual. Haplotypes thus reveal allele origin in data of crossbred individuals.

Results

We introduce a new method for haplotype inference without pedigree that allows nonrandom mating and that can use genotype data of the parental populations and of a crossbred population. The aim of the method is to estimate line origin of alleles. The method has a Bayesian set up with a Dirichlet Process as prior for the haplotypes in the two parental populations. The basic idea is that only a subset of the complete set of possible haplotypes is present in the population.

Conclusion

Line origin of approximately 95% of the alleles at heterozygous sites was assessed correctly in both simulated and real data. Comparing accuracy of haplotype frequencies inferred with the new algorithm to the accuracy of haplotype frequencies inferred with PHASE, an existing algorithm for haplotype inference, showed that the DP algorithm outperformed PHASE in situations of crossbreeding and that PHASE performed better in situations of random mating.     

You can't always get what you want: size assortative mating by mutual mate choice as a resolution of sexual conflict

Background  

Assortative mating patterns for mate quality traits like body size are often observed in nature. However, the underlying mechanisms that cause assortative mating patterns are less well known. Sexual selection is one important explanation for assortment, suggesting that i) one (usually the female) or both sexes could show preferences for mates of similar size or ii) mutual mate choice could resolve sexual conflict over quality traits into assortment. We tested these hypotheses experimentally in the socially monogamous cichlid fish Pelvicachromis taeniatus, in which mate choice is mutual.     

Genotoxic effects of two nickel-compounds in somatic cells of Drosophila melanogaster

In view of the scarcely available information on the in vivo mutagenic and co-mutagenic activity of nickel, the genotoxic potential of two nickel-compounds, nickel chloride (NiCl(2)) and nickel sulphate (NiSO(4)), was assessed in Drosophila melanogaster by measuring two different genetic endpoints. On the one hand, we used the wing-spot assay, which is based on the principle that the loss of heterozygosity of two suitable recessive markers, multiple wing hairs (mwh) and flare-3 (flr(3)), can lead to the formation of mutant clones in the imaginal disks of larval cells. On the other hand, the in vivo comet assay, which detects single- and double-strand DNA breaks, was also used with larval haemocytes. These cells offer several advantages: they are highly sensitive to genotoxic agents, the sampling and processing methodologies are quite simple and the level of basal DNA damage is relatively low. No significant increases in the frequencies of the three categories of mutant spots (i.e. small single spots, large single spots, and twin spots) were observed in the wing-spot assay; however, NiSO(4) induced significant dose-dependent increases in DNA damage in the comet assay. In addition, the combined treatments with gamma-radiation and NiCl(2) and NiSO(4) showed a slight but significant increase in the frequency of the three categories of mutant spots compared with the frequency induced by gamma-radiation alone, indicating that both nickel compounds have a synergistic interaction. These results support the assumption that both nickel compounds could act as co-mutagens interfering with DNA-repair processes and that the in vivo comet assay is a sensitive and effective method for detecting the DNA damage induced by NiSO(4) in haemocytes of D. melanogaster.     

Genome-wide association and genomic prediction for host response to porcine reproductive and respiratory syndrome virus infection

Background

Host genetics has been shown to play a role in porcine reproductive and respiratory syndrome (PRRS), which is the most economically important disease in the swine industry. A region on Sus scrofa chromosome (SSC) 4 has been previously reported to have a strong association with serum viremia and weight gain in pigs experimentally infected with the PRRS virus (PRRSV). The objective here was to identify haplotypes associated with the favorable phenotype, investigate additional genomic regions associated with host response to PRRSV, and to determine the predictive ability of genomic estimated breeding values (GEBV) based on the SSC4 region and based on the rest of the genome. Phenotypic data and 60 K SNP genotypes from eight trials of ~200 pigs from different commercial crosses were used to address these objectives.

Results

Across the eight trials, heritability estimates were 0.44 and 0.29 for viral load (VL, area under the curve of log-transformed serum viremia from 0 to 21 days post infection) and weight gain to 42 days post infection (WG), respectively. Genomic regions associated with VL were identified on chromosomes 4, X, and 1. Genomic regions associated with WG were identified on chromosomes 4, 5, and 7. Apart from the SSC4 region, the regions associated with these two traits each explained less than 3% of the genetic variance. Due to the strong linkage disequilibrium in the SSC4 region, only 19 unique haplotypes were identified across all populations, of which four were associated with the favorable phenotype. Through cross-validation, accuracies of EBV based on the SSC4 region were high (0.55), while the rest of the genome had little predictive ability across populations (0.09).

Conclusions

Traits associated with response to PRRSV infection in growing pigs are largely controlled by genomic regions with relatively small effects, with the exception of SSC4. Accuracies of EBV based on the SSC4 region were high compared to the rest of the genome. These results show that selection for the SSC4 region could potentially reduce the effects of PRRS in growing pigs, ultimately reducing the economic impact of this disease.     

Analysis of genomic damage in the mutagen-sensitive mus-201 mutant of Drosophila melanogaster by arbitrarily primed PCR (AP-PCR) fingerprinting.

DNA repair mechanisms are important to maintain the stability of the genome. In Drosophila melanogaster, the mus-201 gene is required in the excision repair process. To study the contribution of the mus-201 gene in the stability of the Drosophila genome, we have used the arbitrarily primed PCR fingerprinting method (AP-PCR). We have analysed the changes in the genomic DNA fingerprints from the progeny of wild-type males crossed with mus-201 repair-deficient or repair-proficient females. After induction of DNA damage with 2-acetylaminofluorene (2-AAF) in the wild-type parental males, quantitative and qualitative differences in the AP-PCR fingerprints were detected between the two crosses, and the estimate of the genomic damage detected by AP-PCR has clearly shown that the mus-201 repair deficiency is associated with an increase of genomic damage. The predominant type of alterations detected by AP-PCR under the mus-201 repair-deficient conditions agree with the results obtained in microsatellite PCR analysis, suggesting that the role of the mus-201 gene, necessary in excision repair, is not associated to the mismatch repair process. The work reported here demonstrates that the AP-PCR is a suitable technique to analyse genetic alterations in D. melanogaster and, consequently, can be used to compare the susceptibility to genomic damage of different DNA repair mutants.     

Apoptosis,cytokine and chemokine induction by non-structural 1 (NS1) proteins encoded by different influenza subtypes

Background

Influenza pandemic remains a serious threat to human health. Viruses of avian origin, H5N1, H7N7 and H9N2, have repeatedly crossed the species barrier to infect humans. Recently, a novel strain originated from swine has evolved to a pandemic. This study aims at improving our understanding on the pathogenic mechanism of influenza viruses, in particular the role of non-structural (NS1) protein in inducing pro-inflammatory and apoptotic responses.

Methods

Human lung epithelial cells (NCI-H292) was used as an in-vitro model to study cytokine/chemokine production and apoptosis induced by transfection of NS1 mRNA encoded by seven infleunza subtypes (seasonal and pandemic H1, H2, H3, H5, H7, and H9), respectively.

Results

The results showed that CXCL-10/IP10 was most prominently induced (> 1000 folds) and IL-6 was slightly induced (< 10 folds) by all subtypes. A subtype-dependent pattern was observed for CCL-2/MCP-1, CCL3/MIP-1α, CCL-5/RANTES and CXCL-9/MIG; where induction by H5N1 was much higher than all other subtypes examined. All subtypes induced a similar temporal profile of apoptosis following transfection. The level of apoptosis induced by H5N1 was remarkably higher than all others. The cytokine/chemokine and apoptosis inducing ability of the 2009 pandemic H1N1 was similar to previous seasonal strains.

Conclusions

In conclusion, the NS1 protein encoded by H5N1 carries a remarkably different property as compared to other avian and human subtypes, and is one of the keys to its high pathogenicity. NCI-H292 cells system proves to be a good in-vitro model to delineate the property of NS1 proteins.