p66ShcA functions as a contextual promoter of breast cancer metastasis.

BACKGROUND: The p66ShcA redox protein is the longest isoform of the Shc1 gene and is variably expressed in breast cancers. In response to a variety of stress stimuli, p66ShcA becomes phosphorylated on serine 36, which allows it to translocate from the cytoplasm to the mitochondria where it stimulates the formation of reactive oxygen species (ROS). Conflicting studies suggest both pro- and anti-tumorigenic functions for p66ShcA, which prompted us to examine the contribution of tumor cell-intrinsic functions of p66ShcA during breast cancer metastasis. METHODS: We tested whether p66ShcA impacts the lung-metastatic ability of breast cancer cells. Breast cancer cells characteristic of the ErbB2+/luminal (NIC) or basal (4T1) subtypes were engineered to overexpress p66ShcA. In addition, lung-metastatic 4T1 variants (4T1-537) were engineered to lack endogenous p66ShcA via Crispr/Cas9 genomic editing. p66ShcA null cells were then reconstituted with wild-type p66ShcA or a mutant (S36A) that cannot translocate to the mitochondria, thereby lacking the ability to stimulate mitochondrial-dependent ROS production. These cells were tested for their ability to form spontaneous metastases from the primary site or seed and colonize the lung in experimental (tail vein) metastasis assays. These cells were further characterized with respect to their migration rates, focal adhesion dynamics, and resistance to anoikis in vitro. Finally, their ability to survive in circulation and seed the lungs of mice was assessed in vivo. RESULTS: We show that p66ShcA increases the lung-metastatic potential of breast cancer cells by augmenting their ability to navigate each stage of the metastatic cascade. A non-phosphorylatable p66ShcA-S36A mutant, which cannot translocate to the mitochondria, still potentiated breast cancer cell migration, lung colonization, and growth of secondary lung metastases. However, breast cancer cell survival in the circulation uniquely required an intact p66ShcA S36 phosphorylation site. CONCLUSION: This study provides the first evidence that both mitochondrial and non-mitochondrial p66ShcA pools collaborate in breast cancer cells to promote their maximal metastatic fitness.

Dynamic Neutrophil-to-Lymphocyte Ratio: A Novel Prognosis Measure for Triple-Negative Breast Cancer.

BACKGROUND: The neutrophil-to-lymphocyte ratio (NLR) is a measure of systemic inflammation and a prognostic factor for multiple malignancies. This study assesses the value of the NLR as an independent prognostic marker in triple-negative breast cancer (TNBC) and explores the association between dynamic NLR changes and patient outcomes. METHODS: The study retrospectively analyzed a prospectively maintained database including patients 18 to 80 years old with TNBC treated at the authors' institution between 2006 to 2016. Clinical and demographic data were collected, including blood test results and treatments received. Age at diagnosis, stage of disease, and NLR scores were tested for association with overall and disease-free survival in uni- and multivariate Cox models. RESULTS: The inclusion criteria were met by 329 women with a median age of 58. Most of the patients had early-stage disease (30.1% with stage 1 and 47% with stage 2 malignancy). An NLR higher than 2.84 at diagnosis was associated with decreased overall survival (hazard ratio [HR], 1.8; 95% confidence interval [CI], 1.023-3.176), whereas an NLR higher than 7.82 at any time during the follow-up period was a strong predictor of 5-year mortality (HR, 10.76; 95% CI, 4.193-26.58), independent of age or stage of disease. Patients who experienced recurrence had a higher NLR than their counterparts during the 6 months before recurrence. The NLR also significantly rose during the final 18 months of life (p < 0.01). CONCLUSION: The NLR is an important prognostic marker in TNBC, both at diagnosis and during the course of the disease. Moreover, dynamic changes in NLR strongly correlate with disease recurrence and the time of death.

Promiscuous DNA-binding of a mutant zinc finger protein corrupts the transcriptome and diminishes cell viability.

The rules of engagement between zinc finger transcription factors and DNA have been partly defined by in vitro DNA-binding and structural studies, but less is known about how these rules apply in vivo. Here, we demonstrate how a missense mutation in the second zinc finger of Krüppel-like factor-1 (KLF1) leads to degenerate DNA-binding specificity in vivo, resulting in ectopic transcription and anemia in the Nan mouse model. We employed ChIP-seq and 4sU-RNA-seq to identify aberrant DNA-binding events genome wide and ectopic transcriptional consequences of this binding. We confirmed novel sequence specificity of the mutant recombinant zinc finger domain by performing biophysical measurements of in vitro DNA-binding affinity. Together, these results shed new light on the mechanisms by which missense mutations in DNA-binding domains of transcription factors can lead to autosomal dominant diseases.

