Regeneration and characterization of somatic hybrids combining sweet orange and mandarin/mandarin hybrid cultivars for citrus scion improvement

Protoplast fusion between sweet orange and mandarin/mandarin hybrids scion cultivars was performed following the model "diploid embryogenic callus protoplast + diploid mesophyll-derived protoplast". Protoplasts were isolated from embryogenic calli of 'Pera' and 'Westin' sweet orange cultivars (Citrus sinensis) and from young leaves of 'Fremont', Nules', and 'Thomas' mandarins (C. reticulata), and 'Nova' tangelo [C. reticulata x (C. paradisi x C. reticulata)]. The regenerated plants were characterized based on their leaf morphology (thickness), ploidy level, and simple sequence repeat (SSR) molecular markers. Plants were successfully generated only when 'Pera' sweet orange was used as the embryogenic parent. Fifteen plants were regenerated being 7 tetraploid and 8 diploid. Based on SSR molecular markers analyses all 7 tetraploid regenerated plants revealed to be allotetraploids (somatic hybrids), including 2 from the combination of 'Pera' sweet orange + 'Fremont' mandarin, 3 'Pera' sweet orange + 'Nules' mandarin, and 2 'Pera' sweet orange + 'Nova' tangelo, and all the diploid regenerated plants showed the 'Pera' sweet orange marker profile. Somatic hybrids were inoculated with Alternaria alternata and no disease symptoms were detected 96 h post-inoculation. This hybrid material has the potential to be used as a tetraploid parent in interploid crosses for citrus scion breeding.

The effect of ploidy on the fatty acid profile during the reproductive cycle of female rainbow trout (Oncorhynchus mykiss)

Rainbow trout Oncorhynchus mykiss triploids are regularly produced in fish farms to improve growth because the triploid females do not develop ovaries during the reproductive cycle. In this study, the tissue fatty acid allocations in triploid and diploid females were compared during the reproductive cycle to determine whether the ploidy influences the fatty acid profile of fish produced in aquaculture. The ovaries, liver, and white muscle fatty acid contents of diploid and triploid females were analyzed during the reproductive cycle. Diploid females tend to accumulate more polyunsaturated fatty acids than triploids during some phases of the reproductive cycle, and this profile was compensated by an increase in saturated and monounsaturated fatty acids in triploid females. Arachidonic acid (C20:4n6) was the main n6 polyunsaturated fatty acid in the ovaries of diploid females during the most advanced phases of the reproductive cycle, and docosahexaenoic acid (C22:6n3) was the main n3 polyunsaturated fatty acid. In triploid females, the percentage of both of these polyunsaturated fatty acids was lower than in diploid females during the most advanced phases of the reproductive cycle. In general, the lack of ovary development altered the hepatic synthesis of some fatty acids, mainly monounsaturated and polyunsaturated fatty acids, decreasing the content of the main fatty acids in the white muscle and, consequently, the mobilization of these specific fatty acids to the ovaries.

Cytogenetic and molecular analyses in troglobitic and epigean species of Pimelodella (Siluriformes: Heptapteridae) from Brazil

Samples from seven different locations of the genus Pimelodella were genetically examined, two caves (exclusively subterranean, upper Tocantins River and Sao Francisco River) and five epigean (from upper Parana River basin). Cytogenetic analyses revealed the same diploid number (2n=46) for all species besides similarities in both number and location of nucleolar organizer regions and C bands. FISH with 5S rDNA probes and CMA(3) staining indicated significant differences among the studied species. Application of PCR-RFLP in ATPase 6 and 8 mitochondrial genes allowed building a minimum evolution phenogram identifying the close evolutionary relationship among groups. Both chromosomal and molecular data were useful to infer the relationships among studied Pimelodella species.

Genetic diversity and population structure of Musa accessions in ex situ conservation

