A conserved NR5A1-responsive enhancer regulates SRY in testis-determination

The Y-linked SRY gene initiates mammalian testis-determination. However, how the expression of SRY is regulated remains elusive. Here, we demonstrate that a conserved steroidogenic factor-1 (SF-1)/NR5A1 binding enhancer is required for appropriate SRY expression to initiate testis-determination in humans. Comparative sequence analysis of SRY 5' regions in mammals identified an evolutionary conserved SF-1/NR5A1-binding motif within a 250 bp region of open chromatin located 5 kilobases upstream of the SRY transcription start site. Genomic analysis of 46,XY individuals with disrupted testis-determination, including a large multigenerational family, identified unique single-base substitutions of highly conserved residues within the SF-1/NR5A1-binding element. In silico modelling and in vitro assays demonstrate the enhancer properties of the NR5A1 motif. Deletion of this hemizygous element by genome-editing, in a novel in vitro cellular model recapitulating human Sertoli cell formation, resulted in a significant reduction in expression of SRY. Therefore, human NR5A1 acts as a regulatory switch between testis and ovary development by upregulating SRY expression, a role that may predate the eutherian radiation. We show that disruption of an enhancer can phenocopy variants in the coding regions of SRY that cause human testis dysgenesis. Since disease causing variants in enhancers are currently rare, the regulation of gene expression in testis-determination offers a paradigm to define enhancer activity in a key developmental process.

Extraocular muscle stem cells exhibit distinct cellular properties associated with non-muscle molecular signatures

Skeletal muscle stem cells (MuSCs) are recognised as functionally heterogeneous. Cranial MuSCs are reported to have greater proliferative and regenerative capacity when compared with those in the limb. A comprehensive understanding of the mechanisms underlying this functional heterogeneity is lacking. Here, we have used clonal analysis, live imaging and single cell transcriptomic analysis to identify crucial features that distinguish extraocular muscle (EOM) from limb muscle stem cell populations. A MyogeninntdTom reporter showed that the increased proliferation capacity of EOM MuSCs correlates with deferred differentiation and lower expression of the myogenic commitment gene Myod. Unexpectedly, EOM MuSCs activated in vitro expressed a large array of extracellular matrix components typical of mesenchymal non-muscle cells. Computational analysis underscored a distinct co-regulatory module, which is absent in limb MuSCs, as driver of these features. The EOM transcription factor network, with Foxc1 as key player, appears to be hardwired to EOM identity as it persists during growth, disease and in vitro after several passages. Our findings shed light on how high-performing MuSCs regulate myogenic commitment by remodelling their local environment and adopting properties not generally associated with myogenic cells.

Multi-omics insights into the biological mechanisms underlying statistical gene-by-lifestyle interactions with smoking and alcohol consumption

Though both genetic and lifestyle factors are known to influence cardiometabolic outcomes, less attention has been given to whether lifestyle exposures can alter the association between a genetic variant and these outcomes. The Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Consortium's Gene-Lifestyle Interactions Working Group has recently published investigations of genome-wide gene-environment interactions in large multi-ancestry meta-analyses with a focus on cigarette smoking and alcohol consumption as lifestyle factors and blood pressure and serum lipids as outcomes. Further description of the biological mechanisms underlying these statistical interactions would represent a significant advance in our understanding of gene-environment interactions, yet accessing and harmonizing individual-level genetic and 'omics data is challenging. Here, we demonstrate the coordinated use of summary-level data for gene-lifestyle interaction associations on up to 600,000 individuals, differential methylation data, and gene expression data for the characterization and prioritization of loci for future follow-up analyses. Using this approach, we identify 48 genes for which there are multiple sources of functional support for the identified gene-lifestyle interaction. We also identified five genes for which differential expression was observed by the same lifestyle factor for which a gene-lifestyle interaction was found. For instance, in gene-lifestyle interaction analysis, the T allele of rs6490056 (ALDH2) was associated with higher systolic blood pressure, and a larger effect was observed in smokers compared to non-smokers. In gene expression studies, this allele is associated with decreased expression of ALDH2, which is part of a major oxidative pathway. Other results show increased expression of ALDH2 among smokers. Oxidative stress is known to contribute to worsening blood pressure. Together these data support the hypothesis that rs6490056 reduces expression of ALDH2, which raises oxidative stress, leading to an increase in blood pressure, with a stronger effect among smokers, in whom the burden of oxidative stress is greater. Other genes for which the aggregation of data types suggest a potential mechanism include: GCNT4×current smoking (HDL), PTPRZ1×ever-smoking (HDL), SYN2×current smoking (pulse pressure), and TMEM116×ever-smoking (mean arterial pressure). This work demonstrates the utility of careful curation of summary-level data from a variety of sources to prioritize gene-lifestyle interaction loci for follow-up analyses.

