Kallikrein genes are associated with lupus and glomerular basement membrane-specific antibody-induced nephritis in mice and humans.

Immune-mediated nephritis contributes to disease in systemic lupus erythematosus, Goodpasture syndrome (caused by antibodies specific for glomerular basement membrane [anti-GBM antibodies]), and spontaneous lupus nephritis. Inbred mouse strains differ in susceptibility to anti-GBM antibody-induced and spontaneous lupus nephritis. This study sought to clarify the genetic and molecular factors that maybe responsible for enhanced immune-mediated renal disease in these models. When the kidneys of 3 mouse strains sensitive to anti-GBM antibody-induced nephritis were compared with those of 2 control strains using microarray analysis, one-fifth of the underexpressed genes belonged to the kallikrein gene family,which encodes serine esterases. Mouse strains that upregulated renal and urinary kallikreins exhibited less evidence of disease. Antagonizing the kallikrein pathway augmented disease, while agonists dampened the severity of anti-GBM antibody-induced nephritis. In addition, nephritis-sensitive mouse strains had kallikrein haplotypes that were distinct from those of control strains, including several regulatory polymorphisms,some of which were associated with functional consequences. Indeed, increased susceptibility to anti-GBM antibody-induced nephritis and spontaneous lupus nephritis was achieved by breeding mice with a genetic interval harboring the kallikrein genes onto a disease-resistant background. Finally, both human SLE and spontaneous lupus nephritis were found to be associated with kallikrein genes, particularly KLK1 and the KLK3 promoter, when DNA SNPs from independent cohorts of SLE patients and controls were compared. Collectively, these studies suggest that kallikreins are protective disease-associated genes in anti-GBM antibody-induced nephritis and lupus.

An unexpected role of semaphorin3a-neuropilin-1 signaling in lymphatic vessel maturation and valve formation.

RATIONALE: Lymphatic vasculature plays important roles in tissue fluid homeostasis maintenance and in the pathology of human diseases. Yet, the molecular mechanisms that control lymphatic vessel maturation remain largely unknown. OBJECTIVE: We analyzed the gene expression profiles of ex vivo isolated lymphatic endothelial cells to identify novel lymphatic vessel expressed genes and we investigated the role of semaphorin 3A (Sema3A) and neuropilin-1 (Nrp-1) in lymphatic vessel maturation and function. METHODS AND RESULTS: Lymphatic and blood vascular endothelial cells from mouse intestine were isolated using fluorescence-activated cell sorting, and transcriptional profiling was performed. We found that the axonal guidance molecules Sema3A and Sema3D were highly expressed by lymphatic vessels. Importantly, we found that the semaphorin receptor Nrp-1 is expressed on the perivascular cells of the collecting lymphatic vessels. Treatment of mice in utero (E12.5-E16.5) with an antibody that blocks Sema3A binding to Nrp-1 but not with an antibody that blocks VEGF-A binding to Nrp-1 resulted in a complex phenotype of impaired lymphatic vessel function, enhanced perivascular cell coverage, and abnormal lymphatic vessel and valve morphology. CONCLUSIONS: Together, these results reveal an unanticipated role of Sema3A-Nrp-1 signaling in the maturation of the lymphatic vascular network likely via regulating the perivascular cell coverage of the vessels thus affecting lymphatic vessel function and lymphatic valve development.

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.

MsrR contributes to cell surface characteristics and virulence in Staphylococcus aureus.

MsrR, a factor contributing to methicillin resistance in Staphylococcus aureus, belongs to the LytR-CpsA-Psr family of cell envelope-associated proteins. Deletion of msrR increased cell size and aggregation, and altered envelope properties, leading to a temporary reduction in cell surface hydrophobicity, diminished colony-spreading ability, and an increased susceptibility to Congo red. The reduced phosphorus content of purified cell walls of the msrR mutant suggested a reduction in wall teichoic acids, which may explain some of the observed phenotypes. Microarray analysis of the msrR deletion mutant revealed only minor changes in the global transcriptome, suggesting that MsrR has structural rather than regulatory functions. Importantly, virulence of the msrR mutant was decreased in a nematode-killing assay as well as in rat experimental endocarditis. MsrR is therefore likely to play a role in cell envelope maintenance, cell separation, and pathogenicity of S. aureus.

