Functional evaluation of genetic and environmental regulators of p450 mRNA levels

Variations in the activities of Cytochrome P450s are one of the major factors responsible for inter-individual differences in drug clearance rates, which may cause serious toxicity or inefficacy of therapeutic drugs. Various mRNA level is one of the key factors for different activity of the major P450 genes. Although both genetic and environmental regulators of P450 gene expression have been widely investigated, few studies have evaluated the functional importance of cis- and trans-regulatory factors and environmental factors in the modulation of inter-individual expression variations of the P450 genes. In this study, we measured the mRNA levels of seven major P450 genes (CYP1A1, CYP1A2, CYP2C9, CYP2C19, CYP2D6, CYP3A4 and CYP3A5) in 96 liver biopsy samples from Chinese population. Both trans-acting (mRNA levels and non-synonymous SNPs of putative regulator genes) and cis-acting (gene copy number and functional SNPs) factors were investigated to identify the determinants of the expression variations of these seven P450 genes. We found that expression variations of most P450 genes, regulator genes and housekeeping genes were positively correlated at the mRNA level. After partial correlation analysis using ACTB and GAPDH expression to eliminate the effect of global regulators, a UPGMA (Unweighted Pair Group Method with Arithmetic Mean) tree was constructed to reveal the effects of specific regulation networks potentially masked by global regulators. Combined with the functional analysis of regulators, our results suggested that expression variation at the mRNA level was mediated by several factors in a gene-specific manner. Cis-acting genetic variants might play key roles in the expression variation of CYP2D6 and CYP3A5, environmental inducers might play key roles in CYP1A1 and CYP1A2 variation and global regulators might play key roles in CYP2C9 variation. In addition, the functions of regulators that play less important roles in controlling expression variation for each P450 gene were determined.     

Survey of polymorphic sequence variation in the immediate 5' region of human DNA repair genes

Systematic screens have revealed extensive DNA sequence variation existing in the human population. Studies of the role of polymorphic genetic variants in explaining the association of family history with risk of common disease have generally focused on variants predicted to disrupt protein structure and activity. Recent studies have identified genetic variation in the level of expression of many genes, variation that is potentially biologically relevant in explaining individual variation in disease risk. In a survey of data available for 108 DNA repair genes that have been systematically screened for sequence variation, an average of 3.3 SNPs per gene were found to exist at a variant allele frequency of at least 0.02 in the region 2kb upstream from the 5'-untranslated region. One-third of the genes harbored a SNP with an allele frequency of at least 0.02 within a predicted promotor element. These variants are distributed among promoter elements that average 20 elements per gene. The frequency of polymorphic SNPs in CpG islands was 0.8 per gene, while the frequency of SNPs in the 5'-UTR was 0.7 per gene. The recognition of extensive genetic variation with potential to impact levels of gene expression, and thereby exacerbate the impact of amino acid substitution variants on the activity of proteins, increases the complexity of analyses required to explain the molecular genetic basis for the familial contribution to the sporadic incidence of common disease.     

Large-scale profiling and identification of potential regulatory mechanisms for allelic gene expression in colorectal cancer cells

Allelic variation in gene expression is common in humans and this variation is associated with phenotypic variation. In this study, we employed high-density single nucleotide polymorphism (SNP) chips containing 13,900 exonic SNPs to identify genes with allelic gene expression in cells from colorectal cancer cell lines. We found 2 monoallelically expressed genes (ERAP2 and MYLK4), 32 genes with an allelic imbalance in their expression, and 13 genes showing allele substitution by RNA editing. Among a total of 34 allelically expressed genes in colorectal cancer cells, 15 genes (44.1%) were associated with cis-acting eQTL, indicating that large portions of allelically expressed genes are regulated by cis-acting mechanisms of gene expression. In addition, potential regulatory variants present in the proximal promoter regions of genes showing either monoallelic expression or allelic imbalance were not tightly linked with coding SNPs, which were detected with allelic gene expression. These results suggest that multiple rare variants could be involved in the cis-acting regulatory mechanism of allelic gene expression. In the comparison with allelic gene expression data from Centre d'Etude du Polymorphisme Humain (CEPH) family B cells, 12 genes showed B-cell specific allelic imbalance and 1 noncoding SNP showed colorectal cancer cell-specific allelic imbalance. In addition, different patterns of allele substitution were observed between B cells and colorectal cancer cells. Overall, our study not only indicates that allelic gene expression is common in colorectal cancer cells, but our study also provides a better understanding of allele-specific gene expression in colorectal cancer cells.     

Genetics of sputum gene expression in chronic obstructive pulmonary disease

