Molecular ecology of the sleeper shark subgenus Somniosus (Somniosus) reveals genetic homogeneity within species and lack of support for S. antarcticus.

Inferences made from molecular data support regional stock assessment goals by providing insights into the genetic population dynamics of enigmatic species. Population genomics metrics, such as genetic diversity and population connectivity, serve as useful proxies for species health and stability. Sleeper sharks (genus Somniosus) are ecologically important deep-sea predators, estimated to reach ages of 250 to 300 yr and taking decades to reach sexual maturity. The subgenus Somniosus (Somniosus) is comprised of 3 species: S. pacificus, S. microcephalus, and S. antarcticus. Given the life history strategy of somniosids, they are vulnerable to overfishing and population declines. Further, data to assess the stocks of these species are limited. To address this deficiency, we used the reduced representation library method Restriction-site Associated DNA sequencing (RADseq) to conduct phylogenomic and population genomics analyses, providing novel information for use in stock assessments. Our results strongly support the species status of S. microcephalus (N = 79), but recover S. antarcticus (N = 2) intermixed within the S. pacificus (N = 170) clade. Population genomics analyses reveal genetic homogeneity within S. pacificus and S. microcephalus, and estimates of effective population size were in the hundreds for both species. Kinship analysis identified 2 first-degree relative pairs within our dataset (1 within each species). Our results contribute new information for stock assessments of these uniquely long-lived species by providing the strongest molecular evidence to date for the synonymization of S. antarcticus and S. pacificus, as well as estimating population genomic metrics for each supported species within the Somniosus (Somniosus) subgenus.

Reduced genetic diversity associated with the northern expansion of an amphibian species with high habitat specialization, Ascaphus truei, resolved using two types of genetic markers.

Reconstruction of historical relationships between geographic regions within a species' range can indicate dispersal patterns and help predict future responses to shifts in climate. Ascaphus truei (coastal tailed frog) is an indicator species of the health of forests and perennial streams in the Coastal and Cascade Mountains of the Pacific Northwest of North America. We used two genetic techniques-microsatellite and genotype-by-sequencing (GBS)-to compare the within-region genetic diversity of populations near the northern extent of the species' range (British Columbia, Canada) to two geographic regions in British Columbia and two in Washington, USA, moving toward the core of the range. Allelic richness and heterozygosity declined substantially as latitude increased. The northernmost region had the lowest mean expected heterozygosities for both techniques (microsatellite, M = 0.20, SE = 0.080; GBS, M = 0.025, SE = 0.0010) and the southernmost region had the highest (microsatellite, M = 0.88, SE = 0.054; GBS, M = 0.20, SE = 0.0029). The northernmost regions (NC and MC) clustered together in population structure models for both genetic techniques. Our discovery of reduced diversity may have important conservation and management implications for population connectivity and the response of A. truei to climate change.

Autumn shifts in cold tolerance metabolites in overwintering adult mountain pine beetles.

The mountain pine beetle, Dendroctonus ponderosae (Coleoptera: Curculionidae) is a major forest pest of pines in western North America. Beetles typically undergo a one-year life cycle with larval cold hardening in preparation for overwintering. Two-year life cycle beetles have been observed but not closely studied. This study tracks cold-hardening and preparation for overwintering by adult mountain pine beetles in their natal galleries. Adults were collected in situ between September and December 2016 for a total of nine time points during 91 days. Concentrations of 41 metabolites in these pooled samples were assessed using quantitative nuclear magnetic resonance (NMR). Levels of glycerol and proline increased significantly with lowering temperature during the autumn. Newly eclosed mountain pine beetles appear to prepare for winter by generating the same cold-tolerance compounds found in other insect larvae including mountain pine beetle, but high on-site mortality suggested that two-year life cycle adults have a less efficacious acclimation process. This is the first documentation of cold acclimation metabolite production in overwintering new adult beetles and is evidence of physiological plasticity that would allow evolution by natural selection of alternate life cycles (shortened or lengthened) under a changing climate or during expansion into new geoclimatic areas.

Selection of the sex-linked inhibitor of apoptosis in mountain pine beetle (Dendroctonus ponderosae) driven by enhanced expression during early overwintering.

