Genetic variation affecting host-parasite interactions: different genes affect different aspects of sigma virus replication and transmission in Drosophila melanogaster

In natural populations, genetic variation affects resistance to disease. Knowing how much variation exists, and understanding the genetic architecture of this variation, is important for medicine, for agriculture, and for understanding evolutionary processes. To investigate the extent and nature of genetic variation affecting resistance to pathogens, we are studying a tractable model system: Drosophila melanogaster and its natural pathogen the vertically transmitted sigma virus. We show that considerable genetic variation affects transmission of the virus from parent to offspring. However, maternal and paternal transmission of the virus is affected by different genes. Maternal transmission is a simple Mendelian trait: most of the genetic variation is explained by a polymorphism in ref(2)P, a gene already well known to affect resistance to sigma. In contrast, there is considerable genetic variation in paternal transmission that cannot be explained by ref(2)P and is caused by other loci on chromosome 2. Furthermore, we found no genetic correlation between paternal transmission of the virus and resistance to infection by the sigma virus following injection. This suggests that different loci affect viral replication and paternal transmission.     

The genetic architecture of resistance to virus infection in Drosophila

Variation in susceptibility to infection has a substantial genetic component in natural populations, and it has been argued that selection by pathogens may result in it having a simpler genetic architecture than many other quantitative traits. This is important as models of host–pathogen co‐evolution typically assume resistance is controlled by a small number of genes. Using the Drosophila melanogaster multiparent advanced intercross, we investigated the genetic architecture of resistance to two naturally occurring viruses, the sigma virus and DCV (Drosophila C virus). We found extensive genetic variation in resistance to both viruses. For DCV resistance, this variation is largely caused by two major‐effect loci. Sigma virus resistance involves more genes – we mapped five loci, and together these explained less than half the genetic variance. Nonetheless, several of these had a large effect on resistance. Models of co‐evolution typically assume strong epistatic interactions between polymorphisms controlling resistance, but we were only able to detect one locus that altered the effect of the main effect loci we had mapped. Most of the loci we mapped were probably at an intermediate frequency in natural populations. Overall, our results are consistent with major‐effect genes commonly affecting susceptibility to infectious diseases, with DCV resistance being a near‐Mendelian trait.     

The recent spread of a vertically transmitted virus through populations of Drosophila melanogaster

The sigma virus is a vertically transmitted pathogen that commonly infects natural populations of Drosophila melanogaster. This virus is the only known host-specific pathogen of D. melanogaster, and so offers a unique opportunity to study the genetics of Drosophila-viral interactions in a natural system. To elucidate the population genetic processes that operate in sigma virus populations, we collected D. melanogaster from 10 populations across three continents. We found that the sigma virus had a prevalence of 0-15% in these populations. Compared to other RNA viruses, we found that levels of viral genetic diversity are very low across Europe and North America. Based on laboratory measurements of the viral substitution rate, we estimate that most European and North American viral isolates shared a common ancestor approximately 200 years ago. We suggest two explanations for this: the first is that D. melanogaster has recently acquired the sigma virus; the second is that a single viral type has recently swept through D. melanogaster populations. Furthermore, in contrast to Drosophila populations, we find that the sigma viral populations are highly structured. This is surprising for a vertically transmitted pathogen that has a similar migration rate to its host. We suggest that the low structure in the viral populations can be explained by the smaller effective population size of the virus.     

Phenotypic variation and natural selection at catsup, a pleiotropic quantitative trait gene in Drosophila

Quantitative traits are shaped by networks of pleiotropic genes . To understand the mechanisms that maintain genetic variation for quantitative traits in natural populations and to predict responses to artificial and natural selection, we must evaluate pleiotropic effects of underlying quantitative trait genes and define functional allelic variation at the level of quantitative trait nucleotides (QTNs). Catecholamines up (Catsup), which encodes a negative regulator of tyrosine hydroxylase , the rate-limiting step in the synthesis of the neurotransmitter dopamine, is a pleiotropic quantitative trait gene in Drosophila melanogaster. We used association mapping to determine whether the same or different QTNs at Catsup are associated with naturally occurring variation in multiple quantitative traits. We sequenced 169 Catsup alleles from a single population and detected 33 polymorphisms with little linkage disequilibrium (LD). Different molecular polymorphisms in Catsup are independently associated with variation in longevity, locomotor behavior, and sensory bristle number. Most of these polymorphisms are potentially functional variants in protein coding regions, have large effects, and are not common. Thus, Catsup is a pleiotropic quantitative trait gene, but individual QTNs do not have pleiotropic effects. Molecular population genetic analyses of Catsup sequences are consistent with balancing selection maintaining multiple functional polymorphisms.     

