Y-Linked Copy Number Polymorphism of Target of Rapamycin Is Associated with Sexual Size Dimorphism in Seed Beetles.

The Y chromosome is theorized to facilitate evolution of sexual dimorphism by accumulating sexually antagonistic loci, but empirical support is scarce. Due to the lack of recombination, Y chromosomes are prone to degenerative processes, which poses a constraint on their adaptive potential. Yet, in the seed beetle, Callosobruchus maculatus segregating Y linked variation affects male body size and thereby sexual size dimorphism (SSD). Here, we assemble C. maculatus sex chromosome sequences and identify molecular differences associated with Y-linked SSD variation. The assembled Y chromosome is largely euchromatic and contains over 400 genes, many of which are ampliconic with a mixed autosomal and X chromosome ancestry. Functional annotation suggests that the Y chromosome plays important roles in males beyond primary reproductive functions. Crucially, we find that, besides an autosomal copy of the gene target of rapamycin (TOR), males carry an additional TOR copy on the Y chromosome. TOR is a conserved regulator of growth across taxa, and our results suggest that a Y-linked TOR provides a male specific opportunity to alter body size. A comparison of Y haplotypes associated with male size difference uncovers a copy number variation for TOR, where the haplotype associated with decreased male size, and thereby increased sexual dimorphism, has two additional TOR copies. This suggests that sexual conflict over growth has been mitigated by autosome to Y translocation of TOR followed by gene duplications. Our results reveal that despite of suppressed recombination, the Y chromosome can harbor adaptive potential as a male-limited supergene.

Selection on Accessible Chromatin Regions in Capsella grandiflora.

Accurate estimates of genome-wide rates and fitness effects of new mutations are essential for an improved understanding of molecular evolutionary processes. Although eukaryotic genomes generally contain a large noncoding fraction, functional noncoding regions and fitness effects of mutations in such regions are still incompletely characterized. A promising approach to characterize functional noncoding regions relies on identifying accessible chromatin regions (ACRs) tightly associated with regulatory DNA. Here, we applied this approach to identify and estimate selection on ACRs in Capsella grandiflora, a crucifer species ideal for population genomic quantification of selection due to its favorable population demography. We describe a population-wide ACR distribution based on ATAC-seq data for leaf samples of 16 individuals from a natural population. We use population genomic methods to estimate fitness effects and proportions of positively selected fixations (α) in ACRs and find that intergenic ACRs harbor a considerable fraction of weakly deleterious new mutations, as well as a significantly higher proportion of strongly deleterious mutations than comparable inaccessible intergenic regions. ACRs are enriched for expression quantitative trait loci (eQTL) and depleted of transposable element insertions, as expected if intergenic ACRs are under selection because they harbor regulatory regions. By integrating empirical identification of intergenic ACRs with analyses of eQTL and population genomic analyses of selection, we demonstrate that intergenic regulatory regions are an important source of nearly neutral mutations. These results improve our understanding of selection on noncoding regions and the role of nearly neutral mutations for evolutionary processes in outcrossing Brassicaceae species.

Khoe-San Genomes Reveal Unique Variation and Confirm the Deepest Population Divergence in Homo sapiens.

The southern African indigenous Khoe-San populations harbor the most divergent lineages of all living peoples. Exploring their genomes is key to understanding deep human history. We sequenced 25 full genomes from five Khoe-San populations, revealing many novel variants, that 25% of variants are unique to the Khoe-San, and that the Khoe-San group harbors the greatest level of diversity across the globe. In line with previous studies, we found several gene regions with extreme values in genome-wide scans for selection, potentially caused by natural selection in the lineage leading to Homo sapiens and more recent in time. These gene regions included immunity-, sperm-, brain-, diet-, and muscle-related genes. When accounting for recent admixture, all Khoe-San groups display genetic diversity approaching the levels in other African groups and a reduction in effective population size starting around 100,000 years ago. Hence, all human groups show a reduction in effective population size commencing around the time of the Out-of-Africa migrations, which coincides with changes in the paleoclimate records, changes that potentially impacted all humans at the time.

Functional and evolutionary genomic inferences in Populus through genome and population sequencing of American and European aspen.

