Theoretical models of the influence of genomic architecture on the dynamics of speciation

A long‐standing problem in evolutionary biology has been determining whether and how gradual, incremental changes at the gene level can account for rapid speciation and bursts of adaptive radiation. Using genome‐scale computer simulations, we extend previous theory showing how gradual adaptive change can generate nonlinear population transitions, resulting in the rapid formation of new, reproductively isolated species. We show that these transitions occur via a mechanism rooted in a basic property of biological heredity: the organization of genes in genomes. Genomic organization of genes facilitates two processes: (i) the build‐up of statistical associations among large numbers of genes and (ii) the action of divergent selection on persistent combinations of alleles. When a population has accumulated a critical amount of standing, divergently selected variation, the combination of these two processes allows many mutations of small effect to act synergistically and precipitously split one population into two discontinuous, reproductively isolated groups. Periods of allopatry, chromosomal linkage among loci, and large‐effect alleles can facilitate this process under some conditions, but are not required for it. Our results complement and extend existing theory on alternative stable states during population divergence, distinct phases of speciation and the rapid emergence of multilocus barriers to gene flow. The results are thus a step towards aligning population genomic theory with modern empirical studies.     

Selection on Wing Allometry in Drosophila Melanogaster

Five bivariate distributions of wing dimensions of Drosophila melanogaster were measured, in flies 1) subjected to four defined environmental regimes during development, 2) taken directly from nature in seven U.S. states, 3) selected in ten populations for change in wing form, and 4) sampled from 21 long inbred wild-type lines. Environmental stresses during development altered both wing size and the ratios of wing dimensions, but regardless of treatment all wing dimensions fell near a common allometric baseline in each bivariate distribution. The wings of wild-caught flies from seven widely separated localities, and of their laboratory-reared offspring, also fell along the same baselines. However, when flies were selected divergently for lateral offset from these developmental baselines, response to selection was rapid in every case. The mean divergence in offset between oppositely selected lines was 14.68 SD of the base population offset, after only 15 generations of selection at 20%. Measurements of 21 isofemale lines, founded from wild-caught flies and maintained in small populations for at least 22 years, showed large reductions in phenotypic variance of offsets within lines, but a large increase in the variance among lines. The variance of means of isofemale lines within collection localities was ten times the variance of means among localities of newly established wild lines. These observations show that much additive genetic variance exists for individual dimensions within the wing, such that bivariate developmental patterns can be changed in any direction by selection or by drift. The relative invariance of the allometric baselines of wing morphology in nature is most easily explained as the result of continuous natural selection around a local optimum of functional design.     

Longevity differences among lines artificially selected for developmental time and wing length in Drosophila buzzatii

We assessed the indirect response of longevity in lines selected for wing length (WL) and developmental time (DT). Longevity in selection lines was compared to laboratory control lines and the offspring of recently collected females. Wild flies (W lines), flies from lines selected for fast development (F lines), and for fast development and large wing length (L lines) outlived control laboratory lines (C lines) and lines selected for fast development and short wing (S lines). The decline in longevity in S lines is in line with the idea that body size and longevity are correlated and may be the result of the fixation of alleles at loci affecting pleiotropically the two traits under selection and longevity. In addition, inbreeding and artificial selection affected the correlation between wing length and longevity that occurs in natural populations of Drosophila buzzatii, suggesting that correlations between traits are not a perdurable feature in a population.     

A Genetic Map of Quantitative Trait Loci for Body Weight in the Mouse

The genetic basis of body weight in the mouse was investigated by measuring frequency changes of microsatellite marker alleles in lines divergently selected for body weight from a base population of a cross between two inbred strains. In several regions of the genome, sharp peaks of frequency change at linked markers were detected, which suggested the presence of single genes of moderate effect, although in several other regions, significant frequency changes occurred over large portions of chromosomes. A method based on maximum likelihood was used to infer effects and map positions of quantitative trait loci (QTLs) based on genotype frequencies at one or more marker loci. Eleven QTLs with effects in the range 0.17-0.28 phenotypic standard deviations were detected; but under an additive model, these did not fully account for the observed selection response. Tests for the presence of more than one QTL in regions where there were large changes of marker allele frequency were mostly inconclusive.     

An Experimental Test of Adaptive Introgression in Locally Adapted Populations of Splash Pool Copepods

