Genetic variation in chromatin state across multiple tissues in Drosophila melanogaster.

We use ATAC-seq to examine chromatin accessibility for four different tissues in Drosophila melanogaster: adult female brain, ovaries, and both wing and eye-antennal imaginal discs from males. Each tissue is assayed in eight different inbred strain genetic backgrounds, seven associated with a reference quality genome assembly. We develop a method for the quantile normalization of ATAC-seq fragments and test for differences in coverage among genotypes, tissues, and their interaction at 44099 peaks throughout the euchromatic genome. For the strains with reference quality genome assemblies, we correct ATAC-seq profiles for read mis-mapping due to nearby polymorphic structural variants (SVs). Comparing coverage among genotypes without accounting for SVs results in a highly elevated rate (55%) of identifying false positive differences in chromatin state between genotypes. After SV correction, we identify 1050, 30383, and 4508 regions whose peak heights are polymorphic among genotypes, among tissues, or exhibit genotype-by-tissue interactions, respectively. Finally, we identify 3988 candidate causative variants that explain at least 80% of the variance in chromatin state at nearby ATAC-seq peaks.

Genetic variation in P-element dysgenic sterility is associated with double-strand break repair and alternative splicing of TE transcripts.

The germline mobilization of transposable elements (TEs) by small RNA mediated silencing pathways is conserved across eukaryotes and critical for ensuring the integrity of gamete genomes. However, genomes are recurrently invaded by novel TEs through horizontal transfer. These invading TEs are not targeted by host small RNAs, and their unregulated activity can cause DNA damage in germline cells and ultimately lead to sterility. Here we use hybrid dysgenesis-a sterility syndrome of Drosophila caused by transposition of invading P-element DNA transposons-to uncover host genetic variants that modulate dysgenic sterility. Using a panel of highly recombinant inbred lines of Drosophila melanogaster, we identified two linked quantitative trait loci (QTL) that determine the severity of dysgenic sterility in young and old females, respectively. We show that ovaries of fertile genotypes exhibit increased expression of splicing factors that suppress the production of transposase encoding transcripts, which likely reduces the transposition rate and associated DNA damage. We also show that fertile alleles are associated with decreased sensitivity to double-stranded breaks and enhanced DNA repair, explaining their ability to withstand high germline transposition rates. Together, our work reveals a diversity of mechanisms whereby host genotype modulates the cost of an invading TE, and points to genetic variants that were likely beneficial during the P-element invasion.

Predicting evolutionary change at the DNA level in a natural Mimulus population.

Evolution by natural selection occurs when the frequencies of genetic variants change because individuals differ in Darwinian fitness components such as survival or reproductive success. Differential fitness has been demonstrated in field studies of many organisms, but it remains unclear how well we can quantitatively predict allele frequency changes from fitness measurements. Here, we characterize natural selection on millions of Single Nucleotide Polymorphisms (SNPs) across the genome of the annual plant Mimulus guttatus. We use fitness estimates to calibrate population genetic models that effectively predict allele frequency changes into the next generation. Hundreds of SNPs experienced "male selection" in 2013 with one allele at each SNP elevated in frequency among successful male gametes relative to the entire population of adults. In the following generation, allele frequencies at these SNPs consistently shifted in the predicted direction. A second year of study revealed that SNPs had effects on both viability and reproductive success with pervasive trade-offs between fitness components. SNPs favored by male selection were, on average, detrimental to survival. These trade-offs (antagonistic pleiotropy and temporal fluctuations in fitness) may be essential to the long-term maintenance of alleles. Despite the challenges of measuring selection in the wild, the strong correlation between predicted and observed allele frequency changes suggests that population genetic models have a much greater role to play in forward-time prediction of evolutionary change.

Herpes Simplex Virus 1 Mutant with Point Mutations in UL39 Is Impaired for Acute Viral Replication in Mice, Establishment of Latency, and Explant-Induced Reactivation.

