Tissue-specific genetic variation suggests distinct molecular pathways between body shape phenotypes and colorectal cancer

It remains unknown whether adiposity subtypes are differentially associated with colorectal cancer (CRC). To move beyond single-trait anthropometric indicators, we derived four multi-trait body shape phenotypes reflecting adiposity subtypes from principal components analysis on body mass index, height, weight, waist-to-hip ratio, and waist and hip circumference. A generally obese (PC1) and a tall, centrally obese (PC3) body shape were both positively associated with CRC risk in observational analyses in 329,828 UK Biobank participants (3728 cases). In genome-wide association studies in 460,198 UK Biobank participants, we identified 3414 genetic variants across four body shapes and Mendelian randomization analyses confirmed positive associations of PC1 and PC3 with CRC risk (52,775 cases/45,940 controls from GECCO/CORECT/CCFR). Brain tissue-specific genetic instruments, mapped to PC1 through enrichment analysis, were responsible for the relationship between PC1 and CRC, while the relationship between PC3 and CRC was predominantly driven by adipose tissue-specific genetic instruments. This study suggests distinct putative causal pathways between adiposity subtypes and CRC.

Identification of taxonomic changes in the fecal bacteriome associated with colorectal polyps and cancer: potential biomarkers for early diagnosis

Colorectal cancer (CRC) commonly arises in individuals with premalignant colon lesions known as polyps, with both conditions being influenced by gut microbiota. Host-related factors and inherent characteristics of polyps and tumors may contribute to microbiome variability, potentially acting as confounding factors in the discovery of taxonomic biomarkers for both conditions. In this study we employed shotgun metagenomics to analyze the taxonomic diversity of bacteria present in fecal samples of 90 clinical subjects (comprising 30 CRC patients, 30 with polyps and 30 controls). Our findings revealed a decrease in taxonomic richness among individuals with polyps and CRC, with significant dissimilarities observed among the study groups. We identified significant alterations in the abundance of specific taxa associated with polyps (Streptococcaceae, Lachnoclostridium, and Ralstonia) and CRC (Lactobacillales, Clostridiaceae, Desulfovibrio, SFB, Ruminococcus, and Faecalibacterium). Clostridiaceae exhibited significantly lower abundance in the early stages of CRC. Additionally, our study revealed a positive co-occurrence among underrepresented genera in CRC, while demonstrating a negative co-occurrence between Faecalibacterium and Desulfovibrio, suggesting potential antagonistic relationships. Moreover, we observed variations in taxonomic richness and/or abundance within the polyp and CRC bacteriome linked to polyp size, tumor stage, dyslipidemia, diabetes with metformin use, sex, age, and family history of CRC. These findings provide potential new biomarkers to enhance early CRC diagnosis while also demonstrating how intrinsic host factors contribute to establishing a heterogeneous microbiome in patients with CRC and polyps.

A 1 year course of starch- and sucrose-reduced diet used by irritable bowel syndrome patients with diarrhoea and the effect of genetic variants

Background

Irritable bowel syndrome is a heterogeneous syndrome and it is difficult to find an effective treatment. Previously, a starch- and sucrose-reduced diet (SSRD) demonstrated promising short-term outcomes. It was proposed that genetic variants in the sucrose-isomaltase gene might influence this success. Our aim in this work was to extend the follow-up study to 1 year and to analyse the effect of the genetic variants of genes involved in starch and sucrose metabolism.

Methods

IBS-SSS questionnaire, IBS-QoL questionnaire and questionnaires about adherence, difficulty and food assessment were sent to 34 patients after 6 months and 1 year after the end of the dietary intervention. In addition, 11 genes involved in sucrose and starch metabolism were sequenced.

Results

Twenty-three participants responded to the 6 months follow-up and 16 to the 1 year follow-up. IBS-SSS total value increased 59.71% in the 6 months follow-up compared with the end of the intervention (p = 0.0018), and 55.39% in the 1 year follow-up (p = 0.0166); while IBS-QoL score decreased 24.09% (p = 0.0002) and 18.07% (p = 0.0022), respectively. The adherence decreased by 29.11% (p = 4.8 × 10-5) and 27.21% (p = 0.0054), respectively. In addition, carriers of pathogenic variants on the SI gene showed a slightly better performance than non-carriers. Finally, the participants showed less satisfaction over time with 18 allowed foods in the diet.

Conclusion

Over time the SSRD is difficult to follow and the genotype might affect the performance of the diet. Since this diet could be a promising therapeutic option, a larger cohort needs to be analysed to validate the results and to compare it with other diets.

Genetic Variants as Predictors of the Success of Colorectal Cancer Treatments

BACKGROUND

Some genetic polymorphisms (SNPs) have been proposed as predictors for different colorectal cancer (CRC) outcomes. This work aims to assess their performance in our cohort and find new SNPs associated with them.

METHODS

A total of 833 CRC cases were analyzed for seven outcomes, including the use of chemotherapy, and stratified by tumor location and stage. The performance of 63 SNPs was assessed using a generalized linear model and area under the receiver operating characteristic curve, and local SNPs were detected using logistic regressions.

RESULTS

In total 26 of the SNPs showed an AUC > 0.6 and a significant association (p < 0.05) with one or more outcomes. However, clinical variables outperformed some of them, and the combination of genetic and clinical data showed better performance. In addition, 49 suggestive (p < 5 × 10-6) SNPs associated with one or more CRC outcomes were detected, and those SNPs were located at or near genes involved in biological mechanisms associated with CRC.

