Lessons from model organisms: phenotypic robustness and missing heritability in complex disease

Understanding genetic variants in context

Over the last three decades, human genetics has gone from dissecting high-penetrance Mendelian diseases to discovering the vast and complex genetic etiology of common human diseases. In tackling this complexity, scientists have discovered the importance of numerous genetic processes - most notably functional regulatory elements - in the development and progression of these diseases. Simultaneously, scientists have increasingly used multiplex assays of variant effect to systematically phenotype the cellular consequences of millions of genetic variants. In this article, we argue that the context of genetic variants - at all scales, from other genetic variants and gene regulation to cell biology to organismal environment - are critical components of how we can employ genomics to interpret these variants, and ultimately treat these diseases. We describe approaches to extend existing experimental assays and computational approaches to examine and quantify the importance of this context, including through causal analytic approaches. Having a unified understanding of the molecular, physiological, and environmental processes governing the interpretation of genetic variants is sorely needed for the field, and this perspective argues for feasible approaches by which the combined interpretation of cellular, animal, and epidemiological data can yield that knowledge.

LTP2 hypomorphs show genotype-by-environment interaction in early seedling traits in Arabidopsis thaliana

Isogenic individuals can display seemingly stochastic phenotypic differences, limiting the accuracy of genotype-to-phenotype predictions. The extent of this phenotypic variation depends in part on genetic background, raising questions about the genes involved in controlling stochastic phenotypic variation. Focusing on early seedling traits in Arabidopsis thaliana, we found that hypomorphs of the cuticle-related gene LIPID TRANSFER PROTEIN 2 (LTP2) greatly increased variation in seedling phenotypes, including hypocotyl length, gravitropism and cuticle permeability. Many ltp2 hypocotyls were significantly shorter than wild-type hypocotyls while others resembled the wild-type. Differences in epidermal properties and gene expression between ltp2 seedlings with long and short hypocotyls suggest a loss of cuticle integrity as the primary determinant of the observed phenotypic variation. We identified environmental conditions that reveal or mask the increased variation in ltp2 hypomorphs and found that increased expression of its closest paralog LTP1 is necessary for ltp2 phenotypes. Our results illustrate how decreased expression of a single gene can generate starkly increased phenotypic variation in isogenic individuals in response to an environmental challenge.

LTP2 hypomorphs show genotype-by-environment interaction in early seedling traits in Arabidopsis thaliana

Isogenic individuals can display seemingly stochastic phenotypic differences, limiting the accuracy of genotype-to-phenotype predictions. The extent of this phenotypic variation depends in part on genetic background, raising questions about the genes involved in controlling stochastic phenotypic variation. Focusing on early seedling traits in Arabidopsis thaliana, we found that hypomorphs of the cuticle-related gene LTP2 greatly increased variation in seedling phenotypes, including hypocotyl length, gravitropism and cuticle permeability. Many ltp2 hypocotyls were significantly shorter than wild-type hypocotyls while others resembled the wild type. Differences in epidermal properties and gene expression between ltp2 seedlings with long and short hypocotyls suggest a loss of cuticle integrity as the primary determinant of the observed phenotypic variation. We identified environmental conditions that reveal or mask the increased variation in ltp2 hypomorphs, and found that increased expression of its closest paralog LTP1 is necessary for ltp2 phenotypes. Our results illustrate how decreased expression of a single gene can generate starkly increased phenotypic variation in isogenic individuals in response to an environmental challenge.

Proportionality of single nucleotide causation

Assessing the proportionality of causal relationships help us pick out the most relevant causes of an effect. In this paper, I nuance and apply the concept of proportionality from the philosophy of causation debate to the mapping of associations between variations in single nucleotides in the DNA and complex phenotypic traits, such as cancer, bipolar disorder, and happiness. I discuss under what circumstances single nucleotides, as far as these can be understood as causes, may satisfy the criterion of proportionality in an interventionist understanding of causality. I suggest that the overall relatively low stability and explanatory power of such variants may indicate that, in the causation of complex phenotypic traits, there are alternative causal levels that are more proportional than the level of nucleotides. I suggest network modules as candidates for more proportional organism-internal causes of complex phenotypes. Additionally, I address the broadness of many phenotypic traits investigated in GWAS, as well as the selection between several different proportional causes of an effect.

A natural genetic variation screen identifies insulin signaling, neuronal communication, and innate immunity as modifiers of hyperglycemia in the absence of Sirt1

Variation in the onset, progression, and severity of symptoms associated with metabolic disorders such as diabetes impairs the diagnosis and treatment of at-risk patients. Diabetes symptoms, and patient variation in these symptoms, are attributed to a combination of genetic and environmental factors, but identifying the genes and pathways that modify diabetes in humans has proven difficult. A greater understanding of genetic modifiers and the ways in which they interact with metabolic pathways could improve the ability to predict a patient's risk for severe symptoms, as well as enhance the development of individualized therapeutic approaches. In this study, we use the Drosophila Genetic Reference Panel to identify genetic variation influencing hyperglycemia associated with loss of Sirt1 function. Through analysis of individual candidate functions, physical interaction networks, and gene set enrichment analysis, we identify not only modifiers involved in canonical glucose metabolism and insulin signaling, but also genes important for neuronal signaling and the innate immune response. Furthermore, reducing the expression of several of these candidates suppressed hyperglycemia, making them potential candidate therapeutic targets. These analyses showcase the diverse processes contributing to glucose homeostasis and open up several avenues of future investigation.

Functional assessment of the "two-hit" model for neurodevelopmental defects in Drosophila and X. laevis

We previously identified a deletion on chromosome 16p12.1 that is mostly inherited and associated with multiple neurodevelopmental outcomes, where severely affected probands carried an excess of rare pathogenic variants compared to mildly affected carrier parents. We hypothesized that the 16p12.1 deletion sensitizes the genome for disease, while "second-hits" in the genetic background modulate the phenotypic trajectory. To test this model, we examined how neurodevelopmental defects conferred by knockdown of individual 16p12.1 homologs are modulated by simultaneous knockdown of homologs of "second-hit" genes in Drosophila melanogaster and Xenopus laevis. We observed that knockdown of 16p12.1 homologs affect multiple phenotypic domains, leading to delayed developmental timing, seizure susceptibility, brain alterations, abnormal dendrite and axonal morphology, and cellular proliferation defects. Compared to genes within the 16p11.2 deletion, which has higher de novo occurrence, 16p12.1 homologs were less likely to interact with each other in Drosophila models or a human brain-specific interaction network, suggesting that interactions with "second-hit" genes may confer higher impact towards neurodevelopmental phenotypes. Assessment of 212 pairwise interactions in Drosophila between 16p12.1 homologs and 76 homologs of patient-specific "second-hit" genes (such as ARID1B and CACNA1A), genes within neurodevelopmental pathways (such as PTEN and UBE3A), and transcriptomic targets (such as DSCAM and TRRAP) identified genetic interactions in 63% of the tested pairs. In 11 out of 15 families, patient-specific "second-hits" enhanced or suppressed the phenotypic effects of one or many 16p12.1 homologs in 32/96 pairwise combinations tested. In fact, homologs of SETD5 synergistically interacted with homologs of MOSMO in both Drosophila and X. laevis, leading to modified cellular and brain phenotypes, as well as axon outgrowth defects that were not observed with knockdown of either individual homolog. Our results suggest that several 16p12.1 genes sensitize the genome towards neurodevelopmental defects, and complex interactions with "second-hit" genes determine the ultimate phenotypic manifestation.

