Aberrant MNX1 expression associated with t(7;12)(q36;p13) pediatric acute myeloid leukemia induces the disease through altering histone methylation

Certain subtypes of acute myeloid leukemia (AML) in children have inferior outcome, such as AML with translocation t(7;12)(q36;p13) leading to an MNX1::ETV6 fusion along with high expression of MNX1. We have identified the transforming event in this AML and possible ways of treatment. Retroviral expression of MNX1 was able to induce AML in mice, with similar gene expression and pathway enrichment to t(7;12) AML patient data. Importantly, this leukemia was only induced in immune incompetent mice using fetal but not adult hematopoietic stem and progenitor cells. The restriction in transforming capacity to cells from fetal liver is in alignment with t(7;12)(q36;p13) AML being mostly seen in infants. Expression of MNX1 led to increased histone 3 lysine 4 mono-, di- and trimethylation, reduction in H3K27me3, accompanied with changes in genome-wide chromatin accessibility and genome expression, likely mediated through MNX1 interaction with the methionine cycle and methyltransferases. MNX1 expression increased DNA damage, depletion of the Lin-/Sca1+/c-Kit+ population and skewing toward the myeloid lineage. These effects, together with leukemia development, were prevented by pre-treatment with the S-adenosylmethionine analog Sinefungin. In conclusion, we have shown the importance of MNX1 in development of AML with t(7;12), supporting a rationale for targeting MNX1 and downstream pathways.

DNMT and HDAC inhibition induces immunogenic neoantigens from human endogenous retroviral element-derived transcripts

Immunotherapies targeting cancer-specific neoantigens have revolutionized the treatment of cancer patients. Recent evidence suggests that epigenetic therapies synergize with immunotherapies, mediated by the de-repression of endogenous retroviral element (ERV)-encoded promoters, and the initiation of transcription. Here, we use deep RNA sequencing from cancer cell lines treated with DNA methyltransferase inhibitor (DNMTi) and/or Histone deacetylase inhibitor (HDACi), to assemble a de novo transcriptome and identify several thousand ERV-derived, treatment-induced novel polyadenylated transcripts (TINPATs). Using immunopeptidomics, we demonstrate the human leukocyte antigen (HLA) presentation of 45 spectra-validated treatment-induced neopeptides (t-neopeptides) arising from TINPATs. We illustrate the potential of the identified t-neopeptides to elicit a T-cell response to effectively target cancer cells. We further verify the presence of t-neopeptides in AML patient samples after in vivo treatment with the DNMT inhibitor Decitabine. Our findings highlight the potential of ERV-derived neoantigens in epigenetic and immune therapies.

High-resolution transcriptomic and epigenetic profiling identifies novel regulators of COPD

Patients with chronic obstructive pulmonary disease (COPD) are still waiting for curative treatments. Considering its environmental cause, we hypothesized that COPD will be associated with altered epigenetic signaling in lung cells. We generated genome-wide DNA methylation maps at single CpG resolution of primary human lung fibroblasts (HLFs) across COPD stages. We show that the epigenetic landscape is changed early in COPD, with DNA methylation changes occurring predominantly in regulatory regions. RNA sequencing of matched fibroblasts demonstrated dysregulation of genes involved in proliferation, DNA repair, and extracellular matrix organization. Data integration identified 110 candidate regulators of disease phenotypes that were linked to fibroblast repair processes using phenotypic screens. Our study provides high-resolution multi-omic maps of HLFs across COPD stages. We reveal novel transcriptomic and epigenetic signatures associated with COPD onset and progression and identify new candidate regulators involved in the pathogenesis of chronic lung diseases. The presence of various epigenetic factors among the candidates demonstrates that epigenetic regulation in COPD is an exciting research field that holds promise for novel therapeutic avenues for patients.

DNA methylation landscape of tumor-associated macrophages reveals pathways, transcription factors and prognostic value relevant to triple-negative breast cancer patients

The accumulation of myeloid cells, particularly tumor-associated macrophages (TAMs), characterizes the tumor microenvironment (TME) of many solid cancers, including breast cancer. Compared to healthy tissue-resident macrophages, TAMs acquire distinct transcriptomes and tumor-promoting functions by largely unknown mechanisms. Here, we hypothesize the involvement of TME signaling and subsequent epigenetic reprogramming of TAMs. Using the 4T1 mouse model of triple-negative breast cancer, we demonstrate that the presence of cancer cells significantly alters the DNA methylation landscape of macrophages and, to a lesser extent, bone marrow-derived monocytes (BMDMs). TAM methylomes, dissected into BMDM-originating and TAM-specific epigenetic programs, implicated transcription factors (TFs) and signaling pathways involved in TAM reprogramming, correlated with cancer-specific gene expression patterns. Utilizing published single-cell gene expression data, we linked microenvironmentally-derived signals to the cancer-specific DNA methylation landscape of TAMs. These integrative analyses highlighted the role of altered cytokine production in the TME (eg, TGF-β, IFN-γ and CSF1) on the induction of specific TFs (eg, FOSL2, STAT1 and RUNX3) responsible for the epigenetic reprogramming of TAMs. DNA methylation deconvolution identified a TAM-specific signature associated with the identified signaling pathways and TFs, corresponding with severe tumor grade and poor prognosis of breast cancer patients. Similarly, immunosuppressive TAM functions were identified, such as induction of the immune inhibitory receptor-ligand PD-L1 by DNA hypomethylation of Cd274. Collectively, these results provide strong evidence that the epigenetic landscapes of macrophages and monocytes are perturbed by the presence of breast cancer, pointing to molecular mechanisms of TAM reprogramming, impacting patient outcomes.

