Host variation in type I interferon signaling genes (MX1), C-C chemokine receptor type 5 gene, and major histocompatibility complex class I alleles in treated HIV+ noncontrollers predict viral reservoir size.

OBJECTIVE: Prior genomewide association studies have identified variation in major histocompatibility complex (MHC) class I alleles and C-C chemokine receptor type 5 gene (CCR5Δ32) as genetic predictors of viral control, especially in 'elite' controllers, individuals who remain virally suppressed in the absence of therapy. DESIGN: Cross-sectional genomewide association study. METHODS: We analyzed custom whole exome sequencing and direct human leukocyte antigen (HLA) typing from 202 antiretroviral therapy (ART)-suppressed HIV+ noncontrollers in relation to four measures of the peripheral CD4+ T-cell reservoir: HIV intact DNA, total (t)DNA, unspliced (us)RNA, and RNA/DNA. Linear mixed models were adjusted for potential covariates including age, sex, nadir CD4+ T-cell count, pre-ART HIV RNA, timing of ART initiation, and duration of ART suppression. RESULTS: Previously reported 'protective' host genetic mutations related to viral setpoint (e.g. among elite controllers) were found to predict smaller HIV reservoir size. The HLA 'protective' B∗57:01 was associated with significantly lower HIV usRNA (q = 3.3 × 10-3), and among the largest subgroup, European ancestry individuals, the CCR5Δ32 deletion was associated with smaller HIV tDNA (P = 4.3 × 10-3) and usRNA (P = 8.7 × 10-3). In addition, genomewide analysis identified several single nucleotide polymorphisms in MX1 (an interferon stimulated gene) that were significantly associated with HIV tDNA (q = 0.02), and the direction of these associations paralleled MX1 gene eQTL expression. CONCLUSIONS: We observed a significant association between previously reported 'protective' MHC class I alleles and CCR5Δ32 with the HIV reservoir size in noncontrollers. We also found a novel association between MX1 and HIV total DNA (in addition to other interferon signaling relevant genes, PPP1CB, DDX3X). These findings warrant further investigation in future validation studies.

A Large-Scale Association Study Detects Novel Rare Variants, Risk Genes, Functional Elements, and Polygenic Architecture of Prostate Cancer Susceptibility.

To identify rare variants associated with prostate cancer susceptibility and better characterize the mechanisms and cumulative disease risk associated with common risk variants, we conducted an integrated study of prostate cancer genetic etiology in two cohorts using custom genotyping microarrays, large imputation reference panels, and functional annotation approaches. Specifically, 11,984 men (6,196 prostate cancer cases and 5,788 controls) of European ancestry from Northern California Kaiser Permanente were genotyped and meta-analyzed with 196,269 men of European ancestry (7,917 prostate cancer cases and 188,352 controls) from the UK Biobank. Three novel loci, including two rare variants (European ancestry minor allele frequency < 0.01, at 3p21.31 and 8p12), were significant genome wide in a meta-analysis. Gene-based rare variant tests implicated a known prostate cancer gene (HOXB13), as well as a novel candidate gene (ILDR1), which encodes a receptor highly expressed in prostate tissue and is related to the B7/CD28 family of T-cell immune checkpoint markers. Haplotypic patterns of long-range linkage disequilibrium were observed for rare genetic variants at HOXB13 and other loci, reflecting their evolutionary history. In addition, a polygenic risk score (PRS) of 188 prostate cancer variants was strongly associated with risk (90th vs. 40th-60th percentile OR = 2.62, P = 2.55 × 10-191). Many of the 188 variants exhibited functional signatures of gene expression regulation or transcription factor binding, including a 6-fold difference in log-probability of androgen receptor binding at the variant rs2680708 (17q22). Rare variant and PRS associations, with concomitant functional interpretation of risk mechanisms, can help clarify the full genetic architecture of prostate cancer and other complex traits. SIGNIFICANCE: This study maps the biological relationships between diverse risk factors for prostate cancer, integrating different functional datasets to interpret and model genome-wide data from over 200,000 men with and without prostate cancer.See related commentary by Lachance, p. 1637.

Author Correction: Multiplexed droplet single-cell RNA-sequencing using natural genetic variation.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Whole exome and whole genome sequencing with dried blood spot DNA without whole genome amplification.

Newborn screening (NBS) for rare conditions is performed in all 50 states in the USA. We have partnered with the California Department of Public Health Genetic Disease Laboratory to determine whether sufficient DNA can be extracted from archived dried blood spots (DBS) for next-generation sequencing in the hopes that next-generation sequencing can play a role in NBS. We optimized the DNA extraction and sequencing library preparation protocols for residual infant DBS archived over 20 years ago and successfully obtained acceptable whole exome and whole genome sequencing data. This sequencing study using DBS DNA without whole genome amplification prior to sequencing library preparation provides evidence that properly stored residual newborn DBS are a satisfactory source of DNA for genetic studies.