The mountain pine beetle (Dendroctonus ponderosae) is an insect native to western North America; however, its geographical range has recently expanded north in BC and east into Alberta. To understand the population structure in the areas of expansion, 16 gene-linked microsatellites were screened and compared to neutral microsatellites using outlier analyses of Fst and Fct values. One sex-linked gene, inhibitor of apoptosis (IAP), showed a strong signature of positive selection for neo-X alleles and was analyzed for evidence of adaptive variation. Alleles of IAP were sequenced, and differences between the neo-X and neo-Y alleles were consistent with neutral evolution suggesting that the neo-Y allele may not be under functional constraints. Neo-Y alleles were amplified from gDNA, but not effectively from cDNA, suggesting that there was little IAP expression from neo-Y alleles. There were no differences in overall IAP expression between males and females with the common northern neo-X allele suggesting that the neo-X allele in males compensates for the reduced expression of neo-Y alleles. However, males lacking the most common northern neo-X allele thought to be selected for in northern populations had reduced overall IAP expression in early October-at a time when beetles are preparing for overwintering. This suggests that the most common allele may have more rapid upregulation. The reduced function of neo-Y alleles of IAP suggested by both sequence differences and lower levels of expression may foster a highly selective environment for neo-X alleles such as the common northern allele with more efficient upregulation.

Adaptive and neutral markers both show continent-wide population structure of mountain pine beetle (Dendroctonus ponderosae).

Assessments of population genetic structure and demographic history have traditionally been based on neutral markers while explicitly excluding adaptive markers. In this study, we compared the utility of putatively adaptive and neutral single-nucleotide polymorphisms (SNPs) for inferring mountain pine beetle population structure across its geographic range. Both adaptive and neutral SNPs, and their combination, allowed range-wide structure to be distinguished and delimited a population that has recently undergone range expansion across northern British Columbia and Alberta. Using an equal number of both adaptive and neutral SNPs revealed that adaptive SNPs resulted in a stronger correlation between sampled populations and inferred clustering. Our results suggest that adaptive SNPs should not be excluded prior to analysis from neutral SNPs as a combination of both marker sets resulted in better resolution of genetic differentiation between populations than either marker set alone. These results demonstrate the utility of adaptive loci for resolving population genetic structure in a nonmodel organism.

How the mountain pine beetle (Dendroctonus ponderosae) breached the Canadian Rocky Mountains.

The mountain pine beetle (MPB; Dendroctonus ponderosae Hopkins), a major pine forest pest native to western North America, has extended its range north and eastward during an ongoing outbreak. Determining how the MPB has expanded its range to breach putative barriers, whether physical (nonforested prairie and high elevation of the Rocky Mountains) or climatic (extreme continental climate where temperatures can be below -40 °C), may contribute to our general understanding of range changes as well as management of the current epidemic. Here, we use a panel of 1,536 single nucleotide polymorphisms (SNPs) to assess population genetic structure, connectivity, and signals of selection within this MPB range expansion. Biallelic SNPs in MPB from southwestern Canada revealed higher genetic differentiation and lower genetic connectivity than in the northern part of its range. A total of 208 unique SNPs were identified using different outlier detection tests, of which 32 returned annotations for products with putative functions in cholesterol synthesis, actin filament contraction, and membrane transport. We suggest that MPB has been able to spread beyond its previous range by adjusting its cellular and metabolic functions, with genome scale differentiation enabling populations to better withstand cooler climates and facilitate longer dispersal distances. Our study is the first to assess landscape-wide selective adaptation in an insect. We have shown that interrogation of genomic resources can identify shifts in genetic diversity and putative adaptive signals in this forest pest species.

Characterization of class I- and class II-like major histocompatibility complex loci in pedigrees of North Atlantic right whales.

North Atlantic right whales have one of the lowest levels of genetic variation at minisatellite loci, microsatellite loci, and mitochondrial control region haplotypes among mammals. Here, adaptive variation at the peptide binding region of class I and class II DRB-like genes of the major histocompatibility complex was assessed. Amplification of a duplicated region in 222 individuals revealed at least 11 class II alleles. Six alleles were assigned to the locus Eugl-DRB1 and 5 alleles were assigned to the locus Eugl-DRB2 by assessing segregation patterns of alleles from 81 parent/offspring pedigrees. Pedigree analysis indicated that these alleles segregated into 12 distinct haplotypes. Genotyping a smaller subset of unrelated individuals (n = 5 and 10, respectively) using different primer sets revealed at least 2 class II pseudogenes (with ≥ 4 alleles) and at least 3 class I loci (with ≥ 6 alleles). Class II sequences were significantly different from neutrality at peptide binding sites suggesting loci may be under the influence of balancing selection. Trans-species sharing of alleles was apparent for class I and class II sequences. Characterization of class II loci represents the first step in determining the relationship between major histocompatibility complex variability and factors affecting health and reproduction in this species.