CONCLUSIONS

Some SNPs with clinical data can be used in our population as predictors of some CRC outcomes, and the local SNPs detected in our study could be feasible markers that need further validation as predictors.

From Omic Layers to Personalized Medicine in Colorectal Cancer: The Road Ahead

Colorectal cancer is a major health concern since it is a highly diagnosed cancer and the second cause of death among cancers. Thus, the most suitable biomarkers for its diagnosis, prognosis, and treatment have been studied to improve and personalize the prevention and clinical management of colorectal cancer. The emergence of omic techniques has provided a great opportunity to better study CRC and make personalized medicine feasible. In this review, we will try to summarize how the analysis of the omic layers can be useful for personalized medicine and the existing difficulties. We will discuss how single and multiple omic layer analyses have been used to improve the prediction of the risk of CRC and its outcomes and how to overcome the challenges in the use of omic layers in personalized medicine.

The effect of starch- and sucrose-reduced diet accompanied by nutritional and culinary recommendations on the symptoms of irritable bowel syndrome patients with diarrhoea

Background

Irritable bowel syndrome (IBS) is a common gastrointestinal condition which entails a high burden in the quality of life (QoL) of patients. Nutritional interventions have been proposed to alleviate symptoms, since still no effective treatments exist for IBS.

Objectives

Our aim is to analyse the feasibility of the use of starch- and sucrose-reduced diet (SSRD).

Design

In this study, we used a SSRD accompanied by nutritional and culinary recommendations to measure the effects in IBS patients with diarrhoea.

Methods

In all, 34 participants completed a 4-week nutritional intervention based on SSRD. Symptoms, QoL and dietary habits were assessed by several questionnaires that were completed at the beginning, daily, after 2 weeks, at the end, and after 2 months.

Results

85.29% of the participants reached the primary endpoint [reduction of 50 points or more in IBS-symptom severity scale (SSS)], and 58.82% the secondary endpoint (reduction of 50% or more in IBS-SSS). The relief of symptoms and improvement of the QoL were significant after 2 weeks of intervention, at the end and after 2 months. Dietary habits were consistent with the diet and high adherence was achieved.

Conclusions

SSRD and individualized nutritional and culinary guidance improved symptoms and QoL of IBS patients with diarrhoea, with a high adherence.

Host Genetics and Microbiota Interactions in Colorectal Cancer: Shared or Independent Risk?

The role of microbiota in colorectal cancer has been studied since alterations in its composition were observed. In addition, there are more and more pieces of evidence that microbiota could be implicated in colorectal cancer progression. Thus, the components of the microbiota could be biomarkers for the diagnosis and prognosis of colorectal cancer. In addition, it is important to address how the microbiota interacts with the host and how the host shapes the microbiota, in order to understand the biological pathways and mechanisms involved in their relationship and the consequences of their interactions in colorectal cancer. Thereby, it could be possible to find feasible measures and treatments to prevent or better diagnose colorectal cancer. In this review, we will try to summarize the role of the microbiota in colorectal cancer and its interactions with the host and the host genetics, coming to some conclusions that could be useful to find the gaps in our knowledge and propose future steps in this field.

Detailed stratified GWAS analysis for severe COVID-19 in four European populations

Given the highly variable clinical phenotype of Coronavirus disease 2019 (COVID-19), a deeper analysis of the host genetic contribution to severe COVID-19 is important to improve our understanding of underlying disease mechanisms. Here, we describe an extended genome-wide association meta-analysis of a well-characterized cohort of 3255 COVID-19 patients with respiratory failure and 12 488 population controls from Italy, Spain, Norway and Germany/Austria, including stratified analyses based on age, sex and disease severity, as well as targeted analyses of chromosome Y haplotypes, the human leukocyte antigen region and the SARS-CoV-2 peptidome. By inversion imputation, we traced a reported association at 17q21.31 to a ~0.9-Mb inversion polymorphism that creates two highly differentiated haplotypes and characterized the potential effects of the inversion in detail. Our data, together with the 5th release of summary statistics from the COVID-19 Host Genetics Initiative including non-Caucasian individuals, also identified a new locus at 19q13.33, including NAPSA, a gene which is expressed primarily in alveolar cells responsible for gas exchange in the lung.

Applicability of polygenic risk scores in endometriosis clinical presentation

Background

Risk prediction is an essential part of preventative medicine and in recent years genomic information has become an interesting factor in risk models. Polygenic risk scores (PRS) combine the effect of many genetic variations into a single score which has been shown to have predictive value for many diseases. This study aimed to investigate the association between PRS for endometriosis and the clinical presentation of the disease.

Methods

Women with endometriosis (N = 172) were identified at the Department of Gynecology. All participants answered questionnaires regarding sociodemographic factors, lifestyle habits and medical history, registered bowel symptoms on the Visual Analog Scale for Irritable Bowel Syndrome and passed blood samples. DNA was extracted and samples were genotyped, and a PRS was calculated based on previous genome-wide association studies of endometriosis. Inflammatory proteins and TSH receptor antibodies (TRAb) in serum were analyzed.

Results

