Summary from the authors: A linked-read approach to museomics: Higher quality de novo genome assemblies from degraded tissues

We aimed to sequence and compare all the DNA (eg., the genome) of a bunch of different deer mice (genus Peromyscus) species to understand how some deer mice survive in hot deserts with little to no water. A number of deer mice tissue samples were available through natural history museums, which house the raw materials for genetic and biodiversity investigations, but the samples had been collected many years earlier. Older samples produce lower quality DNA that has been broken into many pieces over time. Our normal sequencing procedure selectively removes small fragments of DNA, which would essentially throw away all the DNA we wanted to sequence for these older samples! To circumvent this, we were able to use a different DNA library preparation method called linked-read sequencing (LRS). LRS uses standard short-read sequencing technology, but adds additional information about the location of DNA fragments within the genome by bundling and barcoding DNA fragments that are located near each other prior to sequencing (eg., ‘links’ DNA fragments together in ‘genome-space’). We found that this method improves the overall quality and completeness of genome assemblies from historical tissue samples, in less time and with less effort than traditional shot-gun sequencing methods. This alternative method may be particularly valuable for building high-quality genome assemblies for extinct species for which there are no new samples being collected for or endangered species that are difficult to sample or collect. LRS adds to the suite of genomic methods that continue to unlock the secrets of natural history collections and enable fine-scale genetic measurement of change through time.

This summary was written by the study’s first author, Jocelyn Colella.

Read the full text here.

Video credit: Jocelyn Colella. Peromyscus in the field.

Full Text: Colella JP, Tigano A, MacManes MD. A linked-read approach to museomics: Higher quality de novo genome assemblies from degraded tissues. Mol Ecol Resour. 2020;20:871–881.

Summary from the authors: A metagenomic assessment of microbial eukaryotic diversity in the global ocean

Marine microbial eukaryotes are key components of planktonic ecosystems in all ocean biomes. They are, along with cyanobacteria, responsible for nearly half of the global primary production, and play important roles in food-web dynamics as grazers and parasites, carbon export to the deep ocean, and nutrient remineralization. Currently, one of the most common approaches to survey their diversity is sequencing marker genes amplified from genomic DNA extracted from microbial assemblages. However, this approach requires a PCR step, which is known to introduce biases in microbial diversity estimates. One alternative to overcome this issue involves exploiting the taxonomic information contained in metagenomes, which use massive shotgun sequencing of the same DNA extracts with the goal of assessing the putative functions of environmental microbes.

In this study we investigated the potential of metagenomics to provide taxonomic reports of marine microbial eukaryotes. The overall diversity reported by this approach was similar to that obtained by amplicon sequencing, although the latter performed poorly for some taxonomic groups. We then studied the diversity of picoeukaryotes and nanoeukaryotes using 91 metagenomes from surface down to bathypelagic layers in different oceans, unveiling a clear separation of taxonomic groups between size fractions and depth layers.

Overall, this study shows metagenomics as an excellent resource for taxonomic exploration of marine microbial eukaryotes.

Summary of the relevance of main eukaryotic taxonomic groups within two size fractions of marine plankton (picoeukaryotes [0.2-3 µm] and nanoeukaryotes [3-20µm]) and in two different layers of the global ocean (photic [0-200 m] and aphotic [200-4000m]) as seen by metagenomics. The median of the relative abundance was calculated for each taxonomic group with samples from the 4 categories (pico-photic, pico-aphotic, nano-photic, nano-aphotic) and dots represent these median values transformed to a 0-100 scale. Dots are then colored based on the category where the taxonomic group is most relevant.

This summary was written by the study’s first author, Aleix Obiol.

Full article:
Obiol, A., Giner, C. R., Sánchez, P., Duarte, C. M., Acinas, S. G., & Massana, R. (2020). A metagenomic assessment of microbial eukaryotic diversity in the global ocean. Molecular Ecology Resources. https://doi.org/10.1111/1755-0998.13147

Summary from the authors: Latent Dirichlet Allocation reveals spatial and taxonomic structure in a DNA-based census of soil biodiversity from a tropical forest