A model of developmental canalization, applied to human cranial form

Developmental mechanisms that canalize or compensate perturbations of organismal development (targeted or compensatory growth) are widely considered a prerequisite of individual health and the evolution of complex life, but little is known about the nature of these mechanisms. It is even unclear if and how a “target trajectory” of individual development is encoded in the organism’s genetic-developmental system or, instead, emerges as an epiphenomenon. Here we develop a statistical model of developmental canalization based on an extended autoregressive model. We show that under certain assumptions the strength of canalization and the amount of canalized variance in a population can be estimated, or at least approximated, from longitudinal phenotypic measurements, even if the target trajectories are unobserved. We extend this model to multivariate measures and discuss reifications of the ensuing parameter matrix. We apply these approaches to longitudinal geometric morphometric data on human postnatal craniofacial size and shape as well as to the size of the frontal sinuses. Craniofacial size showed strong developmental canalization during the first 5 years of life, leading to a 50% reduction of cross-sectional size variance, followed by a continual increase in variance during puberty. Frontal sinus size, by contrast, did not show any signs of canalization. Total variance of craniofacial shape decreased slightly until about 5 years of age and increased thereafter. However, different features of craniofacial shape showed very different developmental dynamics. Whereas the relative dimensions of the nasopharynx showed strong canalization and a reduction of variance throughout postnatal development, facial orientation continually increased in variance. Some of the signals of canalization may owe to independent variation in developmental timing of cranial components, but our results indicate evolved, partly mechanically induced mechanisms of canalization that ensure properly sized upper airways and facial dimensions.

Developmental mechanisms that canalize or compensate perturbations of organismal development are a prerequisite of individual health and the evolution of complex life. However, surprisingly little is known about these mechanisms, partly because the “target trajectories” of individual development cannot be directly observed. Here we develop a statistical model of developmental canalization that allows one to estimate the strength of canalization and the amount of canalized variance in a population even if the target trajectories are unobserved. We applied these approaches to data on human postnatal craniofacial growth. Whereas overall craniofacial size was strongly canalized during the first 5 years of age, frontal sinus size did not show any signs of canalization. The relative dimensions of the nasopharynx showed strong canalization and a reduction of variance throughout postnatal development, whereas other shape features, such as facial orientation, continually increased in variance. Our results indicate evolved, partly mechanically induced mechanisms of canalization that ensure properly sized upper airways and facial dimensions.

Genomic Variation, Evolvability, and the Paradox of Mental Illness

Twentieth-century genetics was hard put to explain the irregular behavior of neuropsychiatric disorders. Autism and schizophrenia defy a principle of natural selection; they are highly heritable but associated with low reproductive success. Nevertheless, they persist. The genetic origins of such conditions are confounded by the problem of variable expression, that is, when a given genetic aberration can lead to any one of several distinct disorders. Also, autism and schizophrenia occur on a spectrum of severity, from mild and subclinical cases to the overt and disabling. Such irregularities reflect the problem of missing heritability; although hundreds of genes may be associated with autism or schizophrenia, together they account for only a small proportion of cases. Techniques for higher resolution, genomewide analysis have begun to illuminate the irregular and unpredictable behavior of the human genome. Thus, the origins of neuropsychiatric disorders in particular and complex disease in general have been illuminated. The human genome is characterized by a high degree of structural and behavioral variability: DNA content variation, epistasis, stochasticity in gene expression, and epigenetic changes. These elements have grown more complex as evolution scaled the phylogenetic tree. They are especially pertinent to brain development and function. Genomic variability is a window on the origins of complex disease, neuropsychiatric disorders, and neurodevelopmental disorders in particular. Genomic variability, as it happens, is also the fuel of evolvability. The genomic events that presided over the evolution of the primate and hominid lineages are over-represented in patients with autism and schizophrenia, as well as intellectual disability and epilepsy. That the special qualities of the human genome that drove evolution might, in some way, contribute to neuropsychiatric disorders is a matter of no little interest.

The Malmö Offspring Study (MOS): design, methods and first results

As cardio metabolic disease manifestations tend to cluster in families there is a need to better understand the underlying mechanisms in order to further develop preventive strategies. In fact, genetic markers used in genetic risk scores, important as they are, will not be able alone to explain these family clusters. Therefore, the search goes on for the so called missing heritability to better explain these associations. Shared lifestyle and social conditions in families, but also early life influences may be of importance. Gene-environmental interactions should be explored. In recent years interest has grown for the role of diet-microbiota associations, as microbiota patterns may be shared by family members. In the Malmö Offspring Study that started in 2013, we have so far been able to examine about 4700 subjects (18-71 years) representing children and grandchildren of index subjects from the first generation, examined in the Malmö Diet Cancer Study during 1991 to 1996. This will provide rich data and opportunities to analyse family traits of chronic disease across three generations. We will provide extensive genotyping and phenotyping including cardiovascular and respiratory function, as well as markers of glucose metabolism. In addition, also cognitive function will be assessed. A 4-day online dietary recall will be conducted and gut as well as oral microbiota analysed. The ambition is to provide one of the first large-scale European family studies with individual data across three generations, which could deepen our knowledge about the role of family traits for chronic disease and its underlying mechanisms.

Chaperone biomarkers of lifespan and penetrance track the dosages of many other proteins

Many traits vary among isogenic individuals in homogeneous environments. In microbes, plants and animals, variation in the protein chaperone system affects many such traits. In the animal model C. elegans, the expression level of hsp-16.2 chaperone biomarkers correlates with or predicts the penetrance of mutations and lifespan after heat shock. But the physiological mechanisms causing cells to express different amounts of the biomarker were unknown. Here, we used an in vivo microscopy approach to dissect different contributions to cell-to-cell variation in hsp-16.2 expression in the intestines of young adult animals, which generate the most lifespan predicting signal. While we detected both cell autonomous intrinsic noise and signaling noise, we found both contributions were relatively unimportant. The major contributor to cell-to-cell variation in biomarker expression was general differences in protein dosage. The hsp-16.2 biomarker reveals states of high or low effective dosage for many genes.

Abundance of ethnically biased microsatellites in human gene regions

