Gene expression profiling in breast cancer: understanding the molecular basis of histologic grade to improve prognosis.

BACKGROUND: Histologic grade in breast cancer provides clinically important prognostic information. However, 30%-60% of tumors are classified as histologic grade 2. This grade is associated with an intermediate risk of recurrence and is thus not informative for clinical decision making. We examined whether histologic grade was associated with gene expression profiles of breast cancers and whether such profiles could be used to improve histologic grading. METHODS: We analyzed microarray data from 189 invasive breast carcinomas and from three published gene expression datasets from breast carcinomas. We identified differentially expressed genes in a training set of 64 estrogen receptor (ER)-positive tumor samples by comparing expression profiles between histologic grade 3 tumors and histologic grade 1 tumors and used the expression of these genes to define the gene expression grade index. Data from 597 independent tumors were used to evaluate the association between relapse-free survival and the gene expression grade index in a Kaplan-Meier analysis. All statistical tests were two-sided. RESULTS: We identified 97 genes in our training set that were associated with histologic grade; most of these genes were involved in cell cycle regulation and proliferation. In validation datasets, the gene expression grade index was strongly associated with histologic grade 1 and 3 status; however, among histologic grade 2 tumors, the index spanned the values for histologic grade 1-3 tumors. Among patients with histologic grade 2 tumors, a high gene expression grade index was associated with a higher risk of recurrence than a low gene expression grade index (hazard ratio = 3.61, 95% confidence interval = 2.25 to 5.78; P < .001, log-rank test). CONCLUSIONS: Gene expression grade index appeared to reclassify patients with histologic grade 2 tumors into two groups with high versus low risks of recurrence. This approach may improve the accuracy of tumor grading and thus its prognostic value.

The future is genome-wide.

A report of the annual meeting of the European Society of Human Genetics, Amsterdam, 6-9 May 2006.

Blood pressure loci identified with a gene-centric array.

Raised blood pressure (BP) is a major risk factor for cardiovascular disease. Previous studies have identified 47 distinct genetic variants robustly associated with BP, but collectively these explain only a few percent of the heritability for BP phenotypes. To find additional BP loci, we used a bespoke gene-centric array to genotype an independent discovery sample of 25,118 individuals that combined hypertensive case-control and general population samples. We followed up four SNPs associated with BP at our p < 8.56 × 10(-7) study-specific significance threshold and six suggestively associated SNPs in a further 59,349 individuals. We identified and replicated a SNP at LSP1/TNNT3, a SNP at MTHFR-NPPB independent (r(2) = 0.33) of previous reports, and replicated SNPs at AGT and ATP2B1 reported previously. An analysis of combined discovery and follow-up data identified SNPs significantly associated with BP at p < 8.56 × 10(-7) at four further loci (NPR3, HFE, NOS3, and SOX6). The high number of discoveries made with modest genotyping effort can be attributed to using a large-scale yet targeted genotyping array and to the development of a weighting scheme that maximized power when meta-analyzing results from samples ascertained with extreme phenotypes, in combination with results from nonascertained or population samples. Chromatin immunoprecipitation and transcript expression data highlight potential gene regulatory mechanisms at the MTHFR and NOS3 loci. These results provide candidates for further study to help dissect mechanisms affecting BP and highlight the utility of studying SNPs and samples that are independent of those studied previously even when the sample size is smaller than that in previous studies.

Mutation prevalence of cerebral cavernous malformation genes in Spanish patients.

OBJECTIVE To study the molecular genetic and clinical features of cerebral cavernous malformations (CCM) in a cohort of Spanish patients. METHODS We analyzed the CCM1, CCM2, and CCM3 genes by MLPA and direct sequencing of exons and intronic boundaries in 94 familial forms and 41 sporadic cases of CCM patients of Spanish extraction. When available, RNA studies were performed seeking for alternative or cryptic splicing. RESULTS A total of 26 pathogenic mutations, 22 of which predict truncated proteins, were identified in 29 familial forms and in three sporadic cases. The repertoire includes six novel non-sense and frameshift mutations in CCM1 and CCM3. We also found four missense mutations, one of them located at the third NPXY motif of CCM1 and another one that leads to cryptic splicing of CCM1 exon 6. We found four genomic deletions with the loss of the whole CCM2 gene in one patient and a partial loss of CCM1and CCM2 genes in three other patients. Four families had mutations in CCM3. The results include a high frequency of intronic variants, although most of them localize out of consensus splicing sequences. The main symptoms associated to clinical debut consisted of cerebral haemorrhage, migraines and epileptic seizures. The rare co-occurrence of CCM with Noonan and Chiari syndromes and delayed menarche is reported. CONCLUSIONS Analysis of CCM genes by sequencing and MLPA has detected mutations in almost 35% of a Spanish cohort (36% of familial cases and 10% of sporadic patients). The results include 13 new mutations of CCM genes and the main clinical symptoms that deserves consideration in molecular diagnosis and genetic counselling of cerebral cavernous malformations.

