Upstream Single Nucleotide Polymorphisms Research Articles

Identification of Salmonella Pullorum Factors Affecting Immune Reaction in Macrophages from the Avian Host.

The host-specific Salmonella serovar S. Pullorum (SP) modulates the chicken immune response to a Th2-biased response associated with persistent infection. This is different from the Th1-biased immune response induced by the genetically close serovar, S. Enteritidis (SE). Based on core genome differences between SP and SE, we used three complementary bioinformatics approaches to identify SP genes, which may be important for stimulation of the immune response. Defined mutants were constructed in selected genes, and the infection potential and ability of mutants to stimulate cytokine production in avian derived HD11 macrophages were determined. Deletion of large genomic regions unique to SP did not change infection potential nor immune stimulation significantly. Mutants in genes with conserved single nucleotide polymorphisms (SNPs) between the two serovars in the region 100 bp upstream of the start codon (conserved upstream SNPs [CuSNPs]) such as sseE, osmB, tolQ, a putative immune antigen, and a putative persistent infection factor, exhibited differences in induction of inflammatory cytokines compared to wild-type SP, suggesting a possible role of these CuSNPs in immune regulation. Single nucleotide SP mutants correcting for the CuSNP difference were constructed in the upstream region of sifA and pipA. The SNP corrected pipA mutant expressed pipA at a higher level than the wild-type SP strain, and the mutant differentially caused upregulation of proinflammatory cytokines. It suggests that this CuSNP is important for the suppression of proinflammatory responses. In conclusion, this study has identified putative immune stimulating factors of relevance to the difference in infection dynamics between SP and SE in avian macrophages. IMPORTANCE Salmonella Pullorum is host specific to avian species, where it causes life-threatening infection in young birds. It is unknown why it is host restricted and causes systemic disease, rather than gastroenteritis normally seen with Salmonella. In the present study, we identified genes and single nucleotide polymorphisms (SNPs; relative to the broad-host-range type Salmonella Enteritidis), which affected survival and immune induction in macrophages from hens suggesting a role in development of the host specific infection. Further studies of such genes may enable understanding of which genetic factors determine the development of host specific infection by S. Pullorum. In this study, we developed an in silico approach to predict candidate genes and SNPs for development of the host-specific infection and the specific induction of immunity associated with this infection. This study flow can be used in similar studies in other clades of bacteria.

Analysis of Association Between MGMT and p53 Gene Single Nucleotide Polymorphisms and Laryngeal Cancer.

To investigate the p53 and O6-methylguanine DNA methyltransferase (MGMT)5' upstream sequence gene promoter regions for single nucleotide polymorphisms and explore the p53 gene 5' upstream sequence consisting of two haplotypes to provide a genetic marker for the incidence of laryngeal squamous cell carcinoma. We included 96 cases of laryngeal squamous cell carcinoma and 102 controls. We used SNaPshot micro-sequencing analysis of the MGMT promoter region for four single nucleotide polymorphisms and p53 gene 5' upstream sequence loci (rs1625649, rs2287499, rs2287498, rs228749) genotypes. We calculated and compared two groups for genotypic and allelic frequencies, applied HaploView4.2 for computing rs2287499, rs2287498, rs228749 values and haplotype frequencies and tested control loci and Hardy-Weinberg equilibrium. All the experimental data were statistically evaluated using SPSS17.0. The Chi-square test was used for statistical analysis with p<0.05 indicating statistical significance. 5'Upstream single nucleotide polymorphisms rs1625649, rs2287499, rs2287498, rs228749 of p53 were polymorphic in both patient and control groups. There was no statistical significance in frequency distributions for the four loci genotypes when comparing patients and healthy controls (Chi-square values were 4.47, 0.98, 1.67, 4.68, respectively; p>0.05). However, allelic frequencies of the MGMT promoter region locus rs1625649 between patients and healthy control groups were statistically significantly different (chi-square value of 5.77; p<0.05). Differences between allelic frequencies for the p53 gene 5' upstream sequence loci rs2287499, rs2287498 and rs228749 between patients and the healthy control group were not statistically significant (Chi-square values were 1.11,1.56,3.36; p>0.05). Nor were those for the two haplotypes of rs2287499, rs2287498 and rs228749 between patients and the healthy control group were not statistically significant (Chi-square value 1.46, p>0.05). MGMT gene polymorphism appears to be associated with the incidence of laryngeal cancer.

