Two missense mutations in SALL4 in a patient with microphthalmia, coloboma, and optic nerve hypoplasia.

To investigate the genetic etiology of anophthalmia and microphthalmia, we used exome sequencing in a Caucasian female with unilateral microphthalmia and coloboma, bilateral optic nerve hypoplasia, ventricular and atrial septal defects, and growth delays. We found two sequence variants in SALL4 - c.[575C>A], predicting p.(Ala192Glu), that was paternally inherited, and c.[2053G>C], predicting p.(Asp685His), that was maternally inherited. Haploinsufficiency for SALL4 due to nonsense or frameshift mutations has been associated with acro-renal ocular syndrome that is characterized by eye defects including Duane anomaly and coloboma, in addition to radial ray malformations and renal abnormalities. Our report is the first description of structural eye defects associated with two missense variants in SALL4 inherited in trans; the absence of reported findings in both parents suggests that both sequence variants are hypomorphic mutations and that both are needed for the ocular phenotype. SALL4 is expressed in the developing lens and regulates BMP4, leading us to speculate that altered BMP4 expression was responsible for the eye defects, but we could not demonstrate altered BMP4 expression in vitro after using small interfering RNAs (siRNAs) to reduce SALL4 expression. We conclude that SALL4 hypomorphic variants may influence eye development.

Genomewide Association Study of Tacrolimus Concentrations in African American Kidney Transplant Recipients Identifies Multiple CYP3A5 Alleles.

We previously reported that tacrolimus (TAC) trough blood concentrations for African American (AA) kidney allograft recipients were lower than those observed in white patients. Subtherapeutic TAC troughs may be associated with acute rejection (AR) and AR-associated allograft failure. This variation in TAC troughs is due, in part, to differences in the frequency of the cytochrome P450 CYP3A5*3 allele (rs776746, expresses nonfunctional enzyme) between white and AA recipients; however, even after accounting for this variant, variability in AA-associated troughs is significant. We conducted a genomewide association study of TAC troughs in AA kidney allograft recipients to search for additional genetic variation. We identified two additional CYP3A5 variants in AA recipients independently associated with TAC troughs: CYP3A5*6 (rs10264272) and CYP3A5*7 (rs41303343). All three variants and clinical factors account for 53.9% of the observed variance in troughs, with 19.8% of the variance coming from demographic and clinical factors including recipient age, glomerular filtration rate, anticytomegalovirus drug use, simultaneous pancreas-kidney transplant and antibody induction. There was no evidence of common genetic variants in AA recipients significantly influencing TAC troughs aside from the CYP3A gene. These results reveal that additional and possibly rare functional variants exist that account for the additional variation.

Exome sequencing in 32 patients with anophthalmia/microphthalmia and developmental eye defects.

Anophthalmia/microphthalmia (A/M) is a genetically heterogeneous birth defect for which the etiology is unknown in more than 50% of patients. We used exome sequencing with the ACE Exome(TM) (Personalis, Inc; 18 cases) and UCSF Genomics Core (21 cases) to sequence 28 patients with A/M and four patients with varied developmental eye defects. In the 28 patients with A/M, we identified de novo mutations in three patients (OTX2, p.(Gln91His), RARB, p.Arg387Cys and GDF6, p.Ala249Glu) and inherited mutations in STRA6 in two patients. In patients with developmental eye defects, a female with cataracts and cardiomyopathy had a de novo COL4A1 mutation, p.(Gly773Arg), expanding the phenotype associated with COL4A1 to include cardiomyopathy. A male with a chorioretinal defect, microcephaly, seizures and sensorineural deafness had two PNPT1 mutations, p.(Ala507Ser) and c.401-1G>A, and we describe eye defects associated with this gene for the first time. Exome sequencing was efficient for identifying mutations in pathogenic genes for which there is no clinical testing available and for identifying cases that expand phenotypic spectra, such as the PNPT1 and COL4A1-associated disorders described here.

Gene-based association identifies SPATA13-AS1 as a pharmacogenomic predictor of inhaled short-acting beta-agonist response in multiple population groups.

