Strong and lasting impacts of past global warming on baleen whales and their prey.

Global warming is affecting the population dynamics and trophic interactions across a wide range of ecosystems and habitats. Translating these real-time effects into their long-term consequences remains a challenge. The rapid and extreme warming period that occurred after the Last Glacial Maximum (LGM) during the Pleistocene-Holocene transition (7-12 thousand years ago) provides an opportunity to gain insights into the long-term responses of natural populations to periods with global warming. The effects of this post-LGM warming period have been assessed in many terrestrial taxa, whereas insights into the impacts of rapid global warming on marine taxa remain limited, especially for megafauna. In order to understand how large-scale climate fluctuations during the post-LGM affected baleen whales and their prey, we conducted an extensive, large-scale analysis of the long-term effects of the post-LGM warming on abundance and inter-ocean connectivity in eight baleen whale and seven prey (fish and invertebrates) species across the Southern and the North Atlantic Ocean; two ocean basins that differ in key oceanographic features. The analysis was based upon 7032 mitochondrial DNA sequences as well as genome-wide DNA sequence variation in 100 individuals. The estimated temporal changes in genetic diversity during the last 30,000 years indicated that most baleen whale populations underwent post-LGM expansions in both ocean basins. The increase in baleen whale abundance during the Holocene was associated with simultaneous changes in their prey and climate. Highly correlated, synchronized and exponential increases in abundance in both baleen whales and their prey in the Southern Ocean were indicative of a dramatic increase in ocean productivity. In contrast, the demographic fluctuations observed in baleen whales and their prey in the North Atlantic Ocean were subtle, varying across taxa and time. Perhaps most important was the observation that the ocean-wide expansions and decreases in abundance that were initiated by the post-LGM global warming, continued for millennia after global temperatures stabilized, reflecting persistent, long-lasting impacts of global warming on marine fauna.

Return to the Sea, Get Huge, Beat Cancer: An Analysis of Cetacean Genomes Including an Assembly for the Humpback Whale (Megaptera novaeangliae).

Cetaceans are a clade of highly specialized aquatic mammals that include the largest animals that have ever lived. The largest whales can have ∼1,000× more cells than a human, with long lifespans, leaving them theoretically susceptible to cancer. However, large-bodied and long-lived animals do not suffer higher risks of cancer mortality than humans-an observation known as Peto's Paradox. To investigate the genomic bases of gigantism and other cetacean adaptations, we generated a de novo genome assembly for the humpback whale (Megaptera novaeangliae) and incorporated the genomes of ten cetacean species in a comparative analysis. We found further evidence that rorquals (family Balaenopteridae) radiated during the Miocene or earlier, and inferred that perturbations in abundance and/or the interocean connectivity of North Atlantic humpback whale populations likely occurred throughout the Pleistocene. Our comparative genomic results suggest that the evolution of cetacean gigantism was accompanied by strong selection on pathways that are directly linked to cancer. Large segmental duplications in whale genomes contained genes controlling the apoptotic pathway, and genes inferred to be under accelerated evolution and positive selection in cetaceans were enriched for biological processes such as cell cycle checkpoint, cell signaling, and proliferation. We also inferred positive selection on genes controlling the mammalian appendicular and cranial skeletal elements in the cetacean lineage, which are relevant to extensive anatomical changes during cetacean evolution. Genomic analyses shed light on the molecular mechanisms underlying cetacean traits, including gigantism, and will contribute to the development of future targets for human cancer therapies.

Fin whale (Balaenoptera physalus) mitogenomics: A cautionary tale of defining sub-species from mitochondrial sequence monophyly.

