Combining RT-PCR-seq and RNA-seq to catalog all genic elements encoded in the human genome.

Within the ENCODE Consortium, GENCODE aimed to accurately annotate all protein-coding genes, pseudogenes, and noncoding transcribed loci in the human genome through manual curation and computational methods. Annotated transcript structures were assessed, and less well-supported loci were systematically, experimentally validated. Predicted exon-exon junctions were evaluated by RT-PCR amplification followed by highly multiplexed sequencing readout, a method we called RT-PCR-seq. Seventy-nine percent of all assessed junctions are confirmed by this evaluation procedure, demonstrating the high quality of the GENCODE gene set. RT-PCR-seq was also efficient to screen gene models predicted using the Human Body Map (HBM) RNA-seq data. We validated 73% of these predictions, thus confirming 1168 novel genes, mostly noncoding, which will further complement the GENCODE annotation. Our novel experimental validation pipeline is extremely sensitive, far more than unbiased transcriptome profiling through RNA sequencing, which is becoming the norm. For example, exon-exon junctions unique to GENCODE annotated transcripts are five times more likely to be corroborated with our targeted approach than with extensive large human transcriptome profiling. Data sets such as the HBM and ENCODE RNA-seq data fail sampling of low-expressed transcripts. Our RT-PCR-seq targeted approach also has the advantage of identifying novel exons of known genes, as we discovered unannotated exons in ~11% of assessed introns. We thus estimate that at least 18% of known loci have yet-unannotated exons. Our work demonstrates that the cataloging of all of the genic elements encoded in the human genome will necessitate a coordinated effort between unbiased and targeted approaches, like RNA-seq and RT-PCR-seq.

Regulation of human class I alcohol dehydrogenases by bile acids

Class I alcohol dehydrogenases (ADH1s) are the rate-limiting enzymes for ethanol and vitamin A (retinol) metabolism in the liver . Because previous studies have shown that human ADH1 enzymes may participate in bile acid metabolism, we investigated whether the bile acid-activated nuclear receptor farnesoid X receptor (FXR) regulates ADH1 genes. In human hepatocytes, both the endogenous FXR ligand chenodeoxycholic acid and synthetic FXR-specific agonist GW4064 increased ADH1 mRNA, protein, and activity. Moreover, overexpression of a constitutively active form of FXR induced ADH1A and ADH1B expression, whereas silencing of FXR abolished the effects of FXR agonists on ADH1 expression and activity. Transient transfection studies and electrophoretic mobility shift assays revealed functional FXR response elements in the ADH1A and ADH1B proximal promoters, thus indicating that both genes are direct targets of FXR. These findings provide the first evidence for direct connection of bile acid signaling and alcohol metabolism.

Revisiting spatial distribution and biochemical composition of calcium-containing crystals in human osteoarthritic articular cartilage.

INTRODUCTION: Calcium-containing (CaC) crystals, including basic calcium phosphate (BCP) and calcium pyrophosphate dihydrate (CPP), are associated with destructive forms of osteoarthritis (OA). We assessed their distribution and biochemical and morphologic features in human knee OA cartilage. METHODS: We prospectively included 20 patients who underwent total knee replacement (TKR) for primary OA. CaC crystal characterization and identification involved Fourier-transform infra-red spectrometry and scanning electron microscopy of 8 to 10 cartilage zones of each knee, including medial and lateral femoral condyles and tibial plateaux and the intercondyle zone. Differential expression of genes involved in the mineralization process between cartilage with and without calcification was assessed in samples from 8 different patients by RT-PCR. Immunohistochemistry and histology studies were performed in 6 different patients. RESULTS: Mean (SEM) age and body mass index of patients at the time of TKR was 74.6 (1.7) years and 28.1 (1.6) kg/m², respectively. Preoperative X-rays showed joint calcifications (chondrocalcinosis) in 4 cases only. The medial femoro-tibial compartment was the most severely affected in all cases, and mean (SEM) Kellgren-Lawrence score was 3.8 (0.1). All 20 OA cartilages showed CaC crystals. The mineral content represented 7.7% (8.1%) of the cartilage weight. All patients showed BCP crystals, which were associated with CPP crystals for 8 joints. CaC crystals were present in all knee joint compartments and in a mean of 4.6 (1.7) of the 8 studied areas. Crystal content was similar between superficial and deep layers and between medial and femoral compartments. BCP samples showed spherical structures, typical of biological apatite, and CPP samples showed rod-shaped or cubic structures. The expression of several genes involved in mineralization, including human homolog of progressive ankylosis, plasma-cell-membrane glycoprotein 1 and tissue-nonspecific alkaline phosphatase, was upregulated in OA chondrocytes isolated from CaC crystal-containing cartilages. CONCLUSIONS: CaC crystal deposition is a widespread phenomenon in human OA articular cartilage involving the entire knee cartilage including macroscopically normal and less weight-bearing zones. Cartilage calcification is associated with altered expression of genes involved in the mineralisation process.