Abstract Background Banana cultivars are mostly derived from hybridization between wild diploid subspecies of Musa acuminata (A genome) and M. balbisiana (B genome), and they exhibit various levels of ploidy and genomic constitution. The Embrapa ex situ Musa collection contains over 220 accessions, of which only a few have been genetically characterized. Knowledge regarding the genetic relationships and diversity between modern cultivars and wild relatives would assist in conservation and breeding strategies. Our objectives were to determine the genomic constitution based on Internal Transcribed Spacer (ITS) regions polymorphism and the ploidy of all accessions by flow cytometry and to investigate the population structure of the collection using Simple Sequence Repeat (SSR) loci as co-dominant markers based on Structure software, not previously performed in Musa. Results From the 221 accessions analyzed by flow cytometry, the correct ploidy was confirmed or established for 212 (95.9%), whereas digestion of the ITS region confirmed the genomic constitution of 209 (94.6%). Neighbor-joining clustering analysis derived from SSR binary data allowed the detection of two major groups, essentially distinguished by the presence or absence of the B genome, while subgroups were formed according to the genomic composition and commercial classification. The co-dominant nature of SSR was explored to analyze the structure of the population based on a Bayesian approach, detecting 21 subpopulations. Most of the subpopulations were in agreement with the clustering analysis. Conclusions The data generated by flow cytometry, ITS and SSR supported the hypothesis about the occurrence of homeologue recombination between A and B genomes, leading to discrepancies in the number of sets or portions from each parental genome. These phenomenons have been largely disregarded in the evolution of banana, as the “single-step domestication” hypothesis had long predominated. These findings will have an impact in future breeding approaches. Structure analysis enabled the efficient detection of ancestry of recently developed tetraploid hybrids by breeding programs, and for some triploids. However, for the main commercial subgroups, Structure appeared to be less efficient to detect the ancestry in diploid groups, possibly due to sampling restrictions. The possibility of inferring the membership among accessions to correct the effects of genetic structure opens possibilities for its use in marker-assisted selection by association mapping.

Entstehung und Diversifizierung der mittel- und westeuropäischen Salicornia-Arten

HintergrundrnDie hygrohalophytische Gattung Salicornia ist in Mittel- und Westeuropa durch vier nah verwandte, sympatrisch vorkommende Arten vertreten. Es handelt sich um die zwei tetraploiden Arten S. procumbens und S. stricta und die diploiden Arten S. europaea und S. ramosissima. Morphologisch lassen sich die Arten zwar nur schwer voneinander unterscheiden, die morphologische Variation ist aber wiederum so hoch, dass mehrere distinkte Arten/Morphotypen unterschieden werden können. Bezüglich ihrer Verteilung im hochdynamischen Lebensraum Salzwiese findet man die verschiedenen Arten/Morphotypen in überlappenden Bereichen des Habitats. Ihr relativ vorhersagbares Auftreten entlang eines ökologischen Gradienten innerhalb ihres Lebensraumes scheint jedoch für eine ökologische Differenzierung der verschiedenen Arten/Morphotypen zu sprechen. Aufgrund des sympatrischen Vorkommens der scheinbar ökologisch und morphologisch differenzierten Morphotypen stellt sich die Frage, durch welche Prozesse diese entstanden sein könnten (genetische und ökologische Differenzierung) aber auch welche Prozesse die dauerhafte Koexistenz der Arten (reproduktive Isolationsmechanismen) aufrechterhalten.rnZielsetzungrnZiel dieser Arbeit war es, die Entstehung und Diversifizierung der mittel- und westeuropäischen Salicornia-Arten anhand von molekulargenetischen, ökologischen und reproduktionsbiologischen Methoden zu untersuchen.rnMethodenrnAnhand einer AFLP-Fragmentanalyse mit 89 Herkünften aus Großbritannien, Frankreich und Deutschland wurden molekulare Phylogenien erstellt sowie eine Hauptkomponenten- und Clusteranalyse durchgeführt. Um die ökologische Differenzierung und phänotypische Plastizität der vier Arten/Morphotypen zu untersuchen wurde ein reziprokes Transplantationsexperiment durchgeführt. Um die reproduktiven Isolationsmechanismen der Arten/Morphotypen zu untersuchen, wurden verschiedene Beobachtungen und Experimente durchgeführt.rnErgebnissernDie molekularen Analysen konnten zwar die beiden Artengruppen (Ploidiestufen) trennen, lieferten aber innerhalb dieser weder ein taxonomisches noch ein geographisches Signal. Akzessionen mit identischer Morphologie aus der gleichen Population verteilten sich in den Analysen in verschiedene genetische Cluster. Identische Morphotypen aus verschiedenen geographischen Regionen gruppieren teilweise zusammen. Das Transplantationsexperiment zeigte für die beiden tetraploiden Arten S. procumbens und S. stricta eine deutliche ökologische Differenzierung, bei S. procumbens in Form von verminderter Fitness und einer beschleunigten Phänologie, bei S. stricta nur in Form einer veränderten Phänologie. Bezüglich der Plastizität zeigten beide tetraploiden Arten eine konstante Morphologie. Die beiden diploiden Taxa S. europaea und S. ramosissima zeigten weder eine klare ökologische Differenzierung noch eine konstante Morphologie. Bezüglich der Reproduktionsbiologie konnte bestätigt werden, dass Selbstung bei allen Taxa der hauptsächliche Reproduktionsmodus ist. Bei den tetraploiden Taxa zeigte sich zwar ein geringes Maß an Fremdbefruchtung, bei den diploiden Taxa führen dagegen morphologische Besonderheiten zu hochgradiger Selbstung.rnRésumérnDie in Mittel- und Westeuropa vorkommenden Salicornia-Arten stellen keine evolutionären Einheiten dar. Die beiden tetraploiden Taxa sollten auf Grund ihrer parallelen Entstehung und ökologischen Differenzierung als Ökotypen angesprochen werden. Beide Ökotypen weisen ein hohes Ausbreitungspotential aus und persistieren als Inzuchtlinien mit geringem Anteil an Fremdbestäubung. Die diploiden Taxa sind weder ökologisch differenziert noch morphologisch stabil und sollten deshalb als nur ein morphologisch sehr variables, aus zahlreichen weitverbreiteten Inzuchtlinien bestehendes Taxon angesehen werden.

