High genetic load without purging in caribou, a diverse species at risk.

High intra-specific genetic diversity is associated with adaptive potential, which is key for resilience to global change. However, high variation may also support deleterious alleles through genetic load, thereby increasing the risk of inbreeding depression if population sizes decrease. Purging of deleterious variation has been demonstrated in some threatened species. However, less is known about the costs of declines and inbreeding in species with large population sizes and high genetic diversity even though this encompasses many species globally that are expected to undergo population declines. Caribou is a species of ecological and cultural significance in North America with a wide distribution supporting extensive phenotypic variation but with some populations undergoing significant declines resulting in their at-risk status in Canada. We assessed intra-specific genetic variation, adaptive divergence, inbreeding, and genetic load across populations with different demographic histories using an annotated chromosome-scale reference genome and 66 whole-genome sequences. We found high genetic diversity and nine phylogenomic lineages across the continent with adaptive diversification of genes, but also high genetic load among lineages. We found highly divergent levels of inbreeding across individuals, including the loss of alleles by drift but not increased purging in inbred individuals, which had more homozygous deleterious alleles. We also found comparable frequencies of homozygous deleterious alleles between lineages regardless of nucleotide diversity. Thus, further inbreeding may need to be mitigated through conservation efforts. Our results highlight the "double-edged sword" of genetic diversity that may be representative of other species atrisk affected by anthropogenic activities.

Circadian gene variation in relation to breeding season and latitude in allochronic populations of two pelagic seabird species complexes.

Annual cues in the environment result in physiological changes that allow organisms to time reproduction during periods of optimal resource availability. Understanding how circadian rhythm genes sense these environmental cues and stimulate the appropriate physiological changes in response is important for determining the adaptability of species, especially in the advent of changing climate. A first step involves characterizing the environmental correlates of natural variation in these genes. Band-rumped and Leach's storm-petrels (Hydrobates spp.) are pelagic seabirds that breed across a wide range of latitudes. Importantly, some populations have undergone allochronic divergence, in which sympatric populations use the same breeding sites at different times of year. We investigated the relationship between variation in key functional regions of four genes that play an integral role in the cellular clock mechanism-Clock, Bmal1, Cry2 and Per2-with both breeding season and absolute latitude in these two species complexes. We discovered that allele frequencies in two genes, Clock and Bmal1, differed between seasonal populations in one archipelago, and also correlated with absolute latitude of breeding colonies. These results indicate that variation in these circadian rhythm genes may be involved in allochronic speciation, as well as adaptation to photoperiod at breeding locations.

Genomic population structure and inbreeding history of Lake Superior caribou.

Caribou (Rangifer tarandus) have experienced dramatic declines in both range and population size across Canada over the past century. Boreal caribou (R. t. caribou), 1 of the 12 Designatable Units, has lost approximately half of its historic range in the last 150 years, particularly along the southern edge of its distribution. Despite this overall northward contraction, some populations have persisted at the trailing range edge, over 150 km south of the continuous boreal caribou range in Ontario, along the coast and nearshore islands of Lake Superior. The population history of caribou along Lake Superior remains unclear. It appears that these caribou likely represent a remnant distribution at the trailing edge of the receding population of boreal caribou, but they may also exhibit local adaptation to the coastal environment. A better understanding of the population structure and history of caribou along Lake Superior is important for their conservation and management. Here, we use high-coverage whole genomes (N = 20) from boreal, eastern migratory, and barren-ground caribou sampled in Manitoba, Ontario, and Quebec to investigate population structure and inbreeding histories. We discovered that caribou from the Lake Superior range form a distinct group but also found some evidence of gene flow with the continuous boreal caribou range. Notably, caribou along Lake Superior demonstrated relatively high levels of inbreeding (measured as runs of homozygosity; ROH) and genetic drift, which may contribute to the differentiation observed between ranges. Despite inbreeding, caribou along Lake Superior retained high heterozygosity, particularly in genomic regions without ROH. These results suggest that they present distinct genomic characteristics but also some level of gene flow with the continuous range. Our study provides key insights into the genomics of the southernmost range of caribou in Ontario, beginning to unravel the evolutionary history of these small, isolated caribou populations.