In-vitro function of upstream visfatin polymorphisms that are associated with adverse cardiometabolic parameters in obese children.

BackgroundVisfatin is an adipokine associated with glucose and lipid metabolism. We previously reported two visfatin upstream single nucleotide polymorphisms (SNPs), c.-3187G > A (rs11977021) and c.-1537C > T (rs61330082), which were in perfect linkage disequilibrium, in a Singaporean cohort of severely obese children and are associated with visfatin level and adverse cardiometabolic parameters. We aim to functionally characterize the effect of c.-3187G > A and c.-1537C > T SNPs on basal transcriptional activity.MethodsA 1.6 kb and 3.7 kb upstream promoter region of the visfatin gene was amplified by polymerase chain reaction and separately cloned into luciferase reporter vectors. Successful clones were transfected into human embryonic kidney (HEK293T) and human breast carcinoma (MCF7) cells and in-vitro dual-luciferase assay was performed. Electrophoretic mobility shift assay (EMSA) was also conducted to examine the binding affinity between transcription factors and visfatin promoter sequences.ResultsVariant promoter with only c.-1537C > T SNP did not show a change in transcriptional activity as compared to the wild type. However, variant promoter with both c.-3187G > A and c.-1537C > T SNPs showed a statistically significant increase of 1.41 fold (p < 0.01) in transcriptional activity. The longer 3.7kbp visfatin promoter sequence was also shown to have significantly higher transcriptional activity (p < 0.05) as compared to the shorter 1.6kbp visfatin promoter. Both c.-3187G > A and c.-1537C > T variants showed an increased binding with nuclear protein.Discussion and conclusionsWe have demonstrated for the first time that visfatin variant promoter with both c.-3187G > A and c.-1537C > T SNPs result in an increase in transcriptional activity. This supports our previous finding and postulation that these SNPs contribute to elevated visfatin levels which may mediate higher triglyceride levels, severe systolic blood pressure and severe hypertension in obese children. These SNPs may co-operatively affect enhancer or silencer function to regulate transcriptional activity. In conclusion, this study shows that upstream visfatin SNPs could potentially affect phenotypic outcome in obese children through alteration of circulating visfatin level.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-016-3315-9) contains supplementary material, which is available to authorized users.

Variants in the Regulatory Region of WNT5A Reduced Risk of Cardiac Conotruncal Malformations in the Chinese Population.

WNT5A is one of the most highly investigated non-canonical Wnt ligands and is involved in the embryonic heart development, especially in formation of the cardiac conotruncal region by regulating the migration and differentiation of cardiac neural crest (CNC) and second heart field (SHF) cells. No study to date has comprehensively characterized the WNT5A regulatory variants in patients with congenital heart malformations (CHMs). The association between regulatory variants of the WNT5A gene and CHMs was examined in case-control association study in 1,210 CHMs and 798 controls. Individuals carrying a homozygous genotype CC (rs524153) or GG (rs504849) had a similarly reduced risk of conotruncal malformations. The homozygous genotypes (CC for rs524153 and GG for rs504849) were associated with a lower WNT5A transcriptional level compared with the transcriptional level of those with wild-type genotypes. Further functional analysis revealed that an additional upstream single nucleotide polymorphisms (SNP) rs371954924 (–5244GCCA > CC) in a linkage disequilibrium (LD) block with the above genotyped SNPs decreased WNT5A expression through the attenuated binding affinity with the transcription factor SOX9. This is the first demonstration that genetic variants in the regulatory regions of WNT5A play a vital role in sporadic conotruncal malformations susceptibility through the changeable expression of the WNT5A gene.

Epigenetics in health and disease: heralding the EWAS era

Epigenetics in health and disease: heralding the EWAS era

Visfatin and its genetic variants are associated with obesity‐related morbidities and cardiometabolic risk in severely obese children