Permanent genetic resources added to molecular ecology resources database 1 April 2013-31 May 2013.

This article documents the addition of 234 microsatellite marker loci to the Molecular Ecology Resources Database. Loci were developed for the following species: Acipenser sinensis, Aleochara bilineata, Aleochara bipustulata, Barbus meridionalis, Colossoma macropomum, Delia radicum, Drosophila nigrosparsa, Fontainea picrosperma, Helianthemum cinereum, Liomys pictus, Megabalanus azoricus, Pelteobagrus vachelli, Pleuragramma antarcticum, Podarcis hispanica type 1A, Sardinella brasiliensis and Sclerotinia homoeocarpa. These loci were cross-tested on the following species: Acipenser dabryanus, Barbus balcanicus, Barbus barbus, Barbus cyclolepis, Drosophila hydei, Drosophila melanogaster, Drosophila obscura, Drosophila subobscura, Fontainea australis, Fontainea fugax, Fontainea oraria, Fontainea rostrata, Fontainea venosa, Podarcis bocagei, Podarcis carbonelli, Podarcis liolepis, Podarcis muralis and Podarcis vaucheri.

Spatial genetic structure of the mountain pine beetle (Dendroctonus ponderosae) outbreak in western Canada: historical patterns and contemporary dispersal.

Environmental change has a wide range of ecological consequences, including species extinction and range expansion. Many studies have shown that insect species respond rapidly to climatic change. A mountain pine beetle epidemic of record size in North America has led to unprecedented mortality of lodgepole pine, and a significant range expansion to the northeast of its historic range. Our goal was to determine the spatial genetic variation found among outbreak population from which genetic structure, and dispersal patterns may be inferred. Beetles from 49 sampling locations throughout the outbreak area in western Canada were analysed at 13 microsatellite loci. We found significant north-south population structure as evidenced by: (i) Bayesian-based analyses, (ii) north-south genetic relationships and diversity gradients; and (iii) a lack of isolation-by-distance in the northernmost cluster. The north-south structure is proposed to have arisen from the processes of postglacial colonization as well as recent climate-driven changes in population dynamics. Our data support the hypothesis of multiple sources of origin for the outbreak and point to the need for population specific information to improve our understanding and management of outbreaks. The recent range expansion across the Rocky Mountains into the jack/lodgepole hybrid and pure jack pine zones of northern Alberta is consistent with a northern British Columbia origin. We detected no loss of genetic variability in these populations, indicating that the evolutionary potential of mountain pine beetle to adapt has not been reduced by founder events. This study illustrates a rapid range-wide response to the removal of climatic constraints, and the potential for range expansion of a regional population.

Spatial genetic structure of a symbiotic beetle-fungal system: toward multi-taxa integrated landscape genetics.

Spatial patterns of genetic variation in interacting species can identify shared features that are important to gene flow and can elucidate co-evolutionary relationships. We assessed concordance in spatial genetic variation between the mountain pine beetle (Dendroctonus ponderosae) and one of its fungal symbionts, Grosmanniaclavigera, in western Canada using neutral genetic markers. We examined how spatial heterogeneity affects genetic variation within beetles and fungi and developed a novel integrated landscape genetics approach to assess reciprocal genetic influences between species using constrained ordination. We also compared landscape genetic models built using Euclidean distances based on allele frequencies to traditional pair-wise Fst. Both beetles and fungi exhibited moderate levels of genetic structure over the total study area, low levels of structure in the south, and more pronounced fungal structure in the north. Beetle genetic variation was associated with geographic location while that of the fungus was not. Pinevolume and climate explained beetle genetic variation in the northern region of recent outbreak expansion. Reciprocal genetic relationships were only detectedin the south where there has been alonger history of beetle infestations. The Euclidean distance and Fst-based analyses resulted in similar models in the north and over the entire study area, but differences between methods in the south suggest that genetic distances measures should be selected based on ecological and evolutionary contexts. The integrated landscape genetics framework we present is powerful, general, and can be applied to other systems to quantify the biotic and abiotic determinants of spatial genetic variation within and among taxa.

