Robert Kraaij

Best Practices and Joint Calling of the HumanExome BeadChip: The CHARGE Consortium

Genotyping arrays are a cost effective approach when typing previously-identified genetic polymorphisms in large numbers of samples. One limitation of genotyping arrays with rare variants (e.g., minor allele frequency [MAF] <0.01) is the difficulty that automated clustering algorithms have to accurately detect and assign genotype calls. Combining intensity data from large numbers of samples may increase the ability to accurately call the genotypes of rare variants. Approximately 62,000 ethnically diverse samples from eleven Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Consortium cohorts were genotyped with the Illumina HumanExome BeadChip across seven genotyping centers. The raw data files for the samples were assembled into a single project for joint calling. To assess the quality of the joint calling, concordance of genotypes in a subset of individuals having both exome chip and exome sequence data was analyzed. After exclusion of low performing SNPs on the exome chip and non-overlap of SNPs derived from sequence data, genotypes of 185,119 variants (11,356 were monomorphic) were compared in 530 individuals that had whole exome sequence data. A total of 98,113,070 pairs of genotypes were tested and 99.77% were concordant, 0.14% had missing data, and 0.09% were discordant. We report that joint calling allows the ability to accurately genotype rare variation using array technology when large sample sizes are available and best practices are followed. The cluster file from this experiment is available at www.chargeconsortium.com/main/exomechip.

An inactivating mutation of the luteinizing hormone receptor causes amenorrhea in a 46,XX female.

Hypergonadotropic hypogonadism is characterized by decreased gonadal function due to the inability of the gonads to respond to pituitary gonadotropins. Hypergonadotropic hypogonadism in females has many causes, among which are ovarian dysgenesis and abnormalities of the ovarian receptors for the pituitary gonadotropins. We evaluated a woman who presented with amenorrhea due to hypergonadotropic hypogonadism, but who had structurally normal ovaries. She is a sister of two previously identified 46,XY male pseudohermaphrodites with Leydig cell hypoplasia. Injection of hCG did not cause any change in plasma levels of estradiol or progesterone, suggesting complete ovarian resistance to LH. Analysis of the DNA sequence of the LH receptor gene revealed that the patient is homozygous for the same single base change as her two brothers. This mutation causes substitution of an alanine residue by a proline at position 593. In vitro analysis of the mutant LH receptor in cultured human embryonic kidney 293 cells documented that the receptor is unable to stimulate adenylyl cyclase in response to hCG. Plasma levels of estradiol and progesterone were low, whereas LH and FSH levels were increased. On histological analysis of the ovary, follicles were seen at all developmental stages. Nonetheless, primary amenorrhea had been present for 5 yr, and repeated measurements of plasma estradiol and progesterone indicate that ovulation does not occur. These results document the existence of inherited LH resistance as a cause of primary amenorrhea in women. The combined clinical and molecular observations are consistent with previous experimental data suggesting that in humans, LH is necessary for ovulation but follicular maturation can occur in the presence of FSH alone.

Quantifying prion disease penetrance using large population control cohorts

Large genomic reference data sets reveal a spectrum of pathogenicity in the prion protein gene and provide genetic validation for a therapeutic strategy in prion disease. Share trumps rare No longer just buzz words, “patient empowerment” and “data sharing” are enabling breakthrough research on rare genetic diseases. Although more than 100,000 genetic variants are believed to drive disease in humans, little is known about penetrance—the probability that a mutation will actually cause disease in the carrier. This conundrum persists because small sample sizes breed imperfect alliance estimates between mutations and disease risk. Now, a patient-turned-scientist joined with a large bioinformatics team to analyze vast amounts of shared data—from the Exome Aggregation Consortium and the 23andMe database—to provide insights into genetic-variant penetrance and possible treatment approaches for a rare, fatal genetic prion disease. More than 100,000 genetic variants are reported to cause Mendelian disease in humans, but the penetrance—the probability that a carrier of the purported disease-causing genotype will indeed develop the disease—is generally unknown. We assess the impact of variants in the prion protein gene (PRNP) on the risk of prion disease by analyzing 16,025 prion disease cases, 60,706 population control exomes, and 531,575 individuals genotyped by 23andMe Inc. We show that missense variants in PRNP previously reported to be pathogenic are at least 30 times more common in the population than expected on the basis of genetic prion disease prevalence. Although some of this excess can be attributed to benign variants falsely assigned as pathogenic, other variants have genuine effects on disease susceptibility but confer lifetime risks ranging from <0.1 to ~100%. We also show that truncating variants in PRNP have position-dependent effects, with true loss-of-function alleles found in healthy older individuals, a finding that supports the safety of therapeutic suppression of prion protein expression.

