New target genes in endometrial tumors show a role for the estrogen-receptor pathway in microsatellite-unstable cancers.

Microsatellite instability (MSI) in tumors results in an accumulation of mutations in (target) genes. Previous studies suggest that the profile of target genes differs according to tumor type. This paper describes the first genome-wide search for target genes for mismatch repair-deficient endometrial cancers. Genes expressed in normal endometrium containing coding repeats were analyzed for mutations in tumors. We identified 44 possible genes of which seven are highly mutated (>15%). Some candidates were also found mutated in colorectal and gastric tumors. The most frequently mutated gene, NRIP1 encoding nuclear receptor-interacting protein 1, was silenced in an endometrial tumor cell line and expression microarray experiments were performed. Silencing of NRIP1 was associated with differences in the expression of several genes in the estrogen-receptor network. Furthermore, an enrichment of genes related to cell cycle (regulation) and replication was observed. We present a new profile of target genes, some of them tissue specific, whereas others seem to play a more general role in MSI tumors. The high-mutation frequency combined with the expression data suggest, for the first time, an involvement of NRIP1 in endometrial cancer development.

Organization of the variable region of the immunoglobulin heavy-chain gene locus of the rat.

We have mapped and annotated the variable region of the immunoglobulin heavy (IGH) gene locus of the Brown Norway (BN) rat (assembly V3.4; Rat Genomic Sequence Consortium). In addition to known variable region genes, we found 12 novel previously unidentified functional IGHV genes and 1 novel functional IGHD gene. In total, the variable region of the rat IGH locus is composed of at least 353 unique IGHV genes, 21 IGHD genes, and 5 IGHJ genes, of which 131, 14, and 4 are potentially functional genes, respectively. Of all species studied so far, the rat seems to have the highest number of functional IGHV genes in the genome. Rat IGHV genes can be classified into 13 IGHV families based on nucleotide sequence identity. The variable region of the BN rat spans a total length of approximately 4.9 Mb and is organized in a typical translocon organization. Like the mouse, members of the various IGHV gene families are more or less grouped together on the genome, albeit some members of IGHV gene families are found intermingled with each other. In the rat, the largest IGHV gene families are IGHV1, IGHV2, and IGHV5. The overall conclusion is that the genomic organization of the variable region of the rat IGH locus is strikingly similar to that of the mouse, illustrating the close evolutionary relationship between these two species.

PROS1 analysis in 87 pedigrees with hereditary protein S deficiency demonstrates striking genotype-phenotype associations.

Hereditary protein S (PS) deficiency predisposes to venous thrombosis. Previously, we demonstrated a difference in risk of venous thrombosis between PS deficiency type I and type III. We used direct sequencing, multiplex ligation-dependent probe amplification (MLPA), and linkage analysis to study whether this difference could be explained by molecular heterogeneity. The study contained two sets of families with PS deficiency type I (cohort 1; 35 probands, 155 relatives) or type III (cohort 2; 52 probands, 241 relatives). In cohort 1, a mixed type I/type III PS-deficient phenotype was observed in 66% of the pedigrees. A total of 34 probands carried a mutant PROS1 allele, compared to one proband in cohort 2 (P<10(-10)). The proband's mutation was identified in all type I, but only in 57% of type III PS deficient relatives. MLPA-analysis in the mutation negative families did not reveal PROS1 deletions or insertions. Linkage analysis in 16 families showed cosegregation of PROS1 markers in the family with type I deficiency, but not in the 15 families with type III deficiency. The genotype-phenotype associations point to differences in genetic architecture. Whereas PS deficiency type I is a monogenic disease due to PROS1 allelic heterozygosity, PS deficiency type III is most likely a more complex or heterogeneous disorder.

A novel lipid-based drug carrier targeted to the non-parenchymal cells, including hepatic stellate cells, in the fibrotic livers of bile duct ligated rats.