As species struggle to keep pace with the rapidly warming climate, adaptive introgression of beneficial alleles from closely related species or populations provides a possible avenue for rapid adaptation. We investigate the potential for adaptive introgression in the copepod, Tigriopus californicus, by hybridizing two populations with divergent heat tolerance limits. We subjected hybrids to strong heat selection for 15 generations followed by whole-genome resequencing. Utilizing a hybridize evolve and resequence (HER) technique, we can identify loci responding to heat selection via a change in allele frequency. We successfully increased the heat tolerance (measured as LT₅₀) in selected lines, which was coupled with higher frequencies of alleles from the southern (heat tolerant) population. These repeatable changes in allele frequencies occurred on all 12 chromosomes across all independent selected lines, providing evidence that heat tolerance is polygenic. These loci contained genes with lower protein-coding sequence divergence than the genome-wide average, indicating that these loci are highly conserved between the two populations. In addition, these loci were enriched in genes that changed expression patterns between selected and control lines in response to a nonlethal heat shock. Therefore, we hypothesize that the mechanism of heat tolerance divergence is explained by differential gene expression of highly conserved genes. The HER approach offers a unique solution to identifying genetic variants contributing to polygenic traits, especially variants that might be missed through other population genomic approaches.     

Investigating the molecular basis of local adaptation to thermal stress: population differences in gene expression across the transcriptome of the copepod Tigriopus californicus

ABSTRACT: BACKGROUND: Geographic variation in the thermal environment impacts a broad range of biochemical and physiological processes and can be a major selective force leading to local population adaptation. In the intertidal copepod Tigriopus californicus, populations along the coast of California show differences in thermal tolerance that are consistent with adaptation, i.e., southern populations withstand thermal stresses that are lethal to northern populations. To understand the genetic basis of these physiological differences, we use an RNA-seq approach to compare genome-wide patterns of gene expression in two populations known to differ in thermal tolerance. RESULTS: Observed differences in gene expression between the southern (San Diego) and the northern (Santa Cruz) populations included both the number of affected loci as well as the identity of these loci. However, the most pronounced differences concerned the amplitude of upregulation of genes producing heat shock proteins (Hsps) and genes involved in ubiquitination and proteolysis. Among the hsp genes, orthologous pairs show markedly different thermal responses as the amplitude of hsp response was greatly elevated in the San Diego population, most notably in members of the hsp70 gene family. There was no evidence of accelerated evolution at the sequence level for hsp genes. Among other sets of genes, cuticle genes were up-regulated in SD but down-regulated in SC, and mitochondrial genes were down-regulated in both populations. CONCLUSIONS: Marked changes in gene expression were observed in response to acute sub-lethal thermal stress in the copepod T. californicus. Although some qualitative differences were observed between populations, the most pronounced differences involved the magnitude of induction of numerous hsp and ubiquitin genes. These differences in gene expression suggest that evolutionary divergence in the regulatory pathway(s) involved in acute temperature stress may offer at least a partial explanation of population differences in thermal tolerance observed in Tigriopus.     

Rflps for Somatotropic Genes Identify Quantitative Trait Loci for Growth in Mice

Restriction fragment length polymorphisms for somatotropic genes were tested for associations with body weight and postweaning growth rate in mice. Polymorphisms for growth hormone (GH) and insulin-like growth factor 2 (IGF-2) genes were identified in stock population lines which had been subjected to long-term selection for high 42-day body weight (H lines) or randomly mated (FP and C lines). Two F2 populations of mice (5F2 and MF2) were generated from crosses between a single H line of mice and two unselected control lines and subsequently, two divergently weight selected sublines were generated from each F2 population. The GHh allele which had originally been fixed in three of four H lines and absent from all FP and C lines was found to have a significant (P less than 0.01) effect on 42-day weight and postweaning growth rate in the F2 populations. However, GHh was associated with lower 42-day weight in the F2 populations, suggesting that the positive association between GHh and weight in the stock population was unique to the high weight selected genetic background of those lines. In agreement with this, the frequency of GHh increased in sublines selected for high 42-day weight and decreased in sublines selected for low 42-day weight. The IGF-2H5 allele was associated with higher weights in a sex-dependent manner in 5F2. In the high selected subline derived from 5F2, a significant increase in the frequency of IGF-2H5 was observed. Therefore this allele, in contrast to GHh, appears to be a positive indicator of growth irrespective of genetic background.     

Effects of cis and trans Genetic Ancestry on Gene Expression in African Americans

Variation in gene expression is a fundamental aspect of human phenotypic variation. Several recent studies have analyzed gene expression levels in populations of different continental ancestry and reported population differences at a large number of genes. However, these differences could largely be due to non-genetic (e.g., environmental) effects. Here, we analyze gene expression levels in African American cell lines, which differ from previously analyzed cell lines in that individuals from this population inherit variable proportions of two continental ancestries. We first relate gene expression levels in individual African Americans to their genome-wide proportion of European ancestry. The results provide strong evidence of a genetic contribution to expression differences between European and African populations, validating previous findings. Second, we infer local ancestry (0, 1, or 2 European chromosomes) at each location in the genome and investigate the effects of ancestry proximal to the expressed gene (cis) versus ancestry elsewhere in the genome (trans). Both effects are highly significant, and we estimate that 12±3% of all heritable variation in human gene expression is due to cis variants.     

