Search

International audience; The proper tissue-specific regulation of gene expression is essential for development and homeostasis in metazoans. However, the illegitimate expression of normally tissue-restricted genes-like testis- or placenta-specific genes-is frequently observed in tumors; this promotes transformation, but also allows immunotherapy. Two important questions are: how is the expression of these genes controlled in healthy cells? And how is this altered in cancer? To address these questions, we used an unbiased approach to test the ability of 350 distinct genetic or epigenetic perturbations to induce the illegitimate expression of over 40 tissue-restricted genes in primary human cells. We find that almost all of these genes are remarkably resistant to reactivation by a single alteration in signaling pathways or chromatin regulation. However, a few genes differ and are more readily activated; one is the placenta-expressed gene ADAM12, which promotes invasion. Using cellular systems, an animal model, and bioinformatics, we find that a non-canonical but druggable TGF-β/KAT2A/TAK1 axis controls ADAM12 induction in normal and cancer cells. More broadly, our data show that illegitimate gene expression in cancer is an heterogeneous phenomenon, with a few genes activatable by simple events, and most genes likely requiring a combination of events to become reactivated.

International audience; To better study and manage chestnut trees and species, we identified nuclear single nucleotide polymorphism (SNP) markers using restriction-associated DNA sequencing. Out of 343 loci tested, 68 SNP markers were selected that withhold stringent quality criteria such as quasi-systematic amplification across species and Mendelian segregation in both purebred and hybrid individuals. They provide sufficient power for species, hybrids and backcross characterization as well as for clonal identification, as shown by a comparison with single sequenced repeat (SSR) loci.

DNA barcoding and metabarcoding is increasingly used to effectively and precisely assess and monitor biodiversity in aquatic ecosystems. As these methods rely on data availability and quality of barcode reference libraries, it is important to develop and follow best practices to ensure optimal quality and traceability of the metadata associated with the reference barcodes used for identification. Sufficient metadata, as well as vouchers, corresponding to each reference barcode must be available to ensure reliable barcode library curation and, thereby, provide trustworthy baselines for downstream molecular species identification. This document (1) specifies the data and metadata required to ensure the relevance, the accessibility and traceability of DNA barcodes and (2) specifies the recommendations for DNA harvesting and for the storage of both voucher specimens/samples and barcode data, (undefined)

Background Multidrug-resistant (MDR) Mycobacterium tuberculosis complex strains not detected by commercial molecular drug susceptibility testing (mDST) assays due to the RpoB I491F resistance mutation are threatening the control of MDR tuberculosis (MDR-TB) in Eswatini. Methods We investigate the evolution and spread of MDR strains in Eswatini with a focus on bedaquiline (BDQ) and clofazimine (CFZ) resistance using whole-genome sequencing in two collections ((1) national drug resistance survey, 2009–2010; (2) MDR strains from the Nhlangano region, 2014–2017). Results MDR strains in collection 1 had a high cluster rate (95%, 117/123 MDR strains) with 55% grouped into the two largest clusters (gCL3, n = 28; gCL10, n = 40). All gCL10 isolates, which likely emerged around 1993 (95% highest posterior density 1987–1998), carried the mutation RpoB I491F that is missed by commercial mDST assays. In addition, 21 (53%) gCL10 isolates shared a Rv0678 M146T mutation that correlated with elevated minimum inhibitory concentrations (MICs) to BDQ and CFZ compared to wild type isolates. gCL10 isolates with the Rv0678 M146T mutation were also detected in collection 2. Conclusion The high clustering rate suggests that transmission has been driving the MDR-TB epidemic in Eswatini for three decades. The presence of MDR strains in Eswatini that are not detected by commercial mDST assays and have elevated MICs to BDQ and CFZ potentially jeopardizes the successful implementation of new MDR-TB treatment guidelines. Measures to limit the spread of these outbreak isolates need to be implemented urgently.

