Richard F. Lyman

The Drosophila melanogaster Genetic Reference Panel

A major challenge of biology is understanding the relationship between molecular genetic variation and variation in quantitative traits, including fitness. This relationship determines our ability to predict phenotypes from genotypes and to understand how evolutionary forces shape variation within and between species. Previous efforts to dissect the genotype–phenotype map were based on incomplete genotypic information. Here, we describe the Drosophila melanogaster Genetic Reference Panel (DGRP), a community resource for analysis of population genomics and quantitative traits. The DGRP consists of fully sequenced inbred lines derived from a natural population. Population genomic analyses reveal reduced polymorphism in centromeric autosomal regions and the X chromosome, evidence for positive and negative selection, and rapid evolution of the X chromosome. Many variants in novel genes, most at low frequency, are associated with quantitative traits and explain a large fraction of the phenotypic variance. The DGRP facilitates genotype–phenotype mapping using the power of Drosophila genetics.

Systems genetics of complex traits in Drosophila melanogaster.

Determining the genetic architecture of complex traits is challenging because phenotypic variation arises from interactions between multiple, environmentally sensitive alleles. We quantified genome-wide transcript abundance and phenotypes for six ecologically relevant traits in D. melanogaster wild-derived inbred lines. We observed 10,096 genetically variable transcripts and high heritabilities for all organismal phenotypes. The transcriptome is highly genetically intercorrelated, forming 241 transcriptional modules. Modules are enriched for transcripts in common pathways, gene ontology categories, tissue-specific expression and transcription factor binding sites. The high degree of transcriptional connectivity allows us to infer genetic networks and the function of predicted genes from annotations of other genes in the network. Regressions of organismal phenotypes on transcript abundance implicate several hundred candidate genes that form modules of biologically meaningful correlated transcripts affecting each phenotype. Overlapping transcripts in modules associated with different traits provide insight into the molecular basis of pleiotropy between complex traits.

Natural variation in genome architecture among 205 Drosophila melanogaster Genetic Reference Panel lines

The Drosophila melanogaster Genetic Reference Panel (DGRP) is a community resource of 205 sequenced inbred lines, derived to improve our understanding of the effects of naturally occurring genetic variation on molecular and organismal phenotypes. We used an integrated genotyping strategy to identify 4,853,802 single nucleotide polymorphisms (SNPs) and 1,296,080 non-SNP variants. Our molecular population genomic analyses show higher deletion than insertion mutation rates and stronger purifying selection on deletions. Weaker selection on insertions than deletions is consistent with our observed distribution of genome size determined by flow cytometry, which is skewed toward larger genomes. Insertion/deletion and single nucleotide polymorphisms are positively correlated with each other and with local recombination, suggesting that their nonrandom distributions are due to hitchhiking and background selection. Our cytogenetic analysis identified 16 polymorphic inversions in the DGRP. Common inverted and standard karyotypes are genetically divergent and account for most of the variation in relatedness among the DGRP lines. Intriguingly, variation in genome size and many quantitative traits are significantly associated with inversions. Approximately 50% of the DGRP lines are infected with Wolbachia, and four lines have germline insertions of Wolbachia sequences, but effects of Wolbachia infection on quantitative traits are rarely significant. The DGRP complements ongoing efforts to functionally annotate the Drosophila genome. Indeed, 15% of all D. melanogaster genes segregate for potentially damaged proteins in the DGRP, and genome-wide analyses of quantitative traits identify novel candidate genes. The DGRP lines, sequence data, genotypes, quality scores, phenotypes, and analysis and visualization tools are publicly available.

Dopa decarboxylase (Ddc) affects variation in Drosophila longevity

Mutational analyses in model organisms have shown that genes affecting metabolism and stress resistance regulate life span, but the genes responsible for variation in longevity in natural populations are largely unidentified. Previously, we mapped quantitative trait loci (QTLs) affecting variation in longevity between two Drosophila melanogaster strains. Here, we show that the longevity QTL in the 36E;38B cytogenetic interval on chromosome 2 contains multiple closely linked QTLs, including the Dopa decarboxylase (Ddc) locus. Complementation tests to mutations show that Ddc is a positional candidate gene for life span in these strains. Linkage disequilibrium (LD) mapping in a sample of 173 alleles from a single population shows that three common molecular polymorphisms in Ddc account for 15.5% of the genetic contribution to variance in life span from chromosome 2. The polymorphisms are in strong LD, and the effects of the haplotypes on longevity suggest that the polymorphisms are maintained by balancing selection. DDC catalyzes the final step in the synthesis of the neurotransmitters, dopamine and serotonin. Thus, these data implicate variation in the synthesis of bioamines as a factor contributing to natural variation in individual life span.

Drosophila bristles and the nature of quantitative genetic variation

Numbers of Drosophila sensory bristles present an ideal model system to elucidate the genetic basis of variation for quantitative traits. Here, we review recent evidence that the genetic architecture of variation for bristle numbers is surprisingly complex. A substantial fraction of the Drosophila genome affects bristle number, indicating pervasive pleiotropy of genes that affect quantitative traits. Further, a large number of loci, often with sex- and environment-specific effects that are also conditional on background genotype, affect natural variation in bristle number. Despite this complexity, an understanding of the molecular basis of natural variation in bristle number is emerging from linkage disequilibrium mapping studies of individual candidate genes that affect the development of sensory bristles. We show that there is naturally segregating genetic variance for environmental plasticity of abdominal and sternopleural bristle number. For abdominal bristle number this variance can be attributed in part to an abnormal abdomen-like phenotype that resembles the phenotype of mutants defective in catecholamine biosynthesis. Dopa decarboxylase (Ddc) encodes the enzyme that catalyses the final step in the synthesis of dopamine, a major Drosophila catecholamine and neurotransmitter. We found that molecular polymorphisms at Ddc are indeed associated with variation in environmental plasticity of abdominal bristle number.

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.

