Jennifer K. Grenier

Behavioral idiosyncrasy reveals genetic control of phenotypic variability

Significance If we could rear genetically identical individuals from a variety of genetic backgrounds and rear them in the same environment, how much phenotypic variation between individuals of the same genotype would we see? Would different genetic backgrounds differ in their degree of variability? What would account for these differences? We used Drosophila inbred lines to address these questions focusing on variability in locomotor handedness. We show that different genotypes vary dramatically in their propensity for variability, that phenotypic variability itself, as a trait, can be heritable, and that loci affecting variability can be mapped. The genetic control of variability has received little attention in quantitative genetics despite the important role variability plays in explaining phenotypic variation between individuals. Quantitative genetics has primarily focused on describing genetic effects on trait means and largely ignored the effect of alternative alleles on trait variability, potentially missing an important axis of genetic variation contributing to phenotypic differences among individuals. To study the genetic effects on individual-to-individual phenotypic variability (or intragenotypic variability), we used Drosophila inbred lines and measured the spontaneous locomotor behavior of flies walking individually in Y-shaped mazes, focusing on variability in locomotor handedness, an assay optimized to measure variability. We discovered that some lines had consistently high levels of intragenotypic variability among individuals, whereas lines with low variability behaved as although they tossed a coin at each left/right turn decision. We demonstrate that the degree of variability is itself heritable. Using a genome-wide association study (GWAS) for the degree of intragenotypic variability as the phenotype across lines, we identified several genes expressed in the brain that affect variability in handedness without affecting the mean. One of these genes, Ten-a, implicates a neuropil in the central complex of the fly brain as influencing the magnitude of behavioral variability, a brain region involved in sensory integration and locomotor coordination. We validated these results using genetic deficiencies, null alleles, and inducible RNAi transgenes. Our study reveals the constellation of phenotypes that can arise from a single genotype and shows that different genetic backgrounds differ dramatically in their propensity for phenotypic variabililty. Because traditional mean-focused GWASs ignore the contribution of variability to overall phenotypic variation, current methods may miss important links between genotype and phenotype.

Global Diversity Lines–A Five-Continent Reference Panel of Sequenced Drosophila melanogaster Strains

Reference collections of multiple Drosophila lines with accumulating collections of “omics” data have proven especially valuable for the study of population genetics and complex trait genetics. Here we present a description of a resource collection of 84 strains of Drosophila melanogaster whose genome sequences were obtained after 12 generations of full-sib inbreeding. The initial rationale for this resource was to foster development of a systems biology platform for modeling metabolic regulation by the use of natural polymorphisms as perturbations. As reference lines, they are amenable to repeated phenotypic measurements, and already a large collection of metabolic traits have been assayed. Another key feature of these strains is their widespread geographic origin, coming from Beijing, Ithaca, Netherlands, Tasmania, and Zimbabwe. After obtaining 12.5× coverage of paired-end Illumina sequence reads, SNP and indel calls were made with the GATK platform. Thorough quality control was enabled by deep sequencing one line to >100×, and single-nucleotide polymorphisms and indels were validated using ddRAD-sequencing as an orthogonal platform. In addition, a series of preliminary population genetic tests were performed with these single-nucleotide polymorphism data for assessment of data quality. We found 83 segregating inversions among the lines, and as expected these were especially abundant in the African sample. We anticipate that this will make a useful addition to the set of reference D. melanogaster strains, thanks to its geographic structuring and unusually high level of genetic diversity.

Correlated variation and population differentiation in satellite DNA abundance among lines of Drosophila melanogaster.

Significance Most eukaryotic genomes harbor large amounts of highly repetitive satellite DNA primarily in centromeric regions. Closely related Drosophila species have nearly complete turnover of the types and quantities of simple sequence repeats. However, the detailed dynamics of turnover remains unclear, in part due to technical challenges in examining these highly repetitive sequences. We present a method (k-Seek) that identifies and quantifies simple sequence repeats from whole genome sequences. By characterizing natural variation in tandem repeats within Drosophila melanogaster, we identified many novel repeats and found that geographically isolated populations show differentiation patterns that are, unexpectedly, incongruous with demographic history. Moreover, repeats undergo correlated change in abundance, providing additional insight into the dynamics of satellite DNA and genome evolution. Tandemly repeating satellite DNA elements in heterochromatin occupy a substantial portion of many eukaryotic genomes. Although often characterized as genomic parasites deleterious to the host, they also can be crucial for essential processes such as chromosome segregation. Adding to their interest, satellite DNA elements evolve at high rates; among Drosophila, closely related species often differ drastically in both the types and abundances of satellite repeats. However, due to technical challenges, the evolutionary mechanisms driving this rapid turnover remain unclear. Here we characterize natural variation in simple-sequence repeats of 2–10 bp from inbred Drosophila melanogaster lines derived from multiple populations, using a method we developed called k-Seek that analyzes unassembled Illumina sequence reads. In addition to quantifying all previously described satellite repeats, we identified many novel repeats of low to medium abundance. Many of the repeats show population differentiation, including two that are present in only some populations. Interestingly, the population structure inferred from overall satellite quantities does not recapitulate the expected population relationships based on the demographic history of D. melanogaster. We also find that some satellites of similar sequence composition are correlated across lines, revealing concerted evolution. Moreover, correlated satellites tend to be interspersed with each other, further suggesting that concerted change is partially driven by higher order structure. Surprisingly, we identified negative correlations among some satellites, suggesting antagonistic interactions. Our study demonstrates that current genome assemblies vastly underestimate the complexity, abundance, and variation of highly repetitive satellite DNA and presents approaches to understand their rapid evolutionary divergence.