Yellow fever is a vector-borne disease endemic in tropical regions of Africa, where 90% of the global burden occurs, and Latin America. It is notoriously under-reported with uncertainty arising from a complex transmission cycle including a sylvatic reservoir and non-specific symptom set. Resulting estimates of burden, particularly in Africa, are highly uncertain. We examine two established models of yellow fever transmission within a Bayesian model averaging framework in order to assess the relative evidence for each model’s assumptions and to highlight possible data gaps. Our models assume contrasting scenarios of the yellow fever transmission cycle in Africa. The first takes the force of infection in each province to be static across the observation period; this is synonymous with a constant infection pressure from the sylvatic reservoir. The second model assumes the majority of transmission results from the urban cycle; in this case, the force of infection is dynamic and defined through a fixed value of R0 in each province. Both models are coupled to a generalised linear model of yellow fever occurrence which uses environmental covariates to allow us to estimate transmission intensity in areas where data is sparse. We compare these contrasting descriptions of transmission through a Bayesian framework and trans-dimensional Markov chain Monte Carlo sampling in order to assess each model’s evidence given the range of uncertainty in parameter values. The resulting estimates allow us to produce Bayesian model averaged predictions of yellow fever burden across the African endemic region. We find strong support for the static force of infection model which suggests a higher proportion of yellow fever transmission occurs as a result of infection from an external source such as the sylvatic reservoir. However, the model comparison highlights key data gaps in serological surveys across the African endemic region. As such, conclusions concerning the most prevalent transmission routes for yellow fever will be limited by the sparsity of data which is particularly evident in the areas with highest predicted transmission intensity. Our model and estimation approach provides a robust framework for model comparison and predicting yellow fever burden in Africa. However, key data gaps increase uncertainty surrounding estimates of model parameters and evidence. As more mathematical models are developed to address new research questions, it is increasingly important to compare them with established modelling approaches to highlight uncertainty in structures and data., Author summary Yellow fever (YF) is notoriously under reported due to non-specific symptom spectrum and the true burden is highly uncertain as a result of a complex transmission cycle. As such, estimates surrounding YF burden are highly uncertain and the mechanisms behind transmission are often unclear. We assess these mechanisms and the resulting uncertainty by estimating two existing models of YF transmission within a product space framework. This allows us to produce updated estimates of transmission intensity and to compare the relative support for each model given the data. We find strong support for a model assuming a static force of infection, approximating the constant infection pressure from the sylvatic reservoir of YF. We also highlight areas where data is sparse, often the same areas estimated to have especially high transmission intensity. This is the first robust multi-model approach to applied YF modelling and provides a framework that could be extended to other disease models.

International audience; Macrophages are key cells of innate immune response and serve as the first line of defense against bacteria. Transcription profiling of bacteria-infected macrophages could provide important insights on the pathogenicity and host defense mechanisms during infection. We have examined transcription profiles of bovine monocyte-derived macrophages (bMDMs) isolated from the blood of 12 animals and infected in vitro with two strains of Streptococcus agalactiae. Illumina sequencing of RNA from 36 bMDMs cultures exposed in vitro to either one of two sequence types of S. agalactiae (ST103 or ST12) for 6 h and unchallenged controls was performed. Analyses of over 1,656 million high-quality paired-end sequence reads revealed 5,936 and 6,443 differentially expressed genes (p < 0.05) in bMDMs infected with ST103 and ST12, respectively, versus unchallenged controls. Moreover, 588 genes differentially expressed between bMDMs infected with ST103 versus ST12 were identified. Ingenuity pathway analysis of the differentially up-regulated genes in the bMDMs infected with ST103 revealed significant enrichment for granulocyte adhesion and diapedesis, while significant enrichment for the phagosome formation pathway was found among down-regulated genes. Moreover, Ingenuity pathway analysis of the differentially up-regulated genes in the bMDMs infected with ST12 showed significant enrichment for type 1/type 2 T helper cell activation, while the complement activation pathway was overrepresented in the down-regulated genes. Our study identified pathogen-induced regulation of key genes and pathways involved in the immune response of macrophages against infection but also likely involved in bacterial evasion of the host immune system. These results may contribute to better understanding of the mechanisms underlying subclinical infection such as bovine streptococcal mastitis.