Alternative reproductive tactics, sexual selection and the effects of inbreeding in the ant Hypoponera opacior

Inbreeding can lead to a fitness reduction due to the unmasking of deleterious recessive alleles and the loss of heterosis. Therefore, most sexually reproducing organisms avoid inbreeding, often by disperal. Besides the avoidance of inbreeding, dispersal lowers intraspecific competition on a local scale and leads to a spreading of genotypes into new habitats. In social insects, winged reproductives disperse and mate during nuptial flights. Therafter, queens independently found a new colony. However, some species also produce wingless sexuals as an alternative reproductive tactic. Wingless sexuals mate within or close to their colony and queens either stay in the nest or they found a new colony by budding. During this dependent colony foundation, wingless queens are accompanied by a fraction of nestmate workers. The production of wingless reproductives therefore circumvents the risks associated with dispersal and independent colony foundation. However, the absence of dispersal can lead to inbreeding and local competition.rnIn my PhD-project, I investigated the mating biology of Hypoponera opacior, an ant that produces winged and wingless reproductives in a population in Arizona. Besides the investigation of the annual reproductive cycle, I particularly focused on the consequences of wingless reproduction. An analysis of sex ratios in wingless sexuals should reveal the relative importance of local resource competition among queens (that mainly compete for the help of workers) and local mate competition among males. Further, sexual selection was expected to act on wingless males that were previously found to mate with and mate-guard pupal queens in response to local mate competition. We studied whether males are able to adapt their mating behaviour to the current competitive situation in the nest and which traits are under selection in this mating situation. Last, we investigated the extent and effects of inbreeding. As the species appeared to produce non-dispersive males and queens quite frequently, we assumed to find no or only weak negative effects of inbreeding and potentially mechanisms that moderate inbreeding levels despite frequent nest-matings.rnWe found that winged and wingless males and queens are produced during two separate seasons of the year. Winged sexuals emerge in early summer and conduct nuptial flights in July, when climate conditions due to frequent rainfalls lower the risks of dispersal and independent colony foundation. In fall, wingless sexuals are produced that reproduce within the colonies leading to an expansion on the local scale. The absence of dispersal during this second reproductive season resulted in a local genetic population viscosity and high levels of inbreeding within the colonies. Male-biased sex ratios in fall indicated a greater importance of local resource competition among queens than local mate competition among males. Males were observed to adjust mate-guarding durations to the competitive situation (i.e. the number of competing males and pupae) in the nest, an adaptation that helps maximising their reproductive success. Further, sexual selection was found to act on the timing of emergence as well as on body size in these males, i.e. earlier emerging and larger males show a higher mating success. Genetic analyses revealed that wingless males do not actively avoid inbreeding by choosing less related queens as mating partners. Further, we detected diploid males, a male type that is produced instead of diploid females if close relatives mate. In contrast to many other Hymenopteran species, diploid males were here viable and able to sire sterile triploid offspring. They did not differ in lifespan, body size and mating success from “normal” haploid males. Hence, diploid male production in H. opacior is less costly than in other social Hymenopteran species. No evidence of inbreeding depression was found on the colony level but more inbred colonies invested more resources into the production of sexuals. This effect was more pronounced in the dispersive summer generation. The increased investment in outbreeding sexuals can be regarded as an active strategy to moderate the extent and effects of inbreeding. rnIn summary, my thesis describes an ant species that has evolved alternative reproductive tactics as an adaptation to seasonal environmental variations. Hereby, the species is able to maintain its adaptive mating system without suffering from negative effects due to the absence of dispersal flights in fall.rn

