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.

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.

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 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.

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.

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.

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.

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.

Genetic architecture of isolation between two species of Silene with sex chromosomes and Haldane's rule

Examination of the genetic architecture of hybrid breakdown can provide insight into the genetic mechanisms of commonly observed isolating phenomena such as Haldane's rule. We used line-cross analysis to dissect the genetic architecture of divergence between two plant species that exhibit Haldane's rule for male sterility and rarity, Silene latifolia and Silene diclinis. We made 15 types of crosses, including reciprocal F1, F2, backcrosses, and later-generation crosses, grew the seeds to flowering, and measured the number of viable ovules, proportion of viable pollen, and sex ratio. Typically, Haldane's rule for male rarity in XY animal hybrids is explained by interactions involving recessive X-linked alleles that are deleterious when hemizygous (dominance theory), whereas sterility is explained by rapid evolution of spermatogenesis genes (faster-male evolution). In contrast, we found that the genetic mechanisms underlying Haldane's rule between the two Silene species did not follow these conventions. Dominance theory was sufficient to explain male sterility, but male rarity likely involved faster-male evolution. We also found an effect of the neo-sex chromosomes of S. diclinis on the extreme rarity of some hybrid males. Our findings suggest that the genetic architecture of Haldane's rule in dioecious plants may differ from those commonly found in animals.

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.

Sex, war, and disease: the role of parasite infection on weapon development and mating success in a horned beetle (Gnatocerus cornutus)

While parasites and immunity are widely believed to play important roles in the evolution of male ornaments, their potential influence on systems where male weaponry is the object of sexual selection is poorly understood. We experimentally infect larval broad-horned flour beetles with a tapeworm and study the consequent effects on: 1) adult male morphology 2) male-male contests for mating opportunities, and 3) induction of the innate immune system. We find that infection significantly reduces adult male size in ways that are expected to reduce mating opportunities in nature. The sum of our morphological, competition, and immunological data indicate that during a life history stage where no new resources are acquired, males allocate their finite resources in a way that increases future mating potential.

Hyperexpression of the X chromosome in both sexes results in extensive female bias of X-linked genes in the flour beetle

A genome's ability to produce two separate sexually dimorphic phenotypes is an intriguing biological mystery. Microarray-based studies of a handful of model systems suggest that much of the mystery can be explained by sex-biased gene expression evolved in response to sexually antagonistic selection. We present the first whole-genome study of sex-biased expression in the red flour beetle, Tribolium castaneum. Tribolium is a model for the largest eukaryotic order, Coleoptera, and we show that in whole-body adults, ∼20% of the transcriptome is differentially regulated between the sexes. Among T. castaneum, Drosophila melanogaster, and Anopheles gambiae, we identify 416 1:1:1 orthologs with conserved sex-biased expression. Overrepresented functional categories among sex-biased genes are primarily those involved in gamete production and development. The genomic distribution of sex-biased genes in T. castaneum is distinctly nonrandom, with the strongest deficit of male-biased genes on the X chromosome (9 of 793) of any species studied to date. Tribolium also shows a significant enrichment of X-linked female-biased genes (408 of 793). Our analyses suggest that the extensive female bias of Tribolium X chromosome gene expression is due to hyperexpression of X-linked genes in both males and females. We propose that the overexpression of X chromosomes in females is an evolutionary side effect of the need to dosage compensate in males and that mechanisms to reduce female X chromosome gene expression to autosomal levels are sufficient but imperfect.

Adaptive evolution of young gene duplicates in mammals

Duplicate genes act as a source of genetic material from which new functions arise. They exist in large numbers in every sequenced eukaryotic genome and may be responsible for many differences in phenotypes between species. However, recent work searching for the targets of positive selection in humans has largely ignored duplicated genes due to complications in orthology assignment. Here we find that a high proportion of young gene duplicates in the human, macaque, mouse, and rat genomes have experienced adaptive natural selection. Approximately 10% of all lineage-specific duplicates show evidence for positive selection on their protein sequences, larger than any reported amount of selection among single-copy genes in these lineages using similar methods. We also find that newly duplicated genes that have been transposed to new chromosomal locations are significantly more likely to have undergone positive selection than the ancestral copy. Human-specific duplicates evolving under adaptive natural selection include a surprising number of genes involved in neuronal and cognitive functions. Our results imply that genome scans for selection that ignore duplicated loci are missing a large fraction of all adaptive substitutions. The results are also in agreement with the classical model of evolution by gene duplication, supporting a common role for neofunctionalization in the long-term maintenance of gene duplicates.

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.

The genome of the model beetle and pest Tribolium castaneum

