High-Resolution Taxonomic Profiling of the Subgingival Microbiome for Biomarker Discovery and Periodontitis Diagnosis

The oral microbiome plays a key role for caries, periodontitis, and systemic diseases. A method for rapid, high-resolution, robust taxonomic profiling of subgingival bacterial communities for early detection of periodontitis biomarkers would therefore be a useful tool for individualized medicine. Here, we used Illumina sequencing of the V1-V2 and V5-V6 hypervariable regions of the 16S rRNA gene. A sample stratification pipeline was developed in a pilot study of 19 individuals, 9 of whom had been diagnosed with chronic periodontitis. Five hundred twenty-three operational taxonomic units (OTUs) were obtained from the V1-V2 region and 432 from the V5-V6 region. Key periodontal pathogens like Porphyromonas gingivalis, Treponema denticola, and Tannerella forsythia could be identified at the species level with both primer sets. Principal coordinate analysis identified two outliers that were consistently independent of the hypervariable region and method of DNA extraction used. The linear discriminant analysis (LDA) effect size algorithm (LEfSe) identified 80 OTU-level biomarkers of periodontitis and 17 of health. Health- and periodontitis-related clusters of OTUs were identified using a connectivity analysis, and the results confirmed previous studies with several thousands of samples. A machine learning algorithm was developed which was trained on all but one sample and then predicted the diagnosis of the left-out sample (jackknife method). Using a combination of the 10 best biomarkers, 15 of 17 samples were correctly diagnosed. Training the algorithm on time-resolved community profiles might provide a highly sensitive tool to detect the onset of periodontitis.     

Sex-driven differences in Egyptian mongoose’s (<Emphasis Type="Italic">Herpestes ichneumon</Emphasis>) diet in its northwestern European range

Sexual dimorphism in carnivores can result from, or induce, variations in diet and foraging behaviour between individuals of different sexes. Sex-driven behavioural changes in feeding habits may also result in a reduction in intraspecific competition for resources, avoiding dietary overlap by concentrating on different prey sizes/types. We therefore evaluated the variation in feeding habits of both males and females of an invasive and range expanding southern European carnivore, the Egyptian mongoose (Herpestes ichneumon), through gut content analysis of road or predator-control killed animals. The analysis of 59 gut samples, collected in Portugal, revealed that northwestern populations of mongooses prey mostly upon mammals [especially lagomorphs, percentage of biomass (PB) = 44%], reptiles (PB = 28%) and arthropods (PB = 2% but percentage of occurrence = 35%). However, females seem to focus mainly on reptiles (PB = 51%) and mammals (PB = 38%), whilst males consume mostly mammals (PB = 75%). The results suggest that this variation maybe a result of sexual dimorphism (since some variation is documented in a few Egyptian mongoose’s skull measures, namely condylobasal length and canine diameter), with males shifting their diets to prey from which they can obtain more energy. Females seem to be focused on predating lighter animals, probably females or juveniles, a fact that has important game management implications.     

Docking and scoring with alternative side-chain conformations

We describe a scoring and modeling procedure for docking ligands into protein models that have either modeled or flexible side-chain conformations. Our methodical contribution comprises a procedure for generating new potentials of mean force for the ROTA scoring function which we have introduced previously for optimizing side-chain conformations with the tool IRECS. The ROTA potentials are specially trained to tolerate small-scale positional errors of atoms that are characteristic of (i) side-chain conformations that are modeled using a sparse rotamer library and (ii) ligand conformations that are generated using a docking program. We generated both rigid and flexible protein models with our side-chain prediction tool IRECS and docked ligands to proteins using the scoring function ROTA and the docking programs FlexX (for rigid side chains) and FlexE (for flexible side chains). We validated our approach on the forty screening targets of the DUD database. The validation shows that the ROTA potentials are especially well suited for estimating the binding affinity of ligands to proteins. The results also show that our procedure can compensate for the performance decrease in screening that occurs when using protein models with side chains modeled with a rotamer library instead of using X-ray structures. The average runtime per ligand of our method is 168 seconds on an Opteron V20z, which is fast enough to allow virtual screening of compound libraries for drug candidates.     

Heavy-Metal Stress and Developmental Patterns of Arbuscular Mycorrhizal Fungi

