MSPminer: abundance-based reconstitution of microbial pan-genomes from shotgun metagenomic data.

MOTIVATION: Analysis toolkits for shotgun metagenomic data achieve strain-level characterization of complex microbial communities by capturing intra-species gene content variation. Yet, these tools are hampered by the extent of reference genomes that are far from covering all microbial variability, as many species are still not sequenced or have only few strains available. Binning co-abundant genes obtained from de novo assembly is a powerful reference-free technique to discover and reconstitute gene repertoire of microbial species. While current methods accurately identify species core parts, they miss many accessory genes or split them into small gene groups that remain unassociated to core clusters. RESULTS: We introduce MSPminer, a computationally efficient software tool that reconstitutes Metagenomic Species Pan-genomes (MSPs) by binning co-abundant genes across metagenomic samples. MSPminer relies on a new robust measure of proportionality coupled with an empirical classifier to group and distinguish not only species core genes but accessory genes also. Applied to a large scale metagenomic dataset, MSPminer successfully delineates in a few hours the gene repertoires of 1661 microbial species with similar specificity and higher sensitivity than existing tools. The taxonomic annotation of MSPs reveals microorganisms hitherto unknown and brings coherence in the nomenclature of the species of the human gut microbiota. The provided MSPs can be readily used for taxonomic profiling and biomarkers discovery in human gut metagenomic samples. In addition, MSPminer can be applied on gene count tables from other ecosystems to perform similar analyses. AVAILABILITY AND IMPLEMENTATION: The binary is freely available for non-commercial users at www.enterome.com/downloads. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Gut microbiome alpha-diversity is not a marker of Parkinson's disease and multiple sclerosis.

The gut-brain axis may play a central role in the pathogenesis of neurological disorders. Dozens of case-control studies have been carried out to identify bacterial markers by the use of targeted metagenomics. Alterations of several taxonomic profiles have been confirmed across several populations, however, no consensus has been made regarding alpha-diversity. A recent publication has described and validated a novel method based on richness and evenness measures of the gut microbiome in order to reduce the complexity and multiplicity of alpha-diversity indices. We used these recently described richness and evenness composite measures to investigate the potential link between gut microbiome alpha-diversity and neurological disorders and to determine to what extent it could be used as a marker to diagnose neurological disorders from stool samples. We performed an exhaustive review of the literature to identify original published clinical studies including 16S rRNA gene sequencing on Parkinson's disease, multiple Sclerosis and Alzheimer's disease. Richness and evenness factors loadings were quantified from sequencing files in addition with the Shannon diversity index. For each disease, we performed a meta-analysis comparing the indices between patients and healthy controls. Seven studies were meta-analysed for Parkinson's disease, corresponding to 1067 subjects (631 Parkinson's Disease/436 healthy controls). Five studies were meta-analysed for multiple sclerosis, corresponding to 303 subjects (164 Multiple Sclerosis/139 healthy controls). For Alzheimer's disease, the meta-analysis was not done as only two studies matched our criteria. Neither richness nor evenness was significantly altered in Parkinson's disease and multiple sclerosis patients in comparison to healthy controls (P-value > 0.05). Shannon index was neither associated with neurological disorders (P-value > 0.05). After adjusting for age and sex, none of the alpha-diversity measures were associated with Parkinson's Disease. This is the first report investigating systematically alpha-diversity and its potential link to neurological disorders. Our study has demonstrated that unlike in other gastro-intestinal, immune and metabolic disorders, loss of bacterial diversity is not associated with Parkinson's disease and multiple sclerosis.

Impact of laparoscopic Roux-en-Y gastric bypass and sleeve gastrectomy on gut microbiota: a metagenomic comparative analysis.

