Gancho Trifonu Slavov

Comparative genomics reveals insights into avian genome evolution and adaptation

Birds are the most species-rich class of tetrapod vertebrates and have wide relevance across many research fields. We explored bird macroevolution using full genomes from 48 avian species representing all major extant clades. The avian genome is principally characterized by its constrained size, which predominantly arose because of lineage-specific erosion of repetitive elements, large segmental deletions, and gene loss. Avian genomes furthermore show a remarkably high degree of evolutionary stasis at the levels of nucleotide sequence, gene synteny, and chromosomal structure. Despite this pattern of conservation, we detected many non-neutral evolutionary changes in protein-coding genes and noncoding regions. These analyses reveal that pan-avian genomic diversity covaries with adaptations to different lifestyles and convergent evolution of traits.

Population genomics of Populus trichocarpa identifies signatures of selection and adaptive trait associations

Forest trees are dominant components of terrestrial ecosystems that have global ecological and economic importance. Despite distributions that span wide environmental gradients, many tree populations are locally adapted, and mechanisms underlying this adaptation are poorly understood. Here we use a combination of whole-genome selection scans and association analyses of 544 Populus trichocarpa trees to reveal genomic bases of adaptive variation across a wide latitudinal range. Three hundred ninety-seven genomic regions showed evidence of recent positive and/or divergent selection and enrichment for associations with adaptive traits that also displayed patterns consistent with natural selection. These regions also provide unexpected insights into the evolutionary dynamics of duplicated genes and their roles in adaptive trait variation.

Genome resequencing reveals multiscale geographic structure and extensive linkage disequilibrium in the forest tree Populus trichocarpa

• Plant population genomics informs evolutionary biology, breeding, conservation and bioenergy feedstock development. For example, the detection of reliable phenotype-genotype associations and molecular signatures of selection requires a detailed knowledge about genome-wide patterns of allele frequency variation, linkage disequilibrium and recombination. • We resequenced 16 genomes of the model tree Populus trichocarpa and genotyped 120 trees from 10 subpopulations using 29,213 single-nucleotide polymorphisms. • Significant geographic differentiation was present at multiple spatial scales, and range-wide latitudinal allele frequency gradients were strikingly common across the genome. The decay of linkage disequilibrium with physical distance was slower than expected from previous studies in Populus, with r(2) dropping below 0.2 within 3-6 kb. Consistent with this, estimates of recent effective population size from linkage disequilibrium (N(e) ≈ 4000-6000) were remarkably low relative to the large census sizes of P. trichocarpa stands. Fine-scale rates of recombination varied widely across the genome, but were largely predictable on the basis of DNA sequence and methylation features. • Our results suggest that genetic drift has played a significant role in the recent evolutionary history of P. trichocarpa. Most importantly, the extensive linkage disequilibrium detected suggests that genome-wide association studies and genomic selection in undomesticated populations may be more feasible in Populus than previously assumed.

Contrasting patterns of evolution following whole genome versus tandem duplication events in Populus

Comparative analysis of multiple angiosperm genomes has implicated gene duplication in the expansion and diversification of many gene families. However, empirical data and theory suggest that whole-genome and small-scale duplication events differ with respect to the types of genes preserved as duplicate pairs. We compared gene duplicates resulting from a recent whole genome duplication to a set of tandemly duplicated genes in the model forest tree Populus trichocarpa. We used a combination of microarray expression analyses of a diverse set of tissues and functional annotation to assess factors related to the preservation of duplicate genes of both types. Whole genome duplicates are 700 bp longer and are expressed in 20% more tissues than tandem duplicates. Furthermore, certain functional categories are over-represented in each class of duplicates. In particular, disease resistance genes and receptor-like kinases commonly occur in tandem but are significantly under-retained following whole genome duplication, while whole genome duplicate pairs are enriched for members of signal transduction cascades and transcription factors. The shape of the distribution of expression divergence for duplicated pairs suggests that nearly half of the whole genome duplicates have diverged in expression by a random degeneration process. The remaining pairs have more conserved gene expression than expected by chance, consistent with a role for selection under the constraints of gene balance. We hypothesize that duplicate gene preservation in Populus is driven by a combination of subfunctionalization of duplicate pairs and purifying selection favoring retention of genes encoding proteins with large numbers of interactions.

