ABERRATIONS OF A PUTATIVE TUMOR SUPPRESSOR GENE SEL1L IN PANCREATIC DUCTAL ADENOCARCINOMA

Introduction: Pancreatic cancer is the fourth leading cause of cancer-related death among males and females in the United States. Sel-1-like (SEL1L) is a putative tumor suppressor gene that is downregulated in a significant proportion of human pancreatic ductal adenocarcinoma (PDAC). It was hypothesized that SEL1L expression could be down-modulated by somatic mutation, loss of heterozygosity (LOH), CpG island hypermethylation and/or aberrantly expressed microRNAs (miRNAs). Material and methods: In 42 PDAC tumors, the SEL1L coding region was amplified using reverse transcription polymerase chain reaction (RT-PCR), and analyzed by agarose gel electrophoresis and sequenced to search for mutations. Using fluorescent fragment analysis, two intragenic microsatellites in the SEL1L gene region were examined to detect LOH in a total of 73 pairs of PDAC tumors and normal-appearing adjacent tissues. Bisulfite DNA sequencing was performed to determine the methylation status of the SEL1L promoter in 41 PDAC tumors and 6 PDAC cell lines. Using real-time quantitative PCR, the expression levels of SEL1L mRNA and 7 aberrantly upregulated miRNAs that potentially target SEL1L were assessed in 42 PDAC tumor and normal pairs. Statistical methods were applied to evaluate the correlation between SEL1L mRNA and the miRNAs. Further the interaction was determined by functional analysis using a molecular biological approach. Results: No mutations were detected in the SEL1L coding region. More than 50% of the samples displayed abnormally alternate or aberrant spliced transcripts of SEL1L. About 14.5% of the tumors displayed LOH at the CAR/CAL microsatellite locus and 10.7% at the RepIN20 microsatellite locus. However, the presence of LOH did not show significant association with SEL1L downregulation. No methylation was observed in the SEL1L promoter. Statistical analysis showed that SEL1L mRNA expression levels significantly and inversely correlated with the expression of hsa-mir-143, hsa-mir-155, and hsa-mir-223. Functional analysis indicated that hsa-mir-155 acted as a suppressor of SEL1L in PL18 and MDAPanc3 PDAC cell lines. Discussion: Evidence from these studies suggested that SEL1L was possibly downregulated by aberrantly upregulated miRNAs in PDAC. Future studies should be directed towards developing a better understanding of the mechanisms for generation of aberrant SEL1L transcripts, and further analysis of miRNAs that may downregulate SEL1L.

Molecular events associated with trisomy of chromosome 7 in mouse skin chemical carcinogenesis

The primary objective of this study has been to investigate the effects at the molecular level of trisomy of mouse chromosome 7 in chemically induced skin tumors. It was previously proposed that the initiation event in the mouse skin carcinogenesis model is a heterozygous mutation of the Ha-ras-1 gene, mapped to chromosome 7. Previous studies in this laboratory identified trisomy 7 as one of the primary nonrandom cytogenetic abnormalities found in the majority of severely dysplastic papillomas and squamous cell carcinomas induced in SENCAR mice by an initiation-promotion protocol. Therefore, the first hypothesis tested was that trisomy 7 occurs by specific duplication of the chromosome carrying a mutated Ha-ras-1 allele. Results of a quantitative analysis of normal/mutated allelic ratios of the Ha-ras-1 gene confirmed this hypothesis, showing that most of the tumors exhibited overrepresentation of the mutated allele in the form of 1/2, 0/3, and 0/2 (normal/mutated) ratios. In addition, histopathological analysis of the tumors showed an apparent association between the degree of malignancy and the dosage of the mutated Ha-ras-1 allele. To determine the mechanism for loss of the normal Ha-ras-1 allele, found in 30% of the tumors, a comparison of constitutional and tumor genotypes was performed at different informative loci of chromosome 7. By combining Southern blot and polymerase chain reaction fragment length polymorphism analyses of DNAs extracted from squamous cell carcinomas, complete loss of heterozygosity was detected in 15 of 20 tumors at the Hbb locus, and in 5 of 5 tumors at the int-2 locus, both distal to Ha-ras-1. In addition, polymerase chain reaction analysis of DNA extracted from papillomas indicated that loss of heterozygosity occurs in late-stage lesions exhibiting a high degree of dysplasia and areas of microinvasion, suggesting that this event may be associated to the acquisition of the malignant phenotype. Allelic dosage analysis of tumors that had become homozygous at Hbb but retained heterozygosis at Ha-ras-1, indicated that loss of heterozygosity on mouse chromosome 7 occurs by a mitotic recombination mechanism. Overall, these findings suggest the presence of a putative tumor suppressor locus on the 7F1-ter region of mouse chromosome 7. Thus, loss of function by homozygosis at this putative suppressor locus may complement activation of the Ha-ras-1 gene during tumor progression, and might be associated with the malignant conversion stage of mouse skin carcinogenesis. ^

