Expression of defensin genes across house fly (Musca domestica) life history gives insight into immune system subfunctionalization.

Animals encounter diverse microbial communities throughout their lifetime, which exert varying selection pressures. Antimicrobial peptides (AMPs), which lyse or inhibit microbial growth, are a first line of defense against some of these microbes. Here we examine how developmental variation in microbial exposure has affected the evolution of expression and amino acid sequences of Defensins (an ancient class of AMPs) in the house fly (Musca domestica). The house fly is a well-suited model for this work because it trophically associates with varying microbial communities throughout its life history and its genome contains expanded families of AMPs, including Defensins. We identified two subsets of house fly Defensins: one expressed in larvae or pupae, and the other expressed in adults. The amino acid sequences of these two Defensin subsets form distinct monophyletic clades, and they are located in separate gene clusters in the genome. The adult-expressed Defensins evolve faster than larval/pupal Defensins, consistent with different selection pressures across developmental stages. Our results therefore suggest that varied microbial communities encountered across life history can shape the evolutionary trajectories of immune genes.

New insights into immune genes and other expanded gene families of the house fly, Musca domestica, from an improved whole genome sequence.

The house fly, Musca domestica, is a pest of livestock, transmits pathogens of human diseases, and is a model organism in multiple biological research areas. The first house fly genome assembly was published in 2014 and has been of tremendous use to the community of house fly biologists, but that genome is discontiguous and incomplete by contemporary standards. To improve the house fly reference genome, we sequenced, assembled, and annotated the house fly genome using improved techniques and technologies that were not available at the time of the original genome sequencing project. The new genome assembly is substantially more contiguous and complete than the previous genome. The new genome assembly has a scaffold N50 of 12.46 Mb, which is a 50-fold improvement over the previous assembly. In addition, the new genome assembly is within 1% of the estimated genome size based on flow cytometry, whereas the previous assembly was missing nearly one-third of the predicted genome sequence. The improved genome assembly has much more contiguous scaffolds containing large gene families. To provide an example of the benefit of the new genome, we used it to investigate tandemly arrayed immune gene families. The new contiguous assembly of these loci provides a clearer picture of the regulation of the expression of immune genes, and it leads to new insights into the selection pressures that shape their evolution.

Gene-Level, but Not Chromosome-Wide, Divergence between a Very Young House Fly Proto-Y Chromosome and Its Homologous Proto-X Chromosome.

X and Y chromosomes are usually derived from a pair of homologous autosomes, which then diverge from each other over time. Although Y-specific features have been characterized in sex chromosomes of various ages, the earliest stages of Y chromosome evolution remain elusive. In particular, we do not know whether early stages of Y chromosome evolution consist of changes to individual genes or happen via chromosome-scale divergence from the X. To address this question, we quantified divergence between young proto-X and proto-Y chromosomes in the house fly, Musca domestica. We compared proto-sex chromosome sequence and gene expression between genotypic (XY) and sex-reversed (XX) males. We find evidence for sequence divergence between genes on the proto-X and proto-Y, including five genes with mitochondrial functions. There is also an excess of genes with divergent expression between the proto-X and proto-Y, but the number of genes is small. This suggests that individual proto-Y genes, but not the entire proto-Y chromosome, have diverged from the proto-X. We identified one gene, encoding an axonemal dynein assembly factor (which functions in sperm motility), that has higher expression in XY males than XX males because of a disproportionate contribution of the proto-Y allele to gene expression. The upregulation of the proto-Y allele may be favored in males because of this gene's function in spermatogenesis. The evolutionary divergence between proto-X and proto-Y copies of this gene, as well as the mitochondrial genes, is consistent with selection in males affecting the evolution of individual genes during early Y chromosome evolution.

Population Genomics Reveals Incipient Speciation, Introgression, and Adaptation in the African Mona Monkey (Cercopithecus mona).

