Leprosy in wild chimpanzees.

Humans are considered as the main host for Mycobacterium leprae1, the aetiological agent of leprosy, but spillover has occurred to other mammals that are now maintenance hosts, such as nine-banded armadillos and red squirrels2,3. Although naturally acquired leprosy has also been described in captive nonhuman primates4-7, the exact origins of infection remain unclear. Here we describe leprosy-like lesions in two wild populations of western chimpanzees (Pan troglodytes verus) in Cantanhez National Park, Guinea-Bissau and Taï National Park, Côte d'Ivoire, West Africa. Longitudinal monitoring of both populations revealed the progression of disease symptoms compatible with advanced leprosy. Screening of faecal and necropsy samples confirmed the presence of M. leprae as the causative agent at each site and phylogenomic comparisons with other strains from humans and other animals show that the chimpanzee strains belong to different and rare genotypes (4N/O and 2F). These findings suggest that M. leprae may be circulating in more wild animals than suspected, either as a result of exposure to humans or other unknown environmental sources.

PipeIT2: A tumor-only somatic variant calling workflow for molecular diagnostic Ion Torrent sequencing data.

Precision oncology relies on the accurate identification of somatic mutations in cancer patients. While the sequencing of the tumoral tissue is frequently part of routine clinical care, the healthy counterparts are rarely sequenced. We previously published PipeIT, a somatic variant calling workflow specific for Ion Torrent sequencing data enclosed in a Singularity container. PipeIT combines user-friendly execution, reproducibility and reliable mutation identification, but relies on matched germline sequencing data to exclude germline variants. Expanding on the original PipeIT, here we describe PipeIT2 to address the clinical need to define somatic mutations in the absence of germline control. We show that PipeIT2 achieves a > 95% recall for variants with variant allele fraction >10%, reliably detects driver and actionable mutations and filters out most of the germline mutations and sequencing artifacts. With its performance, reproducibility, and ease of execution, PipeIT2 is a valuable addition to molecular diagnostics laboratories.

A new paradigm for leprosy diagnosis based on host gene expression.

Transcriptional profiling is a powerful tool to investigate and detect human diseases. In this study, we used bulk RNA-sequencing (RNA-Seq) to compare the transcriptomes in skin lesions of leprosy patients or controls affected by other dermal conditions such as granuloma annulare, a confounder for paucibacillary leprosy. We identified five genes capable of accurately distinguishing multibacillary and paucibacillary leprosy from other skin conditions. Indoleamine 2,3-dioxygenase 1 (IDO1) expression alone was highly discriminatory, followed by TLR10, BLK, CD38, and SLAMF7, whereas the HS3ST2 and CD40LG mRNA separated multi- and paucibacillary leprosy. Finally, from the main differentially expressed genes (DEG) and enriched pathways, we conclude that paucibacillary disease is characterized by epithelioid transformation and granuloma formation, with an exacerbated cellular immune response, while multibacillary leprosy features epithelial-mesenchymal transition with phagocytic and lipid biogenesis patterns in the skin. These findings will help catalyze the development of better diagnostic tools and potential host-based therapeutic interventions. Finally, our data may help elucidate host-pathogen interplay driving disease clinical manifestations.

New 2-Ethylthio-4-methylaminoquinazoline derivatives inhibiting two subunits of cytochrome bc1 in Mycobacterium tuberculosis.