Global DNA Methylation Analysis of Cancer-Associated Fibroblasts Reveals Extensive Epigenetic Rewiring Linked with RUNX1 Upregulation in Breast Cancer Stroma

Cancer cells recruit and rewire normal fibroblasts in their microenvironment to become protumorigenic cancer-associated fibroblasts (CAF). These CAFs are genomically stable, yet their transcriptional programs are distinct from those of their normal counterparts. Transcriptional regulation plays a major role in this reprogramming, but the extent to which epigenetic modifications of DNA also contribute to the rewiring of CAF transcription is not clear. Here we address this question by dissecting the epigenetic landscape of breast CAFs. Applying tagmentation-based whole-genome bisulfite sequencing in a mouse model of breast cancer, we found that fibroblasts undergo massive DNA methylation changes as they transition into CAFs. Transcriptional and epigenetic analyses revealed RUNX1 as a potential mediator of this process and identified a RUNX1-dependent stromal gene signature. Coculture and mouse models showed that both RUNX1 and its stromal signature are induced as normal fibroblasts transition into CAFs. In breast cancer patients, RUNX1 was upregulated in CAFs, and expression of the RUNX1 signature was associated with poor disease outcome, highlighting the relevance of these findings to human disease. This work presents a comprehensive genome-wide map of DNA methylation in CAFs and reveals a previously unknown facet of the dynamic plasticity of the stroma.

SIGNIFICANCE

The first genome-wide map of DNA methylation in breast cancer-associated fibroblasts unravels a previously unknown facet of the dynamic plasticity of the stroma, with far-reaching therapeutic implications.

Transcriptional regulation of the thymus master regulator Foxn1

FOXN1 is a transcription factor critical for the development of both thymic epithelial cell (TEC) and hair follicle cell (HFC) compartments. However, mechanisms controlling its expression remain poorly understood. To address this question, we performed thorough analyses of the evolutionary conservation and chromatin status of the Foxn1 locus in different tissues and states and identified several putative cis-regulatory regions unique to TECs versus HFCs. Furthermore, experiments using genetically modified mice with specific deletions in the Foxn1 locus and additional bioinformatic analyses helped us identify key regions and transcription factors involved in either positive or negative regulation of Foxn1 in both TECs and HFCs. Specifically, we identified SIX1 and FOXN1 itself as key factors inducing Foxn1 expression in embryonic and neonatal TECs. Together, our data provide important mechanistic insights into the transcriptional regulation of the Foxn1 gene in TEC versus HFC and highlight the role of FOXN1 in its autoregulation.

The antileukemic activity of decitabine upon PML/RARA-negative AML blasts is supported by all-trans retinoic acid: in vitro and in vivo evidence for cooperation

The prognosis of AML patients with adverse genetics, such as a complex, monosomal karyotype and TP53 lesions, is still dismal even with standard chemotherapy. DNA-hypomethylating agent monotherapy induces an encouraging response rate in these patients. When combined with decitabine (DAC), all-trans retinoic acid (ATRA) resulted in an improved response rate and longer overall survival in a randomized phase II trial (DECIDER; NCT00867672). The molecular mechanisms governing this in vivo synergism are unclear. We now demonstrate cooperative antileukemic effects of DAC and ATRA on AML cell lines U937 and MOLM-13. By RNA-sequencing, derepression of >1200 commonly regulated transcripts following the dual treatment was observed. Overall chromatin accessibility (interrogated by ATAC-seq) and, in particular, at motifs of retinoic acid response elements were affected by both single-agent DAC and ATRA, and enhanced by the dual treatment. Cooperativity regarding transcriptional induction and chromatin remodeling was demonstrated by interrogating the HIC1, CYP26A1, GBP4, and LYZ genes, in vivo gene derepression by expression studies on peripheral blood blasts from AML patients receiving DAC + ATRA. The two drugs also cooperated in derepression of transposable elements, more effectively in U937 (mutated TP53) than MOLM-13 (intact TP53), resulting in a “viral mimicry” response. In conclusion, we demonstrate that in vitro and in vivo, the antileukemic and gene-derepressive epigenetic activity of DAC is enhanced by ATRA.

Club cells employ regeneration mechanisms during lung tumorigenesis

The high plasticity of lung epithelial cells, has for many years, confounded the correct identification of the cell-of-origin of lung adenocarcinoma (LUAD), one of the deadliest malignancies worldwide. Here, we employ lineage-tracing mouse models to investigate the cell of origin of Eml4-Alk LUAD, and show that Club and Alveolar type 2 (AT2) cells give rise to tumours. We focus on Club cell originated tumours and find that Club cells experience an epigenetic switch by which they lose their lineage fidelity and gain an AT2-like phenotype after oncogenic transformation. Single-cell transcriptomic analyses identified two trajectories of Club cell evolution which are similar to the ones used during lung regeneration, suggesting that lung epithelial cells leverage on their plasticity and intrinsic regeneration mechanisms to give rise to a tumour. Together, this study highlights the role of Club cells in LUAD initiation, identifies the mechanism of Club cell lineage infidelity, confirms the presence of these features in human tumours, and unveils key mechanisms conferring LUAD heterogeneity.

