Natural variation in gene expression and viral susceptibility revealed by neural progenitor cell villages.

Human genome variation contributes to diversity in neurodevelopmental outcomes and vulnerabilities; recognizing the underlying molecular and cellular mechanisms will require scalable approaches. Here, we describe a "cell village" experimental platform we used to analyze genetic, molecular, and phenotypic heterogeneity across neural progenitor cells from 44 human donors cultured in a shared in vitro environment using algorithms (Dropulation and Census-seq) to assign cells and phenotypes to individual donors. Through rapid induction of human stem cell-derived neural progenitor cells, measurements of natural genetic variation, and CRISPR-Cas9 genetic perturbations, we identified a common variant that regulates antiviral IFITM3 expression and explains most inter-individual variation in susceptibility to the Zika virus. We also detected expression QTLs corresponding to GWAS loci for brain traits and discovered novel disease-relevant regulators of progenitor proliferation and differentiation such as CACHD1. This approach provides scalable ways to elucidate the effects of genes and genetic variation on cellular phenotypes.

Pervasive Regulatory Functions of mRNA Structure Revealed by High-Resolution SHAPE Probing.

mRNAs can fold into complex structures that regulate gene expression. Resolving such structures de novo has remained challenging and has limited our understanding of the prevalence and functions of mRNA structure. We use SHAPE-MaP experiments in living E. coli cells to derive quantitative, nucleotide-resolution structure models for 194 endogenous transcripts encompassing approximately 400 genes. Individual mRNAs have exceptionally diverse architectures, and most contain well-defined structures. Active translation destabilizes mRNA structure in cells. Nevertheless, mRNA structure remains similar between in-cell and cell-free environments, indicating broad potential for structure-mediated gene regulation. We find that the translation efficiency of endogenous genes is regulated by unfolding kinetics of structures overlapping the ribosome binding site. We discover conserved structured elements in 35% of UTRs, several of which we validate as novel protein binding motifs. RNA structure regulates every gene studied here in a meaningful way, implying that most functional structures remain to be discovered.

The genetic basis of parental care evolution in monogamous mice.

Parental care is essential for the survival of mammals, yet the mechanisms underlying its evolution remain largely unknown. Here we show that two sister species of mice, Peromyscus polionotus and Peromyscus maniculatus, have large and heritable differences in parental behaviour. Using quantitative genetics, we identify 12 genomic regions that affect parental care, 8 of which have sex-specific effects, suggesting that parental care can evolve independently in males and females. Furthermore, some regions affect parental care broadly, whereas others affect specific behaviours, such as nest building. Of the genes linked to differences in nest-building behaviour, vasopressin is differentially expressed in the hypothalamus of the two species, with increased levels associated with less nest building. Using pharmacology in Peromyscus and chemogenetics in Mus, we show that vasopressin inhibits nest building but not other parental behaviours. Together, our results indicate that variation in an ancient neuropeptide contributes to interspecific differences in parental care.

Loss of schooling behavior in cavefish through sight-dependent and sight-independent mechanisms.

BACKGROUND: Surface populations of Astyanax mexicanus, living in rivers like their common ancestors, school, while several, independently derived cave populations of the same species have lost schooling behavior. RESULTS: We quantify schooling behavior in individual A. mexicanus and identify quantitative trait loci (QTL) for this trait. We find that the evolutionary modulation of schooling has both vision-dependent and -independent components. We also quantify differences in the lateral line and vision between cavefish and surface fish and relate these differences to the evolutionary loss of schooling behavior. We provide evidence that a monoamine neurotransmitter may have played a role in the evolution of schooling behavior. CONCLUSIONS: We find that vision is essential for schooling tendency in A. mexicanus, while the lateral line has a small effect on this behavior. Schooling behavior in A. mexicanus has evolved both through changes in sensory systems and through changes in genetic loci that likely act downstream of sensory inputs.

Adaptive evolution of multiple traits through multiple mutations at a single gene.