Extensive pollen flow in two ecologically contrasting populations of Populus trichocarpa

Pollen‐mediated gene flow was measured in two populations of black cottonwood using direct (paternity analysis) and indirect (correlated paternity) methods. The Marchel site was an area with an approximate radius of 250 m in a large continuous stand growing in a mesic habitat in western Oregon. In contrast, the Vinson site was an area with a radius of approximately 10 km and consisted of small, disjunct and isolated stands in the high desert of eastern Oregon. Pollen immigration was extensive in both populations, and was higher in the Marchel site (0.54 ± 0.02) than in the substantially larger and more isolated Vinson site (0.32 ± 0.02). Pollen pool differentiation among mothers was approximately five times stronger in the Vinson population (ΦFT = 0.253, N = 27 mothers) than in the Marchel population (ΦFT = 0.052, N = 5 mothers). Pollen dispersal was modelled using a mixed dispersal curve that incorporated pollen immigration. Predicted pollination frequencies generated based on this curve were substantially more accurate than those based on the widely used exponential power dispersal curve. Male neighbourhood sizes (sensu Wright 1946 ) estimated using paternity analysis and pollen pool differentiation were remarkably similar. They were three to five times smaller in the Vinson population, which reflected the substantial ecological and demographic differences between the two populations. When the same mathematical function was used, applying direct and indirect methods resulted in similar pollen dispersal curves, thus confirming the value of indirect methods as a viable lower‐cost alternative to paternity analysis.

A 34K SNP genotyping array for Populus trichocarpa: design, application to the study of natural populations and transferability to other Populus species.

Genetic mapping of quantitative traits requires genotypic data for large numbers of markers in many individuals. For such studies, the use of large single nucleotide polymorphism (SNP) genotyping arrays still offers the most cost‐effective solution. Herein we report on the design and performance of a SNP genotyping array for Populus trichocarpa (black cottonwood). This genotyping array was designed with SNPs pre‐ascertained in 34 wild accessions covering most of the species latitudinal range. We adopted a candidate gene approach to the array design that resulted in the selection of 34 131 SNPs, the majority of which are located in, or within 2 kb of, 3543 candidate genes. A subset of the SNPs on the array (539) was selected based on patterns of variation among the SNP discovery accessions. We show that more than 95% of the loci produce high quality genotypes and that the genotyping error rate for these is likely below 2%. We demonstrate that even among small numbers of samples (n = 10) from local populations over 84% of loci are polymorphic. We also tested the applicability of the array to other species in the genus and found that the number of polymorphic loci decreases rapidly with genetic distance, with the largest numbers detected in other species in section Tacamahaca. Finally, we provide evidence for the utility of the array to address evolutionary questions such as intraspecific studies of genetic differentiation, species assignment and the detection of natural hybrids.

Genome-wide association studies and prediction of 17 traits related to phenology, biomass and cell wall composition in the energy grass Miscanthus sinensis

Increasing demands for food and energy require a step change in the effectiveness, speed and flexibility of crop breeding. Therefore, the aim of this study was to assess the potential of genome-wide association studies (GWASs) and genomic selection (i.e. phenotype prediction from a genome-wide set of markers) to guide fundamental plant science and to accelerate breeding in the energy grass Miscanthus. We generated over 100 000 single-nucleotide variants (SNVs) by sequencing restriction site-associated DNA (RAD) tags in 138 Micanthus sinensis genotypes, and related SNVs to phenotypic data for 17 traits measured in a field trial. Confounding by population structure and relatedness was severe in naïve GWAS analyses, but mixed-linear models robustly controlled for these effects and allowed us to detect multiple associations that reached genome-wide significance. Genome-wide prediction accuracies tended to be moderate to high (average of 0.57), but varied dramatically across traits. As expected, predictive abilities increased linearly with the size of the mapping population, but reached a plateau when the number of markers used for prediction exceeded 10 000–20 000, and tended to decline, but remain significant, when cross-validations were performed across subpopulations. Our results suggest that the immediate implementation of genomic selection in Miscanthus breeding programs may be feasible.