Host–parasite coevolution: genetic variation in a virus population and the interaction with a host gene

Host–parasite coevolution is considered to be an important factor in maintaining genetic variation in resistance to pathogens. Drosophila melanogaster is naturally infected by the sigma virus, a vertically transmitted and host‐specific pathogen. In fly populations, there is a large amount of genetic variation in the transmission rate from parent to offspring, much of which is caused by major‐effect resistance polymorphisms. We have found that there are similarly high levels of genetic variation in the rate of paternal transmission among 95 different isolates of the virus as in the host. However, when we examined a transmission‐blocking gene in the host, we found that it was effective across virus isolates. Therefore, the high levels of genetic variation observed in this system do not appear to be maintained because of coevolution resulting from interactions between this host gene and parasite genes.     

Population-based resequencing of experimentally evolved populations reveals the genetic basis of body size variation in Drosophila melanogaster

Body size is a classic quantitative trait with evolutionarily significant variation within many species. Locating the alleles responsible for this variation would help understand the maintenance of variation in body size in particular, as well as quantitative traits in general. However, successful genome-wide association of genotype and phenotype may require very large sample sizes if alleles have low population frequencies or modest effects. As a complementary approach, we propose that population-based resequencing of experimentally evolved populations allows for considerable power to map functional variation. Here, we use this technique to investigate the genetic basis of natural variation in body size in Drosophila melanogaster. Significant differentiation of hundreds of loci in replicate selection populations supports the hypothesis that the genetic basis of body size variation is very polygenic in D. melanogaster. Significantly differentiated variants are limited to single genes at some loci, allowing precise hypotheses to be formed regarding causal polymorphisms, while other significant regions are large and contain many genes. By using significantly associated polymorphisms as a priori candidates in follow-up studies, these data are expected to provide considerable power to determine the genetic basis of natural variation in body size.     

Genetic variation for dorsal-ventral patterning of the Drosophila melanogaster eggshell

Patterning of the insect eggshell is an excellent system for exploring the molecular basis of phenotypic variation. In Drosophila melanogaster, two dorsal-anterior respiratory appendages are produced in response to signaling through the Epidermal growth factor receptor (Egfr). Previous work implicates Egfr pathway function in both intraspecific variation for dorsal appendage spacing (DAS) on the eggshell, as well as interspecific differences in dorsal appendage number and location. To test the hypothesis that genetic variation in Egfr contributes to variation in eggshell patterning, we have made use of naturally occurring intraspecific variation for DAS as a model quantitative trait. We found that there is substantial segregating genetic variation for DAS in D. melanogaster, and have tested for associations with 289 common polymorphisms in the Egfr locus. A marginal association was seen with two polymorphic sites in Egfr; however, we failed to replicate these findings in a second population, or in a modified quantitative complementation test designed to specifically test the effects of the putative polymorphisms. Therefore, we conclude that the polymorphisms we have identified in Egfr do not contribute to variation in DAS, and further work is required to understand the genetic architecture of this trait.     

Rhabdoviruses in two species of Drosophila: vertical transmission and a recent sweep

Insects are host to a diverse range of vertically transmitted micro-organisms, but while their bacterial symbionts are well-studied, little is known about their vertically transmitted viruses. We have found that two sigma viruses (Rhabdoviridae) recently discovered in Drosophila affinis and Drosophila obscura are both vertically transmitted. As is the case for the sigma virus of Drosophila melanogaster, we find that both males and females can transmit these viruses to their offspring. Males transmit lower viral titers through sperm than females transmit through eggs, and a lower proportion of their offspring become infected. In natural populations of D. obscura in the United Kingdom, we found that 39% of flies were infected and that the viral population shows clear evidence of a recent expansion, with extremely low genetic diversity and a large excess of rare polymorphisms. Using sequence data we estimate that the virus has swept across the United Kingdom within the past ～11 years, during which time the viral population size doubled approximately every 9 months. Using simulations based on our lab estimates of transmission rates, we show that the biparental mode of transmission allows the virus to invade and rapidly spread through populations at rates consistent with those measured in the field. Therefore, as predicted by our simulations, the virus has undergone an extremely rapid and recent increase in population size. In light of this and earlier studies of a related virus in D. melanogaster, we conclude that vertically transmitted rhabdoviruses may be common in insects and that these host-parasite interactions can be highly dynamic.     

Quantitative trait loci for locomotor behavior in Drosophila melanogaster

