Temporal variability is a personalized feature of the human microbiome

Background

It is now apparent that the complex microbial communities found on and in the human body vary across individuals. What has largely been missing from previous studies is an understanding of how these communities vary over time within individuals. To the extent to which it has been considered, it is often assumed that temporal variability is negligible for healthy adults. Here we address this gap in understanding by profiling the forehead, gut (fecal), palm, and tongue microbial communities in 85 adults, weekly over 3 months.

Results

We found that skin (forehead and palm) varied most in the number of taxa present, whereas gut and tongue communities varied more in the relative abundances of taxa. Within each body habitat, there was a wide range of temporal variability across the study population, with some individuals harboring more variable communities than others. The best predictor of these differences in variability across individuals was microbial diversity; individuals with more diverse gut or tongue communities were more stable in composition than individuals with less diverse communities.

Conclusions

Longitudinal sampling of a relatively large number of individuals allowed us to observe high levels of temporal variability in both diversity and community structure in all body habitats studied. These findings suggest that temporal dynamics may need to be considered when attempting to link changes in microbiome structure to changes in health status. Furthermore, our findings show that, not only is the composition of an individual’s microbiome highly personalized, but their degree of temporal variability is also a personalized feature.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0531-y) contains supplementary material, which is available to authorized users.     

Efficacy of Distortion Correction on Diffusion Imaging: Comparison of FSL Eddy and Eddy_Correct Using 30 and 60 Directions Diffusion Encoding

Diffusion imaging is a unique noninvasive tool to detect brain white matter trajectory and integrity in vivo. However, this technique suffers from spatial distortion and signal pileup or dropout originating from local susceptibility gradients and eddy currents. Although there are several methods to mitigate these problems, most techniques can be applicable either to susceptibility or eddy-current induced distortion alone with a few exceptions. The present study compared the correction efficiency of FSL tools, “eddy_correct” and the combination of “eddy” and “topup” in terms of diffusion-derived fractional anisotropy (FA). The brain diffusion images were acquired from 10 healthy subjects using 30 and 60 directions encoding schemes based on the electrostatic repulsive forces. For the 30 directions encoding, 2 sets of diffusion images were acquired with the same parameters, except for the phase-encode blips which had opposing polarities along the anteroposterior direction. For the 60 directions encoding, non–diffusion-weighted and diffusion-weighted images were obtained with forward phase-encoding blips and non–diffusion-weighted images with the same parameter, except for the phase-encode blips, which had opposing polarities. FA images without and with distortion correction were compared in a voxel-wise manner with tract-based spatial statistics. We showed that images corrected with eddy and topup possessed higher FA values than images uncorrected and corrected with eddy_correct with trilinear (FSL default setting) or spline interpolation in most white matter skeletons, using both encoding schemes. Furthermore, the 60 directions encoding scheme was superior as measured by increased FA values to the 30 directions encoding scheme, despite comparable acquisition time. This study supports the combination of eddy and topup as a superior correction tool in diffusion imaging rather than the eddy_correct tool, especially with trilinear interpolation, using 60 directions encoding scheme.     

A human-specific allelic group of the MHC DRB1 gene in primates

Background

Diversity among human leukocyte antigen (HLA) molecules has been maintained by host-pathogen coevolution over a long period of time. Reflecting this diversity, the HLA loci are the most polymorphic in the human genome. One characteristic of HLA diversity is long-term persistence of allelic lineages, which causes trans-species polymorphisms to be shared among closely related species. Modern humans have disseminated across the world after their exodus from Africa, while chimpanzees have remained in Africa since the speciation event between humans and chimpanzees. It is thought that modern humans have recently acquired resistance to novel pathogens outside Africa. In the present study, we investigated HLA alleles that could contribute to this local adaptation in humans and also studied the contribution of natural selection to human evolution by using molecular data.

Results

Phylogenetic analysis of HLA-DRB1 genes identified two major groups, HLA Groups A and B. Group A formed a monophyletic clade distinct from DRB1 alleles in other Catarrhini, suggesting that Group A is a human-specific allelic group. Our estimates of divergence time suggested that seven HLA-DRB1 Group A allelic lineages in humans have been maintained since before the speciation event between humans and chimpanzees, while chimpanzees possess only one DRB1 allelic lineage (Patr-DRB1*03), which is a sister group to Group A. Experimental data showed that some Group A alleles bound to peptides derived from human-specific pathogens. Of the Group A alleles, three exist at high frequencies in several local populations outside Africa.

Conclusions