The rate of global deposition of Cd, Pb, and Zn has decreased over the past few decades, but heavy metals already in the soil may be mobilized by local and global changes in soil conditions and exert toxic effects on soil microorganisms. We examined in vitro effects of Cd, Pb, and Zn on critical life stages in metal-sensitive ecotypes of arbuscular mycorrhizal (AM) fungi, including spore germination, presymbiotic hyphal extension, presymbiotic sporulation, symbiotic extraradical mycelium expansion, and symbiotic sporulation. Despite long-term culturing under the same low-metal conditions, two species, Glomus etunicatum and Glomus intraradices, had different levels of sensitivity to metal stress. G. etunicatum was more sensitive to all three metals than was G. intraradices. A unique response of increased presymbiotic hyphal extension occurred in G. intraradices exposed to Cd and Pb. Presymbiotic hyphae of G. intraradices formed presymbiotic spores, whose initiation was more affected by heavy metals than was presymbiotic hyphal extension. In G. intraradices grown in compartmentalized habitats with only a portion of the extraradical mycelium exposed to metal stress, inhibitory effects of elevated metal concentrations on symbiotic mycelial expansion and symbiotic sporulation were limited to the metal-enriched compartment. Symbiotic sporulation was more sensitive to metal exposure than symbiotic mycelium expansion. Patterns exhibited by G. intraradices spore germination, presymbiotic hyphal extension, symbiotic extraradical mycelium expansion, and sporulation under elevated metal concentrations suggest that AM fungi may be able to survive in heavy metal-contaminated environments by using a metal avoidance strategy.     

Dynamic modulation of prohormone convertase 2 (PC2)-mediated precursor processing by 7B2 protein: preferential effect on glucagon synthesis

The small neuroendocrine protein 7B2 is required for the production of active prohormone convertase 2 (PC2), an enzyme involved in the synthesis of peptide hormones, such as glucagon and proopiomelanocortin-derived α-melanocyte-stimulating hormone. However, whether 7B2 can dynamically modulate peptide production through regulation of PC2 activity remains unclear. Infection of the pancreatic alpha cell line α-TC6 with 7B2-encoding adenovirus efficiently increased production of glucagon, whereas siRNA-mediated knockdown of 7B2 significantly decreased stored glucagon. Furthermore, rescue of 7B2 expression in primary pituitary cultures prepared from 7B2 null mice restored melanocyte-stimulating hormone production, substantiating the role of 7B2 as a regulatory factor in peptide biosynthesis. In anterior pituitary and pancreatic beta cell lines, however, overexpression of 7B2 affected neither production nor secretion of peptides despite increased release of active PC2. In direct contrast, 7B2 overexpression decreased the secretion and increased the activity of PC2 within α-TC6 cells; the increased intracellular concentration of active PC2 within these cells may therefore account for the enhanced production of glucagon. In line with these findings, we found elevated circulating glucagon levels in 7B2-overexpressing cast/cast mice in vivo. Surprisingly, when proopiomelanocortin and proglucagon were co-expressed in either pituitary or pancreatic alpha cell lines, proglucagon processing was preferentially decreased when 7B2 was knocked down. Taken together, these results suggest that proglucagon cleavage has a greater dependence on PC2 activity than other precursors and moreover that 7B2-dependent routing of PC2 to secretory granules is cell line-specific. The manipulation of 7B2 could therefore represent an effective way to selectively regulate synthesis of certain PC2-dependent peptides.     

Habitat traits and food availability determine the response of marine invertebrates to ocean acidification

Energy availability and local adaptation are major components in mediating the effects of ocean acidification (OA) on marine species. In a long‐term study, we investigated the effects of food availability and elevated pCO₂ (ca. 400, 1000 and 3000 μatm) on growth of newly settled Amphibalanus (Balanus) improvisus to reproduction, and on their offspring. We also compared two different populations, which were presumed to differ in their sensitivity to pCO₂ due to differing habitat conditions: Kiel Fjord, Germany (Western Baltic Sea) with naturally strong pCO₂ fluctuations, and the Tjärnö Archipelago, Sweden (Skagerrak) with far lower fluctuations. Over 20 weeks, survival, growth, reproduction and shell strength of Kiel barnacles were all unaffected by elevated pCO₂, regardless of food availability. Moulting frequency and shell corrosion increased with increasing pCO₂ in adults. Larval development and juvenile growth of the F1 generation were tolerant to increased pCO₂, irrespective of parental treatment. In contrast, elevated pCO₂ had a strong negative impact on survival of Tjärnö barnacles. Specimens from this population were able to withstand moderate levels of elevated pCO₂ over 5 weeks when food was plentiful but showed reduced growth under food limitation. Severe levels of elevated pCO₂ negatively impacted growth of Tjärnö barnacles in both food treatments. We demonstrate a conspicuously higher tolerance to elevated pCO₂ in Kiel barnacles than in Tjärnö barnacles. This tolerance was carried over from adults to their offspring. Our findings indicate that populations from fluctuating pCO₂ environments are more tolerant to elevated pCO₂ than populations from more stable pCO₂ habitats. We furthermore provide evidence that energy availability can mediate the ability of barnacles to withstand moderate CO₂ stress. Considering the high tolerance of Kiel specimens and the possibility to adapt over many generations, near future OA alone does not seem to present a major threat for A. improvisus.     