Visfatin is an adipokine, associated with obesity and possibly glucose regulation. The aim of this study was to examine the association of visfatin and its genetic variants with adiposity, cardiometabolic risk factors and obesity-related morbidities in obese children. Anthropometric measurements, dual energy X-ray absorptiometry scan, fasting blood samples and oral glucose tolerance tests were performed for 243 obese children. We screened the visfatin gene of 24 obese subjects and then performed genotyping of identified genetic variants in other 219 obese children through direct DNA sequencing. Fasting serum visfatin correlated with measures of obesity and liver enzymes and was elevated in obese children with abnormal glucose tolerance and non-alcoholic fatty liver disease. The two upstream single nucleotide polymorphisms, -3187G>A (rs11977021) and -1537C>T (rs61330082), were at complete linkage disequilibrium. The AA genotype of -3187G>A was associated with higher serum visfatin (6.17 ± 0.76 ng mL(-1) vs. 3.92 ± 0.44 ng mL(-1)) and higher triglyceride (1.39 ± 0.08 mmol L(-1) vs. 1.19 ± 0.07 mmol L(-1)) as compared with the GG genotype. There was also a significant linear increase in serum visfatin across GG to GA to AA genotype of -3187G>A, indicating possible additive effect of A allele. The dominant GA + AA genotype model of +21426G>A (rs2302559) was associated with lower serum visfatin (3.83 ± 0.56 ng mL(-1) vs. 5.13 ± 0.34 ng mL(-1)) and lower plasma glucose (4.37 ± 0.08 mmol L(-1) vs. 4.77 ± 0.12 mmol L(-1)) as compared with the GG genotype. Visfatin and its genetic variants were associated with adiposity, obesity-related morbidities and adverse cardiometabolic parameters. This supported our hypothesis that visfatin plays a significant role in the development of obesity-related morbidities and cardiometabolic risk.

Targeting the association of calgranulin B (S100A9) with insulin resistance and type 2 diabetes

Calgranulin B (S100A9) was recognized as a candidate type 2 diabetes (T2D) gene in the genomic profiling of muscle from a rodent model of T2D and identifying the human orthologs of genes localized in T2D susceptibility regions. Circulating and S100A9 expressions in muscle and adipose tissue, isolated fat cells, and mouse models were evaluated. A common 5'-upstream single-nucleotide polymorphism (SNP; rs3014866) for S100A9 was analyzed, as well as the effects of weight loss and treatments in vitro with recombinant S100A9. S100a9 expression was increased in muscle of diabetic mice (1.6-fold, p = 0.002), and in muscle from subjects with impaired glucose tolerance (∼4-fold, p = 0.028; n = 34). The rs3014866 SNP was associated with circulating S100A9 and the risk of T2D, having TT carriers at 28 % (p = 0.03) lower risk (n = 1,450). Indeed, increased circulating S100A9 (∼4-fold, p = 0.03; n = 206) and subcutaneous (2-fold, p = 0.01) and omental (1.4-fold, p = 0.04) S100A9 gene expressions (n = 83) in TT carriers run in parallel to decreased fasting glucose and glycated hemoglobin. Accordingly, metformin led to increased S100A9 mRNA in ex vivo-treated adipose tissue explants (n = 5/treatment). Otherwise, obese subjects showed a compensatory increase in circulating and S100A9 expressions in adipose (n = 126), as further demonstrated by decreased levels after diet- (-34 %, p = 0.002; n = 20) and surgery-induced (-58 %, p = 0.02; n = 8) weight loss. Lipopolysaccharide led to increased S100A9 in adipose from mice (n = 5/treatment) while recombinant S100A9 downregulated inflammation in adipocytes (n = 3/treatment). Current findings support the strategy of testing differentially expressed genes in mice and human orthologs associated with T2D. The increased S100A9 reported for obesity and insulin resistance may be envisioned as a compensatory mechanism for inflammation.

GDF5 single-nucleotide polymorphism rs143383 is associated with lumbar disc degeneration in Northern European women

ObjectiveLumbar disc degeneration (LDD) is a serious social and medical problem which has been shown to be highly heritable. It has similarities with peripheral joint osteoarthritis (OA) in terms of both epidemiology and pathologic processes. A few known genetic variants have been identified using a candidate gene approach, but many more are thought to exist. GDF5 is a gene whose variants have been shown to play a role in skeletal height as well as predisposing to peripheral joint OA. In vitro, the gene product growth differentiation factor 5 has been shown to promote growth and repair of animal disc. This study was undertaken to investigate whether the GDF5 gene plays a role in LDD.MethodsWe investigated whether the 5′ upstream single-nucleotide polymorphism (SNP) variant rs143383 was associated with LDD, using plain radiography and magnetic resonance imaging to identify disc space narrowing and osteophytes, in 5 population cohorts from Northern Europe.ResultsAn association between LDD and the SNP rs143383 was identified in women, with the same risk allele as in knee and hip OA (odds ratio 1.72 [95% confidence interval 1.15–2.57], P = 0.008).ConclusionOur findings in 5 population cohorts from Northern Europe indicate that a variant in the GDF5 gene is a risk factor for LDD in women. Many more such variants are predicted to exist, but this result highlights the growth and differentiation cellular pathway as a possible route to a better understanding of the process behind lumbar disc degeneration.