Identification of New Biological Pathways Involved in Skin Aging From the Analysis of French Women Genome-Wide Data

Skin aging is an ineluctable process leading to the progressive loss of tissue integrity and is characterized by various outcomes such as wrinkling and sagging. Researchers have identified impacting environmental factors (sun exposure, smoking, etc.) and several molecular mechanisms leading to skin aging. We have previously performed genome-wide association studies (GWAS) in 502 very-well characterized French women, looking for associations with four major outcomes of skin aging, namely, photoaging, solar lentigines, wrinkling, and sagging, and this has led to new insights into the molecular mechanisms of skin aging. Since individual SNP associations in GWAS explain only a small fraction of the genetic impact in complex polygenic phenotypes, we have made the integration of these genotypes into the reference Kegg biological pathways and looked for associations by the gene set enrichment analysis (GSEA) approach. 106 pathways were tested for association with the four outcomes of skin aging. This biological pathway analysis revealed new relevant pathways and genes, some likely specific of skin aging such as the WNT7B and PRKCA genes in the “melanogenesis” pathway and some likely involved in global aging such as the DDB1 gene in the “nucleotide excision repair” pathway, not picked up in the previously published GWAS. Overall, our results suggest that the four outcomes of skin aging possess specific molecular mechanisms such as the “proteasome” and “mTOR signaling pathway” but may also share common molecular mechanisms such as “nucleotide excision repair.”

Multitrait GWAS to connect disease variants and biological mechanisms

Genome-wide association studies (GWASs) have uncovered a wealth of associations between common variants and human phenotypes. Here, we present an integrative analysis of GWAS summary statistics from 36 phenotypes to decipher multitrait genetic architecture and its link with biological mechanisms. Our framework incorporates multitrait association mapping along with an investigation of the breakdown of genetic associations into clusters of variants harboring similar multitrait association profiles. Focusing on two subsets of immunity and metabolism phenotypes, we then demonstrate how genetic variants within clusters can be mapped to biological pathways and disease mechanisms. Finally, for the metabolism set, we investigate the link between gene cluster assignment and the success of drug targets in randomized controlled trials.

Deriving stratified effects from joint models investigating gene-environment interactions

Background

Models including an interaction term and performing a joint test of SNP and/or interaction effect are often used to discover Gene-Environment (GxE) interactions. When the environmental exposure is a binary variable, analyses from exposure-stratified models which consist of estimating genetic effect in unexposed and exposed individuals separately can be of interest. In large-scale consortia focusing on GxE interactions in which only the joint test has been performed, it may be challenging to get summary statistics from both exposure-stratified and marginal (i.e not accounting for interaction) models.

Results

In this work, we developed a simple framework to estimate summary statistics in each stratum of a binary exposure and in the marginal model using summary statistics from the “joint” model. We performed simulation studies to assess our estimators’ accuracy and examined potential sources of bias, such as correlation between genotype and exposure and differing phenotypic variances within exposure strata. Results from these simulations highlight the high theoretical accuracy of our estimators and yield insights into the impact of potential sources of bias. We then applied our methods to real data and demonstrate our estimators’ retained accuracy after filtering SNPs by sample size to mitigate potential bias.