CLIC5: a novel ETV6 target gene in childhood acute lymphoblastic leukemia.

The most common rearrangement in childhood precursor B-cell acute lymphoblastic leukemia is the t(12;21)(p13;q22) translocation resulting in the ETV6-AML1 fusion gene. A frequent concomitant event is the loss of the residual ETV6 allele suggesting a critical role for the ETV6 transcriptional repressor in the etiology of this cancer. However, the precise mechanism through which loss of functional ETV6 contributes to disease pathogenesis is still unclear. To investigate the impact of ETV6 loss on the transcriptional network and to identify new transcriptional targets of ETV6, we used whole transcriptome analysis of both pre-B leukemic cell lines and patients combined with chromatin immunoprecipitation. Using this integrative approach, we identified 4 novel direct ETV6 target genes: CLIC5, BIRC7, ANGPTL2 and WBP1L To further evaluate the role of chloride intracellular channel protein CLIC5 in leukemogenesis, we generated cell lines overexpressing CLIC5 and demonstrated an increased resistance to hydrogen peroxide-induced apoptosis. We further described the implications of CLIC5's ion channel activity in lysosomal-mediated cell death, possibly by modulating the function of the transferrin receptor with which it colocalizes intracellularly. For the first time, we showed that loss of ETV6 leads to significant overexpression of CLIC5, which in turn leads to decreased lysosome-mediated apoptosis. Our data suggest that heightened CLIC5 activity could promote a permissive environment for oxidative stress-induced DNA damage accumulation, and thereby contribute to leukemogenesis.

An overlapping set of genes is regulated by both NFIB and the glucocorticoid receptor during lung maturation.

BACKGROUND: Lung maturation is a late fetal developmental event in both mice and humans. Because of this, lung immaturity is a serious problem in premature infants. Disruption of genes for either the glucocorticoid receptor (Nr3c1) or the NFIB transcription factors results in perinatal lethality due to lung immaturity. In both knockouts, the phenotype includes excess cell proliferation, failure of saccularization and reduced expression of markers of epithelial differentiation. This similarity suggests that the two genes may co-regulate a specific set of genes essential for lung maturation. RESULTS: We analyzed the roles of these two transcription factors in regulating transcription using ChIP-seq data for NFIB, and RNA expression data and motif analysis for both. Our new ChIP-seq data for NFIB in lung at E16.5 shows that NFIB binds to a NFI motif. This motif is over-represented in the promoters of genes that are under-expressed in Nfib-KO mice at E18.5, suggesting an activator role for NFIB. Using available microarray data from Nr3c1-KO mice, we further identified 52 genes that are under-expressed in both Nfib and Nr3c1 knockouts, an overlap which is 13.1 times larger than what would be expected by chance. Finally, we looked for enrichment of 738 recently published transcription factor motifs in the promoters of these putative target genes and found that the NFIB and glucocorticoid receptor motifs were among the most enriched, suggesting that a subset of these genes may be directly activated by Nfib and Nr3c1. CONCLUSIONS: Our data provide the first evidence for Nfib and Nr3c1 co-regulating genes related to lung maturation. They also establish that the in vivo DNA-binding specificity of NFIB is the same as previously seen in vitro, and highly similar to that of the other NFI-family members NFIA, NFIC and NFIX.

Mapping the single cell spatial immune landscapes of the melanoma microenvironment.

