Vocal learning-associated convergent evolution in mammalian proteins and regulatory elements

Vocal production learning ("vocal learning") is a convergently evolved trait in vertebrates. To identify brain genomic elements associated with mammalian vocal learning, we integrated genomic, anatomical, and neurophysiological data from the Egyptian fruit bat (Rousettus aegyptiacus) with analyses of the genomes of 215 placental mammals. First, we identified a set of proteins evolving more slowly in vocal learners. Then, we discovered a vocal motor cortical region in the Egyptian fruit bat, an emergent vocal learner, and leveraged that knowledge to identify active cis-regulatory elements in the motor cortex of vocal learners. Machine learning methods applied to motor cortex open chromatin revealed 50 enhancers robustly associated with vocal learning whose activity tended to be lower in vocal learners. Our research implicates convergent losses of motor cortex regulatory elements in mammalian vocal learning evolution.

Identification of constrained sequence elements across 239 primate genomes

Noncoding DNA is central to our understanding of human gene regulation and complex diseases1,2, and measuring the evolutionary sequence constraint can establish the functional relevance of putative regulatory elements in the human genome3-9. Identifying the genomic elements that have become constrained specifically in primates has been hampered by the faster evolution of noncoding DNA compared to protein-coding DNA10, the relatively short timescales separating primate species11, and the previously limited availability of whole-genome sequences12. Here we construct a whole-genome alignment of 239 species, representing nearly half of all extant species in the primate order. Using this resource, we identified human regulatory elements that are under selective constraint across primates and other mammals at a 5% false discovery rate. We detected 111,318 DNase I hypersensitivity sites and 267,410 transcription factor binding sites that are constrained specifically in primates but not across other placental mammals and validate their cis-regulatory effects on gene expression. These regulatory elements are enriched for human genetic variants that affect gene expression and complex traits and diseases. Our results highlight the important role of recent evolution in regulatory sequence elements differentiating primates, including humans, from other placental mammals.

The Genomic Basis of Adaptation to High Elevations in Africanized Honey Bees

A range of different genetic architectures underpin local adaptation in nature. Honey bees (Apis mellifera) in the Eastern African Mountains harbor high frequencies of two chromosomal inversions that likely govern adaptation to this high-elevation habitat. In the Americas, honey bees are hybrids of European and African ancestries and adaptation to latitudinal variation in climate correlates with the proportion of these ancestries across the genome. It is unknown which, if either, of these forms of genetic variation governs adaptation in honey bees living at high elevations in the Americas. Here, we performed whole-genome sequencing of 29 honey bees from both high- and low-elevation populations in Colombia. Analysis of genetic ancestry indicated that both populations were predominantly of African ancestry, but the East African inversions were not detected. However, individuals in the higher elevation population had significantly higher proportions of European ancestry, likely reflecting local adaptation. Several genomic regions exhibited particularly high differentiation between highland and lowland bees, containing candidate loci for local adaptation. Genes that were highly differentiated between highland and lowland populations were enriched for functions related to reproduction and sperm competition. Furthermore, variation in levels of European ancestry across the genome was correlated between populations of honey bees in the highland population and populations at higher latitudes in South America. The results are consistent with the hypothesis that adaptation to both latitude and elevation in these hybrid honey bees are mediated by variation in ancestry at many loci across the genome.

Genome sequencing of 2000 canids by the Dog10K consortium advances the understanding of demography, genome function and architecture

BACKGROUND

The international Dog10K project aims to sequence and analyze several thousand canine genomes. Incorporating 20 × data from 1987 individuals, including 1611 dogs (321 breeds), 309 village dogs, 63 wolves, and four coyotes, we identify genomic variation across the canid family, setting the stage for detailed studies of domestication, behavior, morphology, disease susceptibility, and genome architecture and function.

RESULTS

