RESUMEN
Methods of transcriptional profiling have made it possible to compare gene expression between females and males on a genome-wide scale. Such studies have revealed that sex-biased gene expression is abundant in many species, although its extent may vary greatly among tissues or developmental stages. In species with genetic sex determination, sex chromosome-specific processes, such as dosage compensation, also may influence sex-biased gene expression. Sex-biased genes, especially those with male-biased expression, often show elevated rates of both protein sequence and gene expression divergence between species, which could have a number of causes, including sexual selection, sexual antagonism, and relaxed selective constraint. Here, we review our current knowledge of sex-biased gene expression in both model and nonmodel organisms, as well as the biological and technical factors that should be considered when analyzing sex-biased expression. We also discuss current approaches to uncover the evolutionary forces that govern the evolution of sex-biased genes.
Asunto(s)
Compensación de Dosificación (Genética) , Regulación de la Expresión Génica , Cromosomas Sexuales , Animales , Evolución Biológica , Epigénesis Genética , Evolución Molecular , Femenino , Humanos , Masculino , Selección Genética , Caracteres Sexuales , Procesos de Determinación del SexoRESUMEN
RNA sequencing (RNAseq) methodology has experienced a burst of technological developments in the last decade, which has opened up opportunities for studying the mechanisms of adaptation to environmental factors at both the organismal and cellular level. Selecting the most suitable experimental approach for specific research questions and model systems can, however, be a challenge and researchers in ecology and evolution are commonly faced with the choice of whether to study gene expression variation in whole bodies, specific tissues, and/or single cells. A wide range of sometimes polarised opinions exists over which approach is best. Here, we highlight the advantages and disadvantages of each of these approaches to provide a guide to help researchers make informed decisions and maximise the power of their study. Using illustrative examples of various ecological and evolutionary research questions, we guide the readers through the different RNAseq approaches and help them identify the most suitable design for their own projects.
RESUMEN
Oak gall wasps typically exhibit a life cycle with one sexual and one asexual generation each year. These wasps can carry various endosymbionts, one of which is the maternally inherited bacterium Wolbachia that can induce several reproductive manipulations on its host. Cytoplasmic incompatibility (CI) has been described as the most prominent of these manipulations. CI leads to embryonic mortality in the hosts' offspring when infected males mate with either uninfected females or with females that harbour different Wolbachia strains. It has been hypothesized that Wolbachia can induce CI in oak gall wasps. To address this hypothesis, we derived a mathematical model to investigate the spread of a bacterial infection in naive populations and to determine the plausibility of CI occurrence. To validate our model, we used published data from Wolbachia-infected Belonocnema kinseyi populations in two approaches. Our first approach uses measurements of infection frequencies and maternal transmission in the sexual generation. For the second approach, we extended the model to compare predictions to estimates of mtDNA-haplotypes, which, like Wolbachia, are maternally inherited, and can therefore be associated with the infection. Both approaches indicate that CI is present in these populations. Our model can be generalized to investigate the occurrence of CI not only for oak gall wasps but also for other species.
Asunto(s)
Partenogénesis , Avispas , Wolbachia , Wolbachia/fisiología , Animales , Avispas/microbiología , Avispas/fisiología , Femenino , Masculino , Quercus/microbiología , Modelos Biológicos , Simbiosis , CitoplasmaRESUMEN
Drosophila melanogaster is a leading model in population genetics and genomics, and a growing number of whole-genome data sets from natural populations of this species have been published over the last years. A major challenge is the integration of disparate data sets, often generated using different sequencing technologies and bioinformatic pipelines, which hampers our ability to address questions about the evolution of this species. Here we address these issues by developing a bioinformatics pipeline that maps pooled sequencing (Pool-Seq) reads from D. melanogaster to a hologenome consisting of fly and symbiont genomes and estimates allele frequencies using either a heuristic (PoolSNP) or a probabilistic variant caller (SNAPE-pooled). We use this pipeline to generate the largest data repository of genomic data available for D. melanogaster to date, encompassing 271 previously published and unpublished population samples from over 100 locations in >20 countries on four continents. Several of these locations have been sampled at different seasons across multiple years. This data set, which we call Drosophila Evolution over Space and Time (DEST), is coupled with sampling and environmental metadata. A web-based genome browser and web portal provide easy access to the SNP data set. We further provide guidelines on how to use Pool-Seq data for model-based demographic inference. Our aim is to provide this scalable platform as a community resource which can be easily extended via future efforts for an even more extensive cosmopolitan data set. Our resource will enable population geneticists to analyze spatiotemporal genetic patterns and evolutionary dynamics of D. melanogaster populations in unprecedented detail.
Asunto(s)
Drosophila melanogaster , Metagenómica , Animales , Drosophila melanogaster/genética , Frecuencia de los Genes , Genética de Población , GenómicaRESUMEN
Genetic variation is the fuel of evolution, with standing genetic variation especially important for short-term evolution and local adaptation. To date, studies of spatiotemporal patterns of genetic variation in natural populations have been challenging, as comprehensive sampling is logistically difficult, and sequencing of entire populations costly. Here, we address these issues using a collaborative approach, sequencing 48 pooled population samples from 32 locations, and perform the first continent-wide genomic analysis of genetic variation in European Drosophila melanogaster. Our analyses uncover longitudinal population structure, provide evidence for continent-wide selective sweeps, identify candidate genes for local climate adaptation, and document clines in chromosomal inversion and transposable element frequencies. We also characterize variation among populations in the composition of the fly microbiome, and identify five new DNA viruses in our samples.
Asunto(s)
Drosophila melanogaster/genética , Genoma de los Insectos , Variación Estructural del Genoma , Microbiota , Selección Genética , Aclimatación/genética , Altitud , Animales , Virus ADN , Drosophila melanogaster/virología , Europa (Continente) , Genoma Mitocondrial , Haplotipos , Virus de Insectos , Masculino , Filogeografía , Polimorfismo de Nucleótido SimpleRESUMEN
MOTIVATION: Homology search methods are dominated by the central paradigm that sequence similarity is a proxy for common ancestry and, by extension, functional similarity. For determining sequence similarity in proteins, most widely used methods use models of sequence evolution and compare amino-acid strings in search for conserved linear stretches. Probabilistic models or sequence profiles capture the position-specific variation in an alignment of homologous sequences and can identify conserved motifs or domains. While profile-based search methods are generally more accurate than simple sequence comparison methods, they tend to be computationally more demanding. In recent years, several methods have emerged that perform protein similarity searches based on domain composition. However, few methods have considered the linear arrangements of domains when conducting similarity searches, despite strong evidence that domain order can harbour considerable functional and evolutionary signal. RESULTS: Here, we introduce an alignment scheme that uses a classical dynamic programming approach to the global alignment of domains. We illustrate that representing proteins as strings of domains (domain arrangements) and comparing these strings globally allows for a both fast and sensitive homology search. Further, we demonstrate that the presented methods complement existing methods by finding similar proteins missed by popular amino-acid-based comparison methods. AVAILABILITY: An implementation of the presented algorithms, a web-based interface as well as a command-line program for batch searching against the UniProt database can be found at http://rads.uni-muenster.de. Furthermore, we provide a JAVA API for programmatic access to domain-stringbased search methods.
Asunto(s)
Algoritmos , Biología Computacional , Dominios y Motivos de Interacción de Proteínas , Proteínas/química , Alineación de Secuencia/métodos , Análisis de Secuencia de Proteína/métodos , Programas Informáticos , Proteínas Activadoras de GTPasa/química , Modelos Estadísticos , Fosfotransferasas (Fosfomutasas)/química , Curva ROC , Saccharomyces cerevisiaeRESUMEN
Modularity is a hallmark of molecular evolution. Whether considering gene regulation, the components of metabolic pathways or signaling cascades, the ability to reuse autonomous modules in different molecular contexts can expedite evolutionary innovation. Similarly, protein domains are the modules of proteins, and modular domain rearrangements can create diversity with seemingly few operations in turn allowing for swift changes to an organism's functional repertoire. Here, we assess the patterns and functional effects of modular rearrangements at high resolution. Using a well resolved and diverse group of pancrustaceans, we illustrate arrangement diversity within closely related organisms, estimate arrangement turnover frequency and establish, for the first time, branch-specific rate estimates for fusion, fission, domain addition and terminal loss. Our results show that roughly 16 new arrangements arise per million years and that between 64% and 81% of these can be explained by simple, single-step modular rearrangement events. We find evidence that the frequencies of fission and terminal deletion events increase over time, and that modular rearrangements impact all levels of the cellular signaling apparatus and thus may have strong adaptive potential. Novel arrangements that cannot be explained by simple modular rearrangements contain a significant amount of repeat domains that occur in complex patterns which we term "supra-repeats". Furthermore, these arrangements are significantly longer than those with a single-step rearrangement solution, suggesting that such arrangements may result from multi-step events. In summary, our analysis provides an integrated view and initial quantification of the patterns and functional impact of modular protein evolution in a well resolved phylogenetic tree. This article is part of a Special Issue entitled: The emerging dynamic view of proteins: Protein plasticity in allostery, evolution and self-assembly.
Asunto(s)
Evolución Biológica , Filogenia , Proteínas/genéticaRESUMEN
Ectotherms such as insects are animals whose body temperature largely depends on ambient temperature and temperature variations provide a selection pressure affecting the geographical distribution of these species. However, over the course of evolution, some insect species managed to colonize environments characterized by various temperature ranges. Therefore, insects provide an excellent study system to investigate the basis of adaptation to temperature changes and extremes. We are generally using the vinegar fly Drosophila ananassae as a model system to investigate the genetic basis of cold tolerance. This species has expanded from its tropical ancestral range to more temperate regions resulting in a cosmopolitan, domestic distribution. Previously, we identified candidate genes significantly associated with cold tolerance in this species. We now established molecular genetic tools to assess the function of these genes. Using CRISPR/Cas9 methodology for genome editing and the PiggyBac system, the Cas9 enzyme was successfully integrated into the genome of three fly strains with different levels of cold tolerance. We further report on preliminary findings that the Cas9 integration itself did not have a consistent effect on tolerance to cold. In conclusion, we offer with our study the molecular tools that allow studying stress-related candidate genes in D. ananassae in the future. In addition, we point out and provide guidance on the challenges that come with genome editing in a non-model species.
Asunto(s)
Sistemas CRISPR-Cas , Drosophila , Animales , Drosophila/genética , Temperatura , Frío , AclimataciónRESUMEN
Broadly distributed species must cope with diverse and changing environmental conditions, including various forms of stress. Cosmopolitan populations of Drosophila melanogaster are more tolerant to oxidative stress than those from the species' ancestral range in sub-Saharan Africa, and the degree of tolerance is associated with an insertion/deletion polymorphism in the 3' untranslated region of the Metallothionein A (MtnA) gene that varies clinally in frequency. We examined oxidative stress tolerance and the transcriptional response to oxidative stress in cosmopolitan and sub-Saharan African populations of D. melanogaster, including paired samples with allelic differences at the MtnA locus. We found that the effect of the MtnA polymorphism on oxidative stress tolerance was dependent on the genomic background, with the deletion allele increasing tolerance only in a northern, temperate population. Genes that were differentially expressed under oxidative stress included MtnA and other metallothioneins, as well as those involved in glutathione metabolism and other genes known to be part of the oxidative stress response or the general stress response. A gene coexpression analysis revealed further genes and pathways that respond to oxidative stress including those involved in additional metabolic processes, autophagy, and apoptosis. There was a significant overlap among the genes induced by oxidative and cold stress, which suggests a shared response pathway to these two stresses. Interestingly, the MtnA deletion was associated with consistent changes in the expression of many genes across all genomic backgrounds, regardless of the expression level of the MtnA gene itself. We hypothesize that this is an indirect effect driven by the loss of microRNA binding sites within the MtnA 3' untranslated region.
Asunto(s)
Drosophila melanogaster , Metalotioneína , Animales , Drosophila melanogaster/metabolismo , Mutación INDEL , Metalotioneína/genética , Metalotioneína/metabolismo , Metalotioneína/farmacología , Estrés Oxidativo/genética , Polimorfismo GenéticoRESUMEN
Activation of CD40-signaling contributes to the initiation, progression and drug resistance of B cell lymphomas. We contributed to this knowledge by showing that constitutive CD40-signaling in B cells induces B cell hyperplasia and finally B cell lymphoma development in transgenic mice. CD40 activates, among others, the non-canonical NF-ĸB signaling, which is constitutively activated in several human B cell lymphomas and is therefore presumed to contribute to lymphopathogenesis. This prompted us to study the regulatory role of the non-canonical NF-ĸB transcription factor RelB in lymphomagenesis. To this end, we crossed mice expressing a constitutively active CD40 receptor in B cells with conditional RelB-KO mice. Ablation of RelB attenuated pre-malignant B cell expansion, and resulted in an impaired survival and activation of long-term CD40-stimulated B cells. Furthermore, we found that hyperactivation of non-canonical NF-кB signaling enhances the retention of B cells in the follicles of secondary lymphoid organs. RNA-Seq-analysis revealed that several genes involved in B-cell migration, survival, proliferation and cytokine signaling govern the transcriptional differences modulated by the ablation of RelB in long-term CD40-stimulated B cells. Inactivation of RelB did not abrogate lymphoma development. However, lymphomas occurred with a lower incidence and had a longer latency period. In summary, our data suggest that RelB, although it is not strictly required for malignant transformation, accelerates the lymphomagenesis of long-term CD40-stimulated B cells by regulating genes involved in migration, survival and cytokine signaling.
Asunto(s)
Linfoma de Células B , Linfoma , Factor de Transcripción ReIB , Animales , Linfocitos B , Antígenos CD40/genética , Citocinas , Humanos , Linfoma de Células B/genética , Ratones , Ratones Transgénicos , FN-kappa B , Factor de Transcripción ReIB/genéticaRESUMEN
Drosophila melanogaster is an important model for antiviral immunity in arthropods, but very few DNA viruses have been described from the family Drosophilidae. This deficiency limits our opportunity to use natural host-pathogen combinations in experimental studies, and may bias our understanding of the Drosophila virome. Here, we report fourteen DNA viruses detected in a metagenomic analysis of 6668 pool-sequenced Drosophila, sampled from forty-seven European locations between 2014 and 2016. These include three new nudiviruses, a new and divergent entomopoxvirus, a virus related to Leptopilina boulardi filamentous virus, and a virus related to Musca domestica salivary gland hypertrophy virus. We also find an endogenous genomic copy of galbut virus, a double-stranded RNA partitivirus, segregating at very low frequency. Remarkably, we find that Drosophila Vesanto virus, a small DNA virus previously described as a bidnavirus, may be composed of up to twelve segments and thus represent a new lineage of segmented DNA viruses. Two of the DNA viruses, Drosophila Kallithea nudivirus and Drosophila Vesanto virus are relatively common, found in 2 per cent or more of wild flies. The others are rare, with many likely to be represented by a single infected fly. We find that virus prevalence in Europe reflects the prevalence seen in publicly available datasets, with Drosophila Kallithea nudivirus and Drosophila Vesanto virus the only ones commonly detectable in public data from wild-caught flies and large population cages, and the other viruses being rare or absent. These analyses suggest that DNA viruses are at lower prevalence than RNA viruses in D.melanogaster, and may be less likely to persist in laboratory cultures. Our findings go some way to redressing an earlier bias toward RNA virus studies in Drosophila, and lay the foundation needed to harness the power of Drosophila as a model system for the study of DNA viruses.
RESUMEN
The tremendous diversity of Hymenoptera is commonly attributed to the evolution of parasitoidism in the last common ancestor of parasitoid sawflies (Orussidae) and wasp-waisted Hymenoptera (Apocrita). However, Apocrita and Orussidae differ dramatically in their species richness, indicating that the diversification of Apocrita was promoted by additional traits. These traits have remained elusive due to a paucity of sawfly genome sequences, in particular those of parasitoid sawflies. Here, we present comparative analyses of draft genomes of the primarily phytophagous sawfly Athalia rosae and the parasitoid sawfly Orussus abietinus. Our analyses revealed that the ancestral hymenopteran genome exhibited traits that were previously considered unique to eusocial Apocrita (e.g., low transposable element content and activity) and a wider gene repertoire than previously thought (e.g., genes for CO2 detection). Moreover, we discovered that Apocrita evolved a significantly larger array of odorant receptors than sawflies, which could be relevant to the remarkable diversification of Apocrita by enabling efficient detection and reliable identification of hosts.
Asunto(s)
Especiación Genética , Genoma de los Insectos , Interacciones Huésped-Parásitos/genética , Himenópteros/genética , Secuencia de Aminoácidos , Animales , Secuencia Conservada , Elementos Transponibles de ADN , Femenino , Dosificación de Gen , Glicoproteínas/genética , Herbivoria/genética , Inmunidad/genética , Proteínas de Insectos/genética , Masculino , Familia de Multigenes , Receptores Odorantes/genética , Conducta Social , Visión Ocular/genéticaRESUMEN
BACKGROUND: Genes with sex-biased expression often show rapid molecular evolution between species. Previous population genetic and comparative genomic studies of Drosophila melanogaster and D. simulans revealed that male-biased genes have especially high rates of adaptive evolution. To test if this is also the case for other lineages within the melanogaster group, we investigated gene expression in D. ananassae, a species that occurs in structured populations in tropical and subtropical regions. We used custom-made microarrays and published microarray data to characterize the sex-biased expression of 129 D. ananassae genes whose D. melanogaster orthologs had been classified previously as male-biased, female-biased, or unbiased in their expression and had been studied extensively at the population-genetic level. For 43 of these genes we surveyed DNA sequence polymorphism in a natural population of D. ananassae and determined divergence to the sister species D. atripex and D. phaeopleura. RESULTS: Sex-biased expression is generally conserved between D. melanogaster and D. ananassae, with the majority of genes exhibiting the same bias in the two species. However, about one-third of the genes have either gained or lost sex-biased expression in one of the species and a small proportion of genes (approximately 4%) have changed bias from one sex to the other. The male-biased genes of D. ananassae show evidence of positive selection acting at the protein level. However, the signal of adaptive protein evolution for male-biased genes is not as strong in D. ananassae as it is in D. melanogaster and is limited to genes with conserved male-biased expression in both species. Within D. ananassae, a significant signal of adaptive evolution is also detected for female-biased and unbiased genes. CONCLUSIONS: Our findings extend previous observations of widespread adaptive protein evolution to an independent Drosophila lineage, the D. ananassae subgroup. However, the rate of adaptive evolution is not greater for male-biased genes than for female-biased or unbiased genes, which suggests that there are differences in sex-biased gene evolution between the two lineages.
Asunto(s)
Drosophila melanogaster/genética , Drosophila/genética , Evolución Molecular , Genes de Insecto , Animales , Drosophila/metabolismo , Drosophila melanogaster/metabolismo , Femenino , Masculino , Análisis de Secuencia por Matrices de Oligonucleótidos , Filogenia , Polimorfismo Genético , Caracteres Sexuales , Especificidad de la EspecieRESUMEN
Ectothermic species such as insects are particularly vulnerable to climatic fluctuations. Nevertheless, many insects that evolved and diversified in the tropics have successfully colonized temperate regions all over the globe. To shed light on the genetic basis of cold tolerance in such species, we conducted a quantitative trait locus (QTL) mapping experiment for chill coma recovery time (CCRT) in Drosophila ananassae, a cosmopolitan species that has expanded its range from tropical to temperate regions. We created a mapping population of recombinant inbred advanced intercross lines (RIAILs) from two founder strains with diverging CCRT phenotypes. The RIAILs were phenotyped for their CCRT and, together with the founder strains, genotyped for polymorphic markers with double-digest restriction site-associated DNA (ddRAD) sequencing. Using a hierarchical mapping approach that combined standard interval mapping and a multiple-QTL model, we mapped three QTL which altogether explained 64% of the phenotypic variance. For two of the identified QTL, we found evidence of epistasis. To narrow down the list of cold tolerance candidate genes, we cross-referenced the QTL intervals with genes that we previously identified as differentially expressed in response to cold in D. ananassae, and with thermotolerance candidate genes of D. melanogaster Among the 58 differentially expressed genes that were contained within the QTL, GF15058 showed a significant interaction of the CCRT phenotype and gene expression. Further, we identified the orthologs of four D. melanogaster thermotolerance candidate genes, MtnA, klarsicht, CG5246 (D.ana/GF17132) and CG10383 (D.ana/GF14829) as candidates for cold tolerance in D. ananassae.
Asunto(s)
Respuesta al Choque por Frío/genética , Drosophila/genética , Sitios de Carácter Cuantitativo , Animales , Mapeo Cromosómico , Frío , Genes de Insecto , Variación Genética , FenotipoRESUMEN
Prezygotic mating isolation has been a major interest of evolutionary biologists during the past several decades because it is likely to represent one of the first stages in the transition from populations to species. Mate discrimination is one of the most commonly measured forms of prezygotic isolation and appears to be relatively common among closely related species. In some cases, it has been used as a measure to distinguish populations from subspecies, races, and sister species, yet the influences of various evolutionary mechanisms that may generate mate discrimination are largely unknown. In this study, we measured the level and pattern of mate discrimination among 18 populations of a cosmopolitan drosophilid species, Drosophila ananassae, from throughout its geographical range and its sister species, Drosophila pallidosa, which has a restricted geographical distribution in the South Pacific Islands. In addition, we measured genetic differentiation between all 18 populations using mitochondrial DNA polymorphism data. Mate discrimination varies considerably throughout the species range, being higher among populations outside the ancestral Indonesian range, and highest in the South Pacific. Our results suggest that colonization and genetic differentiation may have an influence on the evolutionary origin of mate discrimination. Our phylogeographical approach clarifies the ancestral relationships of several populations from the South Pacific that show particularly strong mate discrimination and suggests that they may be in the early stages of speciation. Furthermore, both the genetic and behavioral results cast doubt on the status of D. pallidosa as a good species.
Asunto(s)
Drosophila/genética , Polimorfismo Genético , Conducta Sexual Animal , Animales , Asia , Australia , Brasil , ADN Mitocondrial/genética , Drosophila/fisiología , Evolución Molecular , Femenino , Geografía , Masculino , Datos de Secuencia Molecular , Islas del Pacífico , Filogenia , Análisis de Secuencia de ADNRESUMEN
Coping with daily and seasonal temperature fluctuations is a key adaptive process for species to colonize temperate regions all over the globe. Over the past 18,000 years, the tropical species Drosophila ananassae expanded its home range from tropical regions in Southeast Asia to more temperate regions. Phenotypic assays of chill coma recovery time (CCRT) together with previously published population genetic data suggest that only a small number of genes underlie improved cold hardiness in the cold-adapted populations. We used high-throughput RNA sequencing to analyze differential gene expression before and after exposure to a cold shock in coldtolerant lines (those with fast chill coma recovery, CCR) and cold-sensitive lines (slow CCR) from a population originating from Bangkok, Thailand (the ancestral species range). We identified two candidate genes with a significant interaction between cold tolerance and cold shock treatment: GF14647 and GF15058. Further, our data suggest that selection for increased cold tolerance did not operate through the increased activity of heat shock proteins, but more likely through the stabilization of the actin cytoskeleton and a delayed onset of apoptosis.
RESUMEN
It has been experimentally shown that DNA methylation is involved in the regulation of gene expression and the silencing of transposable element activity in eukaryotes. The variable levels of DNA methylation among different insect species indicate an evolutionarily flexible role of DNA methylation in insects, which due to a lack of comparative data is not yet well-substantiated. Here, we use computational methods to trace signatures of DNA methylation across insects by analyzing transcriptomic and genomic sequence data from all currently recognized insect orders. We conclude that: 1) a functional methylation system relying exclusively on DNA methyltransferase 1 is widespread across insects. 2) DNA methylation has potentially been lost or extremely reduced in species belonging to springtails (Collembola), flies and relatives (Diptera), and twisted-winged parasites (Strepsiptera). 3) Holometabolous insects display signs of reduced DNA methylation levels in protein-coding sequences compared with hemimetabolous insects. 4) Evolutionarily conserved insect genes associated with housekeeping functions tend to display signs of heavier DNA methylation in comparison to the genomic/transcriptomic background. With this comparative study, we provide the much needed basis for experimental and detailed comparative analyses required to gain a deeper understanding on the evolution and function of DNA methylation in insects.
Asunto(s)
Metilación de ADN/genética , Evolución Molecular , Holometabola/genética , Insectos/genética , Animales , Artrópodos/genética , Genoma/genética , Holometabola/metabolismo , Insectos/metabolismo , Filogenia , Alineación de SecuenciaRESUMEN
Mesenchymal stromal cells (MSCs) are crucial components of the bone marrow (BM) microenvironment essential for regulating self-renewal, survival, and differentiation of hematopoietic stem/progenitor cells (HSPCs) in the stem cell niche. MSCs are functionally altered in myelodysplastic syndromes (MDS) and exhibit an altered methylome compared with MSCs from healthy controls, thus contributing to disease progression. To determine whether MSCs are amenable to epigenetic therapy and if this affects their function, we examined growth, differentiation, and HSPC-supporting capacity of ex vivo-expanded MSCs from MDS patients in comparison with age-matched healthy controls after direct treatment in vitro with the hypomethylating agent azacitidine (AZA). Strikingly, we find that AZA exerts a direct effect on healthy as well as MDS-derived MSCs such that they favor support of healthy over malignant clonal HSPC expansion in coculture experiments. RNA-sequencing analyses of MSCs identified stromal networks regulated by AZA. Notably, these comprise distinct molecular pathways crucial for HSPC support, foremost extracellular matrix molecules (including collagens) and interferon pathway components. Our study demonstrates that the hypomethylating agent AZA exerts its antileukemic activity in part through a direct effect on the HSPC-supporting BM niche and provides proof of concept for the therapeutic potential of epigenetic treatment of diseased MSCs. In addition, our comprehensive data set of AZA-sensitive gene networks represents a valuable framework to guide future development of targeted epigenetic niche therapy in myeloid malignancies such as MDS and acute myeloid leukemia.
Asunto(s)
Azacitidina/farmacología , Hematopoyesis/efectos de los fármacos , Adipogénesis/efectos de los fármacos , Adulto , Anciano , Anciano de 80 o más Años , Células de la Médula Ósea/citología , Diferenciación Celular/efectos de los fármacos , Supervivencia Celular/efectos de los fármacos , Femenino , Redes Reguladoras de Genes/efectos de los fármacos , Humanos , Inmunofenotipificación , Masculino , Células Madre Mesenquimatosas/citología , Células Madre Mesenquimatosas/efectos de los fármacos , Células Madre Mesenquimatosas/metabolismo , Persona de Mediana Edad , Síndromes Mielodisplásicos , Osteogénesis/efectos de los fármacos , Factor de Crecimiento Transformador beta1/genética , Factor de Crecimiento Transformador beta1/metabolismo , Adulto JovenRESUMEN
Effectively targeting leukemia-initiating cells (LIC) in FLT3-ITD-mutated acute myeloid leukemia (AML) is crucial for cure. Tyrosine kinase inhibitors (TKI) have limited impact as single agents, failing to eradicate LIC in the bone marrow. Using primary AML samples and a patient-derived xenograft model, we investigated whether combining the FLT3-selective TKI crenolanib with the hypomethylating agent azacitidine (AZA) eliminates FLT3-ITD LIC and whether efficacy of this combination depends on co-existing mutations. Using multiparameter flow cytometry, we show FLT3-ITD occurs within the most primitive Lin-/CD33(+)/CD45dim/CD34+CD38- LIC compartment. Crenolanib alone could not target FLT3-ITD LIC in contact with niche cells while addition of AZA overcame stromal protection resulting in dramatically reduced clonogenic capacity of LIC in vitro and severely impaired engraftment in NSG mice. Strikingly, FLT3-mutated samples harboring TET2 mutations were completely resistant to crenolanib whereas neither NPM1 nor DNMT3A mutations influenced response. Conversely, primary AML LIC harboring either TET2, DNMT3A or NPM1 mutations did not show increased sensitivity to AZA. In summary, resistance of FLT3-ITD LIC to TKI depends on co-existing epigenetic mutations. However, AZA + crenolanib effectively abrogates stromal protection and inhibits survival of FLT3-ITD LIC irrespective of mutations, providing evidence for this combination as a means to suppress residual LIC.
RESUMEN
Sex-biased gene expression (i.e., the differential expression of genes between males and females) is common among sexually reproducing species. However, genes often differ in their sex-bias classification or degree of sex bias between species. There is also an unequal distribution of sex-biased genes (especially male-biased genes) between the X chromosome and the autosomes. We used whole-genome expression data and evolutionary rate estimates for two different Drosophilid lineages, melanogaster and obscura, spanning an evolutionary time scale of around 50 Myr to investigate the influence of sex-biased gene expression and chromosomal location on the rate of molecular evolution. In both lineages, the rate of protein evolution correlated positively with the male/female expression ratio. Genes with highly male-biased expression, genes expressed specifically in male reproductive tissues, and genes with conserved male-biased expression over long evolutionary time scales showed the fastest rates of evolution. An analysis of sex-biased gene evolution in both lineages revealed evidence for a "fast-X" effect in which the rate of evolution was greater for X-linked than for autosomal genes. This pattern was particularly pronounced for male-biased genes. Genes located on the obscura "neo-X" chromosome, which originated from a recent X-autosome fusion, showed rates of evolution that were intermediate between genes located on the ancestral X-chromosome and the autosomes. This suggests that the shift to X-linkage led to an increase in the rate of molecular evolution.