Functional metagenomic analysis of the human gut microbiome to identify novel salt tolerance genes

The ability to adapt to and respond to increases in external osmolarity is an important characteristic that enables bacteria to survive and proliferate in different environmental niches. When challenged with increased osmolarity, due to sodium chloride (NaCl) for example, bacteria elicit a phased response; firstly via uptake of potassium (K+), which is known as the primary response. This primary response is followed by the secondary response which is characterised by the synthesis or uptake of compatible solutes (osmoprotectants). The overall osmotic stress response is much broader however, involving many diverse cellular systems and processes. These ancillary mechanisms are arguably more interesting and give a more complete view of the osmotic stress response. The aim of this thesis was to identify novel genetic loci from the human gut microbiota that confer increased tolerance to osmotic stress using a functional metagenomic approach. Functional metagenomics is a powerful tool that enables the identification of novel genes from as yet uncultured bacteria from diverse environments through cloning, heterologous expression and phenotypic identification of a desired trait. Functional metagenomics does not rely on any previous sequence information to known genes and can therefore enable the discovery of completely novel genes and assign functions to new or known genes. Using a functional metagenomic approach, we have assigned a novel function to previously annotated genes; murB, mazG and galE, as well as a putative brp/blh family beta-carotene 15,15’-monooxygenase. Finally, we report the identification of a completely novel salt tolerance determinant with no current known homologues in the databases. Overall the genes identified originate from diverse taxonomic and phylogenetic groups commonly found in the human gastrointestinal (GI) tract, such as Collinsella and Eggerthella, Akkermansia and Bacteroides from the phyla Actinobacteria, Verrucomicrobia and Bacteroidetes, respectively. In addition, a number of the genes appear to have been acquired via lateral gene transfer and/or encoded on a prophage. To our knowledge, this thesis represents the first investigation to identify novel genes from the human gut microbiota involved in the bacterial osmotic stress response.

A genome-wide association study identifies susceptibility variants for type 2 diabetes in Han Chinese.

To investigate the underlying mechanisms of T2D pathogenesis, we looked for diabetes susceptibility genes that increase the risk of type 2 diabetes (T2D) in a Han Chinese population. A two-stage genome-wide association (GWA) study was conducted, in which 995 patients and 894 controls were genotyped using the Illumina HumanHap550-Duo BeadChip for the first genome scan stage. This was further replicated in 1,803 patients and 1,473 controls in stage 2. We found two loci not previously associated with diabetes susceptibility in and around the genes protein tyrosine phosphatase receptor type D (PTPRD) (P = 8.54x10(-10); odds ratio [OR] = 1.57; 95% confidence interval [CI] = 1.36-1.82), and serine racemase (SRR) (P = 3.06x10(-9); OR = 1.28; 95% CI = 1.18-1.39). We also confirmed that variants in KCNQ1 were associated with T2D risk, with the strongest signal at rs2237895 (P = 9.65x10(-10); OR = 1.29, 95% CI = 1.19-1.40). By identifying two novel genetic susceptibility loci in a Han Chinese population and confirming the involvement of KCNQ1, which was previously reported to be associated with T2D in Japanese and European descent populations, our results may lead to a better understanding of differences in the molecular pathogenesis of T2D among various populations.

Evolution of the sex-related locus and genomic features shared in microsporidia and fungi.

BACKGROUND: Microsporidia are obligate intracellular, eukaryotic pathogens that infect a wide range of animals from nematodes to humans, and in some cases, protists. The preponderance of evidence as to the origin of the microsporidia reveals a close relationship with the fungi, either within the kingdom or as a sister group to it. Recent phylogenetic studies and gene order analysis suggest that microsporidia share a particularly close evolutionary relationship with the zygomycetes. METHODOLOGY/PRINCIPAL FINDINGS: Here we expanded this analysis and also examined a putative sex-locus for variability between microsporidian populations. Whole genome inspection reveals a unique syntenic gene pair (RPS9-RPL21) present in the vast majority of fungi and the microsporidians but not in other eukaryotic lineages. Two other unique gene fusions (glutamyl-prolyl tRNA synthetase and ubiquitin-ribosomal subunit S30) that are present in metazoans, choanoflagellates, and filasterean opisthokonts are unfused in the fungi and microsporidians. One locus previously found to be conserved in many microsporidian genomes is similar to the sex locus of zygomycetes in gene order and architecture. Both sex-related and sex loci harbor TPT, HMG, and RNA helicase genes forming a syntenic gene cluster. We sequenced and analyzed the sex-related locus in 11 different Encephalitozoon cuniculi isolates and the sibling species E. intestinalis (3 isolates) and E. hellem (1 isolate). There was no evidence for an idiomorphic sex-related locus in this Encephalitozoon species sample. According to sequence-based phylogenetic analyses, the TPT and RNA helicase genes flanking the HMG genes are paralogous rather than orthologous between zygomycetes and microsporidians. CONCLUSION/SIGNIFICANCE: The unique genomic hallmarks between microsporidia and fungi are independent of sequence based phylogenetic comparisons and further contribute to define the borders of the fungal kingdom and support the classification of microsporidia as unusual derived fungi. And the sex/sex-related loci appear to have been subject to frequent gene conversion and translocations in microsporidia and zygomycetes.