The genetic mosaic suggests a new role for hitchhiking in ecological speciation

Early in ecological speciation, the genomically localized effects of divergent selection cause heterogeneity among loci in divergence between incipient species. We call this pattern of genomic variability in divergence the 'genetic mosaic of speciation'. Previous studies have used F(ST) outliers as a way to identify divergently selected genomic regions, but the nature of the relationship between outlier loci and quantitative trait loci (QTL) involved in reproductive isolation has not yet been quantified. Here, we show that F(ST) outliers between a pair of incipient species are significantly clustered around QTL for traits that cause ecologically based reproductive isolation. Around these key QTL, extensive 'divergence hitchhiking' occurs because reduced inter-race mating and negative selection decrease the opportunity for recombination between chromosomes bearing different locally adapted QTL alleles. Divergence hitchhiking is likely to greatly increase the opportunity for speciation in populations that are sympatric, regardless of whether initial divergence was sympatric or allopatric. Early in ecological speciation, analyses of population structure, gene flow or phylogeography based on different random or arbitrarily chosen neutral markers should be expected to conflict--only markers in divergently selected genomic regions will reveal the evolutionary history of adaptive divergence and ecologically based reproductive isolation. Species retain mosaic genomes for a very long time, and gene exchange in hybrid zones can vary dramatically among loci. However, in hybridizing species, the genomic regions that affect ecologically based reproductive isolation are difficult to distinguish from regions that have diverged for other reasons.     

Using genome scans of DNA polymorphism to infer adaptive population divergence

Elucidating the genetic basis of adaptive population divergence is a goal of central importance in evolutionary biology. In principle, it should be possible to identify chromosomal regions involved in adaptive divergence by screening genome-wide patterns of DNA polymorphism to detect the locus-specific signature of positive directional selection. In the case of spatially separated populations that inhabit different environments or sympatric populations that exploit different ecological niches, it is possible to identify loci that underlie divergently selected traits by comparing relative levels of differentiation among large numbers of unlinked markers. In this review I first address the question of whether diversifying selection on polygenic traits can be expected to produce predictable patterns of allelic variation at the underlying quantitative trait loci (QTL), and whether the locus-specific effects of selection can be reliably detected against the genome-wide backdrop of stochastic variability. I then review different approaches that have been developed to identify loci involved in adaptive population divergence and I discuss the relative merits of model-based approaches that rely on assumptions about population structure vs. model-free approaches that are based on empirical distributions of summary statistics. Finally, I consider the evolutionary and functional insights that might be gained by conducting genome scans for loci involved in adaptive population divergence.     

Knockdown resistance to heat stress and slow recovery from chill coma are genetically associated in a quantitative trait locus region of chromosome 2 in Drosophila melanogaster

In insects, two ecologically relevant traits of thermal adaptation are knockdown resistance to high temperature (KRHT) and chill-coma recovery (CCR). Chromosome 2 of Drosophila melanogaster was tested for quantitative trait loci (QTL) affecting both CCR and KRHT in backcrosses between homosequential lines that are fixed for the standard (noninverted) sequence of this autosome. These lines were obtained by artificial selection on KRHT and subsequent inbreeding from a stock that was derived from a single wild population. Heat-induced expression of the 70KD heat-shock protein (Hsp70) was also examined for variation between the lines. Composite interval mapping was performed for each trait on each reciprocal backcross, identifying one QTL region in the middle of chromosome 2 for both KRHT and CCR. The largest estimates of additive effects were found in pericentromeric regions of chromosome 2, accounting for 10–14% (CCR) and 10–17% (KRHT) of the phenotypic variance in BC populations. No QTL was found in the region of the heat-shock factor ( hsf ) gene. However, the two parental lines have diverged in the heat-induced Hsp70 expression. Distribution of KRHT QTL on chromosome 2 was similar between this study based on crosses between lines selected from a single wild population and previous work based on crosses between selection lines from different continents. Colocalized QTL showed a trade–off association between CCR and KRHT, which should be the result of either multiple, tightly linked trait-specific genes or a single gene with pleiotropic effects on the traits. We discuss candidate loci contained within the QTL regions.     

How Small Are the Smallest Selectable Domains of Form?

Two lines of Drosophila melanogaster from the same base population were selected in opposite directions to produce simultaneous antagonistic changes in two very small (less than 0.2 mm) and closely adjacent (less than 0.3 mm) dimensions within the base of the wing. Wing dimensions near the targeted area became differentiated by large positive and negative percentage differences, while only small homogeneous percentage changes occurred in the remainder of the wing. If very small regions of morphology (less than 100 cells across) can respond to selection almost independently, even in small population samples, then the control of developmental detail must involve many genes, and the diversity of possible outcomes in development and adaptation must be large.     