In the process of generating herpes simplex virus 1 (HSV-1) mutations in the viral regulatory gene encoding infected cell protein 0 (ICP0), we isolated a viral mutant, termed KOS-NA, that was severely impaired for acute replication in the eyes and trigeminal ganglia (TG) of mice, defective in establishing a latent infection, and reactivated poorly from explanted TG. To identify the secondary mutation(s) responsible for the impaired phenotypes of this mutant, we sequenced the KOS-NA genome and noted that it contained two nonsynonymous mutations in UL39, which encodes the large subunit of ribonucleotide reductase, ICP6. These mutations resulted in lysine-to-proline (residue 393) and arginine-to-histidine (residue 950) substitutions in ICP6. To determine whether alteration of these amino acids was responsible for the KOS-NA phenotypes in vivo, we recombined the wild-type UL39 gene into the KOS-NA genome and rescued its acute replication phenotypes in mice. To further establish the role of UL39 in KOS-NA's decreased pathogenicity, the UL39 mutations were recombined into HSV-1 (generating UL39mut), and this mutant virus showed reduced ocular and TG replication in mice comparable to that of KOS-NA. Interestingly, ICP6 protein levels were reduced in KOS-NA-infected cells relative to the wild-type protein. Moreover, we observed that KOS-NA does not counteract caspase 8-induced apoptosis, unlike wild-type strain KOS. Based on alignment studies with other HSV-1 ICP6 homologs, our data suggest that amino acid 950 of ICP6 likely plays an important role in ICP6 accumulation and inhibition of apoptosis, consequently impairing HSV-1 pathogenesis in a mouse model of HSV-1 infection.IMPORTANCE HSV-1 is a major human pathogen that infects ∼80% of the human population and can be life threatening to infected neonates or immunocompromised individuals. Effective therapies for treatment of recurrent HSV-1 infections are limited, which emphasizes a critical need to understand in greater detail the events that modulate HSV-1 replication and pathogenesis. In the current study, we identified a neuroattenuated HSV-1 mutant (i.e., KOS-NA) that contains novel mutations in the UL39 gene, which codes for the large subunit of ribonucleotide reductase (also known as ICP6). This mutant form of ICP6 was responsible for the attenuation of KOS-NA in vivo and resulted in diminished ICP6 protein levels and antiapoptotic effect. Thus, we have determined that subtle alteration of the UL39 gene regulates expression and functions of ICP6 and severely impacts HSV-1 pathogenesis, potentially making KOS-NA a promising vaccine candidate against HSV-1.

Identifying Loci Contributing to Natural Variation in Xenobiotic Resistance in Drosophila.

Natural populations exhibit a great deal of interindividual genetic variation in the response to toxins, exemplified by the variable clinical efficacy of pharmaceutical drugs in humans, and the evolution of pesticide resistant insects. Such variation can result from several phenomena, including variable metabolic detoxification of the xenobiotic, and differential sensitivity of the molecular target of the toxin. Our goal is to genetically dissect variation in the response to xenobiotics, and characterize naturally-segregating polymorphisms that modulate toxicity. Here, we use the Drosophila Synthetic Population Resource (DSPR), a multiparent advanced intercross panel of recombinant inbred lines, to identify QTL (Quantitative Trait Loci) underlying xenobiotic resistance, and employ caffeine as a model toxic compound. Phenotyping over 1,700 genotypes led to the identification of ten QTL, each explaining 4.5-14.4% of the broad-sense heritability for caffeine resistance. Four QTL harbor members of the cytochrome P450 family of detoxification enzymes, which represent strong a priori candidate genes. The case is especially strong for Cyp12d1, with multiple lines of evidence indicating the gene causally impacts caffeine resistance. Cyp12d1 is implicated by QTL mapped in both panels of DSPR RILs, is significantly upregulated in the presence of caffeine, and RNAi knockdown robustly decreases caffeine tolerance. Furthermore, copy number variation at Cyp12d1 is strongly associated with phenotype in the DSPR, with a trend in the same direction observed in the DGRP (Drosophila Genetic Reference Panel). No additional plausible causative polymorphisms were observed in a full genomewide association study in the DGRP, or in analyses restricted to QTL regions mapped in the DSPR. Just as in human populations, replicating modest-effect, naturally-segregating causative variants in an association study framework in flies will likely require very large sample sizes.

Dissecting complex traits using the Drosophila Synthetic Population Resource.

