Gene expression differences consistent with water loss reduction underlie desiccation tolerance of natural Drosophila populations.

BACKGROUND: Climate change is one of the main factors shaping the distribution and biodiversity of organisms, among others by greatly altering water availability, thus exposing species and ecosystems to harsh desiccation conditions. However, most of the studies so far have focused on the effects of increased temperature. Integrating transcriptomics and physiology is key to advancing our knowledge on how species cope with desiccation stress, and these studies are still best accomplished in model organisms. RESULTS: Here, we characterized the natural variation of European D. melanogaster populations across climate zones and found that strains from arid regions were similar or more tolerant to desiccation compared with strains from temperate regions. Tolerant and sensitive strains differed not only in their transcriptomic response to stress but also in their basal expression levels. We further showed that gene expression changes in tolerant strains correlated with their physiological response to desiccation stress and with their cuticular hydrocarbon composition, and functionally validated three of the candidate genes identified. Transposable elements, which are known to influence stress response across organisms, were not found to be enriched nearby differentially expressed genes. Finally, we identified several tRNA-derived small RNA fragments that differentially targeted genes in response to desiccation stress. CONCLUSIONS: Overall, our results showed that basal gene expression differences across individuals should be analyzed if we are to understand the genetic basis of differential stress survival. Moreover, tRNA-derived small RNA fragments appear to be relevant across stress responses and allow for the identification of stress-response genes not detected at the transcriptional level.

The genomic basis of copper tolerance in Drosophila is shaped by a complex interplay of regulatory and environmental factors.

BACKGROUND: Escalation in industrialization and anthropogenic activity have resulted in an increase of pollutants released into the environment. Of these pollutants, heavy metals such as copper are particularly concerning due to their bio-accumulative nature. Due to its highly heterogeneous distribution and its dual nature as an essential micronutrient and toxic element, the genetic basis of copper tolerance is likely shaped by a complex interplay of genetic and environmental factors. RESULTS: In this study, we utilized the natural variation present in multiple populations of Drosophila melanogaster collected across Europe to screen for variation in copper tolerance. We found that latitude and the degree of urbanization at the collection sites, rather than any other combination of environmental factors, were linked to copper tolerance. While previously identified copper-related genes were not differentially expressed in tolerant vs. sensitive strains, genes involved in metabolism, reproduction, and protease induction contributed to the differential stress response. Additionally, the greatest transcriptomic and physiological responses to copper toxicity were seen in the midgut, where we found that preservation of gut acidity is strongly linked to greater tolerance. Finally, we identified transposable element insertions likely to play a role in copper stress response. CONCLUSIONS: Overall, by combining genome-wide approaches with environmental association analysis, and functional analysis of candidate genes, our study provides a unique perspective on the genetic and environmental factors that shape copper tolerance in natural D. melanogaster populations and identifies new genes, transposable elements, and physiological traits involved in this complex phenotype.

Population-scale long-read sequencing uncovers transposable elements associated with gene expression variation and adaptive signatures in Drosophila.

High quality reference genomes are crucial to understanding genome function, structure and evolution. The availability of reference genomes has allowed us to start inferring the role of genetic variation in biology, disease, and biodiversity conservation. However, analyses across organisms demonstrate that a single reference genome is not enough to capture the global genetic diversity present in populations. In this work, we generate 32 high-quality reference genomes for the well-known model species D. melanogaster and focus on the identification and analysis of transposable element variation as they are the most common type of structural variant. We show that integrating the genetic variation across natural populations from five climatic regions increases the number of detected insertions by 58%. Moreover, 26% to 57% of the insertions identified using long-reads were missed by short-reads methods. We also identify hundreds of transposable elements associated with gene expression variation and new TE variants likely to contribute to adaptive evolution in this species. Our results highlight the importance of incorporating the genetic variation present in natural populations to genomic studies, which is essential if we are to understand how genomes function and evolve.

Loss of HDAC11 accelerates skeletal muscle regeneration in mice.