The emergence of multi-drug (MDR-TB) and extensively-drug resistant tuberculosis (XDR-TB) is a major threat to the global management of tuberculosis (TB) worldwide. New chemical entities are of need to treat drug-resistant TB. In this study, the mode of action of new, potent quinazoline derivatives was investigated against Mycobacterium tuberculosis (M. tb). Four derivatives 11626141, 11626142, 11626252 and 11726148 showed good activity (MIC ranging from 0.02-0.09 µg/mL) and low toxicity (TD50 ≥ 5µg/mL) in vitro against M. tb strain H37Rv and HepG2 cells, respectively. 11626252 was the most selective compound from this series. Quinazoline derivatives were found to target cytochrome bc1 by whole-genome sequencing of mutants selected with 11626142. Two resistant mutants harboured the transversion T943G (Trp312Gly) and the transition G523A (Gly175Ser) in the cytochrome bc1 complex cytochrome b subunit (QcrB). Interestingly, a third mutant QuinR-M1 contained a mutation in the Rieske iron-sulphur protein (QcrA) leading to resistance to quinazoline and other QcrB inhibitors, the first report of cross-resistance involving QcrA. Modelling of both QcrA and QcrB revealed that all three resistance mutations are located in the stigmatellin pocket, as previously observed for other QcrB inhibitors such as Q203, AX-35, and lansoprazole sulfide (LPZs). Further analysis of the mode of action in vitro revealed that 11626252 exposure leads to ATP depletion, a decrease in the oxygen consumption rate and also overexpression of the cytochrome bd oxidase in M. tb. Our findings suggest that quinazoline-derived compounds are a new and attractive chemical entity for M. tb drug development targeting two separate subunits of the cytochrome bc1 complex.

Emergence of Mycobacterium leprae Rifampin Resistance Evaluated by Whole-Genome Sequencing after 48 Years of Irregular Treatment.

A case of Mycobacterium leprae rifampin resistance after irregular antileprosy treatments since 1971 is reported. Whole-genome sequencing from four longitudinal samples indicated relapse due to acquired rifampin resistance and not to reinfection with another strain. A putative compensatory mutation in rpoC was also detected. Clinical improvement was achieved using an alternative therapy.

EspL is essential for virulence and stabilizes EspE, EspF and EspH levels in Mycobacterium tuberculosis.

The ESX-1, type VII, secretion system represents the major virulence determinant of Mycobacterium tuberculosis, one of the most successful intracellular pathogens. Here, by combining genetic and high-throughput approaches, we show that EspL, a protein of 115 amino acids, is essential for mediating ESX-1-dependent virulence and for stabilization of EspE, EspF and EspH protein levels. Indeed, an espL knock-out mutant was unable to replicate intracellularly, secrete ESX-1 substrates or stimulate innate cytokine production. Moreover, proteomic studies detected greatly reduced amounts of EspE, EspF and EspH in the espL mutant as compared to the wild type strain, suggesting a role for EspL as a chaperone. The latter conclusion was further supported by discovering that EspL interacts with EspD, which was previously demonstrated to stabilize the ESX-1 substrates and effector proteins, EspA and EspC. Loss of EspL also leads to downregulation in M. tuberculosis of WhiB6, a redox-sensitive transcriptional activator of ESX-1 genes. Overall, our data highlight the importance of a so-far overlooked, though conserved, component of the ESX-1 secretion system and begin to delineate the role played by EspE, EspF and EspH in virulence and host-pathogen interaction.

Ancient genomes reveal a high diversity of Mycobacterium leprae in medieval Europe.

Studying ancient DNA allows us to retrace the evolutionary history of human pathogens, such as Mycobacterium leprae, the main causative agent of leprosy. Leprosy is one of the oldest recorded and most stigmatizing diseases in human history. The disease was prevalent in Europe until the 16th century and is still endemic in many countries with over 200,000 new cases reported annually. Previous worldwide studies on modern and European medieval M. leprae genomes revealed that they cluster into several distinct branches of which two were present in medieval Northwestern Europe. In this study, we analyzed 10 new medieval M. leprae genomes including the so far oldest M. leprae genome from one of the earliest known cases of leprosy in the United Kingdom-a skeleton from the Great Chesterford cemetery with a calibrated age of 415-545 C.E. This dataset provides a genetic time transect of M. leprae diversity in Europe over the past 1500 years. We find M. leprae strains from four distinct branches to be present in the Early Medieval Period, and strains from three different branches were detected within a single cemetery from the High Medieval Period. Altogether these findings suggest a higher genetic diversity of M. leprae strains in medieval Europe at various time points than previously assumed. The resulting more complex picture of the past phylogeography of leprosy in Europe impacts current phylogeographical models of M. leprae dissemination. It suggests alternative models for the past spread of leprosy such as a wide spread prevalence of strains from different branches in Eurasia already in Antiquity or maybe even an origin in Western Eurasia. Furthermore, these results highlight how studying ancient M. leprae strains improves understanding the history of leprosy worldwide.