BACKGROUND: Bariatric surgery is an effective therapeutic procedure for morbidly obese patients. The 2 most common interventions are sleeve gastrectomy (SG) and laparoscopic Roux-en-Y gastric bypass (LRYGB). OBJECTIVES: The aim of this study was to compare microbiome long-term microbiome after SG and LRYGB surgery in obese patients. SETTING: University Hospital, France; University Hospital, United States; and University Hospital, Switzerland. METHODS: Eighty-nine and 108 patients who underwent SG and LRYGB, respectively, were recruited. Stools were collected before and 6 months after surgery. Microbial DNA was analyzed with shotgun metagenomic sequencing (SOLiD 5500 xl Wildfire). MSPminer, a novel innovative tool to characterize new in silico biological entities, was used to identify 715 Metagenomic Species Pan-genome. One hundred forty-eight functional modules were analyzed using GOmixer and KEGG database. RESULTS: Both interventions resulted in a similar increase of Shannon's diversity index and gene richness of gut microbiota, in parallel with weight loss, but the changes of microbial composition were different. LRYGB led to higher relative abundance of aero-tolerant bacteria, such as Escherichia coli and buccal species, such as Streptococcus and Veillonella spp. In contrast, anaerobes, such as Clostridium, were more abundant after SG, suggesting better conservation of anaerobic conditions in the gut. Enrichment of Akkermansia muciniphila was also observed after both surgeries. Function-level changes included higher potential for bacterial use of supplements, such as vitamin B12, B1, and iron upon LRYGB. CONCLUSION: Microbiota changes after bariatric surgery depend on the nature of the intervention. LRYGB induces greater taxonomic and functional changes in gut microbiota than SG. Possible long-term health consequences of these alterations remain to be established.

An integrated in silico gene mapping strategy in inbred mice.

In recent years in silico analysis of common laboratory mice has been introduced and subsequently applied, in slightly different ways, as a methodology for gene mapping. Previously we have demonstrated some limitation of the methodology due to sporadic genetic correlations across the genome. Here, we revisit the three main aspects that affect in silico analysis. First, we report on the use of marker maps: we compared our existing 20,000 SNP map to the newly released 140,000 SNP map. Second, we investigated the effect of varying strain numbers on power to map QTL. Third, we introduced a novel statistical approach: a cladistic analysis, which is well suited for mouse genetics and has increased flexibility over existing in silico approaches. We have found that in our examples of complex traits, in silico analysis by itself does fail to uniquely identify quantitative trait gene (QTG)-containing regions. However, when combined with additional information, it may significantly help to prioritize candidate genes. We therefore recommend using an integrated work flow that uses other genomic information such as linkage regions, regions of shared ancestry, and gene expression information to obtain a list of candidate genes from the genome.

A genome-wide study of lupus: preliminary analysis and data release.

Systemic lupus erythematosus (SLE) is a complex autoimmune disease in which genetics is one of the underlying risk factors. Complicating the analysis of SLE is the fact that the phenotype is heterogeneous, with variations in clinical features, and subgroups are associated with a variety of different autoantibodies. Many association studies of candidate genes as well as linkage studies have generated significant results, highlighting the multigenetic component of the disease. Those findings need to be replicated by independent laboratories and many other genes still remain to be identified. Here, we present our findings on the first genome-wide association study performed in 105 Caucasian patients and controls using the Affymetrix NspI array. The SLE patients studied here were all characterized by the presence of autoantibodies against the U1-6 small nuclear ribonucleoprotein antigen spliceosomal complex. Testing each SNP individually for association with disease, we identified 7 markers with an uncorrected P-value < 0.0001. Of those, three were in reported regions of linkage, namely 20p12, 2q37, and 4p15. To increase the power of our study, we included 60 controls corresponding to the parents from the publicly available CEPH families. The analysis of the 165 cases and controls leads to 28 SNPs being associated with an uncorrected P-value < 0.0001.

Ultrafine mapping of SNPs from mouse strains C57BL/6J, DBA/2J, and C57BLKS/J for loci contributing to diabetes and atherosclerosis susceptibility.