Biodiversity inventories can now be built by collecting and sequencing DNA from the environment, which is not only easier, faster and cheaper than direct observation, but also much more comprehensive and systematic. This gives in particular unprecedented access to little-known microbial diversity. Tapping these data to answer community ecology questions, however, can prove a daunting task, as classical statistical approaches often fall short of the size and complexity of molecular datasets. To uncover the spatial structure of soil biodiversity over 12 ha of primary tropical forest in French Guiana, we borrowed a probabilistic model from text analysis. After demonstrating the performance of the method on simulated data, we used it to capture the co-occurrence and covariance patterns of more than 25,000 taxa of bacteria, protists, fungi and metazoans across 1,131 soil samples, collected every 10 m – a dataset that led to a previous publication in Mol. Ecol. (Zinger et al. 2019). We find that, even though the forest plot is at first sight rather uniform, bacteria, protists and fungi are all clearly structured into three assemblages matching the environmental heterogeneity of the plot, whereas metazoans are unstructured at that scale. We then work though the practical problems ecologists may encounter using this approach, such as whether to use presence-absence or read-count data, how to choose the number of assemblages and how to assess the robustness of the results. Finally, we discuss the potential use of related methods in community ecology and biogeography, and argue that probabilistic models are a way forward for analyzing the ever-expanding amount of data generated by the field.

Left: Primary tropical forest understory on the plot where the data were collected, Nouragues ecological research station, French Guiana. Right: Spatial distribution of assemblages of co-occurring soil taxa (OTUs), obtained by Latent Dirichlet Allocation from OTU presence-absence information only, over a 12-ha plot of plateau forest sampled every 10 m (top); and two main axes of environmental variation over the same forest plot, derived from Airborne Laser Scanning (bottom). Bacteria, protists and fungi exhibit a spatial pattern matching the environmental heterogeneity of the plot: the blue, green and red assemblages match the terra firme, hydromorphic and exposed rock parts of the plot, respectively. In contrast, metazoans such as annelids can be shown to be spatially unstructured at that scale. Sampled locations are indicated by dark dots, and values have been interpolated between samples using kriging.

Full article:
Sommeria-Klein G, Zinger L, Coissac E, et al. (2020). Latent Dirichlet Allocation reveals spatial and taxonomic structure in a DNA-based census of soil biodiversity from a tropical forest. Molecular Ecology Resour. 20:371–386. https://doi.org/10.1111/1755-0998.13109

References:
Zinger, L., Taberlet, P., et al. (2019). Body size determines soil community assembly
in a tropical forest. Molecular Ecology, 28(3), 528–543. https://doi.org/10.1111/mec.14919

Summary from the authors: inbreeding and management in captive populations

Pacific salmon hatcheries aim to supplement declining wild populations and support commercial and recreational fisheries. However, there are also risks associated with hatcheries because the captive and wild environments are inherently different. It is important to understand these risks in order to maximize the success of hatcheries. Inbreeding, which occurs when related individuals interbreed, is one risk that may inadvertently be higher in hatcheries due to space limitations and other factors. Inbred fish may have reduced fitness and survival compared to non-inbred fish. We quantified inbreeding and its effect on key fitness traits across four generations in two hatchery populations of adult Chinook salmon that were derived from the same source. We utilized recent advancements in DNA sequencing technology, which provide much more precise estimates of inbreeding and its potential effects on fitness. Our results indicate that inbreeding may not be severe in salmon hatcheries, even small ones, provided that appropriate management practices are followed. However, we documented an influence of inbreeding on the phenology of adult spawners, which could have biological implications for individual fitness and population productivity. Our findings provide a better understanding of changes that may occur in hatchery salmon and will further inform research on “best” hatchery practices to minimize potential risks. 

Article: Waters CD, Hard JJ, Fast DE, Knudsen CM, Bosch WJ, Naish KA. 2020. Genomic and phenotypic effects of inbreeding across two different hatchery management regimes in Chinook salmon. Molecular Ecology https://doi.org/10.1111/mec.15356.

Summary from the authors: Individualized mating system estimation using genomic data

Mimulus guttatus