The advent of massive parallel sequencing technologies has resulted in an increase of studies based upon complete mitochondrial genome DNA sequences that revisit the taxonomic status within and among species. Spatially distinct monophyly in such mitogenomic genealogies, i.e., the sharing of a recent common ancestor among con-specific samples collected in the same region has been viewed as evidence for subspecies. Several recent studies in cetaceans have employed this criterion to suggest subsequent intraspecific taxonomic revisions. We reason that employing intra-specific, spatially distinct monophyly at non-recombining, clonally inherited genomes is an unsatisfactory criterion for defining subspecies based upon theoretical (genetic drift) and practical (sampling effort) arguments. This point was illustrated by a re-analysis of a global mitogenomic assessment of fin whales, Balaenoptera physalus spp., published by Archer et al. (2013), which proposed to further subdivide the Northern Hemisphere fin whale subspecies, B. p. physalus. The proposed revision was based upon the detection of spatially distinct monophyly among North Atlantic and North Pacific fin whales in a genealogy based upon complete mitochondrial genome DNA sequences. The extended analysis conducted in this study (1676 mitochondrial control region, 162 complete mitochondrial genome DNA sequences and 20 microsatellite loci genotyped in 380 samples) revealed that the apparent monophyly among North Atlantic fin whales reported by Archer et al. (2013) to be due to low sample sizes. In conclusion, defining sub-species from monophyly (i.e., the absence of para- or polyphyly) can lead to erroneous conclusions due to relatively "trivial" aspects, such as sampling. Basic population genetic processes (i.e., genetic drift and migration) also affect the time to the most recent common ancestor and hence the probability that individuals in a sample are monophyletic.

Whales in warming water: Assessing breeding habitat diversity and adaptability in Oceania's changing climate.

In the context of a changing climate, understanding the environmental drivers of marine megafauna distribution is important for conservation success. The extent of humpback whale breeding habitats and the impact of temperature variation on their availability are both unknown. We used 19 years of dedicated survey data from seven countries and territories of Oceania (1,376 survey days), to investigate humpback whale breeding habitat diversity and adaptability to climate change. At a fine scale (1 km resolution), seabed topography was identified as an important influence on humpback whale distribution. The shallowest waters close to shore or in lagoons were favored, although humpback whales also showed flexible habitat use patterns with respect to shallow offshore features such as seamounts. At a coarse scale (1° resolution), humpback whale breeding habitats in Oceania spanned a thermal range of 22.3-27.8°C in August, with interannual variation up to 2.0°C. Within this range, both fine and coarse scale analyses of humpback whale distribution suggested local responses to temperature. Notably, the most detailed dataset was available from New Caledonia (774 survey days, 1996-2017), where encounter rates showed a negative relationship to sea surface temperature, but were not related to the El Niño Southern Oscillation or the Antarctic Oscillation from previous summer, a proxy for feeding conditions that may impact breeding patterns. Many breeding sites that are currently occupied are predicted to become unsuitably warm for this species (>28°C) by the end of the 21st century. Based on modeled ecological relationships, there are suitable habitats for relocation in archipelagos and seamounts of southern Oceania. Although distribution shifts might be restrained by philopatry, the apparent plasticity of humpback whale habitat use patterns and the extent of suitable habitats support an adaptive capacity to ocean warming in Oceania breeding grounds.

Evaluation of fecal hormones for noninvasive research on reproduction and stress in humpback whales (Megaptera novaeangliae).

Fecal hormone analysis shows high potential for noninvasive assessment of population-level patterns in stress and reproduction of marine mammals. However, the marine environment presents unique challenges for fecal sample collection. Data are still lacking on collection methodology and assay validations for most species, particularly for those mysticete whales that have variable diets. In this study we tested collection techniques for fecal samples of free-swimming humpback whales (Megaptera novaeangliae), and validated immunoassays for five steroid and thyroid hormones. Resulting data were used for preliminary physiological validations, i.e., comparisons to independently confirmed sex and reproductive state. Pregnant females had significantly higher fecal progestins and glucocorticoids than did other demographic categories of whales. Two possible cases of previously undetected pregnancies were noted. Males had significantly higher fecal testosterone metabolites than nonpregnant females. Fecal glucocorticoids were significantly elevated in pregnant females and mature males compared to nonpregnant females. Calf fecal samples had elevated concentrations of all fecal hormones. Fecal thyroid hormones showed a significant seasonal decline from spring to summer. Though sample sizes were small, and sampling was necessarily opportunistic, these patterns indicate that noninvasive fecal hormone analysis may facilitate studies of reproduction, stress and potentially energetics in humpback whales.

