Divergence thresholds and divergent biodiversity estimates: can metabarcoding reliably describe zooplankton communities?

DNA metabarcoding is a promising method for describing communities and estimating biodiversity. This approach uses high‐throughput sequencing of targeted markers to identify species in a complex sample. By convention, sequences are clustered at a predefined sequence divergence threshold (often 3%) into operational taxonomic units (OTUs) that serve as a proxy for species. However, variable levels of interspecific marker variation across taxonomic groups make clustering sequences from a phylogenetically diverse dataset into OTUs at a uniform threshold problematic. In this study, we use mock zooplankton communities to evaluate the accuracy of species richness estimates when following conventional protocols to cluster hypervariable sequences of the V4 region of the small subunit ribosomal RNA gene (18S) into OTUs. By including individually tagged single specimens and “populations” of various species in our communities, we examine the impact of intra‐ and interspecific diversity on OTU clustering. Communities consisting of single individuals per species generated a correspondence of 59–84% between OTU number and species richness at a 3% divergence threshold. However, when multiple individuals per species were included, the correspondence between OTU number and species richness dropped to 31–63%. Our results suggest that intraspecific variation in this marker can often exceed 3%, such that a single species does not always correspond to one OTU. We advocate the need to apply group‐specific divergence thresholds when analyzing complex and taxonomically diverse communities, but also encourage the development of additional filtering steps that allow identification of artifactual rRNA gene sequences or pseudogenes that may generate spurious OTUs.     

In situ hatching of invertebrate diapausing eggs from ships' ballast sediment

Ships that enter the Great Lakes laden with cargo carry only residual ballast water and sediment in ballast tanks. These ships are designated ‘no ballast on board’ (NOBOB) and constitute > 90% of inbound traffic. We conducted in situ experiments using emergence traps to assess the viability and the introduction potential of invertebrate diapausing stages present in ships’ ballast sediment. All trials commenced while vessels operated on the lower lakes (Erie, Ontario) and were completed 6–11 days later at ports on the upper lakes (Michigan, Lake Superior). Eight trials were conducted on four ships using five different ballast sediments. Hatching was observed on every ship, although not from all sediments on all ships. Overall hatch rates were very low (0.5 individuals per 500 g sediment), typically involving activation of < 0.05% of total eggs present. Five species of rotifers and copepod nauplii were hatched from ballast sediments, although only one or two species typically hatched from any one sediment. Results of this study indicate that hatching of diapausing eggs contained in ballast sediment of NOBOB ships poses a relatively low risk of invasion to the Great Lakes. However, as reproduction may occur in tanks, and non‐indigenous species may be involved in numerous introduction events, the risk posed by this vector is small but potentially important. While dormancy is a characteristic enabling enhanced survival during transportation in ballast tanks, it becomes a hindrance for introduction.     

Drosophilid Assemblages at Different Urbanization Levels in the City of Porto Alegre, State of Rio Grande do Sul, Southern Brazil

The present study analyzed the drosophilid assemblages in different levels of urbanization in the city of Porto Alegre, Rio Grande do Sul, Brazil. Collections were carried out in 2008 in three different environments: a highly urbanized area????Jardim Botanico,?? a forested area with intermediary urbanization????Parque Gabriel Knijnik,?? and in a relatively well-preserved forested area, although threatened by the urban growth????Morro Santana.?? In Jardim Botanico, 36 species belonging to four genera were found, with high abundance of exotic species as Drosophila simulans Sturtevant and Zaprionus indianus (Gupta). In Parque Gabriel Knijnik, 33 species that belonged to four genera were found, with higher abundances of native species belonging to the Drosophila tripunctata species group and Drosophila willistoni species subgroup, and lower abundance of exotic species. As for Morro Santana, 32 species and three genera were found, with higher abundances of native groups, low representativeness of exotic species, and absence of Zaprionus indianus. The analysis of the Jaccard index showed higher similarity in the species composition between samples collected in summer and autumn, and between samples collected in winter and spring. On the other hand, the Morisita index differentiated Jardim Botanico from the other two studied sites. Our results show that Morro Santana is an important area of native biodiversity, reinforcing, therefore, the inclusion of this area in the project for the creation of an ecological corridor as proposed by the Ministry of the Environment of Brazil.     

