A new multiplexed magnetic capture-Droplet digital PCR tool for monitoring wildlife population health and pathogen surveillance.

Anthropogenic stressors are exacerbating the emergence and spread of pathogens worldwide. In regions like the Arctic, where ecosystems are particularly susceptible, marked changes are predicted in regional diversity, intensity, and patterns of infectious diseases. To understand such rapidly changing host-pathogen dynamics and mitigate the impacts of novel pathogens, we need sensitive disease surveillance tools. We developed and validated a novel multiplexed, magnetic capture, and ddPCR tool for the surveillance of multiple pathogens in polar bears, a sentinel species that is considered susceptible to climate change and other stressors with a pan-Arctic distribution. Through sequence-specific magnetic capture, we concentrated five target template sequences from three zoonotic bacteria (Erysipelothrix rhusiopathiae, Francisella tularensis, and Mycobacterium tuberculosis complex) and two parasitic (Toxoplasma gondii and Trichinella spp.) pathogens from large quantities (<100 g) of host tissue. We then designed and validated two multiplexed probe-based ddPCR assays for the amplification and detection of the low-concentration target DNA. Validations used 48 polar bear tissues (muscle and liver). We detected 14, 1, 3, 4, and 22 tissue positives for E. rhusiopathiae, F. tularensis, M. tuberculosis complex, T. gondii, and Trichinella spp., respectively. These multiplexed assays offer a rapid, specific tool for quantifying and monitoring the changing geographical and host distributions of pathogens relevant to human and animal health.

No evidence of inbreeding depression despite a historical severe bottleneck in the endangered Bermuda petrel (Pterodroma cahow).

The Bermuda petrel Pterodroma cahow is an island endemic seabird that belongs to the Procellariiformes, one of the most endangered orders of birds. Historical records suggest a significant population size decline following human settlement in Bermuda, bringing the species to near extinction. Since the 1950s, the population has been recovering aided by the implementation of an ongoing conservation plan. However, it still faces several threats, and negative genetic effects resulting from that drastic decline are to be expected, including inbreeding and genetic drift. We studied genetic diversity and levels of inbreeding, and their effects on individual fitness and mating choice. We also tested for a genetic signature of the recent demographic bottleneck. For this, we analyzed variation in thousands of nuclear single-nucleotide polymorphisms derived from double digest restriction site-associated DNA sequencing and 1 mitochondrial gene (cytochrome oxidase I). The results revealed that the Bermuda petrel suffered a recent genetic bottleneck and shows low mitochondrial diversity compared with other petrel species. Conversely, nuclear diversity was similar to that of other endangered petrels. Inbreeding levels were not high overall, although some individuals were highly inbred. However, we found no evidence that individual inbreeding or relatedness between mates affected hatching success, or that mate choice is influenced by kinship in this very small population.

Green credit guideline and enterprise export green-sophistication.

Green credit is a major policy innovation to guide enterprises to participate in environmental governance actively. This study uses the data of Chinese A-share listed companies from 2007 to 2016, takes the green credit guideline (GCG) issued in 2012 as a quasi-natural experiment, and uses a difference in difference (DID) model to test the effect of GCG on the enterprises' export green-sophistication (EGS) and its internal and external mechanisms. The study finds that GCG improves enterprises' EGS and research and development (R&D) investment is the intermediation channel for GCG to affect EGS. Results of heterogeneity analysis show that the role of GCG in promoting EGS is significantly reflected in enterprises that the government does not subsidize, enterprises in areas with a low degree of financial marketization development, state-owned enterprises, and enterprises with a high degree of equity incentive.

Genotyping-in-thousands by sequencing (GT-seq) of noninvasive faecal and degraded samples: A new panel to enable ongoing monitoring of Canadian polar bear populations.

Genetic monitoring using noninvasive samples provides a complement or alternative to traditional population monitoring methods. However, next-generation sequencing approaches to monitoring typically require high quality DNA and the use of noninvasive samples (e.g., scat) is often challenged by poor DNA quality and contamination by nontarget species. One promising solution is a highly multiplexed sequencing approach called genotyping-in-thousands by sequencing (GT-seq), which can enable cost-efficient genomics-based monitoring for populations based on noninvasively collected samples. Here, we develop and validate a GT-seq panel of 324 single nucleotide polymorphisms (SNPs) optimized for genotyping of polar bears based on DNA from noninvasively collected faecal samples. We demonstrate (1) successful GT-seq genotyping of DNA from a range of sample sources, including successful genotyping (>50% loci) of 62.9% of noninvasively collected faecal samples determined to contain polar bear DNA; and (2) that we can reliably differentiate individuals, ascertain sex, assess relatedness, and resolve population structure of Canadian polar bear subpopulations based on a GT-seq panel of 324 SNPs. Our GT-seq data reveal spatial-genetic patterns similar to previous polar bear studies but at lesser cost per sample and through use of noninvasively collected samples, indicating the potential of this approach for population monitoring. This GT-seq panel provides the foundation for a noninvasive toolkit for polar bear monitoring and can contribute to community-based programmes - a framework which may serve as a model for wildlife conservation and management for species worldwide.