Effects of fishing rope strength on the severity of large whale entanglements.

Entanglement in fixed fishing gear affects whales worldwide. In the United States, deaths of North Atlantic right (Eubalaena glacialis) and humpback whales (Megaptera novaeangliae) have exceeded management limits for decades. We examined live and dead whales entangled in fishing gear along the U.S. East Coast and the Canadian Maritimes from 1994 to 2010. We recorded whale species, age, and injury severity and determined rope polymer type, breaking strength, and diameter of the fishing gear. For the 132 retrieved ropes from 70 cases, tested breaking strength range was 0.80-39.63 kN (kiloNewtons) and the mean was 11.64 kN (SD 8.29), which is 26% lower than strength at manufacture (range 2.89-53.38 kN, mean = 15.70 kN [9.89]). Median rope diameter was 9.5 mm. Right and humpback whales were found in ropes with significantly stronger breaking strengths at time of manufacture than minke whales (Balaenoptera acuturostrata) (19.30, 17.13, and 10.47 mean kN, respectively). Adult right whales were found in stronger ropes (mean 34.09 kN) than juvenile right whales (mean 15.33 kN) and than all humpback whale age classes (mean 17.37 kN). For right whales, severity of injuries increased since the mid 1980s, possibly due to changes in rope manufacturing in the mid 1990s that resulted in production of stronger ropes at the same diameter. Our results suggest that broad adoption of ropes with breaking strengths of ≤ 7.56 kN (≤ 1700 lbsf) could reduce the number of life-threatening entanglements for large whales by at least 72%, and yet could provide sufficient strength to withstand the routine forces involved in many fishing operations. A reduction of this magnitude would achieve nearly all the mitigation legally required for U.S. stocks of North Atlantic right and humpback whales. Ropes with reduced breaking strength should be developed and tested to determine the feasibility of their use in a variety of fisheries.

Molecular biomarkers for chronological age in animal ecology.

The chronological age of an individual animal predicts many of its biological characteristics, and these in turn influence population-level ecological processes. Animal age information can therefore be valuable in ecological research, but many species have no external features that allow age to be reliably determined. Molecular age biomarkers provide a potential solution to this problem. Research in this area of molecular ecology has so far focused on a limited range of age biomarkers. The most commonly tested molecular age biomarker is change in average telomere length, which predicts age well in a small number of species and tissues, but performs poorly in many other situations. Epigenetic regulation of gene expression has recently been shown to cause age-related modifications to DNA and to cause changes in abundance of several RNA types throughout animal lifespans. Age biomarkers based on these epigenetic changes, and other new DNA-based assays, have already been applied to model organisms, humans and a limited number of wild animals. There is clear potential to apply these marker types more widely in ecological studies. For many species, these new approaches will produce age estimates where this was previously impractical. They will also enable age information to be gathered in cross-sectional studies and expand the range of demographic characteristics that can be quantified with molecular methods. We describe the range of molecular age biomarkers that have been investigated to date and suggest approaches for developing the newer marker types as age assays in nonmodel animal species.

Population structure of humpback whales in the western and central South Pacific Ocean as determined by vocal exchange among populations.

For cetaceans, population structure is traditionally determined by molecular genetics or photographically identified individuals. Acoustic data, however, has provided information on movement and population structure with less effort and cost than traditional methods in an array of taxa. Male humpback whales (Megaptera novaeangliae) produce a continually evolving vocal sexual display, or song, that is similar among all males in a population. The rapid cultural transmission (the transfer of information or behavior between conspecifics through social learning) of different versions of this display between distinct but interconnected populations in the western and central South Pacific region presents a unique way to investigate population structure based on the movement dynamics of a song (acoustic) display. Using 11 years of data, we investigated an acoustically based population structure for the region by comparing stereotyped song sequences among populations and years. We used the Levenshtein distance technique to group previously defined populations into (vocally based) clusters based on the overall similarity of their song display in space and time. We identified the following distinct vocal clusters: western cluster, 1 population off eastern Australia; central cluster, populations around New Caledonia, Tonga, and American Samoa; and eastern region, either a single cluster or 2 clusters, one around the Cook Islands and the other off French Polynesia. These results are consistent with the hypothesis that each breeding aggregation represents a distinct population (each occupied a single, terminal node) in a metapopulation, similar to the current understanding of population structure based on genetic and photo-identification studies. However, the central vocal cluster had higher levels of song-sharing among populations than the other clusters, indicating that levels of vocal connectivity varied within the region. Our results demonstrate the utility and value of using culturally transmitted vocal patterns as a way of defining connectivity to infer population structure. We suggest vocal patterns be incorporated by the International Whaling Commission in conjunction with traditional methods in the assessment of structure.