For most complex traits we have a poor understanding of the positions, phenotypic effects, and population frequencies of the underlying genetic variants contributing to their variation. Recently, several groups have developed multi-parent advanced intercross mapping panels in different model organisms in an attempt to improve our ability to characterize causative genetic variants. These panels are powerful and are particularly well suited to the dissection of phenotypic variation generated by rare alleles and loci segregating multiple functional alleles. We describe studies using one such panel, the Drosophila Synthetic Population Resource (DSPR), and the implications for our understanding of the genetic basis of complex traits. In particular, we note that many loci of large effect appear to be multiallelic. If multiallelism is a general rule, analytical approaches designed to identify multiallelic variants should be a priority for both genome-wide association studies (GWASs) and multi-parental panels.

Can asexuality confer a short-term advantage? Investigating apparent biogeographic success in the apomictic triploid fern Myriopteris gracilis.

PREMISE OF THE STUDY: Although asexual taxa are generally seen as evolutionary dead ends, asexuality appears to provide a short-term benefit in some taxa, including a wider geographic distribution compared to sexual relatives. However, this may be an illusion created by multiple, morphologically cryptic, asexual lineages, each occupying a relatively small area. In this study we investigate the role of multiple lineages in the biogeography of Myriopteris gracilis Fée (Pteridaceae), a North American apomictic triploid fern species with a particularly large range. METHODS: Range-wide asexuality was assessed by counting spores/sporangium in 606 Myriopteris gracilis specimens from across the species range, and lineage structure was assessed with both plastid DNA sequence and Genotyping By Sequencing (GBS) SNP datasets. KEY RESULTS: Spore counting of >600 specimens identified no sexual populations, establishing that Myriopteris gracilis is exclusively asexual. The plastid data estimated the crown age of M. gracilis at ca. 2.5 mya and identified two lineages, each largely confined to the eastern or western portions of the range. These groups were further subdivided by the GBS data, revealing at least seven asexual lineages of varying geographic distributions, each occupying a relatively small portion of the total range of M. gracilis. CONCLUSIONS: Although maintained exclusively through asexual reproduction, the broad distribution of Myriopteris gracilis is a compilation of numerous, independently formed asexual lineages. Since no single asexual lineage occupies the full extent of the species distribution, recurrent lineage formation should be considered when evaluating the short-term benefit of asexuality in this taxon and others.

Genetic dissection of the Drosophila melanogaster female head transcriptome reveals widespread allelic heterogeneity.

Modern genetic mapping is plagued by the "missing heritability" problem, which refers to the discordance between the estimated heritabilities of quantitative traits and the variance accounted for by mapped causative variants. One major potential explanation for the missing heritability is allelic heterogeneity, in which there are multiple causative variants at each causative gene with only a fraction having been identified. The majority of genome-wide association studies (GWAS) implicitly assume that a single SNP can explain all the variance for a causative locus. However, if allelic heterogeneity is prevalent, a substantial amount of genetic variance will remain unexplained. In this paper, we take a haplotype-based mapping approach and quantify the number of alleles segregating at each locus using a large set of 7922 eQTL contributing to regulatory variation in the Drosophila melanogaster female head. Not only does this study provide a comprehensive eQTL map for a major community genetic resource, the Drosophila Synthetic Population Resource, but it also provides a direct test of the allelic heterogeneity hypothesis. We find that 95% of cis-eQTLs and 78% of trans-eQTLs are due to multiple alleles, demonstrating that allelic heterogeneity is widespread in Drosophila eQTL. Allelic heterogeneity likely contributes significantly to the missing heritability problem common in GWAS studies.

Genetic analysis of variation in lifespan using a multiparental advanced intercross Drosophila mapping population.