Insight into the evolution and origin of leprosy bacilli from the genome sequence of Mycobacterium lepromatosis.

Mycobacterium lepromatosis is an uncultured human pathogen associated with diffuse lepromatous leprosy and a reactional state known as Lucio's phenomenon. By using deep sequencing with and without DNA enrichment, we obtained the near-complete genome sequence of M. lepromatosis present in a skin biopsy from a Mexican patient, and compared it with that of Mycobacterium leprae, which has undergone extensive reductive evolution. The genomes display extensive synteny and are similar in size (∼3.27 Mb). Protein-coding genes share 93% nucleotide sequence identity, whereas pseudogenes are only 82% identical. The events that led to pseudogenization of 50% of the genome likely occurred before divergence from their most recent common ancestor (MRCA), and both M. lepromatosis and M. leprae have since accumulated new pseudogenes or acquired specific deletions. Functional comparisons suggest that M. lepromatosis has lost several enzymes required for amino acid synthesis whereas M. leprae has a defective heme pathway. M. lepromatosis has retained all functions required to infect the Schwann cells of the peripheral nervous system and therefore may also be neuropathogenic. A phylogeographic survey of 227 leprosy biopsies by differential PCR revealed that 221 contained M. leprae whereas only six, all from Mexico, harbored M. lepromatosis. Phylogenetic comparisons indicate that M. lepromatosis is closer than M. leprae to the MRCA, and a Bayesian dating analysis suggests that they diverged from their MRCA approximately 13.9 Mya. Thus, despite their ancient separation, the two leprosy bacilli are remarkably conserved and still cause similar pathologic conditions.

Rv3852 (H-NS) of Mycobacterium tuberculosis Is Not Involved in Nucleoid Compaction and Virulence Regulation.

A handful of nucleoid-associated proteins (NAPs) regulate the vast majority of genes in a bacterial cell. H-NS, the histone-like nucleoid-structuring protein, is one of these NAPs and protects Escherichia coli from foreign gene expression. Though lacking any sequence similarity with E. coli H-NS, Rv3852 was annotated as the H-NS ortholog in Mycobacterium tuberculosis, as it resembles human histone H1. The role of Rv3852 was thoroughly investigated by immunoblotting, subcellular localization, construction of an unmarked rv3852 deletion in the M. tuberculosis genome, and subsequent analysis of the resulting Δrv3852 strain. We found that Rv3852 was predominantly present in the logarithmic growth phase with a decrease in protein abundance in stationary phase. Furthermore, it was strongly associated with the cell membrane and not detected in the cytosolic fraction, nor was it secreted. The Δrv3852 strain displayed no growth defect or morphological abnormalities. Quantitative measurement of nucleoid localization in the Δrv3852 mutant strain compared to that in the parental H37Rv strain showed no difference in nucleoid position or spread. Infection of macrophages as well as severe combined immunodeficient (SCID) mice demonstrated that loss of Rv3852 had no detected influence on the virulence of M. tuberculosis We thus conclude that M. tuberculosis Rv3852 is not involved in pathogenesis and is not a typical NAP. The existence of an as yet undiscovered Rv3852 ortholog cannot be excluded, although this role is likely played by the well-characterized Lsr2 protein.IMPORTANCEMycobacterium tuberculosis is the causative agent of the lung infection tuberculosis, claiming more than 1.5 million lives each year. To understand the mechanisms of latent infection, where M. tuberculosis can stay dormant inside the human host, we require deeper knowledge of the basic biology and of the regulatory networks. In our work, we show that Rv3852, previously annotated as H-NS, is not a typical nucleoid-associated protein (NAP) as expected from its initial annotation. Rv3852 from M. tuberculosis has neither influence on nucleoid shape or compaction nor a role in virulence. Our findings reduce the repertoire of identified nucleoid-associated proteins in M. tuberculosis to four transcription regulators and underline the importance of genetic studies to assign a function to bacterial genes.