Microsatellites-a type of short tandem repeat (STR)-have been used for decades as putatively neutral markers to study the genetic structure of diverse human populations. However, recent studies have demonstrated that some microsatellites contribute to gene expression, cis heritability, and phenotype. As a corollary, some microsatellites may contribute to differential gene expression and RNA/protein structure stability in distinct human populations. To test this hypothesis, we investigate genotype frequencies, functional relevance, and adaptive potential of microsatellites in five super-populations (ethnicities) drawn from the 1000 Genomes Project. We discover 3,984 ethnically-biased microsatellite loci (EBML); for each EBML at least one ethnicity has genotype frequencies statistically different from the remaining four. South Asian, East Asian, European, and American EBML show significant overlap; on the contrary, the set of African EBML is mostly unique. We cross-reference the 3,984 EBML with 2,060 previously identified expression STRs (eSTRs); repeats known to affect gene expression (64 total) are over-represented. The most significant pathway enrichments are those associated with the matrisome: a broad collection of genes encoding the extracellular matrix and its associated proteins. At least 14 of the EBML have established links to human disease. Analysis of the 3,984 EBML with respect to known selective sweep regions in the genome shows that allelic variation in some of them is likely associated with adaptive evolution.

Natural Genetic Variation Screen in Drosophila Identifies Wnt Signaling, Mitochondrial Metabolism, and Redox Homeostasis Genes as Modifiers of Apoptosis

Apoptosis is the primary cause of degeneration in a number of neuronal, muscular, and metabolic disorders. These diseases are subject to a great deal of phenotypic heterogeneity in patient populations, primarily due to differences in genetic variation between individuals. This creates a barrier to effective diagnosis and treatment. Understanding how genetic variation influences apoptosis could lead to the development of new therapeutics and better personalized treatment approaches. In this study, we examine the impact of the natural genetic variation in the Drosophila Genetic Reference Panel (DGRP) on two models of apoptosis-induced retinal degeneration: overexpression of p53 or reaper (rpr). We identify a number of known apoptotic, neural, and developmental genes as candidate modifiers of degeneration. We also use Gene Set Enrichment Analysis (GSEA) to identify pathways that harbor genetic variation that impact these apoptosis models, including Wnt signaling, mitochondrial metabolism, and redox homeostasis. Finally, we demonstrate that many of these candidates have a functional effect on apoptosis and degeneration. These studies provide a number of avenues for modifying genes and pathways of apoptosis-related disease.

Epistasis regulates the developmental stability of the mouse craniofacial shape

Fluctuating asymmetry is a classic concept linked to organismal development. It has traditionally been used as a measure of developmental instability, which is the inability of an organism to buffer environmental fluctuations during development. Developmental stability has a genetic component that influences the final phenotype of the organism and can lead to congenital disorders. According to alternative hypotheses, this genetic component might be either the result of additive genetic effects or a by-product of developmental gene networks. Here we present a genome-wide association study of the genetic architecture of fluctuating asymmetry of the skull shape in mice. Geometric morphometric methods were applied to quantify fluctuating asymmetry: we estimated fluctuating asymmetry as Mahalanobis distances to the mean asymmetry, correcting first for genetic directional asymmetry. We applied the marginal epistasis test to study epistasis among genomic regions. Results showed no evidence of additive effects but several interacting regions significantly associated with fluctuating asymmetry. Among the candidate genes overlapping these interacting regions we found an over-representation of genes involved in craniofacial development. A gene network is likely to be associated with skull developmental stability, and genes originally described as buffering genes (e.g., Hspa2) might occupy central positions within these networks, where regulatory elements may also play an important role. Our results constitute an important step in the exploration of the molecular roots of developmental stability and the first empirical evidence about its genetic architecture.

Translesion Synthesis in Plants: Ultraviolet Resistance and Beyond

Plant genomes sustain various forms of DNA damage that stall replication forks. Translesion synthesis (TLS) is one of the pathways to overcome stalled replication in which specific polymerases (TLS polymerase) perform bypass synthesis across DNA damage. This article gives a brief overview of plant TLS polymerases. In Arabidopsis, DNA polymerase (Pol) ζ, η, κ, θ, and λ and Reversionless1 (Rev1) are shown to be involved in the TLS. For example, AtPolη bypasses ultraviolet (UV)-induced cyclobutane pyrimidine dimers in vitro. Disruption of AtPolζ or AtPolη increases root stem cell death after UV irradiation. These results suggest that AtPolζ and ATPolη bypass UV-induced damage, prevent replication arrest, and allow damaged cells to survive and grow. In general, TLS polymerases have low fidelity and often induce mutations. Accordingly, disruption of AtPolζ or AtRev1 reduces somatic mutation frequency, whereas disruption of AtPolη elevates it, suggesting that plants have both mutagenic and less mutagenic TLS activities. The stalled replication fork can be resolved by a strand switch pathway involving a DNA helicase Rad5. Disruption of both AtPolζ and AtRAD5a shows synergistic or additive effects in the sensitivity to DNA-damaging agents. Moreover, AtPolζ or AtRev1 disruption elevates homologous recombination frequencies in somatic tissues. These results suggest that the Rad5-dependent pathway and TLS are parallel. Plants grown in the presence of heat shock protein 90 (HSP90) inhibitor showed lower mutation frequencies, suggesting that HSP90 regulates mutagenic TLS in plants. Hypersensitivities of TLS-deficient plants to γ-ray and/or crosslink damage suggest that plant TLS polymerases have multiple roles, as reported in other organisms.

Mutations, protein homeostasis, and epigenetic control of genome integrity

From bacteria to humans, ancient stress responses enable organisms to contend with damage to both the genome and the proteome. These pathways have long been viewed as fundamentally separate responses. Yet recent discoveries from multiple fields have revealed surprising links between the two. Many DNA-damaging agents also target proteins, and mutagenesis induced by DNA damage produces variant proteins that are prone to misfolding, degradation, and aggregation. Likewise, recent studies have observed pervasive engagement of a p53-mediated response, and other factors linked to maintenance of genomic integrity, in response to misfolded protein stress. Perhaps most remarkably, protein aggregation and self-assembly has now been observed in multiple proteins that regulate the DNA damage response. The importance of these connections is highlighted by disease models of both cancer and neurodegeneration, in which compromised DNA repair machinery leads to profound defects in protein quality control, and vice versa.

Baldspot/ELOVL6 is a conserved modifier of disease and the ER stress response

Endoplasmic reticulum (ER) stress is an important modifier of human disease. Genetic variation in response genes is linked to inter-individual differences in the ER stress response. However, the mechanisms and pathways by which genetic modifiers are acting on the ER stress response remain unclear. In this study, we characterize the role of the long chain fatty acid elongase Baldspot (ELOVL6) in modifying the ER stress response and disease. We demonstrate that loss of Baldspot rescues degeneration and reduces IRE1 and PERK signaling and cell death in a Drosophila model of retinitis pigmentosa and ER stress (Rh1^G69D). Dietary supplementation of stearate bypasses the need for Baldspot activity. Finally, we demonstrate that Baldspot regulates the ER stress response across different tissues and induction methods. Our findings suggest that ELOVL6 is a promising target in the treatment of not only retinitis pigmentosa, but a number of different ER stress-related disorders.