We report the analysis of > 48 M single-nucleotide, indel, and structural variants spanning the autosomes, X chromosome, and mitochondria. We discover more than 75% of variation for 239 sampled breeds. Allele sharing analysis indicates that 94.9% of breeds form monophyletic clusters and 25 major clades. German Shepherd Dogs and related breeds show the highest allele sharing with independent breeds from multiple clades. On average, each breed dog differs from the UU_Cfam_GSD_1.0 reference at 26,960 deletions and 14,034 insertions greater than 50 bp, with wolves having 14% more variants. Discovered variants include retrogene insertions from 926 parent genes. To aid functional prioritization, single-nucleotide variants were annotated with SnpEff and Zoonomia phyloP constraint scores. Constrained positions were negatively correlated with allele frequency. Finally, the utility of the Dog10K data as an imputation reference panel is assessed, generating high-confidence calls across varied genotyping platform densities including for breeds not included in the Dog10K collection.

CONCLUSIONS

We have developed a dense dataset of 1987 sequenced canids that reveals patterns of allele sharing, identifies likely functional variants, informs breed structure, and enables accurate imputation. Dog10K data are publicly available.

Evolutionary constraint and innovation across hundreds of placental mammals

Zoonomia is the largest comparative genomics resource for mammals produced to date. By aligning genomes for 240 species, we identify bases that, when mutated, are likely to affect fitness and alter disease risk. At least 332 million bases (~10.7%) in the human genome are unusually conserved across species (evolutionarily constrained) relative to neutrally evolving repeats, and 4552 ultraconserved elements are nearly perfectly conserved. Of 101 million significantly constrained single bases, 80% are outside protein-coding exons and half have no functional annotations in the Encyclopedia of DNA Elements (ENCODE) resource. Changes in genes and regulatory elements are associated with exceptional mammalian traits, such as hibernation, that could inform therapeutic development. Earth's vast and imperiled biodiversity offers distinctive power for identifying genetic variants that affect genome function and organismal phenotypes.

Relating enhancer genetic variation across mammals to complex phenotypes using machine learning

Protein-coding differences between species often fail to explain phenotypic diversity, suggesting the involvement of genomic elements that regulate gene expression such as enhancers. Identifying associations between enhancers and phenotypes is challenging because enhancer activity can be tissue-dependent and functionally conserved despite low sequence conservation. We developed the Tissue-Aware Conservation Inference Toolkit (TACIT) to associate candidate enhancers with species' phenotypes using predictions from machine learning models trained on specific tissues. Applying TACIT to associate motor cortex and parvalbumin-positive interneuron enhancers with neurological phenotypes revealed dozens of enhancer-phenotype associations, including brain size-associated enhancers that interact with genes implicated in microcephaly or macrocephaly. TACIT provides a foundation for identifying enhancers associated with the evolution of any convergently evolved phenotype in any large group of species with aligned genomes.

Evolutionary constraint and innovation across hundreds of placental mammals

Zoonomia is the largest comparative genomics resource for mammals produced to date. By aligning genomes for 240 species, we identify bases that, when mutated, are likely to affect fitness and alter disease risk. At least 332 million bases (~10.7%) in the human genome are unusually conserved across species (evolutionarily constrained) relative to neutrally evolving repeats, and 4552 ultraconserved elements are nearly perfectly conserved. Of 101 million significantly constrained single bases, 80% are outside protein-coding exons and half have no functional annotations in the Encyclopedia of DNA Elements (ENCODE) resource. Changes in genes and regulatory elements are associated with exceptional mammalian traits, such as hibernation, that could inform therapeutic development. Earth's vast and imperiled biodiversity offers distinctive power for identifying genetic variants that affect genome function and organismal phenotypes.

Integrating gene annotation with orthology inference at scale

Annotating coding genes and inferring orthologs are two classical challenges in genomics and evolutionary biology that have traditionally been approached separately, limiting scalability. We present TOGA (Tool to infer Orthologs from Genome Alignments), a method that integrates structural gene annotation and orthology inference. TOGA implements a different paradigm to infer orthologous loci, improves ortholog detection and annotation of conserved genes compared with state-of-the-art methods, and handles even highly fragmented assemblies. TOGA scales to hundreds of genomes, which we demonstrate by applying it to 488 placental mammal and 501 bird assemblies, creating the largest comparative gene resources so far. Additionally, TOGA detects gene losses, enables selection screens, and automatically provides a superior measure of mammalian genome quality. TOGA is a powerful and scalable method to annotate and compare genes in the genomic era.

Leveraging base-pair mammalian constraint to understand genetic variation and human disease