Melanoma is a highly immunogenic malignancy with an elevated mutational burden, diffuse lymphocytic infiltration, and one of the highest response rates to immune checkpoint inhibitors (ICIs). However, over half of all late-stage patients treated with ICIs will either not respond or develop progressive disease. Spatial imaging technologies are being increasingly used to study the melanoma tumor microenvironment (TME). The goal of such studies is to understand the complex interplay between the stroma, melanoma cells, and immune cell-types as well as their association with treatment response. Investigators seeking a better understanding of the role of cell location within the TME and the importance of spatial expression of biomarkers are increasingly turning to highly multiplexed imaging approaches to more accurately measure immune infiltration as well as to quantify receptor-ligand interactions (such as PD-1 and PD-L1) and cell-cell contacts. CyTOF-IMC (Cytometry by Time of Flight - Imaging Mass Cytometry) has enabled high-dimensional profiling of melanomas, allowing researchers to identify complex cellular subpopulations and immune cell interactions with unprecedented resolution. Other spatial imaging technologies, such as multiplexed immunofluorescence and spatial transcriptomics, have revealed distinct patterns of immune cell infiltration, highlighting the importance of spatial relationships, and their impact in modulating immunotherapy responses. Overall, spatial imaging technologies are just beginning to transform our understanding of melanoma biology, providing new avenues for biomarker discovery and therapeutic development. These technologies hold great promise for advancing personalized medicine to improve patient outcomes in melanoma and other solid malignancies.

Melanomas with concurrent BRAF non-p.V600 and NF1 loss-of-function mutations are targetable by BRAF/MEK inhibitor combination therapy.

Although combination BRAF/MEK inhibition has produced significant survival benefits for BRAF p.V600 mutant melanomas, targeted therapies approved for BRAF non-p.V600 mutant melanomas remain limited. Through the analysis of 772 cutaneous melanoma exomes, we reveal that BRAF non-p.V600 mutations co-occurs more frequently with NF1 loss, but not with oncogenic NRAS mutations, than expected by chance. We present cell signaling data, which demonstrate that BRAF non-p.V600 mutants can signal as monomers and dimers within an NF1 loss context. Concordantly, BRAF inhibitors that inhibit both monomeric and dimeric BRAF synergize with MEK inhibition to significantly reduce cell viability in vitro and tumor growth in vivo in BRAF non-p.V600 mutant melanomas with co-occurring NF1 loss-of-function mutations. Our data suggest that patients harboring BRAF non-p.V600 mutant melanomas may benefit from current FDA-approved BRAF/MEK inhibitor combination therapy currently reserved for BRAF p.V600 mutant patients.

Spatially mapping the immune landscape of melanoma using imaging mass cytometry.

Melanoma is an immunogenic cancer with a high response rate to immune checkpoint inhibitors (ICIs). It harbors a high mutation burden compared with other cancers and, as a result, has abundant tumor-infiltrating lymphocytes (TILs) within its microenvironment. However, understanding the complex interplay between the stroma, tumor cells, and distinct TIL subsets remains a substantial challenge in immune oncology. To properly study this interplay, quantifying spatial relationships of multiple cell types within the tumor microenvironment is crucial. To address this, we used cytometry time-of-flight (CyTOF) imaging mass cytometry (IMC) to simultaneously quantify the expression of 35 protein markers, characterizing the microenvironment of 5 benign nevi and 67 melanomas. We profiled more than 220,000 individual cells to identify melanoma, lymphocyte subsets, macrophage/monocyte, and stromal cell populations, allowing for in-depth spatial quantification of the melanoma microenvironment. We found that within pretreatment melanomas, the abundance of proliferating antigen-experienced cytotoxic T cells (CD8+CD45RO+Ki67+) and the proximity of antigen-experienced cytotoxic T cells to melanoma cells were associated with positive response to ICIs. Our study highlights the potential of multiplexed single-cell technology to quantify spatial cell-cell interactions within the tumor microenvironment to understand immune therapy responses.

Inferring the evolutionary history of gene clusters from phylogenetic and gene order data.

Gene duplication is frequent within gene clusters and plays a fundamental role in evolution by providing a source of new genetic material upon which natural selection can act. Although classical phylogenetic inference methods provide some insight into the evolutionary history of a gene cluster, they are not sufficient alone to differentiate single- from multiple gene duplication events and to answer other questions regarding the nature and size of evolutionary events. In this paper, we present an algorithm allowing to infer a set of optimal evolutionary histories for a gene cluster in a single species, according to a general cost model involving variable length duplications (in tandem or inverted), deletions, and inversions. We applied our algorithm to the human olfactory receptor and protocadherin gene clusters, showing that the duplication size distribution differs significantly between the two gene families. The algorithm is available through a web interface at http://www-lbit.iro.umontreal.ca/DILTAG/.

Mutations in the IFNγ-JAK-STAT Pathway Causing Resistance to Immune Checkpoint Inhibitors in Melanoma Increase Sensitivity to Oncolytic Virus Treatment.