A thorough understanding of the heritability, genetic correlations and additive and non-additive variance components of tree growth and wood properties is a requisite for effective tree breeding. This knowledge is essential to maximize genetic gain, that is, the amount of increase in trait performance achieved annually through directional selection. Understanding the genetic attributes of traits targeted by breeding is also important to sustain decade-long genetic progress, that is, the progress made by increasing the average genetic value of the offspring as compared to that of the parental generation. In this study, we report quantitative genetic parameters for fifteen growth, wood chemical and physical traits for the world-famous Eucalyptus urograndis hybrid (E. grandis × E. urophylla). These traits directly impact the optimal use of wood for cellulose pulp, paper, and energy production. A population of 1,000 trees sampled in a progeny trial was phenotyped directly or following the development and use of near-infrared spectroscopy calibration models. Trees were genotyped with 33,398 SNPs and 24,001 DArT-seq genome-wide markers and genomic realized relationship matrices (GRM) were used for parameter estimation with an individual-tree additive-dominant mixed model. Wood chemical properties and wood density showed stronger genetic control than growth, cellulose and fiber traits. Additive effects are the main drivers of genetic variation for all traits, but dominance plays an equally or more important role for growth, singularly in this hybrid. GRM´s with >10,000 markers provided stable relationships estimates and more accurate parameters than pedigrees by capturing the full genetic relationships among individuals and disentangling the non-additive from the additive genetic component. Low correlations between growth and wood properties indicate that simultaneous selection for wood traits can be applied with minor effects on genetic gain for growth. Conversely, moderate to strong correlations between wood density and chemical traits exist, likely due to their interdependency on cell wall structure such that responses to selection will be connected for these traits. Our results illustrate the advantage of using genomewide marker data to inform tree breeding in general and have important consequences for operational breeding of eucalypt urograndis hybrids. Fil: de Lima, Bruno Marco. Universidade de Sao Paulo; Brasil. Ministerio da Agricultura Pecuaria e Abastecimento de Brasil. Empresa Brasileira de Pesquisa Agropecuaria; Brasil Fil: Cappa, Eduardo Pablo. Instituto Nacional de Tecnología Agropecuaria. Centro de Investigación de Recursos Naturales. Instituto de Recursos Biológicos; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas; Argentina Fil: Silva Junior, Orzenil B.. Ministerio da Agricultura Pecuaria e Abastecimento de Brasil. Empresa Brasileira de Pesquisa Agropecuaria; Brasil. Universidade Católica de Brasília; Brasil Fil: Garcia, Carla. No especifíca; Fil: Mansfield, Shawn D.. University of British Columbia; Canadá Fil: Grattapaglia, Dario. Ministerio da Agricultura Pecuaria e Abastecimento de Brasil. Empresa Brasileira de Pesquisa Agropecuaria; Brasil. Universidade Católica de Brasília; Brasil

Erratum in : Correction for Adam et al., Evolutionary history of carbon monoxide dehydrogenase/acetyl-CoA synthase, one of the oldest enzymatic complexes. [Proc Natl Acad Sci U S A. 2018]; International audience; Carbon monoxide dehydrogenase/acetyl-CoA synthase (CODH/ ACS) is a five-subunit enzyme complex responsible for the carbonyl branch of the Wood-Ljungdahl (WL) pathway, considered one of the most ancient metabolisms for anaerobic carbon fixation, but its origin and evolutionary history have been unclear. While traditionally associated with methanogens and acetogens, the presence of CODH/ACS homologs has been reported in a large number of un-cultured anaerobic lineages. Here, we have carried out an exhaustive phylogenomic study of CODH/ACS in over 6,400 archaeal and bacterial genomes. The identification of complete and likely functional CODH/ACS complexes in these genomes significantly expands its distribution in microbial lineages. The CODH/ACS complex displays astounding conservation and vertical inheritance over geological times. Rare intradomain and interdomain transfer events might tie into important functional transitions, including the acquisition of CODH/ACS in some archaeal methanogens not known to fix carbon, the tinkering of the complex in a clade of model bacterial acetogens, or emergence of archaeal-bacterial hybrid complexes. Once these transfers were clearly identified, our results allowed us to infer the presence of a CODH/ACS complex with at least four subunits in the last universal common ancestor (LUCA). Different scenarios on the possible role of ancestral CODH/ACS are discussed. Despite common assumptions, all are equally compatible with an autotrophic, mixotrophic, or heterotrophic LUCA. Functional characterization of CODH/ACS from a larger spectrum of bacterial and archaeal lineages and detailed evolutionary analysis of the WL methyl branch will help resolve this issue.