Thousands of genomic regions have been associated with heritable human diseases, but attempts to elucidate biological mechanisms are impeded by an inability to discern which genomic positions are functionally important. Evolutionary constraint is a powerful predictor of function, agnostic to cell type or disease mechanism. Single-base phyloP scores from 240 mammals identified 3.3% of the human genome as significantly constrained and likely functional. We compared phyloP scores to genome annotation, association studies, copy-number variation, clinical genetics findings, and cancer data. Constrained positions are enriched for variants that explain common disease heritability more than other functional annotations. Our results improve variant annotation but also highlight that the regulatory landscape of the human genome still needs to be further explored and linked to disease.

The contribution of historical processes to contemporary extinction risk in placental mammals

Species persistence can be influenced by the amount, type, and distribution of diversity across the genome, suggesting a potential relationship between historical demography and resilience. In this study, we surveyed genetic variation across single genomes of 240 mammals that compose the Zoonomia alignment to evaluate how historical effective population size (N_e) affects heterozygosity and deleterious genetic load and how these factors may contribute to extinction risk. We find that species with smaller historical N_e carry a proportionally larger burden of deleterious alleles owing to long-term accumulation and fixation of genetic load and have a higher risk of extinction. This suggests that historical demography can inform contemporary resilience. Models that included genomic data were predictive of species' conservation status, suggesting that, in the absence of adequate census or ecological data, genomic information may provide an initial risk assessment.

Mammalian evolution of human cis-regulatory elements and transcription factor binding sites

Understanding the regulatory landscape of the human genome is a long-standing objective of modern biology. Using the reference-free alignment across 241 mammalian genomes produced by the Zoonomia Consortium, we charted evolutionary trajectories for 0.92 million human candidate cis-regulatory elements (cCREs) and 15.6 million human transcription factor binding sites (TFBSs). We identified 439,461 cCREs and 2,024,062 TFBSs under evolutionary constraint. Genes near constrained elements perform fundamental cellular processes, whereas genes near primate-specific elements are involved in environmental interaction, including odor perception and immune response. About 20% of TFBSs are transposable element-derived and exhibit intricate patterns of gains and losses during primate evolution whereas sequence variants associated with complex traits are enriched in constrained TFBSs. Our annotations illuminate the regulatory functions of the human genome.

A genomic timescale for placental mammal evolution

The precise pattern and timing of speciation events that gave rise to all living placental mammals remain controversial. We provide a comprehensive phylogenetic analysis of genetic variation across an alignment of 241 placental mammal genome assemblies, addressing prior concerns regarding limited genomic sampling across species. We compared neutral genome-wide phylogenomic signals using concatenation and coalescent-based approaches, interrogated phylogenetic variation across chromosomes, and analyzed extensive catalogs of structural variants. Interordinal relationships exhibit relatively low rates of phylogenomic conflict across diverse datasets and analytical methods. Conversely, X-chromosome versus autosome conflicts characterize multiple independent clades that radiated during the Cenozoic. Genomic time trees reveal an accumulation of cladogenic events before and immediately after the Cretaceous-Paleogene (K-Pg) boundary, implying important roles for Cretaceous continental vicariance and the K-Pg extinction in the placental radiation.

The functional and evolutionary impacts of human-specific deletions in conserved elements

Conserved genomic sequences disrupted in humans may underlie uniquely human phenotypic traits. We identified and characterized 10,032 human-specific conserved deletions (hCONDELs). These short (average 2.56 base pairs) deletions are enriched for human brain functions across genetic, epigenomic, and transcriptomic datasets. Using massively parallel reporter assays in six cell types, we discovered 800 hCONDELs conferring significant differences in regulatory activity, half of which enhance rather than disrupt regulatory function. We highlight several hCONDELs with putative human-specific effects on brain development, including HDAC5, CPEB4, and PPP2CA. Reverting an hCONDEL to the ancestral sequence alters the expression of LOXL2 and developmental genes involved in myelination and synaptic function. Our data provide a rich resource to investigate the evolutionary mechanisms driving new traits in humans and other species.

Leveraging Base Pair Mammalian Constraint to Understand Genetic Variation and Human Disease