Inverse associations were identified between PRS and spread of endometriosis, involvement of the gastrointestinal tract and hormone treatment. However, significance was lost when calculated as p for trend and the specificity and sensitivity were low. There were no correlations between PRS and TRAb or inflammatory proteins.

Conclusion

The findings indicate that specific PRS should be developed to predict clinical presentations in patient with endometriosis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12905-022-01788-w.

A survey of functional dyspepsia in 361,360 individuals: Phenotypic and genetic cross-disease analyses

BACKGROUND

Functional dyspepsia (FD) is a common gastrointestinal condition of poorly understood pathophysiology. While symptoms' overlap with other conditions may indicate common pathogenetic mechanisms, genetic predisposition is suspected but has not been adequately investigated.

METHODS

Using healthcare, questionnaire, and genetic data from three large population-based biobanks (UK Biobank, EGCUT, and MGI), we surveyed FD comorbidities, heritability, and genetic correlations across a wide spectrum of conditions and traits in 10,078 cases and 351,282 non-FD controls of European ancestry.

KEY RESULTS

In UK Biobank, 281 diagnoses were detected at increased prevalence in FD, based on healthcare records. Among these, gastrointestinal conditions (OR = 4.0, p < 1.0 × 10^-300 ), anxiety disorders (OR = 2.3, p < 1.4 × 10^-27 ), ischemic heart disease (OR = 2.2, p < 2.3 × 10^-76 ), and infectious and parasitic diseases (OR = 2.1, p = 1.5 × 10^-73 ) showed strongest association with FD. Similar results were obtained in an analysis of self-reported conditions and use of medications from questionnaire data. Based on a genome-wide association meta-analysis of genotypes across all cohorts, FD heritability was estimated close to 5% ( hSNP2  = 0.047, p = 0.014). Genetic correlations indicate FD predisposition is shared with several other diseases and traits (r_g  > 0.344), mostly overlapping with those also enriched in FD patients. Suggestive (p < 5.0 × 10^-6 ) association with FD risk was detected for 13 loci, with 2 showing nominal replication (p < 0.05) in an independent cohort of 192 FD patients.

CONCLUSIONS & INFERENCES

FD has a weak heritable component that shows commonalities with multiple conditions across a wide spectrum of pathophysiological domains. This new knowledge contributes to a better understanding of FD etiology and may have implications for improving its treatment.

Interplay between Genome, Metabolome and Microbiome in Colorectal Cancer

Simple Summary

The development of colorectal cancer (CRC) is influenced by the environment, genetics and microbiota. Microbiome and metabolome analyses allowed for the finding of factors and markers associated with adenoma and CRC risk, but the interaction of host genomics with those omic layers remains unclear. Thus, our aim is to add host genome information to find new factors and markers associated with adenoma and CRC risk or to propose biological mechanisms involved in the risk. We found interactions between different omic layers that could be biologically relevant, and the three layers gave complementary information to predict adenoma and CRC risk. Our findings will help to find new markers for adenoma and CRC risk and to analyze biological mechanisms involved in adenoma and CRC development.

Abstract

Background: Colorectal cancer (CRC), a major health concern, is developed depending on environmental, genetic and microbial factors. The microbiome and metabolome have been analyzed to study their role in CRC. However, the interplay of host genetics with those layers in CRC remains unclear. Methods: 120 individuals were sequenced and association analyses were carried out for adenoma and CRC risk, and for selected components of the microbiome and metabolome. The epistasis between genes located in cholesterol pathways was analyzed; modifiable risk factors were studied using Mendelian randomization; and the three omic layers were used to integrate their data and to build risk prediction models. Results: We detected genetic variants that were associated to components of metabolome or microbiome and adenoma or CRC risk (e.g., in LINC01605, PROKR2 and CCSER1 genes). In addition, we found interactions between genes of cholesterol metabolism, and HDL cholesterol levels affected adenoma (p = 0.0448) and CRC (p = 0.0148) risk. The combination of the three omic layers to build risk prediction models reached high AUC values (>0.91). Conclusions: The use of the three omic layers allowed for the finding of biological mechanisms related to the development of adenoma and CRC, and each layer provided complementary information to build risk prediction models.

GWAS of stool frequency provides insights into gastrointestinal motility and irritable bowel syndrome

Gut dysmotility is associated with constipation, diarrhea, and functional gastrointestinal disorders like irritable bowel syndrome (IBS), although its molecular underpinnings are poorly characterized. We studied stool frequency (defined by the number of bowel movements per day, based on questionnaire data) as a proxy for gut motility in a GWAS meta-analysis including 167,875 individuals from UK Biobank and four smaller population-based cohorts. We identify 14 loci associated with stool frequency (p ≤ 5.0 × 10^-8). Gene set and pathway analyses detected enrichment for genes involved in neurotransmitter/neuropeptide signaling and preferentially expressed in enteric motor neurons controlling peristalsis. PheWAS identified pleiotropic associations with dysmotility syndromes and the response to their pharmacological treatment. The genetic architecture of stool frequency correlates with that of IBS, and UK Biobank participants from the top 1% of stool frequency polygenic score distribution were associated with 5× higher risk of IBS with diarrhea. These findings pave the way for the identification of actionable pathological mechanisms in IBS and the dysmotility syndromes.

Genomewide Association Study of Severe Covid-19 with Respiratory Failure

BACKGROUND

