A global metagenomic map of urban microbiomes and antimicrobial resistance.

We present a global atlas of 4,728 metagenomic samples from mass-transit systems in 60 cities over 3 years, representing the first systematic, worldwide catalog of the urban microbial ecosystem. This atlas provides an annotated, geospatial profile of microbial strains, functional characteristics, antimicrobial resistance (AMR) markers, and genetic elements, including 10,928 viruses, 1,302 bacteria, 2 archaea, and 838,532 CRISPR arrays not found in reference databases. We identified 4,246 known species of urban microorganisms and a consistent set of 31 species found in 97% of samples that were distinct from human commensal organisms. Profiles of AMR genes varied widely in type and density across cities. Cities showed distinct microbial taxonomic signatures that were driven by climate and geographic differences. These results constitute a high-resolution global metagenomic atlas that enables discovery of organisms and genes, highlights potential public health and forensic applications, and provides a culture-independent view of AMR burden in cities.

Ectopic gut colonization: a metagenomic study of the oral and gut microbiome in Crohn's disease.

BACKGROUND: This study aims to characterize, the gut and oral microbiome in Asian subjects with Crohn's disease (CD) using whole genome shotgun sequencing, thereby allowing for strain-level comparison. METHODS: A case-control study with age, sex and ethnicity matched healthy controls was conducted. CD subjects were limited to well-controlled patients without oral manifestations. Fecal and saliva samples were collected for characterization of gut and oral microbiome respectively. Microbial DNA were extracted, libraries prepared and sequenced reads profiled. Taxonomic diversity, taxonomic association, strain typing and microbial gene pathway analyses were conducted. RESULTS: The study recruited 25 subjects with CD and 25 healthy controls. The oral microbe Streptococcus salivarius was found to be enriched and of concordant strains in the gut and oral microbiome of Crohn's disease subjects. This was more likely in CD subjects with higher Crohn's Disease Activity Index (184.3 ± 2.9 vs 67.1 ± 82.5, p = 0.012) and active disease status (Diarrhoea/abdominal pain/blood-in-stool/fever and fatigue) (p = 0.016). Gut species found to be significantly depleted in CD compared to control (Relative abundance: Median[Range]) include: Faecalibacterium prausnitzii (0.03[0.00-4.56] vs 13.69[5.32-18.71], p = 0.010), Roseburia inulinivorans (0.00[0.00-0.03] vs 0.21[0.01-0.53], p = 0.010) and Alistipes senegalensis (0.00[0.00-0.00] vs 0.00[0.00-0.02], p = 0.029). While Clostridium nexile (0.00[0.00-0.12] vs 0.00[0.00-0.00], p = 0.038) and Ruminococcus gnavus (0.43[0.02-0.33] vs 0.00[0.00-0.13], p = 0.043) were found to be enriched. C. nexile enrichment was not found in CD subjects of European descent. Microbial arginine (Linear-discriminant-analysis: 3.162, p = 0.001) and isoprene (Linear-discriminant-analysis: 3.058, p < 0.001) pathways were found at a higher relative abundance level in gut microbiome of Crohn's disease. CONCLUSIONS: There was evidence of ectopic gut colonization by oral bacteria, especially during the active phase of CD. Previously studied gut microbial differences were detected, in addition to novel associations which could have resulted from geographical/ethnic differences to subjects of European descent. Differences in microbial pathways provide possible targets for microbiome modification.

Common Genetic Variations Associated with the Persistence of Immunity following Childhood Immunization.

Vaccines have revolutionized public health, preventing millions of deaths each year, particularly in childhood. Yet, there is considerable variability in the magnitude and persistence of vaccine-induced immunity. Maintenance of specific antibody is essential for continuity of vaccine-induced serological protection. We conducted a genome-wide association study into the persistence of immunity to three childhood vaccines: capsular group C meningococcal (MenC), Haemophilus influenzae type b, and tetanus toxoid (TT) vaccines. We detail associations between variants in a locus containing a family of signal-regulatory proteins and the persistence MenC immunity. We postulate a regulatory role for the lead SNP, with supporting epigenetic and expression quantitative trait loci data. Furthermore, we define associations between SNPs in the human leukocyte antigen (HLA) locus and the persistence of TT-specific immunity. Moreover, we describe four classical HLA alleles, HLA DRB1∗0301, HLA DQB1∗0201, HLA DQB1∗0602, and HLA DRB1∗1501, associated with TT-specific immunity, independent of the lead SNP association.