The inbred mouse strain C57BLKS/J (BKS) carrying a mutation of the leptin receptor lepr(-/-) (BKS-db) is a classic mouse model of type 2 diabetes. While BKS was originally presumed to be a substrain of C57BL/6J (B6), it has become apparent that its genome contains introgressed regions from a DBA/2 (DBA)-like strain and perhaps other unidentified sources. It has been hypothesized that the strikingly enhanced diabetes susceptibility of BKS-db compared with B6-db is conferred by this introgressed DNA. Using high-density single nucleotide polymorphisms, we have mapped the DBA and other contaminating DNA regions present in BKS. Thus, approximately 70% of its genome appears to derive from B6, with approximately 20% from DBA and another 9% from an unidentified donor. Comparison with 56 diverse inbred strains suggests that this donor may be a less common inbred strain or an outbred or wild strain. Using expression data from a B6 x DBA cross, we identified differentially regulated genes between these two strains. Those cis-regulated genes located on DBA-like blocks in BKS constitute primary candidates for genes contributing to diabetes susceptibility in the BKS-db strain. To further prioritize these candidates, we identified those cis-acting expression quantitative trait loci whose expression significantly correlates with diabetes-related phenotypes.

Gut Microbial Diversity Is Reduced in Smokers with Crohn's Disease.

BACKGROUND: Smoking has a negative impact on Crohn's disease (CD), but the mechanisms underlying this association are unclear. We compared the gut microbiota composition of smoking with nonsmoking patients with CD using a metagenomic approach. METHODS: Stool samples and clinical data were collected from current smokers and nonsmokers with CD from France and the Netherlands, matched for country, gender, age, disease activity, and body mass index. Fecal DNA was sequenced on an Illumina HiSeq 2500. On average, 40 million paired-end reads were generated per sample. Gene richness and the Shannon index were computed to assess microbial diversity. Wilcoxon's signed-rank tests for paired samples were performed to detect differences between the 2 groups. RESULTS: In total, 21 smoking and 21 nonsmoking patients with CD were included. Compared with nonsmoking patients, gut microbial gene richness (P = 0.01), genus diversity (P < 0.01), and species diversity (P = 0.01) were decreased in smoking patients. This was accompanied by a reduced relative abundance of the genera Collinsella (P = 0.02), Enterorhabdus (P = 0.02), and Gordonibacter (P = 0.02) in smokers. No statistically significant differences at the species level were observed, although smokers had lower proportions of Faecalibacterium prausnitzii (P = 0.10). CONCLUSIONS: Gut microbial diversity is reduced in smokers with CD compared with nonsmokers with CD. The microbial profile differs between these groups at the genus level. Future studies should evaluate whether intestinal microbes mediate the adverse effects of smoking in CD.

Common sequence variation in FLNB regulates bone structure in women in the general population and FLNB mRNA expression in osteoblasts in vitro.