There is considerable variation in disease behavior among patients infected with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the virus that causes coronavirus disease 2019 (Covid-19). Genomewide association analysis may allow for the identification of potential genetic factors involved in the development of Covid-19.

METHODS

We conducted a genomewide association study involving 1980 patients with Covid-19 and severe disease (defined as respiratory failure) at seven hospitals in the Italian and Spanish epicenters of the SARS-CoV-2 pandemic in Europe. After quality control and the exclusion of population outliers, 835 patients and 1255 control participants from Italy and 775 patients and 950 control participants from Spain were included in the final analysis. In total, we analyzed 8,582,968 single-nucleotide polymorphisms and conducted a meta-analysis of the two case-control panels.

RESULTS

We detected cross-replicating associations with rs11385942 at locus 3p21.31 and with rs657152 at locus 9q34.2, which were significant at the genomewide level (P<5×10-8) in the meta-analysis of the two case-control panels (odds ratio, 1.77; 95% confidence interval [CI], 1.48 to 2.11; P = 1.15×10-10; and odds ratio, 1.32; 95% CI, 1.20 to 1.47; P = 4.95×10-8, respectively). At locus 3p21.31, the association signal spanned the genes SLC6A20, LZTFL1, CCR9, FYCO1, CXCR6 and XCR1. The association signal at locus 9q34.2 coincided with the ABO blood group locus; in this cohort, a blood-group-specific analysis showed a higher risk in blood group A than in other blood groups (odds ratio, 1.45; 95% CI, 1.20 to 1.75; P = 1.48×10-4) and a protective effect in blood group O as compared with other blood groups (odds ratio, 0.65; 95% CI, 0.53 to 0.79; P = 1.06×10-5).

CONCLUSIONS

We identified a 3p21.31 gene cluster as a genetic susceptibility locus in patients with Covid-19 with respiratory failure and confirmed a potential involvement of the ABO blood-group system. (Funded by Stein Erik Hagen and others.).

Gene-Diet Interactions in Colorectal Cancer: Survey Design, Instruments, Participants and Descriptive Data of a Case-Control Study in the Basque Country

Epidemiologic studies have revealed inconsistent evidence of gene-diet interaction in relation to colorectal cancer (CRC). The aim of this study was to analyze them in a sample of cases and controls from the population-based bowel cancer screening program of the Osakidetza/Basque Health Service. This study analyzed dietetic, genetic, demographic, socioeconomic factors and lifestyles. In the present manuscript, the survey design, sampling, instruments, measurements and related quality management were presented. Moreover, we analyze differences between cases and controls in some data, especially those related to diet. The participants were 308 cases and 308 age- and sex-matched subjects as controls. Cases were more likely than controls to have overweight/obesity (67.5% vs. 58.1%, p < 0.05), a lower intake of vitamin B₂ (0.86 ± 0.23 vs. 0.92 ± 0.23 mg/1000 kcal, p < 0.01) and calcium:phosphorus ratio (0.62 ± 0.12 vs. 0.65 ± 0.13, p < 0.01). A higher proportion of cases than controls did not meet the Nutritional Objectives for saturated fatty acids (85.7% vs. 67.5%, p < 0.001) or cholesterol (35.4% vs. 25.0%, p < 0.01). In conclusion, the present study provides valuable data for analyzing the complexity of gene-diet interaction in relation to CRC. The results presented here suggest that overweight/obesity and a high intake of certain dietary components, especially saturated fatty acids and cholesterol, are more frequent in cases than in controls.

Implication of m6A mRNA Methylation in Susceptibility to Inflammatory Bowel Disease

Inflammatory bowel disease (IBD) is a chronic inflammatory condition of the gastrointestinal tract that develops due to the interaction between genetic and environmental factors. More than 160 loci have been associated with IBD, but the functional implication of many of the associated genes remains unclear. N6-Methyladenosine (m6A) is the most abundant internal modification in mRNA. m6A methylation regulates many aspects of mRNA metabolism, playing important roles in the development of several pathologies. Interestingly, SNPs located near or within m6A motifs have been proposed as possible contributors to disease pathogenesis. We hypothesized that certain IBD-associated SNPs could regulate the function of genes involved in IBD development via m6A-dependent mechanisms. We used online available GWAS, m6A and transcriptome data to find differentially expressed genes that harbored m6A-SNPs associated with IBD. Our analysis resulted in five candidate genes corresponding to two of the major IBD subtypes: UBE2L3 and SLC22A4 for Crohn's Disease and TCF19, C6orf47 and SNAPC4 for Ulcerative Colitis. Further analysis using in silico predictions and co-expression analyses in combination with in vitro functional studies showed that our candidate genes seem to be regulated by m6A-dependent mechanisms. These findings provide the first indication of the implication of RNA methylation events in IBD pathogenesis.

Food groups, diet quality and colorectal cancer risk in the Basque Country

BACKGROUND

The results obtained to date concerning food groups, diet quality and colorectal cancer (CRC) risk vary according to criteria used and the study populations.

AIM

To study the relationships between food groups, diet quality and CRC risk, in an adult population of the Basque Country (North of Spain).

METHODS