Locomotion is an integral component of most animal behaviors and many human diseases and disorders are associated with locomotor deficits, but little is known about the genetic basis of natural variation in locomotor behavior. Locomotion is a complex trait, with variation attributable to the joint segregation of multiple interacting quantitative trait loci (QTL), with effects that are sensitive to the environment. We assessed variation in a component of locomotor behavior (locomotor reactivity) in a population of 98 recombinant inbred lines of Drosophila melanogaster and mapped four QTL affecting locomotor reactivity by linkage to polymorphic roo transposable element insertion sites. We used complementation tests of deficiencies to fine map these QTL to 12 chromosomal regions and complementation tests of mutations to identify 13 positional candidate genes affecting locomotor reactivity, including Dopa decarboxylase (Ddc), which catalyzes the final step in the synthesis of serotonin and dopamine. Linkage disequilibrium mapping in a population of 164 second chromosome substitution lines derived from a single natural population showed that polymorphisms at Ddc were associated with naturally occurring genetic variation in locomotor behavior. These data implicate variation in the synthesis of bioamines as a factor contributing to natural variation in locomotor reactivity.     

Genetic variation in rates of nondisjunction: association of two naturally occurring polymorphisms in the chromokinesin nod with increased rates of nondisjunction in Drosophila melanogaster.

Genetic variation in nondisjunction frequency among X chromosomes from two Drosophila melanogaster natural populations is examined in a sensitized assay. A high level of genetic variation is observed (a range of 0.006-0.241). Two naturally occurring variants at the nod locus, a chromokinesin required for proper achiasmate chromosome segregation, are significantly associated with an increased frequency of nondisjunction. Both of these polymorphisms are found at intermediate frequency in widely distributed natural populations. To account for these observations, we propose a general model incorporating unique opportunities for meiotic drive during female meiosis. The oötid competition model can account for both high mean rates of female-specific nondisjunction in Drosophila and humans as well as the standing genetic variation in this critical fitness character in natural populations.     

QTL analysis of floral traits in Louisiana iris hybrids

The formation of hybrid zones between nascent species is a widespread phenomenon. The evolutionary consequences of hybridization are influenced by numerous factors, including the action of natural selection on quantitative trait variation. Here we examine how the genetic basis of floral traits of two species of Louisiana Irises affects the extent of quantitative trait variation in their hybrids. Quantitative trait locus (QTL) mapping was used to assess the size (magnitude) of phenotypic effects of individual QTL, the degree to which QTL for different floral traits are colocalized, and the occurrence of mixed QTL effects. These aspects of quantitative genetic variation would be expected to influence (1) the number of genetic steps (in terms of QTL substitutions) separating the parental species phenotypes; (2) trait correlations; and (3) the potential for transgressive segregation in hybrid populations. Results indicate that some Louisiana Iris floral trait QTL have large effects and QTL for different traits tend to colocalize. Transgressive variation was observed for six of nine traits, despite the fact that mixed QTL effects influence few traits. Overall, our QTL results imply that the genetic basis of floral morphology and color traits might facilitate the maintenance of phenotypic divergence between Iris fulva and Iris brevicaulis, although a great deal of phenotypic variation was observed among hybrids.     

Quantitative genetic analysis of natural variation in body size in Drosophila melanogaster

Latitudinal, genetic variation in body size is a commonly observed phenomenon in many invertebrate species and is shaped by natural selection. In this study, we use a chromosome substitution and a quantitative trait locus (QTL) mapping approach to identify chromosomes and genomic regions associated with adaptive variation in body size in natural populations of Drosophila melanogaster from the extreme ends of clines in South America and Australia. Chromosome substitution revealed the largest effects on chromosome three in both continents, and minor effects on the X and second chromosome. Similarly, QTL analysis of the Australian cline identified QTL with largest effects on the third chromosome, with smaller effects on the second. However, no QTL were found on the X chromosome. We also compared the coincidence of locations of QTL with the locations of five microsatellite loci previously shown to vary clinally in Australia. Permutation tests using both the sum of the LOD scores and the sum distance to nearest QTL peak revealed there were no significant associations between locations of clinal markers and QTL's. The lack of significance may, in part, be due to broad QTL peaks identified in this study. Future studies using higher resolution QTL maps should reveal whether the degree of clinality in microsatellite allele frequencies can be used to identify QTL in traits that vary along an environmental gradient.     

Genetic Load in Natural Populations: Is It Compatible with the Hypothesis That Many Polymorphisms Are Maintained by Natural Selection?