A regulatory SNP in AKAP13 is associated with blood pressure in Koreans

High blood pressure contributes to more than 10 million deaths per year worldwide through stroke and ischemic heart disease. Yet, genome-wide association studies (GWASs) have identified a small fraction of its underlying genetic factors. To identify biologically important single-nucleotide polymorphisms (SNPs) that regulate variations in blood pressure, we analyzed SNPs in a genome-wide association study. Genome-wide genotype data (original study n = 7551, SNP = 352,228; replication study n = 3703, SNP = 20) were obtained from the Korea National Institute of Health, wherein 29,921 of 352,228 SNPs lay within 5 kbp upstream of genes. Linear regression analysis was performed for systolic and diastolic blood pressure (DBP) by controlling for cohort, age, sex and body mass index. For the 20 SNPs that were associated with both blood pressure values, a replication study was performed in an independent population. A total of 20 SNPs were significantly associated with both blood pressure values in the original study, 13 of which lay in a conserved transcription factor-binding site. One SNP (rs11638762), in the GATA-3 binding site upstream of the AKAP13 gene, was significantly replicated in another cohort (P-value of the meta-analysis = 1.4 × 10(-5) for systolic blood pressure and 6.3 × 10(-4) for DBP). A functional GWAS was performed using upstream SNPs, and a novel genetic factor (AKAP13), which is essential for cardiac myocyte development in mice, was identified as a regulator of blood pressure.

Identification of Novel SNPs by Next-Generation Sequencing of the Genomic Region Containing theAPCGene in Colorectal Cancer Patients in China

We described an approach of identifying single nucleotide polymorphisms (SNPs) in complete genomic regions of key genes including promoters, exons, introns, and downstream sequences by combining long-range polymerase chain reaction (PCR) or NimbleGen sequence capture with next-generation sequencing. Using the adenomatous polyposis coli (APC) gene as an example, we identified 210 highly reliable SNPs by next-generation sequencing analysis program MAQ and Samtools, of which 69 were novel ones, in the 123-kb APC genomic region in 27 pair of colorectal cancers and normal adjacent tissues. We confirmed all of the eight randomly selected high-quality SNPs by allele-specific PCR, suggesting that our false discovery rate is negligible. We identified 11 SNPs in the exonic region, including one novel SNP that was not previously reported. Although 10 of them are synonymous, they were predicted to affect splicing by creating or removing exonic splicing enhancers or exonic splicing silencers. We also identified seven SNPs in the upstream region of the APC gene, three of which were only identified in the cancer tissues. Six of these upstream SNPs were predicted to affect transcription factor binding. We also observed that long-range PCR was better in capturing GC-rich regions than the NimbleGen sequence capture technique.

Intron‐1 rs3761548 is related to the defective transcription of Foxp3 in psoriasis through abrogating E47/c‐Myb binding

Foxp3 is a master transcription factor (TF) for development and function of CD4(+)CD25(+)Foxp3(+) regulatory T cells (Treg cells) and is critical for the transcription of target genes. But the transcriptional regulation of Foxp3 itself has not been fully understood until now. Here, we aimed to demonstrate the hypothesis that upstream single nucleotide polymorphism(s) (SNPs) of Foxp3 was/were responsible for the defective transcription of Foxp3 in psoriasis and to explore the mechanism behind this hypothesis. In this study, SNP of large sample was investigated for risk analysis. Mature algorithms, electrophoretic mobility shift and chromatin immunoprecipitation assays were used to identify TF binding site variations. Loss-of-function and overexpression assays and cell cycle blocker assay were performed to identify when and what kind of possible roles the candidate factors play. Our results showed that intron-1 rs3761548 was correlated with a significant susceptibility to psoriasis. The rs3761548 contributed to the decreased resting Foxp3 transcription and impaired acceleration of Foxp3 transcription levels after stimulation in psoriatic patients with genotype AA. We analysed and demonstrated potent new E47/c-Myb -dependent regulation elements in rs3761548, oppositely controlling Foxp3 gene transcription at G1 and G2/M phases of Treg cells in psoriatic patients. For patients with rs3761548 AA, the polymorphism causes loss of bindings to the E47 and c-Myb factors, leading to defective transcription of Foxp3 gene. Further identification of the networks and molecular mechanisms underlying Foxp3 transcription may provide new insights into Foxp3 transcriptional regulation and alternative therapeutic strategies to improve characteristics of autoimmune disorders.