Using the realized relationship matrix to disentangle confounding factors for the estimation of genetic variance components of complex traits

Background

In the analysis of complex traits, genetic effects can be confounded with non-genetic effects, especially when using full-sib families. Dominance and epistatic effects are typically confounded with additive genetic and non-genetic effects. This confounding may cause the estimated genetic variance components to be inaccurate and biased.

Methods

In this study, we constructed genetic covariance structures from whole-genome marker data, and thus used realized relationship matrices to estimate variance components in a heterogenous population of ~ 2200 mice for which four complex traits were investigated. These mice were genotyped for more than 10,000 single nucleotide polymorphisms (SNP) and the variances due to family, cage and genetic effects were estimated by models based on pedigree information only, aggregate SNP information, and model selection for specific SNP effects.

Results and conclusions

We show that the use of genome-wide SNP information can disentangle confounding factors to estimate genetic variances by separating genetic and non-genetic effects. The estimated variance components using realized relationship were more accurate and less biased, compared to those based on pedigree information only. Models that allow the selection of individual SNP in addition to fitting a relationship matrix are more efficient for traits with a significant dominance variance.     

The effect of genomic information on optimal contribution selection in livestock breeding programs

Background

Long-term benefits in animal breeding programs require that increases in genetic merit be balanced with the need to maintain diversity (lost due to inbreeding). This can be achieved by using optimal contribution selection. The availability of high-density DNA marker information enables the incorporation of genomic data into optimal contribution selection but this raises the question about how this information affects the balance between genetic merit and diversity.

Methods

The effect of using genomic information in optimal contribution selection was examined based on simulated and real data on dairy bulls. We compared the genetic merit of selected animals at various levels of co-ancestry restrictions when using estimated breeding values based on parent average, genomic or progeny test information. Furthermore, we estimated the proportion of variation in estimated breeding values that is due to within-family differences.

Results

Optimal selection on genomic estimated breeding values increased genetic gain. Genetic merit was further increased using genomic rather than pedigree-based measures of co-ancestry under an inbreeding restriction policy. Using genomic instead of pedigree relationships to restrict inbreeding had a significant effect only when the population consisted of many large full-sib families; with a half-sib family structure, no difference was observed. In real data from dairy bulls, optimal contribution selection based on genomic estimated breeding values allowed for additional improvements in genetic merit at low to moderate inbreeding levels. Genomic estimated breeding values were more accurate and showed more within-family variation than parent average breeding values; for genomic estimated breeding values, 30 to 40% of the variation was due to within-family differences. Finally, there was no difference between constraining inbreeding via pedigree or genomic relationships in the real data.

Conclusions

The use of genomic estimated breeding values increased genetic gain in optimal contribution selection. Genomic estimated breeding values were more accurate and showed more within-family variation, which led to higher genetic gains for the same restriction on inbreeding. Using genomic relationships to restrict inbreeding provided no additional gain, except in the case of very large full-sib families.     

Viability of invertebrate diapausing eggs exposed to saltwater: implications for Great Lakes’ ship ballast management