Previous data from our group indicate that BMD is linked to chromosome 3p14-p21. Because the filamin B (FLNB gene resides in this region, is the cause of skeletal dysplasias, and was identified among the top genes in our bioinformatics analysis, we hypothesized a role for FLNB in the regulation of bone structure in the general population. Using a tag single nucleotide polymorphism (SNP) approach, a family study of 767 female sibs in which the 3p14-p21 linkage with BMD was previously shown was examined. FLNB variants showing a BMD association were tested in two additional data sets, a study of 1085 UK female twins and a population study (CAIFOS) of 1315 Australian women. Genotype-expression studies were performed in 96 human osteoblast lines to examine the variants in vitro. rs7637505, rs9822918, rs2177153, and rs2001972 showed association with femoral neck (p = 0.0002-0.02) in the family-based study. The twin study provided further support for an association between rs7637505 and femoral neck and spine BMD (p = 0.02-0.03). The CAIFOS study further suggested an association between rs2177153 and rs9822918 and femoral neck BMD (p = 0.004-0.03). Prevalent fractures were increased in carriers of the A allele of rs2177153 (p = 0.009). In vitro studies showed association between rs11130605, itself in strong LD with rs7637505, and FLNB mRNA expression. These findings suggest common variants in FLNB have effects on bone structure in women. Although the location of variants having effects is not entirely consistent, variation at the 5' end of the gene may reflect effects on levels of FLNB transcription efficiency.

Glucocorticoid receptor polymorphisms and intraocular pressure response to intravitreal triamcinolone acetonide.

BACKGROUND: Elevation of intraocular pressure (IOP) following injection of intravitreal triamcinolone acetonide (IVTA) is an important clinical problem. The etiology of the steroid response is poorly understood, although a genetic determinant has long been suspected. We performed a pharmacogenomic association study with glucocorticoid receptor polymorphisms. MATERIALS AND METHODS: Fifty-two patients (56 eyes) who underwent treatment with IVTA for various retinal diseases were genotyped for six well-studied glucocorticoid receptor polymorphisms (ER22/23EK, N363S, BclI, N766N, and single nucleotide polymorphisms (SNPs) within introns 3 and 4). RESULTS: Three polymorphisms (ER22/23EK, N363S, and the intron 3 SNP) were essentially nonpolymorphic within this population sample and excluded from further analysis. The remaining three polymorphisms (BclI, N766N, and within intron 4) passed the Hardy-Weinberg Equilibrium test, indicating good genotyping quality and normal population distribution of allelic frequency. No statistically significant correlations were found between these three polymorphisms and magnitude of IOP elevation following IVTA, using single point association and haplotype analyses. CONCLUSIONS: In this small, pilot study, we found no statistically significant relationship between glucocorticoid receptor polymorphisms and IOP elevation following IVTA. The precise etiology of the steroid response remains obscure. To our knowledge, this is the first published pharmacogenomic study of this common clinical entity.

A comprehensive mouse IBD database for the efficient localization of quantitative trait loci.

Traditional fine-mapping approaches in mouse genetics that go from a linkage region to a candidate gene are very costly and time consuming. Shared ancestry regions, along with the combination of genetics and genomics approaches, provide a powerful tool to shorten the time and effort required to identify a causative gene. In this article we present a novel methodology that predicts IBD (identical by descent) regions between pairs of inbred strains using single nucleotide polymorphism (SNP) maps. We have validated this approach by comparing the IBD regions, estimated using different algorithms, to the results derived using the sequence information in the strains present in the Celera Mouse Database. We showed that based on the current publicly available SNP genotypes, large IBD regions (>1 Mb) can be identified successfully. By assembling a list of 21,514 SNPs in 61 common inbred strains, we inferred IBD regions between all pairs of strains and confirmed, for the first time, that existing quantitative trait genes (QTG) and susceptibility genes all lie outside of IBD regions. We also illustrated how knowledge of IBD structures can be applied to strain selection for future crosses. We have made our results available for data mining and download through a public website ( http://www.mouseibd.florida.scripps.edu ).

A genome-wide set of congenic mouse strains derived from DBA/2J on a C57BL/6J background.

In the analysis of complex traits, congenic strains are powerful tools because they allow characterization of a single locus in the absence of genetic variation throughout the remainder of the genome. Here, we report the construction and initial characterization of a genome-wide panel of congenic strains derived from the donor strain DBA/2J on the background strain C57BL/6J. For many strains, we have carried out high-density SNP genotyping to precisely map the congenic interval and to identify any contaminating regions. Certain strains exhibit striking variation in litter size and in the ratio of females to males. We illustrate the utility of the set by "Mendelizing" the complex trait of myocardial calcification. These 65 strains cover more than 95% of the autosomal genome and should facilitate the analysis of the many genetic trait differences that have been reported between these parental strains.

Fine mapping of a putative tuberculosis-susceptibility locus on chromosome 15q11-13 in African families.

Host genetics plays an important role in individual susceptibility and resistance to infectious diseases, but no genes have yet been identified using genome-wide screens. Twin studies have indicated that tuberculosis susceptibility has a significant host genetic component, and several genes appear to be involved. Recently, a genome-wide linkage analysis of 136 African families identified chromosome 15q11-13 as a region with suggestive evidence of linkage, with a LOD score of 2.0. We tested 10 microsatellite markers and 5 positional candidate genes in this chromosomal region for deviation from random transmission from parents to affected offspring. The polymorphisms, lying in a region of 14 cM, were initially typed in the same 79 Gambian families used in the genome screen. A borderline significant association with a 7 bp deletion in UBE3A (P = 0.01) was found. This polymorphism was then evaluated further in a larger series of families with tuberculosis, including 44 Guinean families and 57 families from South Africa. Testing for association between the deletion and tuberculosis across all the families using the exact symmetry test further supported the association (overall P = 0.002). These fine-mapping data suggest that UBE3A or a closely flanking gene may be a tuberculosis-susceptibility locus.