Inhaled short-acting beta-agonist (SABA) medication is commonly used in asthma patients to rapidly reverse airway obstruction and improve acute symptoms. We performed a genome-wide association study of SABA medication response using gene-based association tests. A linear mixed model approach was first used for single-nucleotide polymorphism associations, and the results were later combined using GATES to generate gene-based associations. Our results identified SPATA13-AS1 as being significantly associated with SABA bronchodilator response in 328 healthy African Americans. In replication, this gene was associated with SABA response among the two separate groups of African Americans with asthma (n=1073, P=0.011 and n=1968, P=0.014), 149 healthy African Americans (P=0.003) and 556 European Americans with asthma (P=0.041). SPATA13-AS1 was also associated with longitudinal SABA medication usage in the two separate groups of African Americans with asthma (n=658, P=0.047 and n=1968, P=0.025). Future studies are needed to delineate the precise mechanism by which SPATA13-AS1 may influence SABA response.

A rare WNT1 missense variant overrepresented in ASD leads to increased Wnt signal pathway activation.

Wnt signaling, which encompasses multiple biochemical pathways that regulate neural development downstream of extracellular Wnt glycoprotein ligands, has been suggested to contribute to major psychiatric disorders including autism spectrum disorders (ASD). We used next-generation sequencing and Sequenom genotyping technologies to resequence 10 Wnt signaling pathway genes in 198 ASD patients and 240 matched controls. Results for single-nucleotide polymorphisms (SNPs) of interest were confirmed in a second set of 91 ASD and 144 control samples. We found a significantly increased burden of extremely rare missense variants predicted to be deleterious by PolyPhen-2, distributed across seven genes in the ASD sample (3.5% in ASD vs 0.8% in controls; Fisher's exact test, odds ratio (OR)=4.37, P=0.04). We also found a missense variant in WNT1 (S88R) that was overrepresented in the ASD sample (8 A/T in 267 ASD (minor allele frequency (MAF)=1.69%) vs 1 A/T in 377 controls (MAF=0.13%), OR=13.0, Fisher's exact test, P=0.0048; OR=8.2 and P=0.053 after correction for population stratification). Functional analysis revealed that WNT1-S88R is more active than wild-type WNT1 in assays for the Wnt/ß-catenin signaling pathway. Our findings of a higher burden in ASD of rare missense variants distributed across 7 of 10 Wnt signaling pathway genes tested, and of a functional variant at the WNT1 locus associated with ASD, support that dysfunction of this pathway contributes to ASD susceptibility. Given recent findings of common molecular mechanisms in ASD, schizophrenia and affective disorders, these loci merit scrutiny in other psychiatric conditions as well.

Beta-adrenergic receptor polymorphisms and cardiac graft function in potential organ donors.

Prior studies have demonstrated associations between beta-adrenergic receptor (ßAR) polymorphisms and left ventricular dysfunction-an important cause of allograft nonutilization for transplantation. We hypothesized that ßAR polymorphisms predispose donor hearts to LV dysfunction after brain death. A total of 1043 organ donors managed from 2001-2006 were initially studied. The following ßAR single nucleotide polymorphisms were genotyped: ß1AR 1165C/G (Arg389Gly), ß1AR 145A/G (Ser49Gly), ß2AR 46G/A (Gly16Arg) and ß2AR 79C/G (Gln27Glu). In multivariable regression analyses, the ß2AR46 SNP was significantly associated with LV systolic dysfunction, with each minor allele additively decreasing the odds for LV ejection fraction <50%. The ß1AR1165 and ß2AR46 SNPs were associated with higher dopamine requirement during the donor management period: donors with the GG and AA genotypes had ORs of 2.64 (95% CI 1.52-4.57) and 2.70 (1.07-2.74) respectively for requiring >10 µg/kg/min of dopamine compared to those with the CC and GG genotypes. However, no significant associations were found between ßAR SNPs and cardiac dysfunction in 364 donors managed from 2007-2008, perhaps due to changes in donor management, lack of power in this validation cohort, or the absence of a true association. ßAR polymorphisms may be associated with cardiac dysfunction after brain death, but these relationships require further study in independent donor cohorts.

A screening study of drug-drug interactions in cerivastatin users: an adverse effect of clopidogrel.

