Assessing microplastic exposure of large marine filter-feeders.

Large filter-feeding animals are potential sentinels for understanding the extent of microplastic pollution, as their mode of foraging and prey mean they are continuously sampling the environment. However, there is considerable uncertainty about the total and mode of exposure (environmental vs trophic). Here, we explore microplastic exposure and ingestion by baleen whales feeding year-round in coastal Auckland waters, New Zealand. Plastic and DNA were extracted concurrently from whale scat, with 32 ± 24 (mean ± SD, n = 21) microplastics per 6 g scat sample detected. Using a novel stochastic simulation modeling incorporating new and previously published DNA diet information, we extrapolate this to total microplastic exposure levels of 24,028 (95% CI: 2119, 69,270) microplastics per mouthful of prey, or 3,408,002 microplastics (95% CI: 295,810, 10,031,370) per day, substantially higher than previous estimates for large filter-feeding animals. Critically, we find that the total exposure is four orders of magnitude more than expected from microplastic measurements of local coastal surface waters. This suggests that trophic transfer, rather than environmental exposure, is the predominant mode of exposure of large filter feeders for microplastic pollution. Measuring plastic concentration from the environment alone significantly underestimates exposure levels, an important consideration for future risk assessment studies.

Multi-locus DNA metabarcoding of zooplankton communities and scat reveal trophic interactions of a generalist predator.

To understand the ecosystem dynamics that underpin the year-round presence of a large generalist consumer, the Bryde's whale (Balaenoptera edeni brydei), we use a DNA metabarcoding approach and systematic zooplankton surveys to investigate seasonal and regional changes in zooplankton communities and if whale diet reflects such changes. Twenty-four zooplankton community samples were collected from three regions throughout the Hauraki Gulf, New Zealand, over two temperature regimes (warm and cool seasons), as well as 20 samples of opportunistically collected Bryde's whale scat. Multi-locus DNA barcode libraries were constructed from 18S and COI gene fragments, representing a trade-off between identification and resolution of metazoan taxa. Zooplankton community OTU occurrence and relative read abundance showed regional and seasonal differences based on permutational analyses of variance in both DNA barcodes, with significant changes in biodiversity indices linked to season in COI only. In contrast, we did not find evidence that Bryde's whale diet shows seasonal or regional trends, but instead indicated clear prey preferences for krill-like crustaceans, copepods, salps and ray-finned fishes independent of prey availability. The year-round presence of Bryde's whales in the Hauraki Gulf is likely associated with the patterns of distribution and abundance of these key prey items.

Incorporating non-equilibrium dynamics into demographic history inferences of a migratory marine species.

Understanding how dispersal and gene flow link geographically separated the populations over evolutionary history is challenging, particularly in migratory marine species. In southern right whales (SRWs, Eubalaena australis), patterns of genetic diversity are likely influenced by the glacial climate cycle and recent history of whaling. Here we use a dataset of mitochondrial DNA (mtDNA) sequences (n = 1327) and nuclear markers (17 microsatellite loci, n = 222) from major wintering grounds to investigate circumpolar population structure, historical demography and effective population size. Analyses of nuclear genetic variation identify two population clusters that correspond to the South Atlantic and Indo-Pacific ocean basins that have similar effective breeder estimates. In contrast, all wintering grounds show significant differentiation for mtDNA, but no sex-biased dispersal was detected using the microsatellite genotypes. An approximate Bayesian computation (ABC) approach with microsatellite markers compared the scenarios with gene flow through time, or isolation and secondary contact between ocean basins, while modelling declines in abundance linked to whaling. Secondary-contact scenarios yield the highest posterior probabilities, implying that populations in different ocean basins were largely isolated and came into secondary contact within the last 25,000 years, but the role of whaling in changes in genetic diversity and gene flow over recent generations could not be resolved. We hypothesise that these findings are driven by factors that promote isolation, such as female philopatry, and factors that could promote dispersal, such as oceanographic changes. These findings highlight the application of ABC approaches to infer the connectivity in mobile species with complex population histories and, currently, low levels of differentiation.

First Direct Evidence for Natal Wintering Ground Fidelity and Estimate of Juvenile Survival in the New Zealand Southern Right Whale Eubalaena australis.

