Clinical, structural and functional implications of mutations and polymorphisms in human NADPH P450 oxidoreductase

Cytochrome P450 proteins are involved in metabolism of drugs and xenobiotics. In the endoplasmic reticulum a single nicotinamide adenine dinucleotide phosphate (NADPH) P450 oxidoreductase (POR) supplies electrons to all microsomal P450s for catalytic activity. POR is a flavoprotein that contains both flavin mononucleotide and flavin adenine dinucleotide as cofactors and uses NADPH as the source of electrons. We have recently reported a number of POR mutations in the patients with disordered steroidogenesis. In the first report we had described missense mutations (A287P, R457H, V492E, C569Y, and V608F) identified in four patients with defects in steroid production. Each POR variant was produced as recombinant N-27 form of the enzyme in bacteria and as full-length form in yeast. Membranes from bacteria or yeast expressing normal or variant POR were purified and their activities were characterized in cytochrome c and CYP17A1 assays. Later we have published a larger study that described a whole range of POR mutations and characterized the mutants/polymorphisms A115V, T142A, M263V, Y459H, A503V, G539R, L565P, R616X, V631I, and F646del from the sequencing of patient DNA. We also studied POR variants Y181D, P228L, R316W, G413S, and G504R that were available in public databases or published literature. Three-dimensional structure of rat POR is known and we have used this structure to deduce the structure-function correlation of POR mutations in human. The missense mutations found in patients with disordered steroidogenesis are generally in the co-factor binding and functionally important domains of POR and the apparent polymorphisms are found in regions with lesser structural importance. A variation in POR can alter the activity of all microsomal P450s, and therefore, can affect the metabolism of drugs and xenobiotics even when the P450s involved are otherwise normal. It is important to study the genetic and biochemical basis of POR variants in human population to gain information about possible differences in P450 mediated reactions among the individuals carrying a variant or polymorphic form of POR that could impact their metabolism.

Systemic analysis of microarray data sets towards identifying genes regulating root development

Nitrogen and water are essential for plant growth and development. In this study, we designed experiments to produce gene expression data of poplar roots under nitrogen starvation and water deprivation conditions. We found low concentration of nitrogen led first to increased root elongation followed by lateral root proliferation and eventually increased root biomass. To identify genes regulating root growth and development under nitrogen starvation and water deprivation, we designed a series of data analysis procedures, through which, we have successfully identified biologically important genes. Differentially Expressed Genes (DEGs) analysis identified the genes that are differentially expressed under nitrogen starvation or drought. Protein domain enrichment analysis identified enriched themes (in same domains) that are highly interactive during the treatment. Gene Ontology (GO) enrichment analysis allowed us to identify biological process changed during nitrogen starvation. Based on the above analyses, we examined the local Gene Regulatory Network (GRN) and identified a number of transcription factors. After testing, one of them is a high hierarchically ranked transcription factor that affects root growth under nitrogen starvation. It is very tedious and time-consuming to analyze gene expression data. To avoid doing analysis manually, we attempt to automate a computational pipeline that now can be used for identification of DEGs and protein domain analysis in a single run. It is implemented in scripts of Perl and R.

Mutations and polymorphisms in the human argininosuccinate lyase (ASL) gene

Argininosuccinate lyase deficiency (ASLD) is caused by a defect of the urea cycle enzyme argininosuccinate lyase (ASL) encoded by the ASL gene. Patients often present early after birth with hyperammonemia but can also manifest outside the neonatal period mainly triggered by excessive protein catabolism. Clinical courses comprise asymptomatic individuals who only excrete the biochemical marker, argininosuccinic acid, in urine, and patients who succumb to their first hyperammonemic decompensation. Some patients without any hyperammonemia develop severe neurological disease. Here, we are providing an update on the molecular basis of ASLD by collecting all published (n = 67) as well as novel mutations (n = 67) of the ASL gene. We compile data on all 160 different genotypes ever identified in 223 ASLD patients, including clinical courses whenever available. Finally, we are presenting structural considerations focusing on the relevance of mutations for ASL homotetramer formation. ASLD can be considered as a panethnic disease with only single founder mutations identified in the Finnish (c.299T>C, p.Ile100Thr) and Arab (c.1060C>T, p.Gln354*) population. Most mutations are private with only few genotypes recurring in unrelated patients. The majority of mutations are missense changes including some with more frequent occurrence such as p.Arg12Gln, p.Ile100Thr, p.Val178Met, p.Arg186Trp, p.Glu189Gly, p.Gln286Arg, and p.Arg385Cys.