HLA Group A alleles are likely to have been retained in human lineages for a long period of time and have not expanded since the divergence of humans and chimpanzees. On the other hand, most orthologs of HLA Group A alleles may have been lost in the chimpanzee due to differences in selective pressures. The presence of alleles with high frequency outside of Africa suggests these HLA molecules result from the local adaptations of humans. Our study helps elucidate the mechanism by which the human adaptive immune system has coevolved with pathogens over a long period of time.     

Crystal Structure of Anti-polysialic Acid Antibody Single Chain Fv Fragment Complexed with Octasialic Acid: INSIGHT INTO THE BINDING PREFERENCE FOR POLYSIALIC ACID*

Polysialic acid is a linear homopolymer of α2–8-linked sialic acids attached mainly onto glycoproteins. Cell surface polysialic acid plays roles in cell adhesion and differentiation events in a manner that is often dependent on the degree of polymerization (DP). Anti-oligo/polysialic acid antibodies have DP-dependent antigenic specificity, and such antibodies are widely utilized in biological studies for detecting and distinguishing between different oligo/polysialic acids. A murine monoclonal antibody mAb735 has a unique preference for longer polymers of polysialic acid (DP >10), yet the mechanism of recognition at the atomic level remains unclear. Here, we report the crystal structure of mAb735 single chain variable fragment (scFv735) in complex with octasialic acid at 1.8 Å resolution. In the asymmetric unit, two scFv735 molecules associate with one octasialic acid. In both complexes of the unit, all the complementarity-determining regions except for L3 interact with three consecutive sialic acid residues out of the eight. A striking feature of the complex is that 11 ordered water molecules bridge the gap between antibody and ligand, whereas the direct antibody-ligand interaction is less extensive. The dihedral angles of the trisialic acid unit directly interacting with scFv735 are not uniform, indicating that mAb735 does not strictly favor the previously proposed helical conformation. Importantly, both reducing and nonreducing ends of the bound ligand are completely exposed to solvent. We suggest that mAb735 gains its apparent high affinity for a longer polysialic acid chain by recognizing every three sialic acid units in a paired manner.     

Analysis of Drosophila Glucuronyl C5-Epimerase: IMPLICATIONS FOR DEVELOPMENTAL ROLES OF HEPARAN SULFATE SULFATION COMPENSATION AND 2-O-SULFATED GLUCURONIC ACID*

During the biosynthesis of heparan sulfate (HS), glucuronyl C5-epimerase (Hsepi) catalyzes C5-epimerization of glucuronic acid (GlcA), converting it to iduronic acid (IdoA). Because HS 2-O-sulfotransferase (Hs2st) shows a strong substrate preference for IdoA over GlcA, C5-epimerization is required for normal HS sulfation. However, the physiological significance of C5-epimerization remains elusive. To understand the role of Hsepi in development, we isolated Drosophila Hsepi mutants. Homozygous mutants are viable and fertile with only minor morphological defects, including the formation of an ectopic crossvein in the wing, but they have a short lifespan. We propose that two mechanisms contribute to the mild phenotypes of Hsepi mutants: HS sulfation compensation and possible developmental roles of 2-O-sulfated GlcA (GlcA2S). HS disaccharide analysis showed that loss of Hsepi resulted in a significant impairment of 2-O-sulfation and induced compensatory increases in N- and 6-O-sulfation. Simultaneous block of Hsepi and HS 6-O-sulfotransferase (Hs6st) activity disrupted tracheoblast formation, a well established FGF-dependent process. This result suggests that the increase in 6-O-sulfation in Hsepi mutants is critical for the rescue of FGF signaling. We also found that the ectopic crossvein phenotype can be induced by expression of a mutant form of Hs2st with a strong substrate preference for GlcA-containing units, suggesting that this phenotype is associated with abnormal GlcA 2-O-sulfation. Finally, we show that Hsepi formed a complex with Hs2st and Hs6st in S2 cells, raising the possibility that this complex formation contributes to the close functional relationships between these enzymes.     

Serologic and molecular evidence for circulation of Crimean-Congo hemorrhagic fever virus in ticks and cattle in Zambia

Crimean-Congo hemorrhagic fever (CCHF) is a tick-borne zoonosis with a high case fatality rate in humans. Although the disease is widely found in Africa, Europe, and Asia, the distribution and genetic diversity of CCHF virus (CCHFV) are poorly understood in African countries. To assess the risks of CCHF in Zambia, where CCHF has never been reported, epidemiologic studies in cattle and ticks were conducted. Through an indirect immunofluorescence assay, CCHFV nucleoprotein-specific serum IgG was detected in 8.4% (88/1,047) of cattle. Among 290 Hyalomma ticks, the principal vector of CCHFV, the viral genome was detected in 11 ticks. Phylogenetic analyses of the CCHFV S and M genome segments revealed that one of the detected viruses was a genetic reassortant between African and Asian strains. This study provides compelling evidence for the presence of CCHFV in Zambia and its transmission to vertebrate hosts.     

