Shigeki Nakagome

Dysbiosis of Salivary Microbiota in Inflammatory Bowel Disease and Its Association With Oral Immunological Biomarkers

Analysis of microbiota in various biological and environmental samples under a variety of conditions has recently become more practical due to remarkable advances in next-generation sequencing. Changes leading to specific biological states including some of the more complex diseases can now be characterized with relative ease. It is known that gut microbiota is involved in the pathogenesis of inflammatory bowel disease (IBD), mainly Crohns disease and ulcerative colitis, exhibiting symptoms in the gastrointestinal tract. Recent studies also showed increased frequency of oral manifestations among IBD patients, indicating aberrations in the oral microbiota. Based on these observations, we analyzed the composition of salivary microbiota of 35 IBD patients by 454 pyrosequencing of the bacterial 16S rRNA gene and compared it with that of 24 healthy controls (HCs). The results showed that Bacteroidetes was significantly increased with a concurrent decrease in Proteobacteria in the salivary microbiota of IBD patients. The dominant genera, Streptococcus, Prevotella, Neisseria, Haemophilus, Veillonella, and Gemella, were found to largely contribute to dysbiosis (dysbacteriosis) observed in the salivary microbiota of IBD patients. Analysis of immunological biomarkers in the saliva of IBD patients showed elevated levels of many inflammatory cytokines and immunoglobulin A, and a lower lysozyme level. A strong correlation was shown between lysozyme and IL-1β levels and the relative abundance of Streptococcus, Prevotella, Haemophilus and Veillonella. Our data demonstrate that dysbiosis of salivary microbiota is associated with inflammatory responses in IBD patients, suggesting that it is possibly linked to dysbiosis of their gut microbiota.

Unequal rates of Y chromosome gene divergence during speciation of the family Ursidae.

Evolution of the bear family Ursidae is well investigated in terms of morphological, paleontological, and genetic features. However, several phylogenetic ambiguities occur within the subfamily Ursinae (the family Ursidae excluding the giant panda and spectacled bear), which may correlate with behavioral traits of female philopatry and male-biased dispersal which form the basis of the observed matriarchal population structure in these species. In the process of bear evolution, we investigate the premise that such behavioral traits may be reflected in patterns of variation among genes with different modes of inheritance: matrilineal mitochondrial DNA (mtDNA), patrilineal Y chromosome, biparentally inherited autosomes, and the X chromosome. In the present study, we sequenced 3 Y-linked genes (3,453 bp) and 4 X-linked genes (4,960 bp) and reanalyzed previously published sequences from autosome genes (2,347 bp) in ursid species to investigate differences in evolutionary rates associated with patterns of inheritance. The results describe topological incongruence between sex-linked genes and autosome genes and between nuclear DNA and mtDNA. In more ancestral branches within the bear phylogeny, Y-linked genes evolved faster than autosome and X-linked genes, consistent with expectations based on male-driven evolution. However, this pattern changes among branches leading to each species within the lineage of Ursinae whereby the evolutionary rates of Y-linked genes have fewer than expected substitutions. This inconsistency between more recent nodes of the bear phylogeny with more ancestral nodes may reflect the influences of sex-biased dispersal as well as molecular evolutionary characteristics of the Y chromosome, and stochastic events in species natural history, and phylogeography unique to ursine bears.

Population‐Specific Susceptibility to Crohn's Disease and Ulcerative Colitis; Dominant and Recessive Relative Risks in the Japanese Population

Crohns disease (CD), a type of chronic inflammatory bowel disease (IBD), is commonly found in European and East Asian countries. The calculated heritability of CD appears to be higher than that of ulcerative colitis (UC), another type of IBD. Recent genome‐wide association studies (GWAS) have identified more than thirty CD‐associated genes/regions in the European population. In the East Asian population, however, a clear association between CD and an associated gene has only been detected with TNFSF15. In order to determine if CD susceptibility differs geographically, nine SNPs from seven of the European CD‐associated genomic regions were selected for analysis. The genotype frequencies for these SNPs were compared among the 380 collected Japanese samples, which consisted of 212 IBD cases and 168 controls. We detected a significant association of both CD and UC with only the TNFSF15 gene. Analysis by the modified genotype relative risk test (mGRR) indicated that the risk allele of TNFSF15 is dominant for CD, but is recessive for UC. These results suggest that CD and UC susceptibility differs between the Japanese and European populations. Furthermore, it is also likely that CD and UC share a causative factor which exhibits a different dominant/recessive relative risk in the Japanese population.

