Verification of microRNA expression in human endometrial adenocarcinoma

BACKGROUND

MicroRNAs are small non-coding RNAs that have been implicated in tumor initiation and progression. In a previous study we identified 138 miRNAs as differentially expressed in endometrial adenocarcinoma compared to normal tissues. One of these miRNAs was miRNA-34a, which regulates several genes involved in the Notch pathway, which is frequently altered in endometrial cancer. The aims of this study were to verify the differential expression of a subset of miRNAs and to scrutinize the regulatory role of mir-34a on the target genes NOTCH1 and DLL1.

METHODS

Twenty-five miRNAs that were previously identified as differentially expressed were subjected to further analysis using qPCR. To investigate the regulation of NOTCH1 and DLL1 by mir-34a, we designed gain- and loss-of-function experiments in Ishikawa and HEK293 cell lines by transfection with a synthetic mir-34a mimic and a mir-34a inhibitor.

RESULTS

Of the 25 validated miRNAs, seven were down-regulated and 18 were up-regulated compared to normal endometrium, which was fully consistent with our previous findings. In addition, the up-regulation of mir-34a led to a significant decrease in mRNA levels of NOTCH1 and DLL1, while down-regulation led to a significant increase in mRNA levels of these two genes.

CONCLUSIONS

We verified both up-regulated and down-regulated miRNAs in the tumor samples, indicating various roles of microRNAs during tumor development. Mir-34a functions as a regulator by decreasing the expression of NOTCH1 and DLL1. Our study is the first to identify a correlation between mir-34a and its target genes NOTCH1 and DLL1 in endometrial adenocarcinoma.

Genome-wide discovery of miRNAs using ensembles of machine learning algorithms and logistic regression

In silico prediction of novel miRNAs from genomic sequences remains a challenging problem. This study presents a genome-wide miRNA discovery software package called GenoScan and evaluates two hairpin classification methods. These methods, one ensemble-based and one using logistic regression were benchmarked along with 15 published methods. In addition, the sequence-folding step is addressed by investigating the impact of secondary structure prediction methods and the choice of input sequence length on prediction performance. Both the accuracy of secondary structure predictions and the miRNA prediction are evaluated. In the benchmark of hairpin classification methods, the regression model achieved highest classification accuracy. Of the structure prediction methods evaluated, ContextFold achieved the highest agreement between predicted and experimentally determined structures. However, both the choice of secondary structure prediction method and input sequence length had limited impact on hairpin classification performance.

Expression of miRNAs in Bull Spermatozoa Correlates with Fertility Rates

Bull spermatozoa are rich in active miRNAs, and it has been shown that specific spermborne miRNAs can be linked to fertility. Thus, expression profiling of spermatozoa could be helpful for understanding male fertility and the ability of spermatozoa to initiate and sustain zygotic, embryonic and foetal development. Herein we hypothesized that bulls with moderate to high fertility can be identified by differences in amounts of certain miRNAs between their ejaculates. RNA samples from spermatozoa of eight brother pairs (one bull with high and one with moderate NRR in each pair) of the Holstein breed were prepared. miRNA was isolated, and the expression of 178 miRNAs was determined by RT-qPCR. Important findings were that highly expressed miRNAs, not linked to NRR status, were identified in the bull sperm samples, which indicate that these miRNAs have an important role in early embryogenesis. A large fraction of the targets genes were phosphoproteins and genes involved in the regulation of transcription. Seven miRNAs (mir-502-5p, mir-1249, mir-320a, mir-34c-3p, mir-19b-3p, mir-27a-5p and mir-148b-3p) were differentially expressed between bulls with moderate and high NRR with a strong tendency towards a higher expression of miRNAs in bulls with moderate fertility. Thus, bulls with a moderate NRR negatively regulate the expression of protein-coding genes, which leads to problems during the pregnancy.

miREC: a database of miRNAs involved in the development of endometrial cancer

BACKGROUND

Endometrial cancer (EC) is the most frequently diagnosed gynecological malignancy and the fourth most common cancer diagnosis overall among women. As with many other forms of cancer, it has been shown that certain miRNAs are differentially expressed in EC and these miRNAs are believed to play important roles as regulators of processes involved in the development of the disease. With the rapidly growing number of studies of miRNA expression in EC, there is a need to organize the data, combine the findings from experimental studies of EC with information from various miRNA databases, and make the integrated information easily accessible for the EC research community.

FINDINGS

The miREC database is an organized collection of data and information about miRNAs shown to be differentially expressed in EC. The database can be used to map connections between miRNAs and their target genes in order to identify specific miRNAs that are potentially important for the development of EC. The aim of the miREC database is to integrate all available information about miRNAs and target genes involved in the development of endometrial cancer, and to provide a comprehensive, up-to-date, and easily accessible source of knowledge regarding the role of miRNAs in the development of EC. Database URL: http://www.mirecdb.org .

CONCLUSIONS

Several databases have been published that store information about all miRNA targets that have been predicted or experimentally verified to date. It would be a time-consuming task to navigate between these different data sources and literature to gather information about a specific disease, such as endometrial cancer. The miREC database is a specialized data repository that, in addition to miRNA target information, keeps track of the differential expression of genes and miRNAs potentially involved in endometrial cancer development. By providing flexible search functions it becomes easy to search for EC-associated genes and miRNAs from different starting points, such as differential expression and genomic loci (based on genomic aberrations).

MicroRNA expression in human endometrial adenocarcinoma

Background

MicroRNAs are small non-coding RNAs that play crucial roles in the pathogenesis of different cancer types. The aim of this study was to identify miRNAs that are differentially expressed in endometrial adenocarcinoma compared to healthy endometrium. These miRNAs can potentially be used to develop a panel for classification and prognosis in order to better predict the progression of the disease and facilitate the choice of treatment strategy.

Methods

Formalin fixed paraffin embedded endometrial tissue samples were collected from the Örebro university hospital. QPCR was used to quantify the expression levels of 742 miRNAs in 30 malignant and 20 normal endometrium samples. After normalization of the qPCR data, miRNAs differing significantly in expression between normal and cancer samples were identified, and hierarchical clustering analysis was used to identify groups of miRNAs with coordinated expression profiles.