Multilocus analysis of genetic divergence between outcrossing Arabidopsis species: evidence of genome-wide admixture

P>Outcrossing Arabidopsis species that diverged from their inbreeding relative Arabidopsis thaliana 5 million yr ago and display a biogeographical pattern of interspecific sympatry vs intraspecific allopatry provides an ideal model for studying impacts of gene introgression and polyploidization on species diversification. Flow cytometry analyses detected ploidy polymorphisms of 2x and 4x in Arabidopsis lyrata ssp. kamchatica of Taiwan. Genomic divergence between species/subspecies was estimated based on 98 randomly chosen nuclear genes. Multilocus analyses revealed a mosaic genome in diploid A. l. kamchatica composed of Arabidopsis halleri-like and A. lyrata-like alleles. Coalescent analyses suggest that the segregation of ancestral polymorphisms alone cannot explain the high inconsistency between gene trees across loci, and that gene introgression via diploid A. l. kamchatica likely distorts the molecular phylogenies of Arabidopsis species. However, not all genes migrated across species freely. Gene ontology analyses suggested that some nonmigrating genes were constrained by natural selection. High levels of estimated ancestral polymorphisms between A. halleri and A. lyrata suggest that gene flow between these species has not completely ceased since their initial isolation. Polymorphism data of extant populations also imply recent gene flow between the species. Our study reveals that interspecific gene flow affects the genome evolution in Arabidopsis.

eIF4E is an important determinant of adhesion and pseudohyphal growth of the yeast S. cerevisiae

eIF4E, the cytoplasmatic cap-binding protein, is required for efficient cap-dependent translation. We have studied the influence of mutations that alter the activity and/or expression level of eIF4E on haploid and diploid cells in the yeast S. cerevisiae. Temperature-sensitive eIF4E mutants with reduced levels of expression and reduced cap-binding affinity clearly show a loss in haploid adhesion and diploid pseudohyphenation upon starvation for nitrogen. Some of these mutations affect the interaction of the cap-structure of mRNAs with the cap-binding groove of eIF4E. The observed reduction in adhesive and pseudohyphenating properties is less evident for an eIF4E mutant that shows reduced interaction with p20 (an eIF4E-binding protein) or for a p20-knockout mutant. Loss of adhesive and pseudohyphenating properties was not only observed for eIF4E mutants but also for knockout mutants of components of eIF4F such as eIF4B and eIF4G1. We conclude from these experiments that mutations that affect components of the eIF4F-complex loose properties such as adhesion and pseudohyphal differentiation, most likely due to less effective translation of required mRNAs for such processes.

Detection of nonrandom association of alleles from the distribution of the number of heterozygous loci in a sample.

The distribution of the number of heterozygous loci in two randomly chosen gametes or in a random diploid zygote provides information regarding the nonrandom association of alleles among different genetic loci. Two alternative statistics may be employed for detection of nonrandom association of genes of different loci when observations are made on these distributions: observed variance of the number of heterozygous loci (s2k) and a goodness-of-fit criterion (X2) to contrast the observed distribution with that expected under the hypothesis of random association of genes. It is shown, by simulation, that s2k is statistically more efficient than X2 to detect a given extent of nonrandom association. Asymptotic normality of s2k is justified, and X2 is shown to follow a chi-square (chi 2) distribution with partial loss of degrees of freedom arising because of estimation of parameters from the marginal gene frequency data. Whenever direct evaluations of linkage disequilibrium values are possible, tests based on maximum likelihood estimators of linkage disequilibria require a smaller sample size (number of zygotes or gametes) to detect a given level of nonrandom association in comparison with that required if such tests are conducted on the basis of s2k. Summarization of multilocus genotype (or haplotype) data, into the different number of heterozygous loci classes, thus, amounts to appreciable loss of information.