The Populus genus is one of the major plant model systems, but genomic resources have thus far primarily been available for poplar species, and primarily Populus trichocarpa (Torr. & Gray), which was the first tree with a whole-genome assembly. To further advance evolutionary and functional genomic analyses in Populus, we produced genome assemblies and population genetics resources of two aspen species, Populus tremula L. and Populus tremuloides Michx. The two aspen species have distributions spanning the Northern Hemisphere, where they are keystone species supporting a wide variety of dependent communities and produce a diverse array of secondary metabolites. Our analyses show that the two aspens share a similar genome structure and a highly conserved gene content with P. trichocarpa but display substantially higher levels of heterozygosity. Based on population resequencing data, we observed widespread positive and negative selection acting on both coding and noncoding regions. Furthermore, patterns of genetic diversity and molecular evolution in aspen are influenced by a number of features, such as expression level, coexpression network connectivity, and regulatory variation. To maximize the community utility of these resources, we have integrated all presented data within the PopGenIE web resource (PopGenIE.org).

Genome-wide association study identified novel candidate loci affecting wood formation in Norway spruce.

Norway spruce is a boreal forest tree species of significant ecological and economic importance. Hence there is a strong imperative to dissect the genetics underlying important wood quality traits in the species. We performed a functional genome-wide association study (GWAS) of 17 wood traits in Norway spruce using 178 101 single nucleotide polymorphisms (SNPs) generated from exome genotyping of 517 mother trees. The wood traits were defined using functional modelling of wood properties across annual growth rings. We applied a Least Absolute Shrinkage and Selection Operator (LASSO-based) association mapping method using a functional multilocus mapping approach that utilizes latent traits, with a stability selection probability method as the hypothesis testing approach to determine a significant quantitative trait locus. The analysis provided 52 significant SNPs from 39 candidate genes, including genes previously implicated in wood formation and tree growth in spruce and other species. Our study represents a multilocus GWAS for complex wood traits in Norway spruce. The results advance our understanding of the genetics influencing wood traits and identifies candidate genes for future functional studies.

Genomic analysis reveals major determinants of cis-regulatory variation in Capsella grandiflora.

Understanding the causes of cis-regulatory variation is a long-standing aim in evolutionary biology. Although cis-regulatory variation has long been considered important for adaptation, we still have a limited understanding of the selective importance and genomic determinants of standing cis-regulatory variation. To address these questions, we studied the prevalence, genomic determinants, and selective forces shaping cis-regulatory variation in the outcrossing plant Capsella grandiflora We first identified a set of 1,010 genes with common cis-regulatory variation using analyses of allele-specific expression (ASE). Population genomic analyses of whole-genome sequences from 32 individuals showed that genes with common cis-regulatory variation (i) are under weaker purifying selection and (ii) undergo less frequent positive selection than other genes. We further identified genomic determinants of cis-regulatory variation. Gene body methylation (gbM) was a major factor constraining cis-regulatory variation, whereas presence of nearby transposable elements (TEs) and tissue specificity of expression increased the odds of ASE. Our results suggest that most common cis-regulatory variation in C. grandiflora is under weak purifying selection, and that gene-specific functional constraints are more important for the maintenance of cis-regulatory variation than genome-scale variation in the intensity of selection. Our results agree with previous findings that suggest TE silencing affects nearby gene expression, and provide evidence for a link between gbM and cis-regulatory constraint, possibly reflecting greater dosage sensitivity of body-methylated genes. Given the extensive conservation of gbM in flowering plants, this suggests that gbM could be an important predictor of cis-regulatory variation in a wide range of plant species.

Haploid selection within a single ejaculate increases offspring fitness.