Rapid identification and detection of pine pathogenic fungi associated with mountain pine beetles by padlock probes.

Fifteen million hectares of pine forests in western Canada have been attacked by the mountain pine beetle (Dendroctonus ponderosae; MPB), leading to devastating economic losses. Grosmannia clavigera and Leptographium longiclavatum, are two fungi intimately associated with the beetles, and are crucial components of the epidemic. To detect and discriminate these two closely related pathogens, we utilized a method based on ligase-mediated nucleotide discrimination with padlock probe technology, and signal amplification by hyperbranched rolling circle amplification (HRCA). Two padlock probes were designed to target species-specific single nucleotide polymorphisms (SNPs) located at the inter-generic spacer 2 region and large subunit of the rRNA respectively, which allows discrimination between the two species. Thirty-four strains of G. clavigera and twenty-five strains of L. longiclavatum representing a broad geographic origin were tested with this assay. The HRCA results were largely in agreement with the conventional identification based on morphology or DNA-based methods. Both probes can also efficiently distinguish the two MPB-associated fungi from other fungi in the MPB, as well as other related fungi in the order Ophiostomatales. We also tested this diagnostic method for the direct detection of these fungi from the DNA of MPB. A nested PCR approach was used to enrich amplicons for signal detection. The results confirmed the presence of these two fungi in MPB. Thus, the padlock probe assay coupled with HRCA is a rapid, sensitive and reproducible method for the identification and detection of these ophiostomatoid fungi.

The effect of hypoxia on gill morphology and ionoregulatory status in the Lake Qinghai scaleless carp, Gymnocypris przewalskii.

Goldfish and crucian carp at low temperature exhibit plasticity in gill morphology during exposure to hypoxia to enhance gas exchange. Hypoxia-induced changes in gill morphology and cellular ultrastructure of the high altitude scaleless carp from Lake Qinghai, China, were investigated to determine whether this is a general characteristic of cold water carp species. Fish were exposed to acute hypoxia (0.3 mg O2 l(-1)) for 24 h followed by 12 h recovery in normoxic water (6 mg O2 l(-1) at 3200 m altitude), with no mortality. Dramatic alterations in gill structure were initiated within 8 h of hypoxia and almost complete by 24 h, and included a gradual reduction of filament epithelial thickness (>50%), elongation of respiratory lamellae, expansion of lamellar respiratory surface area (>60%) and reduction in epithelial water-blood diffusion distance (<50%). An increase in caspase 3 activity in gills occurred following 24 h exposure to hypoxia, indicating possible involvement of apoptosis in gill remodeling. Extensive gill mucous production during hypoxia may have been part of a general stress response or may have played a role in ion exchange and water balance. The large increase in lamellar surface area and reduction in diffusion distance presumably enhances gas transfer during hypoxia (especially in the presence of increased mucous production) but comes with an ionoregulatory cost, as indicated by a 10 and 15% reduction in plasma [Na+] and [Cl-], respectively, within 12-24 h of hypoxia. Within 12 h of hypoxia exposure, ;wavy-convex'-mitochondria rich cells (MRCs) with large apical crypts and numerous branched microvilli were transformed into small ;shallow-basin' cells with a flattened surface. As the apical membrane of MRCs is the site for active ion uptake from the water, a reduction in apical crypt surface area may have contributed to the progressive reduction in plasma [Na+] and [Cl-] observed during hypoxia. The changes in the macro- and ultra-structure of fish gills, and plasma [Na+] and [Cl-] during hypoxia were reversible, showing partial recovery by 12 h following return to normoxia. Although the large morphological changes in the gill observed in the scaleless carp support the hypothesis that gill remodeling during hypoxia is a general characteristic of cold water carp species, the reduced magnitude of the response in scaleless carp relative to goldfish and crucian carp may be a reflection of their more active lifestyle or because they reside in a moderately hypoxic environment at altitude.

Przewalski's naked carp (Gymnocypris przewalskii): an endangered species taking a metabolic holiday in Lake Qinghai, China.