Using machine learning algorithms to enhance the diagnostic performance of electrical impedance myography.

INTRODUCTION/AIMS: We assessed the classification performance of machine learning (ML) using multifrequency electrical impedance myography (EIM) values to improve upon diagnostic outcomes as compared to those based on a single EIM value. METHODS: EIM data was obtained from unilateral excised gastrocnemius in eighty diseased mice (26 D2-mdx, Duchenne muscular dystrophy model, 39 SOD1G93A ALS model, and 15 db/db, a model of obesity-induced muscle atrophy) and 33 wild-type (WT) animals. We assessed the classification performance of a ML random forest algorithm incorporating all the data (multifrequency resistance, reactance and phase values) comparing it to the 50 kHz phase value alone. RESULTS: ML outperformed the 50 kHz analysis as based on receiver-operating characteristic curves and measurement of the area under the curve (AUC). For example, comparing all diseases together versus WT from the test set outputs, the AUC was 0.52 for 50 kHz phase, but was 0.94 for the ML model. Similarly, when comparing ALS versus WT, the AUCs were 0.79 for 50 kHz phase and 0.99 for ML. DISCUSSION: Multifrequency EIM using ML improves upon classification compared to that achieved with a single-frequency value. ML approaches should be considered in all future basic and clinical diagnostic applications of EIM.

A de novo genome assembly and annotation of the southern flying squirrel (Glaucomys volans).

Northern (Glaucomys sabrinus) and southern (Glaucomys volans) flying squirrels are widespread species distributed across North America. Northern flying squirrels are common inhabitants of the boreal forest, also occurring in coniferous forest remnants farther south, whereas the southern flying squirrel range is centered in eastern temperate woodlands. These two flying squirrel species exhibit a hybrid zone across a latitudinal gradient in an area of recent secondary contact. Glaucomys hybrid offspring are viable and can successfully backcross with either parental species, however, the fitness implications of such events are currently unknown. Some populations of G. sabrinus are endangered, and thus, interspecific hybridization is a key conservation concern in flying squirrels. To provide a resource for future studies to evaluate hybridization and possible introgression, we sequenced and assembled a de novo long-read genome from a G. volans individual sampled in southern Ontario, Canada, while four short-read genomes (two G. sabrinus and two G. volans, all from Ontario) were resequenced on Illumina platforms. The final genome assembly consisted of approximately 2.40 Gb with a scaffold N50 of 455.26 Kb. Benchmarking Universal Single-Copy Orthologs reconstructed 3,742 (91.2%) complete mammalian genes and genome annotation using RNA-Seq identified the locations of 19,124 protein-coding genes. The four short-read individuals were aligned to our reference genome to investigate the demographic history of the two species. A principal component analysis clearly separated resequenced individuals, while inferring population size history using the Pairwise Sequentially Markovian Coalescent model noted an approximate species split 1 million years ago, and a single, possibly recently introgressed individual.

Relationships between in vivo surface and ex vivo electrical impedance myography measurements in three different neuromuscular disorder mouse models.

