Elucidating the sustainability of 700 y of Inuvialuit beluga whale hunting in the Mackenzie River Delta, Northwest Territories, Canada.

Beluga whales play a critical role in the subsistence economies and cultural heritage of Indigenous communities across the Arctic, yet the effects of Indigenous hunting on beluga whales remain unknown. Here, we integrate paleogenomics, genetic simulations, and stable δ13C and δ15N isotope analysis to investigate 700 y of beluga subsistence hunting in the Mackenzie Delta area of northwestern Canada. Genetic identification of the zooarchaeological remains, which is based on radiocarbon dating, span three time periods (1290 to 1440 CE; 1450 to 1650 CE; 1800 to 1870 CE), indicates shifts across time in the sex ratio of the harvested belugas. The equal number of females and males harvested in 1450 to 1650 CE versus more males harvested in the two other time periods may reflect changes in hunting practices or temporal shifts in beluga availability. We find temporal shifts and sex-based differences in δ13C of the harvested belugas across time, suggesting historical adaptability in the foraging ecology of the whales. We uncovered distinct mitochondrial diversity unique to the Mackenzie Delta belugas, but found no changes in nuclear genomic diversity nor any substructuring across time. Our findings indicate the genomic stability and continuity of the Mackenzie Delta beluga population across the 700 y surveyed, indicating the impact of Inuvialuit subsistence harvests on the genetic diversity of contemporary beluga individuals has been negligible.

Characterization of a novel papillomavirus identified from a whale (Delphinapterus leucas) pharyngeal metagenomic library.

Here, using viral metagenomic method, a novel whale papillomavirus (temporarily named wPV, GenBank accession number OP856597) was discovered in a whale (Delphinapterus leucas) pharyngeal metagenomic library. The complete genome size of wPV is 7179 bp, with GC content of 54.4% and a nucleotide composition of 23.4% A, 22.3% T, 28.4% G, and 25.9% C. The viral genome has a typical papillomavirus organization pattern, and five ORFs were predicted, including two late genes encoding L1 and L2, and three early genes encoding E1, E2, and E6. Pairwise sequence comparison and phylogenetic analysis based on the L1 gene sequence indicated that wPV may be a novel species within genus Dyodeltapapillomavirus. In addition, the E2 region of wPV was predicted to have a potential recombination event. The discovery of this novel papillomavirus increases our understanding of the viral ecology of marine mammals, providing insights into possible future infectious diseases.

Circumpolar phylogeography and demographic history of beluga whales reflect past climatic fluctuations.

Several Arctic marine mammal species are predicted to be negatively impacted by rapid sea ice loss associated with ongoing ocean warming. However, consequences for Arctic whales remain uncertain. To investigate how Arctic whales responded to past climatic fluctuations, we analysed 206 mitochondrial genomes from beluga whales (Delphinapterus leucas) sampled across their circumpolar range, and four nuclear genomes, covering both the Atlantic and the Pacific Arctic region. We found four well-differentiated mitochondrial lineages, which were established before the onset of the last glacial expansion ~110 thousand years ago. Our findings suggested these lineages diverged in allopatry, reflecting isolation of populations during glacial periods when the Arctic sea-shelf was covered by multiyear sea ice. Subsequent population expansion and secondary contact between the Atlantic and Pacific Oceans shaped the current geographic distribution of lineages, and may have facilitated mitochondrial introgression. Our demographic reconstructions based on both mitochondrial and nuclear genomes showed markedly lower population sizes during the Last Glacial Maximum (LGM) compared to the preceding Eemian and current Holocene interglacial periods. Habitat modelling similarly revealed less suitable habitat during the LGM (glacial) than at present (interglacial). Together, our findings suggested the association between climate, population size, and available habitat in belugas. Forecasts for year 2100 showed that beluga habitat will decrease and shift northwards as oceans continue to warm, putatively leading to population declines in some beluga populations. Finally, we identified vulnerable populations which, if extirpated as a consequence of ocean warming, will lead to a substantial decline of species-wide haplotype diversity.

A comparative analysis of cone photoreceptor morphology in bowhead and beluga whales.