Evidence for the fixation of gene duplications by positive selection in Drosophila

Gene duplications play a key role in the emergence of novel traits and in adaptation. But despite their centrality to evolutionary processes, it is still largely unknown how new gene duplicates are initially fixed within populations and later maintained in genomes. Long-standing debates on the evolution of gene duplications could be settled by determining the relative importance of genetic drift vs. positive selection in the fixation of new gene duplicates. Using the Drosophila Global Diversity Lines (GDL), we have combined genome-wide SNP polymorphism data with a novel set of copy number variant calls and gene expression profiles to characterize the polymorphic phase of new genes. We found that approximately half of the roughly 500 new complete gene duplications segregating in the GDL lead to significant increases in the expression levels of the duplicated genes and that these duplications are more likely to be found at lower frequencies, suggesting a negative impact on fitness. However, we also found that six of the nine gene duplications that are fixed or close to fixation in at least one of the five populations in our study show signs of being under positive selection, and that these duplications are likely beneficial because of dosage effects, with a possible role for additional mutations in two duplications. Our work suggests that in Drosophila, theoretical models that posit that gene duplications are immediately beneficial and fixed by positive selection are most relevant to explain the long-term evolution of gene duplications in this species.

Extensive local adaptation within the chemosensory system following Drosophila melanogaster's global expansion

How organisms adapt to new environments is of fundamental biological interest, but poorly understood at the genetic level. Chemosensory systems provide attractive models to address this problem, because they lie between external environmental signals and internal physiological responses. To investigate how selection has shaped the well-characterized chemosensory system of Drosophila melanogaster, we have analysed genome-wide data from five diverse populations. By couching population genomic analyses of chemosensory protein families within parallel analyses of other large families, we demonstrate that chemosensory proteins are not outliers for adaptive divergence between species. However, chemosensory families often display the strongest genome-wide signals of recent selection within D. melanogaster. We show that recent adaptation has operated almost exclusively on standing variation, and that patterns of adaptive mutations predict diverse effects on protein function. Finally, we provide evidence that chemosensory proteins have experienced relaxed constraint, and argue that this has been important for their rapid adaptation over short timescales.

Dgcr8 and Dicer are essential for sex chromosome integrity during meiosis in males

ABSTRACT Small RNAs play crucial roles in regulating gene expression during mammalian meiosis. To investigate the function of microRNAs (miRNAs) and small interfering RNAs (siRNAs) during meiosis in males, we generated germ-cell-specific conditional deletions of Dgcr8 and Dicer in mice. Analysis of spermatocytes from both conditional knockout lines revealed that there were frequent chromosomal fusions during meiosis, always involving one or both sex chromosomes. RNA sequencing indicates upregulation of Atm in spermatocytes from miRNA-deficient mice, and immunofluorescence imaging demonstrates an increased abundance of activated ATM kinase and mislocalization of phosphorylated MDC1, an ATM phosphorylation substrate. The Atm 3′UTR contains many potential microRNA target sites, and, notably, target sites for several miRNAs depleted in both conditional knockout mice were highly effective at promoting repression. RNF8, a telomere-associated protein whose localization is controlled by the MDC1–ATM kinase cascade, normally associates with the sex chromosomes during pachytene, but in both conditional knockouts redistributed to the autosomes. Taken together, these results suggest that Atm dysregulation in microRNA-deficient germ lines contributes to the redistribution of proteins involved in chromosomal stability from the sex chromosomes to the autosomes, resulting in sex chromosome fusions during meiotic prophase I. Highlighted Article: miRNA-deficient spermatocytes display frequent sex chromosome fusions and fail to progress through meiosis in a process that is probably mediated by dysregulation of Atm.

Survey of Global Genetic Diversity Within the Drosophila Immune System.