International audience; Many different environmental factors (nutrients, pollutants, chemicals, physical activity, lifestyle, physical and mental stress) can modulate epigenetic markers in the developing and adult organism. Epigenetics, in turn, can cause and is associated with several neurodegenerative and aging-dependent human diseases. Alzheimer's disease certainly represents one of the most relevant neurodegenerative disorders due to its incidence and its huge socio-economic impact. Therefore, it is easy to understand why recent literature focuses on the epigenetic modifications associated with Alzheimer's disease and other neurodegenerative disorders. One of the most intriguing and, at the same time, worrying evidence is that even "mild" environmental factors (such as behavioral or physical stress) as well as the under-threshold exposure to pollutants and chemicals, can be effective. Finally, even mild nutrients disequilibria can result in long-lasting and functional alterations of many epigenetic markers, although they don't have an immediate acute effect. Therefore, we will probably have to re-define the current risk threshold for many factors, molecules and stresses. Among the many different environmental factors affecting the epigenome, nutrition represents one of the most investigated fields; the reasons are probably that each person interacts with nutrients and that, in turn, nutrients can modulate at molecular level the epigenetic biochemical pathways. The role that nutrition can exert in modulating epigenetic modifications in Alzheimer's disease will be discussed with particular emphasis on the role of B vitamins and DNA methylation.

Background Leptospirosis is a worldwide zoonosis that is endemic in tropical areas, such as Reunion Island. The species Leptospira interrogans is the primary agent in human infections, but other pathogenic species, such as L. kirschner and L. borgpetersenii, are also associated with human leptospirosis. Methods and Findings In this study, a melting curve analysis of the products that were amplified with the primer pairs lfb1 F/R and G1/G2 facilitated an accurate species classification of Leptospira reference strains. Next, we combined an unsupervised high resolution melting (HRM) method with a new statistical approach using primers to amplify a two variable-number tandem-repeat (VNTR) for typing at the subspecies level. The HRM analysis, which was performed with ScreenClust Software, enabled the identification of genotypes at the serovar level with high resolution power (Hunter-Gaston index 0.984). This method was also applied to Leptospira DNA from blood samples that were obtained from Reunion Island after 1998. We were able to identify a unique genotype that is identical to that of the L. interrogans serovars Copenhageni and Icterohaemorrhagiae, suggesting that this genotype is the major cause of leptospirosis on Reunion Island. Conclusions Our simple, rapid, and robust genotyping method enables the identification of Leptospira strains at the species and subspecies levels and supports the direct genotyping of Leptospira in biological samples without requiring cultures.