Salient Biological Features, Systematics, and Genetic Variation of Populus

The genus Populus includes morphologically diverse species of decid- uous, relatively short-lived, and fast-growing trees. Most species have wide ranges of distribution but tend to occur primarily in riparian or mountainous habitats. Trees from this genus are typically dioecious, flower before leaf emer- gence, and produce large amounts of wind-dispersed pollen or seeds. Seedlings are drought- and shade-intolerant, and their establishment depends on distur- bance and high soil moisture. Asexual reproduction is common and occurs via root sprouting and/or rooting of shoots. Fossil records suggest that the genus appeared in the late Paleocene or early Eocene (i.e., 50-60 million years BP). According to one commonly used classification, the genus is comprised of 29 species divided into six sections, but a number of phylogenetic inconsistencies remain. Natural hybridization both within and among sections is extensive and is believed to have played a major role in the evolution of extant species of Populus. Both neutral molecular markers and adaptive traits reveal high levels of genetic variation within populations. Deviations from Hardy-Weinberg equi- librium are commonly detected in molecular marker studies. These deviations typically have small to moderate magnitudes and tend to be caused by heterozy- gote deficiency, indicating the possible existence of population substructure. Genetic differentiation among populations is much stronger for adaptive traits than for neutral markers, which suggests that divergent selection has played a dominant role in shaping patterns of adaptive genetic variation. Molecular and bioinfor- matic resources are actively being developed for multiple species of Populus, which makes this genus an excellent system for studying tree genetics and genomics.

Gene flow and simulation of transgene dispersal from hybrid poplar plantations

Gene flow is a primary determinant of potential ecological impacts of transgenic trees. However, gene flow is a complex process that must be assessed in the context of realistic genetic, management, and environmental conditions. We measured gene flow from hybrid poplar plantations using morphological and genetic markers, and developed a spatially explicit landscape model to simulate pollination, dispersal, establishment, and mortality in the context of historical and projected disturbance and land-use regimes. Most pollination and seed establishment occurred within 450 m of the source, with a very long tail. Modeled transgene flow was highly context-dependent, strongly influenced by the competitive effects of transgenes, transgenic fertility, plantation rotation length, disturbance regime, and spatial and temporal variation in selection. The use of linked infertility genes even if imperfect, substantially reduced transgene flow in a wide range of modeled scenarios. The significance of seed and vegetative dispersal was highly dependent on plantation size. Our empirical and modeling studies suggest that transgene spread can be spatially extensive. However, the amount of spread is highly dependent on ecological and management context, and can be greatly limited or prevented by management or mitigation genes such as those that cause sexual infertility.

The obscure events contributing to the evolution of an incipient sex chromosome in Populus: a retrospective working hypothesis

Genetic determination of gender is a fundamental developmental and evolutionary process in plants. Although it appears that dioecy in Populus is genetically controlled, the precise gender-determining systems remain unclear. The recently released second draft assembly and annotated gene set of the Populus genome provided an opportunity to revisit this topic. We hypothesized that over evolutionary time, selective pressure has reformed the genome structure and gene composition in the peritelomeric region of the chromosome XIX, which has resulted in a distinctive genome structure and cluster of genes contributing to gender determination in Populus trichocarpa. Multiple lines of evidence support this working hypothesis. First, the peritelomeric region of the chromosome XIX contains significantly fewer single nucleotide polymorphisms than the rest of Populus genome and has a distinct evolutionary history. Second, the peritelomeric end of chromosome XIX contains the largest cluster of the nucleotide-binding site–leucine-rich repeat (NBS–LRR) class of disease resistance genes in the entire Populus genome. Third, there is a high occurrence of small microRNAs on chromosome XIX, which is coincident to the region containing the putative gender-determining locus and the major cluster of NBS–LRR genes. Further, by analyzing the metabolomic profiles of floral bud in male and female Populus trees using a gas chromatography-mass spectrometry, we found that there are gender-specific accumulations of phenolic glycosides. Taken together, these findings led to the hypothesis that resistance to and regulation of a floral pathogen and gender determination coevolved, and that these events triggered the emergence of a nascent sex chromosome. Further studies of chromosome XIX will provide new insights into the genetic control of gender determination in Populus.