Estructura Poblacional de las Yubartas en el Centro y Occidente del Sur del Océano Pacífico, Determinada por el Intercambio Vocal entre las Poblaciones Resumen La estructura poblacional de los cetáceos se determina tradicionalmente por medio de genética molecular o individuos identificados fotográficamente. Sin embargo, los datos acústicos han proporcionado información sobre el movimiento y la estructura poblacional de una variedad de taxones con menos esfuerzo y menos costo que los métodos tradicionales. Las yubartas macho (Megaptera novaeangliae) producen una demostración vocal sexual, o canto, en continua evolución que es similar entre todos los machos en una población. La transmisión cultural rápida (la transferencia de información o comportamientos entre co-específicos por medio del aprendizaje social) de diferentes versiones de esta demostración entre poblaciones distintas pero interconectadas en el centro y occidente en la región del sur del Pacífico presenta una manera única para investigar la estructura poblacional con base en las dinámicas de movimiento de la demostración del canto (acústico). Al usar once años de datos, con base en la acústica investigamos a una estructura poblacional de la región al comparar secuencias estereotipadas de cantos entre las poblaciones y los años. Usamos la técnica de la distancia Levenshtein para agrupar poblaciones previamente definidas en grupos (con base vocal) basados en la similitud general de su canto en el espacio y el tiempo. Identificamos a los siguientes grupos vocales perceptibles: grupo occidental, una población al este de Australia; grupo central, poblaciones alrededor de Nueva Caledonia, Tonga y Samoa Americana; y región oriental, ya sea un grupo único o dos grupos, uno alrededor de las Islas Cook y el otro cerca de la Polinesia Francesa. Estos resultados son consistentes con la hipótesis de que cada conjunto reproductivo representa a una población distinta (cada una ocupó un nodo terminal único) en una meta-población, similar al entendimiento actual de la estructura poblacional con base en los estudios genéticos y de foto-identificación. Sin embargo, el grupo vocal central tuvo niveles más altos de cantos compartidos entre las poblaciones que en los otros grupos, lo que indica que los niveles de conectividad vocal variaron dentro de la región. Nuestros resultados demuestran la utilidad y el valor de usar patrones vocales transmitidos culturalmente como una forma de definir la conectividad para inferir la estructura poblacional. Sugerimos que los patrones vocales sean incorporados por la Comisión Internacional de Caza de Ballenas en conjunto con los métodos tradicionales de evaluación de la estructura.

Caller identification and characterization of individual humpback whale acoustic behaviour.

Acoustic recording tags provide fine-scale data linking acoustic signalling with individual behaviour; however, when an animal is in a group, it is challenging to tease apart calls of conspecifics and identify which individuals produce each call. This, in turn, prohibits a robust assessment of individual acoustic behaviour including call rates and silent periods, call bout production within and between individuals, and caller location. To overcome this challenge, we simultaneously instrumented small groups of humpback whales on a western North Atlantic feeding ground with sound and movement recording tags. This approach enabled a comparison of the relative amplitude of each call across individuals to infer caller identity for 97% of calls. We recorded variable call rates across individuals (mean = 23 calls/h) and groups (mean = 55 calls/h). Calls were produced throughout dives, and most calls were produced in bouts with short inter-call intervals of 2.2 s. Most calls received a likely response from a conspecific within 100 s. This caller identification (ID) method facilitates studying both individual- and group-level acoustic behaviour, yielding novel results about the nature of sequence production and vocal exchanges in humpback whale social calls. Future studies can expand on these caller ID methods for understanding intra-group communication across taxa.