Previous expression quantitative trait loci (eQTL) studies have performed genetic association studies for gene expression, but most of these studies examined lymphoblastoid cell lines from non-diseased individuals. We examined the genetics of gene expression in a relevant disease tissue from chronic obstructive pulmonary disease (COPD) patients to identify functional effects of known susceptibility genes and to find novel disease genes. By combining gene expression profiling on induced sputum samples from 131 COPD cases from the ECLIPSE Study with genomewide single nucleotide polymorphism (SNP) data, we found 4315 significant cis-eQTL SNP-probe set associations (3309 unique SNPs). The 3309 SNPs were tested for association with COPD in a genomewide association study (GWAS) dataset, which included 2940 COPD cases and 1380 controls. Adjusting for 3309 tests (p<1.5e-5), the two SNPs which were significantly associated with COPD were located in two separate genes in a known COPD locus on chromosome 15: CHRNA5 and IREB2. Detailed analysis of chromosome 15 demonstrated additional eQTLs for IREB2 mapping to that gene. eQTL SNPs for CHRNA5 mapped to multiple linkage disequilibrium (LD) bins. The eQTLs for IREB2 and CHRNA5 were not in LD. Seventy-four additional eQTL SNPs were associated with COPD at p<0.01. These were genotyped in two COPD populations, finding replicated associations with a SNP in PSORS1C1, in the HLA-C region on chromosome 6. Integrative analysis of GWAS and gene expression data from relevant tissue from diseased subjects has located potential functional variants in two known COPD genes and has identified a novel COPD susceptibility locus.     

A survey of single nucleotide polymorphisms identified from whole‐genome sequencing and their functional effect in the porcine genome,

Genetic variants detected from sequence have been used to successfully identify causal variants and map complex traits in several organisms. High and moderate impact variants, those expected to alter or disrupt the protein coded by a gene and those that regulate protein production, likely have a more significant effect on phenotypic variation than do other types of genetic variants. Hence, a comprehensive list of these functional variants would be of considerable interest in swine genomic studies, particularly those targeting fertility and production traits. Whole‐genome sequence was obtained from 72 of the founders of an intensely phenotyped experimental swine herd at the U.S. Meat Animal Research Center (USMARC). These animals included all 24 of the founding boars (12 Duroc and 12 Landrace) and 48 Yorkshire–Landrace composite sows. Sequence reads were mapped to the Sscrofa10.2 genome build, resulting in a mean of 6.1 fold (×) coverage per genome. A total of 22 342 915 high confidence SNPs were identified from the sequenced genomes. These included 21 million previously reported SNPs and 79% of the 62 163 SNPs on the PorcineSNP60 BeadChip assay. Variation was detected in the coding sequence or untranslated regions (UTRs) of 87.8% of the genes in the porcine genome: loss‐of‐function variants were predicted in 504 genes, 10 202 genes contained nonsynonymous variants, 10 773 had variation in UTRs and 13 010 genes contained synonymous variants. Approximately 139 000 SNPs were classified as loss‐of‐function, nonsynonymous or regulatory, which suggests that over 99% of the variation detected in our pigs could potentially be ignored, allowing us to focus on a much smaller number of functional SNPs during future analyses.     

Novel and functional regulatory SNPs in the promoter region of <Emphasis Type="Italic">FOXP3</Emphasis> gene in a Gabonese population

Parasites exert a selection pressure on their hosts and are accountable for driving diversity within gene families and immune gene polymorphisms in a host population. The overwhelming response of regulatory T cells during infectious challenges directs the host immune system to lose the ability to mount parasite specific T cell responses. The underlying idea of this study is that regulatory single nucleotide polymorphism (SNPs) can cause significant changes in gene expression in functional immune genes. We identified and investigated regulatory SNPs in the promoter region of the FOXP3 gene in a group of Gabonese individuals exposed to a variety of parasitic infections. We identified two novel and one promoter variants in 40 individual subjects. We further validated these promoter variants for their allelic gene expression using transient transfection assays. Two promoter variants, −794 (C/G) and the other variant −734/−540 (C/T) revealed a significant higher expression of the reporter gene at basal level in comparison to the major allele. The identification of SNPs that modify function in the promoter regions could provide a basis for studying parasite susceptibility in association studies.     

Identification of common genetic variation that modulates alternative splicing