Expression profiling of rectal tumors defines response to neoadjuvant treatment related genes.

To date, no effective method exists that predicts the response to preoperative chemoradiation (CRT) in locally advanced rectal cancer (LARC). Nevertheless, identification of patients who have a higher likelihood of responding to preoperative CRT could be crucial in decreasing treatment morbidity and avoiding expensive and time-consuming treatments. The aim of this study was to identify signatures or molecular markers related to response to pre-operative CRT in LARC. We analyzed the gene expression profiles of 26 pre-treatment biopsies of LARC (10 responders and 16 non-responders) without metastasis using Human WG CodeLink microarray platform. Two hundred and fifty seven genes were differentially over-expressed in the responder patient subgroup. Ingenuity Pathway Analysis revealed a significant ratio of differentially expressed genes related to cancer, cellular growth and proliferation pathways, and c-Myc network. We demonstrated that high Gng4, c-Myc, Pola1, and Rrm1 mRNA expression levels was a significant prognostic factor for response to treatment in LARC patients (p<0.05). Using this gene set, we were able to establish a new model for predicting the response to CRT in rectal cancer with a sensitivity of 60% and 100% specificity. Our results reflect the value of gene expression profiling to gain insight about the molecular pathways involved in the response to treatment of LARC patients. These findings could be clinically relevant and support the use of mRNA levels when aiming to identify patients who respond to CRT therapy.

The MicroArray Quality Control (MAQC)-II study of common practices for the development and validation of microarray-based predictive models.

Gene expression data from microarrays are being applied to predict preclinical and clinical endpoints, but the reliability of these predictions has not been established. In the MAQC-II project, 36 independent teams analyzed six microarray data sets to generate predictive models for classifying a sample with respect to one of 13 endpoints indicative of lung or liver toxicity in rodents, or of breast cancer, multiple myeloma or neuroblastoma in humans. In total, >30,000 models were built using many combinations of analytical methods. The teams generated predictive models without knowing the biological meaning of some of the endpoints and, to mimic clinical reality, tested the models on data that had not been used for training. We found that model performance depended largely on the endpoint and team proficiency and that different approaches generated models of similar performance. The conclusions and recommendations from MAQC-II should be useful for regulatory agencies, study committees and independent investigators that evaluate methods for global gene expression analysis.

Characterization of microsatellite loci and reliable genotyping in a polyploid plant, Mercurialis perennis (Euphorbiaceae).

For many applications in population genetics, codominant simple sequence repeats (SSRs) may have substantial advantages over dominant anonymous markers such as amplified fragment length polymorphisms (AFLPs). In high polyploids, however, allele dosage of SSRs cannot easily be determined and alleles are not easily attributable to potentially diploidized loci. Here, we argue that SSRs may nonetheless be better than AFLPs for polyploid taxa if they are analyzed as effectively dominant markers because they are more reliable and more precise. We describe the transfer of SSRs developed for diploid Mercurialis huetii to the clonal dioecious M. perennis. Primers were tested on a set of 54 male and female plants from natural decaploid populations. Eight of 65 tested loci produced polymorphic fragments. Binary profiles from 4 different scoring routines were used to define multilocus lineages (MLLs). Allowing for fragment differences within 1 MLL, all analyses revealed the same 14 MLLs without conflicting with merigenet, sex, or plot assignment. For semiautomatic scoring, a combination of as few as 2 of the 4 most polymorphic loci resulted in unambiguous discrimination of clones. Our study demonstrates that microsatellite fingerprinting of polyploid plants is a cost efficient and reliable alternative to AFLPs, not least because fewer loci are required than for diploids.