ID3 promotes homologous recombination via non-transcriptional and transcriptional mechanisms and its loss confers sensitivity to PARP inhibition

The inhibitor of DNA-binding 3 (ID3) is a transcriptional regulator that limits interaction of basic helix-loop-helix transcription factors with their target DNA sequences. We previously reported that ID3 loss is associated with mutational signatures linked to DNA repair defects. Here we demonstrate that ID3 exhibits a dual role to promote DNA double-strand break (DSB) repair, particularly homologous recombination (HR). ID3 interacts with the MRN complex and RECQL helicase to activate DSB repair and it facilitates RAD51 loading and downstream steps of HR. In addition, ID3 promotes the expression of HR genes in response to ionizing radiation by regulating both chromatin accessibility and activity of the transcription factor E2F1. Consistently, analyses of TCGA cancer patient data demonstrate that low ID3 expression is associated with impaired HR. The loss of ID3 leads to sensitivity of tumor cells to PARP inhibition, offering new therapeutic opportunities in ID3-deficient tumors.

Epigenetic reprogramming of airway macrophages promotes polarization and inflammation in muco-obstructive lung disease

Lung diseases, such as cystic fibrosis and COPD, are characterized by mucus obstruction and chronic airway inflammation, but their mechanistic link remains poorly understood. Here, we focus on the function of the mucostatic airway microenvironment on epigenetic reprogramming of airway macrophages (AM) and resulting transcriptomic and phenotypical changes. Using a mouse model of muco-obstructive lung disease (Scnn1b-transgenic), we identify epigenetically controlled, differentially regulated pathways and transcription factors involved in inflammatory responses and macrophage polarization. Functionally, AMs from Scnn1b-transgenic mice have reduced efferocytosis and phagocytosis, and excessive inflammatory responses upon lipopolysaccharide challenge, mediated through enhanced Irf1 function and expression. Ex vivo stimulation of wild-type AMs with native mucus impairs efferocytosis and phagocytosis capacities. In addition, mucus induces gene expression changes, comparable with those observed in AMs from Scnn1b-transgenic mice. Our data show that mucostasis induces epigenetic reprogramming of AMs, leading to changes favoring tissue damage and disease progression. Targeting these altered AMs may support therapeutic approaches in patients with muco-obstructive lung diseases.

The age of the bone marrow microenvironment influences B-cell acute lymphoblastic leukemia progression via CXCR5-CXCL13

B-cell acute lymphoblastic leukemia (B-ALL) occurs most commonly in children, whereas chronic myeloid leukemia is more frequent in adults. The myeloid bias of hematopoiesis in elderly individuals has been considered causative, but the age of the bone marrow microenvironment (BMM) may be contributory. Using various murine models of B-ALL in young vs old mice, we recapitulated B-ALL preponderance in children vs adults. We showed differential effects of young vs old BM macrophages on B-ALL cell function. Molecular profiling using RNA- and ATAC-sequencing revealed pronounced differences in young vs old BMM-derived macrophages and enrichment for gene sets associated with inflammation. In concordance with the role of C-X-C motif chemokine (CXCL) 13 for disease-associated B-cell chemoattraction, we found CXCL13 to be highly expressed in young macrophages on a translational compared with a transcriptional level. Inhibition of CXCL13 in BM macrophages impaired leukemia cell migration and decreased the proliferation of cocultured B-ALL cells, whereas recombinant CXCL13 increased pAKT and B-ALL cell expansion. Pretreatment of B-ALL-initiating cells with CXCL13 accelerated B-ALL progression. Deficiency of Cxcr5, the receptor for CXCL13, on B-ALL-initiating cells prolonged murine survival, whereas high expression of CXCR5 in pediatric B-ALL may predict central nervous system relapse. CXCL13 staining was increased in bone sections from pediatric compared with adult patients with B-ALL. Taken together, our study shows that the age of the BMM and, in particular, BM macrophages influence the leukemia phenotype. The CXCR5-CXCL13 axis may act as prognostic marker and an attractive novel target for the treatment of B-ALL.

10-Year Clinical Comparative Study of Ceramic and Composite Veneered Metal Crowns

The rejection of composite veneerings in fixed partial dentures is primarily caused by the inadequacy of the bonds between the metal framework and the composite veneering. The development of improved veneering composites necessitates an investigation of their clinical suitability compared with ceramic materials. Nineteen patients with at least two suitable, adjacent natural teeth for crowning were treated with 64 high noble alloy crowns. The adjacent crowns were veneered with ceramic (IPS inline) and composite materials (SR Adoro). Seven follow-up examinations were carried out over a period of 10 years. The crowns were investigated for mechanical defects, periodontal parameters, and discolorations. The survival rates of the ceramic veneered metal crowns (CeMCs) and composite veneered metal crowns (CoMCs) at the 10-year follow-up were 87.1 and 87.9%, respectively. The success rates of the crowns after 10 years were 83.9% for CeMCs and 51.5% for CoMCs (log-rank test, p = 0.009). No significant differences between the groups were found in the periodontal parameters (Kruskal-Wallis one-way analysis of variance, ANOVA p ⟩ 0.05). After 10 years, discoloration patterns of the two materials differed significantly (Mann-Whitney-U-test, p = 0.017). Thus, despite the improvements associated with CoMCs, CeMCs remain the gold standard for veneered metal crowns.