Co-regulated transcriptional networks contribute to natural genetic variation in Drosophila sleep

Sleep disorders are common in humans, and sleep loss increases the risk of obesity and diabetes. Studies in Drosophila have revealed molecular pathways and neural tissues regulating sleep; however, genes that maintain genetic variation for sleep in natural populations are unknown. Here, we characterized sleep in 40 wild-derived Drosophila lines and observed abundant genetic variation in sleep architecture. We associated sleep with genome-wide variation in gene expression to identify candidate genes. We independently confirmed that molecular polymorphisms in Catsup (Catecholamines up) are associated with variation in sleep and that P-element mutations in four candidate genes affect sleep and gene expression. Transcripts associated with sleep grouped into biologically plausible genetically correlated transcriptional modules. We confirmed co-regulated gene expression using P-element mutants. Quantitative genetic analysis of natural phenotypic variation is an efficient method for revealing candidate genes and pathways.

Linkage disequilibrium mapping of molecular polymorphisms at the scabrous locus associated with naturally occurring variation in bristle number in Drosophila melanogaster.

We evaluated the hypothesis that the Drosophila melanogaster second chromosome gene scabrous (sca), a candidate sensory bristle number quantitative trait locus (QTL), contributes to naturally occurring variation in bristle number. Variation in abdominal and sternopleural bristle number was quantified for wild-derived sca alleles in seven genetic backgrounds: as homozygous second chromosomes (C2) in an isogenic background, homozygous lines in which approximately 20 cM including the sca locus had been introgressed into the isogenic background (sca BC), as C2 and sca BC heterozygotes and hemizygotes against a P element insertional sca allele and a P-induced sca deficiency in the same isogenic background, and as sca BC heterozygotes against the wild-type sca allele of isogenic strain. Molecular restriction map variation was determined for a 45 kb region including the sca locus, and single-stranded conformational polymorphism (SSCP) was examined for the third intron and parts of the third and fourth exons. Associations between each of the 27 molecular polymorphisms and bristle number were evaluated within each genotype and on the first principal component score determined from all seven genotypes, separately for each sex and bristle trait. Permutation tests were used to assess the empirical significance thresholds, accounting for multiple, correlated tests, and correlated markers. Three sites in regulatory regions were associated with female-specific variation in abdominal bristle number, one of which was an SSCP site in the region of the gene associated with regulation of sca in embryonic abdominal segments.

Polygenic mutation in Drosophila melanogaster : estimates from divergence among inbred strains

A highly inbred line of Drosophila melanogaster was subdivided into 25 replicate sublines, which were independently maintained for 100 generations with 10 pairs of unselected flies per generation. The polygenic mutation rate (VM) for two quantitative traits, abdominal and sternopleural bristle number, was estimated from divergence among sublines at 10 generation intervals from generations 30‐100, and from response of each line to divergent selection after more than 65 generations of mutation accumulation. Estimates of VM averaged over males and females both from divergence among lines and from response to selection within lines were 3.3 × 10‐3 VE for abdominal bristles and 1.5 × 10‐3 VE for sternopleural bristles, where VE is the environmental variance. The actual rate of production of mutations affecting these traits may be considerably higher if the traits are under stabilizing selection, and if mutations affecting bristle number have deleterious effects on fitness. There was a substantial component of variance for sex × mutant effect interaction and the sublines evolved highly significant mutational variation in sex dimorphism of abdominal bristle number. Pleiotropic effects on sex dimorphism may be a general property of mutations at loci determining bristle number.

Candidate genes affecting Drosophila life span identified by integrating microarray gene expression analysis and QTL mapping

The current increase in life expectancy observed in industrialized societies underscores the need to achieve a better understanding of the aging process that could help the development of effective strategies to achieve healthy aging. This will require not only identifying genes involved in the aging process, but also understanding how their effects are modulated by environmental factors, such as dietary intake and life style. Although the human genome has been sequenced, it may be impractical to study humans or other long-lived organisms to gain a mechanistic understanding about the aging process. Thus, short-lived animal models are essential to identifying the mechanisms and genes that affect the rate and quality of aging as a first step towards identifying genetic variants in humans. In this study, we investigated gene expression changes between two strains of Drosophila (Oregon and 2b) for which quantitative trait loci (QTLs) affecting life span were identified previously. We collected males and females from both strains at young and old ages, and assessed whole genome variation in transcript abundance using Affymetrix GeneChips. We observed 8217 probe sets with detectable transcripts. A total of 2371 probe sets, representing 2220 genes, exhibited significant changes in transcript abundance with age; and 839 probe sets were differentially expressed between Oregon and 2b. We focused on the 359 probe sets (representing 354 genes) that exhibited significant changes in gene expression both with age and between strains. We used these genes to integrate the analysis of microarray gene expression data, bioinformatics, and the results of genetic mapping studies reported previously, to identify 49 candidate genes and four pathways that could potentially be responsible for regulating life span and involved in the process of aging in Drosophila and humans.