Comparative transcriptomics of extreme phenotypes of human HIV-1 infection and SIV infection in sooty mangabey and rhesus macaque.

High levels of HIV-1 replication during the chronic phase of infection usually correlate with rapid progression to severe immunodeficiency. However, a minority of highly viremic individuals remains asymptomatic and maintains high CD4+ T cell counts. This tolerant profile is poorly understood and reminiscent of the widely studied nonprogressive disease model of SIV infection in natural hosts. Here, we identify transcriptome differences between rapid progressors (RPs) and viremic nonprogressors (VNPs) and highlight several genes relevant for the understanding of HIV-1-induced immunosuppression. RPs were characterized by a specific transcriptome profile of CD4+ and CD8+ T cells similar to that observed in pathogenic SIV-infected rhesus macaques. In contrast, VNPs exhibited lower expression of interferon-stimulated genes and shared a common gene regulation profile with nonpathogenic SIV-infected sooty mangabeys. A short list of genes associated with VNP, including CASP1, CD38, LAG3, TNFSF13B, SOCS1, and EEF1D, showed significant correlation with time to disease progression when evaluated in an independent set of CD4+ T cell expression data. This work characterizes 2 minimally studied clinical patterns of progression to AIDS, whose analysis may inform our understanding of HIV pathogenesis.

Tandem chimerism as a means to increase protein complexity in the human genome.

The "one-gene, one-protein" rule, coined by Beadle and Tatum, has been fundamental to molecular biology. The rule implies that the genetic complexity of an organism depends essentially on its gene number. The discovery, however, that alternative gene splicing and transcription are widespread phenomena dramatically altered our understanding of the genetic complexity of higher eukaryotic organisms; in these, a limited number of genes may potentially encode a much larger number of proteins. Here we investigate yet another phenomenon that may contribute to generate additional protein diversity. Indeed, by relying on both computational and experimental analysis, we estimate that at least 4%-5% of the tandem gene pairs in the human genome can be eventually transcribed into a single RNA sequence encoding a putative chimeric protein. While the functional significance of most of these chimeric transcripts remains to be determined, we provide strong evidence that this phenomenon does not correspond to mere technical artifacts and that it is a common mechanism with the potential of generating hundreds of additional proteins in the human genome.

Human disease and drug pharmacology, complex as real life

In the past decades drug discovery practice has escaped from the complexity of the formerly used phenotypic screening in animals to focus on assessing drug effects on isolated protein targets in the search for drugs that exclusively and potently hit one selected target, thought to be critical for a given disease, while not affecting at all any other target to avoid the occurrence of side-effects. However, reality does not conform to these expectations, and, conversely, this approach has been concurrent with increased attrition figures in late-stage clinical trials, precisely due to lack of efficacy and safety. In this context, a network biology perspective of human disease and treatment has burst into the drug discovery scenario to bring it back to the consideration of the complexity of living organisms and particularly of the (patho)physiological environment where protein targets are (mal)functioning and where drugs have to exert their restoring action. Under this perspective, it has been found that usually there is not one but several disease-causing genes and, therefore, not one but several relevant protein targets to be hit, which do not work on isolation but in a highly interconnected manner, and that most known drugs are inherently promiscuous. In this light, the rationale behind the currently prevailing single-target-based drug discovery approach might even seem a Utopia, while, conversely, the notion that the complexity of human disease must be tackled with complex polypharmacological therapeutic interventions constitutes a difficult-torefuse argument that is spurring the development of multitarget therapies.

The human Rad52 protein exists as a heptameric ring.