An analysis of a case-control study of rhabdomyolysis was conducted to screen for previously unrecognized cytochrome P450 enzyme (CYP) 2C8 inhibitors that may cause other clinically important drug-drug interactions. Medication use in cases of rhabdomyolysis using cerivastatin (n = 72) was compared with that in controls using atorvastatin (n = 287) for the period 1998-2001. The use of clopidogrel was strongly associated with rhabdomyolysis (odds ratio (OR) 29.6; 95% confidence interval (CI), 6.1-143). In a replication effort that used the US Food and Drug Administration (FDA) Adverse Event Reporting System (AERS), it was found that clopidogrel was used more commonly in patients with rhabdomyolysis receiving cerivastatin (17%) than in those receiving atorvastatin (0%, OR infinity; 95% CI = 5.2-infinity). Several medications were tested in vitro for their potential to cause drug-drug interactions. Clopidogrel, rosiglitazone, and montelukast were the most potent inhibitors of cerivastatin metabolism. Clopidogrel and its metabolites also inhibited cerivastatin metabolism in human hepatocytes. These epidemiological and in vitro findings suggest that clopidogrel may cause clinically important, dose-dependent drug-drug interactions with other medications metabolized by CYP2C8.

A common 5'-UTR variant in MATE2-K is associated with poor response to metformin.

Multidrug and toxin extrusion 2 (MATE2-K (SLC47A2)), a polyspecific organic cation exporter, facilitates the renal elimination of the antidiabetes drug metformin. In this study, we characterized genetic variants of MATE2-K, determined their association with metformin response, and elucidated their impact by means of a comparative protein structure model. Four nonsynonymous variants and four variants in the MATE2-K basal promoter region were identified from ethnically diverse populations. Two nonsynonymous variants-c.485C>T and c.1177G>A-were shown to be associated with significantly lower metformin uptake and reduction in protein expression levels. MATE2-K basal promoter haplotypes containing the most common variant, g.-130G>A (>26% allele frequency), were associated with a significant increase in luciferase activities and reduced binding to the transcriptional repressor myeloid zinc finger 1 (MZF-1). Patients with diabetes who were homozygous for g.-130A had a significantly poorer response to metformin treatment, assessed as relative change in glycated hemoglobin (HbA1c) (-0.027 (-0.076, 0.033)), as compared with carriers of the reference allele, g.-130G (-0.15 (-0.17, -0.13)) (P=0.002). Our study showed that MATE2-K plays a role in the antidiabetes response to metformin.

Functional characterization of liver enhancers that regulate drug-associated transporters.

Little is known about how genetic variations in enhancers influence drug response. In this study, we investigated whether nucleotide variations in enhancers that regulate drug transporters can alter their expression levels. Using comparative genomics and liver-specific transcription factor binding site (TFBS) analyses, we identified evolutionary conserved regions (ECRs) surrounding nine liver membrane transporters that interact with commonly used pharmaceuticals. The top 50 ECRs were screened for enhancer activity in vivo, of which five--located around ABCB11, SLC10A1, SLCO1B1, SLCO1A2, and SLC47A1--exhibited significant enhancer activity. Common variants identified in a large ethnically diverse cohort (n = 272) were assayed for differential enhancer activity, and three variants were found to have significant effects on reporter activity as compared with the reference allele. In addition, one variant was associated with reduced SLCO1A2 mRNA expression levels in human liver tissues, and another was associated with increased methotrexate (MTX) clearance in patients. This work provides a general model for the rapid characterization of liver enhancers and identifies associations between enhancer variants and drug response.

Sequencing of TNFAIP3 and association of variants with multiple autoimmune diseases.

The TNFAIP3 locus at 6q23, encoding A20, has been associated with multiple autoimmune diseases (AIDs). In this study, we sequence the coding portions of the gene to identify contributing causal polymorphisms that may explain some of the observed associations. A collection of 123 individuals from the Multiple Autoimmune Disease Genetics Consortium (MADGC) collection, each with multiple AIDs (mean=2.2 confirmed diagnoses), and 397 unrelated healthy controls were used for initial sequencing. A total of 32 polymorphisms were identified in the sequencing experiments, including 16 novel and 11 coding variants. Association testing in the entire MADGC collection (1,008 Caucasians with one or more AIDs and 770 unaffected family controls) revealed association of a novel intronic insertion-deletion polymorphism with rheumatoid arthritis (RA) (odds ratio (OR)=2.48, P=0.041). Genotyping of the most common coding polymorphism, rs2230926, in the MADGC collection and additional control individuals revealed a significant association with Sjögren's syndrome (OR=3.38, P=0.038), Crohn's disease (OR=2.25, P=0.041), psoriasis (OR=0.037, P=0.036) and RA (OR=1.9, P=0.025). Finally, haplotype and additional testing of polymorphisms revealed that cases were enriched for 5' and 3' untranslated region variants (one-sided P-value=0.04), but not specifically for common (>2% minor allele frequency), rare, exonic, intronic, non-synonymous or synonymous variants.