Conclusions

These analyses demonstrated the accuracy of our method in estimating both stratified and marginal summary statistics from a joint model of gene-environment interaction. In addition to facilitating the interpretation of GxE screenings, this work could be used to guide further functional analyses. We provide a user-friendly Python script to apply this strategy to real datasets. The Python script and documentation are available at https://gitlab.pasteur.fr/statistical-genetics/j2s.

JASS: command line and web interface for the joint analysis of GWAS results

Genome-wide association study (GWAS) has been the driving force for identifying association between genetic variants and human phenotypes. Thousands of GWAS summary statistics covering a broad range of human traits and diseases are now publicly available. These GWAS have proven their utility for a range of secondary analyses, including in particular the joint analysis of multiple phenotypes to identify new associated genetic variants. However, although several methods have been proposed, there are very few large-scale applications published so far because of challenges in implementing these methods on real data. Here, we present JASS (Joint Analysis of Summary Statistics), a polyvalent Python package that addresses this need. Our package incorporates recently developed joint tests such as the omnibus approach and various weighted sum of Z-score tests while solving all practical and computational barriers for large-scale multivariate analysis of GWAS summary statistics. This includes data cleaning and harmonization tools, an efficient algorithm for fast derivation of joint statistics, an optimized data management process and a web interface for exploration purposes. Both benchmark analyses and real data applications demonstrated the robustness and strong potential of JASS for the detection of new associated genetic variants. Our package is freely available at https://gitlab.pasteur.fr/statistical-genetics/jass.

Multi-ancestry sleep-by-SNP interaction analysis in 126,926 individuals reveals lipid loci stratified by sleep duration

Both short and long sleep are associated with an adverse lipid profile, likely through different biological pathways. To elucidate the biology of sleep-associated adverse lipid profile, we conduct multi-ancestry genome-wide sleep-SNP interaction analyses on three lipid traits (HDL-c, LDL-c and triglycerides). In the total study sample (discovery + replication) of 126,926 individuals from 5 different ancestry groups, when considering either long or short total sleep time interactions in joint analyses, we identify 49 previously unreported lipid loci, and 10 additional previously unreported lipid loci in a restricted sample of European-ancestry cohorts. In addition, we identify new gene-sleep interactions for known lipid loci such as LPL and PCSK9. The previously unreported lipid loci have a modest explained variance in lipid levels: most notable, gene-short-sleep interactions explain 4.25% of the variance in triglyceride level. Collectively, these findings contribute to our understanding of the biological mechanisms involved in sleep-associated adverse lipid profiles.

Sleep duration is associated with an adverse lipid profile. Here, the authors perform genome-wide gene-by-sleep interaction analysis and find 49 previously unreported lipid loci when considering short or long total sleep time.

VarExp: estimating variance explained by genome-wide GxE summary statistics

Summary

Many genome-wide association studies and genome-wide screening for gene-environment (GxE) interactions have been performed to elucidate the underlying mechanisms of human traits and diseases. When the analyzed outcome is quantitative, the overall contribution of identified genetic variants to the outcome is often expressed as the percentage of phenotypic variance explained. This is commonly done using individual-level genotype data but it is challenging when results are derived through meta-analyses. Here, we present R package, 'VarExp', that allows for the estimation of the percentage of phenotypic variance explained using summary statistics only. It allows for a range of models to be evaluated, including marginal genetic effects, GxE interaction effects and both effects jointly. Its implementation integrates all recent methodological developments and does not need external data to be uploaded by users.

Availability and implementation

The R package is available at https://gitlab.pasteur.fr/statistical-genetics/VarExp.git.

Supplementary information