BACKGROUND: Considerable natural variation for lifespan exists within human and animal populations. Genetically dissecting this variation can elucidate the pathways and genes involved in aging, and help uncover the genetic mechanisms underlying risk for age-related diseases. Studying aging in model systems is attractive due to their relatively short lifespan, and the ability to carry out programmed crosses under environmentally-controlled conditions. Here we investigate the genetic architecture of lifespan using the Drosophila Synthetic Population Resource (DSPR), a multiparental advanced intercross mapping population. RESULTS: We measured lifespan in females from 805 DSPR lines, mapping five QTL (Quantitative Trait Loci) that each contribute 4-5 % to among-line lifespan variation in the DSPR. Each of these QTL co-localizes with the position of at least one QTL mapped in 13 previous studies of lifespan variation in flies. However, given that these studies implicate >90 % of the genome in the control of lifespan, this level of overlap is unsurprising. DSPR QTL intervals harbor 11-155 protein-coding genes, and we used RNAseq on samples of young and old flies to help resolve pathways affecting lifespan, and identify potentially causative loci present within mapped QTL intervals. Broad age-related patterns of expression revealed by these data recapitulate results from previous work. For example, we see an increase in antimicrobial defense gene expression with age, and a decrease in expression of genes involved in the electron transport chain. Several genes within QTL intervals are highlighted by our RNAseq data, such as Relish, a critical immune response gene, that shows increased expression with age, and UQCR-14, a gene involved in mitochondrial electron transport, that has reduced expression in older flies. CONCLUSIONS: The five QTL we isolate collectively explain a considerable fraction of the genetic variation for female lifespan in the DSPR, and implicate modest numbers of genes. In several cases the candidate loci we highlight reside in biological pathways already implicated in the control of lifespan variation. Thus, our results provide further evidence that functional genetics tests targeting these genes will be fruitful, lead to the identification of natural sequence variants contributing to lifespan variation, and help uncover the mechanisms of aging.

Genetic dissection of a model complex trait using the Drosophila Synthetic Population Resource.

Genetic dissection of complex, polygenic trait variation is a key goal of medical and evolutionary genetics. Attempts to identify genetic variants underlying complex traits have been plagued by low mapping resolution in traditional linkage studies, and an inability to identify variants that cumulatively explain the bulk of standing genetic variation in genome-wide association studies (GWAS). Thus, much of the heritability remains unexplained for most complex traits. Here we describe a novel, freely available resource for the Drosophila community consisting of two sets of recombinant inbred lines (RILs), each derived from an advanced generation cross between a different set of eight highly inbred, completely resequenced founders. The Drosophila Synthetic Population Resource (DSPR) has been designed to combine the high mapping resolution offered by multiple generations of recombination, with the high statistical power afforded by a linkage-based design. Here, we detail the properties of the mapping panel of >1600 genotyped RILs, and provide an empirical demonstration of the utility of the approach by genetically dissecting alcohol dehydrogenase (ADH) enzyme activity. We confirm that a large fraction of the variation in this classic quantitative trait is due to allelic variation at the Adh locus, and additionally identify several previously unknown modest-effect trans-acting QTL (quantitative trait loci). Using a unique property of multiparental linkage mapping designs, for each QTL we highlight a relatively small set of candidate causative variants for follow-up work. The DSPR represents an important step toward the ultimate goal of a complete understanding of the genetics of complex traits in the Drosophila model system.

Abundance and distribution of transposable elements in two Drosophila QTL mapping resources.

Here we present computational machinery to efficiently and accurately identify transposable element (TE) insertions in 146 next-generation sequenced inbred strains of Drosophila melanogaster. The panel of lines we use in our study is composed of strains from a pair of genetic mapping resources: the Drosophila Genetic Reference Panel (DGRP) and the Drosophila Synthetic Population Resource (DSPR). We identified 23,087 TE insertions in these lines, of which 83.3% are found in only one line. There are marked differences in the distribution of elements over the genome, with TEs found at higher densities on the X chromosome, and in regions of low recombination. We also identified many more TEs per base pair of intronic sequence and fewer TEs per base pair of exonic sequence than expected if TEs are located at random locations in the euchromatic genome. There was substantial variation in TE load across genes. For example, the paralogs derailed and derailed-2 show a significant difference in the number of TE insertions, potentially reflecting differences in the selection acting on these loci. When considering TE families, we find a very weak effect of gene family size on TE insertions per gene, indicating that as gene family size increases the number of TE insertions in a given gene within that family also increases. TEs are known to be associated with certain phenotypes, and our data will allow investigators using the DGRP and DSPR to assess the functional role of TE insertions in complex trait variation more generally. Notably, because most TEs are very rare and often private to a single line, causative TEs resulting in phenotypic differences among individuals may typically fail to replicate across mapping panels since individual elements are unlikely to segregate in both panels. Our data suggest that "burden tests" that test for the effect of TEs as a class may be more fruitful.