Recent studies of genetically controlled enzyme variation lead to an estimation that at least 30 to 60% of the structural genes are polymorphic in natural populations of many vertebrate and invertebrate species. Some authors have argued that a substantial proportion of these polymorphisms cannot be maintained by natural selection because this would result in an unbearable genetic load. If many polymorphisms are maintained by heterotic natural selection, individuals with much greater than average proportion of homozygous loci should have very low fitness. We have measured in Drosophila melanogaster the fitness of flies homozygous for a complete chromosome relative to normal wild flies. A total of 37 chromosomes from a natural population have been tested using 92 experimental populations. The mean fitness of homozygous flies is 0.12 for second chromosomes, and 0.13 for third chromosomes. These estimates are compatible with the hypothesis that many (more than one thousand) loci are maintained by heterotic selection in natural populations of D. melanogaster.     

Molecular Dissection of a Major Gene Effect on a Quantitative Trait: The Level of Alcohol Dehydrogenase Expression in Drosophila Melanogaster

A molecular mapping experiment shows that a major gene effect on a quantitative trait, the level of alcohol dehydrogenase expression in Drosophila melanogaster, is due to multiple polymorphisms within the Adh gene. These polymorphisms are located in an intron, the coding sequence, and the 3' untranslated region. Because of nonrandom associations among polymorphisms at different sites, the individual effects combine (in some cases epistatically) to produce ``superalleles' with large effect. These results have implications for the interpretation of major gene effects detected by quantitative trait locus mapping methods. They show that large effects due to a single locus may be due to multiple associated polymorphisms (or sequential fixations in isolated populations) rather than individual mutations of large effect.     

Climatic selection on genes and traits after a 100 year-old invasion: a critical look at the temperate-tropical clines in Drosophila melanogaster from eastern Australia

Drosophila melanogaster invaded Australia around 100 years ago, most likely through a northern invasion. The wide range of climatic conditions in eastern Australia across which D. melanogaster is now found provides an opportunity for researchers to identify traits and genes that are associated with climatic adaptation. Allozyme studies indicate clinal patterns for at least four loci including a strong linear cline in Adh and a non-linear cline in alpha-Gpdh. Inversion clines were initially established from cytological studies but have now been validated with larger sample sizes using molecular markers for breakpoints. Recent collections indicate that some genetic markers (Adh and In(3R)Payne) have changed over the last 20 years reflecting continuing evolution. Heritable clines have been established for quantitative traits including wing length/area, thorax length and cold and heat resistance. A cline in egg size independent of body size and a weak cline in competitive ability have also been established. Postulated clinal patterns for resistance to desiccation and starvation have not been supported by extensive sampling. Experiments under laboratory and semi-natural conditions have suggested selective factors generating clinal patterns, particularly for reproductive patterns over winter. Attempts are being made to link clinal variation in traits to specific genes using QTL analysis and the candidate locus approach, and to identify the genetic architecture of trait variation along the cline. This is proving difficult because of inversion polymorphisms that generate disequilibrium among genes. Substantial gaps still remain in linking clines to field selection and understanding the genetic and physiological basis of the adaptive shifts. However D. melanogaster populations in eastern Australia remain an excellent resource for understanding past and future evolutionary responses to climate change.     

The quantitative genetic basis of sex ratio variation in Nasonia vitripennis: a QTL study

Our understanding of how natural selection should shape sex allocation is perhaps more developed than for any other trait. However, this understanding is not matched by our knowledge of the genetic basis of sex allocation. Here, we examine the genetic basis of sex ratio variation in the parasitoid wasp Nasonia vitripennis, a species well known for its response to local mate competition (LMC). We identified a quantitative trait locus (QTL) for sex ratio on chromosome 2 and three weaker QTL on chromosomes 3 and 5. We tested predictions that genes associated with sex ratio should be pleiotropic for other traits by seeing if sex ratio QTL co-occurred with clutch size QTL. We found one clutch size QTL on chromosome 1, and six weaker QTL across chromosomes 2, 3 and 5, with some overlap to regions associated with sex ratio. The results suggest rather limited scope for pleiotropy between these traits.     

Detection of a quantitative trait locus associated with resistance to Ascaris suum infection in pigs