A multiplex PCR for the simultaneous detection and genotyping of the Echinococcus granulosus complex.

Echinococcus granulosus is characterized by high intra-specific variability (genotypes G1-G10) and according to the new molecular phylogeny of the genus Echinococcus, the E. granulosus complex has been divided into E. granulosus sensu stricto (G1-G3), E. equinus (G4), E. ortleppi (G5), and E. canadensis (G6-G10). The molecular characterization of E. granulosus isolates is fundamental to understand the spatio-temporal epidemiology of this complex in many endemic areas with the simultaneous occurrence of different Echinococcus species and genotypes. To simplify the genotyping of the E. granulosus complex we developed a single-tube multiplex PCR (mPCR) allowing three levels of discrimination: (i) Echinococcus genus, (ii) E. granulosus complex in common, and (iii) the specific genotype within the E. granulosus complex. The methodology was established with known DNA samples of the different strains/genotypes, confirmed on 42 already genotyped samples (Spain: 22 and Bulgaria: 20) and then successfully applied on 153 unknown samples (Tunisia: 114, Algeria: 26 and Argentina: 13). The sensitivity threshold of the mPCR was found to be 5 ng Echinoccoccus DNA in a mixture of up to 1 µg of foreign DNA and the specificity was 100% when template DNA from closely related members of the genus Taenia was used. Additionally to DNA samples, the mPCR can be carried out directly on boiled hydatid fluid or on alkaline-lysed frozen or fixed protoscoleces, thus avoiding classical DNA extractions. However, when using Echinococcus eggs obtained from fecal samples of infected dogs, the sensitivity of the mPCR was low (<40%). Thus, except for copro analysis, the mPCR described here has a high potential for a worldwide application in large-scale molecular epidemiological studies on the Echinococcus genus.

Polymorphisms of the DNA repair gene MGMT and risk and progression of head and neck cancer.

Methylating agents are involved in carcinogenesis, and the DNA repair protein O(6)-methylguanine-DNA methyltransferase (MGMT) removes methyl group from O(6)-methylguanine. Genetic variation in DNA repair genes has been shown to contribute to susceptibility to squamous cell carcinoma of the head and neck (SCCHN). We hypothesize that MGMT polymorphisms are associated with risk of SCCHN. In a hospital-based case-control study of 721 patients with SCCHN and 1234 cancer-free controls frequency-matched by age, sex and ethnicity, we genotyped four MGMT polymorphisms, two in exon 3, 16195C>T and 16286C>T and two in the promoter region, 45996G>T and 46346C>A. We found that none of these polymorphisms alone had a significant effect on risk of SCCHN. However, when these four polymorphisms were evaluated together by the number of putative risk genotypes (i.e. 16195CC, 16286CC, 45996GT+TT, and 46346CA+AA), a statistically significantly increased risk of SCCHN was associated with the combined genotypes with three to four risk genotypes, compared with those with zero to two risk genotypes (adjusted odds ratio (OR)=1.27; 95% confidence interval (CI)=1.05-1.53). This increased risk was also more pronounced among young subjects (OR=1.81; 95% CI=1.11-2.96), men (OR=1.24; 95% CI=1.00-1.55), ever smokers (OR=1.25; 95%=1.01-1.56), ever drinkers (OR=1.29; 95% CI=1.04-1.60), patients with oropharyngeal cancer (OR=1.45; 95% CI=1.12-1.87), and oropharyngeal cancer with regional lymph node metastasis (OR=1.52; 95% CI=1.16-1.89). In conclusion, our results suggest that any one of MGMT variants may not have a substantial effect on SCCHN risk, but a joint effect of several MGMT variants may contribute to risk and progression of SCCHN, particularly for oropharyngeal cancer, in non-Hispanic whites.