Results

In comparisons between endometrial adenocarcinoma and normal endometrium samples 138 miRNAs were found to be significantly differentially expressed (p < 0.001) among which 112 miRNAs have not been previous reported for endometrial adenocarcinoma.

Conclusion

Our study shows that several miRNAs are differentially expressed in endometrial adenocarcinoma. These identified miRNA hold great potential as target for classification and prognosis of this disease. Further analysis of the differentially expressed miRNA and their target genes will help to derive new biomarkers that can be used for classification and prognosis of endometrial adenocarcinoma.

Electronic supplementary material

The online version of this article (doi:10.1186/s12935-014-0088-6) contains supplementary material, which is available to authorized users.

Differences in microRNA expression during tumor development in the transition and peripheral zones of the prostate

Background

The prostate is divided into three glandular zones, the peripheral zone (PZ), the transition zone (TZ), and the central zone. Most prostate tumors arise in the peripheral zone (70-75%) and in the transition zone (20-25%) while only 10% arise in the central zone. The aim of this study was to investigate if differences in miRNA expression could be a possible explanation for the difference in propensity of tumors in the zones of the prostate.

Methods

Patients with prostate cancer were included in the study if they had a tumor with Gleason grade 3 in the PZ, the TZ, or both (n=16). Normal prostate tissue was collected from men undergoing cystoprostatectomy (n=20). The expression of 667 unique miRNAs was investigated using TaqMan low density arrays for miRNAs. Student’s t-test was used in order to identify differentially expressed miRNAs, followed by hierarchical clustering and principal component analysis (PCA) to study the separation of the tissues. The ADtree algorithm was used to identify markers for classification of tissues and a cross-validation procedure was used to test the generality of the identified miRNA-based classifiers.

Results

The t-tests revealed that the major differences in miRNA expression are found between normal and malignant tissues. Hierarchical clustering and PCA based on differentially expressed miRNAs between normal and malignant tissues showed perfect separation between samples, while the corresponding analyses based on differentially expressed miRNAs between the two zones showed several misplaced samples. A classification and cross-validation procedure confirmed these results and several potential miRNA markers were identified.

Conclusions

The results of this study indicate that the major differences in the transcription program are those arising during tumor development, rather than during normal tissue development. In addition, tumors arising in the TZ have more unique differentially expressed miRNAs compared to the PZ. The results also indicate that separate miRNA expression signatures for diagnosis might be needed for tumors arising in the different zones. MicroRNA signatures that are specific for PZ and TZ tumors could also lead to more accurate prognoses, since tumors arising in the PZ tend to be more aggressive than tumors arising in the TZ.

Validation of suitable endogenous control genes for quantitative PCR analysis of microRNA gene expression in a rat model of endometrial cancer

BACKGROUND

MicroRNAs are small RNA molecules that negatively regulate gene expression by translational inhibition or mRNA cleavage. The discovery that abnormal expression of particular miRNAs contributes to human disease, including cancer, has spurred growing interest in analysing expression profiles of these molecules. Quantitative polymerase chain reaction is frequently used for quantification of miRNA expression due to its sensitivity and specificity. To minimize experimental error in this system an appropriate endogenous control gene must be chosen. An ideal endogenous control gene should be expressed at a constant level across all samples and its expression stability should be unaffected by the experimental procedure.

RESULTS

The expression and validation of candidate control genes (4.5S RNA(H) A, Y1, 4.5S RNA(H) B, snoRNA, U87 and U6) was examined in 21 rat cell lines to establish the most suitable endogenous control for miRNA analysis in a rat model of cancer. The stability of these genes was analysed using geNorm and NormFinder algorithms. U87 and snoRNA were identified as the most stable control genes, while Y1 was least stable.

CONCLUSION

This study identified the control gene that is most suitable for normalizing the miRNA expression data in rat. That reference gene will be useful when miRNAs expression are analyzed in order to find new miRNA markers for endometrial cancer in rat.

Expression patterns of Phf5a/PHF5A and Gja1/GJA1 in rat and human endometrial cancer

Endometrial adenocarcinoma is the most frequently diagnosed cancer of the female genital tract in the western world. Studies of complex diseases can be difficult to perform on human tumor samples due to the high genetic heterogeneity in human. The use of rat models is preferable since rat has similarities in pathogenesis and histopathological properties to that of human.

A genomic region including the highly conserved Phf5a gene associated to development of EAC has previously been identified in an association study. PHF5A has been suggested to acts as a transcription factor or cofactor in the up regulation of expression of Gja1 gene in the presence of estrogen. It has earlier been shown that the Phf5a gene is down regulated in rat EAC derived cell lines by means of expression microarrays.

We analyzed the expression of Phf5a and Gja1 by qPCR, and potential relations between the two genes in EAC tumors and non-malignant cell lines derived from the BDII rat model. In addition, the expression pattern of these genes was compared in rat and human EAC tumor samples.

Changes in expression for Phf5a/PHF5A were found in tumors from both rat and human even though the observed pattern was not completely consistent between the two species. By separating rat EAC cell lines according to the genetic background, a significant lower expression of Phf5a in one of the two cross backgrounds was revealed, but not for the other. In contrast to other studies, Phf5a/PHF5A regulation of Gja1/GJA1 was not revealed in this study.

Classification of tumor samples from expression data using decision trunks

We present a novel machine learning approach for the classification of cancer samples using expression data. We refer to the method as “decision trunks,” since it is loosely based on decision trees, but contains several modifications designed to achieve an algorithm that: (1) produces smaller and more easily interpretable classifiers than decision trees; (2) is more robust in varying application scenarios; and (3) achieves higher classification accuracy. The decision trunk algorithm has been implemented and tested on 26 classification tasks, covering a wide range of cancer forms, experimental methods, and classification scenarios. This comprehensive evaluation indicates that the proposed algorithm performs at least as well as the current state of the art algorithms in terms of accuracy, while producing classifiers that include on average only 2–3 markers. We suggest that the resulting decision trunks have clear advantages over other classifiers due to their transparency, interpretability, and their correspondence with human decision-making and clinical testing practices.