Alternative splicing of genes is an efficient means of generating variation in protein function. Several disease states have been associated with rare genetic variants that affect splicing patterns. Conversely, splicing efficiency of some genes is known to vary between individuals without apparent ill effects. What is not clear is whether commonly observed phenotypic variation in splicing patterns, and hence potential variation in protein function, is to a significant extent determined by naturally occurring DNA sequence variation and in particular by single nucleotide polymorphisms (SNPs). In this study, we surveyed the splicing patterns of 250 exons in 22 individuals who had been previously genotyped by the International HapMap Project. We identified 70 simple cassette exon alternative splicing events in our experimental system; for six of these, we detected consistent differences in splicing pattern between individuals, with a highly significant association between splice phenotype and neighbouring SNPs. Remarkably, for five out of six of these events, the strongest correlation was found with the SNP closest to the intron–exon boundary, although the distance between these SNPs and the intron–exon boundary ranged from 2 bp to greater than 1,000 bp. Two of these SNPs were further investigated using a minigene splicing system, and in each case the SNPs were found to exert cis-acting effects on exon splicing efficiency in vitro. The functional consequences of these SNPs could not be predicted using bioinformatic algorithms. Our findings suggest that phenotypic variation in splicing patterns is determined by the presence of SNPs within flanking introns or exons. Effects on splicing may represent an important mechanism by which SNPs influence gene function.     

Highly interconnected genes in disease-specific networks are enriched for disease-associated polymorphisms

Background

Complex diseases are associated with altered interactions between thousands of genes. We developed a novel method to identify and prioritize disease genes, which was generally applicable to complex diseases.

Results

We identified modules of highly interconnected genes in disease-specific networks derived from integrating gene-expression and protein interaction data. We examined if those modules were enriched for disease-associated SNPs, and could be used to find novel genes for functional studies. First, we analyzed publicly available gene expression microarray and genome-wide association study (GWAS) data from 13, highly diverse, complex diseases. In each disease, highly interconnected genes formed modules, which were significantly enriched for genes harboring disease-associated SNPs. To test if such modules could be used to find novel genes for functional studies, we repeated the analyses using our own gene expression microarray and GWAS data from seasonal allergic rhinitis. We identified a novel gene, FGF2, whose relevance was supported by functional studies using combined small interfering RNA-mediated knock-down and gene expression microarrays. The modules in the 13 complex diseases analyzed here tended to overlap and were enriched for pathways related to oncological, metabolic and inflammatory diseases. This suggested that this union of the modules would be associated with a general increase in susceptibility for complex diseases. Indeed, we found that this union was enriched with GWAS genes for 145 other complex diseases.

Conclusions

Modules of highly interconnected complex disease genes were enriched for disease-associated SNPs, and could be used to find novel genes for functional studies.     

SNP-based large-scale identification of allele-specific gene expression in human B cells

Polymorphism and variations in gene expression provide the genetic basis for human variation. Allelic variation of gene expression, in particular, may play a crucial role in phenotypic variation and disease susceptibility. To identify genes with allelic expression in human cells, we genotyped genomic DNA and cDNA isolated from 31 immortalized B cell lines from three Centre d'Etude du Polymorphisme Humain (CEPH) families using high-density single-nucleotide polymorphism (SNP) chips containing 13,900 exonic SNPs. We identified seven SNPs in five genes with monoallelic expression, 146 SNPs in 125 genes with allelic imbalance in expression with preferentially higher expression of one allele in a heterozygous individual. The monoallelically expressed genes (ERAP2, MDGA1, LOC644422, SDCCAG3P1 and CLTCL1) were regulated by cis-acting, non-imprinted differential allelic control. In addition, all monoallelic gene expression patterns and allelic imbalances in gene expression in B cells were transmitted from parents to offspring in the pedigree, indicating genetic transmission of allelic gene expression. Furthermore, frequent allele substitution, probably due to RNA editing, was also observed in 21 genes in 23 SNPs as well as in 48 SNPs located in regions containing no known genes. In this study, we demonstrated that allelic gene expression is frequently observed in human B cells, and SNP chips are very useful tools for detecting allelic gene expression. Overall, our data provide a valuable framework for better understanding allelic gene expression in human B cells.     

FCGR2C Polymorphisms Associated with HIV-1 Vaccine Protection Are Linked to Altered Gene Expression of Fc-γ Receptors in Human B Cells

The phase III Thai RV144 vaccine trial showed an estimated vaccine efficacy (VE) to prevent HIV-1 infection of 31.2%, which has motivated the search for immune correlates of vaccine protection. In a recent report, several single nucleotide polymorphisms (SNPs) in FCGR2C were identified to associate with the level of VE in the RV144 trial. To investigate the functional significance of these SNPs, we utilized a large scale B cell RNA sequencing database of 462 individuals from the 1000 Genomes Project to examine associations between FCGR2C SNPs and gene expression. We found that the FCGR2C SNPs that associated with vaccine efficacy in RV144 also strongly associated with the expression of FCGR2A/C and one of them also associated with the expression of Fc receptor-like A (FCRLA), another Fc-γ receptor (FcγR) gene that was not examined in the previous report. These results suggest that the expression of FcγR genes is influenced by these SNPs either directly or in linkage with other causal variants. More importantly, these results motivate further investigations into the potential for a causal association of expression and alternative splicing of FCGR2C and other FcγR genes with the HIV-1 vaccine protection in the RV144 trial and other similar studies.