Drift drives the evolution of chromosome number I: The impact of trait transitions on genome evolution in Coleoptera.

Chromosomal mutations such as fusions and fissions are often thought to be deleterious, especially in heterozygotes (underdominant), and consequently are unlikely to become fixed. Yet, many models of chromosomal speciation ascribe an important role to chromosomal mutations. When the effective population size (Ne) is small, the efficacy of selection is weakened, and the likelihood of fixing underdominant mutations by genetic drift is greater. Thus, it is possible that ecological and phenotypic transitions that modulate Ne facilitate the fixation of chromosome changes, increasing the rate of karyotype evolution. We synthesize all available chromosome number data in Coleoptera and estimate the impact of traits expected to change Ne on the rate of karyotype evolution in the family Carabidae and 12 disparate clades from across Coleoptera. Our analysis indicates that in Carabidae, wingless clades have faster rates of chromosome number increase. Additionally, our analysis indicates clades exhibiting multiple traits expected to reduce Ne, including strict inbreeding, oligophagy, winglessness, and island endemism, have high rates of karyotype evolution. Our results suggest that chromosome number changes are likely fixed by genetic drift despite an initial fitness cost and that chromosomal speciation models may be important to consider in clades with very small Ne.

The origins and evolution of chromosomes, dosage compensation, and mechanisms underlying venom regulation in snakes.

Here we use a chromosome-level genome assembly of a prairie rattlesnake (Crotalus viridis), together with Hi-C, RNA-seq, and whole-genome resequencing data, to study key features of genome biology and evolution in reptiles. We identify the rattlesnake Z Chromosome, including the recombining pseudoautosomal region, and find evidence for partial dosage compensation driven by an evolutionary accumulation of a female-biased up-regulation mechanism. Comparative analyses with other amniotes provide new insight into the origins, structure, and function of reptile microchromosomes, which we demonstrate have markedly different structure and function compared to macrochromosomes. Snake microchromosomes are also enriched for venom genes, which we show have evolved through multiple tandem duplication events in multiple gene families. By overlaying chromatin structure information and gene expression data, we find evidence for venom gene-specific chromatin contact domains and identify how chromatin structure guides precise expression of multiple venom gene families. Further, we find evidence for venom gland-specific transcription factor activity and characterize a complement of mechanisms underlying venom production and regulation. Our findings reveal novel and fundamental features of reptile genome biology, provide insight into the regulation of snake venom, and broadly highlight the biological insight enabled by chromosome-level genome assemblies.

Divergence time estimation of genus Tribolium by extensive sampling of highly conserved orthologs.

Tribolium castaneum, the red flour beetle, is among the most well-studied eukaryotic genetic model organisms. Tribolium often serves as a comparative bridge from highly derived Drosophila traits to other organisms. Simultaneously, as a member of the most diverse order of metazoans, Coleoptera, Tribolium informs us about innovations that accompany hyper diversity. However, understanding the tempo and mode of evolutionary innovation requires well-resolved, time-calibrated phylogenies, which are not available for Tribolium. The most recent effort to understand Tribolium phylogenetics used two mitochondrial and three nuclear markers. The study concluded that the genus may be paraphyletic and reported a broad range for divergence time estimates. Here we employ recent advances in Bayesian methods to estimate the relationships and divergence times among Tribolium castaneum, T. brevicornis, T. confusum, T. freemani, and Gnatocerus cornutus using 1368 orthologs conserved across all five species and an independent substitution rate estimate. We find that the most basal split within Tribolium occurred ~86 Mya [95% HPD 85.90-87.04 Mya] and that the most recent split was between T. freemani and T. castaneum at ~14 Mya [95% HPD 13.55-14.00]. Our results are consistent with broader phylogenetic analyses of insects and suggest that Cenozoic climate changes played a role in the Tribolium diversification.

The within-host dynamics of infection in trans-generationally primed flour beetles.