Hermaphroditic species of plants and animal can produce a mixture of outcrossed and self-fertilized offspring. Estimating the relative frequency of these two outcomes, i.e. the outcrossing rate, has been a major focus in the evolutionary study of reproductive strategies. Outcrossing rate is also a key parameter for plant breeding and for conservation efforts. This paper generalizes a Bayesian method to estimate outcrossing rate (BORICE) using genomic data. Application of the program to an experimental study of Mimulus guttatus illustrates estimation (10% of progeny were selfed), and also how inference of mating system parameters can set up “downstream” evolutionary studies. In the Mimulus study, these downstream analyses included pollination biology (the genetic composition of pollen changed over the season) and local adaptation (inversion polymorphisms exhibit unique patterns of micro spatial structure within the population).

-Professor John K Kelly, University of Kansas

Full article: Colicchio, J., Monnahan, P. J., Wessinger, C. A., Brown, K., Kern, J. R., & Kelly, J. K. (2020). Individualized mating system estimation using genomic data. Molecular ecology resources. https://doi.org/10.1111/1755-0998.13094

Summary from the authors: telomere length predicts remaining lifespan

Close-up of an adult common tern with its prey. Photo credit: Andrea Parisi

Telomeres are DNA structures located at the end of chromosomes. They protect the chromosome, but shorten at each cell division. When telomeres get too short, the normal functioning of cells can be impaired. An individual’s telomere length may therefore predict its future lifespan, and understanding individual telomere dynamics could help to understand ageing in general.

Telomere shortening can be accelerated due to stress, thereby acting as a biomarker of an individual’s health status. However, some studies suggest that individual differences in telomere length are already determined at birth, and largely consistent over life.

We investigated individual telomere dynamics in a long-lived seabird, the common tern. The telomere lengths of 387 individuals, aged from 2 to 24 years, were repeatedly sampled across 10 years. We found that an individual’s telomeres shortened as they got older. Telomere shortening was also slightly increased if individuals had produced more chicks in the previous year. However, the correlation between repeated measures of an individual’s telomere length was very high, even with 6 years between measures. Nevertheless, an individual’s telomere length positively predicted its remaining lifespan, leaving the question of whether lifespan is already partly determined at the start of life.

Full article: Bichet C, Bouwhuis S, Bauch C, Verhulst S, Becker PH, Vedder O. 2019. Telomere length is repeatable, shortens with age and reproductive success, and predicts remaining lifespan in a long-lived seabird. Molecular ecology. https://doi.org/10.1111/mec.15331

Summary from the authors: genetic architecture of sexual dimorphism in an interspecific cross

The evolution of differences among females and males or sexual dimorphism (SD) is very common in animals but rare in plants. These differences emerge because there is a conflict of interests between sexes to maximize their reproductive success. Thus,  moving genes of reproductive traits to low recombining regions such as the sex chromosomes might be one way to solve this conflict at the genomic level. Closely related species with young sex chromosomes, which differ in the degree of SD, are ideal systems to explore the underlining genetic architecture of SD. We have crossed a female from Silene latifolia with marked SD with a male from S. dioica with less SD. We performed a QTL analysis of reproductive and vegetative traits in the F2 hybrids to find out if sexually dimorphic traits are located on the sex chromosomes, and how they contribute to species differences. Our results support that evolutionary young sex chromosomes are important for the expression of both SD and species differences. Moreover, transgressive segregation (traits with extreme values) and a reversal of SD in the F2s indicated that SD is constrained within the species but not in the recombinant hybrids. Sexual selection can, thus, contribute to speciation.

Full article: Baena-Díaz F, Zemp N, Widmar A. 2019. Insights into the genetic architecture of sexual dimorphism from an interspecific cross between two diverging Silene (Caryophyllaceae) species. Molecular ecology. https://doi.org/10.1111/mec.15271

Summary from the authors: City life alters the gut microbiome and stable isotope profiling of the eastern water dragon (Intellagama lesueurii)

Whilst urbanisation poses a major threat to many species, there is growing evidence to suggest that some species, labelled ‘urban adapters’, are thriving within the urban landscape. Urban landscapes differ drastically from native habitats, where urban adapters are often exposed to a more diverse range of novel food items compared to their rural counterparts, which frequently includes human subsidised resources. Diet is one of the most important factors influencing the gut microbiome, an extremely influential symbiotic community that plays a critical role in many processes affecting host health and fitness, including metabolism, nutrition, immunology and development. Here, using populations of the eastern water dragon in Queensland, Australia, we explore the link between urbanisation, diet and gut microbial changes. We show that city dragons exhibit a more diverse gut microbiome than their rural counterparts, and display microbial signatures of a diet that is richer in plant-material and higher in fat. Elevated levels of the Nitrogen-15 isotope in the blood of city dragons also suggests their diet may be richer in protein. These results highlight that urbanisation can have pronounced effects on the gut microbial communities of wild animals, but we do not yet know the possible repercussions of these microbial changes.