Histone deacetylase 11 (HDAC11) is the latest identified member of the histone deacetylase family of enzymes. It is highly expressed in brain, heart, testis, kidney, and skeletal muscle, although its role in these tissues is poorly understood. Here, we investigate for the first time the consequences of HDAC11 genetic impairment on skeletal muscle regeneration, a process principally dependent on its resident stem cells (satellite cells) in coordination with infiltrating immune cells and stromal cells. Our results show that HDAC11 is dispensable for adult muscle growth and establishment of the satellite cell population, while HDAC11 deficiency advances the regeneration process in response to muscle injury. This effect is not caused by differences in satellite cell activation or proliferation upon injury, but rather by an enhanced capacity of satellite cells to differentiate at early regeneration stages in the absence of HDAC11. Infiltrating HDAC11-deficient macrophages could also contribute to this accelerated muscle regenerative process by prematurely producing high levels of IL-10, a cytokine known to promote myoblast differentiation. Altogether, our results show that HDAC11 depletion advances skeletal muscle regeneration and this finding may have potential implications for designing new strategies for muscle pathologies coursing with chronic damage. DATABASE: Data were deposited in NCBI's Gene Expression Omnibus accessible through GEO Series accession number GSE147423.

Transposable elements contribute to the genomic response to insecticides in Drosophila melanogaster.

Most of the genotype-phenotype analyses to date have largely centred attention on single nucleotide polymorphisms. However, transposable element (TE) insertions have arisen as a plausible addition to the study of the genotypic-phenotypic link because of to their role in genome function and evolution. In this work, we investigate the contribution of TE insertions to the regulation of gene expression in response to insecticides. We exposed four Drosophila melanogaster strains to malathion, a commonly used organophosphate insecticide. By combining information from different approaches, including RNA-seq and ATAC-seq, we found that TEs can contribute to the regulation of gene expression under insecticide exposure by rewiring cis-regulatory networks. This article is part of a discussion meeting issue 'Crossroads between transposons and gene regulation'.

Stress response, behavior, and development are shaped by transposable element-induced mutations in Drosophila.

Most of the current knowledge on the genetic basis of adaptive evolution is based on the analysis of single nucleotide polymorphisms (SNPs). Despite increasing evidence for their causal role, the contribution of structural variants to adaptive evolution remains largely unexplored. In this work, we analyzed the population frequencies of 1,615 Transposable Element (TE) insertions annotated in the reference genome of Drosophila melanogaster, in 91 samples from 60 worldwide natural populations. We identified a set of 300 polymorphic TEs that are present at high population frequencies, and located in genomic regions with high recombination rate, where the efficiency of natural selection is high. The age and the length of these 300 TEs are consistent with relatively young and long insertions reaching high frequencies due to the action of positive selection. Besides, we identified a set of 21 fixed TEs also likely to be adaptive. Indeed, we, and others, found evidence of selection for 84 of these reference TE insertions. The analysis of the genes located nearby these 84 candidate adaptive insertions suggested that the functional response to selection is related with the GO categories of response to stimulus, behavior, and development. We further showed that a subset of the candidate adaptive TEs affects expression of nearby genes, and five of them have already been linked to an ecologically relevant phenotypic effect. Our results provide a more complete understanding of the genetic variation and the fitness-related traits relevant for adaptive evolution. Similar studies should help uncover the importance of TE-induced adaptive mutations in other species as well.

Long-Lasting Primed State in Maize Plants: Salicylic Acid and Steroid Signaling Pathways as Key Players in the Early Activation of Immune Responses in Silks.

In the present study, we investigated the induced systemic resistance (ISR) activated by the beneficial fungus Trichoderma atroviride in maize plants, and the early immunological responses triggered after challenge with the ear rot pathogen Fusarium verticillioides. By transcriptional analysis, we were able to identify the gene core set specifically modulated in silks of maize plants expressing ISR. Our results showed that the main transcriptional reprogramming falls into genes involved in five main functional categories: cell structure or cell wall, amino acid and protein metabolism, stress responses, signaling, and transport. Among these ISR-related genes, it is important to highlight novel findings regarding hormone metabolism and signaling. The expression of hormone-dependent genes was in good agreement with the abscisic acid, jasmonic acid, and salicylic acid (SA) levels detected in the plants under study. The experimental design allowed the identification of novel regulatory elements related to a heightened state of defense in silks and suggests that steroids and SA are central components of a master regulatory network controlling the immunity of silks during ISR. The results presented also provide evidence about the molecular mechanisms used by maize silks against F. verticillioides to counteract pathogenic development and host invasion, including pathogenesis-related genes, plant cell-wall reinforcement, fungal cell-wall-degrading enzymes and secondary metabolism.