Gene expression variation underlying tissue-specific responses to copper stress in Drosophila melanogaster.

Copper is one of a handful of biologically necessary heavy metals that is also a common environmental pollutant. Under normal conditions, copper ions are required for many key physiological processes. However, in excess, copper results in cell and tissue damage ranging in severity from temporary injury to permanent neurological damage. Because of its biological relevance, and because many conserved copper-responsive genes respond to nonessential heavy metal pollutants, copper resistance in Drosophila melanogaster is a useful model system with which to investigate the genetic control of the heavy metal stress response. Because heavy metal toxicity has the potential to differently impact specific tissues, we genetically characterized the control of the gene expression response to copper stress in a tissue-specific manner in this study. We assessed the copper stress response in head and gut tissue of 96 inbred strains from the Drosophila Synthetic Population Resource using a combination of differential expression analysis and expression quantitative trait locus mapping. Differential expression analysis revealed clear patterns of tissue-specific expression. Tissue and treatment specific responses to copper stress were also detected using expression quantitative trait locus mapping. Expression quantitative trait locus associated with MtnA, Mdr49, Mdr50, and Sod3 exhibited both genotype-by-tissue and genotype-by-treatment effects on gene expression under copper stress, illuminating tissue- and treatment-specific patterns of gene expression control. Together, our data build a nuanced description of the roles and interactions between allelic and expression variation in copper-responsive genes, provide valuable insight into the genomic architecture of susceptibility to metal toxicity, and highlight candidate genes for future functional characterization.

Functional transcriptomic analysis of the role of MAB-5/Hox in Q neuroblast migration in Caenorhabditis elegans.

BACKGROUND: Directed cell migration is a fundamental process in normal development and in tumor metastasis. In C. elegans the MAB-5/Hox transcription factor is a determinant of posterior migration of the Q neuroblast descendants. In this work, mab-5 transcriptional targets that control Q descendant migration are identified by comparing RNA-seq profiles in wild type and mab-5 mutant backgrounds. RESULTS: Transcriptome profiling is a widely-used and potent tool to identify genes involved in developmental and pathological processes, and is most informative when RNA can be isolated from individual cell or tissue types. Cell-specific RNA samples can be difficult to obtain from invertebrate model organisms such as Drosophila and C. elegans. Here we test the utility of combining a whole organism RNA-seq approach with mab-5 loss and gain-of-function mutants and functional validation using RNAi to identify genes regulated by MAB-5 to control Q descendant migration. We identified 22 genes whose expression was controlled by mab-5 and that controlled Q descendant migration. Genes regulated by mab-5 were enriched for secreted and transmembrane molecules involved in basement membrane interaction and modification, and some affected Q descendant migration. CONCLUSIONS: Our results indicate that a whole-organism RNA-seq approach, when combined with mutant analysis and functional validation, can be a powerful method to identify genes involved in a specific developmental process, in this case Q descendant posterior migration. These genes could act either autonomously in the Q cells, or non-autonomously in other cells that express MAB-5. The identities of the genes regulated by MAB-5 indicate that MAB-5 acts by modifying interactions with the basement membrane, resulting in posterior versus anterior migration.

Genome sequence of herpes simplex virus 1 strain KOS.

Herpes simplex virus type 1 (HSV-1) strain KOS has been extensively used in many studies to examine HSV-1 replication, gene expression, and pathogenesis. Notably, strain KOS is known to be less pathogenic than the first sequenced genome of HSV-1, strain 17. To understand the genotypic differences between KOS and other phenotypically distinct strains of HSV-1, we sequenced the viral genome of strain KOS. When comparing strain KOS to strain 17, there are at least 1,024 small nucleotide polymorphisms (SNPs) and 172 insertions/deletions (indels). The polymorphisms observed in the KOS genome will likely provide insights into the genes, their protein products, and the cis elements that regulate the biology of this HSV-1 strain.