Acoustic behavior of humpback whale calves on the feeding ground: Comparisons across age and implications for vocal development.

Studying sound production at different developmental stages can provide insight into the processes involved in vocal ontogeny. Humpback whales (Megaptera novaeangliae) are a known vocal learning species, but their vocal development is poorly understood. While studies of humpback whale calves in the early stages of their lives on the breeding grounds and migration routes exist, little is known about the behavior of these immature, dependent animals by the time they reach the feeding grounds. In this study, we used data from groups of North Atlantic humpback whales in the Gulf of Maine in which all members were simultaneously carrying acoustic recording tags attached with suction cups. This allowed for assignment of likely caller identity using the relative received levels of calls across tags. We analyzed data from 3 calves and 13 adults. There were high levels of call rate variation among these individuals and the results represent preliminary descriptions of calf behavior. Our analysis suggests that, in contrast to the breeding grounds or on migration, calves are no longer acoustically cryptic by the time they reach their feeding ground. Calves and adults both produce calls in bouts, but there may be some differences in bout parameters like inter-call intervals and bout durations. Calves were able to produce most of the adult vocal repertoire but used different call types in different proportions. Finally, we found evidence of immature call types in calves, akin to protosyllables used in babbling in other mammals, including humans. Overall, the sound production of humpback whale calves on the feeding grounds appears to be already similar to that of adults, but with differences in line with ontogenetic changes observed in other vocal learning species.

Quantifying humpback whale song sequences to understand the dynamics of song exchange at the ocean basin scale.

Humpback whales have a continually evolving vocal sexual display, or "song," that appears to undergo both evolutionary and "revolutionary" change. All males within a population adhere to the current content and arrangement of the song. Populations within an ocean basin share similarities in their songs; this sharing is complex as multiple variations of the song (song types) may be present within a region at any one time. To quantitatively investigate the similarity of song types, songs were compared at both the individual singer and population level using the Levenshtein distance technique and cluster analysis. The highly stereotyped sequences of themes from the songs of 211 individuals from populations within the western and central South Pacific region from 1998 through 2008 were grouped together based on the percentage of song similarity, and compared to qualitatively assigned song types. The analysis produced clusters of highly similar songs that agreed with previous qualitative assignments. Each cluster contained songs from multiple populations and years, confirming the eastward spread of song types and their progressive evolution through the study region. Quantifying song similarity and exchange will assist in understanding broader song dynamics and contribute to the use of vocal displays as population identifiers.

Persistent organic pollutants in female humpback whales Megaptera novaeangliae from the Gulf of Maine.

Contaminant studies in cetaceans can provide information about pollutant levels and patterns in a given region. Due to the confounding effects of reproductive status and maternal offloading in females, these studies typically focus on males. However, an improved understanding of contaminant burdens in female cetaceans is needed to better assess potential impacts to populations. The objectives of this study were to characterize concentrations of persistent organic pollutants (POPs) in blubber of female humpback whales across age classes and to also better characterize maternal offloading of these pollutants to their offspring. A total of 36 blubber biopsy samples of female humpback whales (Megaptera novaeangliae) from the Gulf of Maine were analyzed to examine contaminant loads across females of different ages. Sampled individuals were individually-identified from longitudinal studies and assigned to age class (i.e., adult, subadult, juvenile, calf). Analysis was performed using gas chromatography/mass spectrometry (GC/MS) of POPs including polychlorinated biphenyls (PCBs), dichlorodiphenyltrichloroethanes (DDTs), chlordanes (CHLDs), polybrominated diphenyl ethers (PBDEs), hexachlorocyclohexanes (HCHs). The most abundant POPs were PCB congeners, with summed values ranging from 280 to 12,000 ng/g, lipid weight, which is above recent estimates of the threshold for adverse health effects. We found significant differences in mean values between adults and juveniles and between adults and subadults, with the exception of the less persistent HCHs for the latter. We also found significant differences in mean levels of ∑HCHs between the juveniles and subadults. Changes over age are consistent with maternal offloading and potentially important for evaluating population health and viability.

