Large-scale analysis of sheep rumen metagenome profiles captured by reduced representation sequencing reveals individual profiles are influenced by the environment and genetics of the host.

BACKGROUND: Producing animal protein while reducing the animal's impact on the environment, e.g., through improved feed efficiency and lowered methane emissions, has gained interest in recent years. Genetic selection is one possible path to reduce the environmental impact of livestock production, but these traits are difficult and expensive to measure on many animals. The rumen microbiome may serve as a proxy for these traits due to its role in feed digestion. Restriction enzyme-reduced representation sequencing (RE-RRS) is a high-throughput and cost-effective approach to rumen metagenome profiling, but the systematic (e.g., sequencing) and biological factors influencing the resulting reference based (RB) and reference free (RF) profiles need to be explored before widespread industry adoption is possible. RESULTS: Metagenome profiles were generated by RE-RRS of 4,479 rumen samples collected from 1,708 sheep, and assigned to eight groups based on diet, age, time off feed, and country (New Zealand or Australia) at the time of sample collection. Systematic effects were found to have minimal influence on metagenome profiles. Diet was a major driver of differences between samples, followed by time off feed, then age of the sheep. The RF approach resulted in more reads being assigned per sample and afforded greater resolution when distinguishing between groups than the RB approach. Normalizing relative abundances within the sampling Cohort abolished structures related to age, diet, and time off feed, allowing a clear signal based on methane emissions to be elucidated. Genus-level abundances of rumen microbes showed low-to-moderate heritability and repeatability and were consistent between diets. CONCLUSIONS: Variation in rumen metagenomic profiles was influenced by diet, age, time off feed and genetics. Not accounting for environmental factors may limit the ability to associate the profile with traits of interest. However, these differences can be accounted for by adjusting for Cohort effects, revealing robust biological signals. The abundances of some genera were consistently heritable and repeatable across different environments, suggesting that metagenomic profiles could be used to predict an individual's future performance, or performance of its offspring, in a range of environments. These results highlight the potential of using rumen metagenomic profiles for selection purposes in a practical, agricultural setting.

Combining host and rumen metagenome profiling for selection in sheep: prediction of methane, feed efficiency, production, and health traits.

BACKGROUND: Rumen microbes break down complex dietary carbohydrates into energy sources for the host and are increasingly shown to be a key aspect of animal performance. Host genotypes can be combined with microbial DNA sequencing to predict performance traits or traits related to environmental impact, such as enteric methane emissions. Metagenome profiles were generated from 3139 rumen samples, collected from 1200 dual purpose ewes, using restriction enzyme-reduced representation sequencing (RE-RRS). Phenotypes were available for methane (CH4) and carbon dioxide (CO2) emissions, the ratio of CH4 to CH4 plus CO2 (CH4Ratio), feed efficiency (residual feed intake: RFI), liveweight at the time of methane collection (LW), liveweight at 8 months (LW8), fleece weight at 12 months (FW12) and parasite resistance measured by faecal egg count (FEC1). We estimated the proportion of phenotypic variance explained by host genetics and the rumen microbiome, as well as prediction accuracies for each of these traits. RESULTS: Incorporating metagenome profiles increased the variance explained and prediction accuracy compared to fitting only genomics for all traits except for CO2 emissions when animals were on a grass diet. Combining the metagenome profile with host genotype from lambs explained more than 70% of the variation in methane emissions and residual feed intake. Predictions were generally more accurate when incorporating metagenome profiles compared to genetics alone, even when considering profiles collected at different ages (lamb vs adult), or on different feeds (grass vs lucerne pellet). A reference-free approach to metagenome profiling performed better than metagenome profiles that were restricted to capturing genera from a reference database. We hypothesise that our reference-free approach is likely to outperform other reference-based approaches such as 16S rRNA gene sequencing for use in prediction of individual animal performance. CONCLUSIONS: This paper shows the potential of using RE-RRS as a low-cost, high-throughput approach for generating metagenome profiles on thousands of animals for improved prediction of economically and environmentally important traits. A reference-free approach using a microbial relationship matrix from log10 proportions of each tag normalized within cohort (i.e., the group of animals sampled at the same time) is recommended for future predictions using RE-RRS metagenome profiles.

Efficiency of genotyping by sequencing in inferring genomic relatedness and molecular insights into fat tail selection in Tunisian sheep.