Transmission of Drug-Resistant Leprosy in Guinea-Conakry Detected Using Molecular Epidemiological Approaches.

Molecular drug susceptibility testing was performed on skin biopsies from 24 leprosy patients from Guinea-Conakry for the first time. We identified primary drug resistance in 4 cases and a dapsone-resistant cluster caused by the same strain. Primary transmission of drug-resistant Mycobacterium leprae, including a rifampicin-resistant strain, is reported.

Genomic and transcriptomic analysis of the streptomycin-dependent Mycobacterium tuberculosis strain 18b.

BACKGROUND: The ability of Mycobacterium tuberculosis to establish a latent infection (LTBI) in humans confounds the treatment of tuberculosis. Consequently, there is a need to discover new therapeutic agents that can kill M. tuberculosis both during active disease and LTBI. The streptomycin-dependent strain of M. tuberculosis, 18b, provides a useful tool for this purpose since upon removal of streptomycin (STR) it enters a non-replicating state that mimics latency both in vitro and in animal models. RESULTS: The 4.41 Mb genome sequence of M. tuberculosis 18b was determined and this revealed the strain to belong to clade 3 of the ancient ancestral lineage of the Beijing family. STR-dependence was attributable to insertion of a single cytosine in the 530 loop of the 16S rRNA and to a single amino acid insertion in the N-terminal domain of initiation factor 3. RNA-seq was used to understand the genetic programme activated upon STR-withdrawal and hence to gain insight into LTBI. This revealed reconfiguration of gene expression and metabolic pathways showing strong similarities between non-replicating 18b and M. tuberculosis residing within macrophages, and with the core stationary phase and microaerophilic responses. CONCLUSION: The findings of this investigation confirm the validity of 18b as a model for LTBI, and provide insight into both the evolution of tubercle bacilli and the functioning of the ribosome.

The PhoP-dependent ncRNA Mcr7 modulates the TAT secretion system in Mycobacterium tuberculosis.

The PhoPR two-component system is essential for virulence in Mycobacterium tuberculosis where it controls expression of approximately 2% of the genes, including those for the ESX-1 secretion apparatus, a major virulence determinant. Mutations in phoP lead to compromised production of pathogen-specific cell wall components and attenuation both ex vivo and in vivo. Using antibodies against the native protein in ChIP-seq experiments (chromatin immunoprecipitation followed by high-throughput sequencing) we demonstrated that PhoP binds to at least 35 loci on the M. tuberculosis genome. The PhoP regulon comprises several transcriptional regulators as well as genes for polyketide synthases and PE/PPE proteins. Integration of ChIP-seq results with high-resolution transcriptomic analysis (RNA-seq) revealed that PhoP controls 30 genes directly, whilst regulatory cascades are responsible for signal amplification and downstream effects through proteins like EspR, which controls Esx1 function, via regulation of the espACD operon. The most prominent site of PhoP regulation was located in the intergenic region between rv2395 and PE_PGRS41, where the mcr7 gene codes for a small non-coding RNA (ncRNA). Northern blot experiments confirmed the absence of Mcr7 in an M. tuberculosis phoP mutant as well as low-level expression of the ncRNA in M. tuberculosis complex members other than M. tuberculosis. By means of genetic and proteomic analyses we demonstrated that Mcr7 modulates translation of the tatC mRNA thereby impacting the activity of the Twin Arginine Translocation (Tat) protein secretion apparatus. As a result, secretion of the immunodominant Ag85 complex and the beta-lactamase BlaC is affected, among others. Mcr7, the first ncRNA of M. tuberculosis whose function has been established, therefore represents a missing link between the PhoPR two-component system and the downstream functions necessary for successful infection of the host.