A positive FGFR3/FOXN1 feedback loop underlies benign skin keratosis versus squamous cell carcinoma formation in humans.

Seborrheic keratoses (SKs) are common, benign epithelial tumors of the skin that do not, or very rarely, progress into malignancy, for reasons that are not understood. We investigated this by gene expression profiling of human SKs and cutaneous squamous cell carcinomas (SCCs) and found that several genes previously connected with keratinocyte tumor development were similarly modulated in SKs and SCCs, whereas the expression of others differed by only a few fold. In contrast, the tyrosine kinase receptor FGF receptor-3 (FGFR3) and the transcription factor forkhead box N1 (FOXN1) were highly expressed in SKs, and close to undetectable in SCCs. We also showed that increased FGFR3 activity was sufficient to induce FOXN1 expression, counteract the inhibitory effect of EGFR signaling on FOXN1 expression and differentiation, and induce differentiation in a FOXN1-dependent manner. Knockdown of FOXN1 expression in primary human keratinocytes cooperated with oncogenic RAS in the induction of SCC-like tumors, whereas increased FOXN1 expression triggered the SCC cells to shift to a benign SK-like tumor phenotype, which included increased FGFR3 expression. Thus,we have uncovered a positive regulatory loop between FGFR3 and FOXN1 that underlies a benign versus malignant skin tumor phenotype.

Parallel changes in genetic diversity and species diversity following a natural disturbance.

We examined the spatial and temporal variation of species diversity and genetic diversity in a metacommunity comprising 16 species of freshwater gastropods. We monitored species abundance at five localities of the Ain river floodplain in southeastern France, over a period of four years. Using 190 AFLP loci, we monitored the genetic diversity of Radix balthica, one of the most abundant gastropod species of the metacommunity, twice during that period. An exceptionally intense drought occurred during the last two years and differentially affected the study sites. This allowed us to test the effect of natural disturbances on changes in both genetic and species diversity. Overall, local (alpha) diversity declined as reflected by lower values of gene diversity H(S) and evenness. In parallel, the among-sites (beta) diversity increased at both the genetic (F(ST)) and species (F(STC)) levels. These results suggest that disturbances can lead to similar changes in genetic and community structure through the combined effects of selective and neutral processes.

The Life Course Determinants of Vulnerability in Late Careers

Late career is often seen as a more vulnerable life-stage in the labour market, in which workers may experience a deterioration in job quality. Using a life course perspective and longitudinal data, this article analyses the vulnerability associated with late career by focusing on four occupational dimensions: working-time, career continuity, retirement timing and income change. The research is carried out using data from Switzerland, a country where the age profile of the labour force is an increasing issue. The paper also adopts a cumulative disadvantage perspective to examine the impact of previous work and family life experiences on work life vulnerability at older age. Our data come from the Survey of Health, Ageing and Retirement in Europe (SHARELIFE). The paper uses cluster analysis, sequence analysis and ordered logistic regression. Results show that women with previous family responsibilities resulting in long-term unemployment or caring, often with health complications, are more likely to be vulnerable to deterioration in job quality in late career. This suggests that experiences in the last period of the working life may be just as gendered as earlier periods.

Describing ancient horizontal gene transfers at the nucleotide and gene levels by comparative pathogenicity island genometrics.

MOTIVATION: Lateral gene transfer is a major mechanism contributing to bacterial genome dynamics and pathovar emergence via pathogenicity island (PAI) spreading. However, since few of these genomic exchanges are experimentally reproducible, it is difficult to establish evolutionary scenarios for the successive PAI transmissions between bacterial genera. Methods initially developed at the gene and/or nucleotide level for genomics, i.e. comparisons of concatenated sequences, ortholog frequency, gene order or dinucleotide usage, were combined and applied here to homologous PAIs: we call this approach comparative PAI genometrics. RESULTS: YAPI, a Yersinia PAI, and related islands were compared with measure evolutionary relationships between related modules. Through use of our genometric approach designed for tracking codon usage adaptation and gene phylogeny, an ancient inter-genus PAI transfer was oriented for the first time by characterizing the genomic environment in which the ancestral island emerged and its subsequent transfers to other bacterial genera.

Affymetrix Whole-Transcript Human Gene 1.0 ST array is highly concordant with standard 3' expression arrays.