This observational study included 308 patients diagnosed with CRC and 308 age- and sex-matched subjects as controls. During recruitment, dietary, anthropometric, lifestyle, socioeconomic, demographic and health status information was collected. Adherence to the dietary recommendations was evaluated utilizing the Healthy Eating Index for the Spanish Diet and the MedDietScore. Conditional logistic regressions were used to evaluate the associations of food group intakes, diet quality scores, categorized in tertiles, with CRC risk.

RESULTS

The adjusted models for potential confounding factors showed a direct association between milk and dairy products consumption, in particular high-fat cheeses [odds ratio (OR) third tertile vs first tertile = 1.87, 95% confidence intervals (CI): 1.11-3.16], and CRC risk. While the consumption of fiber-containing foods, especially whole grains (OR third tertile vs first tertile = 0.62, 95%CI: 0.39-0.98), and fatty fish (OR third tertile vs first tertile = 0.53, 95%CI: 0.27-0.99) was associated with a lower risk for CRC. Moreover, higher MD adherence was associated with a reduced CRC risk in adjusted models (OR third tertile vs first tertile = 0.40, 95%CI: 0.20-0.80).

CONCLUSION

Direct associations were found for high-fat cheese, whereas an inverse relation was reported for fiber-containing foods and fatty fish, as well as adherence to a Mediterranean dietary pattern.

MAGI2 Gene Region and Celiac Disease

Celiac disease (CD) patients present a loss of intestinal barrier function due to structural alterations in the tight junction (TJ) network, the most apical unions between epithelial cells. The association of TJ-related gene variants points to an implication of this network in disease susceptibility. This work aims to characterize the functional implication of TJ-related, disease-associated loci in CD pathogenesis. We performed an association study of 8 TJ-related gene variants in a cohort of 270 CD and 91 non-CD controls. The expression level of transcripts located in the associated SNP region was analyzed by RT-PCR in several human tissues and in duodenal biopsies of celiac patients and non-CD controls. (si)RNA-driven silencing combined with gliadin in the Caco2 intestinal cell line was used to analyze the implication of transcripts from the associated region in the regulation of TJ genes. We replicated the association of rs6962966^*A variant [p = 0.0029; OR = 1.88 (95%1.24–2.87)], located in an intron of TJ-related MAGI2 coding gene and upstream of RP4-587D13.2 transcript, bioinformatically classified as a long non-coding RNA (lncRNA). The expression of both genes is correlated and constitutively downregulated in CD intestine. Silencing of lncRNA decreases the levels of MAGI2 protein. At the same time, silencing of MAGI2 affects the expression of several TJ-related genes. The associated region is functionally altered in disease, probably affecting CD-related TJ genes.

Combined Analysis of Methylation and Gene Expression Profiles in Separate Compartments of Small Bowel Mucosa Identified Celiac Disease Patients' Signatures

By GWAS studies on celiac disease, gene expression was studied at the level of the whole intestinal mucosa, composed by two different compartments: epithelium and lamina propria. Our aim is to analyse the gene-expression and DNA methylation of candidate genes in each of these compartments. Epithelium was separated from lamina propria in biopsies of CeD patients and CTRs using magnetic beads. Gene-expression was analysed by RT-PC; methylation analysis required bisulfite conversion and NGS. Reverse modulation of gene-expression and methylation in the same cellular compartment was observed for the IL21 and SH2B3 genes in CeD patients relative to CTRs. Bioinformatics analysis highlighted the regulatory elements in the genomic region of SH2B3 that altered methylation levels. The cREL and TNFAIP3 genes showed methylation patterns that were significantly different between CeD patients and CTRs. In CeD, the genes linked to inflammatory processes are up-regulated, whereas the genes involved in the cell adhesion/integrity of the intestinal barrier are down-regulated. These findings suggest a correlation between gene-expression and methylation profile for the IL21 and SH2B3 genes. We identified a “gene-expression phenotype” of CeD and showed that the abnormal response to dietary antigens in CeD might be related not to abnormalities of gene structure but to the regulation of molecular pathways.

The methylome of the celiac intestinal epithelium harbours genotype-independent alterations in the HLA region

The Human Leucocyte Antigen (HLA) locus and other DNA sequence variants identified in Genome-Wide Association (GWA) studies explain around 50% of the heritability of celiac disease (CD). However, the pathogenesis of CD could be driven by other layers of genomic information independent from sequence variation, such as DNA methylation, and it is possible that allele-specific methylation explains part of the SNP associations. Since the DNA methylation landscape is expected to be different among cell types, we analyzed the methylome of the epithelial and immune cell populations of duodenal biopsies in CD patients and controls separately. We found a cell type-specific methylation signature that includes genes mapping to the HLA region, namely TAP1 and HLA-B. We also performed Immunochip SNP genotyping of the same samples and interrogated the expression of some of the affected genes. Our analysis revealed that the epithelial methylome is characterized by the loss of CpG island (CGI) boundaries, often associated to altered gene expression, and by the increased variability of the methylation across the samples. The overlap between differentially methylated positions (DMPs) and CD-associated SNPs or variants contributing to methylation quantitative trait loci (mQTLs) is minimal. In contrast, there is a notable enrichment of mQTLs among the most significant CD-associated SNPs. Our results support the notion that DNA methylation alterations constitute a genotype-independent event and confirm its role in the HLA region (apart from the well-known, DQ allele-specific effect). Finally, we find that a fraction of the CD-associated variants could exert its phenotypic effect through DNA methylation.