Guenons (tribe Cercopithecini) are the most widely distributed nonhuman primate in the tropical forest belt of Africa and show considerable phenotypic, taxonomic, and ecological diversity. However, genomic information for most species within this group is still lacking. Here, we present a high-quality de novo genome (total 2.90 Gb, contig N50 equal to 22.7 Mb) of the mona monkey (Cercopithecus mona), together with genome resequencing data of 13 individuals sampled across Nigeria. Our results showed differentiation between populations from East and West of the Niger River ∼84 ka and potential ancient introgression in the East population from other mona group species. The PTPRK, FRAS1, BNC2, and EDN3 genes related to pigmentation displayed signals of introgression in the East population. Genomic scans suggest that immunity genes such as AKT3 and IL13 (possibly involved in simian immunodeficiency virus defense), and G6PD, a gene involved in malaria resistance, are under positive natural selection. Our study gives insights into differentiation, natural selection, and introgression in guenons.

Constitutively-expressed and induced immune effectors in the house fly (Musca domestica) and the transcription factors that may regulate them.

Insects possess both infection-induced and constitutively expressed innate immune defences. Some effectors, such as lysozymes and antimicrobial peptides (AMPs), are constitutively expressed in flies, but expression patterns vary across tissues and species. The house fly (Musca domestica L.) has an impressive immune repertoire, with more effector genes than any other flies. We used RNA-seq to explore both constitutive and induced expression of immune effectors in flies. House flies were fed either Pseudomonas aeruginosa or Escherichia coli, or sterile control broth, and gene expression in the gut and carcass was analysed 4 h post-feeding. Flies fed either bacterium did not induce AMP expression, but some lysozyme and AMP genes were constitutively expressed. Prior transcriptome data from flies injected with bacteria also were analysed, and these constitutively expressed genes differed from those induced by bacterial injection. Binding sites for the transcription factor Myc were enriched upstream of constitutively expressed AMP genes, while upstream regions of induced AMPs were enriched for NF-κB binding sites resembling those of the Imd-responsive transcription factor Relish. Therefore, we identified at least two expression repertoires for AMPs in the house fly: constitutively expressed genes that may be regulated by Myc, and induced AMPs likely regulated by Relish.

Ecology and the evolution of sex chromosomes.

Sex chromosomes are common features of animal genomes, often carrying a sex determination gene responsible for initiating the development of sexually dimorphic traits. The specific chromosome that serves as the sex chromosome differs across taxa as a result of fusions between sex chromosomes and autosomes, along with sex chromosome turnover-autosomes becoming sex chromosomes and sex chromosomes 'reverting' back to autosomes. In addition, the types of genes on sex chromosomes frequently differ from the autosomes, and genes on sex chromosomes often evolve faster than autosomal genes. Sex-specific selection pressures, such as sexual antagonism and sexual selection, are hypothesized to be responsible for sex chromosome turnovers, the unique gene content of sex chromosomes and the accelerated evolutionary rates of genes on sex chromosomes. Sex-specific selection has pronounced effects on sex chromosomes because their sex-biased inheritance can tilt the balance of selection in favour of one sex. Despite the general consensus that sex-specific selection affects sex chromosome evolution, most population genetic models are agnostic as to the specific sources of these sex-specific selection pressures, and many of the details about the effects of sex-specific selection remain unresolved. Here, I review the evidence that ecological factors, including variable selection across heterogeneous environments and conflicts between sexual and natural selection, can be important determinants of sex-specific selection pressures that shape sex chromosome evolution. I also explain how studying the ecology of sex chromosome evolution can help us understand important and unresolved aspects of both sex chromosome evolution and sex-specific selection.

Evolution of Sex Determination and Sex Chromosomes: A Novel Alternative Paradigm.

Sex chromosomes can differ between species as a result of evolutionary turnover, a process that can be driven by evolution of the sex determination pathway. Canonical models of sex chromosome turnover hypothesize that a new master sex determining gene causes an autosome to become a sex chromosome or an XY chromosome pair to switch to a ZW pair (or vice versa). Here, a novel paradigm for the evolution of sex determination and sex chromosomes is presented, in which there is an evolutionary transition in the master sex determiner, but the X chromosome remains unchanged. There are three documented examples of the novel paradigm, and it is hypothesized that a similar process could happen in a ZW sex chromosome system. Three other taxa are also identified where the novel paradigm may have occurred, and how it could be distinguished from canonical trajectories in these and additional taxa is also described.

The genome of the stable fly, Stomoxys calcitrans, reveals potential mechanisms underlying reproduction, host interactions, and novel targets for pest control.