Spatially explicit models of divergence and genome hitchhiking

Strong barriers to genetic exchange can exist at divergently selected loci, whereas alleles at neutral loci flow more readily between populations, thus impeding divergence and speciation in the face of gene flow. However, ‘divergence hitchhiking’ theory posits that divergent selection can generate large regions of differentiation around selected loci. ‘Genome hitchhiking’ theory suggests that selection can also cause reductions in average genome‐wide rates of gene flow, resulting in widespread genomic divergence (rather than divergence only around specific selected loci). Spatial heterogeneity is ubiquitous in nature, yet previous models of genetic barriers to gene flow have explored limited combinations of spatial and selective scenarios. Using simulations of secondary contact of populations, we explore barriers to gene flow in various selective and spatial contexts in continuous, two‐dimensional, spatially explicit environments. In general, the effects of hitchhiking are strongest in environments with regular spatial patterning of starkly divergent habitat types. When divergent selection is very strong, the absence of intermediate habitat types increases the effects of hitchhiking. However, when selection is moderate or weak, regular (vs. random) spatial arrangement of habitat types becomes more important than the presence of intermediate habitats per se. We also document counterintuitive processes arising from the stochastic interplay between selection, gene flow and drift. Our results indicate that generalization of results from two‐deme models requires caution and increase understanding of the genomic and geographic basis of population divergence.     

T-DNA insertional mutagenesis for functional genomics in rice

We have produced 22 090 primary transgenic rice plants that carry a T-DNA insertion, which has resulted in 18 358 fertile lines. Genomic DNA gel-blot and PCR analyses have shown that approximately 65% of the population contains more than one copy of the inserted T-DNA. Hygromycin resistance tests revealed that transgenic plants contain an average of 1.4 loci of T-DNA inserts. Therefore, it can be estimated that approximately 25 700 taggings have been generated. The binary vector used in the insertion contained the promoterless beta-glucuronidase (GUS) reporter gene with an intron and multiple splicing donors and acceptors immediately next to the right border. Therefore, this gene trap vector is able to detect a gene fusion between GUS and an endogenous gene, which is tagged by T-DNA. Histochemical GUS assays were carried out in the leaves and roots from 5353 lines, mature flowers from 7026 lines, and developing seeds from 1948 lines. The data revealed that 1.6-2.1% of tested organs were GUS-positive in the tested organs, and that their GUS expression patterns were organ- or tissue-specific or ubiquitous in all parts of the plant. The large population of T-DNA-tagged lines will be useful for identifying insertional mutants in various genes and for discovering new genes in rice.     

Genomic regions associated with antibody response to sheep red blood cells in the chicken

F(1) and F(2) populations were generated by crossing two lines of chickens divergently selected from a common founder population for 32 generations for either high or low antibody response 5 days post-injection of a non-pathogenic antigen, sheep red blood cells (SRBCs). The number of loci with major effects on day 5 SRBC titers was estimated to be more than 7 in this population. There was a significant association between MHC haplotype and day 5 antibody titers as well as body weight at sexual maturity. A significant difference between reciprocal F(2) crosses for both 5- and 12-day antibody titers suggests that sex chromosome and/or parent of origin effects on autosomal loci have an important role in immune response. A single marker-trait association analysis on 1024 genetic markers and 128 F(2) individuals detected 11 genomic regions associated with antibody response traits and 17 regions associated with body weight gain. Several of the genomic regions identified as being associated with antibody response have been described previously, while novel regions associated with antibody response were identified on chromosomes 11 and 24. Based on the lack of overlap of the regions associated with body weight and antibody response, we conclude that while these phenotypes are inversely correlated in the selected lines, they are controlled by distinct genetic loci and may be reflective of intense selection pressure on loci affecting the partitioning of nutrients between the immune system and growth pathways.     

Epigenetic heterogeneity at imprinted loci in normal populations

Genomic imprinting is the phenomenon by which the two alleles of certain genes are differentially expressed according to their parental origin. Extensive analysis of allelic expression at multiple imprinted loci in a normal population has not performed so far. In the present study, we examined the allelic expression pattern of three imprinted genes in a panel of 262 Japanese normal individuals. We observed differences in the extent of maintenance of allele-specific expression of the three genes. The allelic expression of small nuclear ribonucleoprotein N (SNRPN) was stringently regulated while that of multimembrane-spanning polyspecific transporter-like gene 1 (IMPT1) showed a large degree of variation. Significant biallelic expression of insulin-like growth factor II (IGF2) was observed in about 10% of normal individuals. Our findings add to the accumulating evidence for variable allelic expression at multiple loci in a normal human population. This epigenetic heterogeneity can be a stable trait and potentially influence individual phenotypes.