Admixture mapping of ankle-arm index: identification of a candidate locus associated with peripheral arterial disease.

BACKGROUND: Peripheral arterial disease (PAD) is associated with significant morbidity and mortality, and has a higher prevalence in African Americans than Caucasians. Ankle-arm index (AAI) is the ratio of systolic blood pressure in the leg to that in the arm, and, when low, is a marker of PAD. METHODS: The authors used an admixture mapping approach to search for genetic loci associated with low AAI. Using data from 1040 African American participants in the observational, population based Health, Aging, and Body Composition Study who were genotyped at 1322 single nucleotide polymorphisms (SNPs) that are informative for African versus European ancestry and span the entire genome, we estimated genetic ancestry in each chromosomal region and then tested the association between AAI and genetic ancestry at each locus. RESULTS: The authors found a region of chromosome 11 that reaches its peak between 80 and 82 Mb associated with low AAI (p<0.001 for rs12289502 and rs9665943, both within this region). 753 African American participants in the observational, population based Cardiovascular Health Study were genotyped at rs9665943 to test the reproducibility of this association, and this association was also statistically significant (odds ratio (OR) for homozygous African genotype 1.59, 95% confidence interval (CI) 1.12 to 2.27). Another candidate SNP (rs1042602) in the same genomic region was tested in both populations, and was also found to be significantly associated with low AAI in both populations (OR for homozygous African genotype 1.89, 95% CI 1.29 to 2.76). CONCLUSION: This study identifies a novel region of chromosome 11 representing an area with a potential candidate gene associated with PAD in African Americans.

Heat-inducible translationally controlled tumor protein of Trichinella pseudospiralis: cloning and regulation of gene expression.

To elucidate the mechanism of inducing translationally controlled tumor protein (TCTP) in stress adaptation of adenophorean nematodes, the complete coding sequence of TCTP of the infective-stage larvae of Trichinella pseudospiralis was characterized. Two cDNA clones with different 3' untranslated region were identified. Tp-TCTP contained an open reading frame of 534 bp encoding 177 residues. The gene with five introns was expressed as histidine-tagged fusion protein having a molecular mass of 17.5 kDa. Quantitative reverse transcriptase polymerase chain reaction (RT-PCR) analysis showed that TCTP RNA was not accumulated when the infective-stage larvae were heat-shocked for 1 h at 45 or 60 degrees C. Using enzyme-linked immunosorbent assay and antiserum against the fusion protein, the expression of TCTP was found to be up-regulated at the translational level. The data suggest that translational regulation of TCTP may play an important role in the early heat-stress adaptation of the trichinellid. Cluster analysis demonstrated that the TCTP sequence of T. pseudospiralis is closely related to that of T. spiralis, but is diverged from the secernentean species.

Linkage disequilibrium maps constructed with common SNPs are useful for first-pass disease association screens.

To develop an efficient strategy for mapping genetic factors associated with common diseases, we constructed linkage disequilibrium (LD) maps of human chromosomes 5, 7, 17, and X. These maps consist of common single nucleotide polymorphisms at an average intermarker distance of 100 kb. The genotype data from these markers in a panel of American samples of European descent were analyzed to produce blocks of markers in strong pair-wise LD. Power calculations were used to guide block definitions and predicted that high-level LD maps would be useful in initial genome scans for susceptibility alleles in case-control association studies of complex diseases. As anticipated, LD blocks on the X chromosome were larger and covered more of the chromosome than those found on the autosomes.

Methods for genotyping single nucleotide polymorphisms.