Identification of regulatory variants associated with genetic susceptibility to meningococcal disease.

Non-coding genetic variants play an important role in driving susceptibility to complex diseases but their characterization remains challenging. Here, we employed a novel approach to interrogate the genetic risk of such polymorphisms in a more systematic way by targeting specific regulatory regions relevant for the phenotype studied. We applied this method to meningococcal disease susceptibility, using the DNA binding pattern of RELA - a NF-kB subunit, master regulator of the response to infection - under bacterial stimuli in nasopharyngeal epithelial cells. We designed a custom panel to cover these RELA binding sites and used it for targeted sequencing in cases and controls. Variant calling and association analysis were performed followed by validation of candidate polymorphisms by genotyping in three independent cohorts. We identified two new polymorphisms, rs4823231 and rs11913168, showing signs of association with meningococcal disease susceptibility. In addition, using our genomic data as well as publicly available resources, we found evidences for these SNPs to have potential regulatory effects on ATXN10 and LIF genes respectively. The variants and related candidate genes are relevant for infectious diseases and may have important contribution for meningococcal disease pathology. Finally, we described a novel genetic association approach that could be applied to other phenotypes.

Extensive Ethnic Variation and Linkage Disequilibrium at the FCGR2/3 Locus: Different Genetic Associations Revealed in Kawasaki Disease.

The human Fc-gamma receptors (FcγRs) link adaptive and innate immunity by binding immunoglobulin G (IgG). All human low-affinity FcγRs are encoded by the FCGR2/3 locus containing functional single nucleotide polymorphisms (SNPs) and gene copy number variants. This locus is notoriously difficult to genotype and high-throughput methods commonly used focus on only a few SNPs. We performed multiplex ligation-dependent probe amplification for all relevant genetic variations at the FCGR2/3 locus in >4,000 individuals to define linkage disequilibrium (LD) and allele frequencies in different populations. Strong LD and extensive ethnic variation in allele frequencies was found across the locus. LD was strongest for the FCGR2C-ORF haplotype (rs759550223+rs76277413), which leads to expression of FcγRIIc. In Europeans, the FCGR2C-ORF haplotype showed strong LD with, among others, rs201218628 (FCGR2A-Q27W, r2 = 0.63). LD between these two variants was weaker (r2 = 0.17) in Africans, whereas the FCGR2C-ORF haplotype was nearly absent in Asians (minor allele frequency <0.005%). The FCGR2C-ORF haplotype and rs1801274 (FCGR2A-H131R) were in weak LD (r2 = 0.08) in Europeans. We evaluated the importance of ethnic variation and LD in Kawasaki Disease (KD), an acute vasculitis in children with increased incidence in Asians. An association of rs1801274 with KD was previously shown in ethnically diverse genome-wide association studies. Now, we show in 1,028 European KD patients that the FCGR2C-ORF haplotype, although nearly absent in Asians, was more strongly associated with susceptibility to KD than rs1801274 in Europeans. Our data illustrate the importance of interpreting findings of association studies concerning the FCGR2/3 locus with knowledge of LD and ethnic variation.

Comparison of self-reported and register-based hospital medical data on comorbidities in women.

Breast cancer patients commonly present with comorbidities which are known to influence treatment decisions and survival. We aim to examine agreement between self-reported and register-based medical records (National Patient Register [NPR]). Ascertainment of nine conditions, using individually-linked data from 64,961 women enrolled in the Swedish KARolinska MAmmography Project for Risk Prediction of Breast Cancer (KARMA) study. Agreement was assessed using observed proportion of agreement (overall agreement), expected proportion of agreement, and Cohen's Kappa statistic. Two-stage logistic regression models taking into account chance agreement were used to identify potential predictors of overall agreement. High levels of overall agreement (i.e. ≥86.6%) were observed for all conditions. Substantial agreement (Cohen's Kappa) was observed for myocardial infarction (0.74), diabetes (0.71) and stroke (0.64) between self-reported and NPR data. Moderate agreement was observed for preeclampsia (0.51) and hypertension (0.46). Fair agreement was observed for heart failure (0.40) and polycystic ovaries or ovarian cysts (0.27). For hyperlipidemia (0.14) and angina (0.10), slight agreement was observed. In most subgroups we observed negative specific agreement of >90%. There is no clear reference data source for ascertainment of conditions. Negative specific agreement between NPR and self-reported data is consistently high across all conditions.