Genome-wide patterns of local adaptation in Western European Drosophila melanogaster natural populations.

Signatures of spatially varying selection have been investigated both at the genomic and transcriptomic level in several organisms. In Drosophila melanogaster, the majority of these studies have analyzed North American and Australian populations, leading to the identification of several loci and traits under selection. However, several studies based mainly in North American populations showed evidence of admixture that likely contributed to the observed population differentiation patterns. Thus, disentangling demography from selection might be challenging when analyzing these populations. European populations could help identify loci under spatially varying selection provided that no recent admixture from African populations would have occurred. In this work, we individually sequence the genome of 42 European strains collected in populations from contrasting environments: Stockholm (Sweden) and Castellana Grotte (Southern Italy). We found low levels of population structure and no evidence of recent African admixture in these two populations. We thus look for patterns of spatially varying selection affecting individual genes and gene sets. Besides single nucleotide polymorphisms, we also investigated the role of transposable elements in local adaptation. We concluded that European populations are a good dataset to identify candidate loci under spatially varying selection. The analysis of the two populations sequenced in this work in the context of all the available D. melanogaster data allowed us to pinpoint genes and biological processes likely to be relevant for local adaptation. Identifying and analyzing populations with low levels of population structure and admixture should help to disentangle selective from non-selective forces underlying patterns of population differentiation in other species as well.

A Fungal Effector With Host Nuclear Localization and DNA-Binding Properties Is Required for Maize Anthracnose Development.

Plant pathogens have the capacity to manipulate the host immune system through the secretion of effectors. We identified 27 putative effector proteins encoded in the genome of the maize anthracnose pathogen Colletotrichum graminicola that are likely to target the host's nucleus, as they simultaneously contain sequence signatures for secretion and nuclear localization. We functionally characterized one protein, identified as CgEP1. This protein is synthesized during the early stages of disease development and is necessary for anthracnose development in maize leaves, stems, and roots. Genetic, molecular, and biochemical studies confirmed that this effector targets the host's nucleus and defines a novel class of double-stranded DNA-binding protein. We show that CgEP1 arose from a gene duplication in an ancestor of a lineage of monocot-infecting Colletotrichum spp. and has undergone an intense evolution process, with evidence for episodes of positive selection. We detected CgEP1 homologs in several species of a grass-infecting lineage of Colletotrichum spp., suggesting that its function may be conserved across a large number of anthracnose pathogens. Our results demonstrate that effectors targeted to the host nucleus may be key elements for disease development and aid in the understanding of the genetic basis of anthracnose development in maize plants.

Draft Genome Sequences of Mycobacterium setense Type Strain DSM-45070 and the Nonpathogenic Strain Manresensis, Isolated from the Bank of the Cardener River in Manresa, Catalonia, Spain.

We present here the draft genome sequences of two Mycobacterium setense strains. One of them corresponds to the M. setense type strain DSM-45070, originally isolated from a patient with a posttraumatic chronic skin abscess. The other one corresponds to the nonpathogenic M. setense strain Manresensis, isolated from the Cardener River crossing Manresa, Catalonia, Spain. A comparative genomic analysis shows a smaller genome size and fewer genes in M. setense strain Manresensis relative to those of the type strain, and it shows the genome segments unique to each strain.

Natural selection on coding and noncoding DNA sequences is associated with virulence genes in a plant pathogenic fungus.

Natural selection leaves imprints on DNA, offering the opportunity to identify functionally important regions of the genome. Identifying the genomic regions affected by natural selection within pathogens can aid in the pursuit of effective strategies to control diseases. In this study, we analyzed genome-wide patterns of selection acting on different classes of sequences in a worldwide sample of eight strains of the model plant-pathogenic fungus Colletotrichum graminicola. We found evidence of selective sweeps, balancing selection, and positive selection affecting both protein-coding and noncoding DNA of pathogenicity-related sequences. Genes encoding putative effector proteins and secondary metabolite biosynthetic enzymes show evidence of positive selection acting on the coding sequence, consistent with an Arms Race model of evolution. The 5' untranslated regions (UTRs) of genes coding for effector proteins and genes upregulated during infection show an excess of high-frequency polymorphisms likely the consequence of balancing selection and consistent with the Red Queen hypothesis of evolution acting on these putative regulatory sequences. Based on the findings of this work, we propose that even though adaptive substitutions on coding sequences are important for proteins that interact directly with the host, polymorphisms in the regulatory sequences may confer flexibility of gene expression in the virulence processes of this important plant pathogen.