Juvenile survival and recruitment can be more sensitive to environmental, ecological and anthropogenic factors than adult survival, influencing population-level processes like recruitment and growth rate in long-lived, iteroparous species such as southern right whales. Conventionally, Southern right whales are individually identified using callosity patterns, which do not stabilise until 6-12 months, by which time the whale has left its natal wintering grounds. Here we use DNA profiling of skin biopsy samples to identify individual Southern right whales from year of birth and document their return to the species' primary wintering ground in New Zealand waters, the Subantarctic Auckland Islands. We find evidence of natal fidelity to the New Zealand wintering ground by the recapture of 15 of 57 whales, first sampled in year of birth and available for subsequent recapture, during winter surveys to the Auckland Islands in 1995-1998 and 2006-2009. Four individuals were recaptured at the ages of 9 to 11, including two females first sampled as calves in 1998 and subsequently resampled as cows with calves in 2007. Using these capture-recapture records of known-age individuals, we estimate changes in survival with age using Cormack-Jolly-Seber models. Survival is modelled using discrete age classes and as a continuous function of age. Using a bootstrap method to account for uncertainty in model selection and fitting, we provide the first direct estimate of juvenile survival for this population. Our analyses indicate a high annual apparent survival for juveniles at between 0.87 (standard error (SE) 0.17, to age 1) and 0.95 (SE 0.05: ages 2-8). Individual identification by DNA profiling is an effective method for long-term demographic and genetic monitoring, particularly in animals that change identifiable features as they develop or experience tag loss over time.

Cultural traditions across a migratory network shape the genetic structure of southern right whales around Australia and New Zealand.

Fidelity to migratory destinations is an important driver of connectivity in marine and avian species. Here we assess the role of maternally directed learning of migratory habitats, or migratory culture, on the population structure of the endangered Australian and New Zealand southern right whale. Using DNA profiles, comprising mitochondrial DNA (mtDNA) haplotypes (500 bp), microsatellite genotypes (17 loci) and sex from 128 individually-identified whales, we find significant differentiation among winter calving grounds based on both mtDNA haplotype (FST = 0.048, ΦST = 0.109, p < 0.01) and microsatellite allele frequencies (FST = 0.008, p < 0.01), consistent with long-term fidelity to calving areas. However, most genetic comparisons of calving grounds and migratory corridors were not significant, supporting the idea that whales from different calving grounds mix in migratory corridors. Furthermore, we find a significant relationship between δ(13)C stable isotope profiles of 66 Australian southern right whales, a proxy for feeding ground location, and both mtDNA haplotypes and kinship inferred from microsatellite-based estimators of relatedness. This indicates migratory culture may influence genetic structure on feeding grounds. This fidelity to migratory destinations is likely to influence population recovery, as long-term estimates of historical abundance derived from estimates of genetic diversity indicate the South Pacific calving grounds remain at <10% of pre-whaling abundance.

Maximum likelihood estimation for model Mt,α for capture-recapture data with misidentification.

We investigate model Mt,α  for abundance estimation in closed-population capture-recapture studies, where animals are identified from natural marks such as DNA profiles or photographs of distinctive individual features. Model Mt,α  extends the classical model Mt  to accommodate errors in identification, by specifying that each sample identification is correct with probability α and false with probability 1-α. Information about misidentification is gained from a surplus of capture histories with only one entry, which arise from false identifications. We derive an exact closed-form expression for the likelihood for model Mt,α  and show that it can be computed efficiently, in contrast to previous studies which have held the likelihood to be computationally intractable. Our fast computation enables us to conduct a thorough investigation of the statistical properties of the maximum likelihood estimates. We find that the indirect approach to error estimation places high demands on data richness, and good statistical properties in terms of precision and bias require high capture probabilities or many capture occasions. When these requirements are not met, abundance is estimated with very low precision and negative bias, and at the extreme better properties can be obtained by the naive approach of ignoring misidentification error. We recommend that model Mt,α  be used with caution and other strategies for handling misidentification error be considered. We illustrate our study with genetic and photographic surveys of the New Zealand population of southern right whale (Eubalaena australis).

Accounting for female reproductive cycles in a superpopulation capture-recapture framework.