Epigenetic analysis of patients with T-ALL identifies poor outcomes and a hypomethylating agent-responsive subgroup

Adult "T cell" acute lymphoblastic leukemia (T-ALL) is an aggressive hematological malignancy that is associated with poor outcomes, requiring additional therapeutic options. The DNA methylation landscapes of adult T-ALL remain undercharacterized. Here, we systematically analyzed the DNA methylation profiles of normal thymic-sorted T cell subpopulations and 143 primary adult T-ALLs as part of the French GRAALL 2003-2005 trial. Our results indicated that T-ALL is epigenetically heterogeneous consisting of five subtypes (C₁-C₅), which were either associated with co-occurring DNA methyltransferase 3 alpha (DNMT3A)/isocitrate dehydrogenase [NADP(+)] 2 (IDH2) mutations (C₁), TAL bHLH transcription factor 1, erythroid differentiation factor (TAL1) deregulation (C₂), T cell leukemia homeobox 3 (TLX3) (C₃), TLX1/in cis-homeobox A9 (HOXA9) (C₄), or in trans-HOXA9 overexpression (C₅). Integrative analysis of DNA methylation and gene expression identified potential cluster-specific oncogenes and tumor suppressor genes. In addition to an aggressive hypomethylated subgroup (C₁), our data identified an unexpected subset of hypermethylated T-ALL (C₅) associated with poor outcome and primary therapeutic response. Using mouse xenografts, we demonstrated that hypermethylated T-ALL samples exhibited therapeutic responses to the DNA hypomethylating agent 5-azacytidine, which significantly (survival probability; P = 0.001 for C₃, 0.01 for C₄, and 0.0253 for C₅) delayed tumor progression. These findings suggest that epigenetic-based therapies may provide an alternative treatment option in hypermethylated T-ALL.

Methrix: an R/bioconductor package for systematic aggregation and analysis of bisulfite sequencing data

MOTIVATION

Whole genome bisulfite sequencing (WGBS), measures DNA methylation at base pair resolution resulting in large bedGraph like coverage files. Current options for processing such files are hindered by discrepancies in file format specification, speed and memory requirements.

RESULTS

We developed methrix, an R package, which provides a toolset for systematic analysis of large datasets. Core functionality of the package includes a comprehensive bedGraph or similar tab-separated text file reader - which summarizes methylation calls based on annotated reference indices, infers and collapses strands, and handles uncovered reference CpG sites while facilitating a flexible input file format specification. Additional optimized functions for quality control filtering, sub-setting, and visualization allow user-friendly and effective processing of WGBS results. Easy integration with tools for differentially methylated region (DMR) calling and annotation further eases the analysis of genome-wide methylation data. Overall, methrix enriches established WGBS workflows by bringing together computational efficiency and versatile functionality.

AVAILABILITY AND IMPLEMENTATION

Methrix is implemented as an R package, made available under MIT license at https://github.com/CompEpigen/methrix and can be installed from the Bioconductor repository.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

AmpliconDesign - an interactive web server for the design of high-throughput targeted DNA methylation assays

Targeted analysis of DNA methylation patterns based on bisulfite-treated genomic DNA (BT-DNA) is considered as a gold-standard for epigenetic biomarker development. Existing software tools facilitate primer design, primer quality control or visualization of primer localization. However, high-throughput design of primers for BT-DNA amplification is hampered by limits in throughput and functionality of existing tools, requiring users to repeatedly perform specific tasks manually. Consequently, the design of PCR primers for BT-DNA remains a tedious and time-consuming process. To bridge this gap, we developed AmpliconDesign, a webserver providing a scalable and user-friendly platform for the design and analysis of targeted DNA methylation studies based on BT-DNA, e.g. deep amplicon bisulfite sequencing (ampBS-seq) or EpiTYPER MassArray. Core functionality of the web server includes high-throughput primer design and binding site validation based on in silico bisulfite-converted DNA sequences, prediction of fragmentation patterns for EpiTYPER MassArray, an interactive quality control as well as a streamlined analysis workflow for ampBS-seq.

Globally altered epigenetic landscape and delayed osteogenic differentiation in H3.3-G34W-mutant giant cell tumor of bone

The neoplastic stromal cells of giant cell tumor of bone (GCTB) carry a mutation in H3F3A, leading to a mutant histone variant, H3.3-G34W, as a sole recurrent genetic alteration. We show that in patient-derived stromal cells H3.3-G34W is incorporated into the chromatin and associates with massive epigenetic alterations on the DNA methylation, chromatin accessibility and histone modification level, that can be partially recapitulated in an orthogonal cell line system by the introduction of H3.3-G34W. These epigenetic alterations affect mainly heterochromatic and bivalent regions and provide possible explanations for the genomic instability, as well as the osteolytic phenotype of GCTB. The mutation occurs in differentiating mesenchymal stem cells and associates with an impaired osteogenic differentiation. We propose that the observed epigenetic alterations reflect distinct differentiation stages of H3.3 WT and H3.3 MUT stromal cells and add to H3.3-G34W-associated changes.