A time transect of exomes from a Native American population before and after European contact

A major factor for the population decline of Native Americans after European contact has been attributed to infectious disease susceptibility. To investigate whether a pre-existing genetic component contributed to this phenomenon, here we analyse 50 exomes of a continuous population from the Northwest Coast of North America, dating from before and after European contact. We model the population collapse after European contact, inferring a 57% reduction in effective population size. We also identify signatures of positive selection on immune-related genes in the ancient but not the modern group, with the strongest signal deriving from the human leucocyte antigen (HLA) gene HLA-DQA1. The modern individuals show a marked frequency decrease in the same alleles, likely due to the environmental change associated with European colonization, whereby negative selection may have acted on the same gene after contact. The evident shift in selection pressures correlates to the regional European-borne epidemics of the 1800s.

Genome-wide SNP analysis reveals population structure and demographic history of the Ryukyu Islanders in the southern part of the Japanese Archipelago

The Ryukyu Islands are located to the southwest of the Japanese archipelago. Archaeological evidence has revealed the existence of prehistoric cultural differentiation between the northern Ryukyu islands of Amami and Okinawa, and the southern Ryukyu islands of Miyako and Yaeyama. To examine a genetic subdivision in the Ryukyu Islands, we conducted genome-wide single nucleotide polymorphism typing of inhabitants from the Okinawa Islands, the Miyako Islands, and the Yaeyama Islands. Principal component and cluster analyses revealed genetic differentiation among the island groups, especially between Okinawa and Miyako. No genetic affinity was observed between aboriginal Taiwanese and any of the Ryukyu populations. The genetic differentiation observed between the inhabitants of the Okinawa Islands and the Miyako Islands is likely to have arisen due to genetic drift rather than admixture with people from neighboring regions. Based on the observed genetic differences, the divergence time between the inhabitants of Okinawa and Miyako islands was dated to the Holocene. These findings suggest that the Pleistocene inhabitants, whose bones have been found on the southern Ryukyu Islands, did not make a major genetic contribution, if any, to the present-day inhabitants of the southern Ryukyu Islands.

Crohn’s disease risk alleles on the NOD2 locus have been maintained by natural selection on standing variation

Risk alleles for complex diseases are widely spread throughout human populations. However, little is known about the geographic distribution and frequencies of risk alleles, which may contribute to differences in disease susceptibility and prevalence among populations. Here, we focus on Crohns disease (CD) as a model for the evolutionary study of complex disease alleles. Recent genome-wide association studies and classical linkage analyses have identified more than 70 susceptible genomic regions for CD in Europeans, but only a few have been confirmed in non-European populations. Our analysis of eight European-specific susceptibility genes using HapMap data shows that at the NOD2 locus the CD-risk alleles are linked with a haplotype specific to CEU at a frequency that is significantly higher compared with the entire genome. We subsequently examined nine global populations and found that the CD-risk alleles spread through hitchhiking with a high-frequency haplotype (H1) exclusive to Europeans. To examine the neutrality of NOD2, we performed phylogenetic network analyses, coalescent simulation, protein structural prediction, characterization of mutation patterns, and estimations of population growth and time to most recent common ancestor (TMRCA). We found that while H1 was significantly prevalent in European populations, the H1 TMRCA predated human migration out of Africa. H1 is likely to have undergone negative selection because 1) the root of H1 genealogy is defined by a preexisting amino acid substitution that causes serious conformational changes to the NOD2 protein, 2) the haplotype has almost become extinct in Africa, and 3) the haplotype has not been affected by the recent European expansion reflected in the other haplotypes. Nevertheless, H1 has survived in European populations, suggesting that the haplotype is advantageous to this group. We propose that several CD-risk alleles, which destabilize and disrupt the NOD2 protein, have been maintained by natural selection on standing variation because the deleterious haplotype of NOD2 is advantageous in diploid individuals due to heterozygote advantage and/or intergenic interactions.

Comment on “Nuclear Genomic Sequences Reveal that Polar Bears Are an Old and Distinct Bear Lineage”

Based on nuclear and mitochondrial DNA, Hailer et al. (Reports, 20 April 2012, p. 344) suggested early divergence of polar bears from a common ancestor with brown bears and subsequent introgression. Our population genetic analysis that traces each of the genealogies in the independent nuclear loci does not support the evolutionary model proposed by the authors.