Multiplexed droplet single-cell RNA-sequencing using natural genetic variation.

Droplet single-cell RNA-sequencing (dscRNA-seq) has enabled rapid, massively parallel profiling of transcriptomes. However, assessing differential expression across multiple individuals has been hampered by inefficient sample processing and technical batch effects. Here we describe a computational tool, demuxlet, that harnesses natural genetic variation to determine the sample identity of each droplet containing a single cell (singlet) and detect droplets containing two cells (doublets). These capabilities enable multiplexed dscRNA-seq experiments in which cells from unrelated individuals are pooled and captured at higher throughput than in standard workflows. Using simulated data, we show that 50 single-nucleotide polymorphisms (SNPs) per cell are sufficient to assign 97% of singlets and identify 92% of doublets in pools of up to 64 individuals. Given genotyping data for each of eight pooled samples, demuxlet correctly recovers the sample identity of >99% of singlets and identifies doublets at rates consistent with previous estimates. We apply demuxlet to assess cell-type-specific changes in gene expression in 8 pooled lupus patient samples treated with interferon (IFN)-ß and perform eQTL analysis on 23 pooled samples.

A recurrent, non-penetrant sequence variant, p.Arg266Cys in Growth/Differentiation Factor 3 (GDF3) in a female with unilateral anophthalmia and skeletal anomalies.

PURPOSE: The genetic causes of anophthalmia, microphthalmia and coloboma remain poorly understood. Missense mutations in Growth/Differentiation Factor 3 (GDF3) gene have previously been reported in patients with microphthalmia, iridial and retinal colobomas, Klippel-Feil anomaly with vertebral fusion, scoliosis, rudimentary 12th ribs and an anomalous right temporal bone. We used whole exome sequencing with a trio approach to study a female with unilateral anophthalmia, kyphoscoliosis and additional skeletal anomalies. OBSERVATIONS: Exome sequencing revealed that the proposita was heterozygous for c.796C > T, predicting p.Arg266Cys, in GDF3. Sanger sequencing confirmed the mutation and showed that the unaffected mother was heterozygous for the same missense substitution. CONCLUSIONS AND IMPORTANCE: Although transfection studies with the p.Arg266Cys mutation have shown that this amino acid substitution is likely to impair function, non-penetrance for the ocular defects was apparent in this family and has been observed in other families with sequence variants in GDF3. We conclude p.Arg266Cys and other GDF3 mutations can be non-penetrant, making pathogenicity more difficult to establish when sequence variants in this gene are present in patients with structural eye defects.

Tissue Sources for Accurate Measurement of Germline DNA Genotypes in Prostate Cancer Patients Treated With Radical Prostatectomy.

BACKGROUND: Benign tissue from a tumor-containing organ is commonly the only available source for obtaining a patient's unmutated genome for use in cancer research. While it is critical to identify histologically normal tissue that is independent of the tumor lineage, few additional considerations are applied to the choice of this material for such measurements. METHODS: Normal formalin-fixed, paraffin-embedded seminal vesicle, and urethral tissues, in addition to whole blood, were collected from 31 prostate cancer patients having undergone radical prostatectomy. Genotype concordance was evaluated for DNA from each tissue source in relation to whole blood. RESULTS: Overall, there was a greater genotype call rate for DNA derived from urethral tissue (97.0%) in comparison with patient-matched seminal vesicle tissues (95.9%, P = 0.0015). Furthermore, with reference to patient-matched whole blood, urethral samples exhibited higher genotype concordance (94.1%) than that of seminal vesicle samples (92.5%, P = 0.035). CONCLUSIONS: These findings highlight the heterogeneity between diverse sources of DNA in genotype measurement and motivate the consideration of normal tissue biases in tumor-normal analyses. Prostate 77: 425-434, 2017. © 2016 Wiley Periodicals, Inc.

Genotyping Informatics and Quality Control for 100,000 Subjects in the Genetic Epidemiology Research on Adult Health and Aging (GERA) Cohort.