Defining the nitrogen regulated transcriptome of Mycobacterium smegmatis using continuous culture.

BACKGROUND: Nitrogen is essential for microbial growth and its importance is demonstrated by the complex regulatory systems used to control the transport, assimilation and metabolism of nitrogen. Recent studies are beginning to shed light on how mycobacteria respond to nitrogen limitation and several regulators (e.g., GlnR, PII) have been characterized at a molecular level. However, despite this progress, our knowledge of the transcriptional response of mycobacteria to nitrogen limitation and its regulation is confined to batch culture. METHODS: To gain further insight into the response of mycobacteria to nitrogen limitation, we developed a nitrogen-limited chemostat. We compared the transcriptional response of nitrogen-limited cells to carbon-limited cells using RNA-seq analysis in a continuous culture model at a constant growth rate. CONCLUSIONS: Our findings revealed significant changes in the expression of 357 genes (208 upregulated, 149 downregulated; >2-fold change, false discovery rate <5 %) in response to nitrogen limitation in continuous culture. The vast majority of the GlnR regulon (68 %) was differentially expressed under nitrogen limitation in continuous culture and approximately 52 % of the 357 genes overlapped with a previously published study investigating the response of M. smegmatis to nitrogen limitation in batch culture, while expression of only 17 % of the genes identified in batch culture were affected in our chemostat model. Moreover, we identified a unique set of 45 genes involved in the uptake and metabolism of nitrogen that were exclusive to our chemostat model. We observed strong downregulation of pathways for amino acid catabolism (i.e., alanine, aspartate, valine, proline and lysine), suggesting preservation of these amino acids for critical cellular function. We found 16 novel transcriptional regulators that were directly or indirectly involved in the global transcriptomic response of M. smegmatis to nitrogen limitation and identified several non-coding RNAs that might be involved in the transcriptional or post-transcriptional regulation of nitrogen-regulated gene expression. RESULTS: Using nitrogen-limited continuous culture we identified the nitrogen-responsive transcriptome of M. smegmatis, including a number of small non-coding RNAs implicated in controlling nitrogen-regulated gene expression.

Mycobacterium lepromatosis Infections in Nuevo León, Mexico.

The frequency of infection caused by the recently described pathogen Mycobacterium lepromatosis is unknown. Here, we describe the demographics, clinical characteristics, and therapeutic outcomes of five lepromatous leprosy patients suffering from M. lepromatosis infection in Nuevo Léon, Mexico. Diagnosis was facilitated by a new highly specific PCR procedure.

The genome of melon (Cucumis melo L.).

We report the genome sequence of melon, an important horticultural crop worldwide. We assembled 375 Mb of the double-haploid line DHL92, representing 83.3% of the estimated melon genome. We predicted 27,427 protein-coding genes, which we analyzed by reconstructing 22,218 phylogenetic trees, allowing mapping of the orthology and paralogy relationships of sequenced plant genomes. We observed the absence of recent whole-genome duplications in the melon lineage since the ancient eudicot triplication, and our data suggest that transposon amplification may in part explain the increased size of the melon genome compared with the close relative cucumber. A low number of nucleotide-binding site-leucine-rich repeat disease resistance genes were annotated, suggesting the existence of specific defense mechanisms in this species. The DHL92 genome was compared with that of its parental lines allowing the quantification of sequence variability in the species. The use of the genome sequence in future investigations will facilitate the understanding of evolution of cucurbits and the improvement of breeding strategies.

DNA Methylation Landscapes of Prostate Cancer Brain Metastasis Are Shaped by Early Driver Genetic Alterations.