The RAD52 epistasis group was identified in yeast as a group of genes required to repair DNA damaged by ionizing radiation [1]. Genetic evidence indicates that Rad52 functions in Rad51-dependent and Rad51-independent recombination pathways [2] [3] [4]. Consistent with this, purified yeast and human Rad52 proteins have been shown to promote single-strand DNA annealing [5] [6] [7] and to stimulate Rad51-mediated homologous pairing [8] [9] [10] [11]. Electron microscopic examinations of the yeast [12] and human [13] Rad52 proteins have revealed their assembly into ring-like structures in vitro. Using both conventional transmission electron microscopy and scanning transmission electron microscopy (STEM), we found that the human Rad52 protein forms heptameric rings. A three-dimensional (3D) reconstruction revealed that the heptamer has a large central channel. Like the hexameric helicases such as Escherichia coli DnaB [14] [15], bacteriophage T7 gp4b [16] [17], simian virus 40 (SV40) large T antigen [18] and papilloma virus E1 [19], the Rad52 rings show a distinctly chiral arrangement of subunits. Thus, the structures formed by the hexameric helicases may be a more general property of other proteins involved in DNA metabolism, including those, such as Rad52, that do not bind and hydrolyze ATP.

Atorvastatin-loaded hydrogel affects the smooth muscle cells of human veins.

Intimal hyperplasia (IH) is the major cause of stenosis of vein grafts. Drugs such as statins prevent stenosis, but their systemic administration has limited effects. We developed a hyaluronic acid hydrogel matrix, which ensures a controlled release of atorvastatin (ATV) at the site of injury. The release kinetics demonstrated that 100% of ATV was released over 10 hours, independent of the loading concentration of the hydrogel. We investigated the effects of such a delivery on primary vascular smooth muscle cells isolated from human veins. ATV decreased the proliferation, migration, and passage of human smooth muscle cells (HSMCs) across a matrix barrier in a similar dose-dependent (5-10 µM) and time-dependent manner (24-72 hours), whether the drug was directly added to the culture medium or released from the hydrogel. Expression analysis of genes known to be involved in the development of IH demonstrated that the transcripts of both the gap junction protein connexin43 (Cx43) and plasminogen activator inhibitor-1 (PAI-1) were decreased after a 24-48-hour exposure to the hydrogel loaded with ATV, whereas the transcripts of the heme oxygenase (HO-1) and the inhibitor of tissue plasminogen activator were increased. At the protein level, Cx43, PAI-1, and metalloproteinase-9 expression were decreased, whereas HO-1 was upregulated in the presence of ATV. The data demonstrate that ATV released from a hydrogel has effects on HSMCs similar to the drug being freely dissolved in the environment.

Gene expression profiling reveals consistent differences between clinical samples of human leukaemias and their model cell lines.

Microarray gene expression profiles of fresh clinical samples of chronic myeloid leukaemia in chronic phase, acute promyelocytic leukaemia and acute monocytic leukaemia were compared with profiles from cell lines representing the corresponding types of leukaemia (K562, NB4, HL60). In a hierarchical clustering analysis, all clinical samples clustered separately from the cell lines, regardless of leukaemic subtype. Gene ontology analysis showed that cell lines chiefly overexpressed genes related to macromolecular metabolism, whereas in clinical samples genes related to the immune response were abundantly expressed. These findings must be taken into consideration when conclusions from cell line-based studies are extrapolated to patients.

The genomic distribution of intraspecific and interspecific sequence divergence of human segmental duplications relative to human/chimpanzee chromosomal rearrangements

Background: It has been suggested that chromosomal rearrangements harbor the molecular footprint of the biological phenomena which they induce, in the form, for instance, of changes in the sequence divergence rates of linked genes. So far, all the studies of these potential associations have focused on the relationship between structural changes and the rates of evolution of single-copy DNA and have tried to exclude segmental duplications (SDs). This is paradoxical, since SDs are one of the primary forces driving the evolution of structure and function in our genomes and have been linked not only with novel genes acquiring new functions, but also with overall higher DNA sequence divergence and major chromosomal rearrangements.Results: Here we take the opposite view and focus on SDs. We analyze several of the features of SDs, including the rates of intraspecific divergence between paralogous copies of human SDs and of interspecific divergence between human SDs and chimpanzee DNA. We study how divergence measures relate to chromosomal rearrangements, while considering other factors that affect evolutionary rates in single copy DNA. Conclusion: We find that interspecific SD divergence behaves similarly to divergence of single-copy DNA. In contrast, old and recent paralogous copies of SDs do present different patterns of intraspecific divergence. Also, we show that some relatively recent SDs accumulate in regions that carry inversions in sister lineages.