PHFinder: assisted detection of point heteroplasmy in Sanger sequencing chromatograms.

Heteroplasmy is the presence of two or more organellar genomes (mitochondrial or plastid DNA) in an organism, tissue, cell or organelle. Heteroplasmy can be detected by visual inspection of Sanger sequencing chromatograms, where it appears as multiple peaks of fluorescence at a single nucleotide position. Visual inspection of chromatograms is both consuming and highly subjective, as heteroplasmy is difficult to differentiate from background noise. Few software solutions are available to automate the detection of point heteroplasmies, and those that are available are typically proprietary, lack customization or are unsuitable for automated heteroplasmy assessment in large datasets. Here, we present PHFinder, a Python-based, open-source tool to assist in the detection of point heteroplasmies in large numbers of Sanger chromatograms. PHFinder automatically identifies point heteroplasmies directly from the chromatogram trace data. The program was tested with Sanger sequencing data from 100 humpback whales (Megaptera novaeangliae) tissue samples with known heteroplasmies. PHFinder detected most (90%) of the known heteroplasmies thereby greatly reducing the amount of visual inspection required. PHFinder is flexible and enables explicit specification of key parameters to infer double peaks (i.e., heteroplasmies).

Wild pedigrees inform mutation rates and historic abundance in baleen whales.

Phylogeny-based estimates suggesting a low germline mutation rate (µ) in baleen whales have influenced research ranging from assessments of whaling impacts to evolutionary cancer biology. We estimated µ directly from pedigrees in four baleen whale species for both the mitochondrial control region and nuclear genome. The results suggest values higher than those obtained through phylogeny-based estimates and similar to pedigree-based values for primates and toothed whales. Applying our estimate of µ reduces previous genetic-based estimates of preexploitation whale abundance by 86% and suggests that µ cannot explain low cancer rates in gigantic mammals. Our study shows that it is feasible to estimate µ directly from pedigrees in natural populations, with wide-ranging implications for ecological and evolutionary research.

Vulnerability to climate change of United States marine mammal stocks in the western North Atlantic, Gulf of Mexico, and Caribbean.

Climate change and climate variability are affecting marine mammal species and these impacts are projected to continue in the coming decades. Vulnerability assessments provide a framework for evaluating climate impacts over a broad range of species using currently available information. We conducted a trait-based climate vulnerability assessment using expert elicitation for 108 marine mammal stocks and stock groups in the western North Atlantic, Gulf of Mexico, and Caribbean Sea. Our approach combined the exposure (projected change in environmental conditions) and sensitivity (ability to tolerate and adapt to changing conditions) of marine mammal stocks to estimate vulnerability to climate change, and categorize stocks with a vulnerability index. The climate vulnerability score was very high for 44% (n = 47) of these stocks, high for 29% (n = 31), moderate for 20% (n = 22), and low for 7% (n = 8). The majority of stocks (n = 78; 72%) scored very high exposure, whereas 24% (n = 26) scored high, and 4% (n = 4) scored moderate. The sensitivity score was very high for 33% (n = 36) of these stocks, high for 18% (n = 19), moderate for 34% (n = 37), and low for 15% (n = 16). Vulnerability results were summarized for stocks in five taxonomic groups: pinnipeds (n = 4; 25% high, 75% moderate), mysticetes (n = 7; 29% very high, 57% high, 14% moderate), ziphiids (n = 8; 13% very high, 50% high, 38% moderate), delphinids (n = 84; 52% very high, 23% high, 15% moderate, 10% low), and other odontocetes (n = 5; 60% high, 40% moderate). Factors including temperature, ocean pH, and dissolved oxygen were the primary drivers of high climate exposure, with effects mediated through prey and habitat parameters. We quantified sources of uncertainty by bootstrapping vulnerability scores, conducting leave-one-out analyses of individual attributes and individual scorers, and through scoring data quality for each attribute. These results provide information for researchers, managers, and the public on marine mammal responses to climate change to enhance the development of more effective marine mammal management, restoration, and conservation activities that address current and future environmental variation and biological responses due to climate change.