Supplementary data are available at Bioinformatics online.

Genetic Correlations Between Diabetes and Glaucoma: An Analysis of Continuous and Dichotomous Phenotypes

PURPOSE

A genetic correlation is the proportion of phenotypic variance between traits that is shared on a genetic basis. Here we explore genetic correlations between diabetes- and glaucoma-related traits.

DESIGN

Cross-sectional study.

METHODS

We assembled genome-wide association study summary statistics from European-derived participants regarding diabetes-related traits like fasting blood sugar (FBS) and type 2 diabetes (T2D) and glaucoma-related traits (intraocular pressure [IOP], central corneal thickness [CCT], corneal hysteresis [CH], corneal resistance factor [CRF], cup-to-disc ratio [CDR], and primary open-angle glaucoma [POAG]). We included data from the National Eye Institute Glaucoma Human Genetics Collaboration Heritable Overall Operational Database, the UK Biobank, and the International Glaucoma Genetics Consortium. We calculated genetic correlation (r_g) between traits using linkage disequilibrium score regression. We also calculated genetic correlations between IOP, CCT, and select diabetes-related traits based on individual level phenotype data in 2 Northern European population-based samples using pedigree information and Sequential Oligogenic Linkage Analysis Routines.

RESULTS

Overall, there was little r_g between diabetes- and glaucoma-related traits. Specifically, we found a nonsignificant negative correlation between T2D and POAG (r_g = -0.14; P = .16). Using Sequential Oligogenic Linkage Analysis Routines, the genetic correlations between measured IOP, CCT, FBS, fasting insulin, and hemoglobin A1c were null. In contrast, genetic correlations between IOP and POAG (r_g ≥ 0.45; P ≤ 3.0 × 10^-4) and between CDR and POAG were high (r_g = 0.57; P = 2.8 × 10^-10). However, genetic correlations between corneal properties (CCT, CRF, and CH) and POAG were low (r_g range -0.18 to 0.11) and nonsignificant (P ≥ .07).

CONCLUSION

These analyses suggest that there is limited genetic correlation between diabetes- and glaucoma-related traits.

Genetic effects on the commensal microbiota in inflammatory bowel disease patients

Several bacteria in the gut microbiota have been shown to be associated with inflammatory bowel disease (IBD), and dozens of IBD genetic variants have been identified in genome-wide association studies. However, the role of the microbiota in the etiology of IBD in terms of host genetic susceptibility remains unclear. Here, we studied the association between four major genetic variants associated with an increased risk of IBD and bacterial taxa in up to 633 IBD cases. We performed systematic screening for associations, identifying and replicating associations between NOD2 variants and two taxa: the Roseburia genus and the Faecalibacterium prausnitzii species. By exploring the overall association patterns between genes and bacteria, we found that IBD risk alleles were significantly enriched for associations concordant with bacteria-IBD associations. To understand the significance of this pattern in terms of the study design and known effects from the literature, we used counterfactual principles to assess the fitness of a few parsimonious gene-bacteria-IBD causal models. Our analyses showed evidence that the disease risk of these genetic variants were likely to be partially mediated by the microbiome. We confirmed these results in extensive simulation studies and sensitivity analyses using the association between NOD2 and F. prausnitzii as a case study.

Joint Analysis of Multiple Interaction Parameters in Genetic Association Studies