Many taxa exhibit plastic immune responses initiated after primary microbial exposure that provide increased protection against disease-induced mortality and the fitness costs of infection. In several arthropod species, this protection can even be passed from parents to offspring through a phenomenon called trans-generational immune priming. Here, we first demonstrate that trans-generational priming is a repeatable phenomenon in flour beetles (Tribolium castaneum) primed and infected with Bacillus thuringiensis (Bt). We then quantify the within-host dynamics of microbes and host physiological responses in infected offspring from primed and unprimed mothers by monitoring bacterial density and using mRNA-seq to profile host gene expression, respectively, over the acute infection period. We find that priming increases inducible resistance against Bt around a critical temporal juncture where host septicaemic trajectories, and consequently survival, may be determined in unprimed individuals. Our results identify a highly differentially expressed biomarker of priming, containing an EIF4-e domain, in uninfected individuals, as well as several other candidate genes. Moreover, the induction and decay dynamics of gene expression over time suggest a metabolic shift in primed individuals. The identified bacterial and gene expression dynamics are likely to influence patterns of bacterial fitness and disease transmission in natural populations.

The fragile Y hypothesis: Y chromosome aneuploidy as a selective pressure in sex chromosome and meiotic mechanism evolution.

Loss of the Y-chromosome is a common feature of species with chromosomal sex determination. However, our understanding of why some lineages frequently lose Y-chromosomes while others do not is limited. The fragile Y hypothesis proposes that in species with chiasmatic meiosis the rate of Y-chromosome aneuploidy and the size of the recombining region have a negative correlation. The fragile Y hypothesis provides a number of novel insights not possible under traditional models. Specifically, increased rates of Y aneuploidy may impose positive selection for (i) gene movement off the Y; (ii) translocations and fusions which expand the recombining region; and (iii) alternative meiotic segregation mechanisms (achiasmatic or asynaptic). These insights as well as existing evidence for the frequency of Y-chromosome aneuploidy raise doubt about the prospects for long-term retention of the human Y-chromosome despite recent evidence for stable gene content in older non-recombining regions.

Haldane's Rule: Genetic Bases and Their Empirical Support.

There are few patterns in evolution that are as rigidly held as Haldane's rule (HR), which states, "When in the first generation between hybrids between 2 species, 1 sex is absent, rare, or sterile, that sex is always the heterogametic sex." Yet despite considerable attention for almost a century, questions persist as to how many independent examples exist and what is (are) the underlying genetic cause(s). Here, we review recent evidence extending HR to plants, where previously it has only been documented in animals. We also discuss recent comparative analyses that show much more variation in sex-chromosome composition than previously recognized, thus increasing the number of potential independent origins of HR dramatically. Finally, we review the standing of genetic theories proposed to explain HR in light of the new examples and new molecular understanding.

Wright was right: leveraging old data and new methods to illustrate the critical role of epistasis in genetics and evolution.

Much of evolutionary theory is predicated on assumptions about the relative importance of simple additive versus complex epistatic genetic architectures. Previous work suggests traits strongly associated with fitness will lack additive genetic variation, whereas traits less strongly associated with fitness are expected to exhibit more additive genetic variation. We use a quantitative genetics method, line cross analysis, to infer genetic architectures that contribute to trait divergence. By parsing over 1,600 datasets by trait type, clade, and cross divergence, we estimated the relative importance of epistasis across the tree of life. In our comparison between life-history traits and morphological traits, we found greater epistatic contributions to life-history traits. Our comparison between plants and animals showed that animals have more epistatic contribution to trait divergence than plants. In our comparison of within-species versus between-species crosses, we found that only animals exhibit a greater epistatic contribution to trait divergence as divergence increases. While many scientists have argued that epistasis is ultimately of little importance, our results show that epistasis underlies much of trait divergence and must be accounted for in theory and practical applications like domestication, conservation breeding design, and understanding complex diseases.

Gene expression levels are correlated with synonymous codon usage, amino acid composition, and gene architecture in the red flour beetle, Tribolium castaneum.