Interpopulation Variation Frequency of Human Inosine 5′-Monophosphate Dehydrogenase Type II (IMPDH2) Genetic Polymorphisms

Inosine monophosphate dehydrogenase type II (IMPDH2) is the target for immunosuppression by mycophenolic acid and has been linked to resistance of tumor cells to chemotherapy. Determining the frequency of IMPDH2 genetic polymorphisms can inform the design of clinical studies investigating the impact of IMPDH2 genetic variability on both cancer therapy and immunosuppression. Frequencies of three IMPDH2 polymorphisms (rs4974081, rs5848860, and rs11557540) in >400 DNA samples from four different racial/ethnic groups (Caucasian, African American, Hispanic, and Asian populations) were characterized by the pyrosequencing genotyping method. For rs5848860 1591 CTT/- (500 G/EG) and rs11557540 1345 A/G (418 D/G), we did not observe any variant alleles in all DNA samples from these populations, which suggests that these variants could simply be sequencing errors rather than real polymorphisms. The observed frequency of the 5'-upstream single-nucleotide polymorphism (SNP) rs4974081 A/G was similar to that previously reported in the NCBI databank (dbSNP). An in silico functional analysis using FASTSNP predicts that this promoter SNP rs4974081 (-3624 A/G) could be a potential transcription factor binding site. This finding suggests that rs4974081 could be a good candidate SNP for association studies with immunosuppressive and chemotherapeutic therapy outcomes.

Genetic association between 5′-upstream single-nucleotide polymorphisms of PDGFRB and schizophrenia in a Korean population

PDGFRB is located on chromosome 5q31–q32, a chromosomal region identified by linkage analyses to contain a susceptibility gene for schizophrenia (SCZ). Recent research has focused on the role of the N-methyl-d-aspartate (NMDA) receptor in the pathogenesis of SCZ. D4 dopamine receptor-mediated transactivation of the gene encoding platelet-derived growth factor receptor beta (PDGFRB) has immediate effects on synaptic neurotransmission via calcium-dependent inactivation of NMDA receptors. In this study, we investigate the association between the PDGFRB gene and SCZ in a Korean population. We screened 6 single-nucleotide polymorphisms (SNPs) in the 5′-upstream region of PDGFRB and conducted a case–control study of 381 SCZ patients and 752 controls. The genotype and haplotype frequencies of 3 of the 6 SNPs [SNP1 (g.−1924T>C, rs3756314), SNP3 (g.−1772A>G, rs3756312) and SNP4 (rs3756311, g.−1658G>A)] were significantly associated with SCZ [SNP1, corrected p=0.012 (co-dominant model), 0.002 (Dominant model), and 0.506 (Recessive model); SNP3 and 4, corrected p=0.003, 0.009, and 0.049]. Haplotype analysis also revealed that ht1 (CGG) and ht2 (TAA) were significantly associated with SCZ (ht1, corrected p=0.018, 0.340, and 0.010; ht2, corrected p=0.002, 0.009, and 0.016). Transient transfection in neuronal cells revealed that ht1 had higher luciferase activity than the vector alone. Furthermore, Pdgfrb expression was increased in the frontal cortex and hippocampus in a mouse model of SCZ induced by MK801. We conclude that SNPs of the 5′-upstream region of PDGFRB are associated with SCZ in a Korean population. These are weak positives that require future studies to confirm these results.

A Variant in XPNPEP2 Is Associated with Angioedema Induced by Angiotensin I–Converting Enzyme Inhibitors