Mutations of the PC2 substrate binding pocket alter enzyme specificity

By taking advantage of the recently published furin structure, whose catalytic domain shares high homology with other proprotein convertases, we designed mutations in the catalytic domain of PC2, altering residues Ser206, Thr271, Asp278, ArgGlu282, AlaSer323, Leu341, Asn365, and Ser380, which are both conserved and specific to this convertase, and substituting residues specific to PC1 and/or furin. In order to investigate the determinants of PC2 specificity, we have tested the mutated enzymes against a set of proenkephalin-derived substrates, as well as substrates representing Arg, Ala, Leu, Phe, and Glu positional scanning variants of a peptide B-derived substrate. We found that the exchange of the Ser206 residue with Arg or Lys led to a total loss of activity. Increased positive charge of the substrate generally resulted in an increased specificity constant. Most intriguingly, the RE281GR mutation, corresponding to a residue placed distantly in the S6 pocket, evoked the largest changes in the specificity pattern. The D278E and N356S mutations resulted in distinct alterations in PC2 substrate preferences. However, when other residues that distinguish PC2 from other convertases were substituted with PC1-like or furin-like equivalents, there was no significant alteration of the PC2 specificity pattern, suggesting that the overall structure of the substrate binding cleft rather than individual residues specifies substrate binding.     

Molecular mechanism of 14-3-3 protein-mediated inhibition of plant nitrate reductase

14-3-3 proteins regulate key processes in eukaryotic cells including nitrogen assimilation in plants by tuning the activity of nitrate reductase (NR), the first and rate-limiting enzyme in this pathway. The homodimeric NR harbors three cofactors, each of which is bound to separate domains, thus forming an electron transfer chain. 14-3-3 proteins inhibit NR by binding to a conserved phosphorylation site localized in the linker between the heme and molybdenum cofactor-containing domains. Here, we have investigated the molecular mechanism of 14-3-3-mediated NR inhibition using a fragment of the enzyme lacking the third domain, allowing us to analyze electron transfer from the heme cofactor via the molybdenum center to nitrate. The kinetic behavior of the inhibited Mo-heme fragment indicates that the principal point at which 14-3-3 acts is the electron transfer from the heme to the molybdenum cofactor. We demonstrate that this is not due to a perturbation of the reduction potentials of either the heme or the molybdenum center and conclude that 14-3-3 most likely inhibits nitrate reductase by inducing a conformational change that significantly increases the distance between the two redox-active sites.     

Evolutionary Dynamics of the Leucine-Rich Repeat Receptor-Like Kinase (LRR-RLK) Subfamily in Angiosperms