Single cell analysis reveals similar functional competence of dominant and nondominant CD8 T-cell clonotypes.

Immune protection from infectious diseases and cancer is mediated by individual T cells of different clonal origin. Their functions are tightly regulated but not yet fully characterized. Understanding the contribution of each T cell will improve the prediction of immune protection based on laboratory assessment of T-cell responses. Here we developed techniques for simultaneous molecular and functional assessment of single CD8 T cells directly ex vivo. We studied two groups of patients with melanoma after vaccination with two closely related tumor antigenic peptides. Vaccination induced T cells with strong memory and effector functions, as found in virtually all T cells of the first patient group, and fractions of T cells in the second group. Interestingly, high functionality was not restricted to dominant clonotypes. Rather, dominant and nondominant clonotypes acquired equal functional competence. In parallel, this was also found for EBV- and CMV-specific T cells. Thus, the nondominant clonotypes may contribute similarly to immunity as their dominant counterparts.

Update of the FANTOM web resource: from mammalian transcriptional landscape to its dynamic regulation.

The international Functional Annotation Of the Mammalian Genomes 4 (FANTOM4) research collaboration set out to better understand the transcriptional network that regulates macrophage differentiation and to uncover novel components of the transcriptome employing a series of high-throughput experiments. The primary and unique technique is cap analysis of gene expression (CAGE), sequencing mRNA 5'-ends with a second-generation sequencer to quantify promoter activities even in the absence of gene annotation. Additional genome-wide experiments complement the setup including short RNA sequencing, microarray gene expression profiling on large-scale perturbation experiments and ChIP-chip for epigenetic marks and transcription factors. All the experiments are performed in a differentiation time course of the THP-1 human leukemic cell line. Furthermore, we performed a large-scale mammalian two-hybrid (M2H) assay between transcription factors and monitored their expression profile across human and mouse tissues with qRT-PCR to address combinatorial effects of regulation by transcription factors. These interdependent data have been analyzed individually and in combination with each other and are published in related but distinct papers. We provide all data together with systematic annotation in an integrated view as resource for the scientific community (http://fantom.gsc.riken.jp/4/). Additionally, we assembled a rich set of derived analysis results including published predicted and validated regulatory interactions. Here we introduce the resource and its update after the initial release.

Differences in osteogenic and apoptotic genes between osteoporotic and osteoarthritic patients.

BACKGROUND Osteoporosis is a metabolic disorder characterized by a reduction in bone mass and deterioration in the microarchitectural structure of the bone, leading to a higher risk for spontaneous and fragility fractures.The main aim was to study the differences between human bone from osteoporotic and osteoarthritic patients about gene expression (osteogenesis and apoptosis), bone mineral density, microstructural and biomechanic parameters. METHODS We analyzed data from 12 subjects: 6 with osteoporotic hip fracture (OP) and 6 with hip osteoarthritis (OA), as the control group. All subjects underwent medical history, analytical determinations, densitometry, histomorphometric and biochemical study. The expression of 86 genes of osteogenesis and 86 genes of apoptosis was studied in pool of bone samples from patients with OP and OA by PCR array. RESULTS We observed that most of the genes of apoptosis and osteogenesis show a decrease in gene expression in the osteoporotic group in comparison with the osteoarthritic group. The histomorphometric study shows a lower bone quality in the group of patients with hip fractures compared to the osteoarthritic group. CONCLUSIONS The bone tissue of osteoporotic fracture patients is more fragile than the bone of OA patients. Our results showed an osteoporotic bone with a lower capacities for differentiation and osteoblastic activity as well as a lower rate of apoptosis than osteoarthritic bone. These results are related with structural and biochemical parameters.