Mitotic mapping of a novel tumor suppressor locus on human chromosome 3q important in osteosarcoma tumorigenesis

Tumor-specific loss of constitutional heterozygosity by deletion, mitotic recombination or nondisjunction is a common mechanism for tumor suppressor allele inactivation. When loss of heterozygosity is the result of mitotic recombination, or a segmental deletion event, only a portion of the chromosome is lost. This information can be used to map the location of new tumor suppressor genes. In osteosarcoma, the highest frequencies of loss of heterozygosity have been reported for chromosomes 3q, 13q, 17p. On chromosomes 13q and 17p, allelic losses are associated with loss of function at the retinoblastoma susceptibility locus (RB1) and the p53 locus, respectively. Chromosome 3q is also of particular interest because the high percent of loss of heterozygosity (62%-75%) suggests the presence of another tumor suppressor important for osteosarcoma tumorigenesis. To localize this putative tumor suppressor gene, we used polymorphic markers on chromosome 3q to find the smallest common region of allele loss. This putative tumor suppressor was localized to a 700 kb region on chromosome 3q26.2 between the polymorphic loci D3S1282 and D3S1246. ^

Characterization and subchromosomal location of a novel tumor suppressor gene on human chromosome 7

Alterations in oncogenes and tumor suppressor genes (TSGs) are considered to be critical steps in oncogenesis. Consistent deletions and loss of heterozygosity (LOH) of polymorphic markers in a determinate chromosomal fragment are known to be indicative of a closely mapping TSG. Deletion of the long arm of chromosome 7 (hchr 7) is a frequent trait in many kinds of human primary tumors. LOH was studied with an extensive set of markers on chromosome 7q in several types of human neoplasias (primary breast, prostate, colon, ovarian and head and neck carcinomas) to determine the location of a putative TSG. The extent of LOH varied depending the type of tumor studied but all the LOH curves we obtained had a peak at (C-A)$\sb{\rm n}$ microsatellite repeat D7S522 at 7q31.1 and showed a Gaussian distribution. The high incidence of LOH in all tumor types studied suggests that a TSG relevant to the development of epithelial cancers is present on the 7q31.1. To investigate whether the putative TSG is conserved in the syntenic mouse locus, we studied LOH of 30 markers along mouse chromosome 6 (mchr 6) in chemically induced squamous cell carcinomas (SCCs). Tumors were obtained from SENCAR and C57BL/6 x DBA/2 F1 females by a two-stage carcinogenesis protocol. The high incidence of LOH in the tumor types studied suggests that a TSG relevant to the development of epithelial cancers is present on mchr 6 A1. Since this segment is syntenic with the hchr 7q31, these data indicate that the putative TSG is conserved in both species. Functional evidence for the existence of a TSG in hchr 7 was obtained by microcell fusion transfer of a single hchr 7 into a murine SCC-derived cell line. Five out of seven hybrids had two to three-fold longer latency periods for in vivo tumorigenicity assays than parental cells. One of the unrepressed hybrids had a deletion in the introduced chromosome 7 involving q31.1-q31.3, confirming the LOH data. ^

Evidence for the involvement of proximal chromosome 3P loci in the progression of Von Hippel-Lindau related tumors