Genome sequence of herpes simplex virus 1 strain McKrae.

The herpes simplex virus 1 (HSV-1) strain McKrae is highly virulent compared to other wild-type strains of HSV-1. To help us better understand the genetic determinants that lead to differences in the pathogenicity of McKrae and other HSV-1 strains, we sequenced its genome. Comparing the sequence of McKrae's genome to that of strain 17 revealed that the genomes differ by at least 752 single nucleotide polymorphisms (SNPs) and 86 insertion/deletion events (indels). Although the majority of these polymorphisms reside in noncoding regions, 241 SNPs and 10 indels alter the protein-coding sequences of 58 open reading frames. Some of these variations are expected to contribute to the pathogenic phenotype of McKrae.

Gene expression variation underlying tissue-specific responses to copper stress in Drosophila melanogaster.

Copper is one of a handful of biologically necessary heavy metals that is also a common environmental pollutant. Under normal conditions, copper ions are required for many key physiological processes. However, in excess, copper quickly results in cell and tissue damage that can range in severity from temporary injury to permanent neurological damage. Because of its biological relevance, and because many conserved copper-responsive genes also respond to other non-essential heavy metal pollutants, copper resistance in Drosophila melanogaster is a useful model system with which to investigate the genetic control of the response to heavy metal stress. Because heavy metal toxicity has the potential to differently impact specific tissues, we genetically characterized the control of the gene expression response to copper stress in a tissue-specific manner in this study. We assessed the copper stress response in head and gut tissue of 96 inbred strains from the Drosophila Synthetic Population Resource (DSPR) using a combination of differential expression analysis and expression quantitative trait locus (eQTL) mapping. Differential expression analysis revealed clear patterns of tissue-specific expression, primarily driven by a more pronounced gene expression response in gut tissue. eQTL mapping of gene expression under control and copper conditions as well as for the change in gene expression following copper exposure (copper response eQTL) revealed hundreds of genes with tissue-specific local cis-eQTL and many distant trans-eQTL. eQTL associated with MtnA, Mdr49, Mdr50, and Sod3 exhibited genotype by environment effects on gene expression under copper stress, illuminating several tissue- and treatment-specific patterns of gene expression control. Together, our data build a nuanced description of the roles and interactions between allelic and expression variation in copper-responsive genes, provide valuable insight into the genomic architecture of susceptibility to metal toxicity, and highlight many candidate genes for future functional characterization.

The genetic basis of adaptation to copper pollution in Drosophila melanogaster.

Introduction: Heavy metal pollutants can have long lasting negative impacts on ecosystem health and can shape the evolution of species. The persistent and ubiquitous nature of heavy metal pollution provides an opportunity to characterize the genetic mechanisms that contribute to metal resistance in natural populations. Methods: We examined variation in resistance to copper, a common heavy metal contaminant, using wild collections of the model organism Drosophila melanogaster. Flies were collected from multiple sites that varied in copper contamination risk. We characterized phenotypic variation in copper resistance within and among populations using bulked segregant analysis to identify regions of the genome that contribute to copper resistance. Results and Discussion: Copper resistance varied among wild populations with a clear correspondence between resistance level and historical exposure to copper. We identified 288 SNPs distributed across the genome associated with copper resistance. Many SNPs had population-specific effects, but some had consistent effects on copper resistance in all populations. Significant SNPs map to several novel candidate genes involved in refolding disrupted proteins, energy production, and mitochondrial function. We also identified one SNP with consistent effects on copper resistance in all populations near CG11825, a gene involved in copper homeostasis and copper resistance. We compared the genetic signatures of copper resistance in the wild-derived populations to genetic control of copper resistance in the Drosophila Synthetic Population Resource (DSPR) and the Drosophila Genetic Reference Panel (DGRP), two copper-naïve laboratory populations. In addition to CG11825, which was identified as a candidate gene in the wild-derived populations and previously in the DSPR, there was modest overlap of copper-associated SNPs between the wild-derived populations and laboratory populations. Thirty-one SNPs associated with copper resistance in wild-derived populations fell within regions of the genome that were associated with copper resistance in the DSPR in a prior study. Collectively, our results demonstrate that the genetic control of copper resistance is highly polygenic, and that several loci can be clearly linked to genes involved in heavy metal toxicity response. The mixture of parallel and population-specific SNPs points to a complex interplay between genetic background and the selection regime that modifies the effects of genetic variation on copper resistance.