The histone variant mutation H3.3-G34W occurs in the majority of giant cell tumor of bone (GCTB). By profiling patient-derived GCTB tumor cells, the authors show that this mutation associates with epigenetic alterations in heterochromatic and bivalent regions that contribute to an impaired osteogenic differentiation and the osteolytic phenotype of GCTB.

Implementing a new curriculum for computer-assisted restorations in prosthetic dentistry

INTRODUCTION

Computer-aided design/computer-aided manufacturing (CAD/CAM) of fixed prosthetic restorations has gained popularity in the last decade. However, this field of dentistry has not been integrated in the dental curriculum at most universities.

MATERIAL AND METHODS

According to the method of Kern, a curriculum was designed and established on a voluntary basis in the prosthetic education of a German dental school. The success of the implementation was measured by evaluation carried out by the participants on a visual analogue scale. Furthermore, the clinical performance of the fabricated restorations was evaluated.

RESULTS

Ninety-four percent of all students participated in the CAD/CAM curriculum indicating considerable interest. Nearly half of all students used the acquired knowledge to design crowns for their patients. All restorations fabricated by participants of the new CAD/CAM programme showed good clinical performance.

DISCUSSION

By phasing-in the CAD/CAM training programme, independent CAD/CAM-based fabrication of all-ceramic crowns increased student's self-confidence in tooth preparation. A tendency was found that students using CAD/CAM technology prepared more teeth than their fellow students who did not use CAD/CAM technology. Further studies are required to investigate the influence of independent CAD/CAM-based single-crown fabrication on the quality of the preparation.

Performance of chromID® CARBA medium for carbapenemases-producing Enterobacteriaceae detection during rectal screening

Chromogenic chromID® CARBA medium was compared with CDC method and MacConkey agar with imipenem for its performance in detecting carbapenemase-producing Enterobacteriaceae (CPE) during a faecal screening surveillance program. Double rectal swabs were collected from patients hospitalized in the ICU. One swab was inoculated onto the solid media chromID® CARBA and MacConkey agar with imipenem, while the other was tested according to CDC protocol. Suspected colonies from all procedures were identified to species level and tested for their susceptibility to carbapenems by phenotypic tests. All carbapenem non-susceptible isolates were tested by the modified Hodge test (MHT) and synergy tests. Positive results were confirmed by PCR testing for carbapenemase gene detection. Performance of all three procedures applied was statistically analyzed as compared to MHT and PCR results for the presence of carbapenemase-encoding genes. Out of 177 rectal samples tested, 86 samples were found to contain one or more CPE verified by molecular detection of carbapenemase encoding genes among isolated Enterobacteriaceae. Sensitivity of chromID® CARBA was 96.5 % in clinical samples. Specificity was 91.2 % at the reading level and 100.0 % after Gram staining. chromID® CARBA performed with high accuracy among the phenotypic methods applied, giving early results.

[Radiation-related caries: etiology and possible preventive strategies. What should the radiotherapist know?]

BACKGROUND

After radiation treatment of head-and-neck cancer, the impairment of patient's quality of life still remains an issue. After completion of the treatment course, a substantial number of patients develop so-called radiation caries. In addition, almost 50% of all cases of infectious osteoradionecrosis (iORN) of the jaws are directly associated with radiation caries. This review addresses our current knowledge on the etiology and pathogenesis of radiation caries including possible preventive strategies.

MATERIALS AND METHODS

A PubMed search using the terms "radiation caries" ("radiation related caries", "radiation related damage to dentition") and "radiogenic caries" ("postradiation caries", "dental complications and radiotherapy") was performed. The analysis of its content focused on the etiology, the pathogenesis, and the available knowledge on prophylaxis as well as treatment of radiation caries.

RESULTS

For this review, 60 publications were selected. As main causal factors for radiogenic caries, either indirect impairment, resulting from alterations in the oral environment (e.g., radiation-induced xerostomia) or direct radiation-induced damage in teeth hard tissues are discussed. Radiation caries remains a lifelong threat and, therefore, requires permanent prevention programs.

CONCLUSION

To enable optimal medical care of the patients during the time course of radiotherapy as well as afterwards, close interdisciplinary cooperation between radiotherapists, oral surgeons, otorhinolaryngologists, and dentists is absolutely essential.

The mind of the species problem

The species problem is the long-standing failure of biologists to agree on how we should identify species and how we should define the word 'species'. The innumerable attacks on the problem have turned the often-repeated question 'what are species?' into a philosophical conundrum. Today, the preferred form of attack is the well-crafted argument, and debaters seem to have stopped inquiring about what new information is needed to solve the problem. However, our knowledge is not complete and we have overlooked something. The species problem can be overcome if we understand our own role, as conflicted investigators, in causing the problem.

Human populations show reduced DNA sequence variation at the Factor IX locus