BACKGROUND: The stable fly, Stomoxys calcitrans, is a major blood-feeding pest of livestock that has near worldwide distribution, causing an annual cost of over $2 billion for control and product loss in the USA alone. Control of these flies has been limited to increased sanitary management practices and insecticide application for suppressing larval stages. Few genetic and molecular resources are available to help in developing novel methods for controlling stable flies. RESULTS: This study examines stable fly biology by utilizing a combination of high-quality genome sequencing and RNA-Seq analyses targeting multiple developmental stages and tissues. In conjunction, 1600 genes were manually curated to characterize genetic features related to stable fly reproduction, vector host interactions, host-microbe dynamics, and putative targets for control. Most notable was characterization of genes associated with reproduction and identification of expanded gene families with functional associations to vision, chemosensation, immunity, and metabolic detoxification pathways. CONCLUSIONS: The combined sequencing, assembly, and curation of the male stable fly genome followed by RNA-Seq and downstream analyses provide insights necessary to understand the biology of this important pest. These resources and new data will provide the groundwork for expanding the tools available to control stable fly infestations. The close relationship of Stomoxys to other blood-feeding (horn flies and Glossina) and non-blood-feeding flies (house flies, medflies, Drosophila) will facilitate understanding of the evolutionary processes associated with development of blood feeding among the Cyclorrhapha.

Temperature-dependent effects of house fly proto-Y chromosomes on gene expression could be responsible for fitness differences that maintain polygenic sex determination.

Sex determination, the developmental process by which sexually dimorphic phenotypes are established, evolves fast. Evolutionary turnover in a sex determination pathway may occur via selection on alleles that are genetically linked to a new master sex determining locus on a newly formed proto-sex chromosome. Species with polygenic sex determination, in which master regulatory genes are found on multiple different proto-sex chromosomes, are informative models to study the evolution of sex determination and sex chromosomes. House flies are such a model system, with male determining loci possible on all six chromosomes and a female-determiner on one of the chromosomes as well. The two most common male-determining proto-Y chromosomes form latitudinal clines on multiple continents, suggesting that temperature variation is an important selection pressure responsible for maintaining polygenic sex determination in this species. Temperature-dependent fitness effects could be manifested through temperature-dependent gene expression differences across proto-Y chromosome genotypes. These gene expression differences may be the result of cis regulatory variants that affect the expression of genes on the proto-sex chromosomes, or trans effects of the proto-Y chromosomes on genes elswhere in the genome. We used RNA-seq to identify genes whose expression depends on proto-Y chromosome genotype and temperature in adult male house flies. We found no evidence for ecologically meaningful temperature-dependent expression differences of sex determining genes between male genotypes, but we were probably not sampling an appropriate developmental time-point to identify such effects. In contrast, we identified many other genes whose expression depends on the interaction between proto-Y chromosome genotype and temperature, including genes that encode proteins involved in reproduction, metabolism, lifespan, stress response, and immunity. Notably, genes with genotype-by-temperature interactions on expression were not enriched on the proto-sex chromosomes. Moreover, there was no evidence that temperature-dependent expression is driven by chromosome-wide cis-regulatory divergence between the proto-Y and proto-X alleles. Therefore, if temperature-dependent gene expression is responsible for differences in phenotypes and fitness of proto-Y genotypes across house fly populations, these effects are driven by a small number of temperature-dependent alleles on the proto-Y chromosomes that may have trans effects on the expression of genes on other chromosomes.

The house fly Y Chromosome is young and minimally differentiated from its ancient X Chromosome partner.