In developing countries, the use of simple and cost-efficient molecular technology is crucial for genetic characterization of local animal resources and better development of conservation strategies. The genotyping by sequencing (GBS) technique, also called restriction enzyme- reduced representational sequencing, is an efficient, cost-effective method for simultaneous discovery and genotyping of many markers. In the present study, we applied a two-enzyme GBS protocol (PstI/MspI) to discover and genotype SNP markers among 197 Tunisian sheep samples. A total of 100 333 bi-allelic SNPs were discovered and genotyped with an SNP call rate of 0.69 and mean sample depth 3.33. The genomic relatedness between 183 samples grouped the samples perfectly to their populations and pointed out a high genetic relatedness of inbred subpopulation reflecting the current adopted reproductive strategies. The genome-wide association study contrasting fat vs. thin-tailed breeds detected 41 significant variants including a peak positioned on OAR20. We identified FOXC1, GMDS, VEGFA, OXCT1, VRTN and BMP2 as the most promising for sheep tail-type trait. The GBS data have been useful to assess the population structure and improve our understanding of the genomic architecture of distinctive characteristics shaped by selection pressure in local sheep breeds. This study successfully investigates a cost-efficient method to discover genotypes, assign populations and understand insights into sheep adaptation to arid area. GBS could be of potential utility in livestock species in developing/emerging countries.

Effects of long-acting, broad spectra anthelmintic treatments on the rumen microbial community compositions of grazing sheep.

Anthelmintic treatment of adult ewes is widely practiced to remove parasite burdens in the expectation of increased ruminant productivity. However, the broad activity spectra of many anthelmintic compounds raises the possibility of impacts on the rumen microbiota. To investigate this, 300 grazing ewes were allocated to treatment groups that included a 100-day controlled release capsule (CRC) containing albendazole and abamectin, a long-acting moxidectin injection (LAI), and a non-treated control group (CON). Rumen bacterial, archaeal and protozoal communities at day 0 were analysed to identify 36 sheep per treatment with similar starting compositions. Microbiota profiles, including those for the rumen fungi, were then generated for the selected sheep at days 0, 35 and 77. The CRC treatment significantly impacted the archaeal community, and was associated with increased relative abundances of Methanobrevibacter ruminantium, Methanosphaera sp. ISO3-F5, and Methanomassiliicoccaceae Group 12 sp. ISO4-H5 compared to the control group. In contrast, the LAI treatment increased the relative abundances of members of the Veillonellaceae and resulted in minor changes to the bacterial and fungal communities by day 77. Overall, the anthelmintic treatments resulted in few, but highly significant, changes to the rumen microbiota composition.

A restriction enzyme reduced representation sequencing approach for low-cost, high-throughput metagenome profiling.

Microbial community profiles have been associated with a variety of traits, including methane emissions in livestock. These profiles can be difficult and expensive to obtain for thousands of samples (e.g. for accurate association of microbial profiles with traits), therefore the objective of this work was to develop a low-cost, high-throughput approach to capture the diversity of the rumen microbiome. Restriction enzyme reduced representation sequencing (RE-RRS) using ApeKI or PstI, and two bioinformatic pipelines (reference-based and reference-free) were compared to bacterial 16S rRNA gene sequencing using repeated samples collected two weeks apart from 118 sheep that were phenotypically extreme (60 high and 58 low) for methane emitted per kg dry matter intake (n = 236). DNA was extracted from freeze-dried rumen samples using a phenol chloroform and bead-beating protocol prior to RE-RRS. The resulting sequences were used to investigate the repeatability of the rumen microbial community profiles, the effect of laboratory and analytical method, and the relationship with methane production. The results suggested that the best method was PstI RE-RRS analyzed with the reference-free approach, which accounted for 53.3±5.9% of reads, and had repeatabilities of 0.49±0.07 and 0.50±0.07 for the first two principal components (PC1 and PC2), phenotypic correlations with methane yield of 0.43±0.06 and 0.46±0.06 for PC1 and PC2, and explained 41±8% of the variation in methane yield. These results were significantly better than for bacterial 16S rRNA gene sequencing of the same samples (p<0.05) except for the correlation between PC2 and methane yield. A Sensitivity study suggested approximately 2000 samples could be sequenced in a single lane on an Illumina HiSeq 2500, meaning the current work using 118 samples/lane and future proposed 384 samples/lane are well within that threshold. With minor adaptations, our approach could be used to obtain microbial profiles from other metagenomic samples.