The identification of precise mutations is required for a complete understanding of the underlying molecular and evolutionary mechanisms driving adaptive phenotypic change. Using plasticine models in the field, we show that the light coat color of deer mice that recently colonized the light-colored soil of the Nebraska Sand Hills provides a strong selective advantage against visually hunting predators. Color variation in an admixed population suggests that this light Sand Hills phenotype is composed of multiple traits. We identified distinct regions within the Agouti locus associated with each color trait and found that only haplotypes associated with light trait values have evidence of selection. Thus, local adaptation is the result of independent selection on many mutations within a single locus, each with a specific effect on an adaptive phenotype, thereby minimizing pleiotropic consequences.

Discrete genetic modules are responsible for complex burrow evolution in Peromyscus mice.

Relative to morphological traits, we know little about how genetics influence the evolution of complex behavioural differences in nature. It is unclear how the environment influences natural variation in heritable behaviour, and whether complex behavioural differences evolve through few genetic changes, each affecting many aspects of behaviour, or through the accumulation of several genetic changes that, when combined, give rise to behavioural complexity. Here we show that in nature, oldfield mice (Peromyscus polionotus) build complex burrows with long entrance and escape tunnels, and that burrow length is consistent across populations, although burrow depth varies with soil composition. This burrow architecture is in contrast with the small, simple burrows of its sister species, deer mice (P. maniculatus). When investigated under laboratory conditions, both species recapitulate their natural burrowing behaviour. Genetic crosses between the two species reveal that the derived burrows of oldfield mice are dominant and evolved through the addition of multiple genetic changes. In burrows built by first-generation backcross mice, entrance-tunnel length and the presence of an escape tunnel can be uncoupled, suggesting that these traits are modular. Quantitative trait locus analysis also indicates that tunnel length segregates as a complex trait, affected by at least three independent genetic regions, whereas the presence of an escape tunnel is associated with only a single locus. Together, these results suggest that complex behaviours--in this case, a classic 'extended phenotype'--can evolve through multiple genetic changes each affecting distinct behaviour modules.

Evidence of adaptation from ancestral variation in young populations of beach mice.

To understand how organisms adapt to novel habitats, which involves both demographic and selective events, we require knowledge of the evolutionary history of populations and also selected alleles. There are still few cases in which the precise mutations (and hence, defined alleles) that contribute to adaptive change have been identified in nature; one exception is the genetic basis of camouflaging pigmentation of oldfield mice (Peromyscus polionotus) that have colonized the sandy dunes of Florida's Gulf Coast. To quantify the genomic impact of colonization as well as the signature of selection, we resequenced 5000 1.5-kb noncoding loci as well as a 160-kb genomic region surrounding the melanocortin-1 receptor (Mc1r), a gene that contributes to pigmentation differences, in beach and mainland populations. Using a genome-wide phylogenetic approach, we recovered a single monophyletic group comprised of beach mice, consistent with a single colonization event of the Gulf Coast. We also found evidence of a severe founder event, estimated to have occurred less than 3000 years ago. In this demographic context, we show that all beach subspecies share a single derived light Mc1r allele, which was likely selected from standing genetic variation that originated in the mainland. Surprisingly, we were unable to identify a clear signature of selection in the Mc1r region, despite independent evidence that this locus contributes to adaptive coloration. Nonetheless, these data allow us to reconstruct and compare the evolutionary history of populations and alleles to better understand how adaptive evolution, following the colonization of a novel habitat, proceeds in nature.

Double digest RADseq: an inexpensive method for de novo SNP discovery and genotyping in model and non-model species.