CYR61 and alphaVbeta5 integrin cooperate to promote invasion and metastasis of tumors growing in preirradiated stroma.

Radiotherapy is widely used to treat human cancer. Patients locally recurring after radiotherapy, however, have increased risk of metastatic progression and poor prognosis. The clinical management of postradiation recurrences remains an unresolved issue. Tumors growing in preirradiated tissues have an increased fraction of hypoxic cells and are more metastatic, a condition known as tumor bed effect. The transcription factor hypoxia inducible factor (HIF)-1 promotes invasion and metastasis of hypoxic tumors, but its role in the tumor bed effect has not been reported. Here, we show that tumor cells derived from SCCVII and HCT116 tumors growing in a preirradiated bed, or selected in vitro through repeated cycles of severe hypoxia, retain invasive and metastatic capacities when returned to normoxia. HIF activity, although facilitating metastatic spreading of tumors growing in a preirradiated bed, is not essential. Through gene expression profiling and gain- and loss-of-function experiments, we identified the matricellular protein CYR61 and alphaVbeta5 integrin as proteins cooperating to mediate these effects. The anti-alphaV integrin monoclonal antibody 17E6 and the small molecular alphaVbeta3/alphaVbeta5 integrin inhibitor EMD121974 suppressed invasion and metastasis induced by CYR61 and attenuated metastasis of tumors growing within a preirradiated field. These results represent a conceptual advance to the understanding of the tumor bed effect and identify CYR61 and alphaVbeta5 integrin as proteins that cooperate to mediate metastasis. They also identify alphaV integrin inhibition as a potential therapeutic approach for preventing metastasis in patients at risk for postradiation recurrences.

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.

Long-range heterogeneity at the 3' ends of human mRNAs.

The publication of a draft of the human genome and of large collections of transcribed sequences has made it possible to study the complex relationship between the transcriptome and the genome. In the work presented here, we have focused on mapping mRNA 3' ends onto the genome by use of the raw data generated by the expressed sequence tag (EST) sequencing projects. We find that at least half of the human genes encode multiple transcripts whose polyadenylation is driven by multiple signals. The corresponding transcript 3' ends are spread over distances in the kilobase range. This finding has profound implications for our understanding of gene expression regulation and of the diversity of human transcripts, for the design of cDNA microarray probes, and for the interpretation of gene expression profiling experiments.

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.

A combined molecular clinical predictor of survival validated with the rtog-0525 cohort

For glioblastoma (GBM), survival classification has primarily relied on clinical criteria, exemplified by the Radiation Therapy Oncology Group (RTOG) recursive partitioning analysis (RPA). We sought to improve tumor classification by combining tumor biomarkers with the clinical RPA data. To accomplish this, we first developed 4 molecular biomarkers derived from gene expression profiling, a glioma CpG island methylator phenotype, a novel MGMT promoter methylation assay, and IDH1 mutations. A molecular predictor (MP) model was created with these 4 biomarkers on a training set of 220 retrospectively collected archival GBMtumors. ThisMPwas further combined with RPA classification to develop a molecular-clinical predictor (MCP). The median survivals for the combined, 4-class MCP were 65 months, 31 months, 13 months, and 9 months, which was significantly improved when compared with the RPA alone. The MCP was then applied to 725 samples from the RTOG-0525 cohort, showing median survival for each risk group of NR, 26 months, 16 months, and 11 months. The MCP was significantly improved over the RPA at outcome prediction in the RTOG 0525 cohort with a 33%increase in explained variation with respect to survival, validating the result obtained in the training set. To illustrate the benefit of the MCP for patient stratification, we examined progression-free survival (PFS) for patients receiving standard-dose temozolomide (SD-TMZ) vs. dose-dense TMZ (DD-TMZ) in RPA and MCP risk groups. A significant difference between DD-TMZ and SD-TMZ was observed in the poorest surviving MCP risk group with a median PFS of 6 months vs. 3 months (p ¼ 0.048, log-rank test). This difference was not seen using the RPA classification alone. In summary, we have developed a combined molecular-clinical predictor that appears to improve outcome prediction when compared with clinical variables alone. This MCP may serve to better identify patients requiring intensive treatments beyond the standard of care.

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.