The cetacean visual system is a product of selection pressures favoring underwater vision, yet relatively little is known about it across taxa. Previous studies report several mutations in the opsin genetic sequence in cetaceans, suggesting the evolutionary complete or partial loss of retinal cone photoreceptor function in mysticete and odontocete lineages, respectively. Despite this, limited anatomical evidence suggests cone structures are partially maintained but with absent outer and inner segments in the bowhead retina. The functional consequence and anatomical distributions associated with these unique cone morphologies remain unclear. The current study further investigates the morphology and distribution of cone photoreceptors in the bowhead whale and beluga retina and evaluates the potential functional capacity of these cells' alternative to photoreception. Refined histological and advanced microscopic techniques revealed two additional cone morphologies in the bowhead and beluga retina that have not been previously described. Two proteins involved in magnetosensation were present in these cone structures suggesting the possibility for an alternative functional role in responding to changes in geomagnetic fields. These findings highlight a revised understanding of the unique evolution of cone and gross retinal anatomy in cetaceans, and provide prefatory evidence of potential functional reassignment of these cells.

IMMUNOMEDIATOR GENE TRANSCRIPTION PROFILING IN BELUGA WHALE (DELPHINAPTERUS LEUCAS) CLINICAL CASES.

There is an unmet need for specific diagnostics of immune perturbations and inflammation in beluga whale (Delphinapterus leucas) clinical care. Quantitative real-time polymerase chain reaction (qPCR) has been used to measure immunomediator gene transcription in beluga whales. The study hypothesis was that a qPCR-based immunomediator assay would supplement routine clinical data with specific and sensitive information on immune status. Two beluga whale clinical cases provided an opportunity to test this hypothesis: a whale with a skin laceration and a whale with gastrointestinal inflammation. Mitogen-stimulated immunomediator gene transcription (MSIGT) was compared between the cases and healthy contact whales. In both case studies, mitogens increased transcription of IL1B, PTGS2 (Cox-2), TNF, HIF1A, and IL2 but decreased IL10 transcription in peripheral blood mononuclear cells (PBMC) from the abnormal whale over the control. Correlations were identified between most immunomediators tested and one or more standard blood clinical values. Considering all 15 immunomediators tested, the whale with gastrointestinal inflammation had a more unique MSIGT signature than the whale with a laceration. These results support further elucidation of beluga whale PBMC cytokine profiles for use as immune biomarkers.

Associations between organohalogen exposure and thyroid- and steroid-related gene responses in St. Lawrence Estuary belugas and minke whales.

Elevated concentrations of persistent organic pollutants (POPs) and emerging halogenated flame retardants (HFRs) have been reported in tissues of the endangered St. Lawrence Estuary (Canada) beluga population as well as in minke whales visiting that same feeding area. This study examined the linkages between blubber concentrations of POPs and emerging HFRs, and transcription in skin of genes involved in the regulation of thyroid and steroid axes in belugas and minke whales from the St. Lawrence Estuary. In belugas, concentrations of PCBs, OCs and hexabromobenzene (HBB) were positively correlated with the transcription of thyroid- and/or steroid-related genes, while Dec-604 CB concentrations were negatively associated with the transcription of glucocorticoid and thyroid genes. In minke whales, PBDE concentrations changed positively with Esrß transcript levels and HBB concentrations negatively with Nr3c1 transcripts. Present results suggest that several biological functions including reproduction and energetic metabolism may represent potential targets for organohalogens in these whales.

De novo transcriptome assembly and RNA-Seq expression analysis in blood from beluga whales of Bristol Bay, AK.

Assessing the health of marine mammal sentinel species is crucial to understanding the impacts of environmental perturbations on marine ecosystems and human health. In Arctic regions, beluga whales, Delphinapterus leucas, are upper level predators that may serve as a sentinel species, potentially forecasting impacts on human health. While gene expression profiling from blood transcriptomes has widely been used to assess health status and environmental exposures in human and veterinary medicine, its use in wildlife has been limited due to the lack of available genomes and baseline data. To this end we constructed the first beluga whale blood transcriptome de novo from samples collected during annual health assessments of the healthy Bristol Bay, AK stock during 2012-2014 to establish baseline information on the content and variation of the beluga whale blood transcriptome. The Trinity transcriptome assembly from beluga was comprised of 91,325 transcripts that represented a wide array of cellular functions and processes and was extremely similar in content to the blood transcriptome of another cetacean, the bottlenose dolphin. Expression of hemoglobin transcripts was much lower in beluga (25.6% of TPM, transcripts per million) than has been observed in many other mammals. A T12A amino acid substitution in the HBB sequence of beluga whales, but not bottlenose dolphins, was identified and may play a role in low temperature adaptation. The beluga blood transcriptome was extremely stable between sex and year, with no apparent clustering of samples by principle components analysis and <4% of genes differentially expressed (EBseq, FDR<0.05). While the impacts of season, sexual maturity, disease, and geography on the beluga blood transcriptome must be established, the presence of transcripts involved in stress, detoxification, and immune functions indicate that blood gene expression analyses may provide information on health status and exposure. This study provides a wealth of transcriptomic data on beluga whales and provides a sizeable pool of preliminary data for comparison with other studies in beluga whale.