An inescapable consequence of sex in eukaryotes is the evolution of a biphasic life cycle with alternating diploid and haploid phases. The occurrence of selection during the haploid phase can have far-reaching consequences for fundamental evolutionary processes including the rate of adaptation, the extent of inbreeding depression, and the load of deleterious mutations, as well as for applied research into fertilization technology. Although haploid selection is well established in plants, current dogma assumes that in animals, intact fertile sperm within a single ejaculate are equivalent at siring viable offspring. Using the zebrafish Danio rerio, we show that selection on phenotypic variation among intact fertile sperm within an ejaculate affects offspring fitness. Longer-lived sperm sired embryos with increased survival and a reduced number of apoptotic cells, and adult male offspring exhibited higher fitness. The effect on embryo viability was carried over into the second generation without further selection and was equally strong in both sexes. Sperm pools selected by motile phenotypes differed genetically at numerous sites throughout the genome. Our findings clearly link within-ejaculate variation in sperm phenotype to offspring fitness and sperm genotype in a vertebrate and have major implications for adaptive evolution.

The Norway spruce genome sequence and conifer genome evolution.

Conifers have dominated forests for more than 200 million years and are of huge ecological and economic importance. Here we present the draft assembly of the 20-gigabase genome of Norway spruce (Picea abies), the first available for any gymnosperm. The number of well-supported genes (28,354) is similar to the >100 times smaller genome of Arabidopsis thaliana, and there is no evidence of a recent whole-genome duplication in the gymnosperm lineage. Instead, the large genome size seems to result from the slow and steady accumulation of a diverse set of long-terminal repeat transposable elements, possibly owing to the lack of an efficient elimination mechanism. Comparative sequencing of Pinus sylvestris, Abies sibirica, Juniperus communis, Taxus baccata and Gnetum gnemon reveals that the transposable element diversity is shared among extant conifers. Expression of 24-nucleotide small RNAs, previously implicated in transposable element silencing, is tissue-specific and much lower than in other plants. We further identify numerous long (>10,000 base pairs) introns, gene-like fragments, uncharacterized long non-coding RNAs and short RNAs. This opens up new genomic avenues for conifer forestry and breeding.

Variation in Linked Selection and Recombination Drive Genomic Divergence during Allopatric Speciation of European and American Aspens.

Despite the global economic and ecological importance of forest trees, the genomic basis of differential adaptation and speciation in tree species is still poorly understood. Populus tremula and Populus tremuloides are two of the most widespread tree species in the Northern Hemisphere. Using whole-genome re-sequencing data of 24 P. tremula and 22 P. tremuloides individuals, we find that the two species diverged ∼2.2-3.1 million years ago, coinciding with the severing of the Bering land bridge and the onset of dramatic climatic oscillations during the Pleistocene. Both species have experienced substantial population expansions following long-term declines after species divergence. We detect widespread and heterogeneous genomic differentiation between species, and in accordance with the expectation of allopatric speciation, coalescent simulations suggest that neutral evolutionary processes can account for most of the observed patterns of genetic differentiation. However, there is an excess of regions exhibiting extreme differentiation relative to those expected under demographic simulations, which is indicative of the action of natural selection. Overall genetic differentiation is negatively associated with recombination rate in both species, providing strong support for a role of linked selection in generating the heterogeneous genomic landscape of differentiation between species. Finally, we identify a number of candidate regions and genes that may have been subject to positive and/or balancing selection during the speciation process.

Cis-Regulatory Changes Associated with a Recent Mating System Shift and Floral Adaptation in Capsella.

The selfing syndrome constitutes a suite of floral and reproductive trait changes that have evolved repeatedly across many evolutionary lineages in response to the shift to selfing. Convergent evolution of the selfing syndrome suggests that these changes are adaptive, yet our understanding of the detailed molecular genetic basis of the selfing syndrome remains limited. Here, we investigate the role of cis-regulatory changes during the recent evolution of the selfing syndrome in Capsella rubella, which split from the outcrosser Capsella grandiflora less than 200 ka. We assess allele-specific expression (ASE) in leaves and flower buds at a total of 18,452 genes in three interspecific F1 C. grandiflora x C. rubella hybrids. Using a hierarchical Bayesian approach that accounts for technical variation using genomic reads, we find evidence for extensive cis-regulatory changes. On average, 44% of the assayed genes show evidence of ASE; however, only 6% show strong allelic expression biases. Flower buds, but not leaves, show an enrichment of cis-regulatory changes in genomic regions responsible for floral and reproductive trait divergence between C. rubella and C. grandiflora. We further detected an excess of heterozygous transposable element (TE) insertions near genes with ASE, and TE insertions targeted by uniquely mapping 24-nt small RNAs were associated with reduced expression of nearby genes. Our results suggest that cis-regulatory changes have been important during the recent adaptive floral evolution in Capsella and that differences in TE dynamics between selfing and outcrossing species could be important for rapid regulatory divergence in association with mating system shifts.