International audience; The prostate-specific antigen (PSA) is the main diagnostic biomarker for prostate cancer in clinical use, but it lacks specificity and sensitivity, particularly in low dosage values 1. 'How to use PSA' remains a current issue, either for diagnosis as a gray zone corresponding to a concentration in serum of 2.5-10 ng/ml which does not allow a clear differentiation to be made between cancer and noncancer 2 or for patient follow-up as analysis of post-operative PSA kinetic parameters can pose considerable challenges for their practical application 3,4. Alternatively, noncoding RNAs (ncRNAs) are emerging as key molecules in human cancer, with the potential to serve as novel markers of disease, e.g. PCA3 in prostate cancer 5,6 and to reveal uncharacterized aspects of tumor biology. Moreover, data from the ENCODE project published in 2012 showed that different RNA types cover about 62% of the genome. It also appears that the amount of transcriptional regulatory motifs is at least 4.5x higher than the one corresponding to protein-coding exons. Thus, long terminal repeats (LTRs) of human endogenous retroviruses (HERVs) constitute a wide range of putative/candidate transcriptional regulatory sequences, as it is their primary function in infectious retroviruses. HERVs, which are spread throughout the human genome, originate from ancestral and independent infections within the germ line, followed by copy-paste propagation processes and leading to multicopy families occupying 8% of the human genome (note that exons span 2% of our genome). Some HERV loci still express proteins that have been associated with several pathologies including cancer 7-10. We have designed a high-density microarray, in Affymetrix format, aiming to optimally characterize individual HERV loci expression, in order to better understand whether they can be active, if they drive ncRNA transcription or modulate coding gene expression. This tool has been applied in the prostate cancer field (Figure 1).

Background Drug susceptible clinical isolates of Candida albicans frequently become highly tolerant to drugs during chemotherapy, with dreadful consequences to patient health. We used RNA sequencing (RNA-seq) to analyze the transcriptomes of a CDR (Candida Drug Resistance) strain and its isogenic drug sensitive counterpart. Results RNA-seq unveiled differential expression of 228 genes including a) genes previously identified as involved in CDR, b) genes not previously associated to the CDR phenotype, and c) novel transcripts whose function as a gene is uncharacterized. In particular, we show for the first time that CDR acquisition is correlated with an overexpression of the transcription factor encoding gene CZF1. CZF1 null mutants were susceptible to many drugs, independently of known multidrug resistance mechanisms. We show that CZF1 acts as a repressor of β-glucan synthesis, thus negatively regulating cell wall integrity. Finally, our RNA-seq data allowed us to identify a new transcribed region, upstream of the TAC1 gene, which encodes the major CDR transcriptional regulator. Conclusion Our results open new perspectives of the role of Czf1 and of our understanding of the transcriptional and post-transcriptional mechanisms that lead to the acquisition of drug resistance in C. albicans, with potential for future improvements of therapeutic strategies.

By transactivating expression of miRNAs that repress expression of the ZEB1 and ZEB2 transcription factors, p53 inhibits the epithelial–mesenchymal transition., p53 suppresses tumor progression and metastasis. Epithelial–mesenchymal transition (EMT) is a key process in tumor progression and metastasis. The transcription factors ZEB1 and ZEB2 promote EMT. Here, we show that p53 suppresses EMT by repressing expression of ZEB1 and ZEB2. By profiling 92 primary hepatocellular carcinomas (HCCs) and 9 HCC cell lines, we found that p53 up-regulates microRNAs (miRNAs), including miR-200 and miR-192 family members. The miR-200 family members transactivated by p53 then repress ZEB1/2 expression. p53-regulated miR-192 family members also repress ZEB2 expression. Inhibition or overexpression of the miRNAs affects p53-regulated EMT by altering ZEB1 and ZEB2 expression. Our findings indicate that p53 can regulate EMT, and that p53-regulated miRNAs are critical mediators of p53-regulated EMT.