The impact of the genetic background on the genome make-up of tumor cells

Endometrial adenocarcinoma (EAC) is the most common form of malignancy in the female genital tract, ranking as the fourth leading form of invasive tumors that affect women. The BDII inbred rat strain has been used as a powerful tumor model in studies of the genetic background of EAC. Females from the BDII strain are prone to develop tumors with an incidence of more than 90%. Development of EAC in BDII female rats has similarities in pathogenesis, histopathological, and molecular properties to that of human, and thus represents a unique model for analysis of EAC tumorigenesis and for comparative studies in human EACs. In a previous study, a set of rat EAC cell lines derived from tumors developed in female crossprogenies between BDII and nonsusceptible rat strains were analyzed by spectral karyotyping (SKY). Here we present an analysis with specific focus on the impact of different genetic backgrounds on the rate and occurrence of genetic aberrations in experimental tumors using data presented in the previous report. We could reveal that the ploidy state, and the abundance and type of structural as well as numerical change differed between the two genetic setups. We have also identified chromosomes harboring aberrations independent of genetic input from the nonsusceptible strains, which provide valuable information for the identification of the genes involved in the development of EAC in the BDII model as well as in human endometrial tumors.

SKY analysis revealed recurrent numerical and structural chromosome changes in BDII rat endometrial carcinomas

Background

Genomic alterations are common features of cancer cells, and some of these changes are proven to be neoplastic-specific. Such alterations may serve as valuable tools for diagnosis and classification of tumors, prediction of clinical outcome, disease monitoring, and choice of therapy as well as for providing clues to the location of crucial cancer-related genes.

Endometrial carcinoma (EC) is the most frequently diagnosed malignancy of the female genital tract, ranking fourth among all invasive tumors affecting women. Cytogenetic studies of human ECs have not produced very conclusive data, since many of these studies are based on karyotyping of limited number of cases and no really specific karyotypic changes have yet been identified. As the majority of the genes are conserved among mammals, the use of inbred animal model systems may serve as a tool for identification of underlying genes and pathways involved in tumorigenesis in humans. In the present work we used spectral karyotyping (SKY) to identify cancer-related aberrations in a well-characterized experimental model for spontaneous endometrial carcinoma in the BDII rat tumor model.

Results

Analysis of 21 experimental ECs revealed specific nonrandom numerical and structural chromosomal changes. The most recurrent numerical alterations were gains in rat chromosome 4 (RNO4) and losses in RNO15. The most commonly structural changes were mainly in form of chromosomal translocations and were detected in RNO3, RNO6, RNO10, RNO11, RNO12, and RNO20. Unbalanced chromosomal translocations involving RNO3p was the most commonly observed structural changes in this material followed by RNO11p and RNO10 translocations.

Conclusion

The non-random nature of these events, as documented by their high frequencies of incidence, is suggesting for dynamic selection of these changes during experimental EC tumorigenesis and therefore for their potential contribution into development of this malignancy. Comparative molecular analysis of the identified genetic changes in this tumor model with those reported in the human ECs may provide new insights into underlying genetic changes involved in EC development and tumorigenesis.

Loss of glutathione peroxidase 3 expression is correlated with epigenetic mechanisms in endometrial adenocarcinoma

Glutathione peroxidase 3 (GPX3) is one of the key enzymes in the cellular defense against oxidative stress and the hepatocyte growth factor receptor, (MET) has been suggested to be influenced by the GPX3 gene expression. In a previous microarray study performed by our group, Gpx3 was identified as a potential biomarker for rat endometrial adenocarcinoma (EAC), since the expression was highly downregulated in rat EAC tumors. Herein, we have investigated the mRNA expression and Gpx3 and Met in rat EAC by real time quantitative PCR (qPCR), and the methylation status of Gpx3. In addition we have examined the expression of GPX3 and MET in 30 human EACs of different FIGO grades and 20 benign endometrial tissues. We found that the expression of GPX3 was uniformly down regulated in both rat and human EAC, regardless of tumor grade or histopathological subtype, implying that the down-regulation is an early event in EAC. The rate of Gpx3 promoter methylation reaches 91%, where biallelic methylation was present in 90% of the methylated tumors. The expression of the Met oncogene was slightly upregulated in EACs that showed loss of expression of Gpx3, but no tumor suppressor activity of Gpx3/GPX3 was detected. Preliminary results also suggest that the production of H₂O₂ is higher in rat endometrial tumors with down-regulated Gpx3 expression. A likely consequence of loss of GPX3 protein function would be a higher amount of ROS in the cancer cell environment. Thus, the results suggest important clinical implications of the GPX3 expression in EAC, both as a molecular biomarker for EAC and as a potential target for therapeutic interventions.

Gene expression profiling predicts a three-gene expression signature of endometrial adenocarcinoma in a rat model

Background

In the Western world, endometrial cancers are the most common gynaecological neoplastic disorders among women. Initial symptoms are often vague and may be confused with several other conditions or disorders. Thus, there is a need for an easy and reliable diagnostic tool. The objective of this work was to identify a gene expression signature specific for endometrial adenocarcinomas to be used for testing potential endometrial biomarkers.

Results

Changes in expression between endometrial adenocarcinomas and non-/pre-malignant endometrium from the BDII EAC rat model were compared in cDNA microarray assays. By employing classification analysis (Weka) on the expression data from approximately 5600 cDNA clones and TDT analysis on genotype data, we identified a three-gene signature (Gpx3, Bgn and Tgfb3). An independent analysis of differential expression, revealed a total of 354 cDNA clones with significant changes in expression. Among the 10 best ranked clones, Gpx3, Bgn and Tgfb3 were found.

Conclusion