Tribolium castaneum is a member of the most species-rich eukaryotic order, a powerful model organism for the study of generalized insect development, and an important pest of stored agricultural products. We describe its genome sequence here. This omnivorous beetle has evolved the ability to interact with a diverse chemical environment, as shown by large expansions in odorant and gustatory receptors, as well as P450 and other detoxification enzymes. Development in Tribolium is more representative of other insects than is Drosophila, a fact reflected in gene content and function. For example, Tribolium has retained more ancestral genes involved in cell-cell communication than Drosophila, some being expressed in the growth zone crucial for axial elongation in short-germ development. Systemic RNA interference in T. castaneum functions differently from that in Caenorhabditis elegans, but nevertheless offers similar power for the elucidation of gene function and identification of targets for selective insect control.

Evolutionary and biomedical insights from the rhesus macaque genome

The rhesus macaque (Macaca mulatta) is an abundant primate species that diverged from the ancestors of Homo sapiens about 25 million years ago. Because they are genetically and physiologically similar to humans, rhesus monkeys are the most widely used nonhuman primate in basic and applied biomedical research. We determined the genome sequence of an Indian-origin Macaca mulatta female and compared the data with chimpanzees and humans to reveal the structure of ancestral primate genomes and to identify evidence for positive selection and lineage-specific expansions and contractions of gene families. A comparison of sequences from individual animals was used to investigate their underlying genetic diversity. The complete description of the macaque genome blueprint enhances the utility of this animal model for biomedical research and improves our understanding of the basic biology of the species.

POPULATION DIFFERENTIATION IN THE BEETLE TRIBOLIUM 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.

The evolution of mammalian gene families

Gene families are groups of homologous genes that are likely to have highly similar functions. Differences in family size due to lineage-specific gene duplication and gene loss may provide clues to the evolutionary forces that have shaped mammalian genomes. Here we analyze the gene families contained within the whole genomes of human, chimpanzee, mouse, rat, and dog. In total we find that more than half of the 9,990 families present in the mammalian common ancestor have either expanded or contracted along at least one lineage. Additionally, we find that a large number of families are completely lost from one or more mammalian genomes, and a similar number of gene families have arisen subsequent to the mammalian common ancestor. Along the lineage leading to modern humans we infer the gain of 689 genes and the loss of 86 genes since the split from chimpanzees, including changes likely driven by adaptive natural selection. Our results imply that humans and chimpanzees differ by at least 6% (1,418 of 22,000 genes) in their complement of genes, which stands in stark contrast to the oft-cited 1.5% difference between orthologous nucleotide sequences. This genomic “revolving door” of gene gain and loss represents a large number of genetic differences separating humans from our closest relatives.

Experimental Methods for Measuring Gene Interactions

The role of epistasis in evolution has long been contentious. Resolving the issue requires empirical measurements that are statistically adequate and evolutionarily relevant. We review experimental methods for measuring epistasis, some that are commonly used but weak and others that are less frequently used but stronger. We review statistical genetic methods based on analyses of variances and means as well as molecular genetic methods for detecting gene interactions. We also highlight relevant empirical studies that illustrate the implementation of particular methods. In spite of the inherent weaknesses of most methods, epistasis is surprisingly common. We conclude with a discussion of how technologies for investigating genome-wide epistasis are bridging the gap between physiological and statistical epistasis for model organisms.

Maternal expression increases the rate of bicoid evolution by relaxing selective constraint

Population genetic theory predicts that maternal effect genes will evolve differently than genes expressed in both sexes because selection is only half as effective on autosomal genes expressed in one sex but not the other. Here, we use sequences of the tandem gene duplicates, bicoid (bcd) and zerknüllt (zen), to test the prediction that, with similar coefficients of purifying selection, a maternal effect gene evolves more rapidly than a zygotic gene because of this reduction in selective constraint. We find that the maternal effect gene, bcd, is evolving more rapidly than zygotically expressed, zen, providing the first direct confirmation of this prediction of maternal effect theory from molecular evidence. Our results extend current explanations for the accelerated rate of bcd evolution by providing an evolutionary mechanism, relaxed selective constraint, that allows bcd the evolutionary flexibility to escape the typical functional constraints of early developmental genes. We discuss general implications of our findings for the role of maternal effect genes in early developmental patterning.

CAFE: a computational tool for the study of gene family evolution

SUMMARY

We present CAFE (Computational Analysis of gene Family Evolution), a tool for the statistical analysis of the evolution of the size of gene families. It uses a stochastic birth and death process to model the evolution of gene family sizes over a phylogeny. For a specified phylogenetic tree, and given the gene family sizes in the extant species, CAFE can estimate the global birth and death rate of gene families, infer the most likely gene family size at all internal nodes, identify gene families that have accelerated rates of gain and loss (quantified by a p-value) and identify which branches cause the p-value to be small for significant families.

AVAILABILITY

Software is available from http://www.bio.indiana.edu/~hahnlab/Software.html

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.