In fibrotic livers, collagen producing hepatic stellate cells (HSC) represent a major target for antifibrotic therapies. We designed liposomes with surface-coupled mannose 6-phosphate (M6P) modified human serum albumin (HSA) to target HSC via the M6P receptor. In this study we determined the pharmacokinetics and target specificity of M6P-HSA-liposomes in a rat model of liver fibrosis. Ten minutes after injection of [(3)H]-M6P-HSA-liposomes 90% of the dose has cleared the circulation. The blood elimination of these liposomes was counteracted by free M6P-HSA and polyinosinic acid, a competitive inhibitor of scavenger receptors. The M6P-HSA-liposomes accumulated in HSC. However, also Kupffer cells and endothelial cells contributed to the uptake of M6P-HSA-liposomes in the fibrotic livers. Polyinosinic acid inhibited the accumulation of the liposomes in Kupffer cells and liver endothelial cells, but not in HSC. PCR analysis revealed that cultured HSC express scavenger receptors. This was confirmed by Western blotting, although activation of HSC diminishes scavenger receptor protein expression. In conclusion, in a rat model for liver fibrosis M6P-HSA-liposomes can be efficiently targeted to non-parenchymal cells, including HSC. M6P receptors and scavenger receptors are involved in the cellular recognition of these liposomes, allowing multiple pharmacological interference in different pathways involved in the fibrosis.

A verification protocol for the probe sequences of Affymetrix genome arrays reveals high probe accuracy for studies in mouse, human and rat.

BACKGROUND: The Affymetrix GeneChip technology uses multiple probes per gene to measure its expression level. Individual probe signals can vary widely, which hampers proper interpretation. This variation can be caused by probes that do not properly match their target gene or that match multiple genes. To determine the accuracy of Affymetrix arrays, we developed an extensive verification protocol, for mouse arrays incorporating the NCBI RefSeq, NCBI UniGene Unique, NIA Mouse Gene Index, and UCSC mouse genome databases. RESULTS: Applying this protocol to Affymetrix Mouse Genome arrays (the earlier U74Av2 and the newer 430 2.0 array), the number of sequence-verified probes with perfect matches was no less than 85% and 95%, respectively; and for 74% and 85% of the probe sets all probes were sequence verified. The latter percentages increased to 80% and 94% after discarding one or two unverifiable probes per probe set, and even further to 84% and 97% when, in addition, allowing for one or two mismatches between probe and target gene. Similar results were obtained for other mouse arrays, as well as for human and rat arrays. Based on these data, refined chip definition files for all arrays are provided online. Researchers can choose the version appropriate for their study to (re)analyze expression data. CONCLUSION: The accuracy of Affymetrix probe sequences is higher than previously reported, particularly on newer arrays. Yet, refined probe set definitions have clear effects on the detection of differentially expressed genes. We demonstrate that the interpretation of the results of Affymetrix arrays is improved when the new chip definition files are used.

A statistical multiprobe model for analyzing cis and trans genes in genetical genomics experiments with short-oligonucleotide arrays.

Short-oligonucleotide arrays typically contain multiple probes per gene. In genetical genomics applications a statistical model for the individual probe signals can help in separating "true" differential mRNA expression from "ghost" effects caused by polymorphisms, misdesigned probes, and batch effects. It can also help in detecting alternative splicing, start, or termination.

Changes in winter depression phenotype correlate with white blood cell gene expression profiles: a combined metagene and gene ontology approach.

In the present study we evaluate the feasibility of gene expression in white blood cells as a peripheral marker for winter depression. Sixteen patients with winter type seasonal affective disorder were included in the study. Blood was taken by venous puncture at three time points; in winter prior and following bright light therapy and in summer. RNA was isolated, converted into cRNA, amplified and hybridized on Illumina® gene expression arrays. The raw optical array data were quantile normalized and thereafter analyzed using a metagene approach, based on previously published Affymetrix gene array data. The raw data were also subjected to a secondary analysis focusing on circadian genes and genes involved in serotonergic neurotransmission. Differences between the conditions were analyzed, using analysis of variance on the principal components of the metagene score matrix. After correction for multiple testing no statistically significant differences were found. Another approach uses the correlation between metagene factor weights and the actual expression values, averaged over conditions. When comparing the correlations of winter vs. summer and bright light therapy vs. summer significant changes for several metagenes were found. Subsequent gene ontology analyses (DAVID and GeneTrail) of 5 major metagenes suggest an interaction between brain and white blood cells. The hypothesis driven analysis with a smaller group of genes failed to demonstrate any significant effects. The results from the combined metagene and gene ontology analyses support the idea of communication between brain and white blood cells. Future studies will need a much larger sample size to obtain information at the level of single genes.

Use of the EGP-2/Ep-CAM promoter for targeted expression of heterologous genes in carcinoma derived cell lines.