Von Hippel-Lindau (VHL) disease is an autosomal dominant disorder characterized by the development of retinal and central nervous system hemangioblastoma, renal cell carcinoma (RCC), pheochromocytoma and pancreatic islet cell tumors (PICT). The VHL gene maps to chromosome 3p25 and has been shown to be mutated in 57% of sporadic cases of RCC, implicating VHL in the genesis of RCC. We report a multigeneration VHL kindred in which four affected female siblings developed PICT at early ages. Analysis of the three coding exons of the VHL gene in this family revealed a single, missense mutation in codon 238. Inheritance of the 238 mutation has been reported to correlate with a 62% risk of pheochromocytoma development. In this kindred, all affected individuals carried the mutation as well as one additional sibling who showed no evidence of disease. Clinical screening of this individual indicated small ($<$1 cm) pancreatic and kidney tumors. Results suggest that inheritance of the codon 238 mutation does not correlate with early onset pheochromocytoma. Rather, the only individual in the pedigree with pheochromocytoma was the proband's mother who developed bilateral pheochromocytoma at the age of 62. Thus, the VHL codon 238 mutation may predispose to late onset pheochromocytoma in this family; however, it does not explain the preponderance of PICT in the third generation since this mutation has not been reported to increase the risk of developing pancreatic lesions. This suggests that inheritance of the codon 238 mutation and subsequent somatic inactivation of the wild type allele of the VHL gene may not be sufficient to explain the initiation and subsequent progression to malignancy in VHL-associated neoplasms. Since the two tumor types that most frequently progress to malignancy are RCC and PICT, we asked whether loss of heterozygosity (LOH) could be detected proximal to the VHL gene on chromosome 3 in distinct regions of 3p previously implicated by LOH and cytogenetic studies to contain tumor suppressor loci for RCC. LOH was performed on high molecular weight DNA isolated from peripheral blood and frozen tumor tissue of family members using microsatellite markers spanning 3p. Results indicated LOH for all informative 3p loci in tumor tissue from affected individuals with PICT. LOH was detected along the entire length of the chromosome arm and included the proximal region of 3p13-14.2 implicated in the hereditary form of renal cell carcinoma.^ If 3p LOH were a critical event in pancreatic islet cell tumorigenesis, then it should be expected that LOH in sporadic islet cell tumors would also be observed. We expanded LOH studies to include sporadic cases of PICT. Consistent LOH was observed on 3p with a highest frequency LOH in the region 3p21.2. This is the first evidence for an association between chromosome 3 loci and pancreatic islet cell tumorigenesis. (Abstract shortened by UMI.) ^

Cost and power analysis of linkage disequilibrium methods for locating genes influencing quantitative traits in humans

Linkage disequilibrium methods can be used to find genes influencing quantitative trait variation in humans. Linkage disequilibrium methods can require smaller sample sizes than linkage equilibrium methods, such as the variance component approach to find loci with a specific effect size. The increase in power is at the expense of requiring more markers to be typed to scan the entire genome. This thesis compares different linkage disequilibrium methods to determine which factors influence the power to detect disequilibrium. The costs of disequilibrium and equilibrium tests were compared to determine whether the savings in phenotyping costs when using disequilibrium methods outweigh the additional genotyping costs.^ Nine linkage disequilibrium tests were examined by simulation. Five tests involve selecting isolated unrelated individuals while four involved the selection of parent child trios (TDT). All nine tests were found to be able to identify disequilibrium with the correct significance level in Hardy-Weinberg populations. Increasing linked genetic variance and trait allele frequency were found to increase the power to detect disequilibrium, while increasing the number of generations and distance between marker and trait loci decreased the power to detect disequilibrium. Discordant sampling was used for several of the tests. It was found that the more stringent the sampling, the greater the power to detect disequilibrium in a sample of given size. The power to detect disequilibrium was not affected by the presence of polygenic effects.^ When the trait locus had more than two trait alleles, the power of the tests maximized to less than one. For the simulation methods used here, when there were more than two-trait alleles there was a probability equal to 1-heterozygosity of the marker locus that both trait alleles were in disequilibrium with the same marker allele, resulting in the marker being uninformative for disequilibrium.^ The five tests using isolated unrelated individuals were found to have excess error rates when there was disequilibrium due to population admixture. Increased error rates also resulted from increased unlinked major gene effects, discordant trait allele frequency, and increased disequilibrium. Polygenic effects did not affect the error rates. The TDT, Transmission Disequilibrium Test, based tests were not liable to any increase in error rates.^ For all sample ascertainment costs, for recent mutations ($<$100 generations) linkage disequilibrium tests were less expensive than the variance component test to carry out. Candidate gene scans saved even more money. The use of recently admixed populations also decreased the cost of performing a linkage disequilibrium test. ^

GENETIC EVENTS PREDISPOSING TO WILMS' TUMOR