Truncated ETV1, Fused to Novel Tissue-Specific Genes, and Full-Length ETV1 in Prostate Cancer

In this study, we describe the properties of novel ETV1 fusion genes, encoding N-truncated ETV1 (dETV1), and of full-length ETV1, overexpressed in clinical prostate cancer. We detected overexpression of novel ETV1 fusion genes or of full-length ETV1 in 10% of prostate cancers. Novel ETV1 fusion partners included FOXP1, an EST (EST14), and an endogenous retroviral repeat sequence (HERVK17). Like TMPRSS2, EST14 and HERVK17 were prostate-specific and androgen-regulated expressed. This unique expression pattern of most ETV1 fusion partners seems an important determinant in prostate cancer development. In transient reporter assays, full-length ETV1 was a strong transactivator, whereas dETV1 was not. However, several of the biological properties of dETV1 and full-length ETV1 were identical. On stable overexpression, both induced migration and invasion of immortalized nontumorigenic PNT2C2 prostate epithelial cells. In contrast to dETV1, full-length ETV1 also induced anchorage-independent growth of these cells. PNT2C2 cells stably transfected with dETV1 or full-length ETV1 expression constructs showed small differences in induced expression of target genes. Many genes involved in tumor invasion/metastasis, including uPA/uPAR and MMPs, were up-regulated in both cell types. Integrin beta3 (ITGB3) was clearly up-regulated by full-length ETV1 but much less by dETV1. Based on the present data and on previous findings, a novel concept of the role of dETV1 and of full-length ETV1 overexpression in prostate cancer is proposed.

Disconnecting the Yin and Yang Relation of Epidermal Growth Factor Receptor (EGFR)-Mediated Delivery: A Fully Synthetic, EGFR-Targeted Gene Transfer System Avoiding Receptor Activation

The epidermal growth factor receptor (EGFR) is upregulated within a high percentage of solid tumors and hence is an attractive target for tumor-targeted therapies including gene therapy. The natural EGFR ligand epidermal growth factor (EGF) has been used for this purpose, despite the risk of mitogenic effects due to EGFR activation. We have developed a fully synthetic, EGFR-targeted gene delivery system based on PEGylated linear polyethylenimine (LPEI), allowing evaluation of different EGFR-binding peptides in terms of transfection efficiency and EGFR activation. Peptide sequences directly derived from the human EGF molecule enhanced transfection efficiency with concomitant EGFR activation. Only the EGFR-binding peptide GE11, which has been identified by phage display technique, showed specific enhancement of transfection on EGFR-overexpressing tumor cells including glioblastoma and hepatoma, but without EGFR activation. EGFR targeting led to high levels of cell association of fluorescently labeled polyplexes after only 30 min of incubation. EGF pretreatment of cells induced enhanced cellular internalization of all polyplex types tested, pointing at generally enhanced macropinocytosis. EGF polyplexes diminished cell surface expression of EGFR for up to 4 hr, whereas GE11 polyplexes did not. In a clinically relevant orthotopic prostate cancer model, intratumorally injected GE11 polyplexes were superior in inducing transgene expression when compared with untargeted polyplexes.

Validation of transrectal ultrasonographic volumetry for orthotopic prostate tumours in mice

Orthotopic human prostate tumour models in athymic nude mice are regarded as being most suitable for fundamental and pre-clinical research on prostate cancer. The anatomic localization of the tumour in the pelvis, however, provides little possibility for monitoring tumour growth or regression. To assess time-related changes in orthotopic tumour volume, we applied transrectal ultrasonography (TRUS) to the murine prostate. This technique has the advantages of allowing accurate monitoring of tumours during therapeutic manipulations and a reduction of animal use due to a reduction of sacrificing endpoints. To validate the TRUS method, the mouse prostate reconstitution model, RM-9, and the prostate-specific antigen (PSA) producing human prostate cancer xenograft PC-346 were used. Volumetric calliper measurements were performed with a 30 MHz ultrasound probe designed for intra-arterial use in humans. Tumour weight, determined at various time-points, was found to be closely related to actual tumour weight (R = 0.99) and, in the PC-346 model, to the level of PSA in the plasma. Furthermore, the interobserver variation for TRUS was low for tumours above 50 mg. Thus, TRUS for murine prostate tumours proves to be an accurate, reproducible and sensitive method.

Regulation of gonadotropin receptor gene expression

The receptors for the gonadotropins differ from the other G protein-coupled receptors by having a large extracellular hormone-binding domain, encoded by nine or ten exons. Alternative splicing of the large pre-mRNA of approximately 100 kb can result in mRNA species that encode truncated receptor proteins. In this review we discuss the regulation of gonadotropin receptor mRNA expression and the possible roles of alternative splicing in gonadotropin receptor function.