One of the fruits of the Human Genome Project is the discovery of millions of DNA sequence variants in the human genome. The majority of these variants are single nucleotide polymorphisms (SNPs). A dense set of SNP markers opens up the possibility of studying the genetic basis of complex diseases by population approaches. In all study designs, a large number of individuals must be genotyped with a large number of markers. In this review, the current status of SNP genotyping is discussed in terms of the mechanisms of allelic discrimination, the reaction formats, and the detection modalities. A number of genotyping methods currently in use are described to illustrate the approaches being taken. Although no single genotyping method is ideally suited for all applications, a number of good genotyping methods are available to meet the needs of many study designs. The challenges for SNP genotyping in the near future include increasing the speed of assay development, reducing the cost of the assays, and performing multiple assays in parallel. Judging from the accelerated pace of new method development, it is hopeful that an ideal SNP genotyping method will be developed soon.

The birth and death of human single-nucleotide polymorphisms: new experimental evidence and implications for human history and medicine.

Extensive, new databases of single-nucleotide polymorphisms (SNPs) provide a powerful resource for disease gene discovery, and they will be even more useful as more frequency data become available. Interesting observed genomic patterns include SNP deserts (regions of low SNP incidence) and lengthy regions of linkage disequilibrium containing only a few haplotypes. A variety of genetic studies will benefit from SNP resources.

Universal SNP genotyping assay with fluorescence polarization detection.

The degree of fluorescence polarization (FP) of a fluorescent molecule is a reflection of its molecular weight (Mr). FP is therefore a useful detection methodfor homogeneous assays in which the starting reagents and products differ significantly in Mr. We have previously shown that FP is a good detection method for the single-base extension and the 5'-nuclease assays. In this report, we describe a universal, optimized single-base extension assay for genotyping single nucleotide polymorphisms (SNPs). This assay, which we named the template-directed dye-terminator incorporation assay with fluorescence polarization detection (FP-TDI), uses four spectrally distinct dye terminators to achieve universal assay conditions. Even without optimization, approximately 70% of all SNP markers tested yielded robust assays. The addition of an E. coli ssDNA-binding protein just before the FP reading significantly increased FP values of the products and brought the success rate of FP-TDI assays up to 90%. Increasing the amount of dye terminators and reducing the number of thermal cycles in the single-base extension step of the assay increased the separation of the FP values benveen the products corresponding to different genotypes and improved the success rate of the assay to 100%. In this study the genomic DNA samples of 90 individuals were typed for a total of 38 FP-TDI assays (using both the sense and antisense TDI primers for 19 SNP markers). With the previously described modifications, the FP-TDI assay gave unambiguous genotyping data for all the samples tested in the 38 FP-TDI assays. When the genotypes determined by the FP-TDI and 5'-nuclease assays were compared, they were in 100% concordance for all experiments (a total of 3420 genotypes). The four-dye-terminator master mixture described here can be used for assaying any SNP marker and greatly simplifies the SNP genotyping assay design.

Genotyping single-nucleotide polymorphisms by the invader assay with dual-color fluorescence polarization detection.

BACKGROUND: The PCR-Invader assay is a robust, homogeneous assay that has been shown to be highly sensitive and specific in genotyping single-nucleotide polymorphism (SNP) markers. In this study, we introduce two changes to improve the assay: (a) we streamline the PCR-Invader method by assaying both alleles for each SNP in one reaction; and (b) we reduce the cost of the method by adopting fluorescence polarization (FP) as the detection method. METHODS: PCR product was incubated with Invader oligonucleotide and two primary probes at 93 degrees C for 5 min. Signal probes corresponding to the cleaved flaps of the primary probes [labeled with fluorescein and 6-carboxytetramethylrhodamine (TAMRA) dye] and Cleavase VIII enzyme (a flap endonuclease) were then added to the mixture. This reaction mixture was incubated at 63 degrees C for 5 min. FP measurements were made with a fluorescence plate reader. RESULTS: Eighty-eight individuals were genotyped across a panel of 10 SNPs, using PCR product as template, for a total of 880 genotypes. An average "no call" rate of 3.2% was observed after first round of experiments. PCR products were remade in those samples that failed to produce any genotype in the first round, and all gave clear-cut genotypes. When the genotypes determined by the PCR-Invader assay and template-directed dye-terminator incorporation assay with FP were compared, they were in 100% concordance for all SNP markers and experiments. CONCLUSIONS: The improvements introduced in this study make PCR-Invader assay simpler and more cost-effective, and therefore more suitable for high-throughput genotyping.