PURPOSE: Next-generation sequencing studies and CRISPR-Cas9 screens have established mutations in the IFNγ-JAK-STAT pathway as an immune checkpoint inhibitor (ICI) resistance mechanism in a subset of patients with melanoma. We hypothesized ICI resistance mutations in the IFNγ pathway would simultaneously render melanomas susceptible to oncolytic virus (OV) therapy. EXPERIMENTAL DESIGN: Cytotoxicity experiments were performed with a number of OVs on a matched melanoma cell line pair generated from a baseline biopsy and a progressing lesion with complete JAK2 loss from a patient that relapsed on anti-PD-1 therapy, in melanoma lines following JAK1/2 RNA interference (RNAi) and pharmacologic inhibition and in Jak2 knockout (KO) B16-F10 mouse melanomas. Furthermore, we estimated the frequency of genetic alterations in the IFNγ-JAK-STAT pathway in human melanomas. RESULTS: The melanoma line from an anti-PD-1 progressing lesion was 7- and 22-fold more sensitive to the modified OVs, herpes simplex virus 1 (HSV1-dICP0) and vesicular stomatitis virus (VSV-Δ51), respectively, compared with the line from the baseline biopsy. RNAi, JAK1/2 inhibitor studies, and in vivo studies of Jak2 KOs B16-F10 melanomas revealed a significant increase in VSV-Δ51 sensitivity with JAK/STAT pathway inhibition. Our analysis of The Cancer Genome Atlas data estimated that approximately 11% of ICI-naïve cutaneous melanomas have alterations in IFNγ pathway genes that may confer OV susceptibility. CONCLUSIONS: We provide mechanistic support for the use of OVs as a precision medicine strategy for both salvage therapy in ICI-resistant and first-line treatment in melanomas with IFNγ-JAK-STAT pathway mutations. Our study also supports JAK inhibitor-OV combination therapy for treatment-naïve melanomas without IFN signaling defects.See related commentary by Kaufman, p. 3278.

Multi-omic analysis reveals significantly mutated genes and DDX3X as a sex-specific tumor suppressor in cutaneous melanoma.

The high background tumor mutation burden in cutaneous melanoma limits the ability to identify significantly mutated genes (SMGs) that drive this cancer. To address this, we performed a mutation significance study of over 1,000 melanoma exomes, combined with a multi-omic analysis of 470 cases from The Cancer Genome Atlas. We discovered several SMGs with co-occurring loss-of-heterozygosity and loss-of-function mutations, including PBRM1, PLXNC1 and PRKAR1A, which encodes a protein kinase A holoenzyme subunit. Deconvolution of bulk tumor transcriptomes into cancer, immune and stromal components revealed a melanoma-intrinsic oxidative phosphorylation signature associated with protein kinase A pathway alterations. We also identified SMGs on the X chromosome, including the RNA helicase DDX3X, whose loss-of-function mutations were exclusively observed in males. Finally, we found that tumor mutation burden and immune infiltration contain complementary information on survival of patients with melanoma. In summary, our multi-omic analysis provides insights into melanoma etiology and supports contribution of specific mutations to the sex bias observed in this cancer.

Recurrent somatic BRAF insertion (p.V504_R506dup): a tumor marker and a potential therapeutic target in pilocytic astrocytoma.

Pilocytic astrocytoma (PA) is emerging as a tumor entity with dysregulated RAS/RAF/MEK/ERK signaling. In this study, we report the identification of a novel recurrent BRAF insertion (p.V504_R506dup) in five PA cases harboring exclusively this somatic tandem duplication. This recurrent alteration leads to an addition of three amino acids in the kinase domain of BRAF and has functional impact on activating MAPK phosphorylation. Importantly, we show that this mutation confers resistance to RAF inhibitors without changing effectiveness while downstream MEK inhibitors remain effective. Our results further emphasize the importance of BRAF alterations in PA and the need to characterize them in a given tumor as this can affect therapeutic strategies and their potential use as tumor marker in molecular diagnostics.

The artiodactyl APOBEC3 innate immune repertoire shows evidence for a multi-functional domain organization that existed in the ancestor of placental mammals.