A chromosome-level assembly of the seed beetle Callosobruchus maculatus genome with annotation of its repetitive elements.

Callosobruchus maculatus is a major agricultural pest of legume crops worldwide and an established model system in ecology and evolution. Yet, current molecular biological resources for this species are limited. Here, we employ Hi-C sequencing to generate a greatly improved genome assembly and we annotate its repetitive elements in a dedicated in-depth effort where we manually curate and classify the most abundant unclassified repeat subfamilies. We present a scaffolded chromosome-level assembly, which is 1.01 Gb in total length with 86% being contained within the 9 autosomes and the X chromosome. Repetitive sequences accounted for 70% of the total assembly. DNA transposons covered 18% of the genome, with the most abundant superfamily being Tc1-Mariner (9.75% of the genome). This new chromosome-level genome assembly of C. maculatus will enable future genetic and evolutionary studies not only of this important species but of beetles more generally.

Calling Structural Variants with Confidence from Short-Read Data in Wild Bird Populations.

Comprehensive characterization of structural variation in natural populations has only become feasible in the last decade. To investigate the population genomic nature of structural variation, reproducible and high-confidence structural variation callsets are first required. We created a population-scale reference of the genome-wide landscape of structural variation across 33 Nordic house sparrows (Passer domesticus). To produce a consensus callset across all samples using short-read data, we compare heuristic-based quality filtering and visual curation (Samplot/PlotCritic and Samplot-ML) approaches. We demonstrate that curation of structural variants is important for reducing putative false positives and that the time invested in this step outweighs the potential costs of analyzing short-read-discovered structural variation data sets that include many potential false positives. We find that even a lenient manual curation strategy (e.g. applied by a single curator) can reduce the proportion of putative false positives by up to 80%, thus enriching the proportion of high-confidence variants. Crucially, in applying a lenient manual curation strategy with a single curator, nearly all (>99%) variants rejected as putative false positives were also classified as such by a more stringent curation strategy using three additional curators. Furthermore, variants rejected by manual curation failed to reflect the expected population structure from SNPs, whereas variants passing curation did. Combining heuristic-based quality filtering with rapid manual curation of structural variants in short-read data can therefore become a time- and cost-effective first step for functional and population genomic studies requiring high-confidence structural variation callsets.

Germline mutation rate is elevated in young and old parents in Caenorhabditis remanei.

The effect of parental age on germline mutation rate across generations is not fully understood. While some studies report a positive linear relationship of mutation rate with increasing age, others suggest that mutation rate varies with age but not in a linear fashion. We investigated the effect of parental age on germline mutations by generating replicated mutation accumulation lines in Caenorhabditis remanei at three parental ages ("Young T1" [Day 1], "Peak T2" [Day 2], and "Old T5" [Day 5] parents). We conducted whole-genome resequencing and variant calling to compare differences in mutation rates after three generations of mutation accumulation. We found that Peak T2 lines had an overall reduced mutation rate compared to Young T1 and Old T5 lines, but this pattern of the effect varied depending on the variant impact. Specifically, we found no high-impact variants in Peak T2 lines, and modifiers and up- and downstream gene variants were less frequent in these lines. These results suggest that animals at the peak of reproduction have better DNA maintenance and repair compared to young and old animals. We propose that C. remanei start to reproduce before they optimize their DNA maintenance and repair, trading the benefits of earlier onset of reproduction against offspring mutation load. The increase in offspring mutation load with age likely represents germline senescence.

Use of alpha, beta, and gamma diversity measures to characterize seed dispersal by animals.