Although thousands of genomic regions have been associated with heritable human diseases, attempts to elucidate biological mechanisms are impeded by a general inability to discern which genomic positions are functionally important. Evolutionary constraint is a powerful predictor of function that is agnostic to cell type or disease mechanism. Here, single base phyloP scores from the whole genome alignment of 240 placental mammals identified 3.5% of the human genome as significantly constrained, and likely functional. We compared these scores to large-scale genome annotation, genome-wide association studies (GWAS), copy number variation, clinical genetics findings, and cancer data sets. Evolutionarily constrained positions are enriched for variants explaining common disease heritability (more than any other functional annotation). Our results improve variant annotation but also highlight that the regulatory landscape of the human genome still needs to be further explored and linked to disease.

A genomic and morphometric analysis of alpine bumblebees: Ongoing reductions in tongue length but no clear genetic component

Over the last six decades, populations of the bumblebees Bombus sylvicola and Bombus balteatus in Colorado have experienced decreases in tongue length, a trait important for plant-pollinator mutualisms. It has been hypothesized that this observation reflects selection resulting from shifts in floral composition under climate change. Here we used morphometrics and population genomics to determine whether morphological change is ongoing, investigate the genetic basis of morphological variation, and analyse population structure in these populations. We generated a genome assembly of B. balteatus. We then analysed whole-genome sequencing data and morphometric measurements of 580 samples of both species from seven high-altitude localities. Out of 281 samples originally identified as B. sylvicola, 67 formed a separate genetic cluster comprising a newly-discovered cryptic species ("incognitus"). However, an absence of genetic structure within species suggests that gene flow is common between mountains. We found a significant decrease in tongue length between bees collected between 2012-2014 and in 2017, indicating that morphological shifts are ongoing. We did not discover any genetic associations with tongue length, but a SNP related to production of a proteolytic digestive enzyme was implicated in body size variation. We identified evidence of covariance between kinship and both tongue length and body size, which is suggestive of a genetic component of these traits, although it is possible that shared environmental effects between colonies are responsible. Our results provide evidence for ongoing modification of a morphological trait important for pollination and indicate that this trait probably has a complex genetic and environmental basis.

Massive genome inversion drives coexistence of divergent morphs in common quails

The presence of population-specific phenotypes often reflects local adaptation or barriers to gene flow. The co-occurrence of phenotypic polymorphisms that are restricted within the range of a highly mobile species is more difficult to explain. An example of such polymorphisms is in the common quail Coturnix coturnix, a small migratory bird that moves widely during the breeding season in search of new mating opportunities, following ephemeral habitats,^1,2 and whose females may lay successive clutches at different locations while migrating.³ In spite of this vagility, previous studies reported a higher frequency of heavier males with darker throat coloration in the southwest of the distribution (I. Jiménez-Blasco et al., 2015, Int. Union Game Biol., conference). We used population genomics and cytogenetics to explore the basis of this polymorphism and discovered a large inversion in the genome of the common quail. This inversion extends 115 Mbp in length and encompasses more than 7,000 genes (about 12% of the genome), producing two very different forms. Birds with the inversion are larger, have darker throat coloration and rounder wings, are inferred to have poorer flight efficiency, and are geographically restricted despite the high mobility of the species. Stable isotope analyses confirmed that birds carrying the inversion have shorter migratory distances or do not migrate. However, we found no evidence of pre- or post-zygotic isolation, indicating the two forms commonly interbreed and that the polymorphism remains locally restricted because of the effect on behavior. This illustrates a genomic mechanism underlying maintenance of geographically structured polymorphisms despite interbreeding with a lineage with high mobility.

Genomic, Habitat, and Leaf Shape Analyses Reveal a Possible Cryptic Species and Vulnerability to Climate Change in a Threatened Daisy