Metastases from primary prostate cancers to rare locations, such as the brain, are becoming more common due to longer life expectancy resulting from improved treatments. Epigenetic dysregulation is a feature of primary prostate cancer, and distinct DNA methylation profiles have been shown to be associated with the mutually exclusive SPOP-mutant or TMPRSS2-ERG fusion genetic backgrounds. Using a cohort of prostate cancer brain metastases (PCBM) from 42 patients, with matched primary tumors for 17 patients, we carried out a DNA methylation analysis to examine the epigenetic distinction between primary prostate cancer and PCBM, the association between epigenetic alterations and mutational background, and particular epigenetic alterations that may be associated with PCBM. Multiregion sampling of PCBM revealed epigenetic stability within metastases. Aberrant methylation in PCBM was associated with mutational background and PRC2 complex activity, an effect that is particularly pronounced in SPOP-mutant PCBM. While PCBM displayed a CpG island hypermethylator phenotype, hypomethylation at the promoters of genes involved in neuroactive ligand-receptor interaction and cell adhesion molecules such as GABRB3, CLDN8, and CLDN4 was also observed, suggesting that cells from primary tumors may require specific reprogramming to form brain metastasis. This study revealed the DNA methylation landscapes of PCBM and the potential mechanisms and effects of PCBM-associated aberrant DNA methylation. SIGNIFICANCE: DNA methylation analysis reveals the molecular characteristics of PCBM and may serve as a starting point for efforts to identify and target susceptibilities of these rare metastases.

Circulating Cell-Free DNA Captures the Intratumor Heterogeneity in Multinodular Hepatocellular Carcinoma.

PURPOSE: Hepatocellular carcinoma (HCC) is a highly heterogeneous disease, with more than 40% of patients initially diagnosed with multinodular HCCs. Although circulating cell-free DNA (cfDNA) has been shown to effectively detect somatic mutations, little is known about its utility to capture intratumor heterogeneity in patients with multinodular HCC undergoing systemic treatment. MATERIALS AND METHODS: Tumor biopsies and plasma were synchronously collected from seven prospectively recruited patients with HCC before and during systemic therapy. Plasma-derived cfDNA and matched germline were subjected to high-depth targeted sequencing with molecular barcoding. The mutational profile of the cfDNA was compared with whole-exome sequencing from matched tumor biopsies. RESULTS: Genomic data revealed that out of the seven patients, five were considered intrahepatic metastasis and two multicentric HCCs. cfDNA captured the majority of mutations in the tumors and detected significantly more mutations than tumor biopsies. Driver mutations such as CTNNB1 S33C, NRAS Q61R, ARID1A R727fs, and NF1 E2368fs as well as standard-of-care biomarkers of response to targeted therapy were detected only in cfDNA. In the two patients with multicentric HCC, cfDNA detected mutations derived from the genetically independent and spatially distinct nodules. Moreover, cfDNA was not only able to capture clonal mutations but also the subclonal mutations detected in only one of the multiple biopsied nodules. Furthermore, serial cfDNA detected variants of tumor origin emerging during treatment. CONCLUSION: This study revealed that the genetic analysis of cfDNA captures the intratumor heterogeneity in multinodular HCC highlighting the potential for cfDNA as a sensitive and noninvasive tool for precision medicine.

Siglec-7 represents a glyco-immune checkpoint for non-exhausted effector memory CD8+ T cells with high functional and metabolic capacities.

While inhibitory Siglec receptors are known to regulate myeloid cells, less is known about their expression and function in lymphocytes subsets. Here we identified Siglec-7 as a glyco-immune checkpoint expressed on non-exhausted effector memory CD8+ T cells that exhibit high functional and metabolic capacities. Seahorse analysis revealed higher basal respiration and glycolysis levels of Siglec-7+ CD8+ T cells in steady state, and particularly upon activation. Siglec-7 polarization into the T cell immune synapse was dependent on sialoglycan interactions in trans and prevented actin polarization and effective T cell responses. Siglec-7 ligands were found to be expressed on both leukemic stem cells and acute myeloid leukemia (AML) cells suggesting the occurrence of glyco-immune checkpoints for Siglec-7+ CD8+ T cells, which were found in patients' peripheral blood and bone marrow. Our findings project Siglec-7 as a glyco-immune checkpoint and therapeutic target for T cell-driven disorders and cancer.