With growing human genetic and epidemiologic data, there has been increased interest for the study of gene-by-environment (G-E) interaction effects. Still, major questions remain on how to test jointly a large number of interactions between multiple SNPs and multiple exposures. In this study, we first compared the relative performance of four fixed-effect joint analysis approaches using simulated data, considering up to 10 exposures and 300 SNPs: (1) omnibus test, (2) multi-exposure and genetic risk score (GRS) test, (3) multi-SNP and environmental risk score (ERS) test, and (4) GRS-ERS test. Our simulations explored both linear and logistic regression while considering three statistics: the Wald test, the Score test, and the likelihood ratio test (LRT). We further applied the approaches to three large sets of human cohort data (n = 37,664), focusing on type 2 diabetes (T2D), obesity, hypertension, and coronary heart disease with smoking, physical activity, diets, and total energy intake. Overall, GRS-based approaches were the most robust, and had the highest power, especially when the G-E interaction effects were correlated with the marginal genetic and environmental effects. We also observed severe miscalibration of joint statistics in logistic models when the number of events per variable was too low when using either the Wald test or LRT test. Finally, our real data application detected nominally significant interaction effects for three outcomes (T2D, obesity, and hypertension), mainly from the GRS-ERS approach. In conclusion, this study provides guidelines for testing multiple interaction parameters in modern human cohorts including extensive genetic and environmental data.

A test for gene-environment interaction in the presence of measurement error in the environmental variable

The identification of gene-environment interactions in relation to risk of human diseases has been challenging. One difficulty has been that measurement error in the exposure can lead to massive reductions in the power of the test, as well as in bias toward the null in the interaction effect estimates. Leveraging previous work on linear discriminant analysis, we develop a new test of interaction between genetic variants and a continuous exposure that mitigates these detrimental impacts of exposure measurement error in ExG testing by reversing the role of exposure and the diseases status in the fitted model, thus transforming the analysis to standard linear regression. Through simulation studies, we show that the proposed approach is valid in the presence of classical exposure measurement error as well as when there is correlation between the exposure and the genetic variant. Simulations also demonstrated that the reverse test has greater power compared to logistic regression. Finally, we confirmed that our approach eliminates bias from exposure measurement error in estimation. Computing times are reduced by as much as fivefold in this new approach. For illustrative purposes, we applied the new approach to an ExGWAS study of interactions with alcohol and body mass index among 1,145 cases with invasive breast cancer and 1,142 controls from the Cancer Genetic Markers of Susceptibility study.

Discriminating Agonist from Antagonist Ligands of the Nuclear Receptors Using Different Chemoinformatics Approaches

Nuclear receptors (NRs) constitute an important class of therapeutic targets. During the last 4 years, we tackled the pharmacological profile assessment of NR ligands for which we constructed the NRLiSt BDB. We evaluated and compared the performance of different virtual screening approaches: mean of molecular descriptor distribution values, molecular docking and 3D pharmacophore models. The simple comparison of the distribution profiles of 4885 molecular descriptors between the agonist and antagonist datasets didn't provide satisfying results. We obtained an overall good performance with the docking method we used, Surflex-Dock which was able to discriminate agonist from antagonist ligands. But the availability of PDB structures in the "pharmacological-profile-to-predict-bound-state" (agonist-bound or antagonist-bound) and the availability of enough ligands of both pharmacological profiles constituted limits to generalize this protocol for all NRs. Finally, the 3D pharmacophore modeling approach, allowed us to generate selective agonist pharmacophores and selective antagonist pharmacophores that covered more than 99 % of the whole NRLiSt BDB. This study allowed a better understanding of the pharmacological modulation of NRs with small molecules and could be extended to other therapeutic classes.

A new 3p25 locus is associated with liver fibrosis progression in human immunodeficiency virus/hepatitis C virus-coinfected patients