Olearia pannosa is a plant species listed as vulnerable in Australia. Two subspecies are currently recognised (O. pannosa subsp. pannosa (silver daisy) and O. pannosa subsp. cardiophylla (velvet daisy)), which have overlapping ranges but distinct leaf shape. Remnant populations face threats from habitat fragmentation and climate change. We analysed range-wide genomic data and leaf shape variation to assess population diversity and divergence and to inform conservation management strategies. We detected three distinct genetic groupings and a likely cryptic species. Samples identified as O. pannosa subsp. cardiophylla from the Flinders Ranges in South Australia were genetically distinct from all other samples and likely form a separate, range-restricted species. Remaining samples formed two genetic clusters, which aligned with leaf shape differences but not fully with current subspecies classifications. Levels of genetic diversity and inbreeding differed between the three genetic groups, suggesting each requires a separate management strategy. Additionally, we tested for associations between genetic and environmental variation and carried out habitat suitability modelling for O. pannosa subsp. pannosa populations. We found mean annual maximum temperature explained a significant proportion of genomic variance. Habitat suitability modelling identified mean summer maximum temperature, precipitation seasonality and mean annual rainfall as constraints on the distribution of O. pannosa subsp. pannosa, highlighting increasing aridity as a threat for populations located near suitability thresholds. Our results suggest maximum temperature is an important agent of selection on O. pannosa subsp. pannosa and should be considered in conservation strategies. We recommend taxonomic revision of O. pannosa and provide conservation management recommendations.

Genetic Barriers to Historical Gene Flow between Cryptic Species of Alpine Bumblebees Revealed by Comparative Population Genomics

Evidence is accumulating that gene flow commonly occurs between recently diverged species, despite the existence of barriers to gene flow in their genomes. However, we still know little about what regions of the genome become barriers to gene flow and how such barriers form. Here, we compare genetic differentiation across the genomes of bumblebee species living in sympatry and allopatry to reveal the potential impact of gene flow during species divergence and uncover genetic barrier loci. We first compared the genomes of the alpine bumblebee Bombus sylvicola and a previously unidentified sister species living in sympatry in the Rocky Mountains, revealing prominent islands of elevated genetic divergence in the genome that colocalize with centromeres and regions of low recombination. This same pattern is observed between the genomes of another pair of closely related species living in allopatry (B. bifarius and B. vancouverensis). Strikingly however, the genomic islands exhibit significantly elevated absolute divergence (d_XY) in the sympatric, but not the allopatric, comparison indicating that they contain loci that have acted as barriers to historical gene flow in sympatry. Our results suggest that intrinsic barriers to gene flow between species may often accumulate in regions of low recombination and near centromeres through processes such as genetic hitchhiking, and that divergence in these regions is accentuated in the presence of gene flow.

Extreme Differences in Recombination Rate between the Genomes of a Solitary and a Social Bee

Social insect genomes exhibit the highest rates of crossing over observed in plants and animals. The evolutionary causes of these extreme rates are unknown. Insight can be gained by comparing recombination rate variation across the genomes of related social and solitary insects. Here, we compare the genomic recombination landscape of the highly social honey bee, Apis mellifera, with the solitary alfalfa leafcutter bee, Megachile rotundata, by analyzing patterns of linkage disequilibrium in population-scale genome sequencing data. We infer that average recombination rates are extremely elevated in A. mellifera compared with M. rotundata. However, our results indicate that similar factors control the distribution of crossovers in the genomes of both species. Recombination rate is significantly reduced in coding regions in both species, with genes inferred to be germline methylated having particularly low rates. Genes with worker-biased patterns of expression in A. mellifera and their orthologs in M. rotundata have higher than average recombination rates in both species, suggesting that selection for higher diversity in genes involved in worker caste functions in social taxa is not the explanation for these elevated rates. Furthermore, we find no evidence that recombination has modulated the efficacy of selection among genes during bee evolution, which does not support the hypothesis that high recombination rates facilitated positive selection for new functions in social insects. Our results indicate that the evolution of sociality in insects likely entailed selection on modifiers that increased recombination rates genome wide, but that the genomic recombination landscape is determined by the same factors.

Social Parasitism in the Honeybee (Apis mellifera) Is Not Controlled by a Single SNP

The Cape bee (Apis mellifera capensis) is a subspecies of the honeybee, in which workers commonly lay diploid unfertilized eggs via a process known as thelytoky. A recent study aimed to map the genetic basis of this trait in the progeny of a single capensis queen where workers laid either diploid (thelytokous) or haploid (arrhenotokous) eggs. A nonsynonymous single nucleotide polymorphism (SNP) in a gene of unknown function was reported to be strongly associated with thelytoky in this colony. Here, we analyze genome sequences from a global sample of A. mellifera and identify populations where the proposed thelytoky allele at this SNP is common but thelytoky is absent. We also analyze genome sequences of three capensis queens produced by thelytoky and find that, contrary to predictions, they do not carry the proposed thelytoky allele. The proposed SNP is therefore neither sufficient nor required to produce thelytoky in A. mellifera.