Wilms' tumor (WT) is a childhood embryonic tumor of the kidney. In some cases, WT has been associated with a chromosome deletion in the region 11p13. The majority of WT cases, however, have normal karyotypes with no discernable deletions or rearrangements of chromosome 11.^ To study the genetic events predisposing to the development of WT, I have used a number of gene markers specific for chromosome 11. Gene probes for human catalase and apolipoprotein A1 were localized to chromosome 11 by in situ hybridization. A number of other probes previously mapped to chromosome 11 were also used. Nine WT patients who were heterozygous for at least one 11p marker were shown to lose heterozygosity in their tumor DNA. Gene dosage experiments demonstrated that two chromosomes 11 were present although loss of heterozygosity had occurred in all but two cases. By using gene probes from the short and long arms of chromosome 11, I discerned that loss of heterozygosity was due to somatic recombination in four cases, chromosome deletion in two cases, and chromosome loss and reduplication or somatic recombination in these cases. Examination of DNAs from the parents of six of these patients indicated that the alleles that were lost in tumor tissues were alleles inherited from the mother. In sporadic WT cases one would expect the loss of alleles to be random. These data suggest that the loss of alleles resulting in the development of WT is not a random event, however, the significance of this is not known. ^

Heterozygous GHR gene mutation in a child with idiopathic short stature.

Abstract Several monogenic defects have been reported to be associated with idiopathic short stature. Focusing on growth hormone receptor (GHR)-gene alterations, the heterozygosity of the same gene defect may be associated with a range of growth deficits. We found a heterozygous mutation (V144I) within exon 6 of the GHR gene in a patient with a low level of insulin-like growth factor I (IGF-I), normal level of GH, and severe short stature. Despite the lack of statistical difference, an overall tendency for reduced wt-GH-induction of GHR activation and Jak/Stat signalling in cells transiently expressing GHR-V144I alone or co-expressing wt-GHR compared to cells expressing only wt-GHR was found when GH doses were increased. Our results suggest that, although GHR sequence variants are responsible for some functional alterations commonly observed in children with idiopathic short stature, these changes may not explain all the height deficits observed in these subjects.

New generation genomic platforms in investigation of complex diseases and BEN.