Differences in genetic background drives disease variability, even among individuals with identical, causative mutations. Identifying and understanding how genetic variation impacts disease expression could improve diagnosis and treatment of patients. Previous work has linked the endoplasmic reticulum (ER) stress response pathway to disease variability. When misfolded proteins accumulate in the ER, the ER stress response returns the cell to its normal state. Chronic ER stress leads to massive amounts of cell death and tissue degeneration. Limiting tissue loss by regulating the ER stress response has been a major focus of therapeutic development. In this study, we characterize a novel regulator of the ER stress response, the long chain fatty acid elongase Baldspot/ELOVL6. In the absence of this enzyme, cells undergoing ER stress display reduced cell death, and degeneration in a Drosophila disease model. Feeding of excess fatty acids increases degeneration to original disease levels, linking the regulatory activity of Baldspot to its enzymatic activity. Finally, we demonstrate that Baldspot can alter the ER stress response under a variety of other ER stress conditions. Our studies demonstrate that Baldspot/ELOVL6 is a ubiquitous regulator of the ER stress response and is a good candidate therapeutic target.

It's not magic - Hsp90 and its effects on genetic and epigenetic variation

Canalization, or phenotypic robustness in the face of environmental and genetic perturbation, is an emergent property of living systems. Although this phenomenon has long been recognized, its molecular underpinnings have remained enigmatic until recently. Here, we review the contributions of the molecular chaperone Hsp90, a protein that facilitates the folding of many key regulators of growth and development, to canalization of phenotype - and de-canalization in times of stress - drawing on studies in eukaryotes as diverse as baker's yeast, mouse ear cress, and blind Mexican cavefish. Hsp90 is a hub of hubs that interacts with many so-called 'client proteins,' which affect virtually every aspect of cell signaling and physiology. As Hsp90 facilitates client folding and stability, it can epistatically suppress or enable the expression of genetic variants in its clients and other proteins that acquire client status through mutation. Hsp90's vast interaction network explains the breadth of its phenotypic reach, including Hsp90-dependent de novo mutations and epigenetic effects on gene regulation. Intrinsic links between environmental stress and Hsp90 function thus endow living systems with phenotypic plasticity in fluctuating environments. As environmental perturbations alter Hsp90 function, they also alter Hsp90's interaction with its client proteins, thereby re-wiring networks that determine the genotype-to-phenotype map. Ensuing de-canalization of phenotype creates phenotypic diversity that is not simply stochastic, but often has an underlying genetic basis. Thus, extreme phenotypes can be selected, and assimilated so that they no longer require environmental stress to manifest. In addition to acting on standing genetic variation, Hsp90 perturbation has also been linked to increased frequency of de novo variation and several epigenetic phenomena, all with the potential to generate heritable phenotypic change. Here, we aim to clarify and discuss the multiple means by which Hsp90 can affect phenotype and possibly evolutionary change, and identify their underlying common feature: at its core, Hsp90 interacts epistatically through its chaperone function with many other genes and their gene products. Its influence on phenotypic diversification is thus not magic but rather a fundamental property of genetics.

Another Round of "Clue" to Uncover the Mystery of Complex Traits

A plethora of genetic association analyses have identified several genetic risk loci. Technological and statistical advancements have now led to the identification of not only common genetic variants, but also low-frequency variants, structural variants, and environmental factors, as well as multi-omics variations that affect the phenotypic variance of complex traits in a population, thus referred to as complex trait architecture. The concept of heritability, or the proportion of phenotypic variance due to genetic inheritance, has been studied for several decades, but its application is mainly in addressing the narrow sense heritability (or additive genetic component) from Genome-Wide Association Studies (GWAS). In this commentary, we reflect on our perspective on the complexity of understanding heritability for human traits in comparison to model organisms, highlighting another round of clues beyond GWAS and an alternative approach, investigating these clues comprehensively to help in elucidating the genetic architecture of complex traits.

Systematic Profiling of Short Tandem Repeats in the Cattle Genome

Short tandem repeats (STRs), or microsatellites, are genetic variants with repetitive 2–6 base pair motifs in many mammalian genomes. Using high-throughput sequencing and experimental validations, we systematically profiled STRs in five Holsteins. We identified a total of 60,106 microsatellites and generated the first high-resolution STR map, representing a substantial pool of polymorphism in dairy cattle. We observed significant STRs overlap with functional genes and quantitative trait loci (QTL). We performed evolutionary and population genetic analyses using over 20,000 common dinucleotide STRs. Besides corroborating the well-established positive correlation between allele size and variance in allele size, these analyses also identified dozens of outlier STRs based on two anomalous relationships that counter expected characteristics of neutral evolution. And one STR locus overlaps with a significant region of a summary statistic designed to detect STR-related selection. Additionally, our results showed that only 57.1% of STRs located within SNP-based linkage disequilibrium (LD) blocks whereas the other 42.9% were out of blocks. Therefore, a substantial number of STRs are not tagged by SNPs in the cattle genome, likely due to STR's distinct mutation mechanism and elevated polymorphism. This study provides the foundation for future STR-based studies of cattle genome evolution and selection.

Deletion of TLS polymerases promotes homologous recombination in Arabidopsis

Unrepaired DNA damage hinders the maintenance of genome integrity because it blocks the catalytic activity of replicase. The stalled replication fork can be processed through either translesion synthesis (TLS) with specific polymerases, or replication using the undamaged template. To investigate how TLS activities are regulated and how the stalled replication fork is processed in plants, reversion frequencies and homologous recombination (HR) frequencies were analyzed using GUS-based substrates. The HR frequencies in plants deficient in DNA polymerase ζ (Pol ζ) or Rev1 were higher than that in wildtype plants under normal conditions, and were significantly increased by ultraviolet light irradiation. Heat shock protein (HSP) 90 is known to be involved in various stress responses. To examine the role of HSP90 in the regulation of damage tolerance, we analyzed reversion frequencies and HR frequencies in plants grown in the presence of a HSP inhibitor, geldanamycin (GDA). Reversion frequency was lower in GDA-treated plants than in mock-treated plants. Though the HR frequency was higher in GDA-treated wildtype plants than in mock-treated plants, no significant difference was detected in Rev1-deficient plants. In yeast, TLS polymerases interacted with each other or with a replication clump component, proliferating cell nuclear antigen (PCNA). HSP90 interacted with REV1 or REV7 in Nicotiana benthamiana cells. These results suggest that HSP90 interacts with TLS polymerase(s), which promotes error-prone TLS in plants.

Dissecting genetic architecture of startle response in Drosophila melanogaster using multi-omics information

Startle behavior is important for survival, and abnormal startle responses are related to several neurological diseases. Drosophila melanogaster provides a powerful system to investigate the genetic underpinnings of variation in startle behavior. Since mechanically induced, startle responses and environmental conditions can be readily quantified and precisely controlled. The 156 wild-derived fully sequenced lines of the Drosophila Genetic Reference Panel (DGRP) were used to identify SNPs and transcripts associated with variation in startle behavior. The results validated highly significant effects of 33 quantitative trait SNPs (QTSs) and 81 quantitative trait transcripts (QTTs) directly associated with phenotypic variation of startle response. We also detected QTT variation controlled by 20 QTSs (tQTSs) and 73 transcripts (tQTTs). Association mapping based on genomic and transcriptomic data enabled us to construct a complex genetic network that underlies variation in startle behavior. Based on principles of evolutionary conservation, human orthologous genes could be superimposed on this network. This study provided both genetic and biological insights into the variation of startle response behavior of Drosophila melanogaster, and highlighted the importance of genetic network to understand the genetic architecture of complex traits.