Gene duplications are an important factor in plant evolution, and lineage-specific expanded (LSE) genes are of particular interest. Receptor-like kinases expanded massively in land plants, and leucine-rich repeat receptor-like kinases (LRR-RLK) constitute the largest receptor-like kinases family. Based on the phylogeny of 7,554 LRR-RLK genes from 31 fully sequenced flowering plant genomes, the complex evolutionary dynamics of this family was characterized in depth. We studied the involvement of selection during the expansion of this family among angiosperms. LRR-RLK subgroups harbor extremely contrasting rates of duplication, retention, or loss, and LSE copies are predominantly found in subgroups involved in environmental interactions. Expansion rates also differ significantly depending on the time when rounds of expansion or loss occurred on the angiosperm phylogenetic tree. Finally, using a d_N/d_S-based test in a phylogenetic framework, we searched for selection footprints on LSE and single-copy LRR-RLK genes. Selective constraint appeared to be globally relaxed at LSE genes, and codons under positive selection were detected in 50% of them. Moreover, the leucine-rich repeat domains, and specifically four amino acids in them, were found to be the main targets of positive selection. Here, we provide an extensive overview of the expansion and evolution of this very large gene family.Receptor-like kinases (RLKs) constitute one of the largest gene families in plants and expanded massively in land plants (Embryophyta; Lehti-Shiu et al., 2009, 2012). For plant RLK gene families, the functions of most members are often not known (especially in recently expanded families), but some described functions include innate immunity (Albert et al., 2010), pathogen response (Dodds and Rathjen, 2010), abiotic stress (Yang et al., 2010), development (De Smet et al., 2009), and sometimes multiple functions (Lehti-Shiu et al., 2012). The RLKs usually consist of three domains: an N-terminal extracellular domain, a transmembrane domain, and a C-terminal kinase domain (KD). In plants, the KD usually has a Ser/Thr specificity (Shiu and Bleecker, 2001), but Tyr-specific RLKs were also described (e.g. BRASSINOSTEROID INSENSITIVE1; Oh et al., 2009). Interestingly, it was estimated that approximately 20% of RLKs contain a catalytically inactive KD (e.g. STRUBBELIG and CORYNE; Chevalier et al., 2005; Castells and Casacuberta, 2007; Gish and Clark, 2011). In Arabidopsis (Arabidopsis thaliana), 44 RLK subgroups (SGs) were defined by inferring the phylogenetic relationships between the KDs (Shiu and Bleecker, 2001). Interestingly, different SGs show different duplication/retention rates (Lehti-Shiu et al., 2009). Specifically, RLKs involved in stress responses show a high number of tandemly duplicated genes whereas those involved in development do not (Shiu et al., 2004), which suggests that some RLK genes are important for the responses of land plants to a changing environment (Lehti-Shiu et al., 2012). There seem to be relatively few RLK pseudogenes compared with other large gene families, and copy retention was argued to be driven by both drift and selection (Zou et al., 2009; Lehti-Shiu et al., 2012). As most SGs are relatively old and RLK subfamilies expanded independently in several plant lineages, duplicate retention cannot be explained by drift alone, and natural selection is expected to be an important driving factor in RLK gene family retention (Lehti-Shiu et al., 2009).Leucine-rich repeat-receptor-like kinases (LRR-RLKs), which contain up to 30 leucine-rich repeat (LRRs) in their extracellular domain, constitute the largest RLK family (Shiu and Bleecker, 2001). Based on the KD, 15 LRR-RLK SGs have been established in Arabidopsis (Shiu et al., 2004; Lehti-Shiu et al., 2009). So far, two major functions have been attributed to them: defense against pathogens and development (Tang et al., 2010b). LRR-RLKs involved in defense are predominantly found in lineage-specific expanded (LSE) gene clusters, whereas LRR-RLKs involved in development are mostly found in nonexpanded groups (Tang et al., 2010b). It was also discovered that the LRR domains are significantly less conserved than the remaining domains of the LRR-RLK genes (Tang et al., 2010b). In addition, a study of four plant genomes (Arabidopsis, grape [Vitis vinifera], poplar [Populus trichocarpa], and rice [Oryza sativa]) showed that LRR-RLK genes from LSE gene clusters show significantly more indications of positive selection or relaxed constraint than LRR-RLKs from nonexpanded groups (Tang et al., 2010b).The genomes of flowering plants (angiosperms) have been shown to be highly dynamic compared with most other groups of land plants (Leitch and Leitch, 2012). This dynamic is mostly caused by the frequent multiplication of genetic material, followed by a complex pattern of differential losses (i.e. the fragmentation process) and chromosomal rearrangements (Langham et al., 2004; Leitch and Leitch, 2012). Most angiosperm genomes sequenced so far show evidence for at least one whole-genome multiplication event during their evolution (Jaillon et al., 2007; D’Hont et al., 2012; Tomato Genome Consortium, 2012). At a smaller scale, tandem and segmental duplications are also very common in angiosperms (Arabidopsis Genome Initiative, 2000; International Rice Genome Sequencing Project, 2005; Rizzon et al., 2006). Although the most common fate of duplicated genes is to be progressively lost, in some cases they can be retained in the genome, and adaptive as well as nonadaptive scenarios have been discussed to play a role in this preservation process (for review, see Moore and Purugganan, 2005; Hahn, 2009; Innan, 2009; Innan and Kondrashov, 2010). Whole-genome sequences also revealed that the same gene may undergo several rounds of duplication and retention. These LSE genes were shown to evolve under positive selection more frequently than single-copy genes in angiosperms (Fischer et al., 2014). That study analyzed general trends over whole genomes. Here, we ask if, and to what extent, this trend is observable at LRR-RLK genes. As this gene family is very dynamic and large, and in accordance with the results of Tang et al. (2010b), we expect the effect of positive selection to be even more pronounced than in the whole-genome average.We analyzed 33 Embryophyta genomes to investigate the evolutionary history of the LRR-RLK gene family in a phylogenetic framework. Twenty LRR-RLK SGs were identified, and from this data set, we deciphered the evolutionary dynamics of this family within angiosperms. The expansion/reduction rates were contrasted between SGs and species as well as in ancestral branches of the angiosperm phylogeny. We then focused on genes whose number increased dramatically in an SG- and/or species-specific manner (i.e. LSE genes). Those genes are likely to be involved in species-specific cellular processes or adaptive interactions and were used as a template to infer the potential occurrence of positive selection. This led to the identification of sites at which positive selection likely acted. We discuss our results in the light of angiosperm genome evolution and current knowledge of LRR-RLK functions. Positive selection footprints identified in LSE genes highlight the importance of combining evolutionary analysis and functional knowledge to guide further investigations.