(Full text is available at http://www.manu.edu.mk/prilozi). New generation genomic platforms enable us to decipher the complex genetic basis of complex diseases and Balkan Endemic Nephropathy (BEN) at a high-throughput basis. They give valuable information about predisposing Single Nucleotide Polymorphisms (SNPs), Copy Number Variations (CNVs) or Loss of Heterozygosity (LOH) (using SNP-array) and about disease-causing mutations along the whole sequence of candidate-genes (using Next Generation Sequencing). This information could be used for screening of individuals in risk families and moving the main medicine stream to the prevention. They also might have an impact on more effective treatment. Here we discuss these genomic platforms and report some applications of SNP-array technology in a case with familial nephrotic syndrome. Key words: complex diseases, genome wide association studies, SNP, genomic arrays, next generation sequ-encing.

Incomplete Distal Renal Tubular Acidosis from a Heterozygous Mutation of the V-ATPase B1 Subunit

Congenital distal renal tubular acidosis (dRTA) from mutations of the B1 subunit of the V-ATPase is considered an autosomal recessive disease. We analyzed a dRTA kindred with a truncation-mutation of B1 (p.Phe468fsX487) previously shown to have failure of assembly into the V1 domain of the V-ATPase. All heterozygous carriers in this kindred have normal plasma bicarbonate concentrations, thus evaded the diagnosis of RTA. However, inappropriately high urine pH, hypocitraturia, and hypercalciuria are present either individually or in combination in the heterozygotes at baseline. Two of the heterozygotes studied also have inappropriate urinary acidification with acute ammonium chloride loading and impaired urine-blood pCO2 gradient during bicarbonaturia indicating presence of H+ gradient and flux defects. In normal human renal papillae, wild type B1 is located primarily on the plasma membrane but papilla from one of the heterozygote who had kidney stones had renal tissue secured from surgery showed B1 in both plasma membrane as well as a diffuse intracellular staining. Titrating increasing amounts of the mutant B1 subunit did not exhibit negative dominance over the expression, cellular distribution, or H+-pump activity of the wild type B1 in mammalian HEK293 cells and in V-ATPase-deficient S. cerevisiae. This is the first demonstration of renal acidification defects and nephrolithiasis in heterozygous carriers of mutant B1 subunit; which cannot be attributable to negative dominance. We propose that heterozygosity may lead to mild real acidification defects due to haploinsufficiency. B1 heterozygosity should be considered in patients with calcium nephrolithiasis and urinary abnormalities such as alkalinuria or hypocitraturia.

An ARHGEF10 deletion is highly associated with a juvenile-onset inherited polyneuropathy in Leonberger and Saint Bernard dogs

An inherited polyneuropathy (PN) observed in Leonberger dogs has clinical similarities to a genetically heterogeneous group of peripheral neuropathies termed Charcot-Marie-Tooth (CMT) disease in humans. The Leonberger disorder is a severe, juvenile-onset, chronic, progressive, and mixed PN, characterized by exercise intolerance, gait abnormalities and muscle atrophy of the pelvic limbs, as well as inspiratory stridor and dyspnea. We mapped a PN locus in Leonbergers to a 250 kb region on canine chromosome 16 (Praw = 1.16×10-10, Pgenome, corrected = 0.006) utilizing a high-density SNP array. Within this interval is the ARHGEF10 gene, a member of the rho family of GTPases known to be involved in neuronal growth and axonal migration, and implicated in human hypomyelination. ARHGEF10 sequencing identified a 10 bp deletion in affected dogs that removes four nucleotides from the 3'-end of exon 17 and six nucleotides from the 5'-end of intron 17 (c.1955_1958+6delCACGGTGAGC). This eliminates the 3'-splice junction of exon 17, creates an alternate splice site immediately downstream in which the processed mRNA contains a frame shift, and generates a premature stop codon predicted to truncate approximately 50% of the protein. Homozygosity for the deletion was highly associated with the severe juvenile-onset PN phenotype in both Leonberger and Saint Bernard dogs. The overall clinical picture of PN in these breeds, and the effects of sex and heterozygosity of the ARHGEF10 deletion, are less clear due to the likely presence of other forms of PN with variable ages of onset and severity of clinical signs. This is the first documented severe polyneuropathy associated with a mutation in ARHGEF10 in any species.

Echinococcus P29 antigen: molecular characterization and implication on post-surgery follow-up of CE patients infected with different species of the Echinococcus granulosus complex.

The protein P29 is a potential serological marker for post-treatment monitoring of cystic echinococcosis (CE) especially in young patients. We now have demonstrated that P29 is encoded in the Echinococcus genus by a single gene consisting of 7 exons spanning 1.2 kb of DNA. Variability of the p29 gene at inter- and intra-species level was assessed with 50 cDNA and 280 genomic DNA clones isolated from different E. granulosus s.l. isolates (E. granulosus sensu stricto (G1), E. equinus (G4), E. ortleppi (G5), E. canadensis (G6), E. canadensis (G7) and E. canadensis (G10)) as well as four E. multilocularis isolates. Scarce interspecies polymorphism at the p29 locus was observed and affected predominantly E. granulosus s.s. (G1), where we identified two alleles (A1 and A2) coding for identical P29 proteins and yielding in three genotypes (A1/A1, A2/A2 and A1/A2). Genotypic frequencies expected under Hardy-Weinberg equilibrium revealed a high rate of heterozygosity (47%) that strongly supports the hypothesis that E. granulosus s.s. (G1) is predominantly outbreeding. Comparative sequence analyses of the complete p29 gene showed that phylogenetic relationships within the genus Echinococcus were in agreement with those of previous nuclear gene studies. At the protein level, the deduced P29 amino acid (AA) sequences exhibited a high level of conservation, ranging from 97.9% AA sequence identity among the whole E. granulosus s.l. group to 99.58% identity among E. multilocularis isolates. We showed that P29 proteins of these two species differ by three AA substitutions without implication for antigenicity. In Western-blot analyses, serum antibodies from a human CE patient infected with E. canadensis (G6) strongly reacted with recombinant P29 from E. granulosus s.s. (G1) (recEg(G1)P29). In the same line, human anti-Eg(G1)P29 antibodies bound to recEcnd(G6)P29. Thus, minor AA sequence variations appear not to impair the prognostic serological use of P29.

Polymorphous oligodendroglioma of Zülch revisited: a genetically heterogeneous group of anaplastic gliomas including tumors of bona fide oligodendroglial differentiation.

A polymorphous variant of oligodendroglioma was described by K.J. Zülch half a century ago, and is only very sporadically referred to in the subsequent literature. In particular, no comprehensive analysis with respect to clinical or genetic features of these tumors is available. From a current perspective, the term polymorphous oligodendroglioma (pO) may appear as contradictory in terms, as nuclear monotony is a histomorphological hallmark of oligodendrogliomas. For the purpose of this study, we defined pO as diffusely infiltrating gliomas felt to be of oligodendroglial rather than astrocytic differentiation and characterized by the presence of multinucleate tumor giant cells and/or nuclear pleomorphism. In a total of nine patients, we identified tumors consistent with this working definition. All tumors were high-grade. We characterized these with respect to clinical, histomorphological and genetic features. Despite clinical and genetic heterogeneity, we identified a subset of tumors of bona fide oligodendroglial differentiation as characterized by combined loss of heterozygosity of chromosome arms 1p and 19q (LOH 1p19q). Those tumors that lacked LOH 1p19q showed a high frequency of IDH1 mutations and loss of alpha thalassemia/mental retardation syndrome X-linked gene (ATRX) immunoreactivity, indicating a possible phenotypic convergence of true oligodendrogliomas and gliomas of the alternative lengthening of telomeres (ALT) pathway. p53 alterations were common irrespective of the 1p19q status. Histomorphologically, the tumors featured interspersed bizarre multinucleate giant tumor cells, while the background population varied from monotonous to significantly pleomorphic. Our findings indicate, that a rare polymorphous - or "giant cell" - variant of oligodendroglioma does indeed exist.

Combined ATRX/IDH1 immunohistochemistry predicts genotype of oligoastrocytomas

AIMS To assess whether in oligoastrocytomas ATRX deficiency - as a surrogate of the alternative lengthening of telomeres (ALT) pathway - has a role in predicting the presence or absence of loss of heterozygosity of 1p and 19q (LOH), the genetic signature of oligodendroglial differentiation and a favourable prognostic marker. METHODS AND RESULTS A series of 54 oligoastrocytomas were investigated by immunohistochemistry as well as microsatellite analysis for LOH 1p19q. Genetic findings were correlated with morphological assessment. CONCLUSIONS ATRX deficiency was mutually exclusive with LOH. Conversely, ATRX-proficient tumours immunoreactive for R132H-mutant isocitrate dehydrogenase 1 (IDH1) showed a high rate (85%) of LOH. A more oligodendroglioma-like morphology was associated with a higher rate of LOH even in the morphologically ambiguous group of oligoastrocytomas. Our findings support the concept that oligoastrocytomas represent a morphological grey zone rather than a group of truly "mixed" or "intermediate" tumours. More precise classification of diffuse gliomas may also improve grading of borderline cases. We propose an immunohistochemical algorithm for classification of morphologically ambiguous diffuse gliomas. This article is protected by copyright. All rights reserved.

Ignoring heterozygous sites biases phylogenomic estimates of divergence times: implications for the evolutionary history of Microtus voles

Phylogenetic reconstruction of the evolutionary history of closely related organisms may be difficult because of the presence of unsorted lineages and of a relatively high proportion of heterozygous sites that are usually not handled well by phylogenetic programs. Genomic data may provide enough fixed polymorphisms to resolve phylogenetic trees, but the diploid nature of sequence data remains analytically challenging. Here, we performed a phylogenomic reconstruction of the evolutionary history of the common vole (Microtus arvalis) with a focus on the influence of heterozygosity on the estimation of intraspecific divergence times. We used genome-wide sequence information from 15 voles distributed across the European range. We provide a novel approach to integrate heterozygous information in existing phylogenetic programs by repeated random haplotype sampling from sequences with multiple unphased heterozygous sites. We evaluated the impact of the use of full, partial, or no heterozygous information for tree reconstructions on divergence time estimates. All results consistently showed four deep and strongly supported evolutionary lineages in the vole data. These lineages undergoing divergence processes split only at the end or after the last glacial maximum based on calibration with radiocarbon-dated paleontological material. However, the incorporation of information from heterozygous sites had a significant impact on absolute and relative branch length estimations. Ignoring heterozygous information led to an overestimation of divergence times between the evolutionary lineages of M. arvalis. We conclude that the exclusion of heterozygous sites from evolutionary analyses may cause biased and misleading divergence time estimates in closely related taxa.