Gene expression patterns associated with multidrug therapy in multibacillary leprosy.

Multidrug therapy (MDT) has been successfully used in the treatment of leprosy. However, although patients are cured after the completion of MDT, leprosy reactions, permanent disability, and occasional relapse/reinfection are frequently observed in patients. The immune system of multibacillary patients (MB) is not able to mount an effective cellular immune response against M. leprae. Consequently, clearance of bacilli from the body is a slow process and after 12 doses of MDT not all MB patients reduce bacillary index (BI). In this context, we recruited MB patients at the uptake and after 12-month of MDT. Patients were stratified according to the level of reduction of the BI after 12 doses MDT. A reduction of at least one log in BI was necessary to be considered a responder patient. We evaluated the pattern of host gene expression in skin samples with RNA sequencing before and after MDT and between samples from patients with or without one log reduction in BI. Our results demonstrated that after 12 doses of MDT there was a reduction in genes associated with lipid metabolism, inflammatory response, and cellular immune response among responders (APOBEC3A, LGALS17A, CXCL13, CXCL9, CALHM6, and IFNG). Also, by comparing MB patients with lower BI reduction versus responder patients, we identified high expression of CDH19, TMPRSS4, PAX3, FA2H, HLA-V, FABP7, and SERPINA11 before MDT. From the most differentially expressed genes, we observed that MDT modulates pathways related to immune response and lipid metabolism in skin cells from MB patients after MDT, with higher expression of genes like CYP11A1, that are associated with cholesterol metabolism in the group with the worst response to treatment. Altogether, the data presented contribute to elucidate gene signatures and identify differentially expressed genes associated with MDT outcomes in MB patients.

Determination of the melon chloroplast and mitochondrial genome sequences reveals that the largest reported mitochondrial genome in plants contains a significant amount of DNA having a nuclear origin.

BACKGROUND: The melon belongs to the Cucurbitaceae family, whose economic importance among vegetable crops is second only to Solanaceae. The melon has a small genome size (454 Mb), which makes it suitable for molecular and genetic studies. Despite similar nuclear and chloroplast genome sizes, cucurbits show great variation when their mitochondrial genomes are compared. The melon possesses the largest plant mitochondrial genome, as much as eight times larger than that of other cucurbits. RESULTS: The nucleotide sequences of the melon chloroplast and mitochondrial genomes were determined. The chloroplast genome (156,017 bp) included 132 genes, with 98 single-copy genes dispersed between the small (SSC) and large (LSC) single-copy regions and 17 duplicated genes in the inverted repeat regions (IRa and IRb). A comparison of the cucumber and melon chloroplast genomes showed differences in only approximately 5% of nucleotides, mainly due to short indels and SNPs. Additionally, 2.74 Mb of mitochondrial sequence, accounting for 95% of the estimated mitochondrial genome size, were assembled into five scaffolds and four additional unscaffolded contigs. An 84% of the mitochondrial genome is contained in a single scaffold. The gene-coding region accounted for 1.7% (45,926 bp) of the total sequence, including 51 protein-coding genes, 4 conserved ORFs, 3 rRNA genes and 24 tRNA genes. Despite the differences observed in the mitochondrial genome sizes of cucurbit species, Citrullus lanatus (379 kb), Cucurbita pepo (983 kb) and Cucumis melo (2,740 kb) share 120 kb of sequence, including the predicted protein-coding regions. Nevertheless, melon contained a high number of repetitive sequences and a high content of DNA of nuclear origin, which represented 42% and 47% of the total sequence, respectively. CONCLUSIONS: Whereas the size and gene organisation of chloroplast genomes are similar among the cucurbit species, mitochondrial genomes show a wide variety of sizes, with a non-conserved structure both in gene number and organisation, as well as in the features of the noncoding DNA. The transfer of nuclear DNA to the melon mitochondrial genome and the high proportion of repetitive DNA appear to explain the size of the largest mitochondrial genome reported so far.