Canonical ancient sex chromosome pairs consist of a gene rich X (or Z) Chromosome and a male-limited (or female-limited) Y (or W) Chromosome that is gene poor. In contrast to highly differentiated sex chromosomes, nascent sex chromosome pairs are homomorphic or very similar in sequence content. Nascent sex chromosomes can arise if an existing sex chromosome fuses to an autosome or an autosome acquires a new sex-determining locus/allele. Sex chromosomes often differ between closely related species and can even be polymorphic within species, suggesting that nascent sex chromosomes arise frequently over the course of evolution. Previously documented sex chromosome transitions involve changes to both members of the sex chromosome pair (X and Y, or Z and W). The house fly has sex chromosomes that resemble the ancestral fly karyotype that originated ∼100 million yr ago; therefore, the house fly is expected to have X and Y Chromosomes with different gene content. We tested this hypothesis using whole-genome sequencing and transcriptomic data, and we discovered little evidence for genetic differentiation between the X and Y in house fly. We propose that the house fly has retained the ancient X Chromosome, but the ancestral Y was replaced by an X Chromosome carrying a new male determining gene. Our proposed hypothesis provides a mechanism for how one member of a sex chromosome pair can experience evolutionary turnover while the other member remains unaffected.

The X chromosome of the German cockroach, Blattella germanica, is homologous to a fly X chromosome despite 400 million years divergence.

BACKGROUND: Sex chromosome evolution is a dynamic process that can proceed at varying rates across lineages. For example, different chromosomes can be sex-linked between closely related species, whereas other sex chromosomes have been conserved for > 100 million years. Cases of long-term sex chromosome conservation could be informative of factors that constrain sex chromosome evolution. Cytological similarities between the X chromosomes of the German cockroach (Blattella germanica) and most flies suggest that they may be homologous-possibly representing an extreme case of long-term conservation. RESULTS: To test the hypothesis that the cockroach and fly X chromosomes are homologous, we analyzed whole-genome sequence data from cockroaches. We found evidence in both sequencing coverage and heterozygosity that a significant excess of the same genes are on both the cockroach and fly X chromosomes. We also present evidence that the candidate X-linked cockroach genes may be dosage compensated in hemizygous males. Consistent with this hypothesis, three regulators of transcription and chromatin on the fly X chromosome are conserved in the cockroach genome. CONCLUSIONS: Our results support our hypothesis that the German cockroach shares the same X chromosome as most flies. This may represent the convergent evolution of the X chromosome in the lineages leading to cockroaches and flies. Alternatively, the common ancestor of most insects may have had an X chromosome that resembled the extant cockroach and fly X. Cockroaches and flies diverged ∼ 400 million years ago, which would be the longest documented conservation of a sex chromosome. Cockroaches and flies have different mechanisms of sex determination, raising the possibility that the X chromosome was conserved despite the evolution of the sex determination pathway.

Genes Relocated Between Drosophila Chromosome Arms Evolve Under Relaxed Selective Constraints Relative to Non-Relocated Genes.

Gene duplication creates a second copy of a gene either in tandem to the ancestral locus or dispersed to another chromosomal location. When the ancestral copy of a dispersed duplicate is lost from the genome, it creates the appearance that the gene was "relocated" from the ancestral locus to the derived location. Gene relocations may be as common as canonical dispersed duplications in which both the ancestral and derived copies are retained. Relocated genes appear to be under more selective constraints than the derived copies of canonical duplications, and they are possibly as conserved as single-copy non-relocated genes. To test this hypothesis, we combined comparative genomics, population genetics, gene expression, and functional analyses to assess the selection pressures acting on relocated, duplicated, and non-relocated single-copy genes in Drosophila genomes. We find that relocated genes evolve faster than single-copy non-relocated genes, and there is no evidence that this faster evolution is driven by positive selection. In addition, relocated genes are less essential for viability and male fertility than single-copy non-relocated genes, suggesting that relocated genes evolve fast because of relaxed selective constraints. However, relocated genes evolve slower than the derived copies of canonical dispersed duplicated genes. We therefore conclude that relocated genes are under more selective constraints than canonical duplicates, but are not as conserved as single-copy non-relocated genes.

Sexual conflict and the maintenance of genetic variation in natural populations.

Understanding the factors that maintain genetic variation in natural populations is a foundational goal of evolutionary biology. To this end, population geneticists have developed a variety of models that can produce stable polymorphisms. In one of the earliest models, Owen () demonstrated that differences in selection pressures acting on males and females could maintain multiple alleles of a gene at a stable equilibrium. If the selection pressures act in opposite directions in males and females, we refer to this as (inter-) sexual conflict or sexual antagonism (Arnqvist & Rowe, ). Testing if sexual conflict maintains genetic variation in natural populations is a tremendous challenge-it requires both identifying loci that harbor sexually antagonistic alleles and determining whether those alleles are maintained as stable polymorphisms (Mank, ). Doing so genome-wide is even harder because it is not tractable to identify sexually antagonistic alleles and test for stable polymorphisms at all loci. Dutoit et al. () confront this challenge in a paper published in this issue of Molecular Ecology. Using gene expression and population genomic data from the collared flycatcher, Dutoit et al. () identify associations and correlations between genomic signatures of balanced polymorphisms and sexual conflict.