Stress-induced targeting of molecular chaperones in the yeast Saccharomyces cerevisiae

The eukaryotic stress response is an essential mechanism that helps protect cells from a variety of environmental stresses. Cell death can result if cells are not able to properly adapt and protect themselves against adverse stress conditions. Failure to properly deal with stress has implications in human diseases including neurodegenerative disorders and distinct cancers, emphasizing the importance of understanding the eukaryotic stress response in detail. As part of this response, expression of a battery of heat shock proteins (HSP) is induced, which act as molecular chaperones to assist in the repair or triage of unfolded proteins. The 90-kDa HSP (Hsp90) operates in the context of a multi-chaperone complex to promote the maturation of nuclear and cytoplasmic clients. I have discovered that Hsp90 and the co-chaperone Sba1 accumulate in the nucleus of quiescent Saccharomyces cerevisiae cells in a karyopherin-dependent manner. I isolated nuclear accumulation- defective HSP82 mutant alleles to probe the nature of this targeting event and identified a mutant with a single amino acid substitution (I578F) sufficient to prevent nuclear accumulation of Hsp90 in quiescent cells. Diploid hsp82-I578F cells exhibited pronounced defects in spore wall construction and maturation, resulting in catastrophic sporulation. The mislocalization and sporulation phenotypes were shared by another previously identified HSP82 mutant allele, further linking localization to Hsp90 functional status. Pharmacological inhibition of Hsp90 with macbecin in sporulating diploid cells also blocked spore formation, underscoring the importance of this chaperone in this developmental program. The yeast molecular chaperone Hsp104 is a member of the Hsp100 superfamily of AAA+ ATPases. Unlike the Hsp90 family of chaperones, Hsp104 is not restricted to a specific set of client proteins, but rather assists in reactivating stress-denatured proteins by solubilizing protein aggregates. I have discovered that Hsp104, along with the Hsp70 chaperone, Ssa1, and the sHSP Hsp26 accumulate into RNA processing bodies (P- bodies) and stress granules, sites of mRNA metabolism. I found that Hsp104 recruits both Ssa1 and Hsp26 to P-bodies and that these three chaperones are required for stress granule formation. These findings suggest a possible role for chaperones in mRNA metabolism by aiding in the assembly, disassembly or conversion of these enigmatic mRNP complexes. Taken together, the work presented in this dissertation serves to better understand the eukaryotic stress response by illustrating the importance of subcellular-chaperone localization in key biological processes.

Hsp90 nuclear accumulation in quiescence is linked to chaperone function and spore development in yeast.

The 90-kDa heat-shock protein (Hsp90) operates in the context of a multichaperone complex to promote maturation of nuclear and cytoplasmic clients. We have discovered that Hsp90 and the cochaperone Sba1/p23 accumulate in the nucleus of quiescent Saccharomyces cerevisiae cells. Hsp90 nuclear accumulation was unaffected in sba1Delta cells, demonstrating that Hsp82 translocates independently of Sba1. Translocation of both chaperones was dependent on the alpha/beta importin SRP1/KAP95. Hsp90 nuclear retention was coincident with glucose exhaustion and seems to be a starvation-specific response, as heat shock or 10% ethanol stress failed to elicit translocation. We generated nuclear accumulation-defective HSP82 mutants to probe the nature of this targeting event and identified a mutant with a single amino acid substitution (I578F) sufficient to retain Hsp90 in the cytoplasm in quiescent cells. Diploid hsp82-I578F cells exhibited pronounced defects in spore wall construction and maturation, resulting in catastrophic sporulation. The mislocalization and sporulation phenotypes were shared by another previously identified HSP82 mutant allele. Pharmacological inhibition of Hsp90 with macbecin in sporulating diploid cells also blocked spore formation, underscoring the importance of this chaperone in this developmental program.