Estimating the ages of selection signals from different epochs in human history

Genetic variation harbors signatures of natural selection driven by selective pressures that are often unknown. Estimating the ages of selection signals may allow reconstructing the history of environmental changes that shaped human phenotypes and diseases. We have developed an approximate Bayesian computation (ABC) approach to estimate allele ages under a model of selection on new mutations and under demographic models appropriate for human populations. We have applied it to two resequencing data sets: An ultra-high depth data set from a relatively small sample of unrelated individuals and a lower depth data set in a larger sample with transmission information. In addition to evaluating the accuracy of our method based on simulations, for each SNP, we assessed the consistency between the posterior probabilities estimated by the ABC approach and the ancient DNA record, finding good agreement between the two types of data and methods. Applying this ABC approach to data for eight single nucleotide polymorphisms (SNPs), we were able to rule out an onset of selection prior to the dispersal out-of-Africa for three of them and more recent than the spread of agriculture for an additional three SNPs.

Finding the Factors of Reduced Genetic Diversity on X Chromosomes of Macaca fascicularis : Male-Driven Evolution, Demography, and Natural Selection

The ratio of genetic diversity on X chromosomes relative to autosomes in organisms with XX/XY sex chromosomes could provide fundamental insight into the process of genome evolution. Here we report this ratio for 24 cynomolgus monkeys (Macaca fascicularis) originating in Indonesia, Malaysia, and the Philippines. The average X/A diversity ratios in these samples was 0.34 and 0.20 in the Indonesian–Malaysian and Philippine populations, respectively, considerably lower than the null expectation of 0.75. A Philippine population supposed to derive from an ancestral population by founding events showed a significantly lower ratio than the parental population, suggesting a demographic effect for the reduction. Taking sex-specific mutation rate bias and demographic effect into account, expected X/A diversity ratios generated by computer simulations roughly agreed with the observed data in the intergenic regions. In contrast, silent sites in genic regions on X chromosomes showed strong reduction in genetic diversity and the observed X/A diversity ratio in the genic regions cannot be explained by mutation rate bias and demography, indicating that natural selection also reduces the level of polymorphism near genes. Whole-genome analysis of a female cynomolgus monkey also supported the notion of stronger reduction of genetic diversity near genes on the X chromosome.

Colonic transcriptional response to 1α,25(OH)2 vitamin D3 in African- and European-Americans

Colorectal cancer (CRC) is a significant health burden especially among African Americans (AA). Epidemiological studies have correlated low serum vitamin D with CRC risk, and, while hypovitaminosis D is more common and more severe in AA, the mechanisms by which vitamin D modulates CRC risk and how these differ by race are not well understood. Active vitamin D (1α,25(OH)2D3) has chemoprotective effects primarily through transcriptional regulation of target genes in the colon. We hypothesized that transcriptional response to 1α,25(OH)2D3 differs between AA and European Americans (EA) irrespective of serum vitamin D and that regulatory variants could impact transcriptional response. We treated ex vivo colon cultures from 34 healthy subjects (16 AA and 18 EA) with 0.1μM 1α,25(OH)2D3 or vehicle control for 6h and performed genome-wide transcriptional profiling. We found 8 genes with significant differences in transcriptional response to 1α,25(OH)2D3 between AA and EA with definitive replication of inter-ethnic differences for uridine phosphorylase 1 (UPP1) and zinc finger-SWIM containing 4 (ZSWIM4). We performed expression quantitative trait loci (eQTL) mapping and identified response cis-eQTLs for ZSWIM4 as well as for histone deacetylase 3 (HDAC3), the latter of which showed a trend toward significant inter-ethnic differences in transcriptional response. Allele frequency differences of eQTLs for ZSWIM4 and HDAC3 accounted for observed transcriptional differences between populations. Taken together, our results demonstrate that transcriptional response to 1α,25(OH)2D3 differs between AA and EA independent of serum 25(OH)D levels. We provide evidence in support of a genetic regulatory mechanism underlying transcriptional differences between populations for ZSWIM4 and HDAC3. Further work is needed to elucidate how response eQTLs modify vitamin D response and whether genotype and/or transcriptional response correlate with chemopreventive effects. Relevant biomarkers, such as tissue-specific 1α,25(OH)2D3 transcriptional response, could identify individuals likely to benefit from vitamin D for CRC prevention as well as elucidate basic mechanisms underlying CRC disparities.