Levels and patterns of human DNA sequence variation vary widely among loci. However, some of this variation may be due to the different populations used in different studies. So far, few studies of diverse human populations have compared different genetic loci for the same samples of populations and individuals. Here, we present new polymorphism data from intron 4 of the Factor IX gene (FIX) sequenced in diverse Old World populations. An explicit comparison is made with another X-linked gene, PDHA1, for which the sampling of individuals was very similar. Despite having a similar amount of divergence from chimpanzees, as do other nuclear genes, FIX has comparatively much less DNA sequence variation among humans. Nucleotide diversity at FIX is the lowest among the existing non-Y chromosome nuclear gene datasets and is less than 10% of the diversity found at PDHA1. Estimates of effective population size based on FIX are 8,558, about half of the value obtained for PDHA1, and the time to the most recent common ancestry among human FIX gene copies (282,000 years) is one of the most recent estimates reported for human genes. Analyses presented here suggest a history for the FIX region that includes recent positive directional selection, or background, selection. The general conclusion emerging is that very large variations can exist between the histories of similar genomic regions, even when sampling differences are minimized.

The population genetics of the origin and divergence of the Drosophila simulans complex species

The origins and divergence of Drosophila simulans and close relatives D. mauritiana and D. sechellia were examined using the patterns of DNA sequence variation found within and between species at 14 different genes. D. sechellia consistently revealed low levels of polymorphism, and genes from D. sechellia have accumulated mutations at a rate that is approximately 50% higher than the same genes from D. simulans. At synonymous sites, D. sechellia has experienced a significant excess of unpreferred codon substitutions. Together these observations suggest that D. sechellia has had a reduced effective population size for some time, and that it is accumulating slightly deleterious mutations as a result. D. simulans and D. mauritiana are both highly polymorphic and the two species share many polymorphisms, probably since the time of common ancestry. A simple isolation speciation model, with zero gene flow following incipient species separation, was fitted to both the simulans/mauritiana divergence and the simulans/sechellia divergence. In both cases the model fit the data quite well, and the analyses revealed little evidence of gene flow between the species. The exception is one gene copy at one locus in D. sechellia, which closely resembled other D. simulans sequences. The overall picture is of two allopatric speciation events that occurred quite near one another in time.

Parasite populations: the puzzle of Plasmodium

The issue of whether the malaria parasite Plasmodium falciparum is effectively clonal, as some argue, or undergoes outcrossing at a high rate, as many others believe, has been controversial. Recent data support the latter view, though no doubt the puzzle has not yet been laid to rest.

The limits of selection during maize domestication

The domestication of all major crop plants occurred during a brief period in human history about 10,000 years ago. During this time, ancient agriculturalists selected seed of preferred forms and culled out seed of undesirable types to produce each subsequent generation. Consequently, favoured alleles at genes controlling traits of interest increased in frequency, ultimately reaching fixation. When selection is strong, domestication has the potential to drastically reduce genetic diversity in a crop. To understand the impact of selection during maize domestication, we examined nucleotide polymorphism in teosinte branched1, a gene involved in maize evolution. Here we show that the effects of selection were limited to the gene's regulatory region and cannot be detected in the protein-coding region. Although selection was apparently strong, high rates of recombination and a prolonged domestication period probably limited its effects. Our results help to explain why maize is such a variable crop. They also suggest that maize domestication required hundreds of years, and confirm previous evidence that maize was domesticated from Balsas teosinte of southwestern Mexico.

X chromosome evidence for ancient human histories

Diverse African and non-African samples of the X-linked PDHA1 (pyruvate dehydrogenase E1 alpha subunit) locus revealed a fixed DNA sequence difference between the two sample groups. The age of onset of population subdivision appears to be about 200 thousand years ago. This predates the earliest modern human fossils, suggesting the transformation to modern humans occurred in a subdivided population. The base of the PDHA1 gene tree is relatively ancient, with an estimated age of 1.86 million years, a late Pliocene time associated with early species of Homo. PDHA1 revealed very low variation among non-Africans, but in other respects the data are consistent with reports from other X-linked and autosomal haplotype data sets. Like these other genes, but in conflict with microsatellite and mitochondrial data, PDHA1 does not show evidence of human population expansion.

4-[(1H-imidazol-4-yl) methyl] benzamidines and benzylamidines: novel antagonists of the histamine H3 receptor

A series of amidine substituted phenyl-, benzyl-, and phenethylimidazoles based on the known H3 agonist SK&F 91606 (4) has been synthesized and tested as ligands for the histamine H3 receptor. Insertion of a phenyl ring between the imidazole ring and the amidine moiety produces antagonists. The benzyl series was found to be the most potent and was further investigated. Compounds 9c and 18 (entries 5 and 12, Table 1) are potent ligands for the H3 receptor with K(i) values of 16 nM and 7.2 nM respectively. In vivo, both compounds were shown to be equipotent to thioperamide (2), the standard H3 antagonist.

Selfish genes, pleiotropy and the origin of recombination

