Probing the Depths of Biological Diversity During the Second Century of<i>GENETICS</i>
Bibliographic record
Abstract
After a century of GENETICS, we understand better than ever the diversity of life and its immense evolutionary history. Nevertheless, we are still gazing at the tip of the iceberg of biological complexity. For virtually every biological rule, an exception lies in some organism on some branch of the tree of life. Meiosis is fair. Mating is random. Chromosomes govern inheritance. The genetic code is universal. All have exceptions (meiotic drive: Buckler et al. 1999 and Didion et al. 2016; mating: Jiang et al. 2013; inheritance: Fang et al. 2012 and Houri-Ze’evi et al. 2016; genetic code: Saccone et al. 2000). In contemplating what is in store for the journal GENETICS in its second century, we argue that our vision of genetics will move increasingly away from trying to understand the general pattern of biology toward grappling with its variability and, in so doing, better reveal the depths of biological complexity. Biological complexity is often discussed as the product of an evolutionary history spanning the ∼4.1 billion years since the origin of life (Bell et al. 2015), but even this number is misleadingly small. Evolution is not linear: it branches into species, which explore, in parallel, different ways of surviving and reproducing. This exploration spans more than 1014 years of evolutionary discovery (Figure 1), a staggering number. By comparison, there are an estimated 1014 letters in total in all of the published books across human history (Urban 2014). It is no wonder that life is so variable. Summing the total evolutionary history of life. The tree of life is not completely resolved, so the sum total evolution that it spans is unknown. A simple birth model would start with a single species 4.1 billion years ago, undergoing x = 22.25 species doublings, to produce the ∼5 million species alive today (Costello et al. 2013). With the length of each branch spanning τ = 4.1 × 109/x billion years, the sum total amount of evolutionary time captured by the tree of life would be 1.5 × 1015 years (∑i=0x2xτ). Account for extinction (using equation 29b in Nee et al. 1994) and assuming that extinction occurs at 90% the rate of speciation (e.g., Rabosky et al. 2013), the span of evolutionary branches would total 4.0 × 1014 years. Contemplating this breadth of evolutionary history is essential if one is to understand the richness of biology. Historically, we have done the reverse. We have stripped out the complexities, alternative contexts, and species interactions. This is a natural and necessary thing to do when first seeking out the impact of what a gene does or how a population evolves, but a full understanding of biology requires that we expand our viewpoint and consider alternative contexts and understudied organisms from across the tree of life. Why do genes not always perform in the same way? Why do populations not always evolve in the same direction? Although biologists, both theoreticians and empiricists, have been moving in this direction for decades, we argue that in the next century, our focus will shift from a search for general rules to a greater appreciation of biological variability. We describe ways in which we expect this shift to impact genetics and evolution. Over the past century, geneticists have identified possible functions of most genes in Drosophila melanogaster, Saccharomyces cerevisiae, Arabidopsis thaliana, and in our own species, Homo sapiens—a tremendous feat. Furthermore, the commitment to sharing this knowledge through online databases has greatly enabled syntheses and has proven the value of community cooperation for scientific progress (e.g., flybase.org, yeastgenome.org, arabidopsis.org, omim.org, and encodeproject.org). These compilations have significantly spurred scientific discovery, as evidenced by the numerous citations mentioning the “Saccharomyces Genome Database” (over 10,000 on Google Scholar) and “Online Mendelian Inheritance in Man” (over 45,200 on Google Scholar). While genetic data are increasingly available, there are several reasons to believe that much remains unknown. Take for example, gene function. Genes are typically named for their inferred function, often based on mutant screens in a homogenous or inbred lineage. Arguably the best understood eukaryote genetically, S. cerevisiae has 5820 open reading frames (ORFs) (not considering dubious ORFs). The vast majority of these ORFs (88.5%) have now been characterized to some extent (Figure 2). Yet our tendency to name and discuss genes as if they had one function is often unfounded. For example, even considering just coarse Gene Ontology (GO)-Slim categories (Figure 2), 35% of genes have more than one molecular function, 58% are involved in more than one biological process, and much remains to be discovered. The map from gene to function is, in many cases, multipronged and pleiotropic. Most yeast genes have been well characterized with some information about the associated cellular component (88.52%), molecular function (68.44%), and biological process (82.15%). The figure charts the number of GO terms associated with each gene, focusing on “GO-Slim” terms, a broad-level categorization in the ontology (e.g., “vitamin metabolic process” or “transmembrane transport”). Data are from the Saccharomyces Genome Database. Available at: http://downloads.yeastgenome.org/curation/literature/go_slim_mapping.tab. As another example, consider protein–protein interaction networks. The network for yeast was initially inferred in two different studies, which exhibited little overlap in the list of interactions (Uetz et al. 2000; Ito et al. 2001), indicating that many interactions remained undetected. Accounting statistically for the overlap, Grigoriev (2003) concluded that the average yeast protein was involved in five different protein–protein interactions. Consequently, for many genes, the most common of its interactions, as well as the best studied of its functions, may only be scratching the surface; the effects of the gene may be substantially broader. Similarly, in humans, Venkatesan et al. (2009) predicted ∼130,000 protein–protein interactions, the vast majority of which remain undocumented. The fact that some genes, and their interactions, are better studied than others can generate biased conclusions about the biological connectedness of the genes that we study, for example, when comparing cancer-causing vs. control genes (Schaefer et al. 2015). Thus, in addition to being multipronged, the map from gene to function is webbed and contingent upon the state of many interacting genes (epistasis). Even this webbed image of interacting genes is misleadingly static, given that gene expression varies among tissues, between the sexes, and over time (Kwekel et al. 2013; Gutierrez-Arcelus et al. 2015). The map from gene to function is dynamic and responsive to environmental (Scheiner and Lyman 1989; Lacaze et al. 2009) and social cues (Evans and Wheeler 2001; Toth et al. 2007). The lack of full information about the various roles of a gene hampers our ability to predict its potential effects. For example, in stickleback fish, the signaling gene Ectodysplasin (Eda) is well known for its impact on armored plates, thought to be an important defense against predators (Colosimo et al. 2005). Yet, when a polymorphic population was tracked over time, substantial selection was found on plated and nonplated variants of Eda among juvenile fish, before the development of armor (Barrett et al. 2008). Exactly what role(s) the protein EDA plays in these juvenile fish is unknown, but there are many possibilities given that Eda is in a gene family known to impact the development of various ectoderm-derived body parts in mammals (e.g., teeth, hair, and sweat glands; Colosimo et al. 2005). Over the past century, we have attached an ever-growing stack of exceptions and anomalies into a core narrative of simple Mendelian inheritance and a one-gene, one-enzyme, one-function paradigm. It is time to recognize the exceptions for what they are: the outcome of an evolutionary process that need not stick to simple rules. A Gestalt switch, already under way, is needed that grapples with the tangled and dynamic web of genes and their interactions. As we move forward over the next century to understand the immensity of what remains unknown, open repositories of data, synthetic analyses, and probabilistic models will become increasingly important means to capture the diversity of genetic narratives. Analyses such as those mentioned above that compare the overlap among studies (Grigoriev 2003) or that quantify the likelihood of missing data based on how much different genes have been studied (Schaefer et al. 2015) are important ways of identifying gaps. Furthermore, predictive models that combine knowledge of gene functions across species with models of trait evolution across phylogenies would allow probabilistic statements to be made about the likelihood that a gene plays a given role in a given species (e.g., Engelhardt et al. 2005). The last century has also seen a flourishing of evolutionary knowledge, with the development of a strong mathematical foundation and the establishment of a strong experimental and comparative toolkit. Yet evolution is often portrayed in a simple, hill-climbing fashion, according to the core premises of evolution by natural selection (Lewontin 1974): given that individuals vary in a trait, that some variants are more likely to survive and reproduce, and that variants can be inherited from parents to offspring, the offspring generation will differ from the parental generation. Evolution will have occurred. In the simplest of models (in asexual populations, in sexual populations with one gene, or in sexual populations with an infinite number of small-effect genes, all within a constant environment), evolution follows the trajectory that maximizes fitness: populations become successively better adapted to their current environment. By focusing on evolution as a hill-climbing mechanism, however, it becomes difficult to fathom why the biological world is so diverse. Surely, one might think, there would be an optimal genome structure or mating system, a hill that evolution climbed long ago. But this singular view belies reality, and much of evolutionary biology today is grappling with the proliferation of evolutionary trajectories. Why do some organisms reproduce sexually, others asexually (Simon et al. 2003; Whitton et al. 2008)? Why do some organisms separate male and female functions into different individuals, while others combine the sexes (Figure 3) (Tree of Sex Consortium 2014a)? Why do some organisms develop directly, while others proceed through larval stages (Gomez-Mestre et al. 2012)? Why do some organisms have small genomes, others large (Lynch 2007; Elliott and Gregory 2015)? Variation in chromosome number (“Diploid Num.”) and mode of reproduction across the vertebrate tree of life. XY/ZW, chromosomal sex determination (XY in blue; ZW in orange). ESD, environmental sex determination; hermaphrodite, individuals can serve as both males and females; asexual, reproduction without sexual mixing of genomes (including parthenogenesis, gynogenesis, and hybridogenesis). Modified from figure 3 of Tree of Sex Consortium (2014a). Consider, as a special case, the diversity of ploidy levels that underlies several of the quantum leaps in genome size taken by evolution. Dramatically, some organisms have recently experienced the addition of an entire genomic copy (polyploidization), including some animals and numerous plants; polyploidy is particularly common in domesticated crops such as coffee, sugar, wheat, and cotton (Otto and Whitton 2000). Another facet of diversity in ploidy arises from the alternation between haploid and diploid phases in sexually reproducing eukaryotes. Some organisms, including many fungi, algae, and bryophytes, spend the majority of their lives in the haploid phase; others, for instance most vascular plants and animals, in the diploid phase. However, this summary barely scratches the surface of ploidy diversity. Some species alternate between free-living haploid and diploid phases (e.g., red algae). Some species have haploid males and diploid females (e.g., Hymenoptera). Some species start their lives at one ploidy level and end at another (e.g., paternal genome elimination; Gardner and Ross 2014). Why? Empirically, evolutionary biologists examine environmental correlates for clues about the conditions that might favor one ploidy level over another. For example, Lewis (1985) proposed that nutrient limitation may favor haploids over diploids, due to their higher surface area-to-volume ratio. For similar reasons, it has been suggested that toxic environments may favor diploidy (Zörgö et al. 2013). One particularly striking pattern is that haploidy is rare among large-bodied organisms, suggesting that diploidy may have evolved to mask somatic mutations and reduce the risk of cancer among multicellular organisms (Orr 1995; Mable and Otto 1998). On the other hand, many organisms are not dominated by one ploidy phase, and several authors have explored the conditions that maintain both phases, including niche differentiation (Hughes and Otto 1999; Rescan et al. 2016) or contrasting dispersal vs. growth advantages (Bell 1997). How does evolutionary theory move beyond the expected trajectory to grapple with the diversity of life that has evolved? One approach is to examine how genes evolve to alter characteristics of the biological system, such as its physical genome, the genetic map, the mutation rate, and the mode of reproduction. Such models can then be used to predict how different initial conditions (either biotic or abiotic) can generate different outcomes. For example, evolutionary biologists have altered standard models of selection to include genes that modify the alternation of generations between haploid and diploid phases. These models can clarify the conditions under which evolution is predicted to take different paths. For example, selection is generally more efficient in haploid populations than in diploid populations, because alleles are fully expressed. On the other hand, individuals that bear deleterious mutations are much more likely to survive to reproduce if they are diploid and can mask the mutation with a second functional gene copy. These two opposing forces generate conflicting pressures on the life cycle. The evolutionary outcome depends, among other aspects, on the rate of sex and recombination within the population (see review by Mable and Otto 1998). Haploid life cycles are favored in facultatively sexual organisms, when sex and recombination are rare, because genes modifying the life cycle gain from the long-term advantages of more efficient selection. By contrast, diploid life cycles are favored when sex and recombination are common because the modifier soon recombines away from its original lineage and benefits from masking. Making further progress, however, requires that we be able to map life cycles across the tree of life and test predictions from such theories. Major progress has been made recently in the construction of large-scale phylogenies (e.g., of birds, Jetz et al. 2012), progress made possible by the rise of genetic sequence data that are openly available from a large variety of taxa (Benson et al. 2013). Indeed, a first draft of the entire tree of life, combining classical taxonomic information and genetic data, has now been proposed (Hinchliff et al. 2015). Many uncertainties remain, and conflict among data sources is rife across this draft tree, but its existence makes it clearer where the gaps lie, allowing us to target our efforts over the next century to improve the resolution of the tree of life. Alongside phylogeny reconstruction, the development of trait databases [e.g., the Chromosome Counts Database, Rice et al. 2015; genome-size databases, Gregory et al. 2007; the Plant Trait Database (TRY), Kattge et al. 2011; the Tree of Sex Database, Tree of Sex Consortium 2014b; the Encyclopedia of Life, Parr et al. 2014], make it increasingly possible to test predictions and evaluate theories about when and why evolutionary transitions occur. For example, in a recent study, we traced the evolution of sex chromosome fusions along phylogenies of fish and reptiles to evaluate different hypotheses for why fusions establish within populations (neutral drift, chance deleterious events, sexually antagonistic selection; Pennell et al. 2015). Future work will increasingly reveal the context in which evolutionary transitions are favored. Where across the tree of algae do life cycle transitions occur and what features of the organism or its environment best predict when these transitions happen? What explains broad-scale variation in the number of chromosomes or in the genetic map length of chromosomes among taxa? For example, the social Hymenoptera exhibit 10 times higher rates of recombination than their solitary cousins (Ross et al. 2015); do such evolutionary transitions coincide best with changes in effective population size, colony size and structure, or mating system (e.g., monogamy or polyandry)? We are only now building the trees and the databases that will allow the next century of evolutionary genetics to synthesize the data to reveal when and where evolutionary transitions occur. This will not only allow evolutionary hypotheses to be tested, but they will undoubtedly reveal new patterns that require explanation. The past century of genetics has focused on the main effects, the common patterns. The next century of genetics will be the century of diversity, of variance in function, and divergence in the paths taken by evolution. When we characterize genomes and their activity from more than one individual, from unusual species, in different environments, and at various points in development, we will increasingly unveil the unusual roles that genes play in different contexts. We will also come to better understand evolution, not as a singular process, but as a contingent process that can and has led to the endless forms most beautiful on earth. We thank Heath Blackmon for sharing the source code to generate Figure 3; Ailene MacPherson for assistance in generating Figure 1; and Michael Scott, Gil Henriques, and Mark Johnson for helpful comments on the manuscript. We gratefully acknowledge funding from the Natural Sciences and Engineering Research Council of Canada (RGPIN-2016-03711 to S.P.O.) and the Department of Zoology at the University of British Columbia (to L.S.).
Fetched live from OpenAlex and de-inverted. Abstracts are not stored in this database: the inverted indexes are 8.6 GB of the frame’s 9.3 GB of text, and the host has 13 GB free.
How this classification was reachedexpand
Full frame machine prediction
Teacher imitationNot calibrated prevalence, not ground truth. Human validation pending. The Gemma side is a direct model label for every work in the frame, read from the title-only record. The Codex side is a classifier learned from the 10,348 direct Codex labels and calibrated to design-weighted sample rates; fields without enough sample support carry no Codex call. Candidate is the union of the two sides; consensus is their intersection. These outputs are machine_predicted_unvalidated and are not human labels.
Distilled classifier scores by category (both heads)
| Category | Codex | Gemma |
|---|---|---|
| Metaresearch | 0.006 | 0.014 |
| Meta-epidemiology (narrow) | 0.001 | 0.001 |
| Meta-epidemiology (broad) | 0.001 | 0.000 |
| Bibliometrics | 0.002 | 0.002 |
| Science and technology studies | 0.004 | 0.025 |
| Scholarly communication | 0.007 | 0.017 |
| Open science | 0.001 | 0.005 |
| Research integrity | 0.003 | 0.008 |
| Insufficient payload (model declined to judge) | 0.003 | 0.000 |
Machine scores (provisional)
The two teacher heads of the student model, read on this work. A score orders the frame for review; it never asserts a category, and the validation status ships verbatim with every row.
Baseline scores from an immature model (maturity gate not passed, 7 training rounds). Scores rank; they never assert a category.
score_only:v0-immature-baseline · verbatim from the scoring run: score_only means the number may rank works, and no category label ships from itClassification
machine, unvalidatedMachine predicted; a candidate call from one source (direct Gemma or distilled Codex), not a consensus.
How this classification was reached, model by model and score by score, is at the end of the page under "How this classification was reached".