How do novel structures and functions originate? This question has proven more difficult to answer than the question of how existing structures are refined to better suit the environment. Evolution by natural selection explains the latter. Ironically, it is the power of natural selection to maintain genes that are suitable for their current roles that also works against the evolution of entirely new traits. The conservation of genes controlling the early stages of embryo development, from worms, to flies, to humans, is a famous example of natural selection's power to constrain evolution and thwart the spread of evolutionary novelties. However, the gene determining which end of the fly embryo becomes the head has a relatively recent origin. How did this gene, bicoid, escape purifying selection and take on a novel function? The authors investigate the hypothesis that because bicoid is expressed only in females it experiences only half as much constraint as a gene expressed in both sexes. Comparing sequences of bicoid with its duplicate gene zerknüllt, which retains expression in both sexes, the authors show that, as expected, the variation in bicoid is twice that of zerknüllt. The findings suggest that relaxed constraint is an important step in the origin of novel function.

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.

Sexual Selection Favors Female‐Biased Sex Ratios: The Balance between the Opposing Forces of Sex‐Ratio Selection and Sexual Selection

In a verbal model, Trivers and Willard proposed that, whenever there is sexual selection among males, natural selection should favor mothers that produce sons when in good condition but daughters when in poor condition. The predictions of this model have been the subject of recent debate. We present an explicit population genetic model for the evolution of a maternal-effect gene that biases offspring sex ratio. We show that, like local mate competition, sexual selection favors female-biased sex ratios whenever maternal condition affects the reproductive competitive ability of sons. However, Fisherian sex-ratio selection, which favors a balanced sex ratio, is an opposing force. We show that the evolution of maternal sex-ratio biasing by these opposing selection forces requires a positive covariance across environments between the sex-ratio bias toward sons (b) and the mating success of sons (r). This covariance alone is not a sufficient condition for the evolution of maternal sex-ratio biasing; it must be sufficiently positive to outweigh the opposing sex-ratio selection. To identify the necessary and sufficient conditions, we partition total evolutionary change into three components: (1) maternal sex-ratio bias, (2) sexual selection on sons, and (3) sex-ratio selection. Because the magnitude of the first component asymmetrically affects the strength of the second, biasing broods toward females in a poor environment evolves faster than the same degree of bias toward males in a good environment. Consequently, female-biased sex ratios, rather than male-biased sex ratios, are more likely to evolve. We discuss our findings in the context of the primary sex-ratio biases observed in strongly sexually selected species and indicate how this perspective can assist the experimental study of sex ratio evolution.

The effects of constant and fluctuating incubation temperatures on sex determination, growth, and performance in the tortoise Gopherus polyphemus

Temperature-dependent sex determination is one of the best documented yet evolutionarily enigmatic sex-determining systems. The classical theoretical framework suggests that temperature-dependent sex determination will be adaptive when males and females benefit differentially from development at certain temperatures. Empirical evidence has not provided convincing support for this "differential-fitness" hypothesis. Furthermore, since most experiments utilize constant temperature incubation treatments to explore phenotypic response to temperature, few studies have addressed the consequences of incubation under natural conditions. In this study I utilized constant-temperature laboratory incubations and natural-nest incubations to determine the effects of temperature on sex, size, growth, and locomotor performance in the tortoise Gopherus polyphemus. Constant-temperature incubations do induce substantial growth and performance variation in these tortoises. However, the data do not clearly support the differential-fitness hypothesis because (i) growth variation does not result in adult size dimorphism, (ii) performance differences are confined to a very short period after hatching, and (iii) natural incubation temperatures do not vary sufficiently to produce significant phenotypic variation in traits other than sex.

Localization of tumor suppressor gene candidates by cytogenetic and short tandem repeat analyses in tumorigenic human bronchial epithelial cells

Radon exposure is associated with increased risk for bronchogenic carcinoma. Mutagenesis analyses have revealed that radon induces mostly multi-locus chromosome deletions. Based on these findings, it was hypothesized that deletion analysis of multiple radon-induced malignant transformants would reveal common mutations in chromosomal regions containing tumor suppressor genes responsible for malignant transformation. This hypothesis was supported by a previous study in which tumorigenic derivatives of the human papillomavirus 18-immortalized human bronchial epithelial cell line BEP2D were established following irradiation with 30 cGy of high linear energy transfer radon-simulated alpha-particles. Herein, we describe the analyses of 10 additional tumorigenic derivative cell lines resulting from the irradiation of five additional independent BEP2D populations. The new transformants have common cytogenetic changes, including the loss of chromosome (ch)Y, one of three copies of ch8, one of two copies of ch11p15-pter and one of three copies of ch14. These changes are the same as those reported previously. Analysis of PCR-amplified short tandem repeats of informative loci confirmed the loss of heterozygosity (LOH) at 12 loci spanning the length of ch8 in cell lines from four of the total of eight irradiation treatments to date and the loss of chY in all cell lines (8 of 8). LOH analysis with a total of 17 informative loci confirmed loss on ch14 in transformants from seven of eight irradiation treatments and indicated a 0.5-1.7 cM region of common involvement centered around locus D14S306. No LOH was detected at any of the informative loci on ch11. The overall results support our stated hypothesis. Further studies are currently in progress to determine whether the ch8 and ch14 regions contain genes with tumor suppressor function in bronchial epithelial cells.