CD30-positive peripheral T-cell lymphomas share molecular and phenotypic features.

Peripheral T-cell lymphoma, not otherwise specified is a heterogeneous group of aggressive neoplasms with indistinct borders. By gene expression profiling we previously reported unsupervised clusters of peripheral T-cell lymphomas, not otherwise specified correlating with CD30 expression. In this work we extended the analysis of peripheral T-cell lymphoma molecular profiles to prototypical CD30(+) peripheral T-cell lymphomas (anaplastic large cell lymphomas), and validated mRNA expression profiles at the protein level. Existing transcriptomic datasets from peripheral T-cell lymphomas, not otherwise specified and anaplastic large cell lymphomas were reanalyzed. Twenty-one markers were selected for immunohistochemical validation on 80 peripheral T-cell lymphoma samples (not otherwise specified, CD30(+) and CD30(-); anaplastic large cell lymphomas, ALK(+) and ALK(-)), and differences between subgroups were assessed. Clinical follow-up was recorded. Compared to CD30(-) tumors, CD30(+) peripheral T-cell lymphomas, not otherwise specified were significantly enriched in ALK(-) anaplastic large cell lymphoma-related genes. By immunohistochemistry, CD30(+) peripheral T-cell lymphomas, not otherwise specified differed significantly from CD30(-) samples [down-regulated expression of T-cell receptor-associated proximal tyrosine kinases (Lck, Fyn, Itk) and of proteins involved in T-cell differentiation/activation (CD69, ICOS, CD52, NFATc2); upregulation of JunB and MUM1], while overlapping with anaplastic large cell lymphomas. CD30(-) peripheral T-cell lymphomas, not otherwise specified tended to have an inferior clinical outcome compared to the CD30(+) subgroups. In conclusion, we show molecular and phenotypic features common to CD30(+) peripheral T-cell lymphomas, and significant differences between CD30(-) and CD30(+) peripheral T-cell lymphomas, not otherwise specified, suggesting that CD30 expression might delineate two biologically distinct subgroups.

Integrative analysis of RUNX1 downstream pathways and target genes.

BACKGROUND: The RUNX1 transcription factor gene is frequently mutated in sporadic myeloid and lymphoid leukemia through translocation, point mutation or amplification. It is also responsible for a familial platelet disorder with predisposition to acute myeloid leukemia (FPD-AML). The disruption of the largely unknown biological pathways controlled by RUNX1 is likely to be responsible for the development of leukemia. We have used multiple microarray platforms and bioinformatic techniques to help identify these biological pathways to aid in the understanding of why RUNX1 mutations lead to leukemia. RESULTS: Here we report genes regulated either directly or indirectly by RUNX1 based on the study of gene expression profiles generated from 3 different human and mouse platforms. The platforms used were global gene expression profiling of: 1) cell lines with RUNX1 mutations from FPD-AML patients, 2) over-expression of RUNX1 and CBFbeta, and 3) Runx1 knockout mouse embryos using either cDNA or Affymetrix microarrays. We observe that our datasets (lists of differentially expressed genes) significantly correlate with published microarray data from sporadic AML patients with mutations in either RUNX1 or its cofactor, CBFbeta. A number of biological processes were identified among the differentially expressed genes and functional assays suggest that heterozygous RUNX1 point mutations in patients with FPD-AML impair cell proliferation, microtubule dynamics and possibly genetic stability. In addition, analysis of the regulatory regions of the differentially expressed genes has for the first time systematically identified numerous potential novel RUNX1 target genes. CONCLUSION: This work is the first large-scale study attempting to identify the genetic networks regulated by RUNX1, a master regulator in the development of the hematopoietic system and leukemia. The biological pathways and target genes controlled by RUNX1 will have considerable importance in disease progression in both familial and sporadic leukemia as well as therapeutic implications