Sequencing of TGF-beta pathway genes in familial cases of intracranial aneurysm.

BACKGROUND AND PURPOSE: Familial aggregation of intracranial aneurysms (IA) strongly suggests a genetic contribution to pathogenesis. However, genetic risk factors have yet to be defined. For families affected by aortic aneurysms, specific gene variants have been identified, many affecting the receptors to transforming growth factor-beta (TGF-beta). In recent work, we found that aortic and intracranial aneurysms may share a common genetic basis in some families. We hypothesized, therefore, that mutations in TGF-beta receptors might also play a role in IA pathogenesis. METHODS: To identify genetic variants in TGF-beta and its receptors, TGFB1, TGFBR1, TGFBR2, ACVR1, TGFBR3, and ENG were directly sequenced in 44 unrelated patients with familial IA. Novel variants were confirmed by restriction digestion analyses, and allele frequencies were analyzed in cases versus individuals without known intracranial disease. Similarly, allele frequencies of a subset of known SNPs in each gene were also analyzed for association with IA. RESULTS: No mutations were found in TGFB1, TGFBR1, TGFBR2, or ACVR1. Novel variants identified in ENG (p.A60E) and TGFBR3 (p.W112R) were not detected in at least 892 reference chromosomes. ENG p.A60E showed significant association with familial IA in case-control studies (P=0.0080). No association with IA could be found for any of the known polymorphisms tested. CONCLUSIONS: Mutations in TGF-beta receptor genes are not a major cause of IA. However, we identified rare variants in ENG and TGFBR3 that may be important for IA pathogenesis in a subset of families.

Genomic screening identifies novel linkages and provides further evidence for a role of MYH9 in nonsyndromic cleft lip and palate.