Gene expression levels correlate with multiple aspects of gene sequence and gene structure in phylogenetically diverse taxa, suggesting an important role of gene expression levels in the evolution of protein-coding genes. Here we present results of a genome-wide study of the influence of gene expression on synonymous codon usage, amino acid composition, and gene structure in the red flour beetle, Tribolium castaneum. Consistent with the action of translational selection, we find that synonymous codon usage bias increases with gene expression. However, the correspondence between tRNA gene copy number and optimal codons is weak. At the amino acid level, translational selection is suggested by the positive correlation between tRNA gene numbers and amino acid usage, which is stronger for highly expressed genes. In addition, there is a clear trend for increased use of metabolically cheaper, less complex amino acids as gene expression increases. tRNA gene numbers also correlate negatively with amino acid size/complexity (S/C) score indicating the coupling between translational selection and selection to minimize the use of large/complex amino acids. Interestingly, the analysis of 10 additional genomes suggests that the correlation between tRNA gene numbers and amino acid S/C score is widespread and might be explained by selection against negative consequences of protein misfolding. At the level of gene structure, three major trends are detected: 1) complete coding region length increases across low and intermediate expression levels but decreases in highly expressed genes; 2) the average intron size shows the opposite trend, first decreasing with expression, followed by a slight increase in highly expressed genes; and 3) intron density remains nearly constant across all expression levels. These changes in gene architecture are only in partial agreement with selection favoring reduced cost of biosynthesis.

Haldane's rule in marsupials: what happens when both sexes are functionally hemizygous?