The naked carp is an endangered cyprinid that migrates annually between freshwater rivers, where it spawns, and Lake Qinghai, where it feeds and grows. Lake Qinghai is a high-altitude lake (3,200 m) in western China that currently exhibits the following composition (in mmol L(-1): [Na(+)] 200, [Cl(-)] 173, [Mg(2+)] 36, [Ca(2+)] 0.23, [K(+)] 5.3, total CO(2) 21, titration alkalinity 29; osmolality 375 mOsm kg(-1); pH 9.3), but concentrations are increasing because of water diversion and climate change. We studied the physiology of river water to lake water transfer. When river fish are transferred to lake water, there is a transitory metabolic acidosis followed by a slight respiratory alkalosis, and hemoconcentration occurs. All plasma electrolytes rise over the initial 48 h, and final levels in lake water-acclimated fish are very close to lake water concentrations for [Na(+)], [Cl(-)], [K(+)], and osmolality, whereas [Ca(2+)] continues to be regulated well above ambient levels. However, [Mg(2+)] rises to a much greater extent (fourfold in 48 h); final plasma levels in lake fish may reach 12 mmol L(-1) but are still much lower than in lake water (36 mmol L(-1)). At the same time, urine flow rate decreases drastically to <5% of river water values; only the renal excretion of Mg(2+) is maintained. Both gill and kidney Na(+),K(+)-ATPase rapidly decline, with final levels in lake water fish only 30% and 70%, respectively, of those in river water fish. Metabolic rate also quickly decreases on exposure to lake water, with O(2) consumption and ammonia-N excretion rates eventually falling to only 60% and 30%, respectively, of those in river fish, while plasma ammonia rises fivefold. The fish appear to be benefiting from a metabolic holiday at present because of decreases in iono- and osmoregulatory costs while in lake water; elevated plasma [Mg(2+)] and ammonia may be additional factors depressing metabolic rate. If the lake continues to dehydrate, these benefits may change to pathology.

[The significance of mitochondria control region (D-Loop) in intraspecific genetic differentiation of fish].

mtDNA has been widely used as a molecular marker to investigate the phylogenesis of many species. D-Loop has important value for investigating intraspecific genetic differentiation because of its high mutation accumulation. Fish is the most primitive vertebrate, but by far is the most dominant in number. The distribution of fish is broad and its origin is very complicated, has always been a very interesting subject for investigation. D-Loop has many important implications for investigating the intraspecific genetic differentiation of fish. Recently, much insight has been gained by the increasing use of D-Loop as a molecular marker for investigating the intraspecific genetic differentiation of fish. Lake Qinghai is the biggest inland saline lake in China. The main fish in the lake is scaleless carp (Gymnocypris przewalskii). Initial analysis of its D-Loop showed no intraspecific genetic differentiation among the populations in the lake and the rivers around.

Expression of four glutamine synthetase genes in the early stages of development of rainbow trout (Oncorhynchus mykiss) in relationship to nitrogen excretion.

The incorporation of ammonia into glutamine, catalyzed by glutamine synthetase, is thought to be important in the detoxification of ammonia in animals. During early fish development, ammonia is continuously formed as yolk proteins and amino acids are catabolized. We followed the changes in ammonia and urea-nitrogen content, ammonia and urea-nitrogen excretion, glutamine synthetase activity, and mRNA expression of four genes coding for glutamine synthetase (Onmy-GS01-GS04) over 3-80 days post fertilization and in adult liver and skeletal muscle of the rainbow trout (Oncorhynchus mykiss). Both ammonia and urea-nitrogen accumulate before hatching, although the rate of ammonia excretion is considerably higher relative to urea-nitrogen excretion. All four genes were expressed during early development, but only Onmy-GS01 and -GS02 were expressed at appreciable levels in adult liver, and expression was very low in muscle tissue. The high level of expression of Onmy-GS01 and -GS03 prior to hatching corresponded to a linear increase in glutamine synthetase activity. We propose that the induction of glutamine synthetase genes early in development and the subsequent formation of the active protein are preparatory for the increased capacity of the embryo to convert the toxic nitrogen end product, ammonia, into glutamine, which may then be utilized in the ornithine-urea cycle or other pathways.

Dogmas and controversies in the handling of nitrogenous wastes: expression of arginase Type I and II genes in rainbow trout: influence of fasting on liver enzyme activity and mRNA levels in juveniles.