Using genomic data to study insecticide resistance in the house fly, Musca domestica.

The house fly, Musca domestica, is a major pest at livestock facilities throughout the world. Insecticides have been the most common control strategy for flies, but many populations have evolved resistance. The speed by which we are able to identify the mutations responsible for resistance has been a major challenge for the development of high throughput resistance monitoring assays as new insecticides are introduced for control. This is particularly true for mutations that cause trans regulation of a gene, which then results in resistance. In this paper we take advantage of the conserved homology of dipteran chromosomes to assign 3069 genes to chromosomes. Of these, 234 were of toxicological interest (CYPs, esterases/hydrolases, glutathione S-transferases (GSTs) and target sites). The chromosomal location of genes known from linkage analysis studies matched the location predicted by homology mapping in ten out of ten cases, indicating a high reliability of our approach. The CYPs, esterases/hydrolases and GSTs were not randomly distributed throughout the genome. They clustered on chromosomes, but the pattern was different between the CYPs, esterases/hydrolases and GSTs. Examples are provided for how the availability of the house fly genome, combined with an ability to assign genes to chromosomes, will help to accelerate research in house flies.

The faster-X effect: integrating theory and data.

Population genetics theory predicts that X (or Z) chromosomes could play disproportionate roles in speciation and evolutionary divergence, and recent genome-wide analyses have identified situations in which X or Z-linked divergence exceeds that on the autosomes (the so-called 'faster-X effect'). Here, we summarize the current state of both the theory and data surrounding the study of faster-X evolution. Our survey indicates that the faster-X effect is pervasive across a taxonomically diverse array of evolutionary lineages. These patterns could be informative of the dominance or recessivity of beneficial mutations and the nature of genetic variation acted upon by natural selection. We also identify several aspects of disagreement between these empirical results and the population genetic models used to interpret them. However, there are clearly delineated aspects of the problem for which additional modeling and collection of genomic data will address these discrepancies and provide novel insights into the population genetics of adaptation.

Is Multifactorial Sex Determination in the House Fly, Musca domestica (L.), Stable Over Time?

Sex determination pathways evolve rapidly, usually because of turnover of master regulatory genes at the top of the developmental pathway. Polygenic sex determination is expected to be a transient state between ancestral and derived conditions. However, polygenic sex determination has been observed in numerous animal species, including the house fly, Musca domestica House fly males carry a male-determining factor (M) that can be located on any chromosome, and an individual male may have multiple M factors. Females lack M and/or have a dominant allele of the Md-tra gene (Md-tra D ) that acts as a female-determining locus even in the presence of multiple copies of M. We found the frequency and linkage of M in house flies collected in Chino, CA (USA) was relatively unchanged between 1982 and 2014. The frequency of females with Md-tra D in the 2014 collection was 33.6% (n = 140). Analysis of these results, plus previously published data, revealed a strong correlation between the frequencies of Md-tra D and multiple M males, and we find that these populations are expected to have balanced sex ratios. We also find that fitness values that allow for the invasion and maintenance of multiple sex determining loci suggest that sexually antagonistic selection could be responsible for maintaining polygenic sex determination in house fly populations. The stability over time and equilibrium frequencies within populations suggest the house fly polygenic sex determination system is not in transition, and provide guidance for future investigations on the factors responsible for the polymorphism.

Disentangling the relationship between sex-biased gene expression and X-linkage.

