Search

Abstract 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 (~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.

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., (© The Author(s) 2024. Published by Oxford University Press on behalf of the European Society of Evolutionary Biology. All rights reserved. For commercial re-use, please contact reprints@oup.com for reprints and translation rights for reprints. All other permissions can be obtained through our RightsLink service via the Permissions link on the article page on our site—for further information please contact journals.permissions@oup.com.)

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., (© 2024 The Author(s). Molecular Ecology published by John Wiley & Sons Ltd.)

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., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 Elsevier Ltd. All rights reserved.)

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., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Kuhn, Valentin, Stojanovic, Strobl, Babushku, Wang, Rambold, Scheffler, Grath, John-Robbert, Blum, Feuchtinger, Blutke, Weih, Kitamura, Rad, Strobl and Zimber-Strobl.)

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., (© The Author(s) 2021. Published by Oxford University Press on behalf of Genetics Society of America.)

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., (© The Author(s) 2021. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

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., (© The Author(s) 2021. Published by Oxford University Press.)

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., (© The Author(s) 2020. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

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., (© The Author(s) 2020. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

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 ., (Copyright © 2019 Koniger et al.)

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., Competing Interests: The authors declare no conflict of interest.

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., (© 2018 by The American Society of Hematology.)

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.

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 (+) /CD45 dim /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., Competing Interests: CONFLICTS OF INTEREST KSG received honoraria from Celgene. All other authors declare no potential conflicts of interest.

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.

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-string–based search methods.

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., (Copyright © 2013 Elsevier B.V. All rights reserved.)

The phylogeny of insects, one of the most spectacular radiations of life on earth, has received considerable attention. However, the evolutionary roots of one intriguing group of insects, the twisted-wing parasites (Strepsiptera), remain unclear despite centuries of study and debate. Strepsiptera exhibit exceptional larval developmental features, consistent with a predicted step from direct (hemimetabolous) larval development to complete metamorphosis that could have set the stage for the spectacular radiation of metamorphic (holometabolous) insects. Here we report the sequencing of a Strepsiptera genome and show that the analysis of sequence-based genomic data (comprising more than 18 million nucleotides from nearly 4,500 genes obtained from a total of 13 insect genomes), along with genomic metacharacters, clarifies the phylogenetic origin of Strepsiptera and sheds light on the evolution of holometabolous insect development. Our results provide overwhelming support for Strepsiptera as the closest living relatives of beetles (Coleoptera). They demonstrate that the larval developmental features of Strepsiptera, reminiscent of those of hemimetabolous insects, are the result of convergence. Our analyses solve the long-standing enigma of the evolutionary roots of Strepsiptera and reveal that the holometabolous mode of insect development is more malleable than previously thought., (Copyright © 2012 Elsevier Ltd. All rights reserved.)

Plant genomes are generally very large, mostly paleopolyploid, and have numerous gene duplicates and complex genomic features such as repeats and transposable elements. Many of these features have been hypothesized to enable plants, which cannot easily escape environmental challenges, to rapidly adapt. Another mechanism, which has recently been well described as a major facilitator of rapid adaptation in bacteria, animals, and fungi but not yet for plants, is modular rearrangement of protein-coding genes. Due to the high precision of profile-based methods, rearrangements can be well captured at the protein level by characterizing the emergence, loss, and rearrangements of protein domains, their structural, functional, and evolutionary building blocks. Here, we study the dynamics of domain rearrangements and explore their adaptive benefit in 27 plant and 3 algal genomes. We use a phylogenomic approach by which we can explain the formation of 88% of all arrangements by single-step events, such as fusion, fission, and terminal loss of domains. We find many domains are lost along every lineage, but at least 500 domains are novel, that is, they are unique to green plants and emerged more or less recently. These novel domains duplicate and rearrange more readily within their genomes than ancient domains and are overproportionally involved in stress response and developmental innovations. Novel domains more often affect regulatory proteins and show a higher degree of structural disorder than ancient domains. Whereas a relatively large and well-conserved core set of single-domain proteins exists, long multi-domain arrangements tend to be species-specific. We find that duplicated genes are more often involved in rearrangements. Although fission events typically impact metabolic proteins, fusion events often create new signaling proteins essential for environmental sensing. Taken together, the high volatility of single domains and complex arrangements in plant genomes demonstrate the importance of modularity for environmental adaptability of plants.

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.

Genes with sexually dimorphic expression (sex-biased genes) often evolve rapidly and are thought to make an important contribution to reproductive isolation between species. We examined the molecular evolution of sex-biased genes in Drosophila melanogaster and D. ananassae, which represent two independent lineages within the melanogaster group. We find that strong purifying selection limits protein sequence variation within species, but that a considerable fraction of divergence between species can be attributed to positive selection. In D. melanogaster, the proportion of adaptive substitutions between species is greatest for male-biased genes and is especially high for those on the X chromosome. In contrast, male-biased genes do not show unusually high variation within or between populations. A similar pattern is seen at the level of gene expression, where sex-biased genes show high expression divergence between species, but low divergence between populations. In D. ananassae, there is no increased rate of adaptation of male-biased genes, suggesting that the type or strength of selection acting on sex-biased genes differs between lineages.

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.

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.

Phages infecting Lactococcus lactis, a Gram-positive bacterium, are a recurrent problem in the dairy industry. Despite their economical importance, the knowledge on these phages, belonging mostly to Siphoviridae, lags behind that accumulated for members of Myoviridae. The three-dimensional structures of the receptor-binding proteins (RBP) of three lactococcal phages have been determined recently, illustrating their modular assembly and assigning the nature of their bacterial receptor. These RBPs are attached to the baseplate, a large phage organelle, located at the tip of the tail. Tuc2009 baseplate is formed by the products of 6 open read frames, including the RBP. Because phage binding to its receptor induces DNA release, it has been postulated that the baseplate might be the trigger for DNA injection. We embarked on a structural study of the lactococcal phages baseplate, ultimately to gain insight into the triggering mechanism following receptor binding. Structural features of the Tuc2009 baseplate were established using size exclusion chromatography coupled to on-line UV-visible absorbance, light scattering, and refractive index detection (MALS/UV/RI). Combining the results of this approach with literature data led us to propose a "low resolution" model of Tuc2009 baseplate. This model will serve as a knowledge base to submit relevant complexes to crystallization trials.

Since the time bacteriophages were first identified as a major cause of fermentation failure in the dairy industry, researchers have been struggling to develop strategies to exclude them from the dairy environment. Over 70 years of research has led to huge improvements in the consistency and quality of fermented dairy products, while also facilitating an appreciation of the beneficial properties of bacteriophages with respect to dairy product development. With specific reference to Lactococcus lactis and cheese production, this review outlines some recently reported novel methods aimed at limiting the bacteriophage infection as well as highlighting some beneficial aspects of bacteriophage activity.

The complete genome sequences of three lactococcal 936-type bacteriophages, 712, jj50 and P008, were determined. Comparative genomic analysis of these phages with the previously sequenced 936-type phages, sk1 and bIL170, reveals a strict conservation of the overall genetic organization of this geographically diverse phage group. Genetic divergence was mainly observed in the early expressed region of the phage genomes, where a number of deletions, exchanges and insertions appear to have occurred. These genetic differences may be responsible for the observed differential sensitivity to the lactococcal DNA injection blocking protein, Sie(2009), and the abortive infection system, AbiA.

Bacteriophages of the Siphoviridae family utilize a long noncontractile tail to recognize, adsorb to, and inject DNA into their bacterial host. The tail anatomy of the archetypal Siphoviridae lambda has been well studied, in contrast to phages infecting gram-positive bacteria. This report outlines a detailed anatomical description of a typical member of the Siphoviridae infecting a gram-positive bacterium. The tail superstructure of the lactococcal phage Tuc2009 was investigated using N-terminal protein sequencing, Western blotting, and immunogold transmission electron microscopy, allowing a tangible path to be followed from gene sequence through encoded protein to specific architectural structures on the Tuc2009 virion. This phage displays a striking parity with lambda with respect to tail structure, which reenforced a model proposed for Tuc2009 tail architecture. Furthermore, comparisons with lambda and other lactococcal phages allowed the specification of a number of genetic submodules likely to encode specific tail structures.

The first step in the infection process of tailed phages is recognition and binding to the host receptor. This interaction is mediated by the phage antireceptor located in the distal tail structure. The temperate Lactococcus lactis phage TP901-1 belongs to the P335 species of the Siphoviridae family, which also includes the related phage Tuc2009. The distal tail structure of TP901-1 is well characterized and contains a double-disk baseplate and a central tail fiber. The structural tail proteins of TP901-1 and Tuc2009 are highly similar, but the phages have different host ranges and must therefore encode different antireceptors. In order to identify the antireceptors of TP901-1 and Tuc2009, a chimeric phage was generated in which the gene encoding the TP901-1 lower baseplate protein (bppL(TP901-1)) was exchanged with the analogous gene (orf53(2009)) of phage Tuc2009. The chimeric phage (TP901-1C) infected the Tuc2009 host strain efficiently and thus displayed an altered host range compared to TP901-1. Genomic analysis and sequencing verified that TP901-1C is a TP901-1 derivative containing the orf53(2009) gene in exchange for bppL(TP901-1); however, a new sequence in the late promoter region was also discovered. Protein analysis confirmed that TP901-1C contains ORF53(2009) and not the lower baseplate protein BppL(TP901-1), and it was concluded that BppL(TP901-1) and ORF53(2009) constitute antireceptor proteins of TP901-1 and Tuc2009, respectively. Electron micrographs revealed altered baseplate morphology of TP901-1C compared to that of the parental phage.

The tail structures of bacteriophages infecting gram-positive bacteria are largely unexplored, although the phage tail mediates the initial interaction with the host cell. The temperate Lactococcus lactis phage TP901-1 of the Siphoviridae family has a long noncontractile tail with a distal baseplate. In the present study, we investigated the distal tail structures and tail assembly of phage TP901-1 by introducing nonsense mutations into the late transcribed genes dit (orf46), tal(TP901-1) (orf47), bppU (orf48), bppL (orf49), and orf50. Transmission electron microscopy examination of mutant and wild-type TP901-1 phages showed that the baseplate consisted of two different disks and that a central tail fiber is protruding below the baseplate. Evaluation of the mutant tail morphologies with protein profiles and Western blots revealed that the upper and lower baseplate disks consist of the proteins BppU and BppL, respectively. Likewise, Dit and Tal(TP901-1) were shown to be structural tail proteins essential for tail formation, and Tal(TP901-1) was furthermore identified as the tail fiber protein by immunogold labeling experiments. Determination of infection efficiencies of the mutant phages showed that the baseplate is fundamental for host infection and the lower disk protein, BppL, is suggested to interact with the host receptor. In contrast, ORF50 was found to be nonessential for tail assembly and host infection. A model for TP901-1 tail assembly, in which the function of eight specific proteins is considered, is presented.

The genome of bacteriophage Tuc2009 consists of 38347 base pairs on which 57 open reading frames (ORFs) were identified, divided in two oppositely transcribed regions. The leftward-transcribed region harbors three ORFs, two of which are involved in the establishment of lysogeny. The rightward-transcribed region contains 54 ORFs, which are assumed to be required for the lytic life cycle. An exception to the above organization is ORF 10, of unknown function, located within the rightward-transcribed region that has an orientation opposite to the ORFs surrounding it. Transcriptional analysis of the Tuc2009 genome following infection of a sensitive host revealed that most ORFs are transcribed in a sequential manner. ORFs that are presumed to form (part of) the genetic switch along with the superinfection exclusion-encoding gene are transcribed immediately after infection, followed by transcription of the presumed replication region. Subsequent to this, several small transcripts could be identified followed by a single 24-kb transcript. This latter transcript was shown to specify most of the identified structural proteins as well as two proteins required for host lysis. Interestingly, the 24-kb mRNA was shown to undergo splicing through the activity of a type I intron whose removal from the mRNA resulted in the formation of an ORF specifying a major structural protein. Primer extension analysis was employed to identify the 5' ends of mRNA transcripts and the genome and transcriptional data are discussed in relation to other lactococcal bacteriophages.