The evolutionary history of ferns inferred from 25 low-copy nuclear genes.

UNLABELLED: • PREMISE OF THE STUDY: Understanding fern (monilophyte) phylogeny and its evolutionary timescale is critical for broad investigations of the evolution of land plants, and for providing the point of comparison necessary for studying the evolution of the fern sister group, seed plants. Molecular phylogenetic investigations have revolutionized our understanding of fern phylogeny, however, to date, these studies have relied almost exclusively on plastid data.• METHODS: Here we take a curated phylogenomics approach to infer the first broad fern phylogeny from multiple nuclear loci, by combining broad taxon sampling (73 ferns and 12 outgroup species) with focused character sampling (25 loci comprising 35877 bp), along with rigorous alignment, orthology inference and model selection.• KEY RESULTS: Our phylogeny corroborates some earlier inferences and provides novel insights; in particular, we find strong support for Equisetales as sister to the rest of ferns, Marattiales as sister to leptosporangiate ferns, and Dennstaedtiaceae as sister to the eupolypods. Our divergence-time analyses reveal that divergences among the extant fern orders all occurred prior to ∼200 MYA. Finally, our species-tree inferences are congruent with analyses of concatenated data, but generally with lower support. Those cases where species-tree support values are higher than expected involve relationships that have been supported by smaller plastid datasets, suggesting that deep coalescence may be reducing support from the concatenated nuclear data.• CONCLUSIONS: Our study demonstrates the utility of a curated phylogenomics approach to inferring fern phylogeny, and highlights the need to consider underlying data characteristics, along with data quantity, in phylogenetic studies.

ENIGMA and the individual: Predicting factors that affect the brain in 35 countries worldwide.