X chromosomes are preferentially transmitted through females, which may favor the accumulation of X-linked alleles/genes with female-beneficial effects. Numerous studies have shown that genes with sex-biased expression are under- or over-represented on the X chromosomes of a wide variety of organisms. The patterns, however, vary between different animal species, and the causes of these differences are unresolved. Additionally, genes with sex-biased expression tend to be narrowly expressed in a limited number of tissues, and narrowly expressed genes are also non-randomly X-linked in a taxon-specific manner. It is therefore unclear whether the unique gene content of the X chromosome is the result of selection on genes with sex-biased expression, narrowly expressed genes, or some combination of the two. To address this problem, we measured sex-biased expression in multiple Drosophila species and at different developmental time points. These data were combined with available expression measurements from Drosophila melanogaster and mouse to reconcile the inconsistencies in X-chromosome content among taxa. Our results suggest that most of the differences between Drosophila and mammals are confounded by disparate data collection/analysis approaches as well as the correlation between sex bias and expression breadth. Both the Drosophila and mouse X chromosomes harbor an excess of genes with female-biased expression after controlling for the confounding factors, suggesting that the asymmetrical transmission of the X chromosome favors the accumulation of female-beneficial mutations in X-linked genes. However, some taxon-specific patterns remain, and we provide evidence that these are in part a consequence of constraints imposed by the dosage compensation mechanism in Drosophila.

Faster-X evolution of gene expression in Drosophila.

DNA sequences on X chromosomes often have a faster rate of evolution when compared to similar loci on the autosomes, and well articulated models provide reasons why the X-linked mode of inheritance may be responsible for the faster evolution of X-linked genes. We analyzed microarray and RNA-seq data collected from females and males of six Drosophila species and found that the expression levels of X-linked genes also diverge faster than autosomal gene expression, similar to the "faster-X" effect often observed in DNA sequence evolution. Faster-X evolution of gene expression was recently described in mammals, but it was limited to the evolutionary lineages shortly following the creation of the therian X chromosome. In contrast, we detect a faster-X effect along both deep lineages and those on the tips of the Drosophila phylogeny. In Drosophila males, the dosage compensation complex (DCC) binds the X chromosome, creating a unique chromatin environment that promotes the hyper-expression of X-linked genes. We find that DCC binding, chromatin environment, and breadth of expression are all predictive of the rate of gene expression evolution. In addition, estimates of the intraspecific genetic polymorphism underlying gene expression variation suggest that X-linked expression levels are not under relaxed selective constraints. We therefore hypothesize that the faster-X evolution of gene expression is the result of the adaptive fixation of beneficial mutations at X-linked loci that change expression level in cis. This adaptive faster-X evolution of gene expression is limited to genes that are narrowly expressed in a single tissue, suggesting that relaxed pleiotropic constraints permit a faster response to selection. Finally, we present a conceptional framework to explain faster-X expression evolution, and we use this framework to examine differences in the faster-X effect between Drosophila and mammals.

The poly(A) polymerase GLD2 is required for spermatogenesis in Drosophila melanogaster.

The DNA of a developing sperm is normally inaccessible for transcription for part of spermatogenesis in many animals. In Drosophila melanogaster, many transcripts needed for late spermatid differentiation are synthesized in pre-meiotic spermatocytes, but are not translated until later stages. Thus, post-transcriptional control mechanisms are required to decouple transcription and translation during spermatogenesis. In the female germline, developing germ cells accomplish similar decoupling through poly(A) tail alterations to ensure that dormant transcripts are not prematurely translated: a transcript with a short poly(A) tail will remain untranslated, whereas elongating the poly(A) tail permits protein production. In Drosophila, the ovary-expressed cytoplasmic poly(A) polymerase WISPY is responsible for stage-specific poly(A) tail extension in the female germline. Here, we examine the possibility that a recently derived testis-expressed WISPY paralog, GLD2, plays a similar role in the Drosophila male germline. We show that knockdown of Gld2 transcripts causes male sterility, as GLD2-deficient males do not produce mature sperm. Spermatogenesis up to and including meiosis appears normal in the absence of GLD2, but post-meiotic spermatid development rapidly becomes abnormal. Nuclear bundling and F-actin assembly are defective in GLD2 knockdown testes and nuclei fail to undergo chromatin reorganization in elongated spermatids. GLD2 also affects the incorporation of protamines and the stability of dynamin and transition protein transcripts. Our results indicate that GLD2 is an important regulator of late spermatogenesis and is the first example of a Gld-2 family member that plays a significant role specifically in male gametogenesis.