Seed dispersal shapes ecological and evolutionary dynamics of plant populations. Here, we extend classical diversity measures to study the impact of disperser behavior on seed dispersal. We begin by extending our previous diversity structure approach, which partitioned seed source diversity within and among dispersal sites, into the more general framework of traditional diversity measures. This statistical approach allows an assessment of the extent to which foraging behavior shapes α and γ diversity, as well as the divergence in seed sources among dispersal sites, which we call δ. We also introduce tests to facilitate comparisons of diversity among dispersal sites and seed vectors and to compare overall diversity among sampled systems. We then apply these tools to investigate the diversity blend of parentage resulting from seed dispersal by two avian seed vectors with very different social and foraging behaviors: (1) acorn woodpeckers, transporting Quercus agrifolia acorns, and (2) long-wattled umbrellabirds, transporting Oenocarpus bataua palm nuts. Using these diversity and divergence measures, we test the hypothesis that different foraging behaviors generate distinctive diversity partitions for the two focal tree species. This approach provides a new tool for assessment of the impact of dispersal agents on the seed source structure of plant populations, which can be extended to include the impact of virtually any propagule vector for a range of systems.

Using seedling and pericarp tissues to determine maternal parentage of dispersed valley oak recruits.

The spatial pattern of established seedlings yields valuable information about variation in fecundity, dispersal, and spatial structure of distributed recruits, but separating maternal and paternal contributions in monoecious species has been hampered by the "2 parent" problem. It is now possible to determine the maternal parentage of established recruits with genetic assay of maternally derived tissues of the seed or fruit, but the DNA of weathered maternal tissues often yields unreliable genotypes, reducing the practical range of such assay. We develop a mixed assay of seedling and seed (pericarp) tissues and illustrate it with distributed recruits of California valley oak (Quercus lobata Née). Detailed analysis indicates correct maternal assignment rates of canopy patch recruits of 56% (seedling assay only) versus 94% (mixed assay). For open patch recruits, maternal assignment rates were less than 50% (seedling assay only) versus 91% (mixed assay). The strategy of choice is to use seedling genotypes to identify a small set of credible parental candidates and then deploy 3-4 well-chosen pericarp/endocarp loci to reduce that list to a single obvious maternal candidate. The increase in the number of recruits available for subsequent analysis is pronounced, increasing precision and statistical power for subsequent inference.

Foraging patterns of acorn woodpeckers (Melanerpes formicivorus) on valley oak (Quercus lobata Née) in two California oak savanna-woodlands.

Landscape characteristics and social behavior can affect the foraging patterns of seed-dependent animals. We examine the movement of acorns from valley oak (Quercus lobata) trees to granaries maintained by acorn woodpeckers (Melanerpes formicivorus) in two California oak savanna-woodlands differing in the distribution of Q. lobata within each site. In 2004, we sampled Q. lobata acorns from 16 granaries at Sedgwick Reserve in Santa Barbara County and 18 granaries at Hastings Reserve in Monterey County. Sedgwick has lower site-wide density of Q. lobata than Hastings as well as different frequencies of other Quercus species common to both sites. We found acorn woodpeckers foraged from fewer Q. lobata seed source trees (K(g) = 4.1 ± 0.5) at Sedgwick than at Hastings (K(g) = 7.6 ± 0.6) and from fewer effective seed sources (N(em)* = 2.00 and 5.78, respectively). The differences between sites are due to a greater number of incidental seed sources used per granary at Hastings than at Sedgwick. We also found very low levels of seed source sharing between adjacent granaries, indicating that territoriality is strong at both sites and that each social group forages on its own subset of trees. We discovered an interesting spatial pattern in the location of granaries. At Sedgwick, acorn woodpeckers situated their granaries within areas of higher-than-average tree density, while at Hastings, they placed them within areas of lower-than-average tree density, with the outcome that granaries at the two sites were located in areas of similar valley oak density. Our results illustrate that landscape characteristics might influence the number of trees visited by acorn woodpeckers and the locations of territories, while woodpecker social behavior, such as territoriality, shapes which trees are visited and whether they are shared with other social groups.

Evolutionary diversification of the Sm family of RNA-associated proteins.