If multiple linked polymorphisms are under natural selection, then conflicts arise and the efficiency of natural selection is hindered relative to the case of no linkage. This simple interaction between linkage and natural selection creates an opportunity for mutations that raise the level of recombination to increase in frequency and have an enhanced chance of fixation. This important finding by S. Otto and N. Barton means that mutations that raise the recombination rate, but are otherwise neutral, will be selectively favored under fairly general circumstances of multilocus selection and linkage. The effect described by Otto and Barton, which was limited to neutral modifiers, can also be extended to include all modifiers of recombination, both beneficial and deleterious. Computer simulations show that beneficial mutations that also increase recombination have an increased chance of fixation. Similarly, deleterious mutations that also decrease recombination have an increased chance of fixation. The results suggest that a simple model of recombination modifiers, including both neutral and pleiotropic modifiers, is a necessary explanation for the evolutionary origin of recombination.

Gene flow and natural selection in the origin of Drosophila pseudoobscura and close relatives

The divergence of Drosophila pseudoobscura and close relatives D. persimilis and D. pseudoobscura bogotana has been studied using comparative DNA sequence data from multiple nuclear loci. New data from the Hsp82 and Adh regions, in conjunction with existing data from Adh and the Period locus, are examined in the light of various models of speciation. The principal finding is that the three loci present very different histories, with Adh indicating large amounts of recent gene flow among the taxa, while little or no gene flow is apparent in the data from the other loci. The data were compared with predictions from several isolation models of divergence. These models include no gene flow, and they were found to be incompatible with the data. Instead the DNA data, taken together with other evidence, seem consistent with divergence models in which natural selection acts against gene flow at some loci more than at others. This family of models includes some sympatric and parapatric speciation models, as well as models of secondary contact and subsequent reinforcement of sexual isolation.

A coalescent estimator of the population recombination rate

Population genetic models often use a population recombination parameter 4Nc, where N is the effective population size and c is the recombination rate per generation. In many ways 4Nc is comparable to 4Nu, the population mutation rate. Both combine genome level and population level processes, and together they describe the rate of production of genetic variation in a population. However, 4Nc is more difficult to estimate. For a population sample of DNA sequences, historical recombination can only be detected if polymorphisms exist, and even then most recombination events are not detectable. This paper describes an estimator of 4Nc, hereafter designated gamma (gamma), that was developed using a coalescent model for a sample of four DNA sequences with recombination. The reliability of gamma was assessed using multiple coalescent simulations. In general gamma has low to moderate bias, and the reliability of gamma is comparable, though less, than that for a widely used estimator of 4Nu. If there exists an independent estimate of the recombination rate (per generation, per base pair), gamma can be used to estimate the effective population size or the neutral mutation rate.

Estimating ancestral population parameters

The expected numbers of different categories of polymorphic sites are derived for two related models of population history the isolation model, in which an ancestral population splits into two descendents, and the size-change model, in which a single population undergoes an instantaneous change in size. For the isolation model, the observed numbers of shared, fixed, and exclusive polymorphic sites are used to estimate the relative sizes of the three populations, ancestral plus two descendent, as well as the time of the split. For the size change model, the numbers of sites segregating at particular frequencies in the sample are used to estimate the relative sizes of the ancestral and descendent populations plus the time the change took place. Parameters are estimated by choosing values that most closely equate expectations with observations. Computer simulations show that current and historical population parameters can be estimated accurately. The methods are applied to DNA data from two species of Drosophila and to some human mitochondrial DNA sequences.

Mitochondrial and nuclear genes present conflicting portraits of human origins

Human mitochondrial DNA (mtDNA) sequences reveal an abundance of polymorphic sites in which the frequencies of the segregating bases are very different. A typical polymorphism involves one base at low frequency and the other base at high frequency. In contrast, nuclear gene data sets tend to show an excess of polymorphisms in which both segregating bases are at intermediate frequencies. A new statistical test of this difference finds significant differences between mtDNA and nuclear gene data sets reported in the literature. However, differences in the polymorphism patterns could be caused by different sample origins for the different data sets. To examine the mtDNA-nuclear difference more closely, DNA sequences were generated from a portion of the X-linked pyruvate dehydrogenase E1 alpha subunit (PDHA1) locus and from a portion of mitochondrial control region I (CRI) from each of eight individuals, four from sub-Saharan Africa. The two genes revealed a significant difference in the site frequency distribution of polymorphic sites. PDHA1 revealed an excess of intermediate-frequency polymorphisms, while CRI showed an excess of sites with the low-high frequency pattern. The discrepancy suggests that mitochondrial variation has been shaped by natural selection, and may not be ideal for some questions on human origins.

The speciation history of Drosophila pseudoobscura and close relatives: inferences from DNA sequence variation at the period locus

Thirty-five period locus sequences from Drosophila pseudoobscura and its siblings species, D. p. bogotana, D. persimilis, and D. miranda, were studied. A large amount of variation was found within D. pseudoobscura and D. persimilis, consistent with histories of large effective population sizes. D. p. bogotana, however, has a severe reduction in diversity. Combined analysis of per with two other loci, in both D. p. bogotana and D. pseudoobscura, strongly suggest this reduction is due to recent directional selection at or near per within D. p. bogotana. Since D. p. bogotana is highly variable and shares variation with D. pseudoobscura at other loci, the low level of variation at per within D. p. bogotana can not be explained by a small effective population size or by speciation via founder event. Both D. pseudoobscura and D. persimilis have considerable intraspecific gene flow. A large portion of one D. persimilis sequence appears to have arisen via introgression from D. pseudoobscura. The time of this event appears to be well after the initial separation of these two species. The estimated times since speciation are one mya for D. pseudoobscura and D. persimilis and 2 mya since the formation of D. miranda.