Human high-density lipoprotein particles prevent activation of the JNK pathway induced by human oxidised low-density lipoprotein particles in pancreatic beta cells.

AIMS/HYPOTHESIS: We explored the potential adverse effects of pro-atherogenic oxidised LDL-cholesterol particles on beta cell function. MATERIALS AND METHODS: Isolated human and rat islets and different insulin-secreting cell lines were incubated with human oxidised LDL with or without HDL particles. The insulin level was monitored by ELISA, real-time PCR and a rat insulin promoter construct linked to luciferase gene reporter. Cell apoptosis was determined by scoring cells displaying pycnotic nuclei. RESULTS: Prolonged incubation with human oxidised LDL particles led to a reduction in preproinsulin expression levels, whereas the insulin level was preserved in the presence of native LDL-cholesterol. The loss of insulin production occurred at the transcriptional levels and was associated with an increase in activator protein-1 transcriptional activity. The rise in activator protein-1 activity resulted from activation of c-Jun N-terminal kinases (JNK, now known as mitogen-activated protein kinase 8 [MAPK8]) due to a subsequent decrease in islet-brain 1 (IB1; now known as MAPK8 interacting protein 1) levels. Consistent with the pro-apoptotic role of the JNK pathway, oxidised LDL also induced a twofold increase in the rate of beta cell apoptosis. Treatment of the cells with JNK inhibitor peptides or HDL countered the effects mediated by oxidised LDL. CONCLUSIONS/INTERPRETATION: These data provide strong evidence that oxidised LDL particles exert deleterious effects in the progression of beta cell failure in diabetes and that these effects can be countered by HDL particles.

Impact of phenotype definition on genome-wide association signals: empirical evaluation in human immunodeficiency virus type 1 infection.

Discussion on improving the power of genome-wide association studies to identify candidate variants and genes is generally centered on issues of maximizing sample size; less attention is given to the role of phenotype definition and ascertainment. The authors used genome-wide data from patients infected with human immunodeficiency virus type 1 (HIV-1) to assess whether differences in type of population (622 seroconverters vs. 636 seroprevalent subjects) or the number of measurements available for defining the phenotype resulted in differences in the effect sizes of associations between single nucleotide polymorphisms and the phenotype, HIV-1 viral load at set point. The effect estimate for the top 100 single nucleotide polymorphisms was 0.092 (95% confidence interval: 0.074, 0.110) log(10) viral load (log(10) copies of HIV-1 per mL of blood) greater in seroconverters than in seroprevalent subjects. The difference was even larger when the authors focused on chromosome 6 variants (0.153 log(10) viral load) or on variants that achieved genome-wide significance (0.232 log(10) viral load). The estimates of the genetic effects tended to be slightly larger when more viral load measurements were available, particularly among seroconverters and for variants that achieved genome-wide significance. Differences in phenotype definition and ascertainment may affect the estimated magnitude of genetic effects and should be considered in optimizing power for discovering new associations.

Identification of cancer/testis-antigen genes by massively parallel signature sequencing.

Massively parallel signature sequencing (MPSS) generates millions of short sequence tags corresponding to transcripts from a single RNA preparation. Most MPSS tags can be unambiguously assigned to genes, thereby generating a comprehensive expression profile of the tissue of origin. From the comparison of MPSS data from 32 normal human tissues, we identified 1,056 genes that are predominantly expressed in the testis. Further evaluation by using MPSS tags from cancer cell lines and EST data from a wide variety of tumors identified 202 of these genes as candidates for encoding cancer/testis (CT) antigens. Of these genes, the expression in normal tissues was assessed by RT-PCR in a subset of 166 intron-containing genes, and those with confirmed testis-predominant expression were further evaluated for their expression in 21 cancer cell lines. Thus, 20 CT or CT-like genes were identified, with several exhibiting expression in five or more of the cancer cell lines examined. One of these genes is a member of a CT gene family that we designated as CT45. The CT45 family comprises six highly similar (>98% cDNA identity) genes that are clustered in tandem within a 125-kb region on Xq26.3. CT45 was found to be frequently expressed in both cancer cell lines and lung cancer specimens. Thus, MPSS analysis has resulted in a significant extension of our knowledge of CT antigens, leading to the discovery of a distinctive X-linked CT-antigen gene family.