Draft Genome Sequence of Colletotrichum sublineola, a Destructive Pathogen of Cultivated Sorghum.

Colletotrichum sublineola is a filamentous fungus that causes anthracnose disease on sorghum. We report a draft whole-genome shotgun sequence and gene annotation of the nuclear genome of this fungus using Illumina sequencing.

Identification of positive selection in disease response genes within members of the Poaceae.

Millions of years of coevolution between plants and pathogens can leave footprints on their genomes and genes involved on this interaction are expected to show patterns of positive selection in which novel, beneficial alleles are rapidly fixed within the population. Using information about upregulated genes in maize during Colletotrichum graminicola infection and resources available in the Phytozome database, we looked for evidence of positive selection in the Poaceae lineage, acting on protein coding sequences related with plant defense. We found six genes with evidence of positive selection and another eight with sites showing episodic selection. Some of them have already been described as evolving under positive selection, but others are reported here for the first time including genes encoding isocitrate lyase, dehydrogenases, a multidrug transporter, a protein containing a putative leucine-rich repeat and other proteins with unknown functions. Mapping positively selected residues onto the predicted 3-D structure of proteins showed that most of them are located on the surface, where proteins are in contact with other molecules. We present here a set of Poaceae genes that are likely to be involved in plant defense mechanisms and have evidence of positive selection. These genes are excellent candidates for future functional validation.

Lifestyle transitions in plant pathogenic Colletotrichum fungi deciphered by genome and transcriptome analyses.

Colletotrichum species are fungal pathogens that devastate crop plants worldwide. Host infection involves the differentiation of specialized cell types that are associated with penetration, growth inside living host cells (biotrophy) and tissue destruction (necrotrophy). We report here genome and transcriptome analyses of Colletotrichum higginsianum infecting Arabidopsis thaliana and Colletotrichum graminicola infecting maize. Comparative genomics showed that both fungi have large sets of pathogenicity-related genes, but families of genes encoding secreted effectors, pectin-degrading enzymes, secondary metabolism enzymes, transporters and peptidases are expanded in C. higginsianum. Genome-wide expression profiling revealed that these genes are transcribed in successive waves that are linked to pathogenic transitions: effectors and secondary metabolism enzymes are induced before penetration and during biotrophy, whereas most hydrolases and transporters are upregulated later, at the switch to necrotrophy. Our findings show that preinvasion perception of plant-derived signals substantially reprograms fungal gene expression and indicate previously unknown functions for particular fungal cell types.

Plant defense mechanisms are activated during biotrophic and necrotrophic development of Colletotricum graminicola in maize.

Hemibiotrophic plant pathogens first establish a biotrophic interaction with the host plant and later switch to a destructive necrotrophic lifestyle. Studies of biotrophic pathogens have shown that they actively suppress plant defenses after an initial microbe-associated molecular pattern-triggered activation. In contrast, studies of the hemibiotrophs suggest that they do not suppress plant defenses during the biotrophic phase, indicating that while there are similarities between the biotrophic phase of hemibiotrophs and biotrophic pathogens, the two lifestyles are not analogous. We performed transcriptomic, histological, and biochemical studies of the early events during the infection of maize (Zea mays) with Colletotrichum graminicola, a model pathosystem for the study of hemibiotrophy. Time-course experiments revealed that mRNAs of several defense-related genes, reactive oxygen species, and antimicrobial compounds all begin to accumulate early in the infection process and continue to accumulate during the biotrophic stage. We also discovered the production of maize-derived vesicular bodies containing hydrogen peroxide targeting the fungal hyphae. We describe the fungal respiratory burst during host infection, paralleled by superoxide ion production in specific fungal cells during the transition from biotrophy to a necrotrophic lifestyle. We also identified several novel putative fungal effectors and studied their expression during anthracnose development in maize. Our results demonstrate a strong induction of defense mechanisms occurring in maize cells during C. graminicola infection, even during the biotrophic development of the pathogen. We hypothesize that the switch to necrotrophic growth enables the fungus to evade the effects of the plant immune system and allows for full fungal pathogenicity.