The ability to efficiently and accurately determine genotypes is a keystone technology in modern genetics, crucial to studies ranging from clinical diagnostics, to genotype-phenotype association, to reconstruction of ancestry and the detection of selection. To date, high capacity, low cost genotyping has been largely achieved via "SNP chip" microarray-based platforms which require substantial prior knowledge of both genome sequence and variability, and once designed are suitable only for those targeted variable nucleotide sites. This method introduces substantial ascertainment bias and inherently precludes detection of rare or population-specific variants, a major source of information for both population history and genotype-phenotype association. Recent developments in reduced-representation genome sequencing experiments on massively parallel sequencers (commonly referred to as RAD-tag or RADseq) have brought direct sequencing to the problem of population genotyping, but increased cost and procedural and analytical complexity have limited their widespread adoption. Here, we describe a complete laboratory protocol, including a custom combinatorial indexing method, and accompanying software tools to facilitate genotyping across large numbers (hundreds or more) of individuals for a range of markers (hundreds to hundreds of thousands). Our method requires no prior genomic knowledge and achieves per-site and per-individual costs below that of current SNP chip technology, while requiring similar hands-on time investment, comparable amounts of input DNA, and downstream analysis times on the order of hours. Finally, we provide empirical results from the application of this method to both genotyping in a laboratory cross and in wild populations. Because of its flexibility, this modified RADseq approach promises to be applicable to a diversity of biological questions in a wide range of organisms.

Big genomes facilitate the comparative identification of regulatory elements.

The identification of regulatory sequences in animal genomes remains a significant challenge. Comparative genomic methods that use patterns of evolutionary conservation to identify non-coding sequences with regulatory function have yielded many new vertebrate enhancers. However, these methods have not contributed significantly to the identification of regulatory sequences in sequenced invertebrate taxa. We demonstrate here that this differential success, which is often attributed to fundamental differences in the nature of vertebrate and invertebrate regulatory sequences, is instead primarily a product of the relatively small size of sequenced invertebrate genomes. We sequenced and compared loci involved in early embryonic patterning from four species of true fruit flies (family Tephritidae) that have genomes four to six times larger than those of Drosophila melanogaster. Unlike in Drosophila, where virtually all non-coding DNA is highly conserved, blocks of conserved non-coding sequence in tephritids are flanked by large stretches of poorly conserved sequence, similar to what is observed in vertebrate genomes. We tested the activities of nine conserved non-coding sequences flanking the even-skipped gene of the teprhitid Ceratis capitata in transgenic D. melanogaster embryos, six of which drove patterns that recapitulate those of known D. melanogaster enhancers. In contrast, none of the three non-conserved tephritid non-coding sequences that we tested drove expression in D. melanogaster embryos. Based on the landscape of non-coding conservation in tephritids, and our initial success in using conservation in tephritids to identify D. melanogaster regulatory sequences, we suggest that comparison of tephritid genomes may provide a systematic means to annotate the non-coding portion of the D. melanogaster genome. We also propose that large genomes be given more consideration in the selection of species for comparative genomics projects, to provide increased power to detect functional non-coding DNAs and to provide a less biased view of the evolution and function of animal genomes.

Sepsid even-skipped enhancers are functionally conserved in Drosophila despite lack of sequence conservation.

The gene expression pattern specified by an animal regulatory sequence is generally viewed as arising from the particular arrangement of transcription factor binding sites it contains. However, we demonstrate here that regulatory sequences whose binding sites have been almost completely rearranged can still produce identical outputs. We sequenced the even-skipped locus from six species of scavenger flies (Sepsidae) that are highly diverged from the model species Drosophila melanogaster, but share its basic patterns of developmental gene expression. Although there is little sequence similarity between the sepsid eve enhancers and their well-characterized D. melanogaster counterparts, the sepsid and Drosophila enhancers drive nearly identical expression patterns in transgenic D. melanogaster embryos. We conclude that the molecular machinery that connects regulatory sequences to the transcription apparatus is more flexible than previously appreciated. In exploring this diverse collection of sequences to identify the shared features that account for their similar functions, we found a small number of short (20-30 bp) sequences nearly perfectly conserved among the species. These highly conserved sequences are strongly enriched for pairs of overlapping or adjacent binding sites. Together, these observations suggest that the local arrangement of binding sites relative to each other is more important than their overall arrangement into larger units of cis-regulatory function.

Opposed growth factor signals control protein degradation in muscles of Caenorhabditis elegans.