EGP-2, also known as Ep-CAM, is expressed at high levels on the surface of most carcinomas and is therefore considered an attractive target for anticancer strategies. To explore the mechanisms regulating the expression of EGP-2, sequences 3.4 kb upstream of the transcription start site were isolated and assayed for their ability to control the expression of the EGP-2 cDNA, the green fluorescent protein, the luciferase reporter gene and the thymidine kinase and cytosine deaminase suicide genes. Expression of these chimeric constructs as assessed in a range of different cell lines was restricted to cell lines expressing EGP-2. In addition, only cells expressing EGP-2 were sensitive for gancyclovir after being transiently transfected with EGP-2 promoter-driven thymidine kinase. Deletion analyses defined 687 bp upstream as the basic proximal promoter region, which could confer epithelial-specific expression to the GFP reporter gene in vitro. As these EGP-2 sequences can confer promoter activity to reporter and suicide genes in an EGP-2 restricted manner, they may be useful for gene therapy of EGP-2 expressing carcinomas.

EpCAM homologues exhibit epithelial-specific but different expression patterns in the kidney.

The Epithelial Cell Adhesion Molecule (EpCAM) is expressed virtually on normal epithelia in vertebrates. Among different species, the amino acid sequence of EpCAM is highly homologous, indicating that EpCAM is an evolutionary conserved protein. However, differences in the expression pattern of EpCAM homologues have been reported. We hypothesized that differences in expression pattern might be related to the promoter organization of the respective EpCAM homologues. Therefore, we here compared the promoter region of the mouse and human EpCAM homologues. In addition, we compared the expression pattern of the human and murine EpCAM homologues in the hEpCAM transgenic mouse. In silico analysis of EpCAM homologues revealed that the amino acid sequence as well as the domain structure is highly conserved throughout different vertebrates. In silico analysis of the promoter region identified that the human and mouse EpCAM promoters share a low homology. In agreement with this low homology, the murine and human EpCAM promoter contains only a few common transcription factor binding sites. Nevertheless, immunohistochemcial analysis of the expression of human and murine EpCAM in lung, colon, and kidney of the hEpCAM transgenic mouse identified that expression of both homologues is restricted to epithelial cells in these organs. Moreover, in lung and colon the human and murine homologues of EpCAM were co-expressed. In contrast, the EpCAM homologues were only sporadically co-expressed in renal epithelia, although they were distributed similarly along the nephronic segments. Together, these findings indicate an overall conserved regulatory mechanism that ensures epithelial expression of EpCAM homologues, despite the low promoter homology. Furthermore, the fact that murine epithelia express the human homologue of EpCAM indicates that the mouse has transcription factors required for human EpCAM expression.

Sequence polymorphisms cause many false cis eQTLs.

Many investigations have reported the successful mapping of quantitative trait loci (QTLs) for gene expression phenotypes (eQTLs). Local eQTLs, where expression phenotypes map to the genes themselves, are of especially great interest, because they are direct candidates for previously mapped physiological QTLs. Here we show that many mapped local eQTLs in genetical genomics experiments do not reflect actual expression differences caused by sequence polymorphisms in cis-acting factors changing mRNA levels. Instead they indicate hybridization differences caused by sequence polymorphisms in the mRNA region that is targeted by the microarray probes. Many such polymorphisms can be detected by a sensitive and novel statistical approach that takes the individual probe signals into account. Applying this approach to recent mouse and human eQTL data, we demonstrate that indeed many local eQTLs are falsely reported as "cis-acting" or "cis" and can be successfully detected and eliminated with this approach.

Functional analysis of lung tumor suppressor activity at 3p21.3.

The early and frequent occurrence of deletions at 3p21.3 in lung cancer has led to the consideration of this chromosomal region as a lung cancer (LUCA) critical region with tumor suppressor activity. We covered this 19 genes-containing region with overlapping P1 artificial chromosomes (PACs), in which genes are likely accompanied by their own promoters or other regulatory sequences. With these PACs we transfected cells from a small cell lung cancer (SCLC) cell line which readily caused tumors in nude mice. Per PAC we selected two cell clones with a low number of PAC copies integrated at a single genomic site. The selected clones were s.c. injected into nude mice to investigate whether the integrated genes suppressed the tumor-inducing capacity of the original SCLC cell line. We could demonstrate PAC-specific gene expression in the transfected cells. All of the PAC integration sites were different. It appeared that introduction of a PAC or even an empty PAC vector causes some chromosomal instability, which in principle may either promote or inhibit cell growth. However, both cell clones with integration of the same PAC from the centromeric part of the LUCA region in different genomic sites were the sole pair of clones that caused smaller tumors than did the original SCLC cell line. This suggests that rather than the induced chromosomal instability, the DNA sequence of that PAC, which in addition to two protein-encoding genes contains at least one potential miRNA gene, is responsible for the tumor suppressor activity.