DNA sequence variation at the period locus reveals the history of species and speciation events in the Drosophila virilis group

The virilis phylad of the Drosophila virilis group consists of five closely related taxa: D. virilis, D. lummei, D. novamexicana, D. americana americana and D. americana texana. DNA sequences from a 2.1-kb pair portion of the period locus were generated in four to eight individuals from each of the five taxa. We found evidence of recombination and high levels of variation within species. We found no evidence of recent natural selection. Surprisingly there was no evidence of divergence between D. a. americana and D. a. texana, and they collectively appear to have had a large historical effective population size. The ranges of these two taxa overlap in a large hybrid zone that has been delineated in the eastern U.S. on the basis of the geographic pattern of a chromosomal fusion. Also surprisingly, D. novamexicana appears to consist of two distinct groups each with low population size and no gene flow between them.

The effects of mutation and natural selection on codon bias in the genes of Drosophila

Codon bias varies widely among the loci of Drosophila melanogaster, and some of this diversity has been explained by variation in the strength of natural selection. A study of correlations between intron and coding region base composition shows that variation in mutation pattern also contributes to codon bias variation. This finding is corroborated by an analysis of variance (ANOVA), which shows a tendency for introns from the same gene to be similar in base composition. The strength of base composition correlations between introns and codon third positions is greater for genes with low codon bias than for genes with high codon bias. This pattern can be explained by an overwhelming effect of natural selection, relative to mutation, in highly biased loci. In particular, this correlation is absent when examining fourfold degenerate sites of highly biased genes. In general, it appears that selection acts more strongly in choosing among fourfold degenerate codons than among twofold degenerate codons. Although the results indicate regional variation in mutational bias, no evidence is found for large scale regions of compositional homogeneity.

Reduced natural selection associated with low recombination in Drosophila melanogaster

Synonymous codons are not used equally in many organisms, and the extent of codon bias varies among loci. Earlier studies have suggested that more highly expressed loci in Drosophila melanogaster are more biased, consistent with findings from several prokaryotes and unicellular eukaryotes that codon bias is partly due to natural selection for translational efficiency. We link this model of varying selection intensity to the population-genetics prediction that the effectiveness of natural selection is decreased under reduced recombination. In analyses of 385 D. melanogaster loci, we find that codon bias is reduced in regions of low recombination (i.e., near centromeres and telomeres and on the fourth chromosome). The effect does not appear to be a linear function of recombination rate; rather, it seems limited to regions with the very lowest levels of recombination. The large majority of the genome apparently experiences recombination at a sufficiently high rate for effective natural selection against suboptimal codons. These findings support models of the Hill-Robertson effect and genetic hitchhiking and are largely consistent with multiple reports of low levels of DNA sequence variation in regions of low recombination.

Population genetics and phylogenetics of DNA sequence variation at multiple loci within the Drosophila melanogaster species complex

Two regions of the genome, a 1-kbp portion of the zeste locus and a 1.1-kbp portion of the yolk protein 2 locus, were sequenced in six individuals from each of four species: Drosophila melanogaster, D. simulans, D. mauritiana, and D. sechellia. The species and strains were the same as those of a previous study of a 1.9-kbp region of the period locus. No evidence was found for recent balancing or directional selection or for the accumulation of selected differences between species. Yolk protein 2 has a high level of amino acid replacement variation and a low level of synonymous variation, while zeste has the opposite pattern. This contrast is consistent with information on gene function and patterns of codon bias. Polymorphism levels are consistent with a ranking of effective population sizes, from low to high, in the following order: D. sechellia, D. melanogaster, D.mauritiana, and D. simulans. The apparent species relationships are very similar to those suggested by the period locus study. In particular, D. simulans appears to be a large population that is still segregating variation that arose before the separation of D. mauritiana and D. sechellia. It is estimated that the separation of ancestral D. melanogaster from the other species occurred 2.5-3.4 Mya. The separations of D. sechellia and D. mauritiana from ancestral D. simulans appear to have occurred 0.58-0.86 Mya, with D. mauritiana having diverged from ancestral D. simulans 0.1 Myr more recently than D. sechellia.

DNA sequence variation at the period locus within and among species of the Drosophila melanogaster complex

A 1.9-kilobase region of the period locus was sequenced in six individuals of Drosophila melanogaster and from six individuals of each of three sibling species: Drosophila simulans, Drosophila sechellia and Drosophila mauritiana. Extensive genealogical analysis of 174 polymorphic sites reveals a complex history. It appears that D. simulans, as a large population still segregating very old lineages, gave rise to the island species D. mauritiana and D. sechellia. Rather than considering these speciation events as having produced "sister" taxa, it seems more appropriate to consider D. simulans a parent species to D. sechellia and D. mauritiana. The order, in time, of these two phylogenetic events remains unclear. D. mauritiana supports a large number of polymorphisms, many of which are shared with D. simulans, and so appears to have begun and persisted as a large population. In contrast, D. sechellia has very little variation and seems to have experienced a severe population bottleneck. Alternatively, the low variation in D. sechellia could be due to recent directional selection and genetic hitchhiking at or near the per locus.