Studies on the {\it in vivo} function of a phosphatidylinositol transfer protein from {\it Saccharomyces cerevisiae}

Phosphatidylinositol transfer proteins (PI-TP's) catalyze the transfer of phosphatidylinositol and phosphatidylcholine between membranes in vitro. However the in vivo function of these proteins is unknown. In this thesis we have used a combined biochemical and genetic approach to determine the importance of PI-TP in vivo. An oligonucleotide based on the amino terminal sequence of the PI-TP from Saccharomyces cerevisiae, was used to screen a yeast genomic library for the gene encoding PI-TP (PIT1 gene). Yeast strains transformed with the positive clones showed overproduction of transfer activities and transfer protein in the 100,000 x g supernatants. The 5$\sp\prime$ terminus of the PIT1 gene correlates with the predicted codons for residues 3-30 of the determined protein sequence. Tetrad analysis of a heterozygous diploid (PIT1/pit1::LEU2) revealed that the PIT1 gene is essential for cell growth. Non-viable spores could be rescued by transformation of the above diploid prior to sporulation, with a plasmid borne copy of the wild type gene. Sequencing of the entire PIT1 gene has revealed that the PIT1 gene is identical to the SEC14 gene. The sec14 ts mutant which exhibits conditional defects at the Golgi stage of protein secretion, is also temperature sensitive for PI-TP activity in vitro. These findings represent the first instance in which a physiological function has been assigned to any phospholipid transfer protein. ^

Semi-dwarfism and lodging tolerance in tef (Eragrostis tef) is linked to a mutation in the α-Tubulin 1 gene

Genetic improvement of native crops is a new and promising strategy to combat hunger in the developing world. Tef is the major staple food crop for approximately 50 million people in Ethiopia. As an indigenous cereal, it is well adapted to diverse climatic and soil conditions; however, its productivity is extremely low mainly due to susceptibility to lodging. Tef has a tall and weak stem, liable to lodge (or fall over), which is aggravated by wind, rain, or application of nitrogen fertilizer. To circumvent this problem, the first semi-dwarf lodging-tolerant tef line, called kegne, was developed from an ethyl methanesulphonate (EMS)-mutagenized population. The response of kegne to microtubule-depolymerizing and -stabilizing drugs, as well as subsequent gene sequencing and segregation analysis, suggests that a defect in the α-Tubulin gene is functionally and genetically tightly linked to the kegne phenotype. In diploid species such as rice, homozygous mutations in α-Tubulin genes result in extreme dwarfism and weak stems. In the allotetraploid tef, only one homeologue is mutated, and the presence of the second intact α-Tubulin gene copy confers the agriculturally beneficial semi-dwarf and lodging-tolerant phenotype. Introgression of kegne into locally adapted and popular tef cultivars in Ethiopia will increase the lodging tolerance in the tef germplasm and, as a result, will improve the productivity of this valuable crop.

Sex-Specific Differences in Pathogen Susceptibility in Honey Bees (Apis mellifera)

Sex-related differences in susceptibility to pathogens are a common phenomenon in animals. In the eusocial Hymenoptera the two female castes, workers and queens, are diploid and males are haploid. The haploid susceptibility hypothesis predicts that haploid males are more susceptible to pathogen infections compared to females. Here we test this hypothesis using adult male (drone) and female (worker) honey bees (Apis mellifera), inoculated with the gut endoparasite Nosema ceranae and/or black queen cell virus (BQCV). These pathogens were chosen due to previously reported synergistic interactions between Nosema apis and BQCV. Our data do not support synergistic interactions between N. ceranae and BQCV and also suggest that BQCV has limited effect on both drone and worker health, regardless of the infection level. However, the data clearly show that, despite lower levels of N. ceranae spores in drones than in workers, Nosema-infected drones had both a higher mortality and a lower body mass than non-infected drones, across all treatment groups, while the mortality and body mass of worker bees were largely unaffected by N. ceranae infection, suggesting that drones are more susceptible to this pathogen than workers. In conclusion, the data reveal considerable sex-specific differences in pathogen susceptibility in honey bees and highlight the importance of ultimate measures for determining susceptibility, such as mortality and body quality, rather than mere infection levels

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.