Electrical impedance myography (EIM) using surface techniques has shown promise as a means of diagnosing and tracking disorders affecting muscle and assessing treatment efficacy. However, the relationship between such surface-obtained impedance values and pure muscle impedance values has not been established. Here we studied three groups of diseased and wild-type (WT) animals, including a Duchenne muscular dystrophy model (the D2-mdx mouse), an amyotrophic lateral sclerosis (ALS) model (the SOD1 G93A mouse), and a model of fat-related atrophy (the db/db diabetic obese mouse), performing hind limb measurements using a standard surface array and ex vivo measurements on freshly excised gastrocnemius muscle. A total of 101 animals (23 D2-mdx, 43 ALS mice, 12 db/db mice, and corresponding 30 WT mice) were studied with EIM across a frequency range of 8 kHz to 1 MHz. For both D2-mdx and ALS models, moderate strength correlations (Spearman rho values generally ranging from 0.3-0.7, depending on the impedance parameter (i.e., resistance, reactance and phase) were obtained. In these groups of animals, there was an offset in frequency with impedance values obtained at higher surface frequencies correlating more strongly to impedance values obtained at lower ex vivo frequencies. For the db/db model, correlations were comparatively weaker and strongest at very high and very low frequencies. When combining impedance data from all three disease models together, moderate correlations persisted (with maximal Spearman rho values of 0.45). These data support that surface EIM data reflect ex vivo muscle tissue EIM values to a moderate degree across several different diseases, with the highest correlations occurring in the 10-200 kHz frequency range. Understanding these relationships will prove useful for future applications of the technique of EIM in the assessment of neuromuscular disorders.

Population dynamics of caribou shaped by glacial cycles before the last glacial maximum.

Pleistocene glacial cycles influenced the diversification of high-latitude wildlife species through recurrent periods of range contraction, isolation, divergence, and expansion from refugia and subsequent admixture of refugial populations. We investigate population size changes and the introgressive history of caribou (Rangifer tarandus) in western Canada using 33 whole genome sequences coupled with larger-scale mitochondrial data. We found that a major population expansion of caribou occurred starting around 110,000 years ago (kya), the start of the last glacial period. Additionally, we found effective population sizes of some caribou reaching ~700,000 to 1,000,000 individuals, one of the highest recorded historical effective population sizes for any mammal species thus far. Mitochondrial analyses dated introgression events prior to the LGM dating to 20-30 kya and even more ancient at 60 kya, coinciding with colder periods with extensive ice coverage, further demonstrating the importance of glacial cycles and events prior to the LGM in shaping demographic history. Reconstructing the origins and differential introgressive history has implications for predictions on species responses under climate change. Our results have implications for other whole genome analyses using pairwise sequentially Markovian coalescent (PSMC) analyses, as well as highlighting the need to investigate pre-LGM demographic patterns to fully reconstruct the origin of species diversity, especially for high-latitude species.

Cytonuclear discordance in the crowned-sparrows, Zonotrichia atricapilla and Zonotrichia leucophrys.

The golden-crowned (Zonotrichia atricapilla) and white-crowned (Z. leucophrys) sparrows have been presented as a compelling case for rapid speciation. They display divergence in song and plumage with overlap in their breeding ranges implying reproductive isolation, but have almost identical mitochondrial genomes. Previous research proposed hybridization and subsequent mitochondrial introgression as an alternate explanation, but lacked robust nuclear gene trees to distinguish between introgression and incomplete lineage sorting. We test for signatures of these processes between Z. atricapilla and Z. leucophrys, and investigate the relationships among Z. leucophrys subspecies, using mitochondrial sequencing and a reduced representation nuclear genomic dataset. Contrary to the paraphyly evident in mitochondrial gene trees, we confirmed the reciprocal monophyly of Z. atricapilla and Z. leucophrys using large panels of single nucleotide polymorphisms (SNPs). The pattern of cytonuclear discordance is consistent with limited, historical hybridization and mitochondrial introgression, rather than a recent origin and incomplete lineage sorting between recent sister species. We found evidence of nuclear phylogeographic structure within Z. leucophrys with two distinct clades. Altogether, our results indicate deeper divergences between Z. atricapilla and Z. leucophrys than inferred using mitochondrial markers. Our results demonstrate the limitations of relying solely on mitochondrial DNA for taxonomy, and raise questions about the possibility of selection on the mitochondrial genome during temperature oscillations (e.g. during the Pleistocene). Historical mitochondrial introgression facilitated by past environmental changes could cause erroneous dating of lineage splitting in other taxa when based on mitochondrial DNA alone.

Estimating myofiber cross-sectional area and connective tissue deposition with electrical impedance myography: A study in D2-mdx mice.

INTRODUCTION: Surface electrical impedance myography (sEIM) has the potential for providing information on muscle composition and structure noninvasively. We sought to evaluate its use to predict myofiber size and connective tissue deposition in the D2-mdx model of Duchenne muscular dystrophy (DMD). METHODS: We applied a prediction algorithm, the least absolute shrinkage and selection operator, to select specific EIM measurements obtained with surface and ex vivo EIM data from D2-mdx and wild-type (WT) mice (analyzed together or separately). We assessed myofiber cross-sectional area histologically and hydroxyproline (HP), a surrogate measure for connective tissue content, biochemically. RESULTS: Using WT and D2-mdx impedance values together in the algorithm, sEIM gave average root-mean-square errors (RMSEs) of 26.6% for CSA and 45.8% for HP, which translate into mean errors of ±363 µm2 for a mean CSA of 1365 µm2 and of ±1.44 µg HP/mg muscle for a mean HP content of 3.15 µg HP/mg muscle. Stronger predictions were obtained by analyzing sEIM data from D2-mdx animals alone (RMSEs of 15.3% for CSA and 34.1% for HP content). Predictions made using ex vivo EIM data from D2-mdx animals alone were nearly equivalent to those obtained with sEIM data (RMSE of 16.59% for CSA), and slightly more accurate for HP (RMSE of 26.7%). DISCUSSION: Surface EIM combined with a predictive algorithm can provide estimates of muscle pathology comparable to values obtained using ex vivo EIM, and can be used as a surrogate measure of disease severity and progression and response to therapy.

Predicting myofiber cross-sectional area and triglyceride content with electrical impedance myography: A study in db/db mice.

BACKGROUND: Electrical impedance myography (EIM) provides insight into muscle composition and structure. We sought to evaluate its use in a mouse obesity model characterized by myofiber atrophy. METHODS: We applied a prediction algorithm, ie, the least absolute shrinkage and selection operator (LASSO), to surface, needle array, and ex vivo EIM data from db/db and wild-type mice and assessed myofiber cross-sectional area (CSA) histologically and triglyceride (TG) content biochemically. RESULTS: EIM data from all three modalities provided acceptable predictions of myofiber CSA with average root mean square error (RMSE) of 15% in CSA (ie, ±209 µm2 for a mean CSA of 1439 µm2 ) and TG content with RMSE of 30% in TG content (ie, ±7.3 nmol TG/mg muscle for a mean TG content of 25.4 nmol TG/mg muscle). CONCLUSIONS: EIM combined with a predictive algorithm provides reasonable estimates of myofiber CSA and TG content without the need for biopsy.

The role of introgression and ecotypic parallelism in delineating intraspecific conservation units.

Parallel evolution can occur through selection on novel mutations, standing genetic variation or adaptive introgression. Uncovering parallelism and introgressed populations can complicate management of threatened species as parallelism may have influenced conservation unit designations and admixed populations are not generally considered under legislations. We examined high coverage whole-genome sequences of 30 caribou (Rangifer tarandus) from across North America and Greenland, representing divergent intraspecific lineages, to investigate parallelism and levels of introgression contributing to the formation of ecotypes. Caribou are split into four subspecies and 11 extant conservation units, known as designatable units (DUs), in Canada. Using genomes from all four subspecies and six DUs, we undertake demographic reconstruction and confirm two previously inferred instances of parallel evolution in the woodland subspecies and uncover an additional instance of parallelism of the eastern migratory ecotype. Detailed investigations reveal introgression in the woodland subspecies, with introgressed regions found spread throughout the genomes encompassing both neutral and functional sites. Our investigations using whole genomes highlight the difficulties in unequivocally demonstrating parallelism through adaptive introgression in nonmodel species with complex demographic histories, with standing variation and introgression both potentially involved. Additionally, the impact of parallelism and introgression on conservation policy for management units needs to be considered in general, and the caribou designations will need amending in light of our results. Uncovering and decoupling parallelism and differential patterns of introgression will become prevalent with the availability of comprehensive genomic data from nonmodel species, and we highlight the need to incorporate this into conservation unit designations.

The Caribou (Rangifer tarandus) Genome.

Rangifer tarandus, known as caribou or reindeer, is a widespread circumpolar species which presents significant variability in their morphology, ecology, and genetics. A genome was sequenced from a male boreal caribou (R. t. caribou) from Manitoba, Canada. Both paired end and Chicago libraries were constructed and sequenced on Illumina platforms. The final assembly consists of approximately 2.205 Gb, and has a scaffold N50 of 11.765 Mb. BUSCO (Benchmarking Universal Single-Copy Orthologs) reconstructed 3820 (93.1%) complete mammalian genes, and genome annotation identified the locations of 33,177 protein-coding genes. An alignment to the bovine genome was carried out, indicating sequence coverage on all bovine chromosomes. A high-quality reference genome will be invaluable for evolutionary research and for conservation efforts for the species. Further information about the genome, including a FASTA file of the assembly and the annotation files, is available on our caribou genome website. Raw sequence data is available at the National Centre for Biotechnology Information (NCBI), under the BioProject accession number PRJNA549927.

Cryptic species and independent origins of allochronic populations within a seabird species complex (Hydrobates spp.).

Humans are inherently biased towards naming species based on morphological differences, which can lead to reproductively isolated species being mistakenly classified as one if they are morphologically similar. Recognising cryptic diversity is needed to understand drivers of speciation fully, and for accurate estimates of global biodiversity and assessments for conservation. We investigated cryptic species across the range of band-rumped storm-petrels (Hydrobates spp.): highly pelagic, nocturnal seabirds that breed on tropical and sub-tropical islands in the Atlantic and Pacific Oceans. In many breeding colonies, band-rumped storm-petrels have sympatric but temporally isolated (allochronic) populations; we sampled all breeding locations and allochronic populations. Using mitochondrial control region sequences from 754 birds, cytochrome b sequences from 69 birds, and reduced representation sequencing of the nuclear genomes of 133 birds, we uncovered high levels of genetic structuring. Population genomic analyses revealed up to seven unique clusters, and phylogenomic reconstruction showed that these represent seven monophyletic groups. We uncovered up to six independent breeding season switches across the phylogeny, spanning the continuum from genetically undifferentiated temporal populations to full allochronic species. Thus, band-rumped storm-petrels encompass multiple cryptic species, with non-geographic barriers potentially comprising strong barriers to gene flow.

Estimating Myofiber Size With Electrical Impedance Myography: a Study In Amyotrophic Lateral Sclerosis MICE.

INTRODUCTION: A method for quantifying myofiber size noninvasively would find wide use, including primary diagnosis and evaluating response to therapy. METHODS: Using prediction algorithms, including the least absolute shrinkage and selection operator, we applied multifrequency electrical impedance myography (EIM) to amyotrophic lateral sclerosis superoxide dismutase 1 G93A mice of different ages and assessed myofiber size histologically. RESULTS: Multifrequency EIM data provided highly accurate predictions of myofiber size, with a root mean squared error (RMSE) of only 14% in mean myofiber area (corresponding to ± 207 µm2 for a mean area of 1,488 µm2 ) and an RMSE of only 8.8% in predicting the coefficient of variation in fiber size distribution. DISCUSSION: This impedance-based approach provides predictive variables to assess myofiber size and distribution with good accuracy, particularly in diseases in which myofiber atrophy is the predominant histological feature, without the requirement for biopsy or burdensome quantification. Muscle Nerve 58: 713-717, 2018.

Use of carbon-based composites to enhance performance of TiO2 for the simultaneous removal of nitrates and organics from aqueous environments.

The simultaneous photocatalytic removal of nitrate from aqueous environment in presence of organic hole scavenger using TiO2 has long been explored. However, the use of unmodified TiO2 in such reaction resulted in non-performance or release of significant amount of undesirable reaction products in the process, a problem that triggered surface modification of TiO2 for enhanced photocatalytic performance. Previous studies focused on decreasing rate of charge carrier recombination and absorption of light in the visible region. Yet, increasing active sites and adsorption capacity by combining TiO2 with a high surface area adsorbent such as activated carbon (AC) remains unexploited. This study reports the potential of such modification in simultaneous removal of nitrates and oxalic acid in aqueous environment. The adsorptive behaviour of nitrate and oxalic acid on TiO2 and TiO2/AC composites were studied. The Langmuir adsorption coefficient for nitrate was four times greater than that of oxalic acid. However, the amount of oxalic acid adsorbed was about 10 times greater than the amount of nitrate taken up. Despite this advantage, the materials did not appear to produce more active photocatalysts for the simultaneous degradation of nitrate and oxalic acid. The photocatalytic activity of TiO2 and its carbon-based composites was improved by combination with Cu2O particles. Consequently, 2.5 Cu2O/TiO2 exhibited the maximum photocatalytic performance with 57.6 and 99.8% removal of nitrate and oxalic acid, respectively, while selectivity stood at 45.7, 12.4 and 41.9% for NH4+, NO2- and N2, respectively. For the carbon based, 2.5 Cu2O/TiO2-20AC showed removal of 12.7% nitrate and 80.3% oxalic acid and achieved 21.6, 0 and 78.4% selectivity for NH4+, NO2- and N2, respectively. Using the optimal AC loading (20 wt%) resulted in significant decrease in the selectivity for NH4+ with no formation of NO2-, which unveils that selectivity for N2 and low/no selectivity for undesirable products can be manipulated by controlling the rate of consumption of oxalic acid. In contract, no nitrate reduction was observed with Cu2O promoted TiO2-T and its TiO2-(T)-20AC, which may be connected to amorphous nature of TiO2-T and perhaps served as charge carrier trapping sites that impeded activity.

Electrical impedance myography as a biomarker of myostatin inhibition with ActRIIB-mFc: a study in wild-type mice.

AIM: We sought to determine the sensitivity of electrical impedance myography (EIM) to myofiber hypertrophy induced by treatment with various doses of ActRIIB-mFc, an inhibitor of myostatin signaling. METHODS: Wild-type C57BL/6 J mice (n = 40, male) were treated with three different doses of ActRIIB-mFc (i.e., RAP-031) or vehicle twice weekly for 5 weeks. End point assessments included gastrocnemius EIM, force measurements, muscle mass and myofiber size quantification. RESULTS: ActRIIB-mFc increased body mass, muscle mass and myofiber size across all doses. Alterations in EIM 50 kHz phase and center frequency (fc) were also present, with trends in a dose-dependent fashion. Significant correlations between EIM parameters and myofiber/functional data were identified. CONCLUSION: EIM outcomes can serve as effective biomarkers of myostatin signaling inhibition, demonstrating a dose sensitivity and correlation to standard assessments.

Predicting myofiber size with electrical impedance myography: A study in immature mice.

INTRODUCTION: Electrical impedance can be used to estimate cellular characteristics. We sought to determine whether it could be used to approximate myofiber size using standard prediction modeling approaches. METHODS: Forty-four C57BL/6J wild-type immature mice of varying ages underwent electrical impedance myography (EIM) with a needle electrode array placed in the gastrocnemius. Animals were then humanely killed and muscle fixed, stained, and myofiber size quantified. Two different statistical prediction models were then applied. RESULTS: Impedance parameters showed major variation with increasing myofiber size. The prediction models based on EIM data alone were able to predict fiber size, with errors in the range of ±69.05-78.44 µm2 (16.19%-18.40% with respect to the average myofiber size). DISCUSSION: By using well-established statistical models, EIM data alone can provide a satisfactory estimate of myofiber size. Additional study of this approach for approximating myofiber size without the requirement of removing tissue for histological analysis is warranted. Muscle Nerve, 2018.