Superpopulation capture-recapture models are useful for estimating the abundance of long-lived, migratory species because they are able to account for the fluid nature of annual residency at migratory destinations. Here we extend the superpopulation POPAN model to explicitly account for heterogeneity in capture probability linked to reproductive cycles (POPAN-tau). This extension has potential application to a range of species that have temporally variable life stages (e.g., non-annual breeders such as albatrosses and baleen whales) and results in a significant reduction in bias over the standard POPAN model. We demonstrate the utility of this model in simultaneously estimating abundance and annual population growth rate (lamda) in the New Zealand (NZ) southern right whale (Eubalaena australis) from 1995 to 2009. DNA profiles were constructed for the individual identification of more than 700 whales, sampled during two sets of winter expeditions in 1995-1998 and 2006-2009. Due to differences in recapture rates between sexes, only sex-specific models were considered. The POPAN-tau models, which explicitly account for a decrease in capture probability in non-calving years, fit the female data set significantly better than do standard superpopulation models (deltaAIC > 25). The best POPAN-tau model (AIC) gave a super-population estimate of 1162 females for 1995-2009 (95% CL 921, 1467) and an estimated annual increase of 5% (95% CL--2%, 13%). The best model (AIC) gave a superpopulation estimate of 1007 males (95% CL 794, 1276) and an estimated annual increase of 7% (95% CL 5%, 9%) for 1995-2009. Combined, the total superpopulation estimate for 1995-2009 was 2169 whales (95% CL 1836, 2563). Simulations suggest that failure to account for the effect of reproductive status on the capture probability would result in a substantial positive bias (+19%) in female abundance estimates.

Peptides specific for Mycobacterium avium subspecies paratuberculosis infection: diagnostic potential.

Mycobacterium avium subspecies paratuberculosis (Map) is the causative agent of Johne's disease (JD). Current serological diagnostic tests for JD are limited by their sensitivity when used in sub-clinical stages of the disease. Our objective was to identify peptides that mimic diagnostically important Map epitopes that might be incorporated into a new-generation JD diagnostic. Four peptides were isolated from a phage-displayed random peptide library by screening on antibodies derived from Map-infected goats. The peptides were recognised by antibodies from Map-infected goats but not by antibodies from uninfected goats. The peptides elicited immune responses in rabbits, which reacted strongly with bona fide Map antigens proving the peptides were true epitope mimics. To assess the diagnostic value a panel of goat sera was screened for reactivity's with peptides. The peptides were recognised by antibodies from a proportion of goats infected with Map compared with control animals with a diagnostic specificity of 100% and the sensitivity ranged from 50 to 75%. Combinations of any two peptides improved sensitivity 62.5-87.5% and 100% sensitivity was achieved with three of the four peptides in combination. These data suggest peptides representing diagnostically important Map epitopes could be incorporated into a sensitive diagnostic test.

How few whales were there after whaling? Inference from contemporary mtDNA diversity.

Reconstructing the history of exploited populations of whales requires fitting a trajectory through at least three points in time: (i) prior to exploitation, when abundance is assumed to be at the maximum allowed by environmental carrying capacity; (ii) the point of minimum abundance or 'bottleneck', usually near the time of protection or the abandonment of the hunt; and (iii) near the present, when protected populations are assumed to have undergone some recovery. As historical abundance is usually unknown, this trajectory must be extrapolated according to a population dynamic model using catch records, an assumed rate of increase and an estimate of current abundance, all of which have received considerable attention by the International Whaling Commission (IWC). Relatively little attention has been given to estimating minimum abundance (N(min)), although it is clear that genetic and demographic forces at this point are critical to the potential for recovery or extinction of a local population. We present a general analytical framework to improve estimates of N(min) using the number of mtDNA haplotypes (maternal lineages) surviving in a contemporary population of whales or other exploited species. We demonstrate the informative potential of this parameter as an a posteriori constraint on Bayesian logistic population dynamic models based on the IWC Comprehensive Assessment of the intensively exploited southern right whales (Eubalaena australis) and published surveys of mtDNA diversity for this species. Estimated historical trajectories from all demographic scenarios suggested a substantial loss of mtDNA haplotype richness as a result of 19th century commercial whaling and 20th century illegal whaling by the Soviet Union. However, the relatively high rates of population increase used by the IWC assessment predicted a bottleneck that was implausibly narrow (median, 67 mature females), given our corrected estimates of N(min). Further, high levels of remnant sequence diversity (theta) suggested that pre-exploitation abundance was larger than predicted by the logistic model given the catch record, which is known to be incomplete. Our results point to a need to better integrate evolutionary processes into population dynamic models to account for uncertainty in catch records, the influence of maternal fidelity on metapopulation dynamics, and the potential for inverse density dependence (an 'Allee effect') in severely depleted populations.