MicroRNAs Establish Uniform Traits during the Architecture of Vertebrate Embryos

Proper functioning of an organism requires cells and tissues to behave in uniform, well-organized ways. How this optimum of phenotypes is achieved during the development of vertebrates is unclear. Here, we carried out a multi-faceted and single-cell resolution screen of zebrafish embryonic blood vessels upon mutagenesis of single and multi-gene microRNA (miRNA) families. We found that embryos lacking particular miRNA-dependent signaling pathways develop a vascular trait similar to wild-type, but with a profound increase in phenotypic heterogeneity. Aberrant trait variance in miRNA mutant embryos uniquely sensitizes their vascular system to environmental perturbations. We discovered a previously unrecognized role for specific vertebrate miRNAs to protect tissue development against phenotypic variability. This discovery marks an important advance in our comprehension of how miRNAs function in the development of higher organisms.

Stress-Induced Mutagenesis: Implications in Cancer and Drug Resistance

Genomic instability underlies many cancers and generates genetic variation that drives cancer initiation, progression, and therapy resistance. In contrast with classical assumptions that mutations occur purely stochastically at constant, gradual rates, microbes, plants, flies, and human cancer cells possess mechanisms of mutagenesis that are upregulated by stress responses. These generate transient, genetic-diversity bursts that can propel evolution, specifically when cells are poorly adapted to their environments-that is, when stressed. We review molecular mechanisms of stress-response-dependent (stress-induced) mutagenesis that occur from bacteria to cancer, and are activated by starvation, drugs, hypoxia, and other stressors. We discuss mutagenic DNA break repair in Escherichia coli as a model for mechanisms in cancers. The temporal regulation of mutagenesis by stress responses and spatial restriction in genomes are common themes across the tree of life. Both can accelerate evolution, including the evolution of cancers. We discuss possible anti-evolvability drugs, aimed at targeting mutagenesis and other variation generators, that could be used to delay the evolution of cancer progression and therapy resistance.

Variability in a Short Tandem Repeat Mediates Complex Epistatic Interactions in Arabidopsis thaliana

Short tandem repeats (STRs) are hypervariable genetic elements that occur frequently in coding regions. Their high mutation rate readily generates genetic variation, contributing to adaptive evolution and human diseases. We previously reported that natural ELF3 polyglutamine variants cause reciprocal genetic incompatibilities in two divergent Arabidopsis thaliana backgrounds. Here, we dissect the genetic architecture of this incompatibility, revealing as many as four loci putatively interacting with ELF3 We were able to specifically identify one such ELF3-interacting gene, LSH9 We further used a yeast two-hybrid strategy to identify proteins whose physical interactions with ELF3 were affected by polyglutamine tract length. We found two proteins for which this was the case, ELF4 and AtGLDP1. Using these two approaches, we identify specific genetic interactions and physical mechanisms by which the ELF3 polyglutamine tract may mediate the observed genetic incompatibilities. Our work elucidates how STR variation, which is generally underascertained in population-scale sequencing, can contribute to phenotypic variation. Furthermore, our results support our proposal that highly variable STR loci can contribute to the epistatic component of heritability.

Personalized learning: From neurogenetics of behaviors to designing optimal language training

Variability in drug responsivity has prompted the development of Personalized Medicine, which has shown great promise in utilizing genotypic information to develop safer and more effective drug regimens for patients. Similarly, individual variability in learning outcomes has puzzled researchers who seek to create optimal learning environments for students. "Personalized Learning" seeks to identify genetic, neural and behavioral predictors of individual differences in learning and aims to use predictors to help create optimal teaching paradigms. Evidence for Personalized Learning can be observed by connecting research in pharmacogenomics, cognitive genetics and behavioral experiments across domains of learning, which provides a framework for conducting empirical studies from the laboratory to the classroom and holds promise for addressing learning effectiveness in the individual learners. Evidence can also be seen in the subdomain of speech learning, thus providing initial support for the applicability of Personalized Learning to language.

Network biology concepts in complex disease comorbidities

The co-occurrence of diseases can inform the underlying network biology of shared and multifunctional genes and pathways. In addition, comorbidities help to elucidate the effects of external exposures, such as diet, lifestyle and patient care. With worldwide health transaction data now often being collected electronically, disease co-occurrences are starting to be quantitatively characterized. Linking network dynamics to the real-life, non-ideal patient in whom diseases co-occur and interact provides a valuable basis for generating hypotheses on molecular disease mechanisms, and provides knowledge that can facilitate drug repurposing and the development of targeted therapeutic strategies.

Genetic structure of phenotypic robustness in the collaborative cross mouse diallel panel

Developmental stability and canalization describe the ability of developmental systems to minimize phenotypic variation in the face of stochastic micro-environmental effects, genetic variation and environmental influences. Canalization is the ability to minimize the effects of genetic or environmental effects, whereas developmental stability is the ability to minimize the effects of micro-environmental effects within individuals. Despite much attention, the mechanisms that underlie these two components of phenotypic robustness remain unknown. We investigated the genetic structure of phenotypic robustness in the collaborative cross (CC) mouse reference population. We analysed the magnitude of fluctuating asymmetry (FA) and among-individual variation of cranial shape in reciprocal crosses among the eight parental strains, using geometric morphometrics and a diallel analysis based on a Bayesian approach. Significant differences among genotypes were found for both measures, although they were poorly correlated at the level of individuals. An overall positive effect of inbreeding was found for both components of variation. The strain CAST/EiJ exerted a positive additive effect on FA and, to a lesser extent, among-individual variance. Sex- and other strain-specific effects were not significant. Neither FA nor among-individual variation was associated with phenotypic extremeness. Our results support the existence of genetic variation for both developmental stability and canalization. This finding is important because robustness is a key feature of developmental systems. Our finding that robustness is not related to phenotypic extremeness is consistent with theoretical work that suggests that its relationship to stabilizing selection is not straightforward.

Cooperativity of Negative Autoregulation Confers Increased Mutational Robustness

Negative autoregulation is universally found across organisms. In the bacterium Escherichia coli, transcription factors often repress their own expression to form a negative feedback network motif that enables robustness to changes in biochemical parameters. Here we present a simple phenomenological model of a negative feedback transcription factor repressing both itself and another target gene. The strength of the negative feedback is characterized by three parameters: the cooperativity in self-repression, the maximal expression rate of the transcription factor, and the apparent dissociation constant of the transcription factor binding to its own promoter. Analysis of the model shows that the target gene levels are robust to mutations in the transcription factor, and that the robustness improves as the degree of cooperativity in self-repression increases. The prediction is tested in the LexA transcriptional network of E. coli by altering cooperativity in self-repression and promoter strength. Indeed, we find robustness is correlated with the former. Considering the proposed importance of gene regulation in speciation, parameters governing a transcription factor's robustness to mutation may have significant influence on a cell or organism's capacity to evolve.