Empirical evaluation of humpback whale telomere length estimates; quality control and factors causing variability in the singleplex and multiplex qPCR methods.

BACKGROUND: Telomeres, the protective cap of chromosomes, have emerged as powerful markers of biological age and life history in model and non-model species. The qPCR method for telomere length estimation is one of the most common methods for telomere length estimation, but has received recent critique for being too error-prone and yielding unreliable results. This critique coincides with an increasing awareness of the potentials and limitations of the qPCR technique in general and the proposal of a general set of guidelines (MIQE) for standardization of experimental, analytical, and reporting steps of qPCR. In order to evaluate the utility of the qPCR method for telomere length estimation in non-model species, we carried out four different qPCR assays directed at humpback whale telomeres, and subsequently performed a rigorous quality control to evaluate the performance of each assay. RESULTS: Performance differed substantially among assays and only one assay was found useful for telomere length estimation in humpback whales. The most notable factors causing these inter-assay differences were primer design and choice of using singleplex or multiplex assays. Inferred amplification efficiencies differed by up to 40% depending on assay and quantification method, however this variation only affected telomere length estimates in the worst performing assays. CONCLUSION: Our results suggest that seemingly well performing qPCR assays may contain biases that will only be detected by extensive quality control. Moreover, we show that the qPCR method for telomere length estimation can be highly precise and accurate, and thus suitable for telomere measurement in non-model species, if effort is devoted to optimization at all experimental and analytical steps. We conclude by highlighting a set of quality controls which may serve for further standardization of the qPCR method for telomere length estimation, and discuss some of the factors that may cause variation in qPCR experiments.

How many genetic markers to tag an individual? An empirical assessment of false matching rates among close relatives.

Genetic identification of individuals is now commonplace, enabling the application of tagging methods to elusive species or species that cannot be tagged by traditional methods. A key aspect is determining the number of loci required to ensure that different individuals have non-matching multi-locus genotypes. Closely related individuals are of particular concern because of elevated matching probabilities caused by their recent co-ancestry. This issue may be addressed by increasing the number of loci to a level where full siblings (the relatedness category with the highest matching probability) are expected to have non-matching multi-locus genotypes. However, increasing the number of loci to meet this "full-sib criterion" greatly increases the laboratory effort, which in turn may increase the genotyping error rate resulting in an upward-biased mark-recapture estimate of abundance as recaptures are missed due to genotyping errors. We assessed the contribution of false matches from close relatives among 425 maternally related humpback whales, each genotyped at 20 microsatellite loci. We observed a very low (0.5-4%) contribution to falsely matching samples from pairs of first-order relatives (i.e., parent and offspring or full siblings). The main contribution to falsely matching individuals from close relatives originated from second-order relatives (e.g., half siblings), which was estimated at 9%. In our study, the total number of observed matches agreed well with expectations based upon the matching probability estimated for unrelated individuals, suggesting that the full-sib criterion is overly conservative, and would have required a 280% relative increase in effort. We suggest that, under most circumstances, the overall contribution to falsely matching samples from close relatives is likely to be low, and hence applying the full-sib criterion is unnecessary. In those cases where close relatives may present a significant issue, such as unrepresentative sampling, we propose three different genotyping strategies requiring only a modest increase in effort, which will greatly reduce the number of false matches due to the presence of related individuals.

Multi-year progesterone profiles during pregnancy in baleen of humpback whales (Megaptera novaeangliae).