Full article: Littleford‐Colquhoun BL, Weyrich LS, Kent N, Frere CH. City life alters the gut microbiome and stable isotope profiling of the eastern water dragon (Intellagama lesueurii). Mol Ecol. 2019;28:4592–4607. https://doi.org/10.1111/mec.15240

Summary from the authors: Host plant associations and geography interact to shape diversification in a specialist insect herbivore

The sexual generation female of the cynipid gall wasp Belonocnema treatae, which forms locally adapted populations on three closely related live oaks across the southeastern United States. Photo credit: Jena Johnson.

In this study, we wanted to know how geography and ecology predicted population genetic structure among 58 populations of the gall wasp Belonocnema treatae, which exhibits regional specialization on three host plant species across the U.S. Gulf Coast. We combined range-wide sampling with a genotype-by-sequencing approach for 40,699 SNPs across 1,217 individuals. Disentangling the processes underlying geographic and environmental patterns of biodiversity is challenging, as such patterns emerge from eco‐evolutionary processes confounded by spatial autocorrelation among sample units. We evaluated this question using a hierarchical Bayesian model (ENTROPY) to assign individuals to genetic clusters and estimate admixture proportions. Using distance-based Moran’s eigenvector mapping, we generated regression variables that represent varying degrees of spatial autocorrelation in genetic variation among sample sites. These spatial variables, along with host association, were incorporated in distance-based redundancy analysis (dbRDA) to partition the relative contributions of host plant and spatial autocorrelation. This novel approach of combining ENTROPY results with dbRDA to analyze SNP data unveiled a complex mosaic of diversification within and among insect populations forming discrete host associated lineages coupled with geographic variation. This demonstrates that geography and ecology play significant roles in explaining patterns of genomic variation in B. treatae – an emerging model of ecological speciation.

Full article: Driscoe AL, Nice CC, Busbee RW,Hood GR, Egan SP, Ott JR. Host plant associations and geography interact to shape diversification in a specialist insect herbivore. Mol Ecol. 2019;28:4197–4211. https://doi.org/10.1111/mec.15220

Summary from the authors: Polygenic selection drives the evolution of convergent transcriptomic landscapes across continents within a Nearctic sister-species complex

Lake Whitefish (Coregonus clupeaformis)
Adult size benthic (normal) and limnetic (dwarf) sympatric whitefish
Sampled from Cliff Lake, Maine, USA.
Illustration by Paul Vecsei.

The Lake Whitefish (Coregonus clupeaformis) is found in Nearctic post-glacial freshwater lakes, and diverged about 500,000 years ago from its sister species, the European Whitefish (Coregonus lavaretus) found throughout Northern Europe, European Alpine lakes and Russia. These two lineages underwent an adaptive radiation following the last glaciation, resulting in the sympatric occurrence of limnetic and benthic species-pairs. Decades of research has aimed to decipher the process of adaptive divergence and ecological speciation in the Whitefish species complex. Here, we compared independent diverging species-pairs from the two continents to elucidate the genomic and transcriptomic bases associated with the benthic-limnetic diversification. We used a statistical framework to detect polygenic targets of selection associated with phenotypic diversification. We identified a subset of genes that showed convergent patterns of differential expression between limnetic and benthic species across both continents. Those adaptive divergent genes retained a higher degree of shared polymorphism among species-pairs, most likely due to balancing selection, and this genetic variation was associated with changes in levels of gene expression between species. As such, our results indicate that standing genetic variation underlying phenotypes involved in the ecological speciation of the whitefish species-pairs has been partly maintained in parallel across both continents for at least half a million years.

-Clément Rougeux – Postdoctoral Researcher, University of Calgary

Full article: Rougeux C, Gagnaire P‐A, Praebel K, Seehausen O, Bernatchez L. Polygenic selection drives the evolution of convergent transcriptomic landscapes across continents within a Nearctic sister species complex. Mol Ecol. 2019;28:4388–4403. https://doi.org/10.1111/mec.15226