Angiotensin I-converting enzyme inhibitors (ACEi), which are used to treat common cardiovascular diseases, are associated with a potentially life-threatening adverse reaction known as angioedema (AE-ACEi). We have previously documented a significant association between AE-ACEi and low plasma aminopeptidase P (APP) activity. With eight large pedigrees, we hereby demonstrate that this quantitative trait is partially regulated by genetic factors. We tested APP activity using a variance-component QTL analysis of a 10-cM genomewide microsatellite scan enriched with seven markers over two candidate regions. We found significant linkage (LOD = 3.75) to a locus that includes the XPNPEP2 candidate gene encoding membrane-bound APP. Mutation screening of this QTL identified a large coding deletion segregating in one pedigree and an upstream single-nucleotide polymorphism (C-2399A SNP), which segregates in the remaining seven pedigrees. Measured genotype analysis strongly suggests that the linkage signal for APP activity at this locus is accounted for predominantly by the SNP association. In a separate case-control study (20 cases and 60 controls), we found significant association of this SNP to ACEi-induced AE (P=.0364). In conclusion, our findings provide supporting evidence that the C-2399A variant in XPNPEP2 is associated with reduced APP activity and a higher incidence of AE-ACEi.

Haplotype‐controlled analysis of the association of a non‐synonymous single nucleotide polymorphism at DBH (+ 1603C → T) with plasma dopamine β‐hydroxylase activity

The DBH locus controls plasma dopamine beta-hydroxylase activity (pDbetaH). A 5'-upstream single nucleotide polymorphism (SNP) at DBH (-1021C --> T) explains approximately 45% of the variance in pDbetaH, and a non-synonymous SNP in exon 11 (+ 1603C --> T) an additional 2%. However, that regression result underestimates the effect of + 1603C --> T because of its low minor allele frequency. We estimated the biological effect of + 1603C --> T on pDbetaH by comparing subjects of identical -1021CgammaT genotype, in a diagnostically heterogeneous group of subjects of European origin (N = 367). + 1603C --> T genotype associated with pDbetaH within groups of identical genotype at -1021 C --> T, accounting for 5%-16% of the variance. There was no significant linkage disequilibrium between -1021C --> T and + 1603C --> T (D = 0.0058, D' = 0.4774, d(2) = 0.0011, P > 0.05), confirming the validity of assessing the two polymorphisms independently. These results suggest that altered homospecific activity of the enzyme can contribute to variation in pDbetaH. This conclusion informs how associations between DBH and psychiatric disorders should be approached. (c) 2005 Wiley-Liss, Inc.

Discovery of Novel Flavin-Containing Monooxygenase 3 (FMO3) Single Nucleotide Polymorphisms and Functional Analysis of Upstream Haplotype Variants

The flavin-containing monooxygenases (FMOs) are important for xenobiotic metabolism. FMO3, the predominant FMO enzyme in human adult liver, exhibits significant interindividual variation that is poorly understood. This study was designed to identify common FMO3 genetic variants and determine their potential for contributing to interindividual differences in FMO3 expression. FMO3 single nucleotide polymorphism (SNP) discovery was accomplished by resequencing DNA samples from the Coriell Polymorphism Discovery Resource. Population-specific SNP frequencies were determined by multiplexed, single-base extension using DNA from 201 Hispanic American (Mexican descent), 201 African American, and 200 White (northern European descent) subjects. Haplotypes were inferred and population frequencies estimated using PHASE version 2.1. Multiple site-directed mutagenesis was used to introduce inferred upstream haplotypes into an FMO3/luciferase construct for functional analysis in HepG2 cells. Sequence analysis revealed seven FMO3 upstream SNPs, 11 exon SNPs, and 22 intron SNPs. Five of the latter fell within consensus splice sites. A g.72G>T variant (E24D) is predicted to impact the structure of the Rossmann fold involved in FAD binding, whereas a g.11177C>A variant (N61K) is predicted to disrupt the secondary structure of a conserved membrane interaction domain. Seven common (>1%) promoter region haplotypes were inferred in one or more of the study populations that differed in estimated frequency among the groups. Haplotype 2 resulted in an 8-fold increase in promoter activity, whereas haplotypes 8 and 15 exhibited a near complete loss of activity. In conclusion, FMO3 promoter haplotype variants modulate gene function and probably contribute to interindividual differences in FMO3 expression.