BACKGROUND: APOBEC3 (A3) proteins deaminate DNA cytosines and block the replication of retroviruses and retrotransposons. Each A3 gene encodes a protein with one or two conserved zinc-coordinating motifs (Z1, Z2 or Z3). The presence of one A3 gene in mice (Z2-Z3) and seven in humans, A3A-H (Z1a, Z2a-Z1b, Z2b, Z2c-Z2d, Z2e-Z2f, Z2g-Z1c, Z3), suggests extraordinary evolutionary flexibility. To gain insights into the mechanism and timing of A3 gene expansion and into the functional modularity of these genes, we analyzed the genomic sequences, expressed cDNAs and activities of the full A3 repertoire of three artiodactyl lineages: sheep, cattle and pigs. RESULTS: Sheep and cattle have three A3 genes, A3Z1, A3Z2 and A3Z3, whereas pigs only have two, A3Z2 and A3Z3. A comparison between domestic and wild pigs indicated that A3Z1 was deleted in the pig lineage. In all three species, read-through transcription and alternative splicing also produced a catalytically active double domain A3Z2-Z3 protein that had a distinct cytoplasmic localization. Thus, the three A3 genes of sheep and cattle encode four conserved and active proteins. These data, together with phylogenetic analyses, indicated that a similar, functionally modular A3 repertoire existed in the common ancestor of artiodactyls and primates (i.e., the ancestor of placental mammals). This mammalian ancestor therefore possessed the minimal A3 gene set, Z1-Z2-Z3, required to evolve through a remarkable series of eight recombination events into the present day eleven Z domain human repertoire. CONCLUSION: The dynamic recombination-filled history of the mammalian A3 genes is consistent with the modular nature of the locus and a model in which most of these events (especially the expansions) were selected by ancient pathogenic retrovirus infections.

Dual MAPK Inhibition Is an Effective Therapeutic Strategy for a Subset of Class II BRAF Mutant Melanomas.

PURPOSE: Dual MAPK pathway inhibition (dMAPKi) with BRAF and MEK inhibitors improves survival in BRAF V600E/K mutant melanoma, but the efficacy of dMAPKi in non-V600 BRAF mutant tumors is poorly understood. We sought to characterize the responsiveness of class II (enhanced kinase activity, dimerization dependent) BRAF mutant melanoma to dMAPKi. EXPERIMENTAL DESIGN: Tumors from patients with BRAF wild-type (WT), V600E (class I), and L597S (class II) metastatic melanoma were used to generate patient-derived xenografts (PDX). We assembled a panel of melanoma cell lines with class IIa (activation segment) or IIb (p-loop) mutations and compared these with WT or V600E/K BRAF mutant cells. Cell lines and PDXs were treated with BRAFi (vemurafenib, dabrafenib, encorafenib, and LY3009120), MEKi (cobimetinib, trametinib, and binimetinib), or the combination. We identified 2 patients with BRAF L597S metastatic melanoma who were treated with dMAPKi. RESULTS: BRAFi impaired MAPK signaling and cell growth in class I and II BRAF mutant cells. dMAPKi was more effective than either single MAPKi at inhibiting cell growth in all class II BRAF mutant cells tested. dMAPKi caused tumor regression in two melanoma PDXs with class II BRAF mutations and prolonged survival of mice with class II BRAF mutant melanoma brain metastases. Two patients with BRAF L597S mutant melanoma clinically responded to dMAPKi. CONCLUSIONS: Class II BRAF mutant melanoma is growth inhibited by dMAPKi. Responses to dMAPKi have been observed in 2 patients with class II BRAF mutant melanoma. These data provide rationale for clinical investigation of dMAPKi in patients with class II BRAF mutant metastatic melanoma.See related commentary by Johnson and Dahlman, p. 6107.

Duplication and inversion history of a tandemly repeated genes family.

Given a phylogenetic tree for a family of tandemly repeated genes and their signed order on the chromosome, we aim to find the minimum number of inversions compatible with an evolutionary history of this family. This is the first attempt to account for inversions in an evolutionary model of tandemly repeated genes. We present a branch-and-bound algorithm that finds the exact solution, and a polynomial-time heuristic based on the breakpoint distance. We show, on simulated data, that those algorithms can be used to improve phylogenetic inference of tandemly repeated gene families. An application on a published phylogeny of KRAB zinc finger genes is presented.

Characterization of the microDNA through the response to chemotherapeutics in lymphoblastoid cell lines.