HSPA1A is upregulated in periodontal ligament at early stage of tooth movement in rats

Heat shock proteins (HSPs) are molecular chaperones that maintain intracellular protein homeostasis and ensure survival of cells. Continuous orthodontic force on the tooth is considered to be a type of physical stress loaded to the periodontal ligament (PDL). However, little is known about the role of HSPs during tooth movement. This study was performed to examine the expression of HSPs in the PDL during tooth movement using laser microdissection, microarray analysis, real-time RT-PCR and immunohistochemistry. Gene expression of HSPA1A in the pressure zone of the PDL was higher during 6 h of tooth movement than in the control group. Expression of HSPA1A decreased with time. HSPA1A was also detected in the pressure zone of the PDL at the protein level 24 h after the initial tissue change. These results strongly suggest that expression of HSPA1A in the PDL during early stages of tooth movement is a critical factor for tissue reaction.     

High Affinity of the 5'-Upstream Region of a Winged Bean Chymotrypsin Inhibitor Gene for Nuclear Matrix

The 5'-upstream region of a winged bean chymotrypsin inhibitorgene (WCI-3b) was found to have a high affinity for nuclearmatrix. The region, named WCI-3b MAR (matrix attachment region),is highly A+T-rich and contains multiple sites interacting withnuclear matrix. A MAR was also found in the corresponding regionof the WCI-x gene, another active gene of the WCI family. SeveralMAR-binding proteins were detected in the wheat nuclear matrix. ⁴Present address: Friedrich Miescher Institute, P.O. Box 2543,CH-4002 Basel, Switzerland. ⁵Present address: Research Institute for Biological Sciences(RIBS), Kayo-cho, Jyobo-gun, Okayama, 716–1241 Japan.     

Canine X-linked muscular dystrophy in Japan (CXMDJ).

The purpose of this study was to develop a strain of canine X-linked muscular dystrophy (CXMD), a model of Duchenne muscular dystrophy, in Japan. A female beagle was artificially inseminated with frozen-thawed spermatozoa derived from an affected golden retriever. Subsequently, two carrier female dogs (G1 carriers) and four normal male littermates were produced. Thereafter, the two G1 carriers were mated with beagle sires. As a result, each bitch whelped three times, and out of 54 pups, 17 affected male descendants, and 11 carrier female descendants (G2 carriers) were detected. One G2 carrier was then mated with a beagle sire and 15 pups in two whelpings were produced, including five affected males and four carrier females (G3 carriers). A total of 10 female beagles were artificially inseminated to evaluate the fertility of the frozen-thawed spermatozoa from the two affected dogs. The whelping rates of the two affected dogs were 4/5 and the litter sizes were 5.0 +/- 1.41 and 6.0 +/- 0.82, respectively. These results indicate that a canine X-linked muscular dystrophy colony has been established in Japan. We called them CXMDJ.     

Role of the autonomic nervous system in the reduced maximal cardiac output at altitude.

After acclimatization to high altitude, maximal exercise cardiac output (QT) is reduced. Possible contributing factors include 1) blood volume depletion, 2) increased blood viscosity, 3) myocardial hypoxia, 4) altered autonomic nervous system (ANS) function affecting maximal heart rate (HR), and 5) reduced flow demand from reduced muscle work capability. We tested the role of the ANS reduction of HR in this phenomenon in five normal subjects by separately blocking the sympathetic and parasympathetic arms of the ANS during maximal exercise after 2-wk acclimatization at 3,800 m to alter maximal HR. We used intravenous doses of 8.0 mg of propranolol and 0.8 mg of glycopyrrolate, respectively. At altitude, peak HR was 170 +/- 6 beats/min, reduced from 186 +/- 3 beats/min (P = 0.012) at sea level. Propranolol further reduced peak HR to 139 +/- 2 beats/min (P = 0.001), whereas glycopyrrolate increased peak HR to sea level values, 184 +/- 3 beats/min, confirming adequate dosing with each drug. In contrast, peak O(2) consumption, work rate, and QT were similar at altitude under all drug treatments [peak QT = 16.2 +/- 1.2 (control), 15.5 +/- 1.3 (propranolol), and 16.2 +/- 1.1 l/min (glycopyrrolate)]. All QT results at altitude were lower than those at sea level (20.0 +/- 1.8 l/min in air). Therefore, this study suggests that, whereas the ANS may affect HR at altitude, peak QT is unaffected by ANS blockade. We conclude that the effect of altered ANS function on HR is not the cause of the reduced maximal QT at altitude.