Bifid T waves on the ECG and genetic variation in calcium channel voltage-dependent beta 2 subunit gene (CACNB2) in acute Kawasaki disease.

BACKGROUND: We previously described the association of genetic variants in calcium channel genes and susceptibility to Kawasaki disease (KD), an acute, self-limited vasculitis, and the most common cause of acquired cardiac disease in children. Abnormal repolarization of cardiomyocytes and changes in T wave morphology have been reported in KD but have not been studied systematically. METHODS: We analyzed acute and convalescent ECG T wave morphology in two independent cohorts of KD subjects and studied the association between bifid T waves and genetic variants in previously reported genes with SNVs associated with cardiac repolarization. RESULTS: Bifid T waves in limb leads were identified in 24% and 27% of two independent cohorts of acute KD subjects. Calcium channel voltage-dependent beta 2 subunit gene (CACNB2) (rs1409207) showed association with bifid T waves in both cohorts (nominal P = .04 and P = .0003, respectively). This CACNB2 polymorphism also showed association with KD susceptibility in a previously published KD genome wide association study data (nominal P = .009). CONCLUSION: This genotype/phenotype association study uncovered a variant in CACNB2 that may be associated with both KD susceptibility and bifid T waves, a novel signature of altered myocardial repolarization. The present study combined with published reports suggests that genetic variants in calcium channels and intracellular calcium signaling play a prominent role in shaping susceptibility to KD.

PRINCESS: Privacy-protecting Rare disease International Network Collaboration via Encryption through Software guard extensionS.

Motivation: We introduce PRINCESS, a privacy-preserving international collaboration framework for analyzing rare disease genetic data that are distributed across different continents. PRINCESS leverages Software Guard Extensions (SGX) and hardware for trustworthy computation. Unlike a traditional international collaboration model, where individual-level patient DNA are physically centralized at a single site, PRINCESS performs a secure and distributed computation over encrypted data, fulfilling institutional policies and regulations for protected health information. Results: To demonstrate PRINCESS' performance and feasibility, we conducted a family-based allelic association study for Kawasaki Disease, with data hosted in three different continents. The experimental results show that PRINCESS provides secure and accurate analyses much faster than alternative solutions, such as homomorphic encryption and garbled circuits (over 40 000× faster). Availability and Implementation: https://github.com/achenfengb/PRINCESS_opensource. Contact: shw070@ucsd.edu. Supplementary information: Supplementary data are available at Bioinformatics online.

Natural resistance to Meningococcal Disease related to CFH loci: Meta-analysis of genome-wide association studies.

Meningococcal disease (MD) remains an important infectious cause of life threatening infection in both industrialized and resource poor countries. Genetic factors influence both occurrence and severity of presentation, but the genes responsible are largely unknown. We performed a genome-wide association study (GWAS) examining 5,440,063 SNPs in 422 Spanish MD patients and 910 controls. We then performed a meta-analysis of the Spanish GWAS with GWAS data from the United Kingdom (combined cohorts: 897 cases and 5,613 controls; 4,898,259 SNPs). The meta-analysis identified strong evidence of association (P-value ≤ 5 × 10-8) in 20 variants located at the CFH gene. SNP rs193053835 showed the most significant protective effect (Odds Ratio (OR) = 0.62, 95% confidence interval (C.I.) = 0.52-0.73; P-value = 9.62 × 10-9). Five other variants had been previously reported to be associated with susceptibility to MD, including the missense SNP rs1065489 (OR = 0.64, 95% C.I.) = 0.55-0.76, P-value = 3.25 × 10-8). Theoretical predictions point to a functional effect of rs1065489, which may be directly responsible for protection against MD. Our study confirms the association of CFH with susceptibility to MD and strengthens the importance of this link in understanding pathogenesis of the disease.

Sensitive detection of chromatin-altering polymorphisms reveals autoimmune disease mechanisms.