Genetic and Environmental Control of Neurodevelopmental Robustness in Drosophila

Interindividual differences in neuronal wiring may contribute to behavioral individuality and affect susceptibility to neurological disorders. To investigate the causes and potential consequences of wiring variation in Drosophila melanogaster, we focused on a hemilineage of ventral nerve cord interneurons that exhibits morphological variability. We find that late-born subclasses of the 12A hemilineage are highly sensitive to genetic and environmental variation. Neurons in the second thoracic segment are particularly variable with regard to two developmental decisions, whereas its segmental homologs are more robust. This variability “hotspot” depends on Ultrabithorax expression in the 12A neurons, indicating variability is cell-intrinsic and under genetic control. 12A development is more variable and sensitive to temperature in long-established laboratory strains than in strains recently derived from the wild. Strains with a high frequency of one of the 12A variants also showed a high frequency of animals with delayed spontaneous flight initiation, whereas other wing-related behaviors did not show such a correlation and were thus not overtly affected by 12A variation. These results show that neurodevelopmental robustness is variable and under genetic control in Drosophila and suggest that the fly may serve as a model for identifying conserved gene pathways that stabilize wiring in stressful developmental environments. Moreover, some neuronal lineages are variation hotspots and thus may be more amenable to evolutionary change.

Genetic assimilation: a review of its potential proximate causes and evolutionary consequences

BACKGROUND

Most, if not all, organisms possess the ability to alter their phenotype in direct response to changes in their environment, a phenomenon known as phenotypic plasticity. Selection can break this environmental sensitivity, however, and cause a formerly environmentally induced trait to evolve to become fixed through a process called genetic assimilation. Essentially, genetic assimilation can be viewed as the evolution of environmental robustness in what was formerly an environmentally sensitive trait. Because genetic assimilation has long been suggested to play a key role in the origins of phenotypic novelty and possibly even new species, identifying and characterizing the proximate mechanisms that underlie genetic assimilation may advance our basic understanding of how novel traits and species evolve.

SCOPE

This review begins by discussing how the evolution of phenotypic plasticity, followed by genetic assimilation, might promote the origins of new traits and possibly fuel speciation and adaptive radiation. The evidence implicating genetic assimilation in evolutionary innovation and diversification is then briefly considered. Next, the potential causes of phenotypic plasticity generally and genetic assimilation specifically are examined at the genetic, molecular and physiological levels and approaches that can improve our understanding of these mechanisms are described. The review concludes by outlining major challenges for future work.

CONCLUSIONS

Identifying and characterizing the proximate mechanisms involved in phenotypic plasticity and genetic assimilation promises to help advance our basic understanding of evolutionary innovation and diversification.

The Mechanistic Underpinnings of an ago1-Mediated, Environmentally Dependent, and Stochastic Phenotype

The crucial role of microRNAs in plant development is exceedingly well supported; their importance in environmental robustness is studied in less detail. Here, we describe a novel, environmentally dependent phenotype in hypomorphic argonaute1 (ago1) mutants and uncover its mechanistic underpinnings in Arabidopsis (Arabidopsis thaliana). AGO1 is a key player in microRNA-mediated gene regulation. We observed transparent lesions on embryonic leaves of ago1 mutant seedlings. These lesions increased in frequency in full-spectrum light. Notably, the lesion phenotype was most environmentally responsive in ago1-27 mutants. This allele is thought to primarily affect translational repression, which has been linked with the response to environmental perturbation. Using several lines of evidence, we found that these lesions represent dead and dying tissues due to an aberrant hypersensitive response. Although all three canonical defense hormone pathways (salicylic acid, jasmonate, and jasmonate/ethylene pathways) were up-regulated in ago1 mutants, we demonstrate that jasmonate perception drives the lesion phenotype. Double mutants of ago1 and coronatine insensitive1, the jasmonate receptor, showed greatly decreased frequency of affected seedlings. The chaperone HEAT SHOCK PROTEIN 90 (HSP90), which maintains phenotypic robustness in the face of environmental perturbations, is known to facilitate AGO1 function. HSP90 perturbation has been shown previously to up-regulate jasmonate signaling and to increase plant resistance to herbivory. Although single HSP90 mutants showed subtly elevated levels of lesions, double mutant analysis disagreed with a simple epistatic model for HSP90 and AGO1 interaction; rather, both appeared to act nonadditively in producing lesions. In summary, our study identifies AGO1 as a major, largely HSP90-independent, factor in providing environmental robustness to plants.

A Genome-Wide Association Analysis Reveals Epistatic Cancellation of Additive Genetic Variance for Root Length in Arabidopsis thaliana

Efforts to identify loci underlying complex traits generally assume that most genetic variance is additive. Here, we examined the genetics of Arabidopsis thaliana root length and found that the genomic narrow-sense heritability for this trait in the examined population was statistically zero. The low amount of additive genetic variance that could be captured by the genome-wide genotypes likely explains why no associations to root length could be found using standard additive-model-based genome-wide association (GWA) approaches. However, as the broad-sense heritability for root length was significantly larger, and primarily due to epistasis, we also performed an epistatic GWA analysis to map loci contributing to the epistatic genetic variance. Four interacting pairs of loci were revealed, involving seven chromosomal loci that passed a standard multiple-testing corrected significance threshold. The genotype-phenotype maps for these pairs revealed epistasis that cancelled out the additive genetic variance, explaining why these loci were not detected in the additive GWA analysis. Small population sizes, such as in our experiment, increase the risk of identifying false epistatic interactions due to testing for associations with very large numbers of multi-marker genotypes in few phenotyped individuals. Therefore, we estimated the false-positive risk using a new statistical approach that suggested half of the associated pairs to be true positive associations. Our experimental evaluation of candidate genes within the seven associated loci suggests that this estimate is conservative; we identified functional candidate genes that affected root development in four loci that were part of three of the pairs. The statistical epistatic analyses were thus indispensable for confirming known, and identifying new, candidate genes for root length in this population of wild-collected A. thaliana accessions. We also illustrate how epistatic cancellation of the additive genetic variance explains the insignificant narrow-sense and significant broad-sense heritability by using a combination of careful statistical epistatic analyses and functional genetic experiments.

Complex traits, such as many human diseases or climate adaptation and production traits in crops, arise through the action and interaction of many genes and environmental factors. Classic approaches to identify contributing genes generally assume that these factors contribute mainly additive genetic variance. Recent methods, such as genome-wide association studies, often adhere to this additive genetics paradigm. However, additive models of complex traits do not reflect that genes can also contribute with non-additive genetic variance. In this study, we use Arabidopsis thaliana to determine the additive and non-additive genetic contributions to the phenotypic variation in root length. Surprisingly, much of the observed phenotypic variation in root length across genetically divergent strains was explained by epistasis. We mapped seven loci contributing to the epistatic genetic variance and validated four genes in these loci with mutant analysis. For three of these genes, this is their first implication in root development. Together, our results emphasize the importance of considering both non-additive and additive genetic variance when dissecting complex trait variation, in order not to lose sensitivity in genetic analyses.