During the process of speciation, diverging taxa often hybridize and produce offspring wherein the heterogametic sex (i.e., XY or ZW) is unfit (Haldane's rule). Dominance theory seeks to explain Haldane's rule in terms of the difference in X-linked dominance regimes experienced by the sexes. However, X inactivation in female mammals extends the effects of hemizygosity to both sexes. Here, we highlight where the assumptions of dominance theory are particularly problematic in marsupials, where X inactivation uniformly results in silencing the paternal X. We then present evidence of Haldane's rule for sterility but not for viability in marsupials, as well as the first violations of Haldane's rule for these traits among all mammals. Marsupials represent a large taxonomic group possessing heteromorphic sex chromosomes, where the dominance theory cannot explain Haldane's rule. In this light, we evaluate alternative explanations for the preponderance of male sterility in interspecific hybrids, including faster male evolution, X-Y interactions, and genomic conflict hypotheses.

The life and death of gene families.

One of the unique insights provided by the growing number of fully sequenced genomes is the pervasiveness of gene duplication and gene loss. Indeed, several metrics now suggest that rates of gene birth and death per gene are only 10-40% lower than nucleotide substitutions per site, and that per nucleotide, the consequent lineage-specific expansion and contraction of gene families may play at least as large a role in adaptation as changes in orthologous sequences. While gene family evolution is pervasive, it may be especially important in our own evolution since it appears that the "revolving door" of gene duplication and loss has undergone multiple accelerations in the lineage leading to humans. In this paper, we review current understanding of gene family evolution including: methods for inferring copy number change, evidence for adaptive expansion and adaptive contraction of gene families, the origins of new families and deaths of previously established ones, and finally we conclude with a perspective on challenges and promising directions for future research.

Accelerated rate of gene gain and loss in primates.

The molecular changes responsible for the evolution of modern humans have primarily been discussed in terms of individual nucleotide substitutions in regulatory or protein coding sequences. However, rates of nucleotide substitution are slowed in primates, and thus humans and chimpanzees are highly similar at the nucleotide level. We find that a third source of molecular evolution, gene gain and loss, is accelerated in primates relative to other mammals. Using a novel method that allows estimation of rate heterogeneity among lineages, we find that the rate of gene turnover in humans is more than 2.5 times faster than in other mammals and may be due to both mutational and selective forces. By reconciling the gene trees for all of the gene families included in the analysis, we are able to independently verify the numbers of inferred duplications. We also use two methods based on the genome assembly of rhesus macaque to further verify our results. Our analyses identify several gene families that have expanded or contracted more rapidly than is expected even after accounting for an overall rate acceleration in primates, including brain-related families that have more than doubled in size in humans. Many of the families showing large expansions also show evidence for positive selection on their nucleotide sequences, suggesting that selection has been important in shaping copy-number differences among mammals. These findings may help explain why humans and chimpanzees show high similarity between orthologous nucleotides yet great morphological and behavioral differences.

Maternal expression relaxes constraint on innovation of the anterior determinant, bicoid.

The origin of evolutionary novelty is believed to involve both positive selection and relaxed developmental constraint. In flies, the redesign of anterior patterning during embryogenesis is a major developmental innovation and the rapidly evolving Hox gene, bicoid (bcd), plays a critical role. We report evidence for relaxation of selective constraint acting on bicoid as a result of its maternal pattern of gene expression. Evolutionary theory predicts 2-fold greater sequence diversity for maternal effect genes than for zygotically expressed genes, because natural selection is only half as effective acting on autosomal genes expressed in one sex as it is on genes expressed in both sexes. We sample an individual from ten populations of Drosophila melanogaster and nine populations of D. simulans for polymorphism in the tandem gene duplicates bcd, which is maternally expressed, and zerknüllt (zen), which is zygotically expressed. In both species, we find the ratio of bcd to zen nucleotide diversity to be two or more in the coding regions but one in the noncoding regions, providing the first quantitative support for the theoretical prediction of relaxed selective constraint on maternal-effect genes resulting from sex-limited expression. Our results suggest that the accelerated rate of evolution observed for bcd is owing, at least partly, to variation generated by relaxed selective constraint.

Molecular Adaptations for Sensing and Securing Prey and Insight into Amniote Genome Diversity from the Garter Snake Genome.

Colubridae represents the most phenotypically diverse and speciose family of snakes, yet no well-assembled and annotated genome exists for this lineage. Here, we report and analyze the genome of the garter snake, Thamnophis sirtalis, a colubrid snake that is an important model species for research in evolutionary biology, physiology, genomics, behavior, and the evolution of toxin resistance. Using the garter snake genome, we show how snakes have evolved numerous adaptations for sensing and securing prey, and identify features of snake genome structure that provide insight into the evolution of amniote genomes. Analyses of the garter snake and other squamate reptile genomes highlight shifts in repeat element abundance and expansion within snakes, uncover evidence of genes under positive selection, and provide revised neutral substitution rate estimates for squamates. Our identification of Z and W sex chromosome-specific scaffolds provides evidence for multiple origins of sex chromosome systems in snakes and demonstrates the value of this genome for studying sex chromosome evolution. Analysis of gene duplication and loss in visual and olfactory gene families supports a dim-light ancestral condition in snakes and indicates that olfactory receptor repertoires underwent an expansion early in snake evolution. Additionally, we provide some of the first links between secreted venom proteins, the genes that encode them, and their evolutionary origins in a rear-fanged colubrid snake, together with new genomic insight into the coevolutionary arms race between garter snakes and highly toxic newt prey that led to toxin resistance in garter snakes.

Population differentiation in the beetleTribolium castaneum. I. Genetic architecture.

We used joint-scaling analyses in conjunction with rearing temperature variation to investigate the contributions of additive, non-additive, and environmental effects to genetic divergence and incipient speciation among 12 populations of the red flour beetle, Tribolium castaneum, with small levels of pairwise nuclear genetic divergence (0.033 < Nei's D < 0.125). For 15 population pairs we created a full spectrum of line crosses (two parental, two reciprocal F1's, four F2's, and eight backcrosses), reared them at multiple temperatures, and analyzed the numbers and developmental defects of offspring. We assayed a total of 219,388 offspring from 5147 families. Failed crosses occurred predominately in F2's, giving evidence of F2 breakdown within this species. In all cases where a significant model could be fit to the data on offspring number, we observed at least one type of digenic epistasis. We also found maternal and cytoplasmic effects to be common components of divergence among T. castaneum populations. In some cases, the most complex model tested (additive, dominance, epistatic, maternal, and cytoplasmic effects) did not provide a significant fit to the data, suggesting that linkage or higher order epistasis is involved in differentiation between some populations. For the limb deformity data, we observed significant genotype-by-environment interaction in most crosses and pure parent crosses tended to have fewer deformities than hybrid crosses. Complexity of genetic architecture was not correlated with either geographic distance or genetic distance. Our results support the view that genetic incompatibilities responsible for postzygotic isolation, an important component of speciation, may be a natural but serendipitous consequence of nonadditive genetic effects and structured populations.

Population differentiation in the beetle Tribolium castaneum. II. Haldane'S rule and incipient speciation.

The heterogametic sex tends to be rare, absent, sterile, or deformed in F1 hybrid crosses between species, a pattern called Haldane's rule (HR). The introgression of single genes or chromosomal regions from one drosophilid species into the genetic background of another have shown that HR is most often associated with fixed genetic differences in inter-specific crosses. However, because such introgression studies have involved species diverged several hundred thousand generations from a common ancestor, it is not clear whether HR attends the speciation process or results from the accumulation of epistatically acting genes postspeciation. We report the first evidence for HR prior to speciation in crosses between two populations of the red flour beetle, Tribolium castaneum, collected 931 km apart in Colombia and Ecuador. In this cross, HR is manifested as an increase in the proportion of deformed males compared to females and the expression of HR is temperature dependent. Neither population, when crossed to a geographically distant population from Japan, exhibits HR at any rearing temperature. Using joint-scaling analysis and additional data from backcrosses and F2's, we find that the hybrid incompatibilities and the emergence of HR are concurrent processes involving interactions between X-linked and autosomal genes. However, we also find many examples of incompatibilities manifest by F2 and backcross hybrids but not by F1 hybrids and most incompatibilities are not sex different in their effects, even when they involve both X-autosomal interactions and genotype-by-environment interactions. We infer that incipient speciation in flour beetles can occur with or without HR and that significant hybrid incompatibilities result from the accumulation of epistatically acting gene differences between populations without differentially affecting the heterogametic sex in F1 hybrids. The temperature dependence of the incompatibilities supports the inference that genotype-by-environment interactions and adaptation to different environments contribute to the genetic divergence important to postzygotic reproductive isolation.

An information-theoretic approach to estimating the composite genetic effects contributing to variation among generation means: Moving beyond the joint-scaling test for line cross analysis.

The pace and direction of evolution in response to selection, drift, and mutation are governed by the genetic architecture that underlies trait variation. Consequently, much of evolutionary theory is predicated on assumptions about whether genes can be considered to act in isolation, or in the context of their genetic background. Evolutionary biologists have disagreed, sometimes heatedly, over which assumptions best describe evolution in nature. Methods for estimating genetic architectures that favor simpler (i.e., additive) models contribute to this debate. Here we address one important source of bias, model selection in line cross analysis (LCA). LCA estimates genetic parameters conditional on the best model chosen from a vast model space using relatively few line means. Current LCA approaches often favor simple models and ignore uncertainty in model choice. To address these issues we introduce Software for Analysis of Genetic Architecture (SAGA), which comprehensively assesses the potential model space, quantifies model selection uncertainty, and uses model weighted averaging to accurately estimate composite genetic effects. Using simulated data and previously published LCA studies, we demonstrate the utility of SAGA to more accurately define the components of complex genetic architectures, and show that traditional approaches have underestimated the importance of epistasis.

Genome of the Asian longhorned beetle (Anoplophora glabripennis), a globally significant invasive species, reveals key functional and evolutionary innovations at the beetle-plant interface.

BACKGROUND: Relatively little is known about the genomic basis and evolution of wood-feeding in beetles. We undertook genome sequencing and annotation, gene expression assays, studies of plant cell wall degrading enzymes, and other functional and comparative studies of the Asian longhorned beetle, Anoplophora glabripennis, a globally significant invasive species capable of inflicting severe feeding damage on many important tree species. Complementary studies of genes encoding enzymes involved in digestion of woody plant tissues or detoxification of plant allelochemicals were undertaken with the genomes of 14 additional insects, including the newly sequenced emerald ash borer and bull-headed dung beetle. RESULTS: The Asian longhorned beetle genome encodes a uniquely diverse arsenal of enzymes that can degrade the main polysaccharide networks in plant cell walls, detoxify plant allelochemicals, and otherwise facilitate feeding on woody plants. It has the metabolic plasticity needed to feed on diverse plant species, contributing to its highly invasive nature. Large expansions of chemosensory genes involved in the reception of pheromones and plant kairomones are consistent with the complexity of chemical cues it uses to find host plants and mates. CONCLUSIONS: Amplification and functional divergence of genes associated with specialized feeding on plants, including genes originally obtained via horizontal gene transfer from fungi and bacteria, contributed to the addition, expansion, and enhancement of the metabolic repertoire of the Asian longhorned beetle, certain other phytophagous beetles, and to a lesser degree, other phytophagous insects. Our results thus begin to establish a genomic basis for the evolutionary success of beetles on plants.

On the theoretical and empirical framework for studying genetic interactions within and among species.

We present a quantitative genetic (QG) interpretation of the Bateson-Dobzhansky-Muller (BDM) genetic model of speciation in order to unify the theoretical framework for understanding how the genetic differentiation of populations is associated with the process of speciation. Specifically, we compare the QG theory of joint scaling with the Turelli-Orr mathematical formulation of the BDM model. By formally linking the two models, we show that a wealth of empirical methods from QG can be brought to bear on the study of the genetic architecture of hybrid phenotypes to better understand the connections, if any, between microevolution within populations and macroevolution in the origin of species. By integrating the two theories, we make additional novel predictions that enrich the opportunities for empirically testing speciation genetic theory or facets of it, such as Haldane's rule. We show that the connection between the two theories is simple and straightforward for autosomal genes but not for sex-linked genes. Differences between the two approaches highlight key conceptual issues concerning the relevance of epistasis to evolution within and among lineages and to differences in the process of speciation in hermaphrodites and in organisms with separate sexes.

Genomic origins of insect sex chromosomes.

Recent efforts to catalog the diversity of sex chromosome systems coupled with genome sequencing projects are adding a new level of resolution to our understanding of insect sex chromosome origins. Y-chromosome degeneration makes sequencing difficult and may erase homology so rapidly that their origins will often remain enigmatic. X-chromosome origins are better understood, but thus far prove to be remarkably labile, often lacking homology even among close relatives. Furthermore, evidence now suggests that differentiated X or Y-chromosomes may both revert to autosomal inheritance. Data for ZW systems is scarcer, but W and Y-chromosomes seem to share many characteristics. Limited evidence suggests that Z-chromosome homology is more conserved than X counterparts, but broader sampling of both sex chromosome systems is needed.

Estimating tempo and mode of Y chromosome turnover: explaining Y chromosome loss with the fragile Y hypothesis.

Chromosomal sex determination is phylogenetically widespread, having arisen independently in many lineages. Decades of theoretical work provide predictions about sex chromosome differentiation that are well supported by observations in both XY and ZW systems. However, the phylogenetic scope of previous work gives us a limited understanding of the pace of sex chromosome gain and loss and why Y or W chromosomes are more often lost in some lineages than others, creating XO or ZO systems. To gain phylogenetic breadth we therefore assembled a database of 4724 beetle species' karyotypes and found substantial variation in sex chromosome systems. We used the data to estimate rates of Y chromosome gain and loss across a phylogeny of 1126 taxa estimated from seven genes. Contrary to our initial expectations, we find that highly degenerated Y chromosomes of many members of the suborder Polyphaga are rarely lost, and that cases of Y chromosome loss are strongly associated with chiasmatic segregation during male meiosis. We propose the "fragile Y" hypothesis, that recurrent selection to reduce recombination between the X and Y chromosome leads to the evolution of a small pseudoautosomal region (PAR), which, in taxa that require XY chiasmata for proper segregation during meiosis, increases the probability of aneuploid gamete production, with Y chromosome loss. This hypothesis predicts that taxa that evolve achiasmatic segregation during male meiosis will rarely lose the Y chromosome. We discuss data from mammals, which are consistent with our prediction.