PCBs are associated with altered gene transcript profiles in arctic Beluga Whales (Delphinapterus leucas).

High trophic level arctic beluga whales (Delphinapterus leucas) are exposed to persistent organic pollutants (POP) originating primarily from southern latitudes. We collected samples from 43 male beluga harvested by Inuvialuit hunters (2008-2010) in the Beaufort Sea to evaluate the effects of POPs on the levels of 13 health-related gene transcripts using quantitative real-time polymerase chain reaction. Consistent with their role in detoxification, the aryl hydrocarbon receptor (Ahr) (r(2) = 0.18, p = 0.045 for 2008 and 2009) and cytochrome P450 1A1 (Cyp1a1) (r(2) = 0.20, p < 0.001 for 2008 and 2009; r(2) = 0.43, p = 0.049 for 2010) transcripts were positively correlated with polychlorinated biphenyls (PCBs), the dominant POP in beluga. Principal Components Analysis distinguished between these two toxicology genes and 11 other genes primarily involved in growth, metabolism, and development. Factor 1 explained 56% of gene profiles, with these latter 11 gene transcripts displaying greater abundance in years coinciding with periods of low sea ice extent (2008 and 2010). δ(13)C results suggested a shift in feeding ecology and/or change in condition of these ice edge-associated beluga whales during these two years. While this provides insight into the legacy of PCBs in a remote environment, the possible impacts of a changing ice climate on the health of beluga underscores the need for long-term studies.

Groups of related belugas (Delphinapterus leucas) travel together during their seasonal migrations in and around Hudson Bay.

Social structure involving long-term associations with relatives should facilitate the learning of complex behaviours such as long-distance migration. In and around Hudson Bay (Canada), three stocks of beluga whales form a panmictic unit, but have different migratory behaviours associated with different summering areas. We analysed genetic variation at 13 microsatellite loci among 1524 belugas, to test hypotheses about social structure in belugas. We found significant proportions of mother-offspring pairs throughout the migratory cycle, but average relatedness extended beyond close kinship only during migration. Average relatedness was significantly above random expectations for pairs caught at the same site but on different days or months of a year, suggesting that belugas maintain associations with a network of relatives during migration. Pairs involving a female (female-female or male-female) were on average more related than pairs of males, and males seemed to disperse from their matrilineal group to associate with other mature males. Altogether, our results indicate that relatives other than strictly parents, and especially females, play a role in maintaining a social structure that could facilitate the learning of migration routes. Cultural conservatism may limit contributions from nearby summer stocks to endangered stocks such as the Eastern Hudson Bay beluga.

FLOCK provides reliable solutions to the "number of populations" problem.

Identifying groups of individuals forming coherent genetic clusters is relevant to many fields of biology. This paper addresses the K-partition problem: given a collection of genotypes, partition those genotypes into K groups, each group being a sample of the K source populations that are represented in the collection of genotypes. This problem involves allocating genotypes to genetic groups while building those groups at the same time without the use of any other a priori information. FLOCK is a non-Markov chain Monte Carlo (MCMC) algorithm that uses an iterative method to partition a collection of genotypes into k groups. Rules to estimate K are formulated and their validity firmly established by running simulations under several migration rates, migration regimes, number of loci, and values of K. FLOCK tended to build clusters largely consistent with the source samples. The performance of FLOCK was also compared with that of STRUCTURE and BAPS. FLOCK provided more accurate allocations to clusters and more reliable estimates of K; it also ran much faster than STRUCTURE. FLOCK is based on an entirely novel approach and provides a true alternative to the existing, MCMC based, algorithms. FLOCK v.2.0 for microsatellites or for AFLP markers can be downloaded from http://www.bio.ulaval.ca/no_cache/departement/professeurs/fiche_des_professeurs/professeur/11/13/.

[Molecular genetic study of the beluga (Delphinapterus leucas: Cetacea, Monodontidae) summering in the southern Sea of Okhotsk as compared to North American populations].

The structure of the left part of the mtDNA control region has been studied in Delphinapterus leucas captured in summer in the southern Sea of Okhotsk. The data have been compared with published results of similar studies on D. leucas at different sites along the American coast. A high specificity of maternal lineages of D. leucas from the southern Sea of Ok.