Bacterial degradation of xenobiotic compounds: evolution and distribution of novel enzyme activities.

Bacterial dehalogenases catalyse the cleavage of carbon-halogen bonds, which is a key step in aerobic mineralization pathways of many halogenated compounds that occur as environmental pollutants. There is a broad range of dehalogenases, which can be classified in different protein superfamilies and have fundamentally different catalytic mechanisms. Identical dehalogenases have repeatedly been detected in organisms that were isolated at different geographical locations, indicating that only a restricted number of sequences are used for a certain dehalogenation reaction in organohalogen-utilizing organisms. At the same time, massive random sequencing of environmental DNA, and microbial genome sequencing projects have shown that there is a large diversity of dehalogenase sequences that is not employed by known catabolic pathways. The corresponding proteins may have novel functions and selectivities that could be valuable for biotransformations in the future. Apparently, traditional enrichment and metagenome approaches explore different segments of sequence space. This is also observed with alkane hydroxylases, a category of proteins that can be detected on basis of conserved sequence motifs and for which a large number of sequences has been found in isolated bacterial cultures and genomic databases. It is likely that ongoing genetic adaptation, with the recruitment of silent sequences into functional catabolic routes and evolution of substrate range by mutations in structural genes, will further enhance the catabolic potential of bacteria toward synthetic organohalogens and ultimately contribute to cleansing the environment of these toxic and recalcitrant chemicals.

The 172 kb prkA-addAB region from 83 degrees to 97 degrees of the Bacillus subtilis chromosome contains several dysfunctional genes, the glyB marker, many genes encoding transporter proteins, and the ubiquitous hit gene.

A 171812 bp nucleotide sequence between prkA and addAB (83 degrees to 97 degrees) on the genetic map of the Bacillus subtilis 168 chromosome was determined and analysed. An accurate physical/genetic map of this previously poorly described chromosomal region was constructed. One hundred and seventy open reading frames (ORFs) were identified on the DNA fragment. These include the previously described genes cspB, glpPFKD, spoVR, phoAIV, papQ, citRA, sspB, prsA, hpr, pbpF, hemEHY, aprE, comK and addAB. ORF yhaF in this region corresponds to the glyB marker. Among the striking features of this region are: an abundance of genes encoding (putative) transporter proteins, several dysfunctional genes, the ubiquitous hit gene, and five multidrug-resistance-like genes. These analyses have also revealed the existence of numerous paralogues of ORFs in this region: about two-thirds of the putative genes seem to have at least one paralogue in the B. subtilis genome.

Structural and functional genomics of the CPT1B gene for muscle-type carnitine palmitoyltransferase I in mammals.

Muscle-type carnitine palmitoyltransferase I (M-CPT I) is a key enzyme in the control of beta-oxidation of long-chain fatty acids in the heart and skeletal muscle. Because knowledge of the mammalian genes encoding M-CPT I may aid in studies of disturbed energy metabolism, we obtained new genomic and cDNA data for M-CPT I for the human, mouse, rat, and sheep. The introns of these compact genes are 80% (mouse versus rat) and 60% (mouse versus human) identical. Sheep and goat, but not cow, pig, rodent, or human promoter sequences contain a short interspersed repeated sequence (SINE) upstream of highly conserved regulatory elements. These elements constitute two promoters in humans, sheep, and mice, and, contrary to previous reports, there is a second promoter in rats as well. Thus, the transcriptional organization of these genes is more uniform than previously supposed, with interspecies differences in the 5'-ends of the mRNAs reflecting differences in splicing; only in humans extensive splicing and splice variation is found in the 5'- and 3'-untranslated regions. In the mouse, intron retention was detected in heart, muscle, and testes and may indicate an additional mechanism of regulation of M-CPT I expression. Splice variation in the coding region was previously proposed to lead to expression of CPT I enzymes with altered malonyl-CoA sensitivity (Yu, G. S., Lu, Y. C., and Gulick, T. (1998) Biochem. J. 334, 225-231). However, when expressed in the yeast Pichia pastoris, none of three earlier described splice variants had CPT I activity. Therefore, the involvement of splice variation of M-CPT I in the modulation of malonyl-CoA inhibition of fatty acid oxidation may be less relevant than hitherto assumed.