There is growing evidence that human genetic variants contribute to liver fibrosis in subjects with hepatitis C virus (HCV) monoinfection, but this aspect has been little investigated in patients coinfected with HCV and human immunodeficiency virus (HIV). We performed the first genome-wide association study of liver fibrosis progression in patients coinfected with HCV and HIV, using the well-characterized French National Agency for Research on AIDS and Viral Hepatitis CO13 HEPAVIH cohort. Liver fibrosis was assessed by elastography (FibroScan), providing a quantitative fibrosis score. After quality control, a genome-wide association study was conducted on 289 Caucasian patients, for a total of 8,426,597 genotyped (Illumina Omni2.5 BeadChip) or reliably imputed single-nucleotide polymorphisms. Single-nucleotide polymorphisms with P values <10^-6 were investigated in two independent replication cohorts of European patients infected with HCV alone. Two signals of genome-wide significance (P < 5 × 10^-8 ) were obtained. The first, on chromosome 3p25 and corresponding to rs61183828 (P = 3.8 × 10^-9 ), was replicated in the two independent cohorts of patients with HCV monoinfection. The cluster of single-nucleotide polymorphisms in linkage disequilibrium with rs61183828 was located close to two genes involved in mechanisms affecting both cell signaling and cell structure (CAV3) or HCV replication (RAD18). The second signal, obtained with rs11790131 (P = 9.3 × 10^-9 ) on chromosome region 9p22, was not replicated.

CONCLUSION

This genome-wide association study identified a new locus associated with liver fibrosis severity in patients with HIV/HCV coinfection, on chromosome 3p25, a finding that was replicated in patients with HCV monoinfection; these results provide new relevant hypotheses for the pathogenesis of liver fibrosis in patients with HIV/HCV coinfection that may help define new targets for drug development or new prognostic tests, to improve patient care. (Hepatology 2016;64:1462-1472).

Identification of Genes Whose Expression Profile Is Associated with Non-Progression towards AIDS Using eQTLs

BACKGROUND

Many genome-wide association studies have been performed on progression towards the acquired immune deficiency syndrome (AIDS) and they mainly identified associations within the HLA loci. In this study, we demonstrate that the integration of biological information, namely gene expression data, can enhance the sensitivity of genetic studies to unravel new genetic associations relevant to AIDS.

METHODS

We collated the biological information compiled from three databases of expression quantitative trait loci (eQTLs) involved in cells of the immune system. We derived a list of single nucleotide polymorphisms (SNPs) that are functional in that they correlate with differential expression of genes in at least two of the databases. We tested the association of those SNPs with AIDS progression in two cohorts, GRIV and ACS. Tests on permuted phenotypes of the GRIV and ACS cohorts or on randomised sets of equivalent SNPs allowed us to assess the statistical robustness of this method and to estimate the true positive rate.

RESULTS

Eight genes were identified with high confidence (p = 0.001, rate of true positives 75%). Some of those genes had previously been linked with HIV infection. Notably, ENTPD4 belongs to the same family as CD39, whose expression has already been associated with AIDS progression; while DNAJB12 is part of the HSP90 pathway, which is involved in the control of HIV latency. Our study also drew our attention to lesser-known functions such as mitochondrial ribosomal proteins and a zinc finger protein, ZFP57, which could be central to the effectiveness of HIV infection. Interestingly, for six out of those eight genes, down-regulation is associated with non-progression, which makes them appealing targets to develop drugs against HIV.

Evidence after imputation for a role of MICA variants in nonprogression and elite control of HIV type 1 infection

Past genome-wide association studies (GWAS) involving individuals with AIDS have mainly identified associations in the HLA region. Using the latest software, we imputed 7 million single-nucleotide polymorphisms (SNPs)/indels of the 1000 Genomes Project from the GWAS-determined genotypes of individuals in the Genomics of Resistance to Immunodeficiency Virus AIDS nonprogression cohort and compared them with those of control cohorts. The strongest signals were in MICA, the gene encoding major histocompatibility class I polypeptide-related sequence A (P = 3.31 × 10(-12)), with a particular exonic deletion (P = 1.59 × 10(-8)) in full linkage disequilibrium with the reference HCP5 rs2395029 SNP. Haplotype analysis also revealed an additive effect between HLA-C, HLA-B, and MICA variants. These data suggest a role for MICA in progression and elite control of human immunodeficiency virus type 1 infection.

NRLiSt BDB, the manually curated nuclear receptors ligands and structures benchmarking database