The Sm family of proteins is closely associated with RNA metabolism throughout all life. These proteins form homomorphic and heteromorphic rings consisting of six or seven subunits with a characteristic central pore, the presence of which is critical for binding U-rich regions of single-stranded RNA. Eubacteria and Archaea typically carry one or two forms of Sm proteins and assemble one homomorphic ring per Sm protein. Eukaryotes typically carry 16 or more Sm proteins that assemble to form heteromorphic rings which lie at the center of a number of critical RNA-associated small nuclear ribonucleoproteins (snRNPs). High Sm protein diversity and heteromorphic Sm rings are features stretching back to the origin of eukaryotes; very deep phylogenetic divisions among existing Sm proteins indicate simultaneous evolution across essentially all existing eukaryotic life. Two basic forms of heteromorphic Sm rings are found in eukaryotes. Fixed Sm rings are highly stable and static and are assembled around an RNA cofactor. Flexible Sm rings also stabilize and chaperone RNA but assemble in the absence of an RNA substrate and, more significantly, associate with and dissociate from RNA substrates more freely than fixed rings. This suggests that the conformation of flexible Sm rings might be modified in some specific manner to facilitate association and dissociation with RNA. Diversification of eukaryotic Sm proteins may have been initiated by gene transfers and/or genome clashes that accompanied the origin of the eukaryotic cell itself, with further diversification driven by a greater need for steric specificity within increasingly complex snRNPs.

Intron presence-absence polymorphisms in Daphnia.

Here, we report 2 novel intron gains segregating in populations of Daphnia pulex endemic to Oregon. These novel introns do not have an obvious source and are not present in any D. pulex populations outside Oregon, other species of Daphnia that we examined, or any other organism for which sequence data are available. Furthermore, the novel introns are both found in the same gene, a Rab GTPase (rab4), and they appear to differ in their insertion site by one base pair, providing some support to the proto-splice site hypothesis. The rarity of intron-gain polymorphisms is questioned as we discovered 2 events in an initial survey of only 6 nuclear loci in 36 Daphnia individuals. Neutrality tests failed to ascertain a clear selective effect for either intron insertion, and a significant difference in recombination rate was not observed in alleles that contain the novel intron insertion versus alleles lacking it. We conclude that one novel intron insertion segregating at high frequencies in Daphnia populations in Oregon is unlikely to be adaptive and may result from the reduced efficacy of selection in isolated populations of small effective size.

Mutation accumulation in real branches: fitness assays for genomic deleterious mutation rate and effect in large-statured plants.

The genomic deleterious mutation rate and mean effect are central to the biology and evolution of all species. Large-statured plants, such as trees, are predicted to have high mutation rates due to mitotic mutation and the absence of a sheltered germ line, but their size and generation time has hindered genetic study. We develop and test approaches for estimating deleterious mutation rates and effects from viability comparisons within the canopy of large-statured plants. Our methods, inspired by E. J. Klekowski, are a modification of the classic Bateman-Mukai mutation-accumulation experiment. Within a canopy, cell lineages accumulate mitotic mutations independently. Gametes or zygotes produced at more distal points by these cell lineages contain more mitotic mutations than those at basal locations, and within-flower selfs contain more homozygous mutations than between-flower selfs. The resulting viability differences allow demonstration of lethal mutation with experiments similar in size to assays of genetic load and allow estimates of the rate and effect of new mutations with moderate precision and bias similar to that of classic mutation-accumulation experiments in small-statured organisms. These methods open up new possibilities with the potential to provide valuable new insights into the evolutionary genetics of plants.

Endogenous mechanisms for the origins of spliceosomal introns.

Over 30 years since their discovery, the origin of spliceosomal introns remains uncertain. One nearly universally accepted hypothesis maintains that spliceosomal introns originated from self-splicing group-II introns that invaded the uninterrupted genes of the last eukaryotic common ancestor (LECA) and proliferated by "insertion" events. Although this is a possible explanation for the original presence of introns and splicing machinery, the emphasis on a high number of insertion events in the genome of the LECA neglects a considerable body of empirical evidence showing that spliceosomal introns can simply arise from coding or, more generally, nonintronic sequences within genes. After presenting a concise overview of some of the most common hypotheses and mechanisms for intron origin, we propose two further hypotheses that are broadly based on central cellular processes: 1) internal gene duplication and 2) the response to aberrant and fortuitously spliced transcripts. These two nonmutually exclusive hypotheses provide a powerful way to explain the establishment of spliceosomal introns in eukaryotes without invoking an exogenous source.