Novel Genetic Variants for Cartilage Thickness and Hip Osteoarthritis.

Osteoarthritis is one of the most frequent and disabling diseases of the elderly. Only few genetic variants have been identified for osteoarthritis, which is partly due to large phenotype heterogeneity. To reduce heterogeneity, we here examined cartilage thickness, one of the structural components of joint health. We conducted a genome-wide association study of minimal joint space width (mJSW), a proxy for cartilage thickness, in a discovery set of 13,013 participants from five different cohorts and replication in 8,227 individuals from seven independent cohorts. We identified five genome-wide significant (GWS, P≤5·0×10−8) SNPs annotated to four distinct loci. In addition, we found two additional loci that were significantly replicated, but results of combined meta-analysis fell just below the genome wide significance threshold. The four novel associated genetic loci were located in/near TGFA (rs2862851), PIK3R1 (rs10471753), SLBP/FGFR3 (rs2236995), and TREH/DDX6 (rs496547), while the other two (DOT1L and SUPT3H/RUNX2) were previously identified. A systematic prioritization for underlying causal genes was performed using diverse lines of evidence. Exome sequencing data (n = 2,050 individuals) indicated that there were no rare exonic variants that could explain the identified associations. In addition, TGFA, FGFR3 and PIK3R1 were differentially expressed in OA cartilage lesions versus non-lesioned cartilage in the same individuals. In conclusion, we identified four novel loci (TGFA, PIK3R1, FGFR3 and TREH) and confirmed two loci known to be associated with cartilage thickness.The identified associations were not caused by rare exonic variants. This is the first report linking TGFA to human OA, which may serve as a new target for future therapies.

A transductionally retargeted adenoviral vector for virotherapy of Her2/neu-expressing prostate cancer.

The efficacy of adenovirus (Ad)-based gene therapy of solid tumors, such as prostate cancer, is limited. One of the many problems is that the virus infects many different cell types in the body, resulting in high toxicity, whereas the target cancer cells are often less prone to wild-type Ad infection. Our aim was to develop genetically de- and retargeted Ad vectors to reduce off-target effects and increase target infection for prostate cancer. We have previously reported an Ad5 vector specific for the cancer-associated receptor Her2/neu, created by inserting Her2/neu-reactive Affibody(®) molecules (ZH) into the HI loop of a coxsackievirus and adenovirus receptor binding-ablated fiber (Ad[ZH/1]). In addition to virus retargeting to Her2/neu, this virus was further modified from wild-type Ad by changing the RGD motif in the penton base to EGD and by substitution of the KKTK motif in the third shaft repeat to RKSK, resulting in the vector Ad[ZH/3]. The ZH-containing vectors could be produced to high titers and were specific for their target, resulting in efficient infection and killing of Her2/neu-positive androgen-dependent PC346C prostate cancer cells in vitro. Here we show that the oncolytic Ad[ZH/3] vector significantly prolonged survival time and reduced serum prostate-specific antigen levels in an orthotopic prostate tumor model in nude mice to the same extent as wild-type Ad5. Our results show that Her2/neu targeting using Ad-based vectors for prostate cancer is feasible and may serve as a basis for the development of gene therapy of human prostate cancer as well as other Her2/neu-expressing cancers.

Excessive burden of lysosomal storage disorder gene variants in Parkinson's disease

Mutations in the glucocerebrosidase gene (GBA), which cause Gaucher disease, are also potent risk factors for Parkinsons disease. We examined whether a genetic burden of variants in other lysosomal storage disorder genes is more broadly associated with Parkinsons disease susceptibility. The sequence kernel association test was used to interrogate variant burden among 54 lysosomal storage disorder genes, leveraging whole exome sequencing data from 1156 Parkinsons disease cases and 1679 control subjects. We discovered a significant burden of rare, likely damaging lysosomal storage disorder gene variants in association with Parkinsons disease risk. The association signal was robust to the exclusion of GBA, and consistent results were obtained in two independent replication cohorts, including 436 cases and 169 controls with whole exome sequencing and an additional 6713 cases and 5964 controls with exome-wide genotyping. In secondary analyses designed to highlight the specific genes driving the aggregate signal, we confirmed associations at the GBA and SMPD1 loci and newly implicate CTSD, SLC17A5, and ASAH1 as candidate Parkinsons disease susceptibility genes. In our discovery cohort, the majority of Parkinsons disease cases (56%) have at least one putative damaging variant in a lysosomal storage disorder gene, and 21% carry multiple alleles. Our results highlight several promising new susceptibility loci and reinforce the importance of lysosomal mechanisms in Parkinsons disease pathogenesis. We suggest that multiple genetic hits may act in combination to degrade lysosomal function, enhancing Parkinsons disease susceptibility.