Northern glacial refugia for the pygmy shrew Sorex minutus in Europe revealed by phylogeographic analyses and species distribution modelling

The southern European peninsulas (Iberian, Italian and Balkan) are traditionally recognized as glacial refugia from where many species colonized central and northern Europe after the Last Glacial Maximum (LGM). However, evidence that some species had more northerly refugia is accumulating from phylogeographic, palaeontological and palynological studies, and more recently from species distribution modelling (SDM), but further studies are needed to test the idea of northern refugia in Europe. Here, we take a rarely implemented multidisciplinary approach to assess if the pygmy shrew Sorex minutus, a widespread Eurasian mammal species, had northern refugia during the LGM, and if these influenced its postglacial geographic distribution. First, we evaluated the phylogeographic and population expansion patterns using mtDNA sequence data from 123 pygmy shrews. Then, we used SDM to predict present and past (LGM) potential distributions using two different training data sets, two different algorithms (Maxent and GARP) and climate reconstructions for the LGM with two different general circulation models. An LGM distribution in the southern peninsulas was predicted by the SDM approaches, in line with the occurrence of lineages of S. minutus in these areas. The phylogeographic analyses also indicated a widespread and strictly northern‐central European lineage, not derived from southern peninsulas, and with a postglacial population expansion signature. This was consistent with the SDM predictions of suitable LGM conditions for S. minutus occurring across central and eastern Europe, from unglaciated parts of the British Isles to much of the eastern European Plain. Hence, S. minutus likely persisted in parts of central and eastern Europe during the LGM, from where it colonized other northern areas during the late‐glacial and postglacial periods. Our results provide new insights into the glacial and postglacial colonization history of the European mammal fauna, notably supporting glacial refugia further north than traditionally recognized.     

Studies of metabolic phenotypic correlates of 15 obesity associated gene variants

Aims

Genome-wide association studies have identified novel BMI/obesity associated susceptibility loci. The purpose of this study is to determine associations with overweight, obesity, morbid obesity and/or general adiposity in a Danish population. Moreover, we want to investigate if these loci associate with type 2 diabetes and to elucidate potential underlying metabolic mechanisms.

Methods

15 gene variants in 14 loci including TMEM18 (rs7561317), SH2B1 (rs7498665), KCTD15 (rs29941), NEGR1 (rs2568958), ETV5 (rs7647305), BDNF (rs4923461, rs925946), SEC16B (rs10913469), FAIM2 (rs7138803), GNPDA2 (rs10938397), MTCH2 (rs10838738), BAT2 (rs2260000), NPC1 (rs1805081), MAF (rs1424233), and PTER (rs10508503) were genotyped in 18,014 middle-aged Danes.

Results

Five of the 15 gene variants associated with overweight, obesity and/or morbid obesity. Per allele ORs ranged from 1.15–1.20 for overweight, 1.10–1.25 for obesity, and 1.41–1.46 for morbid obesity. Five of the 15 variants moreover associated with increased measures of adiposity. BDNF rs4923461 displayed a borderline BMI-dependent protective effect on type 2 diabetes (0.87 (0.78–0.96, p = 0.008)), whereas SH2B1 rs7498665 associated with nominally BMI-independent increased risk of type 2 diabetes (1.16 (1.07–1.27, p = 7.8×10⁻⁴)).

Conclusions

Associations with overweight and/or obesity and measures of obesity were confirmed for seven out of the 15 gene variants. The obesity risk allele of BDNF rs4923461 protected against type 2 diabetes, which could suggest neuronal and peripheral distinctive ways of actions for the protein. SH2B1 rs7498665 associated with type 2 diabetes independently of BMI.     

Further evidence for the role of ENPP1 in obesity: association with morbid obesity in Finns