Phylogenomics reveals viral sources, transmission, and potential superinfection in early-stage COVID-19 patients in Ontario, Canada.

The emergence and rapid global spread of SARS-CoV-2 demonstrates the importance of infectious disease surveillance, particularly during the early stages. Viral genomes can provide key insights into transmission chains and pathogenicity. Nasopharyngeal swabs were obtained from thirty-two of the first SARS-CoV-2 positive cases (March 18-30) in Kingston Ontario, Canada. Viral genomes were sequenced using Ion Torrent (n = 24) and MinION (n = 27) sequencing platforms. SARS-CoV-2 genomes carried forty-six polymorphic sites including two missense and three synonymous variants in the spike protein gene. The D614G point mutation was the predominate viral strain in our cohort (92.6%). A heterozygous variant (C9994A) was detected by both sequencing platforms but filtered by the ARTIC network bioinformatic pipeline suggesting that heterozygous variants may be underreported in the SARS-CoV-2 literature. Phylogenetic analysis with 87,738 genomes in the GISAID database identified global origins and transmission events including multiple, international introductions as well as community spread. Reported travel history validated viral introduction and transmission inferred by phylogenetic analysis. Molecular epidemiology and evolutionary phylogenetics may complement contact tracing and help reconstruct transmission chains of emerging diseases. Earlier detection and screening in this way could improve the effectiveness of regional public health interventions to limit future pandemics.

A real-time PCR assay to accurately quantify polar bear DNA in fecal extracts.

DNA extracted from fecal samples contains DNA from the focal species, food, bacteria and pathogens. Most DNA quantification methods measure total DNA and cannot differentiate among sources. Despite the desirability of noninvasive fecal sampling for studying wildlife populations, low amounts of focal species DNA make it difficult to use for next-generation sequencing (NGS), where accurate DNA quantification is critical for normalization. Two factors are required prior to using fecal samples in NGS libraries: (1) an accurate quantification method for the amount of target DNA and (2) a determination of the relative amount of target DNA needed for successful single nucleotide polymorphism genotyping assays. Here, we address these needs by developing primers to amplify a 101 bp region of the nuclear F2 gene and a quantitative PCR (qPCR) assay that allows the accurate quantification of the amount of polar bear (Ursus maritimus) DNA in fecal extracts. We test the assay on pure polar bear DNA extracted from muscle tissue and find a high correlation between fluorometric and qPCR quantifications. The qPCR assay was also successfully used to quantify the amount of DNA derived from polar bears in fecal extractions. Orthologs of the F2 gene have been identified across vertebrates; thus, similar qPCR assays could be developed for other species to enable noninvasive studies.

Isolation and characterization of microsatellite loci for storm-petrels.

Primers were developed for 10 microsatellite loci for two species of Oceanodroma storm-petrels. Variability was tested in 27 O. castro and 22 O. monteiroi from the Azores, and 24 O. leucorhoa from Norway. At least six loci amplified reliably and were polymorphic in each species. The number of alleles per locus averaged 4.6, and observed heterozygosities averaged 0.41. Most primers also yielded polymerase chain reaction products in O. tethys, O. hornbyi and Pterodroma phaeopygia. These loci are being used to assay population genetic structure in storm-petrels.

Isolation and characterization of microsatellite loci for eastern foxsnakes (Elaphe gloydi).

New statistical techniques and software have improved our ability to quantify fine-scale population structure. We isolated 11 microsatellite markers for eastern foxsnakes (Elaphe gloydi) and surveyed the variability of each using 115-136 individuals from one population in southwestern Ontario. We determined that all loci were variable and conformed to Hardy-Weinberg expectations and that no locus demonstrated linkage disequilibrium. We are using these markers to quantify the scale of gene flow, and to assess the effect of variable landscapes on connectivity and genetic diversity of three disjunct regional populations in the northern portion of the species' range.