The genetic control of avascular area in mouse oxygen-induced retinopathy.

PURPOSE: The C57BL/6ByJ and BALB/cByJ inbred strains of mice are, respectively, susceptible and resistant to oxygen-induced retinopathy (OIR). The purpose of this work was to investigate the genetic control of the retinal avascular area in mouse OIR using a mapping cross. METHODS: The central retinal avascular area was measured on postnatal day 16 (P16) in C57BL/6ByJ, BALB/cByJ, 101 (C57BL/6ByJ x BALB/cByJ)F2, and 116 (BALB/cByJ x C57BL/6ByJ)F2 mice that had been subjected to the OIR protocol. A genome-wide scan was performed of selected albino and non-albino mice to determine quantitative trait loci associated with weight and avascular area. RESULTS: C57BL/6ByJ mice had significantly larger avascular areas than BALB/cByJ ones. Albino mice of the F2 generation had smaller avascular areas than the non-albino mice. Genotyping was performed at 856 informative single nucleotide polymorphisms approximately evenly distributed across the genome from each of 85 selected F2 mice. Weight, sex, and the paternal grandmother were found to act as additive covariates associated with the avascular area on P16; mapping analyses that used a model incorporating these covariates found a quantitative trait locus on chromosome 7 related to avascular area. Mapping analyses that used a model that did not incorporate covariates found a quantitative trait locus on chromosome 9 related to avascular area. A quantitative trait locus for bodyweight on P16 was mapped to chromosome 5. CONCLUSIONS: The retinal avascular area in the mouse OIR model is under genetic control. Revascularization in OIR is related to the weight, strain of paternal grandmother, sex, and albinism. Our data support the existence of a quantitative trait locus on chromosome 5 that influences weight after exposure to hyperoxia, as well as quantitative trait loci on chromosomes 7 and 9 that modify susceptibility to OIR.

Pre-adult exposure to three heavy metals leads to changes in the head transcriptome of adult flies.

We examined the effect of developmental exposure to three heavy metals - cadmium, copper, and lead - on gene expression in adult head tissue in the model organism Drosophila melanogaster . All metals affected development time and/or gene expression level. While variation in the response to each metal was apparent, two differentially-expressed genes were upregulated in response to all three metal treatments, and 11 genes were downregulated in two of the three treatments. Our work reveals that developmental metal exposure has the potential to have long-lasting, metal-specific effects on gene expression in adults, even after the metal stress has been removed.

Discovery of malathion resistance QTL in Drosophila melanogaster using a bulked phenotyping approach.

Drosophila melanogaster has proved an effective system with which to understand the evolutionary genetics and molecular mechanisms of insecticide resistance. Insecticide use has left signatures of selection in the fly genome, and both functional and quantitative genetic studies in the system have identified genes and variants associated with resistance. Here, we use D. melanogaster and leverage a bulk phenotyping and pooled sequencing "extreme quantitative trait loci" approach to genetically dissect variation in resistance to malathion, an organophosphate insecticide. We resolve 2 quantitative trait loci, one of which implicates allelic variation at the cytochrome P450 gene Cyp6g1, a strong candidate based on previous work. The second shows no overlap with hits from a previous genome-wide association study for malathion resistance, recapitulating other studies showing that different strategies for complex trait dissection in flies can yield apparently different architectures. Notably, we see no genetic signal at the Ace gene. Ace encodes the target of organophosphate insecticide inhibition, and genome-wide association studies have identified strong Ace-linked associations with resistance in flies. The absence of quantitative trait locus implicating Ace here is most likely because our mapping population does not segregate for several of the known functional polymorphisms impacting resistance at Ace, perhaps because our population is derived from flies collected prior to the widespread use of organophosphate insecticides. Our fundamental approach can be an efficient, powerful strategy to dissect genetic variation in resistance traits. Nonetheless, studies seeking to interrogate contemporary insecticide resistance variation may benefit from deriving mapping populations from more recently collected strains.