Douglas G. Scofield

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.

Natural Selection and Recombination Rate Variation Shape Nucleotide Polymorphism Across the Genomes of Three Related Populus Species.

A central aim of evolutionary genomics is to identify the relative roles that various evolutionary forces have played in generating and shaping genetic variation within and among species. Here we use whole-genome resequencing data to characterize and compare genome-wide patterns of nucleotide polymorphism, site frequency spectrum, and population-scaled recombination rates in three species of Populus: Populus tremula, P. tremuloides, and P. trichocarpa. We find that P. tremuloides has the highest level of genome-wide variation, skewed allele frequencies, and population-scaled recombination rates, whereas P. trichocarpa harbors the lowest. Our findings highlight multiple lines of evidence suggesting that natural selection, due to both purifying and positive selection, has widely shaped patterns of nucleotide polymorphism at linked neutral sites in all three species. Differences in effective population sizes and rates of recombination largely explain the disparate magnitudes and signatures of linked selection that we observe among species. The present work provides the first phylogenetic comparative study on a genome-wide scale in forest trees. This information will also improve our ability to understand how various evolutionary forces have interacted to influence genome evolution among related species.

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.

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

Significance Despite long-standing interest in the contribution of cis-regulatory changes to adaptation, we still have a limited understanding of the selective importance and genomic determinants of cis-regulatory variation in natural populations. We use a combination of analyses of allele-specific expression and population genomic analyses to investigate the selective forces and genomic determinants of cis-regulatory variation in the outcrossing plant species Capsella grandiflora. We conclude that gene-specific functional constraints shape cis-regulatory variation and that genes with cis-regulatory variation are under relaxed purifying selection compared with other genes. Finally, we identify a link between gene body methylation and the extent of cis-regulatory constraint in natural populations. 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.

The Ty1-copia LTR retroelement family PARTC is highly conserved in conifers over 200 MY of evolution

Long Terminal Repeat retroelements (LTR-RTs) are a major component of many plant genomes. Although well studied and described in angiosperms, their features and dynamics are poorly understood in gymnosperms. Representative complete copies of a Ty1-copia element isolate in Picea abies and named PARTC were identified in six other conifer species (Picea glauca, Pinus sylvestris, Pinus taeda, Abies sibirica, Taxus baccata and Juniperus communis) covering more than 200 million years of evolution. Here we characterized the structure of this element, assessed its abundance across conifers, studied the modes and timing of its amplification, and evaluated the degree of conservation of its extant copies at nucleotide level over distant species. We demonstrated that the element is ancient, abundant, widespread and its paralogous copies are present in the genera Picea, Pinus and Abies as an LTR-RT family. The amplification leading to the extant copies of PARTC occurred over long evolutionary times spanning 10s of MY and mostly took place after the speciation of the conifers analyzed. The level of conservation of PARTC is striking and may be explained by low substitution rates and limited removal mechanisms for LTR-RTs. These PARTC features and dynamics are representative of a more general scenario for LTR-RTs in gymnosperms quite different from that characterizing the vast majority of LTR-RT elements in angiosperms.

Haploid selection within a single ejaculate increases offspring fitness

Significance Diploid organisms produce haploid gametes for sexual reproduction, resulting in a biphasic life cycle. Although selection during the diploid phase is well understood, selection during the haploid gametic stage and its consequences are largely ignored despite its potential importance for fundamental evolutionary processes, including the rate of adaptation and inbreeding depression, as well as for applied research into fertilization technology. A current dogma assumes that in animals selection on the haploid gametic genotype is minimal. We examined the importance of haploid selection in the zebrafish and found strong fitness consequences of selection on sperm phenotype in the resulting offspring. Genomic data support the idea that these effects may well be the consequence of selection on the haploid sperm genotype. 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.

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.

LTR Retrotransposons Show Low Levels of Unequal Recombination and High Rates of Intraelement Gene Conversion in Large Plant Genomes

Abstract The accumulation and removal of transposable elements (TEs) is a major driver of genome size evolution in eukaryotes. In plants, long terminal repeat (LTR) retrotransposons (LTR-RTs) represent the majority of TEs and form most of the nuclear DNA in large genomes. Unequal recombination (UR) between LTRs leads to removal of intervening sequence and formation of solo-LTRs. UR is a major mechanism of LTR-RT removal in many angiosperms, but our understanding of LTR-RT-associated recombination within the large, LTR-RT-rich genomes of conifers is quite limited. We employ a novel read-based methodology to estimate the relative rates of LTR-RT-associated UR within the genomes of four conifer and seven angiosperm species. We found the lowest rates of UR in the largest genomes studied, conifers and the angiosperm maize. Recombination may also resolve as gene conversion, which does not remove sequence, so we analyzed LTR-RT-associated gene conversion events (GCEs) in Norway spruce and six angiosperms. Opposite the trend for UR, we found the highest rates of GCEs in Norway spruce and maize. Unlike previous work in angiosperms, we found no evidence that rates of UR correlate with retroelement structural features in the conifers, suggesting that another process is suppressing UR in these species. Recent results from diverse eukaryotes indicate that heterochromatin affects the resolution of recombination, by favoring gene conversion over crossing-over, similar to our observation of opposed rates of UR and GCEs. Control of LTR-RT proliferation via formation of heterochromatin would be a likely step toward large genomes in eukaryotes carrying high LTR-RT content.

Variant Calling Using NGS Data in European Aspen (Populus tremula)

Analysis of next-generation sequencing (NGS) data is rapidly becoming an important source of information for genetics and genomics studies. The utility of such data does, however, rely crucially on the accuracy and quality of SNP and genotype calling. Identification of genetic variants (SNPs and short indels) from NGS data is an area of active research and many recent statistical methods have been developed to both improve and quantify the large uncertainty associated with genotype calling. The detection of genetic variants from NGS data is prone to errors, due to multiple factors such as base-calling, alignment errors, and read coverage. Here we highlight some of the issues and review and exemplify some of the recent methods that have been developed for genotype calling. We also provide guidelines for their application to whole-genome re-sequencing data using a data set based on a number of European aspen (Populus tremula) individuals each sequenced to a depth of about 20× coverage per individual.