The Kaiser Permanente (KP) Research Program on Genes, Environment and Health (RPGEH), in collaboration with the University of California-San Francisco, undertook genome-wide genotyping of >100,000 subjects that constitute the Genetic Epidemiology Research on Adult Health and Aging (GERA) cohort. The project, which generated >70 billion genotypes, represents the first large-scale use of the Affymetrix Axiom Genotyping Solution. Because genotyping took place over a short 14-month period, creating a near-real-time analysis pipeline for experimental assay quality control and final optimized analyses was critical. Because of the multi-ethnic nature of the cohort, four different ethnic-specific arrays were employed to enhance genome-wide coverage. All assays were performed on DNA extracted from saliva samples. To improve sample call rates and significantly increase genotype concordance, we partitioned the cohort into disjoint packages of plates with similar assay contexts. Using strict QC criteria, the overall genotyping success rate was 103,067 of 109,837 samples assayed (93.8%), with a range of 92.1-95.4% for the four different arrays. Similarly, the SNP genotyping success rate ranged from 98.1 to 99.4% across the four arrays, the variation depending mostly on how many SNPs were included as single copy vs. double copy on a particular array. The high quality and large scale of genotype data created on this cohort, in conjunction with comprehensive longitudinal data from the KP electronic health records of participants, will enable a broad range of highly powered genome-wide association studies on a diversity of traits and conditions.

DLX4 is associated with orofacial clefting and abnormal jaw development.

Cleft lip and/or palate (CL/P) are common structural birth defects in humans. We used exome sequencing to study a patient with bilateral CL/P and identified a single nucleotide deletion in the patient and her similarly affected son­c.546_546delG, predicting p.Gln183Argfs*57 in the Distal-less 4 (DLX4) gene. The sequence variant was absent from databases, predicted to be deleterious and was verified by Sanger sequencing. In mammals, there are three Dlx homeobox clusters with closely located gene pairs (Dlx1/Dlx2, Dlx3/Dlx4, Dlx5/Dlx6). In situ hybridization showed that Dlx4 was expressed in the mesenchyme of the murine palatal shelves at E12.5, prior to palate closure. Wild-type human DLX4, but not mutant DLX4_c.546delG, could activate two murine Dlx conserved regulatory elements, implying that the mutation caused haploinsufficiency. We showed that reduced DLX4 expression after short interfering RNA treatment in a human cell line resulted in significant up-regulation of DLX3, DLX5 and DLX6, with reduced expression of DLX2 and significant up-regulation of BMP4, although the increased BMP4 expression was demonstrated only in HeLa cells. We used antisense morpholino oligonucleotides to target the orthologous Danio rerio gene, dlx4b, and found reduced cranial size and abnormal cartilaginous elements. We sequenced DLX4 in 155 patients with non-syndromic CL/P and CP, but observed no sequence variants. From the published literature, Dlx1/Dlx2 double homozygous null mice and Dlx5 homozygous null mice both have clefts of the secondary palate. This first finding of a DLX4 mutation in a family with CL/P establishes DLX4 as a potential cause of human clefts.

Novel mutations in PXDN cause microphthalmia and anterior segment dysgenesis.

We used exome sequencing to study a non-consanguineous family with two children who had anterior segment dysgenesis, sclerocornea, microphthalmia, hypotonia and developmental delays. Sanger sequencing verified two Peroxidasin (PXDN) mutations in both sibs--a maternally inherited, nonsense mutation, c.1021C>T predicting p.(Arg341*), and a paternally inherited, 23-basepair deletion causing a frameshift and premature protein truncation, c.2375_2397del23, predicting p.(Leu792Hisfs*67). We re-examined exome data from 20 other patients with structural eye defects and identified two additional PXDN mutations in a sporadic male with bilateral microphthalmia, cataracts and anterior segment dysgenesis--a maternally inherited, frameshift mutation, c.1192delT, predicting p.(Tyr398Thrfs*40) and a paternally inherited, missense substitution that was predicted to be deleterious, c.947 A>C, predicting p.(Gln316Pro). Mutations in PXDN were previously reported in three families with congenital cataracts, microcornea, sclerocornea and developmental glaucoma. The gene is expressed in corneal epithelium and is secreted into the extracellular matrix. Defective peroxidasin has been shown to impair sulfilimine bond formation in collagen IV, a constituent of the basement membrane, implying that the eye defects result because of loss of basement membrane integrity in the developing eye. Our finding of a broader phenotype than previously appreciated for PXDN mutations is typical for exome-sequencing studies, which have proven to be highly effective for mutation detection in patients with atypical presentations. We conclude that PXDN sequencing should be considered in microphthalmia with anterior segment dysgenesis.

Mitochondrial DNA sequence variation and risk of pancreatic cancer.