GBS Data Identify Pigmentation-Specific Genes of Potential Role in Skin-Photosensitization in Two Tunisian Sheep Breeds.

The Tunisian Noire de Thibar sheep breed is a composite breed, recently selected to create animals that are uniformly black in order to avoid skin photosensitization after the ingestion of toxic "hypericum perforatum" weeds, which causes a major economic loss to sheep farmers. We assessed genetic differentiation and estimated marker FST using genotyping-by-sequencing (GBS) data in black (Noire de Thibar) and related white-coated (Queue fine de l'ouest) sheep breeds to identify signals of artificial selection. The results revealed the selection signatures within candidate genes related to coat color, which are assumed to be indirectly involved in the mechanism of photosensitization in sheep. The identified genes could provide important information for molecular breeding.

An entropy-reducing data representation approach for bioinformatic data.

Database URL: https://github.com/AgResearch/data_prism.

Sheep genome functional annotation reveals proximal regulatory elements contributed to the evolution of modern breeds.

Domestication fundamentally reshaped animal morphology, physiology and behaviour, offering the opportunity to investigate the molecular processes driving evolutionary change. Here we assess sheep domestication and artificial selection by comparing genome sequence from 43 modern breeds (Ovis aries) and their Asian mouflon ancestor (O. orientalis) to identify selection sweeps. Next, we provide a comparative functional annotation of the sheep genome, validated using experimental ChIP-Seq of sheep tissue. Using these annotations, we evaluate the impact of selection and domestication on regulatory sequences and find that sweeps are significantly enriched for protein coding genes, proximal regulatory elements of genes and genome features associated with active transcription. Finally, we find individual sites displaying strong allele frequency divergence are enriched for the same regulatory features. Our data demonstrate that remodelling of gene expression is likely to have been one of the evolutionary forces that drove phenotypic diversification of this common livestock species.

GAMOLA2, a Comprehensive Software Package for the Annotation and Curation of Draft and Complete Microbial Genomes.

Expert curated annotation remains one of the critical steps in achieving a reliable biological relevant annotation. Here we announce the release of GAMOLA2, a user friendly and comprehensive software package to process, annotate and curate draft and complete bacterial, archaeal, and viral genomes. GAMOLA2 represents a wrapping tool to combine gene model determination, functional Blast, COG, Pfam, and TIGRfam analyses with structural predictions including detection of tRNAs, rRNA genes, non-coding RNAs, signal protein cleavage sites, transmembrane helices, CRISPR repeats and vector sequence contaminations. GAMOLA2 has already been validated in a wide range of bacterial and archaeal genomes, and its modular concept allows easy addition of further functionality in future releases. A modified and adapted version of the Artemis Genome Viewer (Sanger Institute) has been developed to leverage the additional features and underlying information provided by the GAMOLA2 analysis, and is part of the software distribution. In addition to genome annotations, GAMOLA2 features, among others, supplemental modules that assist in the creation of custom Blast databases, annotation transfers between genome versions, and the preparation of Genbank files for submission via the NCBI Sequin tool. GAMOLA2 is intended to be run under a Linux environment, whereas the subsequent visualization and manual curation in Artemis is mobile and platform independent. The development of GAMOLA2 is ongoing and community driven. New functionality can easily be added upon user requests, ensuring that GAMOLA2 provides information relevant to microbiologists. The software is available free of charge for academic use.

Buccal swabbing as a noninvasive method to determine bacterial, archaeal, and eukaryotic microbial community structures in the rumen.

Analysis of rumen microbial community structure based on small-subunit rRNA marker genes in metagenomic DNA samples provides important insights into the dominant taxa present in the rumen and allows assessment of community differences between individuals or in response to treatments applied to ruminants. However, natural animal-to-animal variation in rumen microbial community composition can limit the power of a study considerably, especially when only subtle differences are expected between treatment groups. Thus, trials with large numbers of animals may be necessary to overcome this variation. Because ruminants pass large amounts of rumen material to their oral cavities when they chew their cud, oral samples may contain good representations of the rumen microbiota and be useful in lieu of rumen samples to study rumen microbial communities. We compared bacterial, archaeal, and eukaryotic community structures in DNAs extracted from buccal swabs to those in DNAs from samples collected directly from the rumen by use of a stomach tube for sheep on four different diets. After bioinformatic depletion of potential oral taxa from libraries of samples collected via buccal swabs, bacterial communities showed significant clustering by diet (R = 0.37; analysis of similarity [ANOSIM]) rather than by sampling method (R = 0.07). Archaeal, ciliate protozoal, and anaerobic fungal communities also showed significant clustering by diet rather than by sampling method, even without adjustment for potentially orally associated microorganisms. These findings indicate that buccal swabs may in future allow quick and noninvasive sampling for analysis of rumen microbial communities in large numbers of ruminants.

Analyses of hypomethylated oil palm gene space.

Demand for palm oil has been increasing by an average of ∼8% the past decade and currently accounts for about 59% of the world's vegetable oil market. This drives the need to increase palm oil production. Nevertheless, due to the increasing need for sustainable production, it is imperative to increase productivity rather than the area cultivated. Studies on the oil palm genome are essential to help identify genes or markers that are associated with important processes or traits, such as flowering, yield and disease resistance. To achieve this, 294,115 and 150,744 sequences from the hypomethylated or gene-rich regions of Elaeis guineensis and E. oleifera genome were sequenced and assembled into contigs. An additional 16,427 shot-gun sequences and 176 bacterial artificial chromosomes (BAC) were also generated to check the quality of libraries constructed. Comparison of these sequences revealed that although the methylation-filtered libraries were sequenced at low coverage, they still tagged at least 66% of the RefSeq supported genes in the BAC and had a filtration power of at least 2.0. A total 33,752 microsatellites and 40,820 high-quality single nucleotide polymorphism (SNP) markers were identified. These represent the most comprehensive collection of microsatellites and SNPs to date and would be an important resource for genetic mapping and association studies. The gene models predicted from the assembled contigs were mined for genes of interest, and 242, 65 and 14 oil palm transcription factors, resistance genes and miRNAs were identified respectively. Examples of the transcriptional factors tagged include those associated with floral development and tissue culture, such as homeodomain proteins, MADS, Squamosa and Apetala2. The E. guineensis and E. oleifera hypomethylated sequences provide an important resource to understand the molecular mechanisms associated with important agronomic traits in oil palm.

Screening for hypogonadism in diabetes 2008/9: results from the Cheshire Primary Care cohort.

AIM: There is increasing awareness of hypogonadism in men with type 2 diabetes but limited data from Primary Care. SUBJECTS AND METHODS: The anonymised records of 6457 male patients aged 18-80 years with diabetes were accessed. Within the last 2 years 391 men (6.0% of total) underwent measurement of serum testosterone. Data search was performed through the centralised data facility afforded by EMIS®, the majority GP systems provider in Cheshire. RESULTS: 4.4% of type 2 diabetes men screened were frankly hypogonadal with a serum total testosterone of less than 8.0 nmol/l. For borderline hypogonadism (serum total testosterone 8-11.99 nmol/l) the proportion of type 2 diabetes men rose to 32.1%. Age adjusted mean (geometric) testosterone was lower in men with type 2 diabetes (13.6 nmol/l 95%CI: 13.1-14.2) vs type 1 diabetes (17.9 nmol/l; 95%CI 15.2-21.0), F=10.3; p=0.0014. For those screened age adjusted body mass index (BMI) was greater in type 2 diabetes at 30.7 (30.1-31.3) vs 28.4 (26.1-30.6)kg/m(2) in type 1 diabetes (F=4.3; p=0.04). Multiple linear regression analysis indicated that there was a statistically significant interaction (P=0.014) between BMI and diabetes type in their relation with log testosterone. For persons with type 1 DM and type 2 DM, testosterone can be expected to decrease by 6% (P=0.002) and by 1% (P=0.002) respectively, for every one unit increment in BMI. CONCLUSIONS: There is manifestly a subset of men with diabetes and androgen deficiency who could benefit from testosterone replacement. BMI has an independent influence on androgen status.

Genetic analysis of MDR1 and inflammatory bowel disease reveals protective effect of heterozygous variants for ulcerative colitis.

BACKGROUND: Single nucleotide polymorphisms (SNPs) in the multidrug transporter MDR1 have been associated with inflammatory bowel disease (IBD) in different studies. However, the data are highly controversial. Recently, 6 haplotype tagging SNPs (tSNPs), representing the haplotype variations of the MDR1 gene, were identified. The aims of this study were to genotype these variants and correlate them to disease phenotype in New Zealand IBD patients. MATERIALS AND METHODS: A total of 784 IBD patients and 200 healthy subjects were genotyped for 5 tSNPs and the triallelic MDR1 variant G2677T/A using the Sequenom MassArray platform. Furthermore, the effects of these variants were examined in correlation with phenotypic clinical features. RESULTS: Heterozygous carriers for the variants C1236T, rs2235046 (an SNP in intron 16), and G2677T/A showed a lower risk of developing ulcerative colitis (C1236T: odds ratio [OR] = 0.63, 95% confidence interval [CI] = 0.42-0.93, P = 0.03; G2677T/A: OR = 0.59, CI = 0.39-0.89, P = 0.02; and rs2235046: OR = 0.59, 95% CI = 0.38-0.91, P = 0.009) as compared with homozygotes. None of the analyzed markers were associated with Crohn's disease on a genotypic level. Subgroup analysis revealed an association for 2 variants with IBD when stratified for age of onset (C1236T SNP and rs3789243). The MDR1 variant C3435T was associated with disease behavior in CD (OR = 1.45, 95% CI = 1.01-2.08, P = 0.04), whereas the SNP rs3789243 was found to be associated with pancolitis in UC patients (OR = 1.35, CI = 1.00-1.82, P = 0.05). CONCLUSIONS: The results of our study support the role of MDR1 as a candidate gene for ulcerative colitis. Furthermore, our results suggest the possibility of a heterozygous advantage for certain MDR1 variants for this disease.

Mushroom intolerance: a novel diet-gene interaction in Crohn's disease.

Carrying a functional single nucleotide polymorphism (L503F, c. 1672 C>T) in the gene for the Na-dependent organic cation transporter (OCTN1), increases the risk of Crohn's disease (CD) in some, but not all, populations. Case-control data on New Zealand Caucasians show no differences for CD risk between individuals carrying the L503F OCTN1 C-allele when compared with those carrying the variant T-allele. However, more of the New Zealand CD cases report intolerance to maize and mushrooms than those who report beneficial effects or no differences. The OCTN1 gene encodes a transporter for ergothionine, a fungal metabolite at high levels in mushrooms but not widely common in other dietary items. An inability to tolerate mushrooms showed statistically significant associations with the variant OCTN1 genotype. That is, among those individuals reporting adverse effects from mushrooms, there was a higher frequency of the variant T-allele when compared with the general population, or with CD patients overall. We believe that this is a novel gene-diet association, suggesting that individuals carrying the OCTN1 variant single nucleotide polymorphism may have an enhanced risk of adverse symptoms associated with consuming mushrooms. Nutrigenomic approaches to dietary recommendations may be appropriate in this group.

Tumor necrosis factor receptor superfamily, member 1B haplotypes increase or decrease the risk of inflammatory bowel diseases in a New Zealand caucasian population.

Inflammatory bowel diseases (IBDs) comprising Crohn disease (CD) and ulcerative colitis (UC) are chronic inflammatory conditions with polygenic susceptibility. Interactions between TNF-alpha and TNF-alpha receptor play a fundamental role in inflammatory response. This study investigates the role that selected single nucleotide polymorphisms (SNPs) and haplotypes in the TNF-alpha receptor (TNSFRSF1B) gene play in the risk of IBD in a New Zealand Caucasian population. DNA samples from 388 CD, 405 UC, 27 indeterminate colitis patients, and 293 randomly selected controls, from Canterbury, New Zealand were screened for 3 common SNPs in TNSFRSF1B: rs1061622 (c.676T > C), rs1061624 (c.*1663A > G), and rs3397 (c.*1690T > C), using TaqMan technologies. Carrying the rs1061624 variant decreased the risk of UC in the left colon (OR 0.73, 95% CI = 0.54-1.00) and of being a smoker at diagnosis (OR 0.62; 95% CI = 0.40-0.96). Carrying the rs3397 variant decreased the risk of penetrating CD (OR 0.62, 95% CI = 0.40-0.95). Three marker haplotype analyses revealed highly significant differences between CD patients and control subjects (chi(2) = 29.9, df = 7, P = .0001) and UC cases and controls (chi(2) = 46.3, df = 7, P < .0001). We conclude that carrying a 3-marker haplotype in the TNSFRSF1B gene may increase (e.g., haplotype of GGC was 2.9-fold more in the CD or UCpatients) or decrease (e.g., TGT was 0.47-fold less in UC patients) the risk of IBD in a New Zealand Caucasian population.

Proteomic analysis of a filamentous fungal endophyte using EST datasets.

Proteomic analysis of many species of fungi, particularly filamentous fungi, is difficult due to the lack of publicly available genome sequence data and the problems associated with cross-species comparisons. Furthermore, the detection of fungal proteins in biological systems where there are a greater number of proteins present from other eukaryote species provides additional challenges. We present an EST-based approach for identifying proteins from a fungal endophyte of temperate grasses and demonstrate that this method is well suited for fungi with minimal sequence data.

Molecular evolution of the Bovini tribe (Bovidae, Bovinae): is there evidence of rapid evolution or reduced selective constraint in Domestic cattle?

BACKGROUND: If mutation within the coding region of the genome is largely not adaptive, the ratio of nonsynonymous (dN) to synonymous substitutions (dS) per site (dN/dS) should be approximately equal among closely related species. Furthermore, dN/dS in divergence between species should be equivalent to dN/dS in polymorphisms. This hypothesis is of particular interest in closely related members of the Bovini tribe, because domestication has promoted rapid phenotypic divergence through strong artificial selection of some species while others remain undomesticated. We examined a number of genes that may be involved in milk production in Domestic cattle and a number of their wild relatives for evidence that domestication had affected molecular evolution. Elevated rates of dN/dS were further queried to determine if they were the result of positive selection, low effective population size (N(e)) or reduced selective constraint. RESULTS: We have found that the domestication process has contributed to higher dN/dS ratios in cattle, especially in the lineages leading to the Domestic cow (Bos taurus) and Mithan (Bos frontalis) and within some breeds of Domestic cow. However, the high rates of dN/dS polymorphism within B. taurus when compared to species divergence suggest that positive selection has not elevated evolutionary rates in these genes. Likewise, the low rate of dN/dS in Bison, which has undergone a recent population bottleneck, indicates a reduction in population size alone is not responsible for these observations. CONCLUSION: The effect of selection depends on effective population size and the selection coefficient (N(e)s). Typically under domestication both selection pressure for traits important in fitness in the wild and Ne are reduced. Therefore, reduced selective constraint could be responsible for the observed elevated evolutionary ratios in domesticated species, especially in B. taurus and B. frontalis, which have the highest dN/dS in the Bovini. This may have important implications for tests of selection such as the McDonald-Kreitman test. Surprisingly we have also detected a significant difference in the supposed neutral substitution rate between synonymous and noncoding sites in the Bovine genome, with a 30% higher rate of substitution at synonymous sites. This is due, at least in part, to an excess of the highly mutable CpG dinucleotides at synonymous sites, which will have implications for time of divergence estimates from molecular data.

Nucleotide-binding oligomerization domain containing 1 (NOD1) haplotypes and single nucleotide polymorphisms modify susceptibility to inflammatory bowel diseases in a New Zealand caucasian population: a case-control study.

BACKGROUND: The nucleotide-binding oligomerization domain containing 1 (NOD1) gene encodes a pattern recognition receptor that senses pathogens, leading to downstream responses characteristic of innate immunity. We investigated the role of NOD1 single nucleotide polymorphisms (SNPs) on IBD risk in a New Zealand Caucasian population, and studied Nod1 expression in response to bacterial invasion in the Caco2 cell line. FINDINGS: DNA samples from 388 Crohn's disease (CD), 405 ulcerative colitis (UC), 27 indeterminate colitis patients and 201 randomly selected controls, from Canterbury, New Zealand were screened for 3 common SNPs in NOD1, using the MassARRAY iPLEX Gold assay. Transcriptional activation of the protein produced by NOD1 (Nod1) was studied after infection of Caco2 cells with Escherichia coli LF82. Carrying the rs2075818 G allele decreased the risk of CD (OR = 0.66, 95% CI = 0.50-0.88, p < 0.002) but not UC. There was an increased frequency of the three SNP (rs2075818, rs2075822, rs2907748) haplotype, CTG (p = 0.004) and a decreased frequency of the GTG haplotype (p = 0.02).in CD. The rs2075822 CT or TT genotypes were at an increased frequency (genotype p value = 0.02), while the rs2907748 AA or AG genotypes showed decreased frequencies in UC (p = 0.04), but not in CD. Functional assays showed that Nod1 is produced 6 hours after bacterial invasion of the Caco2 cell line. CONCLUSION: The NOD1 gene is important in signalling invasion of colonic cells by pathogenic bacteria, indicative of its' key role in innate immunity. Carrying specific SNPs in this gene significantly modifies the risk of CD and/or UC in a New Zealand Caucasian population.

Single nucleotide polymorphism in the tumor necrosis factor-alpha gene affects inflammatory bowel diseases risk.

AIM: To investigate the role that single nucleotide polymorphisms (SNPs) in the promoter of the tumour necrosis factor-alpha (TNF-alpha) gene play in the risk of inflammatory bowel diseases (IBDs) in a New Zealand population, in the context of international studies. METHODS: DNA samples from 388 patients with Crohn's disease (CD), 405 ulcerative colitis (UC), 27 indeterminate colitis (IC) and 201 randomly selected controls, from Canterbury, New Zealand were screened for 3 common polymorphisms in the TNF-alpha receptor: -238 G-->A, -308 G-->A and -857C-->T, using a Taqman assay. A meta-analysis was performed on the data obtained on these polymorphisms combined with that from other published studies. RESULTS: Individuals carrying the -308 G/A allele had a significantly (OR = 1.91, c2 = 17.36, P < 0.0001) increased risk of pancolitis, and a 1.57-fold increased risk (OR = 1.57, c2 = 4.34, P = 0.037) of requiring a bowel resection in UC. Carrying the -857 C/T variant decreased the risk of ileocolonic CD (OR = 0.56, c2 = 4.32, P = 0.037), and the need for a bowel resection (OR = 0.59, c2 = 4.85, P = 0.028). The risk of UC was reduced in individuals who were smokers at diagnosis, (OR = 0.48, c2 = 4.86, P = 0.028). CONCLUSION: TNF-alpha is a key cytokine known to play a role in inflammatory response, and the locus for the gene is found in the IBD3 region on chromosome 6p21, known to be associated with an increased risk for IBD. The -308 G/A SNP in the TNF-alpha promoter is functional, and may account in part for the increased UC risk associated with the IBD3 genomic region. The -857 C/T SNP may decrease IBD risk in certain groups. Pharmaco- or nutrigenomic approaches may be desirable for individuals with such affected genotypes.

Association of DLG5 variants with inflammatory bowel disease in the New Zealand Caucasian population and meta-analysis of the DLG5 R30Q variant.

BACKGROUND: Variants in the DLG5 gene have been associated with inflammatory bowel disease (IBD) in samples from some, but not all populations. In particular, 2 nonsynonymous single-nucleotide polymorphisms (SNPs), R30Q (rs1248696) and P1371Q (rs2289310), have been associated with an increased risk of IBD, and a common haplotype (called haplotype "A") has been associated with reduced risk. METHODS: We genotyped R30Q, P1371Q, and a haplotype A tagging SNP (rs2289311) in a New Zealand Caucasian cohort of 389 Crohn's disease (CD) patients, 406 ulcerative colitis (UC) patients, and 416 population controls. Each SNP was tested for association with disease susceptibility and clinical phenotypes. We also performed a meta-analysis of R30Q data from published association studies. RESULTS: The haplotype A tagging SNP was associated with reduced risk of IBD at the 0.05 significance level (P=0.036) with an allelic odds ratio of 0.83 (95% confidence interval [CI]: 0.69-0.99). Association with haplotype A was strongest (odds ratio approximately 0.57) in UC patients with familial IBD or extraintestinal manifestations. The R30Q and P1371Q polymorphisms were not significantly associated with UC, CD, or IBD. Analysis of male and female data did not find any gender-specific associations. Meta-analysis gave no evidence of association of R30Q with IBD. CONCLUSIONS: Meta-analysis demonstrates that the minor allele of R30Q is not a risk factor for IBD across populations. This study provides some evidence that DLG5 haplotype A is associated with reduced risk of IBD in the New Zealand Caucasian population, but this association will need to be replicated in an independent sample.