The genomic basis of adaptation to high-altitude habitats in the eastern honey bee (Apis cerana)

The eastern honey bee (Apis cerana) is of central importance for agriculture in Asia. It has adapted to a wide variety of environmental conditions across its native range in southern and eastern Asia, which includes high-altitude regions. eastern honey bees inhabiting mountains differ morphologically from neighbouring lowland populations and may also exhibit differences in physiology and behaviour. We compared the genomes of 60 eastern honey bees collected from high and low altitudes in Yunnan and Gansu provinces, China, to infer their evolutionary history and to identify candidate genes that may underlie adaptation to high altitude. Using a combination of F_ST -based statistics, long-range haplotype tests and population branch statistics, we identified several regions of the genome that appear to have been under positive selection. These candidate regions were strongly enriched for coding sequences and had high haplotype homozygosity and increased divergence specifically in highland bee populations, suggesting they have been subjected to recent selection in high-altitude habitats. Candidate loci in these genomic regions included genes related to reproduction and feeding behaviour in honey bees. Functional investigation of these candidate loci is necessary to fully understand the mechanisms of adaptation to high-altitude habitats in the eastern honey bee.

Functional acclimation across microgeographic scales in Dodonaea viscosa

Intraspecific plant functional trait variation provides mechanistic insight into persistence and can infer population adaptive capacity. However, most studies explore intraspecific trait variation in systems where geographic and environmental distances co-vary. Such a design reduces the certainty of trait-environment associations, and it is imperative for studies that make trait-environment associations be conducted in systems where environmental distance varies independently of geographic distance. Here we explored trait variation in such a system, and aimed to: (i) quantify trait variation of parent and offspring generations, and associate this variation to parental environments; (ii) determine the traits which best explain population differences; (iii) compare parent and offspring trait-trait relationships. We characterized 15 plant functional traits in eight populations of a shrub with a maximum separation ca. 100 km. Populations differed markedly in aridity and elevation, and environmental distance varied independently of geographic distance. We measured traits in parent populations collected in the field, as well as their offspring reared in greenhouse conditions. Parent traits regularly associated with their environment. These associations were largely lost in the offspring generation, indicating considerable phenotypic plasticity. An ordination of parent traits showed clear structure with strong influence of leaf area, specific leaf area, stomatal traits, isotope δ¹³C and δ¹⁵N ratios, and N_area, whereas the offspring ordination was less structured. Parent trait-trait correlations were in line with expectations from the leaf economic spectrum. We show considerable trait plasticity in the woody shrub over microgeographic scales (<100 km), indicating it has the adaptive potential within a generation to functionally acclimate to a range of abiotic conditions. Since our study shrub is commonly used for restoration in southern Australia and local populations do not show strong genetic differentiation in functional traits, the potential risks of transferring seed across the broad environmental conditions are not likely to be a significant issue.

Bioclimatic transect networks: Powerful observatories of ecological change

Transects that traverse substantial climate gradients are important tools for climate change research and allow questions on the extent to which phenotypic variation associates with climate, the link between climate and species distributions, and variation in sensitivity to climate change among biomes to be addressed. However, the potential limitations of individual transect studies have recently been highlighted. Here, we argue that replicating and networking transects, along with the introduction of experimental treatments, addresses these concerns. Transect networks provide cost‐effective and robust insights into ecological and evolutionary adaptation and improve forecasting of ecosystem change. We draw on the experience and research facilitated by the Australian Transect Network to demonstrate our case, with examples, to clarify how population‐ and community‐level studies can be integrated with observations from multiple transects, manipulative experiments, genomics, and ecological modeling to gain novel insights into how species and systems respond to climate change. This integration can provide a spatiotemporal understanding of past and future climate‐induced changes, which will inform effective management actions for promoting biodiversity resilience.

Targeted capture to assess neutral genomic variation in the narrow-leaf hopbush across a continental biodiversity refugium