International audience; We investigated space-use patterns of the West African savanna buffalo (Syncerus caffer brachyceros), a little-studied subspecies occurring at the northern limit of the African buffalo's geographical range. This buffalo generally ranges in small herds (about 45 individuals) and has a low body mass (approximately 400 kg) relative to the Cape buffalo (S. c. caffer). We monitored the movements of 7 breeding herds in W Regional Park (Burkina Faso, Benin, Niger) using global positioning system collars and activity data loggers. Habitat selection was analyzed at both large (interseasonal) and small (intraseasonal) scales in a context where resources are segregated spatially at some times of year. Both biotic (primary production and vegetation types) and abiotic (timing of rainfall and surface water) covariates, and the extent to which neighboring herds shared space, were considered. In the dry season buffalo herds ranged close (within 5.3 +/- 2.0 km, mean +/- SD) to segments of permanent rivers. At the onset of the monsoon all herds but 1 (which had year-round access to suitable resources) performed a large (35 +/- 10 km) directional movement in response to a large-scale gradient of primary production. Spatiotemporal dynamics of forage and water resources thus jointly stimulated interseasonal directional movements and shaped large (335 +/- 167 km(2)) annual home ranges. Furthermore, the establishment of home ranges in the wet season appears to be conditioned by a threshold (about 10%) in the availability of perennial grasses. Habitat-selection analysis at intraseasonal scale also underlines the key role played by perennial grasses for buffaloes. The spatial arrangements of home ranges of neighboring herds also suggest that interherd behavioral avoidance is a high-level constraint on foraging processes. The ability of the African buffalo to cope with contrasting environmental conditions throughout most sub-Saharan ecosystems highlights the high behavioral plasticity of this species.

The use of morphological characters for the classification of cyanobacteria has often led to ambiguous strain assignment. In the past two decades, the availability of sequences, such as those of the 16S rRNA,nif,cpcandrpoC1genes, and the use of metagenomics, has steadily increased and has made the reconstruction of evolutionary relationships of some cyanobacterial groups possible in addition to improving strain assignment. Conserved indels (insertions/deletions) are present in all cyanobacterial RpoB (βsubunit of RNA polymerase) sequences presently available in public databases. These indels are located in the Rpb2_6 domain of RpoB, which is involved in DNA binding and DNA-directed RNA polymerase activity. They are variable in length (6–44 aa) and sequence, and form part of what appears to be a longer signature sequence (43–81 aa). Indeed, a number of these sequences turn out to be distinctive among several strains of a given genus and even among strains of a given species. These signature sequences can thus be used to identify cyanobacteria at a subgenus level and can be useful molecular markers to establish the taxonomic positions of cyanobacterial isolates in laboratory cultures, and/or to assess cyanobacterial biodiversity in space and time in natural ecosystems.

International audience; Laboratory mice are well known to be highly susceptible to virulent strains of Yersinia pestis in experimental models of bubonic plague. We have found that Mus spretus-derived SEG/Pas (SEG) mice are exceptionally resistant to virulent CO92 and 6/69 wild type strains. Upon subcutaneous injection of 10(2) colony-forming units (CFU), 90% of females and 68% of males survived, compared with only an 8% survival rate for both male and female C57BL/6 mice. Furthermore, half of the SEG mice survived a challenge of up to 10(7) CFU. The time required for mortality was similar between B6 and SEG, suggesting that survival is dependent on early rather than late processes. The analysis of 322 backcross mice identified three significant quantitative trait loci (QTLs) on chromosomes 3, 4 and 6, with dominant SEG protective alleles. Each QTL increased the survival rate by approximately 20%. The three QTLs function additively, thereby accounting for 67% of the difference between the parental phenotypes. Mice heterozygous for the three QTLs were just as resistant as SEG mice to Y. pestis challenge. The SEG strain therefore offers an invaluable opportunity to unravel mechanisms and underlying genetic factors of resistance against Y. pestis infection.