Female-Specific Association Between Variants on Chromosome 9 and Self-Reported Diagnosis of Irritable Bowel Syndrome

BACKGROUND & AIMS

Genetic factors are believed to affect risk for irritable bowel syndrome (IBS), but there have been no sufficiently powered and adequately sized studies. To identify DNA variants associated with IBS risk, we performed a genome-wide association study (GWAS) of the large UK Biobank population-based cohort, which includes genotype and health data from 500,000 participants.

METHODS

We studied 7,287,191 high-quality single nucleotide polymorphisms in individuals who self-reported a doctor's diagnosis of IBS (cases; n = 9576) compared to the remainder of the cohort (controls; n = 336,499) (mean age of study subjects, 40-69 years). Genome-wide significant findings were further investigated in 2045 patients with IBS from tertiary centers and 7955 population controls from Europe and the United States, and a small general population sample from Sweden (n = 249). Functional annotation of GWAS results was carried out by integrating data from multiple biorepositories to obtain biological insights from the observed associations.

RESULTS

We identified a genome-wide significant association on chromosome 9q31.2 (single nucleotide polymorphism rs10512344; P = 3.57 × 10-8) in a region previously linked to age at menarche, and 13 additional loci of suggestive significance (P < 5.0×10-6). Sex-stratified analyses revealed that the variants at 9q31.2 affect risk of IBS in women only (P = 4.29 × 10-10 in UK Biobank) and also associate with constipation-predominant IBS in women (P = .015 in the tertiary cohort) and harder stools in women (P = .0012 in the population-based sample). Functional annotation of the 9q31.2 locus identified 8 candidate genes, including the elongator complex protein 1 gene (ELP1 or IKBKAP), which is mutated in patients with familial dysautonomia.

CONCLUSIONS

In a sufficiently powered GWAS of IBS, we associated variants at the locus 9q31.2 with risk of IBS in women. This observation may provide additional rationale for investigating the role of sex hormones and autonomic dysfunction in IBS.

Transcription Factor Binding Site Enrichment Analysis in Co-Expression Modules in Celiac Disease

The aim of this study was to construct celiac co-expression patterns at a whole genome level and to identify transcription factors (TFs) that could drive the gliadin-related changes in coordination of gene expression observed in celiac disease (CD). Differential co-expression modules were identified in the acute and chronic responses to gliadin using expression data from a previous microarray study in duodenal biopsies. Transcription factor binding site (TFBS) and Gene Ontology (GO) annotation enrichment analyses were performed in differentially co-expressed genes (DCGs) and selection of candidate regulators was performed. Expression of candidates was measured in clinical samples and the activation of the TFs was further characterized in C2BBe1 cells upon gliadin challenge. Enrichment analyses of the DCGs identified 10 TFs and five were selected for further investigation. Expression changes related to active CD were detected in four TFs, as well as in several of their in silico predicted targets. The activation of TFs was further characterized in C2BBe1 cells upon gliadin challenge, and an increase in nuclear translocation of CAMP Responsive Element Binding Protein 1 (CREB1) and IFN regulatory factor-1 (IRF1) in response to gliadin was observed. Using transcriptome-wide co-expression analyses we are able to propose novel genes involved in CD pathogenesis that respond upon gliadin stimulation, also in non-celiac models.

Cytoplasmic Form of Carlr lncRNA Facilitates Inflammatory Gene Expression upon NF-κB Activation

Long noncoding RNAs (lncRNAs) have emerged as critical regulators of inflammation. To further understand the interaction between inflammatory signaling pathways and lncRNAs, we characterized the function of cardiac and apoptosis-related lncRNA (Carlr), an lncRNA expressed in both mouse and human cells of diverse tissues. Carlr expression is increased following NF-κB signaling in macrophages, with concomitant translocation to, and enrichment of, the transcript in the cytoplasm. Knockdown of Carlr results in impaired expression of NF-κB pathway genes and influences the interaction between macrophages and intestinal cells in an inflammatory environment. In human celiac disease patient samples, increased levels of the Carlr transcript were detected in the cytoplasm, alongside elevated expression of NF-κB pathway genes. These findings suggest that increased Carlr expression and/or cytoplasmic localization is required for efficient NF-κB signaling and is associated with the inflamed tissue state observed in human celiac disease.

An Activating Mutation in STAT3 Results in Neonatal Diabetes Through Reduced Insulin Synthesis

Neonatal diabetes mellitus (NDM) is a rare form of diabetes diagnosed within the first 6 months of life. Genetic studies have allowed the identification of several genes linked to the development of NDM; however, genetic causes for ∼20% of the cases remain to be clarified. Most cases of NDM involve isolated diabetes, but sometimes NDM appears in association with other pathological conditions, including autoimmune diseases. Recent reports have linked activating mutations in STAT3 with early-onset autoimmune disorders that include diabetes of autoimmune origin, but the functional impact of STAT3-activating mutations have not been characterized at the pancreatic β-cell level. By using whole-exome sequencing, we identified a novel missense mutation in the binding domain of the STAT3 protein in a patient with NDM. The functional analyses showed that the mutation results in an aberrant activation of STAT3, leading to deleterious downstream effects in pancreatic β-cells. The identified mutation leads to hyperinhibition of the transcription factor Isl-1 and, consequently, to a decrease in insulin expression. These findings represent the first functional indication of a direct link between an NDM-linked activating mutation in STAT3 and pancreatic β-cell dysfunction.