Understanding calving rates of wild whale populations is critically important for management and conservation. Reproduction of humpback whales (Megaptera novaeangliae) is difficult to monitor and, even with long-term sighting studies, basic physiological information such as pregnancy rates and calving intervals remain poorly understood in many populations. We hypothesized that pregnant whales have sustained elevations in baleen progesterone that temporally correlate with gestation. To test this hypothesis, baleen progesterone profiles from two adult female North Pacific humpbacks, both with extensive sighting records and documented pregnancies, were compared to those of a nulliparous female (adult female never seen with a calf) and a juvenile male. Baleen specimens recovered during necropsy were subsampled every 2 cm from the base to the tip of the plate, with each interval representing 30-45 days of growth. Homogenized baleen powder was assayed for progesterone using enzyme immunoassays. The date of growth of each sampling location on the baleen plate was estimated based on stable isotope analysis of annual δ15N cycles. Progesterone profiles from both pregnant whales showed sustained high progesterone content (>350 ng/g) in areas corresponding to known pregnancies, inferred from calf sightings and post-mortem data. The younger female, estimated to be 13 years old, had higher progesterone during pregnancy than the 44.5 year old, but levels during non-pregnancy were similar. The nulliparous female and the male had low progesterone throughout their baleen plates. Baleen hormone analysis can determine how progesterone concentrations change throughout gestation and has potential for estimating age at first reproduction, pregnancy intervals, failed pregnancies and early calf mortality. Understanding rates of calving and current and historic reproductive patterns in humpbacks is vital to continuing conservation measures in this species.

Patterns of cortisol and corticosterone concentrations in humpback whale (Megaptera novaeangliae) baleen are associated with different causes of death.

Baleen whales are subject to a myriad of natural and anthropogenic stressors, but understanding how these stressors affect physiology is difficult. Measurement of adrenal glucocorticoid (GC) hormones involved in the vertebrate stress response (cortisol and corticosterone) in baleen could help fill this data gap. Baleen analysis is a powerful tool, allowing for a retrospective re-creation of multiple years of GC hormone concentrations at approximately a monthly resolution. We hypothesized that whales that died from acute causes (e.g. ship strike) would have lower levels of baleen GCs than whales that died from extended illness or injury (e.g. long-term entanglement in fishing gear). To test this hypothesis, we extracted hormones from baleen plates of four humpback whales (Megaptera novaeangliae) with well-documented deaths including multiple and chronic entanglements (n = 1, female), ship strike (n = 2, male and female) and chronic illness with nutritional stress (n = 1, male). Over ~3 years of baleen growth and during multiple entanglements, the entangled whale had average corticosterone levels of 80-187% higher than the other three whales but cortisol levels were similar to two of the other three whales. The nutritionally stressed and chronically ill whale showed a slow increase in both cortisol and corticosterone spanning ~3 years, followed by a sharp decline in both hormones before death, possibly indicative of adrenal failure in this moribund individual. This whale's correlation between cortisol and corticosterone was significant but there were no correlations in the other three whales. Our results show that cortisol and corticosterone concentrations vary according to the type and duration of illness or injury. Single-point GC concentrations should be interpreted with caution as low values can occur in whales experiencing pronounced stress and individual baselines can be highly variable. Baleen analysis is a promising tissue type for retrospective analyses of physiological responses to various stressors affecting baleen whales.

Characterizing the culturable surface microbiomes of diverse marine animals.

Biofilm-forming bacteria have the potential to contribute to the health, physiology, behavior and ecology of the host and serve as its first line of defense against adverse conditions in the environment. While metabarcoding and metagenomic information furthers our understanding of microbiome composition, fewer studies use cultured samples to study the diverse interactions among the host and its microbiome, as cultured representatives are often lacking. This study examines the surface microbiomes cultured from three shallow-water coral species and two whale species. These unique marine animals place strong selective pressures on their microbial symbionts and contain members under similar environmental and anthropogenic stress. We developed an intense cultivation procedure, utilizing a suite of culture conditions targeting a rich assortment of biofilm-forming microorganisms. We identified 592 microbial isolates contained within 15 bacterial orders representing 50 bacterial genera, and two fungal species. Culturable bacteria from coral and whale samples paralleled taxonomic groups identified in culture-independent surveys, including 29% of all bacterial genera identified in the Megaptera novaeangliae skin microbiome through culture-independent methods. This microbial repository provides raw material and biological input for more nuanced studies which can explore how members of the microbiome both shape their micro-niche and impact host fitness.