Although the mitochondrial genome exhibits high mutation rates, common mitochondrial DNA (mtDNA) variation has not been consistently associated with pancreatic cancer. Here, we comprehensively examined mitochondrial genomic variation by sequencing the mtDNA of participants (cases = 286, controls = 283) in a San Francisco Bay Area pancreatic cancer case-control study. Five common variants were associated with pancreatic cancer at nominal statistical significance (P < 0.05) with the strongest finding for mt5460g in the ND2 gene [OR = 3.9; 95% confidence interval (CI), 1.5-10; P = 0.004] which encodes an A331T substitution. Haplogroup K was nominally associated with reduced pancreatic cancer risk (OR = 0.32; 95% CI, 0.13-0.76; P = 0.01) when compared with the most common haplogroup, H. A total of 19 haplogroup-specific rare variants yielded nominal statistically significant associations (P < 0.05) with pancreatic cancer risk, with the majority observed in genes involved in oxidative phosphorylation. Weighted-sum statistics were used to identify an aggregate effect of variants in the 22 mitochondrial tRNAs on pancreatic cancer risk (P = 0.02). While the burden of singleton variants in the HV2 and 12S RNA regions was three times higher among European haplogroup N cases than controls, the prevalence of singleton variants in ND4 and ND5 was two to three times higher among African haplogroup L cases than in controls. Together, the results of this study provide evidence that aggregated common and rare variants and the accumulation of singleton variants are important contributors to pancreatic cancer risk.

Next generation genome-wide association tool: design and coverage of a high-throughput European-optimized SNP array.

The success of genome-wide association studies has paralleled the development of efficient genotyping technologies. We describe the development of a next-generation microarray based on the new highly-efficient Affymetrix Axiom genotyping technology that we are using to genotype individuals of European ancestry from the Kaiser Permanente Research Program on Genes, Environment and Health (RPGEH). The array contains 674,517 SNPs, and provides excellent genome-wide as well as gene-based and candidate-SNP coverage. Coverage was calculated using an approach based on imputation and cross validation. Preliminary results for the first 80,301 saliva-derived DNA samples from the RPGEH demonstrate very high quality genotypes, with sample success rates above 94% and over 98% of successful samples having SNP call rates exceeding 98%. At steady state, we have produced 462 million genotypes per week for each Axiom system. The new array provides a valuable addition to the repertoire of tools for large scale genome-wide association studies.

Green technologies for room temperature nucleic acid storage.

Maintaining the long-term integrity of nucleic acids in the laboratory has traditionally required the use of freezers. However, novel nucleic acid stabilization technologies may allow for the storage of DNA and RNA at room temperature in a cost-effective, environmentally friendly manner. In this study, we evaluated two novel products for room temperature DNA storage: Biomatrica's DNA SampleMatrix technology and GenVault's GenTegra DNA technology. We compared the integrity and quality of DNA stored using these products against DNA stored in a -20 C freezer by performing downstream testing with short range PCR, long range PCR, DNA sequencing, and SNP microarrays. In addition, we tested Biomatrica's RNAstable product for its ability to preserve RNA at room temperature for use in a quantitative reverse transcription PCR assay.

Direct determination of haplotypes from single DNA molecules.

Determining the long-range haplotypes in a diploid individual is a major technical challenge. Here we report a method of molecular haplotyping by directly imaging multiple polymorphic sites on individual human DNA molecules simultaneously. We demonstrate the utility of this technology by accurately determining the haplotypes consisting of up to 16 single-nucleotide polymorphisms in genomic regions up to 50 kilobases.

Rapid DNA mapping by fluorescent single molecule detection.

DNA mapping is an important analytical tool in genomic sequencing, medical diagnostics and pathogen identification. Here we report an optical DNA mapping strategy based on direct imaging of individual DNA molecules and localization of multiple sequence motifs on the molecules. Individual genomic DNA molecules were labeled with fluorescent dyes at specific sequence motifs by the action of nicking endonuclease followed by the incorporation of dye terminators with DNA polymerase. The labeled DNA molecules were then stretched into linear form on a modified glass surface and imaged using total internal reflection fluorescence (TIRF) microscopy. By determining the positions of the fluorescent labels with respect to the DNA backbone, the distribution of the sequence motif recognized by the nicking endonuclease can be established with good accuracy, in a manner similar to reading a barcode. With this approach, we constructed a specific sequence motif map of lambda-DNA. We further demonstrated the capability of this approach to rapidly type a human adenovirus and several strains of human rhinovirus.

A simple DNA stretching method for fluorescence imaging of single DNA molecules.

Stretching or aligning DNA molecules onto a surface by means of molecular combing techniques is one of the critical steps in single DNA molecule analysis. However, many of the current studies have focused on lambda-DNA, or other large DNA molecules. There are very few studies on stretching methodologies for DNA molecules generated via PCR (typically smaller than 20 kb). Here we describe a simple method of stretching DNA molecules up to 18 kb in size on a modified glass surface. The very low background fluorescence allows efficient detection of single fluorescent dye labels incorporated into the stretched DNA molecules.