Numerous studies across a wide range of taxa have demonstrated that immune genes are routinely among the most rapidly evolving genes in the genome. This observation, however, does not address what proportion of immune genes undergo strong selection during adaptation to novel environments. Here, we determine the extent of very recent divergence in genes with immune function across five populations of Drosophila melanogaster and find that immune genes do not show an overall trend of recent rapid adaptation. Our population-based approach uses a set of carefully matched control genes to account for the effects of demography and local recombination rate, allowing us to identify whether specific immune functions are putative targets of strong selection. We find evidence that viral-defense genes are rapidly evolving in Drosophila at multiple timescales. Local adaptation to bacteria and fungi is less extreme and primarily occurs through changes in recognition and effector genes rather than large-scale changes to the regulation of the immune response. Surprisingly, genes in the Toll pathway, which show a high rate of adaptive substitution between the D. melanogaster and D. simulans lineages, show little population differentiation. Quantifying the flies for resistance to a generalist Gram-positive bacterial pathogen, we found that this genetic pattern of low population differentiation was recapitulated at the phenotypic level. In sum, our results highlight the complexity of immune evolution and suggest that Drosophila immune genes do not follow a uniform trajectory of strong directional selection as flies encounter new environments.

Melanocyte Stem Cell Activation and Translocation Initiate Cutaneous Melanoma in Response to UV Exposure

Melanoma is one of the deadliest cancers, yet the cells of origin and mechanisms of tumor initiation remain unclear. The majority of melanomas emerge from clear skin without a precursor lesion, but it is unknown whether these melanomas can arise from melanocyte stem cells (MCSCs). Here we employ mouse models to define the role of MCSCs as melanoma cells of origin, demonstrate that MCSC quiescence acts as a tumor suppressor, and identify the extrinsic environmental and molecular factors required for the critical early steps of melanoma initiation. Specifically, melanomas originate from melanoma-competent MCSCs upon stimulation by UVB, which induces MCSC activation and translocation via an inflammation-dependent process. Moreover, the chromatin-remodeling factor Hmga2 in the skin plays a critical role in UVB-mediated melanomagenesis. These findings delineate melanoma formation from melanoma-competent MCSCs following extrinsic stimuli, and they suggest that abrogation of Hmga2 function in the microenvironment can suppress MCSC-originating cutaneous melanomas.

Developmental Origin Governs CD8+ T Cell Fate Decisions during Infection

Heterogeneity is a hallmark feature of the adaptive immune system in vertebrates. Following infection, naive Txa0cells differentiate into various subsets of effector and memory Txa0cells, which help to eliminate pathogens and maintain long-term immunity. The current model suggests there is a single lineage of naive Txa0cells that give rise to different populations of effector and memory Txa0cells depending on the type and amounts of stimulation they encounter during infection. Here, we have discovered that multiple sub-populations of cells exist in the naive CD8+ Txa0cell pool that are distinguished by their developmental origin, unique transcriptional profiles, distinct chromatin landscapes, and different kinetics and phenotypes after microbial challenge. These data demonstrate that the naive CD8+ Txa0cell pool is not as homogeneous as previously thought and offers a new framework for explaining the remarkable heterogeneity in the effector and memory Txa0cell subsets that arise after infection.

Transcriptome profiling of alphaherpesvirus-infected cells treated with the HIV-integrase inhibitor raltegravir reveals profound and specific alterations in host transcription

Anti-microbial compounds typically exert their action by directly interfering with one or more stages of the pathogens life cycle. However, some compounds also have secondary effects on the host that aid in pathogen clearance. Raltegravir is a human immunodeficiency virus (HIV)-integrase inhibitor that has been shown to alter the host immune response to HIV in addition to its direct antiviral effect. Interestingly, raltegravir can also directly inhibit the replication of various herpesviruses. However, the host-targeted effects of this drug in the context of a herpesvirus infection have not been explored. Here, we used felid alphaherpesvirus 1 (FHV-1), a close relative of human alphaherpesvirus 1 (HHV-1) that similarly causes ocular herpes, to characterize the host-targeted effects of raltegravir on corneal epithelial cells during an alphaherpesvirus infection. Using RNA deep sequencing, we found that raltegravir specifically boosts the expression of anti-angiogenic factors and promotes metabolic homeostasis in FHV-1-infected cells. In contrast, few changes in host gene transcription were found in uninfected cells. Importantly, we were able to demonstrate that these effects were specific to raltegravir and independent of the direct-acting antiviral effect of the drug, since treatment with the DNA polymerase inhibitor phosphonoacetic acid did not induce these host-targeted effects. Taken together, these results indicate that raltegravir has profound and specific effects on the host transcription profile of herpesvirus-infected cells that may contribute to the overall antiviral activity of the drug and could provide therapeutic benefits in vivo. Furthermore, this study provides a framework for future efforts evaluating the host-targeted effects of anti-microbial compounds.