The aim of this study was to investigate a series of single-nucleotide polymorphisms (SNPs) in the genes MC2R, MC3R, MC4R, MC5R, POMC, and ENPP1 for association with obesity. Twenty-five SNPs (2-7 SNPs/gene) were genotyped in 246 Finns with extreme obesity (BMI > or = 40 kg/m2) and in 481 lean subjects (BMI 20-25 kg/m2). Of the obese subjects, 23% had concomitant type 2 diabetes. SNPs and SNP haplotypes were tested for association with obesity and type 2 diabetes. Allele frequencies differed between obese and lean subjects for two SNPs in the ENPP1 gene, rs1800949 (P = 0.006) and rs943003 (P = 0.0009). These SNPs are part of a haplotype (rs1800949 C-rs943003 A), which was observed more frequently in lean subjects compared to obese subjects (P = 0.0007). Weaker associations were detected between the SNPs rs1541276 in the MC5R, rs1926065 in the MC3R genes and obesity (P = 0.04 and P = 0.03, respectively), and between SNPs rs2236700 in the MC5R, rs2118404 in the POMC, rs943003 in the ENPP1 genes and type 2 diabetes (P = 0.03, P = 0.02 and P = 0.02, respectively); these associations did not, however, remain significant after correction for multiple testing. In conclusion, a previously unexplored ENPP1 haplotype composed of SNPs rs1800949 and rs943003 showed suggestive evidence for association with adult-onset morbid obesity in Finns. In this study, we did not find association between the frequently studied ENPP1 K121Q variant, nor SNPs in the MCR or POMC genes and obesity or type 2 diabetes.     

Mesopredator behavioral response to olfactory signals of an apex predator

Olfactory signals constitute an important mechanism in interspecific interactions, but little is known regarding their role in communication between predator species. We analyzed the behavioral responses of a mesopredator, the red fox (Vulpes vulpes), to an olfactory cue (scat) of an apex predator, the lynx (Lynx lynx) in Bia?owie?a Primeval Forest, Poland, using video camera traps. Red fox visited sites with scats more often than expected and the duration of their visits was longer at scat sites than at control sites (no scat added). Vigilant behavior, sniffing and scent marking (including over-marking) occurred more often at scat sites compared to control sites, where foxes mainly passed by. Vigilance was most pronounced during the first days of the recordings. Red fox behavior was also influenced by foxes previously visiting scat sites. They sniffed and scent marked (multiple over-marking) more frequently when the lynx scat had been over-marked previously by red fox. Fox visits to lynx scats may be seen as a trade-off between obtaining information on a potential food source (prey killed by lynx) and the potential risk of predation by an apex predator.     

Covalent binding of the 13C-labeled skin sensitizers 5-chloro-2-methylisothiazol-3-one (MCI) and 2-methylisothiazol-3-one (MI) to a model peptide and glutathione

The reactivity of 4-[13C]- and 5-[13C]-5-chloro-2-methylisothiazol-3-one (MCI) and 2-methylisothiazol-3-one (MI) towards a model peptide and glutathione was followed by 13C and 1H[13C] NMR spectroscopy. Both molecules were found to react with GSH but in addition MCI was found to react with histidine and lysine to form adducts of a different nature. Reaction with histidine led to stable substitution adducts through an addition-elimination reaction at position 5 while reaction with lysine led to the formation of open adducts of the thioamide or amide type.     

Analysis of kavalactones from Piper methysticum (kava-kava)

The chemical analysis and quality control of both Piper methysticum G. Forster (kava-kava) and extracts obtained by aqueous acetone or aqueous methanol as well as supercritical fluid extraction are reviewed. In the last two decades various procedures concerning the separation and detection of kavalactones have been routinely carried out by gas chromatography (without previous derivatization of kavalactones) and high performance liquid chromatography but most of them are not validated or only partially validated. Recently, analyses by supercritical fluid chromatography and micellar electrokinetic chromatography have also been reported. Both gas chromatography and high performance liquid chromatography can be used for the analysis of kavalactones with some advantages and disadvantages for each method. Using gas chromatography analysis, methysticin and yangonin, which are two of the major components, are generally not separated. In addition, the high temperature of the injection port caused the decomposition of methysticin. Concerning high performance liquid chromatography analyses, the reversed-phase is generally better because highly reproducible with a very low detection limit for all compounds even if the quantitative analysis of the kavalactones by liquid chromatography needs to be carried out in the absence of light to prevent the cis/trans isomerisation of yangonin.     