Nuclear receptors (NRs) constitute an important class of drug targets. We created the most exhaustive NR-focused benchmarking database to date, the NRLiSt BDB (NRs ligands and structures benchmarking database). The 9905 compounds and 339 structures of the NRLiSt BDB are ready for structure-based and ligand-based virtual screening. In the present study, we detail the protocol used to generate the NRLiSt BDB and its features. We also give some examples of the errors that we found in ChEMBL that convinced us to manually review all original papers. Since extensive and manually curated experimental data about NR ligands and structures are provided in the NRLiSt BDB, it should become a powerful tool to assess the performance of virtual screening methods on NRs, to assist the understanding of NR's function and modulation, and to support the discovery of new drugs targeting NRs. NRLiSt BDB is freely available online at http://nrlist.drugdesign.fr .

Genome-wide expression profiling and tissue array analysis for prediction of recurrences in stage II-III colorectal cancer (CRC)

3572

Background: Despite substantial progress in molecular pathogenesis of colon cancer (CC), no reliable biomarkers of outcome have yet been identified in patients with resected stage II-III CC. Methods: We analyzed genome-wide mRNA expression profiles in 20 stage II or III left side CC from 10 patient who developed metastasis (M+) and 10 disease free patients followed up for at least 4 years (M-) using high-density oligonucleotide microarrays (Agilent technology). RNA from tumor tissue (T) was hybridized against normal tissue (NT) from the same patient and each experiment was replicated 4 times (with 2 dye-swaps). The goal was to select genes both differentially expressed between T and NT and between M+ and M-. A tissu-array was constructed using 212 stage II and III resected CRC (164 CC, 64 rectal cancers) and their matched NT. For survival analysis, immunohistochemistry (IHC) data was dichotomized at the median value. Results: Analysis of microarray data yielded 27 genes that had a 2-fold difference between the expression in T and NT in at least 5 out of 20 patients and for which the average expression was significantly different between M+ and M- (p<0.01, t-test). Among the 6 most differentially expressed genes between M+ and M- in T, 4 of them were found to be involved in interferon γ pathway and could be evaluated by IHC: CXCL9, CXCL13, PPARγ, THSD. In order to assess macrophage and natural killer (NK) cell infiltration, CD68 and CD57 were also analyzed. Intensity was measured by semi-quantitative scores for the first 4 genes and by the number of infiltrating cells for the others. CXCL9, PPARGγ, CXCL13, CD57 and CD68 were significantly underexpressed in T as compared to NT (p<0.0001, paired t-test). The logrank test stratified by cancer site indicated that high IHC expression of CD57 possessed a significantly better recurrence-free survival (RFS) than those low expression (p=0.004). Multivariate Cox analysis identified tumor site (p=0.001), node stage (p<0.001) and CD57 (p=0.002) as independent predictors of RFS. Conclusions: NK cell infiltration within colorectal cancers is associated with prolonged recurrence-free survival in stage II and III CRC.

No significant financial relationships to disclose.

Gene expression profiles of bladder cancers: evidence for a striking effect of in vitro cell models on gene patterns

In order to assess the effect of in vitro models on the expression of key genes known to be implicated in the development or progression of cancer, we quantified by real-time quantitative PCR the expression of 28 key genes in three bladder cancer tissue specimens and in their derived cell lines, studied either as one-dimensional single cell suspensions, two-dimensional monolayers or three-dimensional spheroids. Global analysis of gene expression profiles showed that in vitro models had a dramatic impact upon gene expression. Remarkably, quantitative differences in gene expression of 2-63-fold were observed in 24 out of 28 genes among the cell models. In addition, we observed that the in vitro model which most closely mimicked in vivo mRNA phenotype varied with both the gene and the patient. These results provide evidence that mRNA expression databases based on cancer cell lines, which are studied to provide a rationale for selection of therapy on the basis of molecular characteristics of a patient's tumour, must be carefully interpreted.