International audience; Chrysomyxa rusts are fungal pathogens widely present in the boreal forest. Taxonomic delimitation and precise species identification are difficult within this genus because several species display similar morphological features. We applied a DNA barcode system based on the ribosomal internal transcribed spacer region (ITS), large subunit (28S) ribosomal RNA gene, mitochondrial cytochrome oxidase 1 (CO1) and mitochondrial NADH dehydrogenase subunit 6 (NAD6) in 86 strains from 16 different Chrysomyxa species, including members of the Chrysomyxa ledi species complex. The nuclear ITS and 28S loci revealed higher resolving power than the mitochondrial genes. Amplification of the full CO1 barcode region failed due to the presence of introns limiting the dataset obtained with this barcode. In most cases the ITS barcodes were in agreement with taxonomic species based on phenotypic characters. Nevertheless we observed genetically distinct (different DNA barcodes) lineages within Chrysomyxa pyrolae and Chrysomyxa rhododendri, providing some evidence for allopatric speciation within these morphologically defined species. This finding, together with the observed pattern of host specificities of the studied rust fungi, suggest that species diversification within the C. ledi species complex might be governed by a set of factors such as specialisation to certain Ericaceae species as telial hosts and to a lesser extent specialization to different spruce species as aecial hosts. Moreover allopatric speciation by geographic disruption of species also seems to take place. When our data were integrated into a broader phylogenetic framework the Chrysomyxa genus unexpectedly was not resolved as a monophyletic group. Indeed the spruce cone rusts C. pyrolae and C. monesis coalesced with the pine needle rusts belonging to the genus Coleosporium, whereas the microcyclic species Chrysomyxa weirii was embedded within a clade comprising the genus Melampsora.

Phlebotomine sandflies are vectors of phleboviruses that cause sandfly fever or meningitis with significant implications for public health. Although several strains of these viruses had been isolated in Iran in the late 1970's, there was no recent data about the present situation at the outset of this study. Entomological investigations performed in 2009 and 2011 in Iran collected 4,770 sandflies from 10 different regions. Based on morphological identification, they were sorted into 315 pools according to species, sex, trapping station and date of capture. A phlebovirus, provisionally named Dashli virus (DASHV), was isolated from one pool of Sergentomyia spp, and subsequently DASHV RNA was detected in a second pool of Phlebotomus papatasi. Genetic and phylogenetic analyses based on complete coding genomic sequences indicated that (i) DASHV is most closely related to the Iranian isolates of Sandfly fever Sicilian virus [SFSV], (ii) there is a common ancestor to DASHV, Sandfly fever Sicilian- (SFS) and SFS-like viruses isolated in Italy, India, Turkey, and Cyprus (lineage I), (iii) DASHV is more distantly related with Corfou and Toros viruses (lineage II) although common ancestry is supported with 100% bootstrap, (iii) lineage I can be subdivided into sublineage Ia including all SFSV, SFCV and SFTV except those isolated in Iran which forms sublineage Ib (DASHV). Accordingly, we suggest to approve Sandfly fever Sicilian virus species consisting of the all aforementioned viruses. Owing that most of these viruses have been identified in human patients with febrile illness, DASHV should be considered as a potential human pathogen in Iran., Author summary Phlebotomine sandflies are vectors of phleboviruses that cause sandfly fever or meningitis with significant implications for public health. Although several strains of these viruses had been isolated in Iran in the late 1970's, there was no recent data about the present situation at the outset of this study. Entomological investigations performed in 2009 and 2011 in Iran collected 4,770 sandflies from 10 different regions. A phlebovirus, provisionally named Dashli virus (DASHV), was isolated / detected in two pools. DASHV strain was isolated in cell culture and complete genome sequence was determined. Sequence analysis indicated that (i) DASHV is most closely related to the Iranian isolates of Sandfly fever Sicilian virus [SFSV], a virus that is known to cause self-resolutive incapacitating febrile illness in humans, (ii) there is a common ancestor to DASHV and all other variants of SFSV isolated in Italy, India, Turkey, and Cyprus (lineage I), (iii) DASHV is more distantly related with Corfou and Toros viruses (lineage II) although common ancestry is supported with 100% bootstrap, (iii) lineage I can be subdivided into sublineage Ia including all SFSV strains, whereas Iranian viruses are most closely related and should be individualized as DASHV (sublineage Ib). Although discovered first in the 1940's, SFSV is still listed as "tentative species" by the International Committee for Taxonomy of Viruses. Based on the results described in this study, we propose to approve Sandfly fever Sicilian virus species. Owing that most of these viruses have been identified in human patients with febrile illness, DASHV should be considered as a potential human pathogen in Iran.