Ancestry-based stratified analysis of Immunochip data identifies novel associations with celiac disease

To identify candidate genes in celiac disease (CD), we reanalyzed the whole Immunochip CD cohort using a different approach that clusters individuals based on immunoancestry prior to disease association analysis, rather than by geographical origin. We detected 636 new associated SNPs (P<7.02 × 10^-07) and identified 5 novel genomic regions, extended 8 others previously identified and also detected 18 isolated signals defined by one or very few significant SNPs. To test whether we could identify putative candidate genes, we performed expression analyses of several genes from the top novel region (chr2:134533564-136169524), from a previously identified locus that is now extended, and a gene marked by an isolated SNP, in duodenum biopsies of active and treated CD patients, and non-celiac controls. In the largest novel region, CCNT2 and R3HDM1 were constitutively underexpressed in disease, even after gluten removal. Moreover, several genes within this region were coexpressed in patients, but not in controls. Other novel genes like KIF21B, REL and SORD also showed altered expression in active disease. Apart from the identification of novel CD loci, these results suggest that ancestry-based stratified analysis is an efficient strategy for association studies in complex diseases.

High Expression of Endogenous Retroviral Envelope Gene in the Equine Fetal Part of the Placenta

Endogenous retroviruses (ERVs) are proviral phases of exogenous retroviruses that have co-evolved with vertebrate genomes for millions of years. Previous studies have identified the envelope (env) protein genes of retroviral origin preferentially expressed in the placenta which suggests a role in placentation based on their membrane fusogenic capacity and therefore they have been named syncytins. Until now, all the characterized syncytins have been associated with three invasive placentation types: the endotheliochorial (Carnivora), the synepitheliochorial (Ruminantia), and the hemochorial placentation (human, mouse) where they play a role in the syncytiotrophoblast formation. The purpose of the present study was to evaluate whether EqERV env RNA is expressed in horse tissues as well and investigate if the horse, possessing an epitheliochorial placenta, has "captured" a common retroviral env gene with syncytin-like properties in placental tissues. Interestingly, although in the equine placenta there is no syncytiotrophoblast layer at the maternal-fetal interface, our results showed that EqERV env RNA is highly expressed at that level, as expected for a candidate syncytin-like gene but with reduced abundance in the other somatic tissues (nearly 30-fold lower) thus suggesting a possible role in the placental tissue. Although the horse is one of the few domestic animals with a sequenced genome, few studies have been conducted about the EqERV and their expression in placental tissue has never been investigated.

Genome-wide reexamination of endogenous retroviruses in Rattus norvegicus

Endogenous retroviruses (ERVs) are remnants of retroviral infections that are present in a large number of vertebrate genomes. Based on the proposal that the rat could act as a reservoir of retroviruses, rat ERVs were analysed in silico using a whole-genome approach. To enrich the detected ERV groups, we applied an upgraded approach based on the hidden Markov model. We found 2637 elements that were classified into the following groups: 9 groups of Class I; 15 of Class II, 7 of them previously described; 1 of Class III; and 3 groups whose classification was unclear but were distantly related to Class I. Sixteen ERV groups seemed to be specific to rat. The high number of rat-specific groups might be related to the contact of rats with retroviruses and their role as a reservoir. In addition, the env gene of the more extended groups seemed to be undetectable.

Endogenous retroviruses in domestic animals

Endogenous retroviruses (ERVs) are genomic elements that are present in a wide range of vertebrates. Although the study of ERVs has been carried out mainly in humans and model organisms, recently, domestic animals have become important, and some species have begun to be analyzed to gain further insight into ERVs. Due to the availability of complete genomes and the development of new computer tools, ERVs can now be analyzed from a genome-wide viewpoint. In addition, more experimental work is being carried out to analyze the distribution, expression and interplay of ERVs within a host genome. Cats, cattle, chicken, dogs, horses, pigs and sheep have been scrutinized in this manner, all of which are interesting species in health and economic terms. Furthermore, several studies have noted differences in the number of endogenous retroviruses and in the variability of these elements among different breeds, as well as their expression in different tissues and the effects of their locations, which, in some cases, are near genes. These findings suggest a complex, intriguing relationship between ERVs and host genomes. In this review, we summarize the most important in silico and experimental findings, discuss their implications and attempt to predict future directions for the study of these genomic elements.

Consistency of metagenomic assignment programs in simulated and real data

Background

Metagenomics is the genomic study of uncultured environmental samples, which has been greatly facilitated by the advent of shotgun-sequencing technologies. One of the main focuses of metagenomics is the discovery of previously uncultured microorganisms, which makes the assignment of sequences to a particular taxon a challenge and a crucial step. Recently, several methods have been developed to perform this task, based on different methodologies such as sequence composition or sequence similarity. The sequence composition methods have the ability to completely assign the whole dataset. However, their use in metagenomics and the study of their performance with real data is limited. In this work, we assess the consistency of three different methods (BLAST + Lowest Common Ancestor, Phymm, and Naïve Bayesian Classifier) in assigning real and simulated sequence reads.

Results