Genetic Architecture of Natural Variation in Thermal Responses of Arabidopsis

Wild strains of Arabidopsis (Arabidopsis thaliana) exhibit extensive natural variation in a wide variety of traits, including response to environmental changes. Ambient temperature is one of the major external factors that modulates plant growth and development. Here, we analyze the genetic architecture of natural variation in thermal responses of Arabidopsis. Exploiting wild accessions and recombinant inbred lines, we reveal extensive phenotypic variation in response to ambient temperature in distinct developmental traits such as hypocotyl elongation, root elongation, and flowering time. We show that variation in thermal response differs between traits, suggesting that the individual phenotypes do not capture all the variation associated with thermal response. Genome-wide association studies and quantitative trait locus analyses reveal that multiple rare alleles contribute to the genetic architecture of variation in thermal response. We identify at least 20 genomic regions that are associated with variation in thermal response. Further characterizations of temperature sensitivity quantitative trait loci that are shared between traits reveal a role for the blue-light receptor CRYPTOCHROME2 (CRY2) in thermosensory growth responses. We show the accession Cape Verde Islands is less sensitive to changes in ambient temperature, and through transgenic analysis, we demonstrate that allelic variation at CRY2 underlies this temperature insensitivity across several traits. Transgenic analyses suggest that the allelic effects of CRY2 on thermal response are dependent on genetic background suggestive of the presence of modifiers. In addition, our results indicate that complex light and temperature interactions, in a background-dependent manner, govern growth responses in Arabidopsis.

Network Modules of the Cross-Species Genotype-Phenotype Map Reflect the Clinical Severity of Human Diseases

Recent advances in genome sequencing techniques have improved our understanding of the genotype-phenotype relationship between genetic variants and human diseases. However, genetic variations uncovered from patient populations do not provide enough information to understand the mechanisms underlying the progression and clinical severity of human diseases. Moreover, building a high-resolution genotype-phenotype map is difficult due to the diverse genetic backgrounds of the human population. We built a cross-species genotype-phenotype map to explain the clinical severity of human genetic diseases. We developed a data-integrative framework to investigate network modules composed of human diseases mapped with gene essentiality measured from a model organism. Essential and nonessential genes connect diseases of different types which form clusters in the human disease network. In a large patient population study, we found that disease classes enriched with essential genes tended to show a higher mortality rate than disease classes enriched with nonessential genes. Moreover, high disease mortality rates are explained by the multiple comorbid relationships and the high pleiotropy of disease genes found in the essential gene-enriched diseases. Our results reveal that the genotype-phenotype map of a model organism can facilitate the identification of human disease-gene associations and predict human disease progression.

MIPSTR: a method for multiplex genotyping of germline and somatic STR variation across many individuals

Short tandem repeats (STRs) are highly mutable genetic elements that often reside in regulatory and coding DNA. The cumulative evidence of genetic studies on individual STRs suggests that STR variation profoundly affects phenotype and contributes to trait heritability. Despite recent advances in sequencing technology, STR variation has remained largely inaccessible across many individuals compared to single nucleotide variation or copy number variation. STR genotyping with short-read sequence data is confounded by (1) the difficulty of uniquely mapping short, low-complexity reads; and (2) the high rate of STR amplification stutter. Here, we present MIPSTR, a robust, scalable, and affordable method that addresses these challenges. MIPSTR uses targeted capture of STR loci by single-molecule Molecular Inversion Probes (smMIPs) and a unique mapping strategy. Targeted capture and our mapping strategy resolve the first challenge; the use of single molecule information resolves the second challenge. Unlike previous methods, MIPSTR is capable of distinguishing technical error due to amplification stutter from somatic STR mutations. In proof-of-principle experiments, we use MIPSTR to determine germline STR genotypes for 102 STR loci with high accuracy across diverse populations of the plant A. thaliana. We show that putatively functional STRs may be identified by deviation from predicted STR variation and by association with quantitative phenotypes. Using DNA mixing experiments and a mutant deficient in DNA repair, we demonstrate that MIPSTR can detect low-frequency somatic STR variants. MIPSTR is applicable to any organism with a high-quality reference genome and is scalable to genotyping many thousands of STR loci in thousands of individuals.

Partially redundant enhancers cooperatively maintain Mammalian pomc expression above a critical functional threshold

Cell-specific expression of many genes is conveyed by multiple enhancers, with each individual enhancer controlling a particular expression domain. In contrast, multiple enhancers drive similar expression patterns of some genes involved in embryonic development, suggesting regulatory redundancy. Work in Drosophila has indicated that functionally overlapping enhancers canalize development by buffering gene expression against environmental and genetic disturbances. However, little is known about regulatory redundancy in vertebrates and in genes mainly expressed during adulthood. Here we study nPE1 and nPE2, two phylogenetically conserved mammalian enhancers that drive expression of the proopiomelanocortin gene (Pomc) to the same set of hypothalamic neurons. The simultaneous deletion of both enhancers abolished Pomc expression at all ages and induced a profound metabolic dysfunction including early-onset extreme obesity. Targeted inactivation of either nPE1 or nPE2 led to very low levels of Pomc expression during early embryonic development indicating that both enhancers function synergistically. In adult mice, however, Pomc expression is controlled additively by both enhancers, with nPE1 being responsible for ∼80% and nPE2 for ∼20% of Pomc transcription. Consequently, nPE1 knockout mice exhibit mild obesity whereas nPE2-deficient mice maintain a normal body weight. These results suggest that nPE2-driven Pomc expression is compensated by nPE1 at later stages of development, essentially rescuing the earlier phenotype of nPE2 deficiency. Together, these results reveal that cooperative interactions between the enhancers confer robustness of Pomc expression against gene regulatory disturbances and preclude deleterious metabolic phenotypes caused by Pomc deficiency in adulthood. Thus, our study demonstrates that enhancer redundancy can be used by genes that control adult physiology in mammals and underlines the potential significance of regulatory sequence mutations in common diseases.

Twins and endocrinology

Twins are two independent babies delivered during the same pregnancy and are divided as monozygotic or dizygotic based on their origin. Dizygotic twins are similar to two siblings and have different genetic information. In contrary, monozygotic twins have a similar genetic identity and provide a unique opportunity to evaluate the contribution of genetic and environmental factors of the disease. The endocrine and metabolic disorders affect a large number of the population including the twins. Diabetes, obesity, and autoimmune thyroid disease are the most common endocrine disorders in general practice. It is essential to understand the genetic basis of endocrine disorders for therapy, prognostication and risk assessment for future generations. In this article, we review the endocrine disorders in relation to their occurrence in monozygotic twins to highlight the genetic and environmental contribution.