International shipping has been the dominant vector of nonindigenous species introductions to the Laurentian Great Lakes over the past century. Apparent ballast-mediated invasions have been recorded in recent years, despite the implementation of voluntary ballast water exchange regulations in 1989. Since unregulated no-ballast-on-board vessels currently dominate inbound traffic to the Great Lakes, it has been proposed that live or dormant organisms contained in residual ballast of these vessels may be partially responsible for recent invasions. Alternatively, euryhaline species may pose a significant invasion threat because they can potentially survive ballast exchange. In this study, we explored whether exposure to open-ocean water (32) reduced the viability of invertebrate diapausing eggs in ballast sediments. Sediments collected from three transoceanic ships and from three freshwater habitats were exposed to open-ocean seawater. Egg viability, assessed as the abundance of taxa hatched between exposed and unexposed sediments, was not affected by saltwater exposure in any experiment. Species richness of hatched diapausing eggs was reduced by saltwater exposure in only one of seven trials. Our results indicate that oligostenohaline zooplankton may pose an invasion risk because their diapausing eggs are largely resistant to exposure to open-ocean saltwater.     

Evaluation of stochastic gravity model selection for use in estimating non-indigenous species dispersal and establishment

Predicting dispersal of nonindigenous species (NIS) is an essential component of risk analysis and management, as preventative measures are most readily applied at this stage of the invasion sequence. Gravity models provide one of the most useful techniques available to model dispersal of nonindigenous invasive species (NIS) across heterogeneous landscapes, as these models are able to capture transport patterns of recreational boaters who are dominant vectors responsible for aquatic NIS dispersal. Despite the widespread use of gravity models in forecasting biological invasions, different classes of gravity models have not been evaluated regarding their comparative ability to capture recreational transport patterns and subsequent use in predicting NIS establishment. Here we evaluate model selection between unconstrained, total-flow-constrained, production-constrained and doubly-constrained stochastic gravity models to assess dispersal of the spiny waterflea Bythotrephes between Ontario lakes. Differences between the models relate to the amount of data required and constraints under which calculations of source/destination interactions are made. For each class of gravity model, we then estimated the probability of a lake having established Bythotrephes populations by modeling the relationship between empirical presence/absence data and inbound recreational traffic (i.e. propagule pressure) via boosted regression. The unconstrained gravity model provided the best fit to observed traffic patterns of recreational boaters. However, when output from the gravity models was used to predict Bythotrephes establishment, the doubly-constrained gravity model provided the strongest relationship between inbound recreational traffic and observed Bythotrephes presence/absence, followed by the production-constrained model. Our results indicate production-constrained gravity models offer an acceptable balance between modeling recreational boater traffic and their utility to estimate establishment probabilities.     

YIL113w encodes a functional dual-specificity protein phosphatase which specifically interacts with and inactivates the Slt2/Mpk1p MAP kinase in S. cerevisiae

We show here that the YIL113w gene of Saccharomyces cerevisiae encodes a functional protein phosphatase. Yil113p shows no activity in vitro towards either phosphorylated casein or myelin basic protein. However, Yil113p dephosphorylates activated extracellular signal-regulated kinase 2 MAP kinase indicating that it is a dual-specificity MAP kinase phosphatase. In support of this we find that Yil113p specifically interacts with the stress-activated Slt2/Mpk1p MAP kinase of S. cerevisiae. Furthermore, expression of Yil113p causes the dephosphorylation of Slt2/Mpk1p in vivo, while expression of an inactive mutant of Yil113p causes the accumulation of phosphorylated Slt2/Mpk1p. We conclude that the physiological target of YIL113p is Slt2/Mpk1p.     

Low-multiplicity infection of Moloney murine leukemia virus in mouse cells: effect on number of viral DNA copies and virus production in producer cells.

Mouse cells infected with Moloney murine leukemia virus (M-MuLV) were prepared by two methods, and the number of M-MuLV-specific DNA copies in the infected cells was measured. The number of M-MuLV-specific DNA copies detected varied from one to eight per infected cell in different cell lines. Cells in which multiple rounds of viral infection occurred during establishment had on the average more viral DNA copies than cells in which infection at low multiplicity was performed, followed by cloning of the cells. However, even in cells derived by the low multiplicity of infection method, most cell lines carried more than one copy of M-MuLV-specific DNA. Virus production per cell was also measured, and no strict correlation was observed between the number of M-MuLV DNA copies present and the amount of virus produced.