Nonsyndromic cleft lip with or without cleft palate (NSCLP) is a common birth anomaly that requires prolonged multidisciplinary rehabilitation. Although variation in several genes has been identified as contributing to NSCLP, most of the genetic susceptibility loci have yet to be defined. To identify additional contributory genes, a high-throughput genomic scan was performed using the Illumina Linkage IVb Panel platform. We genotyped 6008 SNPs in nine non-Hispanic white NSCLP multiplex families and a single large African-American NSCLP multiplex family. Fourteen chromosomal regions were identified with LOD>1.5, including six regions not previously reported. Analysis of the data from the African-American and non-Hispanic white families revealed two likely chromosomal regions: 8q21.3-24.12 and 22q12.2-12.3 with LOD scores of 2.98 and 2.66, respectively. On the basis of biological function, syndecan 2 (SDC2) and growth differentiation factor 6 (GDF6) in 8q21.3-24.12 and myosin heavy-chain 9, non-muscle (MYH9) in 22q12.2-12.3 were selected as candidate genes. Association analyses from these genes yielded marginally significant P-values for SNPs in SDC2 and GDF6 (0.01

Tracking a tuberculosis outbreak over 21 years: strain-specific single nucleotide polymorphism-typing combined with targeted whole genome sequencing.

BACKGROUND  Whole genome sequencing (WGS) is increasingly used in molecular-epidemiological investigations of bacterial pathogens, despite cost- and time-intensive analyses. We combined strain-specific single nucleotide polymorphism (SNP)-typing and targeted WGS to investigate a tuberculosis cluster spanning 21 years in Bern, Switzerland. METHODS  Based on genome sequences of three historical outbreak Mycobacterium tuberculosis isolates, we developed a strain-specific SNP-typing assay to identify further cases. We screened 1,642 patient isolates, and performed WGS on all identified cluster isolates. We extracted SNPs to construct genomic networks. Clinical and social data were retrospectively collected. RESULTS  We identified 68 patients associated with the outbreak strain. Most were diagnosed in 1991-1995, but cases were observed until 2011. Two thirds belonged to the homeless and substance abuser milieu. Targeted WGS revealed 133 variable SNP positions among outbreak isolates. Genomic network analyses suggested a single origin of the outbreak, with subsequent division into three sub-clusters. Isolates from patients with confirmed epidemiological links differed by 0-11 SNPs. CONCLUSIONS  Strain-specific SNP-genotyping allowed rapid and inexpensive identification of M. tuberculosis outbreak isolates in a population-based strain collection. Subsequent targeted WGS provided detailed insights into transmission dynamics. This combined approach could be applied to track bacterial pathogens in real-time and at high resolution.

8p23 beta-defensin copy number determination by single-locus pseudogene-based paralog ratio tests risk bias due to low-frequency sequence variations.

BACKGROUND The copy number variation (CNV) in beta-defensin genes (DEFB) on human chromosome 8p23 has been proposed to contribute to the phenotypic differences in inflammatory diseases. However, determination of exact DEFB CN is a major challenge in association studies. Quantitative real-time PCR (qPCR), paralog ratio tests (PRT) and multiplex ligation-dependent probe amplification (MLPA) have been extensively used to determine DEFB CN in different laboratories, but inter-method inconsistencies were observed frequently. In this study we asked which one is superior among the three methods for DEFB CN determination. RESULTS We developed a clustering approach for MLPA and PRT to statistically correlate data from a single experiment. Then we compared qPCR, a newly designed PRT and MLPA for DEFB CN determination in 285 DNA samples. We found MLPA had the best convergence and clustering results of the raw data and the highest call rate. In addition, the concordance rates between MLPA or PRT and qPCR (32.12% and 37.99%, respectively) were unacceptably low with underestimated CN by qPCR. Concordance rate between MLPA and PRT (90.52%) was high but PRT systematically underestimated CN by one in a subset of samples. In these samples a sequence variant which caused complete PCR dropout of the respective DEFB cluster copies was found in one primer binding site of one of the targeted paralogous pseudogenes. CONCLUSION MLPA is superior to PRT and even more to qPCR for DEFB CN determination. Although the applied PRT provides in most cases reliable results, such a test is particularly sensitive to low-frequency sequence variations preferably accumulating in loci like pseudogenes which are most likely not under selective pressure. In the light of the superior performance of multiplex assays, the drawbacks of such single PRTs could be overcome by combining more test markers.

Candidate gene identification following linkage: Search for hypertension susceptibility genes

Following up genetic linkage studies to identify the underlying susceptibility gene(s) for complex disease traits is an arduous yet biologically and clinically important task. Complex traits, such as hypertension, are considered polygenic with many genes influencing risk, each with small effects. Chromosome 2 has been consistently identified as a genomic region with genetic linkage evidence suggesting that one or more loci contribute to blood pressure levels and hypertension status. Using combined positional candidate gene methods, the Family Blood Pressure Program has concentrated efforts in investigating this region of chromosome 2 in an effort to identify underlying candidate hypertension susceptibility gene(s). Initial informatics efforts identified the boundaries of the region and the known genes within it. A total of 82 polymorphic sites in eight positional candidate genes were genotyped in a large hypothesis-generating sample consisting of 1640 African Americans, 1339 whites, and 1616 Mexican Americans. To adjust for multiple comparisons, resampling-based false discovery adjustment was applied, extending traditional resampling methods to sibship samples. Following this adjustment for multiple comparisons, SLC4A5, a sodium bicarbonate transporter, was identified as a primary candidate gene for hypertension. Polymorphisms in SLC4A5 were subsequently genotyped and analyzed for validation in two populations of African Americans (N = 461; N = 778) and two of whites (N = 550; N = 967). Again, SNPs within SLC4A5 were significantly associated with blood pressure levels and hypertension status. While not identifying a single causal DNA sequence variation that is significantly associated with blood pressure levels and hypertension status across all samples, the results further implicate SLC4A5 as a candidate hypertension susceptibility gene, validating previous evidence for one or more genes on chromosome 2 that influence hypertension related phenotypes in the population-at-large. The methodology and results reported provide a case study of one approach for following up the results of genetic linkage analyses to identify genes influencing complex traits. ^

Characterization of glucocorticoid receptor (NR3C1) gene polymorphisms and a family-based association study with hypertension

Hypertension is usually defined as having values of systolic blood pressure ≥140 mmHg, diastolic blood pressure ≥90 mmHg. Hypertension is one of the main adverse effects of glucocorticoid on the cardiovascular system. Glucocorticoids are essential hormones, secreted from adrenal glands in circadian fashion. Glucocorticoid's effect on blood pressure is conveyed by the glucocorticoid receptor (NR3C1), an omnipresent nuclear transcription factor. Although polymorphisms in this gene have long been implicated to be a causal factor for cardiovascular diseases such as hypertension, no study has yet thoroughly interrogated the gene's polymorphisms for their effect on blood pressure levels. Therefore, I have first resequenced ∼30 kb of the gene, encompassing all exons, promoter regions, 5'/3' UTRs as well as at least 1.5 kb of the gene's flanking regions from 114 chromosome 5 monosomic cell lines, comprised of three major American ethnic groups—European American, African American and Mexican American. I observed 115 polymorphisms and 14 common molecularly phased haplotypes. A subset of markers was chosen for genotyping study populations of GENOA (Genetic Epidemiology Network of Atherosclerosis; 1022 non-Hispanic whites, 1228 African Americans and 954 Mexican Americans). Since these study populations include sibships, the family-based association test was performed on 4 blood pressure-related quantitative variables—pulse, systolic blood pressure, diastolic blood pressure and mean arterial pressure. Using these analyses, multiple correlated SNPs are significantly protective against high systolic blood pressure in non-Hispanic whites, which includes rsb198, a SNP formerly associated with beneficial body compositions. Haplotype association analysis also supports this finding and all p-values remained significant after permutation tests. I therefore conclude that multiple correlated SNPs on the gene may confer protection against high blood pressure in non-Hispanic whites. ^

Comparison of multiple comparison methods for identifying differential gene expression in simulated and real papillary thyroid cancer microarray data

The difficulty of detecting differential gene expression in microarray data has existed for many years. Several correction procedures try to avoid the family-wise error rate in multiple comparison process, including the Bonferroni and Sidak single-step p-value adjustments, Holm's step-down correction method, and Benjamini and Hochberg's false discovery rate (FDR) correction procedure. Each multiple comparison technique has its advantages and weaknesses. We studied each multiple comparison method through numerical studies (simulations) and applied the methods to the real exploratory DNA microarray data, which detect of molecular signatures in papillary thyroid cancer (PTC) patients. According to our results of simulation studies, Benjamini and Hochberg step-up FDR controlling procedure is the best process among these multiple comparison methods and we discovered 1277 potential biomarkers among 54675 probe sets after applying the Benjamini and Hochberg's method to PTC microarray data.^

Nonlinear dimension reduction in pathway-based genome-wide association analysis

Pathway based genome wide association study evolves from pathway analysis for microarray gene expression and is under rapid development as a complementary for single-SNP based genome wide association study. However, it faces new challenges, such as the summarization of SNP statistics to pathway statistics. The current study applies the ridge regularized Kernel Sliced Inverse Regression (KSIR) to achieve dimension reduction and compared this method to the other two widely used methods, the minimal-p-value (minP) approach of assigning the best test statistics of all SNPs in each pathway as the statistics of the pathway and the principal component analysis (PCA) method of utilizing PCA to calculate the principal components of each pathway. Comparison of the three methods using simulated datasets consisting of 500 cases, 500 controls and100 SNPs demonstrated that KSIR method outperformed the other two methods in terms of causal pathway ranking and the statistical power. PCA method showed similar performance as the minP method. KSIR method also showed a better performance over the other two methods in analyzing a real dataset, the WTCCC Ulcerative Colitis dataset consisting of 1762 cases, 3773 controls as the discovery cohort and 591 cases, 1639 controls as the replication cohort. Several immune and non-immune pathways relevant to ulcerative colitis were identified by these methods. Results from the current study provided a reference for further methodology development and identified novel pathways that may be of importance to the development of ulcerative colitis.^

Is Microarray Analysis Really Useful and Sufficient to Diagnose Nut Allergy in the Mediterranean Area?

Background: Component-based diagnosis on multiplex platforms is widely used in food allergy but its clinical performance has not been evaluated in nut allergy. Objective: To assess the diagnostic performance of a commercial protein microarray in the determination of specific IgE (sIgE) in peanut, hazelnut, and walnut allergy. Methods: sIgE was measured in 36 peanut-allergic, 36 hazelnut-allergic, and 44 walnut-allergic patients by ISAC 112, and subsequently, sIgE against available components was determined by ImmunoCAP in patients with negative ISAC results. ImmunoCAP was also used to measure sIgE to Ara h 9, Cor a 8, and Jug r 3 in a subgroup of lipid transfer protein (LTP)-sensitized nut-allergic patients (positive skin prick test to LTP-enriched extract). sIgE levels by ImmunoCAP were compared with ISAC ranges. Results: Most peanut-, hazelnut-, and walnut-allergic patients were sensitized to the corresponding nut LTP (Ara h 9, 66.7%; Cor a 8, 80.5%; Jug r 3, 84% respectively). However, ISAC did not detect sIgE in 33.3% of peanut-allergic patients, 13.9% of hazelnut-allergic patients, or 13.6% of walnut-allergic patients. sIgE determination by ImmunoCAP detected sensitization to Ara h 9, Cor a 8, and Jug r 3 in, respectively, 61.5% of peanut-allergic patients, 60% of hazelnut-allergic patients, and 88.3% of walnut-allergic patients with negative ISAC results. In the subgroup of peach LTP?sensitized patients, Ara h 9 sIgE was detected in more cases by ImmunoCAP than by ISAC (94.4% vs 72.2%, P<.05). Similar rates of Cor a 8 and Jug r 3 sensitization were detected by both techniques. Conclusions: The diagnostic performance of ISAC was adequate for hazelnut and walnut allergy but not for peanut allergy. sIgE sensitivity against Ara h 9 in ISAC needs to be improved.

All males are not created equal: Fertility differences depend on gamete recognition polymorphisms in sea urchins

Behaviors, morphologies, and genetic loci directly involved in reproduction have been increasingly shown to be polymorphic within populations. Explaining how such variants are maintained by selection is crucial to understanding the genetic basis of fertility differences, but direct tests of how alleles at reproductive loci affect fertility are rare. In the sea urchin genus Echinometra, the protein bindin mediates sperm attachment to eggs, evolves quickly, and is polymorphic within species. Eggs exposed to experimental sperm mixtures show strong discrimination on the basis of the males’ bindin genotype. Different females produce eggs that nonrandomly select sperm from different males, showing that variable egg–sperm interactions determine fertility. Eggs select sperm with a bindin genotype similar to their own, suggesting strong linkage between female choice and male trait loci. These experiments demonstrate that alleles at a single locus can have a strong effect on fertilization and that reproductive loci may retain functional polymorphisms through epistatic interactions between male and female traits. They also suggest that positive selection at gamete recognition loci like bindin involves strong selection within species on mate choice interactions.