Helminths almost invariably have an over-dispersed distribution in the host population. Human and animal studies have provided evidence suggesting that a large part of this variation is due to host genetic factors. Recently, the heritability for roundworm (Ascaris suum) infection levels in pigs was estimated to be 0.45. We used single nucleotide polymorphism markers to perform a whole-genome scan on 195 pigs experimentally infected with A. suum. A putative quantitative trait locus for worm burden on chromosome 4 covering 2.5 Mbp was identified by measured genotype analysis, although none of the SNPs reached genome-wide significance. To validate the putative quantitative trait locus, we genotyped two of the SNPs within the region in unrelated, informative animals exposed to experimental or natural infections and from which we had worm counts and/or faecal egg counts; the validation studies showed that one of the SNPs (TXNIP) was associated with total worm burden (P < 0.001) and adult worm burden(P < 0.0001), whereas the other SNP (ARNT) was associated with adult worm burden (P < 0.025) in these populations. We were thus able to confirm the existence of the quantitative trait locus on chromosome 4.This is to our knowledge the first report of a quantitative trait locus associated with helminth burden in pigs.     

Genetic basis of cross-resistance to three organophosphate insecticides in Drosophila melanogaster (Diptera: Drosophilidae)

To investigate the genetic basis of cross-resistance to insecticides, we conducted genetic analyses of resistance to three organophosphate insecticides, malathion, prothiophos, and fenitrothion. After isofemale lines resistant and susceptible to all of the three organophosphates had been screened from natural populations of Drosophila melanogaster (Meigen), chromosomal analyses were performed by using chromosome-substituted lines between the resistant and the susceptible lines. The chromosomal analyses revealed that both the second and the third chromosomes contributed to resistance to the organophosphates, suggesting that this resistant line possessed at least two factors for organophosphate resistance. However, the relative contribution of each chromosome was different in resistance to different organophosphates. We further carried out genetic mapping of a resistance factor for each organophosphate on each of the two chromosomes. Each resistance factor was mapped to the position of each chromosome, about II-62 and III-50. Results of the chromosomal analyses and the genetic mapping revealed that at least two resistance factors exhibiting different patterns of cross-resistance to the organophosphates existed within a natural population of D. melanogaster. Based on this research, genetic variation in insecticide resistance within natural populations and complex as well as simple aspects of the mechanism of cross-resistance are discussed.     

Successive increases in the resistance of Drosophila to viral infection through a transposon insertion followed by a Duplication

To understand the molecular basis of how hosts evolve resistance to their parasites, we have investigated the genes that cause variation in the susceptibility of Drosophila melanogaster to viral infection. Using a host-specific pathogen of D. melanogaster called the sigma virus (Rhabdoviridae), we mapped a major-effect polymorphism to a region containing two paralogous genes called CHKov1 and CHKov2. In a panel of inbred fly lines, we found that a transposable element insertion in the protein coding sequence of CHKov1 is associated with increased resistance to infection. Previous research has shown that this insertion results in a truncated messenger RNA that encodes a far shorter protein than the susceptible allele. This resistant allele has rapidly increased in frequency under directional selection and is now the commonest form of the gene in natural populations. Using genetic mapping and site-specific recombination, we identified a third genotype with considerably greater resistance that is currently rare in the wild. In these flies there have been two duplications, resulting in three copies of both the truncated allele of CHKov1 and CHKov2 (one of which is also truncated). Remarkably, the truncated allele of CHKov1 has previously been found to confer resistance to organophosphate insecticides. As estimates of the age of this allele predate the use of insecticides, it is likely that this allele initially functioned as a defence against viruses and fortuitously "pre-adapted" flies to insecticides. These results demonstrate that strong selection by parasites for increased host resistance can result in major genetic changes and rapid shifts in allele frequencies; and, contrary to the prevailing view that resistance to pathogens can be a costly trait to evolve, the pleiotropic effects of these changes can have unexpected benefits.     

Variation in genetic architecture of olfactory behaviour among wild-derived populations of Drosophila melanogaster

Odour-guided behaviour is a quantitative trait determined by many genes that are sensitive to gene-environment interactions. Different natural populations are likely to experience different selection pressures on the genetic underpinnings of chemosensory behaviour. However, few studies have reported comparisons of the quantitative genetic basis of olfactory behaviour in geographically distinct populations. We generated isofemale lines of Drosophila melanogaster from six populations in Argentina and measured larval and adult responses to benzaldehyde. There was significant variation within populations for both larval and adult olfactory behaviour and a significant genotype x sex interaction (GSI) for adult olfactory behaviour. However, there is substantial variation in the contribution of GSI to the total phenotypic variance among populations. Estimates of evolvability are orders of magnitude higher for larvae than for adults. Our results suggest that the potential for evolutionary adaptation to the chemosensory environment is greater at the larval feeding stage than at the adult reproductive stage.