Primary sequence determinants responsible for site-selective dephosphorylation of the PDGF beta-receptor by the receptor-like protein tyrosine phosphatase DEP-1

Site-selective dephosphorylation of receptor tyrosine kinases contributes to receptor regulation. The receptor-like protein tyrosine phosphatase DEP-1 site-selectively dephosphorylates the PDGF beta-receptor. DEP-1 dephosphorylation of original and chimeric phospho-peptides spanning the preferred pY1021 and the less preferred pY857 and pY562 sites was analyzed. Double substitutions of basic residues at -4 and +3 of pY857 and pY562 peptides improved affinity. Substitutions of single amino acids indicated preference for an acidic residue at position -1 and a preference against a basic residue at position +3. DEP-1 site-selective dephosphorylation of PDGF beta-receptor is thus determined by the primary sequence surrounding phosphorylation sites and involves interactions with residues spanning at least between positions -1 and +3.     

Suppressive subtractive hybridization detects extensive genomic diversity in Thermotoga maritima

Comparisons between genomes of closely related bacteria often show large variations in gene content, even between strains of the same species. Such studies have focused mainly on pathogens; here, we examined Thermotoga maritima, a free-living hyperthermophilic bacterium, by using suppressive subtractive hybridization. The genome sequence of T. maritima MSB8 is available, and DNA from this strain served as a reference to obtain strain-specific sequences from Thermotoga sp. strain RQ2, a very close relative (approximately 96% identity for orthologous protein-coding genes, 99.7% identity in the small-subunit rRNA sequence). Four hundred twenty-six RQ2 subtractive clones were sequenced. One hundred sixty-six had no DNA match in the MSB8 genome. These differential clones comprise, in sum, 48 kb of RQ2-specific DNA and match 72 genes in the GenBank database. From the number of identical clones, we estimated that RQ2 contains 350 to 400 genes not found in MSB8. Assuming a similar genome size, this corresponds to 20% of the RQ2 genome. A large proportion of the RQ2-specific genes were predicted to be involved in sugar transport and polysaccharide degradation, suggesting that polysaccharides are more important as nutrients for this strain than for MSB8. Several clones encode proteins involved in the production of surface polysaccharides. RQ2 encodes multiple subunits of a V-type ATPase, while MSB8 possesses only an F-type ATPase. Moreover, an RQ2-specific MutS homolog was found among the subtractive clones and appears to belong to a third novel archaeal type MutS lineage. Southern blot analyses showed that some of the RQ2 differential sequences are found in some other members of the order Thermotogales, but the distribution of these variable genes is patchy, suggesting frequent lateral gene transfer within the group.     

Distribution of spontaneous mutants and inferences about the replication mode of the RNA bacteriophage phi6

When a parent virus replicates inside its host, it must first use its own genome as the template for replication. However, once progeny genomes are produced, the progeny can in turn act as templates. Depending on whether the progeny genomes become templates, the distribution of mutants produced by an infection varies greatly. While information on the distribution is important for many population genetic models, it is also useful for inferring the replication mode of a virus. We have analyzed the distribution of mutants emerging from single bursts in the RNA bacteriophage phi6 and find that the distribution closely matches a Poisson distribution. The match suggests that replication in this bacteriophage is effectively by a stamping machine model in which the parental genome is the main template used for replication. However, because the distribution deviates slightly from a Poisson distribution, the stamping machine is not perfect and some progeny genomes must replicate. By fitting our data to a replication model in which the progeny genomes become replicative at a given rate or probability per round of replication, we estimated the rate to be very low and on the on the order of 10(-4). We discuss whether different replication modes may confer an adaptive advantage to viruses.