In this review, we discuss recent work by the ENIGMA Consortium (http://enigma.ini.usc.edu) - a global alliance of over 500 scientists spread across 200 institutions in 35 countries collectively analyzing brain imaging, clinical, and genetic data. Initially formed to detect genetic influences on brain measures, ENIGMA has grown to over 30 working groups studying 12 major brain diseases by pooling and comparing brain data. In some of the largest neuroimaging studies to date - of schizophrenia and major depression - ENIGMA has found replicable disease effects on the brain that are consistent worldwide, as well as factors that modulate disease effects. In partnership with other consortia including ADNI, CHARGE, IMAGEN and others(1), ENIGMA's genomic screens - now numbering over 30,000 MRI scans - have revealed at least 8 genetic loci that affect brain volumes. Downstream of gene findings, ENIGMA has revealed how these individual variants - and genetic variants in general - may affect both the brain and risk for a range of diseases. The ENIGMA consortium is discovering factors that consistently affect brain structure and function that will serve as future predictors linking individual brain scans and genomic data. It is generating vast pools of normative data on brain measures - from tens of thousands of people - that may help detect deviations from normal development or aging in specific groups of subjects. We discuss challenges and opportunities in applying these predictors to individual subjects and new cohorts, as well as lessons we have learned in ENIGMA's efforts so far.

Genetics of Healthy Aging in Europe.The EU-Integrated Project GEHA (GEnetics of Healthy Aging)

The aim of the 5-year European Union (EU)-Integrated Project GEnetics of Healthy Aging (GEHA), constituted by 25 partners (24 from Europe plus the Beijing Genomics Institute from China), is to identify genes involved in healthy aging and longevity, which allow individuals to survive to advanced old age in good cognitive and physical function and in the absence of major age-related diseases. To achieve this aim a coherent, tightly integrated program of research that unites demographers, geriatricians, geneticists, genetic epidemiologists, molecular biologists, bioinfomaticians, and statisticians has been set up. The working plan is to: (a) collect DNA and information on the health status from an unprecedented number of long-lived 90+ sibpairs (n = 2650) and of younger ethnically matched controls (n = 2650) from 11 European countries; (b) perform a genome-wide linkage scannning in all the sibpairs (a total of 5300 individuals); this investigation will be followed by linkage disequilibrium mapping (LD mapping) of the candidate chromosomal regions; (c) study in cases (i.e., the 2650 probands of the sibpairs) and controls (2650 younger people), genomic regions (chromosome 4, D4S1564, chromosome 11, 11.p15.5) which were identified in previous studies as possible candidates to harbor longevity genes; (d) genotype all recruited subjects for apoE polymorphisms; and (e) genotype all recruited subjects for inherited as well as epigenetic variability of the mitochondrial DNA (mtDNA). The genetic analysis will be performed by 9 high-throughput platforms, within the framework of centralized databases for phenotypic, genetic, and mtDNA data. Additional advanced approaches (bioinformatics, advanced statistics, mathematical modeling, functional genomics and proteomics, molecular biology, molecular genetics) are envisaged to identify the gene variant(s) of interest. The experimental design will also allow (a) to identify gender-specific genes involved in healthy aging and longevity in women and men stratified for ethnic and geographic origin and apoE genotype; (b) to perform a longitudinal survival study to assess the impact of the identified genetic loci on 90+ people mortality; and (c) to develop mathematical and statistical models capable of combining genetic data with demographic characteristics, health status, socioeconomic factors, lifestyle habits.

GENOTYPE-PHENOTYPE ASSOCIATIONS IN NEOVASCULAR AGE-RELATED MACULAR DEGENERATION.

PURPOSE:To examine associations between recognized genetic susceptibility loci and angiographic subphenotypes of the neovascular variant of age-related macular degeneration (nvAMD).METHODS:Participants (247 nvAMD, 52 early age-related macular degeneration [AMD], and 103 controls) were genotyped (complement factor H and ARMS2/HTRA1). nvAMD participants were assigned to one of two subcategories: mainly classic or mainly occult (based on the proportions of classic and occult choroidal neovascularization). nvAMD and early AMD were reassigned to two groups based on the extent and severity of drusen (retinal pigment epithelium dysfunction or not). Univariate and multivariate analysis were used to examine for associations between participant characteristics and genetic loci after adjusting for age, smoking status, and history of cardiovascular disease.RESULTS:Univariate analysis confirmed the known significant associations between AMD stage and age, hypertension, and a history of cardiovascular disease. Those with retinal pigment epithelium dysfunction (F = 5.46; P = 0.02) or a positive smoking history (F = 3.89; P = 0.05) were more likely to have been classified as having mainly an occult rather than a mainly classic lesion. Multivariate analysis showed that significant associations were noted with the number of ARMS2/HTRA1 risk alleles (P

Common variants near FRK/COL10A1 and VEGFA are associated with advanced age-related macular degeneration

Despite significant progress in the identification of genetic loci for age-related macular degeneration (AMD), not all of the heritability has been explained. To identify variants which contribute to the remaining genetic susceptibility, we performed the largest meta-analysis of genome-wide association studies to date for advanced AMD. We imputed 6 036 699 single-nucleotide polymorphisms with the 1000 Genomes Project reference genotypes on 2594 cases and 4134 controls with follow-up replication of top signals in 5640 cases and 52 174 controls. We identified two new common susceptibility alleles, rs1999930 on 6q21-q22.3 near FRK/COL10A1 [odds ratio (OR) 0.87; P = 1.1 × 10(-8)] and rs4711751 on 6p12 near VEGFA (OR 1.15; P = 8.7 × 10(-9)). In addition to the two novel loci, 10 previously reported loci in ARMS2/HTRA1 (rs10490924), CFH (rs1061170, and rs1410996), CFB (rs641153), C3 (rs2230199), C2 (rs9332739), CFI (rs10033900), LIPC (rs10468017), TIMP3 (rs9621532) and CETP (rs3764261) were confirmed with genome-wide significant signals in this large study. Loci in the recently reported genes ABCA1 and COL8A1 were also detected with suggestive evidence of association with advanced AMD. The novel variants identified in this study suggest that angiogenesis (VEGFA) and extracellular collagen matrix (FRK/COL10A1) pathways contribute to the development of advanced AMD.

Epigenetics:Time to translate into transplantation

Substantial progress has been made in identifying genetic loci associated with multifactorial disorders, including variants that seem to impact outcomes following solid organ transplantation. Despite these advances, much of the heritability and susceptibility to chronic disease processes remains unexplained. Epigenetic modifications may exert their effect independently or complementary to genetic variants. Epigenetic modifications can change gene expression without altering the DNA sequence. These modifications are dynamic, potentially heritable, and can be induced by environmental stimuli or drugs. The impact of epigenetic phenomena on the outcomes of organ transplantation is currently poorly understood. Epigenetic modifications can occur during periods of illness; these may persist and potentially influence allograft outcomes. Epigenetic mechanisms influence the activation, proliferation, and differentiation of the immune cells involved in allograft rejection. The donor's epigenome may also impact transplant survival, and initial research has demonstrated that peritransplant conditions induce rapid epigenetic modification within the allograft. Further research will help to define the importance of epigenetic modifications in transplantation. This will potentially lead to the identification of useful biomarkers and the development of novel pharmacotherapies. This review explores the nature of epigenetic modification in disease and the emerging evidence for epigenetic influences on allograft survival.

Rice-arsenate interactions in hydroponics: a three-gene model for tolerance

In this study, the genetic mapping of the tolerance of root growth to 13.3 muM arsenate [As(V)] using the BalaxAzucena population is improved, and candidate genes for further study are identified. A remarkable three-gene model of tolerance is advanced, which appears to involve epistatic interaction between three major genes, two on chromosome 6 and one on chromosome 10. Any combination of two of these genes inherited from the tolerant parent leads to the plant having tolerance. Lists of potential positional candidate genes are presented. These are then refined using whole genome transcriptomics data and bioinformatics. Physiological evidence is also provided that genes related to phosphate transport are unlikely to be behind the genetic loci conferring tolerance. These results offer testable hypotheses for genes related to As(V) tolerance that might offer strategies for mitigating arsenic (As) accumulation in consumed rice.

Association Analysis of Dyslipidemia-Related Genes in Diabetic Nephropathy

Type 1 diabetes (T1D) increases risk of the development of microvascular complications and cardiovascular disease (CVD). Dyslipidemia is a common risk factor in the pathogenesis of both CVD and diabetic nephropathy (DN), with CVD identified as the primary cause of death in patients with DN. In light of this commonality, we assessed single nucleotide polymorphisms (SNPs) in thirty-seven key genetic loci previously associated with dyslipidemia in a T1D cohort using a casecontrol design. SNPs (n = 53) were genotyped using Sequenom in 1467 individuals with T1D (718 cases with proteinuric nephropathy and 749 controls without nephropathy i.e. normal albumin excretion). Cases and controls were white and recruited from the UK and Ireland. Association analyses were performed using PLINK to compare allele frequencies in cases and controls. In a sensitivity analysis, samples from control individuals with reduced renal function (estimated glomerular filtration rate,60 ml/min/1.73 m2) were excluded. Correction for multiple testing was performed by permutation testing. A total of 1394 samples passed quality control filters. Following regression analysis adjusted by collection center, gender, duration of diabetes, and average HbA1c, two SNPs were significantly associated with DN. rs4420638 in the APOC1 region (odds ratio [OR] = 1.51; confidence intervals [CI]: 1.19–1.91; P = 0.001) and rs1532624 in CETP (OR = 0.82; CI: 0.69–0.99; P = 0.034); rs4420638 was also significantly associated in a sensitivity analysis (P = 0.016) together with rs7679 (P = 0.027). However, no association was significant following correction for multiple testing. Subgroup analysis of end-stage renal disease status failed to reveal any association. Our results suggest common variants associated with dyslipidemia are not strongly associated with DN in T1D among white individuals. Our findings, cannot entirely exclude these key genes which are central to the process of dyslipidemia, from involvement in DN pathogenesis as our study had limited power to detect variants of small effect size. Analysis in larger independent cohorts is required.

Suggestive association with ocular phoria at chromosome 6p22

PURPOSE. We conducted a genome-wide association study to identify genetic factors that contribute to the etiology of heterophoria.

METHODS. We measured near and far vertical and horizontal phorias in 988 healthy adults aged 16 to 40 using the Keystone telebinocular with plates 5218 and 5219. We regressed degree of phoria against genotype at 642758 genetic loci. To control for false positives, we applied the conservative genome-wide permutation test to our data.

RESULTS. A locus at 6p22.2 was found to be associated with the degree of near horizontal phoria (P = 2.3 × 10 ). The P value resulting from a genome-wide permutation test was 0.014.

CONCLUSIONS. The strongest association signal arose from an intronic region of the gene ALDH5A1, which encodes the mitochondrial enzyme succinic semialdehyde dehydrogenase (SSADH), an enzyme involved in γ-aminobutyric acid metabolism. Succinic semialdehyde dehydrogenase deficiency, resulting from mutations of ALDH5A1, causes a variety of neural and behavioral abnormalities, including strabismus. Variation in ALDH5A1 is likely to contribute to degree of horizontal phoria.

Deficiency of the ferrous iron transporter FeoAB is linked with metronidazole resistance in Bacteroides fragilis

Background Metronidazole is the most commonly used antimicrobial for Bacteroides fragilis infections and is recommended for prophylaxis of colorectal surgery. Metronidazole resistance is increasing and the mechanisms of resistance are not clear.

Methods A transposon mutant library was generated in B. fragilis 638R (BF638R) to identify the genetic loci associated with resistance to metronidazole.

Results Thirty-two independently isolated metronidazole-resistant mutants had a transposon insertion in BF638R_1421 that encodes the ferrous transport fusion protein (feoAB). Deletion of feoAB resulted in a 10-fold increased MIC of metronidazole for the strain. The metronidazole MIC for the feoAB mutant was similar to that for the parent strain when grown on media supplemented with excess iron, suggesting that the increase seen in the MIC of metronidazole was due to reduced cellular iron transport in the feoAB mutant. The furA gene repressed feoAB transcription in an iron-dependent manner and disruption of furA resulted in constitutive transcription of feoAB, regardless of whether or not iron was present. However, disruption of feoAB also diminished the capacity of BF638R to grow in a mouse intraperitoneal abscess model, suggesting that inorganic ferrous iron assimilation is essential for B. fragilis survival in vivo.

Conclusions Selection for feoAB mutations as a result of metronidazole treatment will disable the pathogenic potential of B. fragilis and could contribute to the clinical efficacy of metronidazole. While mutations in feoAB are probably not a direct cause of clinical resistance, this study provides a key insight into intracellular metronidazole activity and the link with intracellular iron homeostasis.

Identification of putative markers of triclabendazole resistance by a genome-wide analysis of genetically recombinant Fasciola hepatica

Despite years of investigation into triclabendazole (TCBZ) resistance in Fasciola hepatica, the genetic mechanisms responsible remain unknown. Extensive analysis of multiple triclabendazole-susceptible and -resistant isolates using a combination of experimental in vivo and in vitro approaches has been carried out, yet few, if any, genes have been demonstrated experimentally to be associated with resistance phenotypes in the field. In this review we summarize the current understanding of TCBZ resistance from the approaches employed to date. We report the current genomic and genetic resources for F. hepatica that are available to facilitate novel functional genomics and genetic experiments for this parasite in the future. Finally, we describe our own non-biased approach to mapping the major genetic loci involved in conferring TCBZ resistance in F. hepatica.

Independent genomewide screens identify the tumor suppressor VTRNA2-1 as a human epiallele responsive to periconceptional environment

Background: Interindividual epigenetic variation that occurs systemically must be established prior to gastrulation in the very early embryo and, because it is systemic, can be assessed in easily biopsiable tissues. We employ two independent genome-wide approaches to search for such variants.

Results: First, we screen for metastable epialleles by performing genomewide bisulfite sequencing in peripheral blood lymphocyte (PBL) and hair follicle DNA from two Caucasian adults. Second, we conduct a genomewide screen for genomic regions at which PBL DNA methylation is affected by season of conception in rural Gambia. Remarkably, both approaches identify the genomically imprinted VTRNA2-1 as a top environmentally responsive epiallele. We demonstrate systemic and stochastic interindividual variation in DNA methylation at the VTRNA2-1 differentially methylated region in healthy Caucasian and Asian adults and show, in rural Gambians, that periconceptional environment affects offspring VTRNA2-1 epigenotype, which is stable over at least 10 years. This unbiased screen also identifies over 100 additional candidate metastable epialleles, and shows that these are associated with cis genomic features including transposable elements.

Conclusions: The non-coding VTRNA2-1 transcript (also called nc886) is a putative tumor suppressor and modulator of innate immunity. Thus, these data indicating environmentally induced loss of imprinting at VTRNA2-1 constitute a plausible causal pathway linking early embryonic environment, epigenetic alteration, and human disease. More broadly, the list of candidate metastable epialleles provides a resource for future studies of epigenetic variation and human disease.