Most disease associations detected by genome-wide association studies (GWAS) lie outside coding genes, but very few have been mapped to causal regulatory variants. Here, we present a method for detecting regulatory quantitative trait loci (QTLs) that does not require genotyping or whole-genome sequencing. The method combines deep, long-read chromatin immunoprecipitation-sequencing (ChIP-seq) with a statistical test that simultaneously scores peak height correlation and allelic imbalance: the genotype-independent signal correlation and imbalance (G-SCI) test. We performed histone acetylation ChIP-seq on 57 human lymphoblastoid cell lines and used the resulting reads to call 500,066 single-nucleotide polymorphisms de novo within regulatory elements. The G-SCI test annotated 8,764 of these as histone acetylation QTLs (haQTLs)­an order of magnitude larger than the set of candidates detected by expression QTL analysis. Lymphoblastoid haQTLs were highly predictive of autoimmune disease mechanisms. Thus, our method facilitates large-scale regulatory variant detection in any moderately sized cohort for which functional profiling data can be generated, thereby simplifying identification of causal variants within GWAS loci.

Polymorphisms in SP110 are not associated with pulmonary tuberculosis in Indonesians.

Despite being high transmissible, Mycobacterium tuberculosis (M. tuberculosis) infection causes active disease in only 5-10% of disease-susceptible individuals. This has instigated interest in studying potentially underlying genetic host factors and mechanisms in tuberculosis (TB). The recent identification of the Intracellular pathogen resistance 1 (Ipr1) gene, which plays a major role in controlling M. tuberculosis susceptibility and infection severity in mice (Pan et al., 2005), has prompted studies on its human homolog; SP110 in humans. Association of SP110 SNPs with pulmonary TB were first reported in a study on West African families (Tosh et al., 2006). Subsequent attempts to replicate these findings in other populations, including another West African (Ghanaian) cohort (Thye et al., 2006), however, were unsuccessful. Here we have genotyped 20 SNPs located in the SP110 gene, including the previously TB associated variants; rs2114592 and rs3948464, for the first time in a South East Asian cohort from Indonesia. Our study did not reveal any statistically significant associations between SP110 SNPs and pulmonary TB. In addition, a meta-analysis of the two previously TB associated SNPs revealed that these are not associated with TB, further confirming the lack of convincing evidence for SP110 to be implicated in TB susceptibility, as yet in humans.

A genome wide association study of pulmonary tuberculosis susceptibility in Indonesians.

BACKGROUND: There is reason to expect strong genetic influences on the risk of developing active pulmonary tuberculosis (TB) among latently infected individuals. Many of the genome wide linkage and association studies (GWAS) to date have been conducted on African populations. In order to identify additional targets in genetically dissimilar populations, and to enhance our understanding of this disease, we performed a multi-stage GWAS in a Southeast Asian cohort from Indonesia. METHODS: In stage 1, we used the Affymetrix 100 K SNP GeneChip marker set to genotype 259 Indonesian samples. After quality control filtering, 108 cases and 115 controls were analyzed for association of 95,207 SNPs. In stage 2, we attempted validation of 2,453 SNPs with promising associations from the first stage, in 1,189 individuals from the same Indonesian cohort, and finally in stage 3 we selected 251 SNPs from this stage to test TB association in an independent Caucasian cohort (n = 3,760) from Russia. RESULTS: Our study suggests evidence of association (P = 0.0004-0.0067) for 8 independent loci (nominal significance P < 0.05), which are located within or near the following genes involved in immune signaling: JAG1, DYNLRB2, EBF1, TMEFF2, CCL17, HAUS6, PENK and TXNDC4. CONCLUSIONS: Mechanisms of immune defense suggested by some of the identified genes exhibit biological plausibility and may suggest novel pathways involved in the host containment of infection with TB.

A genome-wide association study of hepatitis B vaccine response in an Indonesian population reveals multiple independent risk variants in the HLA region.

We performed a two-stage genome-wide association study (GWAS) of antibody titer in 3614 hepatitis B vaccine recipients from Indonesia's Riau Archipelago, leading to the identification of at least three independent signals within the human leukocyte antigen (HLA) complex. These appear to implicate HLA-DR [rs3135363; P= 6.53 × 10(-22); odds ratio (OR) = 1.53, 95% confidence interval (CI) = 1.35-1.74]; HLA-DP, previously associated with the risk of chronic hepatitis B infection (rs9277535; P= 2.91 × 10(-12); OR = 0.72, 95% CI = 0.63-0.81); and a gene rich HLA Class III interval (rs9267665; P = 1.24 × 10(-17); OR = 2.05, CI = 1.64-2.57). The substantial overlap of these variants and those identified by GWAS of chronic hepatitis B infection confirms vaccine response as a model for infection, while suggesting that the vaccine is least effective in those most at risk of lifelong infection, following exposure to the virus.