Taken together, we present a unique data set of genes with different expression patterns between EACs and non-/pre-malignant endometrium, and specifically we found three genes that were confirmed in two independent analyses. These three genes are candidates for an EAC signature and further evaluations of their involvement in EAC tumorigenesis will be undertaken.

Rearrangement and allelic imbalance on chromosome 5 leads to homozygous deletions in the CDKN2A/2B tumor suppressor gene region in rat endometrial cancer

The inbred BDII rat is a valuable experimental model for the genetic analysis of hormone-dependent endometrial adenocarcinoma (EAC). One common aberration detected previously by comparative genomic hybridization in rat EAC is loss affecting mostly the middle part of rat chromosome 5 (RNO5). First, we applied an RNO5-specific painting probe and four region-specific gene probes onto tumor cell metaphases from 21 EACs, and found that rearrangements involving RNO5 were common. The copy numbers of loci situated on RNO5 were found to be reduced, particularly for the CDKN2A/2B locus. Second, polymerase chain reaction analysis was performed with 22 genes and markers and homozygous deletions of the CDKN2A exon 1beta and CDKN2B genes were detected in 13 EACs (62%) and of CDKN2A exon 1alpha in 12 EACs (57%) Third, the occurrence of allelic imbalance in RNO5 was analyzed using 39 microsatellite markers covering the entire chromosome and frequent loss of heterozygosity was detected. Even more intriguing was the repeated finding of allele switching in a narrow region of 7 Mb across the CDKN2A/2B locus. We conclude that genetic events affecting the middle part of RNO5 (including bands 5q31 approximately q33 and the CDKN2A locus) contribute to the development of EAC in rat, with the CDKN2A locus having a primary role.

Recurrent chromosome 10 aberrations and Tp53 mutations in rat endometrial adenocarcinomas

Human genetic heterogeneity and differences in the environment and life style make analysis of complex diseases such as cancer difficult. By using inbred animal strains, the genetic variability can be minimized and the environmental factors can be reasonably controlled. Endometrial adenocarcinoma (EAC) is the most common gynecologic malignancy, ranking fourth in incidence among tumors in women. The inbred BDII rat strain is genetically prone to spontaneously develop hormone-related EAC, and can be used as a tool to investigate and characterize genetic changes in this tumor type. In the present project, BDII females were crossed to males from two nonsusceptible rat strains and F1, F2, and backcross progeny were produced. Genetic and molecular genetic analysis of tumors showed that rat chromosome 10 (RNO10) was frequently involved in genetic changes. Our data indicate that often there was loss of chromosomal material in the proximal to middle part of the chromosome followed by gains in distal RNO10. This suggested that there is a tumor suppressor gene(s) in the proximal to middle part of RNO10 and an oncogene(s) in the distal part of the chromosome with potential significance in EAC development. The Tp53 gene, located at band RNO10q24-q25, was a strong candidate target for the observed aberrations affecting the middle part of the chromosome. However, our Tp53 gene mutation analyses suggested that a second gene situated very close to Tp53 might be the main target for the observed pattern of genetic changes.

Analysis of chromosome 10 aberrations in rat endometrial cancer-evidence for a tumor suppressor locus distal to Tp53

We have recently shown in the BDII rat model of human endometrial adenocarcinoma (EAC), rat chromosome 10 (RNO10) is frequently involved in chromosomal aberrations. In the present study, we investigated the association between RNO10 deletions, allelic imbalance (AI) at RNO10q24 and Tp53 mutation in 27 rat EAC tumors. We detected chromosomal breakage accompanied by loss of proximal and/or gain of distal parts of RNO10 in approximately 2/3 of the tumors. This finding is suggestive of a tumor suppressor activity encoded from the proximal RNO10. Given the fact that Tp53 is located at RNO10q24-q25, we then performed Tp53 mutation analysis. However, we could not find a strong correlation between AI/deletions at RNO10q24 and Tp53 mutation. Instead, the observed patterns for AI, chromosomal breaks and deletions suggest that major selection was directed against a region located close to, but distal of Tp53. In different human malignancies a similar situation of AI at chromosome band 17p13.3 (HSA17p13.3) unassociated with TP53 mutation has been observed. Although RNO10 is largely homologous to HSA17, the conservation with respect to gene order among them is not extensive. We utilized publicly available draft DNA sequences to study intrachromosomal rearrangement during the divergence between HSA17 and RNO10. By using reciprocal comparison of rat and human genome data, we could substantially narrow down the candidate tumor suppressor region in rat from 3 Mb to a chromosomal segment of about 0.5 Mb in size. These results provide scientific groundwork for identification of the putative tumor suppressor gene(s) at 17p13.3 in human tumors.

RatMap--rat genome tools and data