Both in real and in simulated data, BLAST + Lowest Common Ancestor (BLAST + LCA), Phymm, and Naïve Bayesian Classifier consistently assign a larger number of reads in higher taxonomic levels than in lower levels. However, discrepancies increase at lower taxonomic levels. In simulated data, consistent assignments between all three methods showed greater precision than assignments based on Phymm or Bayesian Classifier alone, since the BLAST + LCA algorithm performed best. In addition, assignment consistency in real data increased with sequence read length, in agreement with previously published simulation results.

Conclusions

The use and combination of different approaches is advisable to assign metagenomic reads. Although the sensitivity could be reduced, the reliability can be increased by using the reads consistently assigned to the same taxa by, at least, two methods, and by training the programs using all available information.

Evolutionary history of bovine endogenous retroviruses in the Bovidae family

Background

Endogenous retroviruses (ERVs) are genomic elements of retroviral origin that are present in the genomes of almost all vertebrates. In cattle, more than 13,000 elements related to ERVs have been detected, and based on the pol gene, 24 families or groups of bovine ERVs have been described. However, information about ERVs in other bovids and the presence of families of related bovine ERVs in different species of the Bovidae family is scarce.

Results

The 24 families of bovine ERVs previously detected in cattle (Bos taurus) were also detected in zebus (Bos indicus) and yaks (Bos grunniens). In addition, six new families, named BoERV25 to BoERV30, were detected in the three Bos species. Five more ruminant species were screened for related ERVs: 26 families were detected in these species, but four families (BoERV24, BoERV26, BoERV28 and BoERV29) were specific to cattle, zebus, yaks and buffalo. An analysis of the homology of the ERVs of cattle, zebus and yaks revealed that the level of LTR divergence was similar between ERVs from cattle and zebus but was less similar between with ERVs from cattle and yaks. In addition, purifying selection was detected in the genes and retroviral regions of clusters of ERVs of cattle, zebus and yaks.

Conclusions

In this work, the 24 ERV families previously identified in cattle were also found in two other species in the Bos genus. In addition, six new bovine ERV families were detected. Based on LTR divergence, the most recently inserted families are from Class II. The divergence of the LTR, used as an indirect estimate of the ERV insertion time, seemed to be influenced by the differences in genome evolution since the divergence of the species. In addition, purifying selection could be acting on clusters of ERVs from different species.

Human papillomavirus (HPV) genotype 18 variants in patients with clinical manifestations of HPV related infections in Bilbao, Spain

BACKGROUND

Human papillomavirus (HPV) variants differ in their biological and chemical properties, and therefore, may present differences in pathogenicity. Most authors classified variants based on the phylogenetic analysis of L1 region. Nevertheless, recombination in HPV samples is becoming a usual finding and thus, characterizing genetic variability in other regions should be essential.

OBJECTIVES

We aimed to characterize the genetic variability of HPV 18 in 5 genomic regions: E6, E7, E4, L1 and the Upstream Regulatory Region (URR), working with both single infection and multiple HPV infection samples. Furthermore, we aimed to assess the prevalence of HPV 18 variants in our region and look for possible existence of recombination as well as analyze the relationship between these variants and the type of lesion.

METHODS

From 2007 to 2010, Clinical Microbiology and Infection Control Department analyzed 44 samples which were positive for HPV 18. Genetic variability was determined in PCR products and variants were assigned to European, Asian-amerindian or African lineage. Recombination and association of variants with different types of lesion was studied.

RESULTS

Genetic analysis of the regions revealed a total of 56 nucleotide variations. European, African and Asian-amerindian variants were found in 25/44 (56.8%), 10/44 (22.7%) and 5/44 (11.4%) samples, respectively. We detected the presence of recombinant variants in 2/44 (4.5%) cases. Samples taken from high-grade squamous intraepithelial lesions (H-SIL) only presented variants with specific-african substitutions.

CONCLUSIONS

Multiple HPV infection, non-european HPV variants prevalence and existence of recombination are considered risk factors for HPV persistence and progression of intraepithelial abnormalities, and therefore, should be taken into consideration in order to help to design and optimize diagnostics protocols as well as improve epidemiologic studies.Our study is one of the few studies in Spain which analyses the genetic variability of HPV18 and we showed the importance of characterizing more than one genomic region in order to detect recombination and classify HPV variants properly.

Staphylococcus prevails in the skin microbiota of long-term immunodeficient mice

Host-commensal relationships in the skin are a complex system governed by variables related to the host, the bacteria and the environment. A disruption of this system may lead to new steady states, which, in turn, may lead to disease. We have studied one such disruption by characterizing the skin microbiota in healthy and immunodepressed (ID) mice. A detailed anatomopathological study failed to reveal any difference between the skin of healthy and ID mice. We sequenced the 16S rDNA V1-V2 gene region to saturation in 10 healthy and 10 ID 8 week-old mice, and found than all of the healthy and two of the ID mice had bacterial communities that were similar in composition to that of human skin, although, presumably because of the uniform raising conditions, less interindividual variation was found in mice. However, eight ID mice showed microbiota dominated by Staphylococcus epidermidis. Quantitative PCR amplification of 16S rDNA gene and of the Staphylococcus-specific TstaG region confirmed the previous results and indicated that the quantitative levels of Staphylococcus were similar in both groups while the total number of 16S copies was greater in the healthy mice. Thus, it is possible that, under long-term immunodeficiency, which removes the acquired but not the native immune system, S.epidermidis may inhibit the growth of other bacteria but does not cause a pathogenic state.