In addition to contractile function, muscle provides a metabolic buffer by degrading protein in times of organismal need. Protein is also degraded during adaptive muscle remodeling upon exercise, but extreme degradation in diverse clinical conditions can compromise function(s) and threaten life. Here, we show how two independent signals interact to control protein degradation. In striated muscles of Caenorhabditis elegans, reduction of insulin-like signaling via DAF-2 insulin/IGF receptor or its intramuscular effector PtdIns-3-kinase (PI3K) causes unexpected activation of MAP kinase (MAPK), consequent activation of pre-existing systems for protein degradation, and progressive impairment of mobility. Degradation is prevented by mutations that increase signal downstream of PI3K or by disruption of autocrine signal from fibroblast growth factor (FGF) via the FGF receptor and its effectors in the Ras-MAPK pathway. Thus, the activity of constitutive protein degradation systems in normal muscle is minimized by a balance between directly interacting signaling pathways, implying that physiological, pathological, or therapeutic alteration of this balance may contribute to muscle remodeling or wasting.

Clustered DNA motifs mark X chromosomes for repression by a dosage compensation complex.

Gene expression in metazoans is regulated not only at the level of individual genes but also in a coordinated manner across large chromosomal domains (for example centromeres, telomeres and imprinted gene clusters) and along entire chromosomes (for example X-chromosome dosage compensation). The primary DNA sequence usually specifies the regulation of individual genes, but the nature of cis-acting information that controls genes over large regions has been elusive: higher-order DNA structure, specific histone modifications, subnuclear compartmentalization and primary DNA sequence are possibilities. One paradigm of chromosome-wide gene regulation is Caenorhabditis elegans dosage compensation in which a large dosage compensation complex (DCC) is targeted to both X chromosomes of hermaphrodites to repress transcript levels by half. This essential process equalizes X-linked gene expression between the sexes (XO males and XX hermaphrodites). Here we report the discovery and dissection of cis-acting sites that mark nematode X chromosomes as targets for gene repression by the DCC. These rex (recruitment element on X) sites are widely dispersed along X and reside in promoters, exons and intergenic regions. rex sites share at least two distinct motifs that act in combination to recruit the DCC. Mutating these motifs severely reduces or abolishes DCC binding in vivo, demonstrating the importance of primary DNA sequence in chromosome-wide regulation. Unexpectedly, the motifs are not enriched on X, but altering motif numbers within rex sites demonstrates that motif co-occurrence in unusually high densities is essential for optimal DCC recruitment. Thus, X-specific repression is established through sequences not specific to X. The distribution of common motifs provides the foundation for repression along an entire chromosome.

Activation of Ras and the mitogen-activated protein kinase pathway promotes protein degradation in muscle cells of Caenorhabditis elegans.

To discover and study intracellular signals that regulate proteolysis in muscle, we have employed transgenic strains of Caenorhabditis elegans that produce a soluble LacZ reporter protein limited to body-wall and vulval muscles. This reporter protein is stable in well-fed wild-type animals, but its degradation is triggered upon a shift to 25 degrees C in a strain carrying a temperature-sensitive activating mutation in the Ras oncogene homologue let-60. These mutants are not physiologically starved, inasmuch as growth rates are normal at 25 degrees C. Ras-induced degradation is not prevented by the presence of cycloheximide added at or before the temperature shift and thus uses preexisting proteolytic systems and signaling components. Furthermore, degradation is triggered when adult animals are shifted to conditions of 25 degrees C, confirming that Ras acutely promotes protein degradation in muscles whose developmental history is normal. Reduction-of-function mutations in the downstream protein kinase Raf (lin-45), MEK (mek-2), or mitogen-activated protein kinase (MAPK) (mpk-1) prevent Ras-induced protein degradation, whereas activated MPK-1 is sufficient to trigger degradation, indicating that this kinase cascade is the principal route by which Ras signaling triggers protein degradation in muscle. This pathway is activated in hypodermal cells by the LET-23 epidermal growth factor receptor homologue, but an activating mutation in let-23 does not promote proteolysis in muscle. Starvation-induced LacZ reporter degradation is unaffected by reduction-of-function mutations in Ras, Raf, MEK, or MAPK, implying that Ras activation and starvation trigger proteolysis by mechanisms that are at least partially independent. This is the first evidence that Ras-Raf-MEK-MAPK signaling activates protein degradation in differentiated muscle.