The Adelaide geosyncline, a mountainous region in central southern Australia, is purported to be an important continental refugium for Mediterranean and semi-arid Australian biota, yet few population genetic studies have been conducted to test this theory. Here, we focus on a plant species distributed widely throughout the region, the narrow-leaf hopbush, Dodonaea viscosa ssp. angustissima, and examine its genetic diversity and population structure. We used a hybrid-capture target enrichment technique to selectively sequence over 700 genes from 89 individuals across 17 sampling locations. We compared 815 single nucleotide polymorphisms among individuals and populations to investigate population genetic structure. Three distinct genetic clusters were identified; a Flinders/Gammon ranges cluster, an Eastern cluster, and a Kangaroo Island cluster. Higher genetic diversity was identified in the Flinders/Gammon Ranges cluster, indicating that this area is likely to have acted as a refugium during past climate oscillations. We discuss these findings and consider the historical range dynamics of these populations. We also provide methodological considerations for population genomics studies that aim to use novel genomic approaches (such as target capture methods) on non-model systems. The application of our findings to restoration of this species across the region are also considered.

Transcriptome sequencing, annotation and polymorphism detection in the hop bush, Dodonaea viscosa

Background

The hop bush, Dodonaea viscosa, is a trans-oceanic species distributed oversix continents. It evolved in Australia where it is found over a wide range of habitat types and is an ecologically important species. Limited genomic resources are currently available for this species, thus our understanding of its evolutionary history and ecological adaptation is restricted. Here, we present a comprehensive transcriptome dataset for future genomic studies into this species.

Methods

We performed Illumina sequencing of cDNA prepared from leaf tissue collected from seven populations of D. viscosa ssp. angustissima and spatulata distributed along an environmental gradient in South Australia. Sequenced reads were assembled to provide a transcriptome resource. Contiguous sequences (contigs) were annotated using BLAST searches against the NCBI non-redundant database and gene ontology definitions were assigned. Single nucleotide polymorphisms were detected for the establishment of a genetic marker set. A comparison between the two subspecies was also carried out.

Results

Illumina sequencing returned 268,672,818 sequence reads, which were de novoassembled into 105,125 contigs. Contigs with significant BLAST alignments (E value < 1e^-5)numbered at 44,191, with 38,311 of these having their most significant hits to sequences from land plant species. Gene Ontology terms were assigned to 28,440 contigs and KEGG analysis identified 146 pathways that the gene products from 5,070 contigs are potentially involved in. The subspecies comparison identified 8,494 fixed SNP differences across 3,979 contiguous sequences, indicating a level of genetic differentiation between them. Across all samples, 248,235 SNPs were detected.

Conclusions

We have established a significant genomic data resource for D. viscosa,providing a comprehensive transcriptomic reference. Genetic differences among morphologically distinct subspecies were found. A wide range of putative gene regions were identified along with a large set of variable SNP markers, providing a basis for studies into the evolution and ecological adaptation of D. viscosa.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1987-1) contains supplementary material, which is available to authorized users.

Higher levels of multiple paternities increase seedling survival in the long-lived tree Eucalyptus gracilis

Studying associations between mating system parameters and fitness in natural populations of trees advances our understanding of how local environments affect seed quality, and thereby helps to predict when inbreeding or multiple paternities should impact on fitness. Indeed, for species that demonstrate inbreeding avoidance, multiple paternities (i.e. the number of male parents per half-sib family) should still vary and regulate fitness more than inbreeding--named here as the 'constrained inbreeding hypothesis'. We test this hypothesis in Eucalyptus gracilis, a predominantly insect-pollinated tree. Fifty-eight open-pollinated progeny arrays were collected from trees in three populations. Progeny were planted in a reciprocal transplant trial. Fitness was measured by family establishment rates. We genotyped all trees and their progeny at eight microsatellite loci. Planting site had a strong effect on fitness, but seed provenance and seed provenance × planting site did not. Populations had comparable mating system parameters and were generally outcrossed, experienced low biparental inbreeding and high levels of multiple paternity. As predicted, seed families that had more multiple paternities also had higher fitness, and no fitness-inbreeding correlations were detected. Demonstrating that fitness was most affected by multiple paternities rather than inbreeding, we provide evidence supporting the constrained inbreeding hypothesis; i.e. that multiple paternity may impact on fitness over and above that of inbreeding, particularly for preferentially outcrossing trees at life stages beyond seed development.