The rat genome database RatMap (http://ratmap.org or http://ratmap.gen.gu.se) has been one of the main resources for rat genome information since 1994. The database is maintained by CMB–Genetics at Göteborg University in Sweden and provides information on rat genes, polymorphic rat DNA-markers and rat quantitative trait loci (QTLs), all curated at RatMap. The database is under the supervision of the Rat Gene and Nomenclature Committee (RGNC); thus much attention is paid to rat gene nomenclature. RatMap presents information on rat idiograms, karyotypes and provides a unified presentation of the rat genome sequence and integrated rat linkage maps. A set of tools is also available to facilitate the identification and characterization of rat QTLs, as well as the estimation of exon/intron number and sizes in individual rat genes. Furthermore, comparative gene maps of rat in regard to mouse and human are provided.

Genetic analysis of susceptibility to endometrial adenocarcinoma in the BDII rat model

Most cancers are genetically complex and heterogeneous, a serious obstacle to identifying specific genes underlying the disease. If inbred animal models are used, then both the genetic constitution and environmental influences can be carefully controlled. Females of the BDII inbred rat strain are genetically predisposed to endometrial cancer; more than 90% of virgin BDII females will develop endometrial adenocarcinoma (EAC) during their life span. BDII females were crossed to males from inbred strains with low EAC incidence (SPRD or BN). When F(1) males were backcrossed to BDII females to generate N(1) populations of offspring, about one fourth of the female progeny developed EAC. With transmission disequilibrium test analysis, significant association was detected in three chromosomal regions (on RNO1, RNO11, and RNO17) in the SPRD crosses and in the short arm of RNO20 in the BN crosses. It appears that several susceptibility genes with minor but cooperating effects are responsible for the susceptibility. Furthermore, it seems clear from the interstrain crosses not only that the onset of tumors depends on the presence of susceptibility alleles from the EAC-prone BDII strain, but also that tumor development is affected by the contribution of a genetic component derived from the nonsusceptible strains.

Inconsistent susceptibility to autoimmunity in inbred LEW rats is due to genetic crossbreeding involving segregation of the arthritis-regulating gene Ncf1

We recently identified a single-nucleotide polymorphism in the Ncf1 gene, a component of the NADPH oxidase complex, to be the cause of one of the strongest identified loci for arthritis severity in rats. This polymorphism was found to be naturally occurring in a collection of inbred rat strains as well as in wild rats. Among the inbred strains we found that different LEW substrains (LEW/Ztm and LEW/Mol), originating from different breeders, showed an allelic discrepancy in Ncf1, suggesting an impact on arthritis susceptibility between these substrains. In fact, the LEW/Mol strain was completely resistant to pristane-induced arthritis, in contrast to the LEW/Ztm strain, which was susceptible. Moreover, the LEW/Mol strain had higher production of radical oxygen species in peripheral blood leukocytes, a phenomenon most likely regulated by the polymorphisms in the Ncf1 gene. However, the phenotypic difference between LEW/Mol and LEW/Ztm is most likely a combination of several genes, of which Ncf1 is suggested to be the major regulating gene. This has also been confirmed by previous linkage analyses involving the LEW/Ztm strain which shows that a QTL on chromosome 12, most likely caused by polymorphism of Ncf1, is the major regulatory gene but that other loci are contributing. That more genes are likely to contribute was shown by a complete genome comparison of the LEW/Ztm and the LEW/Mol rat strains that uncovered an introduction of approximately 37% non-LEW genome into the LEW/Mol strain, which probably was caused by past crossbreeding. Therefore, the LEW/Mol should be regarded as a recombinant inbred strain.

Mapping of glutathione transferase (GST) genes in the rat

Glutathione transferases (GST) make up a large group of related enzymes in mammalian tissues. The enzyme molecules are dimeric and at least 13 different subunits occur in the rat. Each subunit appears to be coded for by a distinct gene, and thus there is a large GST gene family in the rat. Recently, there have been several reports of the mapping of rat GST genes. In the present communication we confirm the previous assignments and extend the data with the mapping to rat chromosome 2 of a previously unmapped GST gene (Gstm1), and with the regional mapping of seven Gstp genes. These mappings provide further evidence for conservation of syntenic gene relationships among mammals. The human homologs of Gstm1 map to chromosome 1, and belong to a group of 9 genes that show conserved synteny on rat chromosome 2. The corresponding murine genes in most cases map to mouse chromosome 3. Similarly, the human homolog of Gstp maps to chromosome 11, and is one of 10 genes that exhibit conserved synteny on rat chromosome 1. The corresponding mouse genes map to mouse chromosome 7. Previously only one gene on rat chromosome 8 had a human homolog on chromosome 6, and rat Gsta1 is the second instance. Based on these mappings it appears that a new group of genes will exhibit conserved synteny on rat chromosome 8, human chromosome 6 and mouse chromosome 9. Interestingly, each of the three groups of conserved synteny seems to span the region across the centromeres of the human chromosomes.

Evolutionary aspects of the genomic organization of rat chromosome 10

Using FISH and RH mapping a chromosomal map of rat chromosome 10 (RNO10) was constructed. Our mapping data were complemented by other published data and the final map was compared to maps of mouse and human chromosomes. RNO10 contained segments homologous to mouse chromosomes (MMU) 11, 16 and 17, with evolutionary breakpoints between the three segments situated in the proximal part of RNO10. Near one of these breakpoints (between MMU17 and 11) we found evidence for an inversion ancestral to the mouse that was not ancestral to the condition in the rat. Within each of the chromosome segments identified, the gene order appeared to be largely conserved. This conservation was particularly clear in the long MMU11-homologous segment. RNO10 also contained segments homologous to three human chromosomes (HSA5, 16, 17). However, within each segment of conserved synteny were signs of more extensive rearrangements. At least 13 different evolutionary breakpoints were indicated in the rat-human comparison. In contrast to what was found between rat and mouse, the rat-human evolutionary breaks were distributed along the entire length of RNO10.

Detailed chromosomal and radiation hybrid mapping in the proximal part of rat Chromosome 10 and gene order comparison with mouse and human

The rat provides valuable and sometimes unique models of human complex diseases. To fully exploit the rat models in biomedical research, it is important to have access to detailed knowledge of the rat genome organization as well as its relation to the human genome. Rat Chromosome 10 (RNO10) harbors several important cancer-related genes. Deletions in the proximal part of RNO10 were repeatedly found in a rat model for endometrial cancer. To identify functional and positional candidate genes in the affected region, we used radiation hybrid (RH) mapping and single- and dual-color fluorescence in situ hybridization (FISH) techniques to construct a detailed chromosomal map of the proximal part of RNO10. The regional localization of 14 genes, most of them cancer-related ( Grin2a, Gspt1, Crebbp, Gfer, Tsc2, Tpsb1, Il9r, Il4, Irf1, Csf2, Sparc, Tp53, Thra1, Gh1), and of five microsatellite markers ( D10Mit10, D10Rat42, D10Rat50, D10Rat72, and D10Rat165) was determined on RNO10. For a fifteenth gene, Ppm1b, which had previously been assigned to RNO10, the map position was corrected to RNO6q12-q13.

Blood group A glycosyltransferase occurring as alleles with high sequence difference is transiently induced during a Nippostrongylus brasiliensis parasite infection

Neutral mucin oligosaccharides from the small intestine of control rats and rats infected with the parasite Nippostrongylus brasiliensis were released and analyzed by gas chromatography-mass spectrometry. Infected animals expressed seven blood group A-like structures that were all absent in the control animals. The blood group A nature of these epitopes was confirmed by blood group A reactivity of the prepared mucins, of which Muc2 was one. Transferase assays and Northern blotting on small intestines from infected animals showed that an alpha-N-acetylgalactosaminyltransferase similar to the human blood group A glycosyltransferase had been induced. The expression was a transient event, with a maximum at day 6 of the 13-day-long infection. The rat blood group A glycosyltransferase was cloned, revealing two forms with an amino acid similarity of 95%. Both types had blood group A transferase activity and were probably allelic because none of 12 analyzed inbred strains carried both types. The second type was found in outbred rats and in one inbred strain. First generation offspring of inbred rats of each type were heterozygous, further supporting the allelic hypothesis. The transient induction and the large allelic variation could suggest that glycosyltransferases are part of a dynamic system altering mucins and other glycoconjugates as a protecting mechanism against microbial challenges.

Genetic identification of multiple susceptibility genes involved in the development of endometrial carcinoma in a rat model

There are clear indications that inheritance plays an essential role in certain cases of human endometrial cancer, and there are at least 2 forms of early-onset heritable endometrial adenocarcinomas (EACs). Females of the BDII inbred rat strain are known to be genetically predisposed to endometrial carcinoma, and we have performed a genetic analysis of susceptibility to endometrial cancer in this strain. F(2) populations were generated by crossing BDII females with males from 2 different strains with a low incidence of EAC, and the occurrence of endometrial cancer was studied. Three chromosome regions associated to EAC susceptibility were identified, and the susceptibility genes in these regions were designated Ecs1, Ecs2 and Ecs3. Our results indicate that the genes affecting susceptibility to EAC are different in the 2 crosses, suggesting that the genes behind the susceptibility in BDII animals may interact with different genes in different genetic backgrounds.

High-density marker loss of heterozygosity analysis of rat chromosome 10 in endometrial adenocarcinoma

Endometrial cancer is a disease with serious impact on the human population, but not much is known about genetic factors involved in this complex disease. Female BDII rats are genetically predisposed to spontaneous endometrial carcinoma, and the BDII inbred strain provides an experimental animal model for endometrial carcinoma development. In the present study, BDII females were crossed with males from two nonsusceptible inbred rat strains. Endometrial adenocarcinomas (EACs) developed in a proportion of the F1 and F2 progeny. We screened 18 EAC solid tumors and 9 EAC cell cultures for loss of heterozygosity (LOH) using fluorescent-PCR-based marker allelotyping methodology with 47 microsatellite markers covering the proximal part of rat chromosome 10 (RNO10). Conclusive evidence was obtained for LOH/deletion involving about 56 cM in the proximal part of RNO10 in DNA from six out of seven informative tumor cell cultures. Analysis of the solid tumors confirmed the presence of LOH in this part of RNO10 in 14 of 17 informative tumors. However, from the studies in the solid tumors it appeared that in fact three separate segments in the proximal part of RNO10 were affected. These three LOH/deletion regions were located approximately in cytogenetic bands 10q11-12, 10q22, and 10q24.

Analysis of genetic changes in rat endometrial carcinomas by means of comparative genomic hybridization

Animals of the BDII inbred rat strain are known to be genetically predisposed to endometrial adenocarcinoma (EAC). Using them as models of human EACs, we studied tumors arising in F1 and F2 progeny from BDII animals crossed to animals from two other inbred strains, in which EACs were quite rare. In order to identify chromosomal regions exhibiting DNA copy number changes, comparative genomic hybridization (CGH) was applied in a series corresponding to 27 different solid tumors, most of which were classified as EACs, from these animals. The main findings from the study were that, although many different chromosomes were involved in copy number variation, some of the changes detected were recurrent and quite specific. Among specific changes found were gains in rat chromosome (RNO) regions 4q12 approximately q22, 6q14 approximately q16, and whole chromosome arms in some of the small metacentric chromosomes (e.g., RNO14, 16, and 18). RNO10 was involved in gain in the terminal and proximal regions. Each of these regions contains previously identified cancer-related genes representing possible candidates to be involved in the development of EAC. Furthermore, it was observed that there were clear differences in the pattern of copy number changes between tumors occurring in the two different crosses, and also between solid tumors and cell cultures. Endometrial cancer is the most common human gynecological cancer, but not much is known about specific genetic changes influencing this disease. Two genetic alterations that have been reported from human endometrial cancer are amplification of the ERBB2 gene and mutations in the 12 codon of the KRAS gene. One case of Erbb2 amplification was found but there were no Kras mutations in the rat material studied. We conclude that molecular genetic analysis of the rat BDII model will provide important new information about EAC in mammals.

Genomewide assessment of genetic alterations in DMBA-induced rat sarcomas: cytogenetic, CGH, and allelotype analyses reveal recurrent DNA copy number changes in rat chromosomes 1, 2, 4, and 7

Rat sarcomas, induced by subcutaneous injections of 7,12-dimethylbenz[a]anthracene (DMBA), were studied with the objective of identifying critical chromosome regions associated with tumorigenesis. We employed three genomewide screening techniques-cytogenetics, CGH, and allelotyping-in 19 DMBA-induced sarcomas in F1 (BN/Han x LE/Mol) rats. The most conspicuous finding in the cytogenetic analysis was a high incidence of trisomy for rat chromosome 2 (RNO2). Signs of gene amplification (hsr) were also seen in several tumors. The CGH analysis revealed that gains in copy number were much more common than losses. The gains mostly affected RNO2 (10/19), RNO12q (7/19), and RNO19q (5/19), as well as the proximal part of RNO4 (8/19) and the distal part of RNO7 (7/19). Reduction in copy number was seen in RNO17 (2/19). For the allelotyping, we used 318 polymorphic microsatellite marker loci covering the entire genome. We identified regions of allelic imbalance affecting most of the rat chromosomes. The highest incidences of recurrent allelic imbalance were observed at loci in certain regions in RNO1, 2, 4, and 7 and at lower incidences in parts of RNO12, 16, 18, and 19. The combined results suggested that genetic alterations detected in RNO2 and RNO12 usually corresponded to complete or partial trisomy, whereas those in RNO1 and RNO7 seemed to involve regional deletions and/or gains. Furthermore, we could confirm that copy number gains occur proximally in RNO4, where a previous study showed amplification of the Met oncogene in a subset of these tumors.

A dual-color FISH framework map for the characterization of the Sai1 tumor suppression region on rat chromosome 5

The analysis of cell hybrids between malignant mouse hepatoma cells and normal rat fibroblasts has previously demonstrated the critical role of a deletion in rat chromosome 5 (RNO5) that was related to an anchorage independent phenotype. Those hybrids that were anchorage independent displayed loss of the entire RNO5 or an interstitial deletion in RNO5. These findings suggested that a putative tumor suppressor gene, Sai1 (suppression of anchorage independence 1), was located within the deleted region. To explore the molecular basis of the tumor suppressor activity of the Sai1 region, we analyzed the RNO5q23-q36 region with several genes and microsatellite markers that could be assigned to the region, as well as with new markers derived by representational difference analysis (RDA) or by microdissection. Dual-color FISH was used to construct a detailed physical map of the entire RNO5. These new data can be used to connect the physical and linkage maps in the rat, as well as to identify the details of the comparative map with other mammalian species including humans and mice. Using as FISH reagents genomic YAC, P1, or phage lambda clones corresponding to RNO5 markers, the order and unique positions of 18 markers could be established. The map provided a framework for the detailed characterization of the deletion found in anchorage independent hybrids. All markers within the bands RNO5q31.3-q35 were shown to be lost, including known cancer-related genes such as Ifna (5q32), Cdkn2a, -b (5q32), Jun (5q34), and Cdkn2c (5q35). However, the aberration in the deletion chromosome turned out to be more complex than originally thought in that we detected the presence of a paracentric inversion in addition to a deletion. The inversion led to the juxtaposition of the gene markers Tal2 (5q24.1) and Cd30lg (5q24.3). The framework map will provide the basis for the detailed physical YAC clone contig mapping of this region, and facilitate the identification and characterization of the Sai1 locus.

Amplification and overexpression of the hepatocyte growth factor receptor (HGFR/MET) in rat DMBA sarcomas

In the present study subcutaneous fibrosarcomas were induced by the carcinogen 7,12-dimethylbenz(a)anthracene (DMBA) in rats from F1 generation cross breedings of two different inbred strains. Comparative genomic hybridization (CGH) analysis, which allows detection of DNA sequence copy changes, was applied to one of the tumors and it was found that there were increased copy numbers of sequences at chromosome 4q12-q21 in this tumor. We have previously determined that the loci for the hepatocyte growth factor (Hgf) and hepatocyte growth factor receptor (Hgfr/Met), a protooncogene, are situated in this particular chromosome region. Using probes for the two genes in FISH (fluorescence in situ hybridization) and in Southern blots we found that the Hgfr/Met gene was amplified in five of the 19 sarcomas studied, and that the Hgf gene was coamplified in two of them. Northern and Western blots and tyrosine phosphorylation analysis showed that the HGF receptor was overexpressed and functional in all five tumors, as well as in two additional tumors. In summary, both amplification and overexpression of the Hgfr/Met gene was found in about 25% of DMBA-induced experimental rat sarcomas, and HGF receptor overexpression alone was seen in two additional tumors. Possibly this reflects an involvement in paracrine or autocrine stimulation of growth and invasiveness by HGF. Our finding could provide a rodent model system to increased knowledge about causality and therapy, which may be applicable to the sizeable fraction of human musculoskeletal tumors displaying MET overexpression.

Rat Chromosome 2: assignment of the genes encoding cyclin B1, interleukin 6 signal transducer, and proprotein convertase 1 to the Mcs1-containing region and identification of new microsatellite markers

The rat Chromosome (Chr) 2 harbors several genes controlling tumor growth or development, blood pressure, and non-insulin-dependent diabetes mellitus. We report that the region (2q1) containing the mammary susceptibility cancer gene Mcs1 also harbors the genes encoding cyclin B1, interleukin 6 signal transducer (gp130), and proprotein convertase 1. We also generated 13 new anonymous microsatellite markers from Chr 2-sorted DNA. These markers, as well as a microsatellite marker in the cyclin B1 gene, were genetically mapped in combination with known markers. A cyclin B1-related gene was also cytogenetically assigned to rat Chr 11q22-q23.

Gene-based anchoring of the rat genetic linkage and cytogenetic maps: new regional localizations, orientation of the linkage groups, and insights into mammalian chromosome evolution

In order to generate anchor points connecting the rat cytogenetic and genetic maps, the cytogenetic position of 62 rat markers (including 55 genes) already localized genetically was determined by fluorescence in situ hybridization. Whenever possible, markers located near one end of the linkage groups were included. These new localizations allowed us to unambiguously orient the 20 autosomal and the X chromosome linkage groups. The position of the centromere in the linkage map could also be determined in the case of several metacentric chromosomes. In addition, the regional localization of 15 other rat genes was determined. These new data bring useful information with respect to comparative mapping with the mouse and the human and to mammalian evolution. They illustrate, for instance, that groups of genes can remain syntenic during mammalian evolution while being subjected to intrachromosomal rearrangements in some lineages (synteny is conserved while gene order is not). This analysis also disclosed cases of synteny conservation in one the two rodent species and the human, while the synteny is split in the other rodent species: such configurations are likely examples of lineage-specific interchromosomal rearrangements associated with speciation.

Rat chromosome 1: regional localization of seven genes (Slc9a3, Srd5a1, Esr, Tcp1, Grik5, Tnnt3, Jak2) and anchoring of the genetic linkage map to the cytogenetic map

Seven genes were regionally localized on rat Chromosome (Chr) 1, from 1p11 to 1q42, and two of these genes were also included in a linkage map. This mapping work integrates the genetic linkage map and the cytogenetic map, and allows us to orient the linkage map with respect to the centromere, and to deduce the approximate position of the centromere in the linkage map. These mapping data also indicate that the Slc9a3 gene, encoding the Na+/H+ exchanger 3, is an unlikely candidate for the blood pressure loci assigned to rat Chr 1. These new localizations expand comparative mapping between rat Chr 1 and mouse or human chromosomes.

Isolation of DNA markers for the rat Sai 1 gene for suppression of anchorage independence by using representational difference analysis

We have applied the representational difference analysis (RDA) to isolate genetic markers for a deletion on the rat chromosome RNO5q22-33. This deletion occurred in anchorage independent sublines of a normal rat fibroblast x mouse hepatoma cell hybrid (BS181) (Islam 1989). Normal rat tissue DNA provided the "tester" and the BS181 hybrid DNA the "driver" in the RDA hybridization/selection reactions. Out of twelve RDA derived DNA sequences that were analyzed in detail using a rat X mouse cell hybrid panel for chromosome mapping, nine (75%) were found to represent RNO5 deletions, whereas the other three were new RFLPs mapping to other chromosomes. In two cases, the RDA sequences were also analyzed by fluorescence in situ hybridization (FISH) and found to give distinct signals in the RNOq22-33 region. This result emphasizes teh significance of the previous cytogenetic analysis of this hybrid, which indicated the presence of a gene for the suppression of anchorage independence, Sai 1, in this deletion region. The RDA derived sequences isolated by this work will provide a valuable source of new genetic markers for the further detailed analysis of the Sai 1 deletion region.

Chromosome mapping of rat histone genes H1fv, H1d, H1t, Th2a and Th2b

Chromosome assignment of the rat histone genes H1t, H1d (H1.4), H1fv (H10), Th2a and Th2b is described. The testicularly expressed histone genes H1t, Th2a and Th2b could be assigned to rat chromosome (RNO) 17 by PCR analysis of somatic cell hybrid DNAs. The H1d gene was mapped to RNO17p12-->p11 by FISH. These genes might form a histone gene cluster homologous to that found on HSA6p21.3 in humans and MMU13A2-3 in mice. The rat histone H1fv gene was assigned to RNO7 by PCR. This result allows the inclusion of rat H1fv to an established conserved group of syntenic genes in rat, mouse and human on chromosomes RNO7, MMU15 and HSA22, respectively.

Microsatellite polymorphism analysis allows the individual assignment of the rat 11 beta-hydroxylase gene (Cyp11b1) and the rat aldosterone synthase gene (Cyp11b2) to chromosome 7 using rat x mouse somatic cell hybrids and identifies differences between and within various rat strains

Mouse hepatoma x rat hepatocyte hybrids that segregate rat chromosomes were used to determine the chromosomal location of the rat genes encoding 11 beta-hydroxylase and aldosterone synthase (Cyp11b1 and Cyp11b2 respectively). By means of species-specific restriction fragments and microsatellite markers both genes were mapped to rat chromosome 7. The Cyp11b1 microsatellite marker was subsequently found to vary in length between and within rat strains. Furthermore, we compared the sequences of Cyp11b1 markers in two genetically hypertensive strains of rat with their normotensive counterparts. Previous studies have indicated that 11 beta-hydroxylase activities in Milan and Lyon hypertensive strains are different from their respective genetic controls. The Cyp11b1 microsatellite regions from Lyon hypotensive and normotensive strains of rat were similar and were both shorter by 15 bases than that of the Lyon hypertensive strain. The Cyp11b1 marker in Milan hypertensive (MHS) and normotensive (MNS) strains differ from all the Lyon strains and from each other. The MHS marker is 12 bases shorter than that of MNS rats. These differences in microsatellite length may provide useful polymorphic markers in cosegregation studies of genetic hypertension in rats.

Molecular cloning of a cDNA coding for a region of an apoprotein from the 'insoluble' mucin complex of rat small intestine

The major part of rat small intestinal mucins occurs as an 'insoluble' glycoprotein complex unextractable in 6 M guanidinium chloride unless disulfide bonds are cleaved. One of the trypsin-resistant high-glycosylated domains of this complex (glycopeptide A) was recently isolated. We have now deglycosylated it with HF, injected it into rabbits and the obtained antiserum was used for expression cloning providing a cDNA clone (VR-1A). This clone contained an open reading frame of 235 amino acids composed of two regions. The deduced N-terminal 53 amino acids included seven Cys residues and only one Ser, followed by a region of 182 residues with 64% Ser and Thr but devoid of Cys residues. Analysis of mRNA revealed a transcript of about 12 kb, identical in size to a band labelled with a probe based on the rat mucin-like protein (MLP/Muc2) cDNA. Pulsed-field gel electrophoresis of genomic rat DNA showed identical bands (380 and 500 kb) when blots were sequentially probed with the MLP/Muc2 probe and VR-1A. A panel of mouse x rat hybrids was used to localize the gene corresponding to both VR-1A and Muc2 to rat chromosome 1. The results strongly suggest that the 'insoluble' mucin complex of the rat small intestine is encoded by the Muc2 gene.

The retinoblastoma-related gene, RB2, maps to human chromosome 16q12 and rat chromosome 19

A retinoblastoma-related human gene, referred to as RB2, has been cloned based on sequence homology of the E1A-binding domain of the retinoblastoma gene. Structural homology with the retinoblastoma gene suggests a possible function of RB2 as a tumor suppressor gene. In this study, we have mapped this gene to human chromosome 16q12.2 and rat chromosome 19, using fluorescence in situ hybridization and somatic hybrid cell analysis, respectively. Based on known syntenic relationships among human, rat and mouse, the data suggest that the mouse homolog resides on chromosome 8. Deletions of chromosome 16q have been found in several human neoplasias (including breast, ovarian, hepatic, and prostatic cancers) which is in support of an involvement of RB2 in human cancer as a tumor suppressor gene.