Recently, a new class of extrachromosomal circular DNA, called microDNA, was identified. They are on average 100 to 400 bp long and are derived from unique non-repetitive genomic regions with high gene density. MicroDNAs are thought to arise from DNA breaks associated with RNA metabolism or replication slippage. Given the paucity of information on this entirely novel phenomenon, we aimed to get an additional insight into microDNA features by performing the microDNA analysis in 20 independent human lymphoblastoid cell lines (LCLs) prior and after treatment with chemotherapeutic drugs. The results showed non-random genesis of microDNA clusters from the active regions of the genome. The size periodicity of 190 bp was observed, which matches DNA fragmentation typical for apoptotic cells. The chemotherapeutic drug-induced apoptosis of LCLs increased both number and size of clusters further suggesting that part of microDNAs could result from the programmed cell death. Interestingly, proportion of identified microDNA sequences has common loci of origin when compared between cell line experiments. While compatible with the original observation that microDNAs originate from a normal physiological process, obtained results imply complementary source of its production. Furthermore, non-random genesis of microDNAs depicted by redundancy between samples makes these entities possible candidates for new biomarker generation.

A childhood acute lymphoblastic leukemia-specific lncRNA implicated in prednisolone resistance, cell proliferation, and migration.

Childhood acute lymphoblastic leukemia (cALL) is the most common pediatric cancer and, despite an 85% cure rate, still represents a major cause of disease-related death in children. Recent studies have implicated long non-coding RNAs (lncRNAs) in cALL etiology, progression, and treatment response. However, barring some exceptions little is known about the functional impact of lncRNAs on cancer biology, which limits their potential as potential therapeutic targets. We wanted to investigate the functional role of lncRNAs identified as specifically overexpressed in pre-B cALL by whole-transcriptome sequencing. Here we report five lncRNAs specifically upregulated in pre-B cALL that had significant impacts on cancer hallmark traits such as cell proliferation, migration, apoptosis, and treatment response. In particular, silencing of the RP11-137H2.4 lncRNA effectively restored normal glucocorticoid (GC) response in a GC-resistant pre-B cALL cell line and specifically modulated expression of members of both the NRAS/BRAF/NF-κB MAPK cascade and cell cycle pathways. Since GC form the cornerstone of cALL chemotherapy and resistance in cALL confers a dismal prognosis, characterizing RP11-137H2.4'sexact role and function in this process will be critical to the development of new therapeutic approaches to overcome GC resistance in children treated for cALL.

Specific expression of novel long non-coding RNAs in high-hyperdiploid childhood acute lymphoblastic leukemia.

Pre-B cell childhood acute lymphoblastic leukemia (pre-B cALL) is a heterogeneous disease involving many subtypes typically stratified using a combination of cytogenetic and molecular-based assays. These methods, although widely used, rely on the presence of known chromosomal translocations, which is a limiting factor. There is therefore a need for robust, sensitive, and specific molecular biomarkers unaffected by such limitations that would allow better risk stratification and consequently better clinical outcome. In this study we performed a transcriptome analysis of 56 pre-B cALL patients to identify expression signatures in different subtypes. In both protein-coding and long non-coding RNAs (lncRNA), we identified subtype-specific gene signatures distinguishing pre-B cALL subtypes, particularly in t(12;21) and hyperdiploid cases. The genes up-regulated in pre-B cALL subtypes were enriched in bivalent chromatin marks in their promoters. LncRNAs is a new and under-studied class of transcripts. The subtype-specific nature of lncRNAs suggests they may be suitable clinical biomarkers to guide risk stratification and targeted therapies in pre-B cALL patients.

Recovering haplotype structure through recombination and gene conversion.

MOTIVATION: Understanding haplotype evolution subject to mutation, recombination and gene conversion is fundamental to understand genetic specificities of human populations and hereditary bases of complex disorders. The goal of this project is to develop new algorithmic tools assisting the reconstruction of historical relationships between haplotypes and the inference of haplotypes from genotypes. RESULTS: We present two new algorithms. The first one finds an optimal pathway of mutations, recombinations and gene conversions leading to a given haplotype of size m from a population of h haplotypes. It runs in time O(mhs(2)), where s is the maximum number of contiguous sites that can be exchanged in a single gene conversion. The second one finds an optimal pathway of mutations and recombinations leading to a given genotype, and runs in time O(mh(2)). Both algorithms are based on a penalty score model and use a dynamic programming approach. We apply the second one to the problem of inferring haplotypes from genotypes, and show how it can be used as an independent tool, or to improve the performance of existing methods. AVAILABILITY: The algorithms have been implemented in JAVA and are available on request.