The recently released Affymetrix Human Gene 1.0 ST array has two major differences compared with standard 3' based arrays: (i) it interrogates the entire mRNA transcript, and (ii) it uses DNA targets. To assess the impact of these differences on array performance, we performed a series of comparative hybridizations between the Human Gene 1.0 ST and the Affymetrix HG-U133 Plus 2.0 and the Illumina HumanRef-8 BeadChip arrays. Additionally, both RNA and DNA targets were hybridized on HG-U133 Plus 2.0 arrays. The results show that the overall reproducibility of the Gene 1.0 ST array is best. When looking only at the high intensity probes, the reproducibility of the Gene 1.0 ST array and the Illumina BeadChip array is equally good. Concordance of array results was assessed using different inter-platform mappings. Agreements are best between the two labeling protocols using HG-U133 Plus 2.0 array. The Gene 1.0 ST array is most concordant with the HG-U133 array hybridized with cDNA targets. This may reflect the impact of the target type. Overall, the high degree of correspondence provides strong evidence for the reliability of the Gene 1.0 ST array.

The polymorphic nature of HIV type 1 env V4 affects the patterns of potential N-glycosylation sites in proviral DNA at the intrahost level.

We have previously shown that env V4 from HIV-1 plasma RNA is highly heterogeneous within a single patient, due to indel-associated polymorphism. In this study, we have analyzed the variability of V4 in proviral DNA from unfractionated PBMC and sorted T and non-T cell populations within individual patients. Our data show that the degree of sequence variability and length polymorphism in V4 from HIV provirus is even higher than we previously reported in plasma. The data also show that the sequence of V4 depends largely on the experimental approach chosen. We could observe no clear trend for compartmentalization of V4 variants in specific cell types. Of interest is the fact that some variants that had been found to be predominant in plasma were not detected in any of the cell subsets analyzed. Consistently with our observations in plasma, V3 was found to be relatively conserved at both interpatient and intrapatient level. Our data show that V4 polymorphism involving insertions and deletions in addition to point mutations results in changes in the patterns of sequons in HIV-1 proviral DNA as well as in plasma RNA. These rearrangements may result in the coexistence, within the same individual, of a swarm of different V4 regions, each characterized by a different carbohydrate surface shield. Further studies are needed to investigate the mechanism responsible for the variability observed in V4 and its role in HIV pathogenesis.

Structure-function analyses point to a polynucleotide-accommodating groove essential for APOBEC3A restriction activities.

Members of the human APOBEC3 family of editing enzymes can inhibit various mobile genetic elements. APOBEC3A (A3A) can block the retrotransposon LINE-1 and the parvovirus adeno-associated virus type 2 (AAV-2) but does not inhibit retroviruses. In contrast, APOBEC3G (A3G) can block retroviruses but has only limited effects on AAV-2 or LINE-1. What dictates this differential target specificity remains largely undefined. Here, we modeled the structure of A3A based on its homology with the C-terminal domain of A3G and further compared the sequence of human A3A to those of 11 nonhuman primate orthologues. We then used these data to perform a mutational analysis of A3A, examining its ability to restrict LINE-1, AAV-2, and foreign plasmid DNA and to edit a single-stranded DNA substrate. The results revealed an essential functional role for the predicted single-stranded DNA-docking groove located around the A3A catalytic site. Within this region, amino acid differences between A3A and A3G are predicted to affect the shape of the polynucleotide-binding groove. Correspondingly, transferring some of these A3A residues to A3G endows the latter protein with the ability to block LINE-1 and AAV-2. These results suggest that the target specificity of APOBEC3 family members is partly defined by structural features influencing their interaction with polynucleotide substrates.

Effects of alpha-phosphoglucomutase deficiency on cell wall properties and fitness in Streptococcus gordonii.

Streptococcus gordonii alpha-phosphoglucomutase, which converts glucose 6-phosphate to glucose 1-phosphate, is encoded by pgm. The pgm transcript is monocistronic and is initiated from a sigma(A)-like promoter. Mutants with a gene disruption in pgm exhibited an altered cell wall muropeptide pattern and a lower teichoic acid content, and had reduced fitness both in vitro and in vivo. In vitro, the reduced fitness included reduced growth, reduced viability in the stationary phase and increased autolytic activity. In vivo, the pgm-deficient strain had a lower virulence in a rat model of experimental endocarditis.