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.

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

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.

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.

An Ultra-Dense Haploid Genetic Map for Evaluating the Highly Fragmented Genome Assembly of Norway Spruce (Picea abies).

Norway spruce (Picea abies (L.) Karst.) is a conifer species of substanital economic and ecological importance. In common with most conifers, the P. abies genome is very large (∼20 Gbp) and contains a high fraction of repetitive DNA. The current P. abies genome assembly (v1.0) covers approximately 60% of the total genome size but is highly fragmented, consisting of >10 million scaffolds. The genome annotation contains 66,632 gene models that are at least partially validated (www.congenie.org), however, the fragmented nature of the assembly means that there is currently little information available on how these genes are physically distributed over the 12 P. abies chromosomes. By creating an ultra-dense genetic linkage map, we anchored and ordered scaffolds into linkage groups, which complements the fine-scale information available in assembly contigs. Our ultra-dense haploid consensus genetic map consists of 21,056 markers derived from 14,336 scaffolds that contain 17,079 gene models (25.6% of the validated gene models) that we have anchored to the 12 linkage groups. We used data from three independent component maps, as well as comparisons with previously published Picea maps to evaluate the accuracy and marker ordering of the linkage groups. We demonstrate that approximately 3.8% of the anchored scaffolds and 1.6% of the gene models covered by the consensus map have likely assembly errors as they contain genetic markers that map to different regions within or between linkage groups. We further evaluate the utility of the genetic map for the conifer research community by using an independent data set of unrelated individuals to assess genome-wide variation in genetic diversity using the genomic regions anchored to linkage groups. The results show that our map is sufficiently dense to enable detailed evolutionary analyses across the P. abies genome.

RPASE: Individual-based allele-specific expression detection without prior knowledge of haplotype phase.

Variation in gene expression is believed to make a significant contribution to phenotypic diversity and divergence. The analysis of allele-specific expression (ASE) can reveal important insights into gene expression regulation. We developed a novel method called RPASE (Read-backed Phasing-based ASE detection) to test for genes that show ASE. With mapped RNA-seq data from a single individual and a list of SNPs from the same individual as the only input, RPASE is capable of aggregating information across multiple dependent SNPs and producing individual-based gene-level tests for ASE. RPASE performs well in simulations and comparisons. We applied RPASE to multiple bird species and found a potentially rich landscape of ASE.

Tracking the NGS revolution: managing life science research on shared high-performance computing clusters.

Background: Next-generation sequencing (NGS) has transformed the life sciences, and many research groups are newly dependent upon computer clusters to store and analyze large datasets. This creates challenges for e-infrastructures accustomed to hosting computationally mature research in other sciences. Using data gathered from our own clusters at UPPMAX computing center at Uppsala University, Sweden, where core hour usage of ∼800 NGS and ∼200 non-NGS projects is now similar, we compare and contrast the growth, administrative burden, and cluster usage of NGS projects with projects from other sciences. Results: The number of NGS projects has grown rapidly since 2010, with growth driven by entry of new research groups. Storage used by NGS projects has grown more rapidly since 2013 and is now limited by disk capacity. NGS users submit nearly twice as many support tickets per user, and 11 more tools are installed each month for NGS projects than for non-NGS projects. We developed usage and efficiency metrics and show that computing jobs for NGS projects use more RAM than non-NGS projects, are more variable in core usage, and rarely span multiple nodes. NGS jobs use booked resources less efficiently for a variety of reasons. Active monitoring can improve this somewhat. Conclusions: Hosting NGS projects imposes a large administrative burden at UPPMAX due to large numbers of inexperienced users and diverse and rapidly evolving research areas. We provide a set of recommendations for e-infrastructures that host NGS research projects. We provide anonymized versions of our storage, job, and efficiency databases.

Bacteria colonising Penstemon digitalis show volatile and tissue-specific responses to a natural concentration range of the floral volatile linalool.

Bacteria on floral tissue can have negative effects by consuming resources and affecting nectar quality, which subsequently could reduce pollinator visitation and plant fitness. Plants however can employ chemical defences to reduce bacteria density. In North American, bee-pollinated Penstemon digitalis, the nectar volatile S-(+)-linalool can influence plant fitness, and terpenes such as linalool are known for their antimicrobial properties suggesting that it may also play a role in plant-microbe interactions. Therefore, we hypothesized linalool could affect bacterial growth on P. digitalis plants/flowers. Because P. digitalis emits linalool from nectar and nectary tissue but not petals, we hypothesised that the effects of linalool could depend on tissue of origin due to varying exposure. We isolated bacteria from nectary tissue, petals and leaves, and compared their growth relative to control using two volatile concentrations representing the natural emission range of linalool. To assess whether effects were specific to linalool, we compared results with the co-occurring nectar volatile, methyl nicotinate. We show that response to floral volatiles can be substance and tissue-origin specific. Because linalool could slow growth rate of bacteria across the P. digitalis phyllosphere, floral emission of linalool could play a role in mediating plant-bacteria interactions in this system.