Computational analysis of the genetic and environmental contributions to disease-related human phenotypes: From predicting adverse lipid response of HIV patients receiving antiretroviral therapy to genome-wide association studies of cardiovascular and lipid related disorders

Genetic variants influence the risk to develop certain diseases or give rise to differences in drug response. Recent progresses in cost-effective, high-throughput genome-wide techniques, such as microarrays measuring Single Nucleotide Polymorphisms (SNPs), have facilitated genotyping of large clinical and population cohorts. Combining the massive genotypic data with measurements of phenotypic traits allows for the determination of genetic differences that explain, at least in part, the phenotypic variations within a population. So far, models combining the most significant variants can only explain a small fraction of the variance, indicating the limitations of current models. In particular, researchers have only begun to address the possibility of interactions between genotypes and the environment. Elucidating the contributions of such interactions is a difficult task because of the large number of genetic as well as possible environmental factors.In this thesis, I worked on several projects within this context. My first and main project was the identification of possible SNP-environment interactions, where the phenotypes were serum lipid levels of patients from the Swiss HIV Cohort Study (SHCS) treated with antiretroviral therapy. Here the genotypes consisted of a limited set of SNPs in candidate genes relevant for lipid transport and metabolism. The environmental variables were the specific combinations of drugs given to each patient over the treatment period. My work explored bioinformatic and statistical approaches to relate patients' lipid responses to these SNPs, drugs and, importantly, their interactions. The goal of this project was to improve our understanding and to explore the possibility of predicting dyslipidemia, a well-known adverse drug reaction of antiretroviral therapy. Specifically, I quantified how much of the variance in lipid profiles could be explained by the host genetic variants, the administered drugs and SNP-drug interactions and assessed the predictive power of these features on lipid responses. Using cross-validation stratified by patients, we could not validate our hypothesis that models that select a subset of SNP-drug interactions in a principled way have better predictive power than the control models using "random" subsets. Nevertheless, all models tested containing SNP and/or drug terms, exhibited significant predictive power (as compared to a random predictor) and explained a sizable proportion of variance, in the patient stratified cross-validation context. Importantly, the model containing stepwise selected SNP terms showed higher capacity to predict triglyceride levels than a model containing randomly selected SNPs. Dyslipidemia is a complex trait for which many factors remain to be discovered, thus missing from the data, and possibly explaining the limitations of our analysis. In particular, the interactions of drugs with SNPs selected from the set of candidate genes likely have small effect sizes which we were unable to detect in a sample of the present size (<800 patients).In the second part of my thesis, I performed genome-wide association studies within the Cohorte Lausannoise (CoLaus). I have been involved in several international projects to identify SNPs that are associated with various traits, such as serum calcium, body mass index, two-hour glucose levels, as well as metabolic syndrome and its components. These phenotypes are all related to major human health issues, such as cardiovascular disease. I applied statistical methods to detect new variants associated with these phenotypes, contributing to the identification of new genetic loci that may lead to new insights into the genetic basis of these traits. This kind of research will lead to a better understanding of the mechanisms underlying these pathologies, a better evaluation of disease risk, the identification of new therapeutic leads and may ultimately lead to the realization of "personalized" medicine.

Antigenicity and immunogenicity of Melan-A/MART-1 derived peptides as targets for tumor reactive CTL in human melanoma.

Some cancer patients mount spontaneous T- and B-cell responses against their tumor cells. Autologous tumor reactive CD8 cytolytic T lymphocyte (CTL) and CD4 T-cell clones as well as antibodies from these patients have been used for the identification of genes encoding the target antigens. This knowledge opened the way for new approaches to the immunotherapy of cancer. In this review, we describe the characterization of the structure-function properties of the melanocyte/melanoma tumor antigen Melan-A/MART-1, the assessment of the T-cell repertoire available against this antigen in healthy individuals, and the analysis of naturally acquired and/or vaccine-induced CTL responses to this antigen in patients with metastatic melanoma.