A study on the minimum number of loci required for genetic evaluation using a finite locus model

For a finite locus model, Markov chain Monte Carlo (MCMC) methods can be used to estimate the conditional mean of genotypic values given phenotypes, which is also known as the best predictor (BP). When computationally feasible, this type of genetic prediction provides an elegant solution to the problem of genetic evaluation under non-additive inheritance, especially for crossbred data. Successful application of MCMC methods for genetic evaluation using finite locus models depends, among other factors, on the number of loci assumed in the model. The effect of the assumed number of loci on evaluations obtained by BP was investigated using data simulated with about 100 loci. For several small pedigrees, genetic evaluations obtained by best linear prediction (BLP) were compared to genetic evaluations obtained by BP. For BLP evaluation, used here as the standard of comparison, only the first and second moments of the joint distribution of the genotypic and phenotypic values must be known. These moments were calculated from the gene frequencies and genotypic effects used in the simulation model. BP evaluation requires the complete distribution to be known. For each model used for BP evaluation, the gene frequencies and genotypic effects, which completely specify the required distribution, were derived such that the genotypic mean, the additive variance, and the dominance variance were the same as in the simulation model. For lowly heritable traits, evaluations obtained by BP under models with up to three loci closely matched the evaluations obtained by BLP for both purebred and crossbred data. For highly heritable traits, models with up to six loci were needed to match the evaluations obtained by BLP.     

Persistence of accuracy of genomic estimated breeding values over generations in layer chickens

Background

The predictive ability of genomic estimated breeding values (GEBV) originates both from associations between high-density markers and QTL (Quantitative Trait Loci) and from pedigree information. Thus, GEBV are expected to provide more persistent accuracy over successive generations than breeding values estimated using pedigree-based methods. The objective of this study was to evaluate the accuracy of GEBV in a closed population of layer chickens and to quantify their persistence over five successive generations using marker or pedigree information.

Methods

The training data consisted of 16 traits and 777 genotyped animals from two generations of a brown-egg layer breeding line, 295 of which had individual phenotype records, while others had phenotypes on 2,738 non-genotyped relatives, or similar data accumulated over up to five generations. Validation data included phenotyped and genotyped birds from five subsequent generations (on average 306 birds/generation). Birds were genotyped for 23,356 segregating SNP. Animal models using genomic or pedigree relationship matrices and Bayesian model averaging methods were used for training analyses. Accuracy was evaluated as the correlation between EBV and phenotype in validation divided by the square root of trait heritability.

Results

Pedigree relationships in outbred populations are reduced by 50% at each meiosis, therefore accuracy is expected to decrease by the square root of 0.5 every generation, as observed for pedigree-based EBV (Estimated Breeding Values). In contrast the GEBV accuracy was more persistent, although the drop in accuracy was substantial in the first generation. Traits that were considered to be influenced by fewer QTL and to have a higher heritability maintained a higher GEBV accuracy over generations. In conclusion, GEBV capture information beyond pedigree relationships, but retraining every generation is recommended for genomic selection in closed breeding populations.     

Retracing micro-epidemics of Chagas disease using epicenter regression

Vector-borne transmission of Chagas disease has become an urban problem in the city of Arequipa, Peru, yet the debilitating symptoms that can occur in the chronic stage of the disease are rarely seen in hospitals in the city. The lack of obvious clinical disease in Arequipa has led to speculation that the local strain of the etiologic agent, Trypanosoma cruzi, has low chronic pathogenicity. The long asymptomatic period of Chagas disease leads us to an alternative hypothesis for the absence of clinical cases in Arequipa: transmission in the city may be so recent that most infected individuals have yet to progress to late stage disease. Here we describe a new method, epicenter regression, that allows us to infer the spatial and temporal history of disease transmission from a snapshot of a population's infection status. We show that in a community of Arequipa, transmission of T. cruzi by the insect vector Triatoma infestans occurred as a series of focal micro-epidemics, the oldest of which began only around 20 years ago. These micro-epidemics infected nearly 5% of the community before transmission of the parasite was disrupted through insecticide application in 2004. Most extant human infections in our study community arose over a brief period of time immediately prior to vector control. According to our findings, the symptoms of chronic Chagas disease are expected to be absent, even if the strain is pathogenic in the chronic phase of disease, given the long asymptomatic period of the disease and short history of intense transmission. Traducción al espa?ol disponible en Alternative Language Text S1/A Spanish translation of this article is available in Alternative Language Text S1.     

Arginylated calreticulin at plasma membrane increases susceptibility of cells to apoptosis

Post-translational modifications of proteins are important for the regulation of cell fate and functions; one of these post-translational modifications is arginylation. We have previously established that calreticulin (CRT), an endoplasmic reticulum resident, is also one of the arginylated substrates found in the cytoplasm. In the present study, we describe that arginylated CRT (R-CRT) binds to the cell membrane and identified its role as a preapoptotic signal. We also show that cells lacking arginyl-tRNA protein transferase are less susceptible to apoptosis than wild type cells. Under these conditions R-CRT is present on the cell membrane but at early stages is differently localized in stress granules. Moreover, cells induced to undergo apoptosis by arsenite show increased R-CRT on their cell surface. Exogenously applied R-CRT binds to the cell membrane and is able to both increase the number of cells undergoing apoptosis in wild type cells and overcome apoptosis resistance in cells lacking arginyl-tRNA protein transferase that express R-CRT on the cell surface. Thus, these results demonstrate the importance of surface R-CRT in the apoptotic response of cells, implying that post-translational arginylation of CRT can regulate its intracellular localization, cell function, and survival.     

Stage‐specific effects of bone morphogenetic proteins on the oligodendrocyte lineage

Oligodendrocyte maturation is regulated by multiple secreted factors present in the brain during critical stages of development. Whereas most of these factors promote oligodendrocyte proliferation and survival, members of the bone morphogenetic protein family (BMPs) recently have been shown to inhibit oligodendrocyte differentiation in vitro. Oligodendrocyte precursors treated with BMPs differentiate to the astrocyte lineage. Given that cells at various stages of the oligodendrocyte lineage have distinct responses to growth factors, we hypothesized that the response to BMP would be stage‐specific. Using highly purified, stage‐specific cultures, we found that BMP has distinct effects on cultured oligodendrocyte preprogenitors, precursors, and mature oligodendrocytes. Oligodendrocyte preprogenitors (PSA‐NCAM+, A2B5−) treated with BMP2 or BMP4 developed a novel astrocyte phenotype characterized by a morphological change and expression of glial fibrillary acidic protein (GFAP) but little glutamine synthetase expression and no labeling with A2B5 antibody. In contrast, treating oligodendrocyte precursors with BMPs resulted in the accumulation of cells with the traditional type 2 astrocyte phenotype (GFAP+, A2B5+). However, many of the cells with an astrocytic morphology did not express GFAP or glutamine synthetase unless thyroid hormone was present in the medium. The addition of fibroblast growth factor along with BMP to either oligodendrocyte preprogenitor or the oligodendrocyte precursor cells inhibited the switch to the astrocyte lineage, whereas platelet‐derived growth factor addition had no effect. Treatment of mature oligodendrocytes with BMP elicited no change in morphology or expression of GFAP. These data suggest that as cells progress through the oligodendrocyte lineage, they show developmentally restricted responses to the BMPs. © 2000 John Wiley & Sons, Inc. J Neurobiol 43: 1–17, 2000     

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.