Through analysis of a cDNA library and third-party annotation of available database sequences, we characterized the full-length coding regions of rainbow trout (Oncorhynchus mykiss) Type I, Onmy-ARG01, and Type II, Onmy-ARG02, arginase genes. Two partial related arginase sequences, Onmy-ARG01b and Onmy-ARG02b, and a full-length zebrafish arginase coding region (Danio rerio), Dare-ARG02, are also reported. Comparison of vertebrate arginase sequences shows that both Type I and Type II genes in bony fishes contain a mitochondrial targeting N-terminal domain. This suggests that the cytosolic Type I arginase found in ureotelic vertebrates arose in the common ancestor of amphibia and mammals. Onmy-ARG01 and Onmy-ARG02 mRNA was detected in liver, kidney, gill, intestine, red muscle and heart tissues. Onmy-ARG01 was expressed at a significantly higher level relative to Onmy-ARG02 in liver and red muscle tissue. To investigate whether there was differential regulation of Onmy-ARG01 and Onmy-ARG02, juvenile trout were fasted for 6 weeks and hepatic enzyme activities and mRNA levels were compared with those of fed control fish. There was a 3-fold increase in liver arginase activity and a 2-fold increase in Onmy-ARG02 mRNA levels but no change in Onmy-ARG01 mRNA levels in fasted fish relative to fed fish. These findings indicate that both types of arginase genes are present and expressed in rainbow trout and that the pattern of expression varies between tissues. The increase in liver arginase activity after a 6-week fast is due, in part, to an increase in the expression of Onmy-ARG02 mRNA levels.

Evolution of glutamine synthetase in vertebrates: multiple glutamine synthetase genes expressed in rainbow trout (Oncorhynchus mykiss).

Glutamine synthetase (GSase) is a key enzyme in nitrogen metabolism and encoded by a single gene in mammals. Using PCR cloning techniques, including RT-PCR from total RNA and PCR from a cDNA library, we find evidence of four expressed GSase mRNAs for the tetraploid rainbow trout. For two of these mRNAs (Onmy-GS01, -GS02) we characterize the full-length coding regions, and for two others (Onmy-GS03, -GS04), we describe partial sequences. Northern analysis of Onmy-GS01, -GS02, -GS03 and -GS04 indicates that (1) Onmy-GS02 is expressed at higher levels relative to the other transcripts in most adult tissues, with the exception of brain and gill, where Onmy-GS01 is at the highest level, and (2) the tissue with the highest level of expression of all four transcripts is the brain, with decreasing levels in the intestine, liver, red muscle, gill/kidney, white muscle and heart. Clearly, rainbow trout possess multiple GSase genes with differing levels of tissue expression, implying manifold potential routes of regulation for this octameric enzyme. Our data also indicate that caution should be taken when interpreting mRNA expression data of a single gene, unless multiple genes have been ruled out. Consistent with a southern blot, phylogenetic and intron sequence analyses imply that the trout genes are encoded by at least four separate loci, belonging to two distinct evolutionary branches. Our data on rainbow trout, together with those from two full-length zebrafish Danio rerio GSase genes compiled from GenBank ESTs, support the idea that fish GSases are polyphyletic and that gene duplications have occurred at multiple points and in independent lineages throughout the evolution of bony fishes.

Linkage Relationships and Haplotype Variation of the Major Histocompatibility Complex Class I A Genes in the Cichlid Fish Oreochromis niloticus.

The haplochromine cichlid species flocks of the East African Great Lakes are one of the best examples of adaptive radiation. Analysis of genetic variation among these species provides valuable information on species relationships and timing of speciation events. Although the haplochromine cichlids generally display little genetic variation, the major histocompatibility complex (Mhc) genes have been found to be highly variable. A study of the linkage relationships of the Mhc class I A genes in the cichlid fish Oreochromis niloticus was therefore undertaken. Class I loci were identified, and their segregation in seven mothers and their haploid embryos was determined. In total, 56 class I A sequences were found among the seven families. A strong concordance of segregation was observed in five haplotypes among the embryos, indicating a close linkage of all loci. The number of loci per haplotype varied from 11 to 17, while the total number of distinct loci found among all families was 22. These findings show that all class I A loci are linked in a single genetic cluster in O. niloticus.