The overdue promise of short tandem repeat variation for heritability

Short tandem repeat (STR) variation has been proposed as a major explanatory factor in the heritability of complex traits in humans and model organisms. However, we still struggle to incorporate STR variation into genotype-phenotype maps. We review here the promise of STRs in contributing to complex trait heritability and highlight the challenges that STRs pose due to their repetitive nature. We argue that STR variants are more likely than single-nucleotide variants to have epistatic interactions, reiterate the need for targeted assays to genotype STRs accurately, and call for more appropriate statistical methods in detecting STR-phenotype associations. Lastly, we suggest that somatic STR variation within individuals may serve as a read-out of disease susceptibility, and is thus potentially a valuable covariate for future association studies.

Host genetics and viral load in primary HIV-1 infection: clear evidence for gene by sex interactions

Research in the past two decades has generated unequivocal evidence that host genetic variations substantially account for the heterogeneous outcomes following human immunodeficiency virus type 1 (HIV-1) infection. In particular, genes encoding human leukocyte antigens (HLA) have various alleles, haplotypes, or specific motifs that can dictate the set-point (a relatively steady state) of plasma viral load (VL), although rapid viral evolution driven by innate and acquired immune responses can obscure the long-term relationships between HLA genotypes and HIV-1-related outcomes. In our analyses of VL data from 521 recent HIV-1 seroconverters enrolled from eastern and southern Africa, HLA-A*03:01 was strongly and persistently associated with low VL in women (frequency = 11.3 %, P < 0.0001) but not in men (frequency = 7.7 %, P = 0.66). This novel sex by HLA interaction (P = 0.003, q = 0.090) did not extend to other frequent HLA class I alleles (n = 34), although HLA-C*18:01 also showed a weak association with low VL in women only (frequency = 9.3 %, P = 0.042, q > 0.50). In a reduced multivariable model, age, sex, geography (clinical sites), previously identified HLA factors (HLA-B*18, B*45, B*53, and B*57), and the interaction term for female sex and HLA-A*03:01 collectively explained 17.0 % of the overall variance in geometric mean VL over a 3-year follow-up period (P < 0.0001). Multiple sensitivity analyses of longitudinal and cross-sectional VL data yielded consistent results. These findings can serve as a proof of principle that the gap of "missing heritability" in quantitative genetics can be partially bridged by a systematic evaluation of sex-specific associations.

Linking genetics to structural biology: complex heterozygosity screening with actin alanine scan alleles identifies functionally related surfaces on yeast actin

Previous genome-level genetic interaction screens with the single essential actin gene of yeast identified 238 nonessential genes that upon deletion result in deleterious, digenic complex haploinsufficiences with an actin null allele. Deletion alleles of these 238 genes were tested for complex heterozygous interactions with 32 actin alanine scan alleles, which target clusters of residues on the surface of actin. A total of 891 deleterious digenic combinations were identified with 203 of the 238 genes. Two-dimensional hierarchical cluster analysis of the interactions identified nine distinct groups, and the alleles within clusters tended to affect localized regions on the surface of actin. The mutants in one cluster all affect electrostatic interactions between stacked subunits in the long pitch helix of the actin filament. A second cluster that contains the most highly interactive alleles may disrupt the tropomyosin/myosin system, as one of the mutants in that cluster cannot support Type V myosin-dependent movement of secretory vesicles in haploids and causes processivity defects in heterozygous diploids. These examples suggest the clusters represent mutations with shared protein−protein interaction defects. These results show that complex heterozygous interaction screens have benefit for detecting actin-related genes and suggest that having actin filaments of mixed composition, containing both mutant and wild-type subunits, presents unique challenges to the cell.

Chromatin regulators, phenotypic robustness, and autism risk

Though extensively characterized clinically, the causes of autism spectrum disorder (ASD) remain a mystery. ASD is known to have a strong genetic basis, but it is genetically very heterogeneous. Recent studies have estimated that de novo disruptive mutations in hundreds of genes may contribute to ASD. However, it is unclear how it is possible for mutations in so many different genes to contribute to ASD. Recent findings suggest that many of the mutations disrupt genes involved in transcription regulation that are expressed prenatally in the developing brain. De novo disruptive mutations are also more frequent in girls with ASD, despite the fact that ASD is more prevalent in boys. In this paper, we hypothesize that loss of robustness may contribute to ASD. Loss of phenotypic robustness may be caused by mutations that disrupt capacitors that operate in the developing brain. This may lead to the release of cryptic genetic variation that contributes to ASD. Reduced robustness is consistent with the observed variability in expressivity and incomplete penetrance. It is also consistent with the hypothesis that the development of the female brain is more robust, and it may explain the higher rate and severity of disruptive de novo mutations in girls with ASD.

Natural genetic variation determines susceptibility to aggregation or toxicity in a C. elegans model for polyglutamine disease

Background

Monogenic gain-of-function protein aggregation diseases, including Huntington’s disease, exhibit substantial variability in age of onset, penetrance, and clinical symptoms, even between individuals with similar or identical mutations. This difference in phenotypic expression of proteotoxic mutations is proposed to be due, at least in part, to the variability in genetic background. To address this, we examined the role of natural variation in defining the susceptibility of genetically diverse individuals to protein aggregation and toxicity, using the Caenorhabditis elegans polyglutamine model.

Results

Introgression of polyQ40 into three wild genetic backgrounds uncovered wide variation in onset of aggregation and corresponding toxicity, as well as alteration in the cell-specific susceptibility to aggregation. To further dissect these relationships, we established a panel of 21 recombinant inbred lines that showed a broad range of aggregation phenotypes, independent of differences in expression levels. We found that aggregation is a transgressive trait, and does not always correlate with measures of toxicity, such as early onset of muscle dysfunction, egg-laying deficits, or reduced lifespan. Moreover, distinct measures of proteotoxicity were independently modified by the genetic background.

Conclusions

Resistance to protein aggregation and the ability to restrict its associated cellular dysfunction are independently controlled by the natural variation in genetic background, revealing important new considerations in the search for targets for therapeutic intervention in conformational diseases. Thus, our C. elegans model can serve as a powerful tool to dissect the contribution of natural variation to individual susceptibility to proteotoxicity.

Please see related commentary by Kaeberlein, http://www.biomedcentral.com/1741-7007/11/102.

Genetics of type 2 diabetes

Type 2 diabetes (T2D) is the result of interaction between environmental factors and a strong hereditary component. We review the heritability of T2D as well as the history of genetic and genomic research in this area. Very few T2D risk genes were identified using candidate gene and linkage-based studies, but the advent of genome-wide association studies has led to the identification of multiple genes, including several that were not previously known to play any role in T2D. Highly replicated genes, for example TCF7L2, KCNQ1 and KCNJ11, are discussed in greater detail. Taken together, the genetic loci discovered to date explain only a small proportion of the observed heritability. We discuss possible explanations for this “missing heritability”, including the role of rare variants, gene-environment interactions and epigenetics. The clinical utility of current findings and avenues of future research are also discussed.