High-definition spatial transcriptomics for in situ tissue profiling.

Spatial and molecular characteristics determine tissue function, yet high-resolution methods to capture both concurrently are lacking. Here, we developed high-definition spatial transcriptomics, which captures RNA from histological tissue sections on a dense, spatially barcoded bead array. Each experiment recovers several hundred thousand transcript-coupled spatial barcodes at 2-µm resolution, as demonstrated in mouse brain and primary breast cancer. This opens the way to high-resolution spatial analysis of cells and tissues.

Sequences of 95 human MHC haplotypes reveal extreme coding variation in genes other than highly polymorphic HLA class I and II.

The most polymorphic part of the human genome, the MHC, encodes over 160 proteins of diverse function. Half of them, including the HLA class I and II genes, are directly involved in immune responses. Consequently, the MHC region strongly associates with numerous diseases and clinical therapies. Notoriously, the MHC region has been intractable to high-throughput analysis at complete sequence resolution, and current reference haplotypes are inadequate for large-scale studies. To address these challenges, we developed a method that specifically captures and sequences the 4.8-Mbp MHC region from genomic DNA. For 95 MHC homozygous cell lines we assembled, de novo, a set of high-fidelity contigs and a sequence scaffold, representing a mean 98% of the target region. Included are six alternative MHC reference sequences of the human genome that we completed and refined. Characterization of the sequence and structural diversity of the MHC region shows the approach accurately determines the sequences of the highly polymorphic HLA class I and HLA class II genes and the complex structural diversity of complement factor C4A/C4B It has also uncovered extensive and unexpected diversity in other MHC genes; an example is MUC22, which encodes a lung mucin and exhibits more coding sequence alleles than any HLA class I or II gene studied here. More than 60% of the coding sequence alleles analyzed were previously uncharacterized. We have created a substantial database of robust reference MHC haplotype sequences that will enable future population scale studies of this complicated and clinically important region of the human genome.

Defining KIR and HLA Class I Genotypes at Highest Resolution via High-Throughput Sequencing.

The physiological functions of natural killer (NK) cells in human immunity and reproduction depend upon diverse interactions between killer cell immunoglobulin-like receptors (KIRs) and their HLA class I ligands: HLA-A, HLA-B, and HLA-C. The genomic regions containing the KIR and HLA class I genes are unlinked, structurally complex, and highly polymorphic. They are also strongly associated with a wide spectrum of diseases, including infections, autoimmune disorders, cancers, and pregnancy disorders, as well as the efficacy of transplantation and other immunotherapies. To facilitate study of these extraordinary genes, we developed a method that captures, sequences, and analyzes the 13 KIR genes and HLA-A, HLA-B, and HLA-C from genomic DNA. We also devised a bioinformatics pipeline that attributes sequencing reads to specific KIR genes, determines copy number by read depth, and calls high-resolution genotypes for each KIR gene. We validated this method by using DNA from well-characterized cell lines, comparing it to established methods of HLA and KIR genotyping, and determining KIR genotypes from 1000 Genomes sequence data. This identified 116 previously uncharacterized KIR alleles, which were all demonstrated to be authentic by sequencing from source DNA via standard methods. Analysis of just two KIR genes showed that 22% of the 1000 Genomes individuals have a previously uncharacterized allele or a structural variant. The method we describe is suited to the large-scale analyses that are needed for characterizing human populations and defining the precise HLA and KIR factors associated with disease. The methods are applicable to other highly polymorphic genes.

In vitro, long-range sequence information for de novo genome assembly via transposase contiguity.

We describe a method that exploits contiguity preserving transposase sequencing (CPT-seq) to facilitate the scaffolding of de novo genome assemblies. CPT-seq is an entirely in vitro means of generating libraries comprised of 9216 indexed pools, each of which contains thousands of sparsely sequenced long fragments ranging from 5 kilobases to > 1 megabase. These pools are "subhaploid," in that the lengths of fragments contained in each pool sums to ∼5% to 10% of the full genome. The scaffolding approach described here, termed fragScaff, leverages coincidences between the content of different pools as a source of contiguity information. Specifically, CPT-seq data is mapped to a de novo genome assembly, followed by the identification of pairs of contigs or scaffolds whose ends disproportionately co-occur in the same indexed pools, consistent with true adjacency in the genome. Such candidate "joins" are used to construct a graph, which is then resolved by a minimum spanning tree. As a proof-of-concept, we apply CPT-seq and fragScaff to substantially boost the contiguity of de novo assemblies of the human, mouse, and fly genomes, increasing the scaffold N50 of de novo assemblies by eight- to 57-fold with high accuracy. We also demonstrate that fragScaff is complementary to Hi-C-based contact probability maps, providing midrange contiguity to support robust, accurate chromosome-scale de novo genome assemblies without the need for laborious in vivo cloning steps. Finally, we demonstrate CPT-seq as a means of anchoring unplaced novel human contigs to the reference genome as well as for detecting misassembled sequences.

Whole-genome haplotyping by dilution, amplification, and sequencing.

Standard whole-genome genotyping technologies are unable to determine haplotypes. Here we describe a method for rapid and cost-effective long-range haplotyping. Genomic DNA is diluted and distributed into multiple aliquots such that each aliquot receives a fraction of a haploid copy. The DNA template in each aliquot is amplified by multiple displacement amplification, converted into barcoded sequencing libraries using Nextera technology, and sequenced in multiplexed pools. To assess the performance of our method, we combined two male genomic DNA samples at equal ratios, resulting in a sample with diploid X chromosomes with known haplotypes. Pools of the multiplexed sequencing libraries were subjected to targeted pull-down of a 1-Mb contiguous region of the X-chromosome Duchenne muscular dystrophy gene. We were able to phase the Duchenne muscular dystrophy region into two contiguous haplotype blocks with a mean length of 494 kb. The haplotypes showed 99% agreement with the consensus base calls made by sequencing the individual DNAs. We subsequently used the strategy to haplotype two human genomes. Standard genomic sequencing to identify all heterozygous SNPs in the sample was combined with dilution-amplification-based sequencing data to resolve the phase of identified heterozygous SNPs. Using this procedure, we were able to phase >95% of the heterozygous SNPs from the diploid sequence data. The N50 for a Yoruba male DNA was 702 kb whereas the N50 for a European female DNA was 358 kb. Therefore, the strategy described here is suitable for haplotyping of a set of targeted regions as well as of the entire genome.

Unusual selection on the KIR3DL1/S1 natural killer cell receptor in Africans.

Interactions of killer cell immunoglobulin-like receptors (KIRs) with major histocompatibility complex (MHC) class I ligands diversify natural killer cell responses to infection. By analyzing sequence variation in diverse human populations, we show that the KIR3DL1/S1 locus encodes two lineages of polymorphic inhibitory KIR3DL1 allotypes that recognize Bw4 epitopes of protein">HLA-A and HLA-B and one lineage of conserved activating KIR3DS1 allotypes, also implicated in Bw4 recognition. Balancing selection has maintained these three lineages for over 3 million years. Variation was selected at D1 and D2 domain residues that contact HLA class I and at two sites on D0, the domain that enhances the binding of KIR3D to HLA class I. HLA-B variants that gained Bw4 through interallelic microconversion are also products of selection. A worldwide comparison uncovers unusual KIR3DL1/S1 evolution in modern sub-Saharan Africans. Balancing selection is weak and confined to D0, KIR3DS1 is rare and KIR3DL1 allotypes with similar binding sites predominate. Natural killer cells express the dominant KIR3DL1 at a high frequency and with high surface density, providing strong responses to cells perturbed in Bw4 expression.

Very long haplotype tracts characterized at high resolution from HLA homozygous cell lines.

The HLA region of chromosome 6 contains the most polymorphic genes in humans. Spanning ~5 Mbp the densely packed region encompasses approximately 175 expressed genes including the highly polymorphic HLA class I and II loci. Most of the other genes and functional elements are also polymorphic, and many of them are directly implicated in immune function or immune-related disease. For these reasons, this complex genomic region is subject to intense scrutiny by researchers with the common goal of aiding further understanding and diagnoses of multiple immune-related diseases and syndromes. To aid assay development and characterization of the classical loci, a panel of cell lines partially or fully homozygous for HLA class I and II was assembled over time by the International Histocompatibility Working Group (IHWG). Containing a minimum of 88 unique HLA haplotypes, we show that this panel represents a significant proportion of European HLA allelic and haplotype diversity (60-95 %). Using a high-density whole genome array that includes 13,331 HLA region SNPs, we analyzed 99 IHWG cells to map the coordinates of the homozygous tracts at a fine scale. The mean homozygous tract length within chromosome 6 from these individuals is 21 Mbp. Within HLA, the mean haplotype length is 4.3 Mbp, and 65 % of the cell lines were shown to be homozygous throughout the entire region. In addition, four cell lines are homozygous throughout the complex KIR region of chromosome 19 (~250 kbp). The data we describe will provide a valuable resource for characterizing haplotypes, designing and refining imputation algorithms and developing assay controls.

Highly parallel oligonucleotide purification and functionalization using reversible chemistry.

We have developed a cost-effective, highly parallel method for purification and functionalization of 5'-labeled oligonucleotides. The approach is based on 5'-hexa-His phase tag purification, followed by exchange of the hexa-His tag for a functional group using reversible reaction chemistry. These methods are suitable for large-scale (micromole to millimole) production of oligonucleotides and are amenable to highly parallel processing of many oligonucleotides individually or in high complexity pools. Examples of the preparation of 5'-biotin, 95-mer, oligonucleotide pools of >40K complexity at micromole scale are shown. These pools are prepared in up to ~16% yield and 90-99% purity. Approaches for using this method in other applications are also discussed.

Deep sequencing of the hepatitis B virus in hepatocellular carcinoma patients reveals enriched integration events, structural alterations and sequence variations.

Chronic hepatitis B virus (HBV) infection is epidemiologically associated with hepatocellular carcinoma (HCC), but its role in HCC remains poorly understood due to technological limitations. In this study, we systematically characterize HBV in HCC patients. HBV sequences were enriched from 48 HCC patients using an oligo-bead-based strategy, pooled together and sequenced using the FLX-Genome-Sequencer. In the tumors, preferential integration of HBV into promoters of genes (P < 0.001) and significant enrichment of integration into chromosome 10 (P < 0.01) were observed. Integration into chromosome 10 was significantly associated with poorly differentiated tumors (P < 0.05). Notably, in the tumors, recurrent integration into the promoter of the human telomerase reverse transcriptase (TERT) gene was found to correlate with increased TERT expression. The preferred region within the HBV genome involved in integration and viral structural alteration is at the 3'-end of hepatitis B virus X protein (HBx), where viral replication/transcription initiates. Upon integration, the 3'-end of the HBx is often deleted. HBx-human chimeric transcripts, the most common type of chimeric transcripts, can be expressed as chimeric proteins. Sequence variation resulting in non-conservative amino acid substitutions are commonly observed in HBV genome. This study highlights HBV as highly mutable in HCC patients with preferential regions within the host and virus genome for HBV integration/structural alterations.

Although divergent in residues of the peptide binding site, conserved chimpanzee Patr-AL and polymorphic human HLA-A*02 have overlapping peptide-binding repertoires.

Patr-AL is an expressed, non-polymorphic MHC class I gene carried by ∼50% of chimpanzee MHC haplotypes. Comparing Patr-AL(+) and Patr-AL(-) haplotypes showed Patr-AL defines a unique 125-kb genomic block flanked by blocks containing classical Patr-A and pseudogene Patr-H. Orthologous to Patr-AL are polymorphic orangutan Popy-A and the 5' part of human pseudogene HLA-Y, carried by ∼10% of HLA haplotypes. Thus, the AL gene alternatively evolved in these closely related species to become classical, nonclassical, and nonfunctional. Although differing by 30 aa substitutions in the peptide-binding α(1) and α(2) domains, Patr-AL and HLA-A*0201 bind overlapping repertoires of peptides; the overlap being comparable with that between the A*0201 and A*0207 subtypes differing by one substitution. Patr-AL thus has the A02 supertypic peptide-binding specificity. Patr-AL and HLA-A*0201 have similar three-dimensional structures, binding peptides in similar conformation. Although comparable in size and shape, the B and F specificity pockets of Patr-AL and HLA-A*0201 differ in both their constituent residues and contacts with peptide anchors. Uniquely shared by Patr-AL, HLA-A*0201, and other members of the A02 supertype are the absence of serine at position 9 in the B pocket and the presence of tyrosine at position 116 in the F pocket. Distinguishing Patr-AL from HLA-A*02 is an unusually electropositive upper face on the α(2) helix. Stimulating PBMCs from Patr-AL(-) chimpanzees with B cells expressing Patr-AL produced potent alloreactive CD8 T cells with specificity for Patr-AL and no cross-reactivity toward other MHC class I molecules, including HLA-A*02. In contrast, PBMCs from Patr-AL(+) chimpanzees are tolerant of Patr-AL.

LTBP2 mutations cause Weill-Marchesani and Weill-Marchesani-like syndrome and affect disruptions in the extracellular matrix.

Latent transforming growth factor (TGF) beta-binding protein 2 (LTBP2) is an extracellular matrix (ECM) protein that associates with fibrillin-1 containing microfibrils. Various factors prompted considering LTBP2 in the etiology of isolated ectopia lentis and associated conditions such as Weill-Marchesani syndrome (WMS) and Marfan syndrome (MFS). LTBP2 was screened in 30 unrelated Iranian patients. Mutations were found only in one WMS proband and one MFS proband. Homozygous c.3529G>A (p.Val1177Met) was shown to cause autosomal recessive WMS or WM-like syndrome by several approaches, including homozygosity mapping. Light, fluorescent, and electron microscopy evidenced disruptions of the microfibrillar network in the ECM of the proband's skin. In conjunction with recent findings regarding other ECM proteins, the results presented strongly support the contention that anomalies in WMS patients are due to disruptions in the ECM. Heterozygous c.1642C >T (p.Arg548*) possibly contributed to MFS-related phenotypes, including ocular manifestations, mitral valve prolapse, and pectus excavatum, but was not cause of MFS.

Whole-genome sequencing of the efficient industrial fuel-ethanol fermentative Saccharomyces cerevisiae strain CAT-1.

The Saccharomyces cerevisiae strains widely used for industrial fuel-ethanol production have been developed by selection, but their underlying beneficial genetic polymorphisms remain unknown. Here, we report the draft whole-genome sequence of the S. cerevisiae strain CAT-1, which is a dominant fuel-ethanol fermentative strain from the sugarcane industry in Brazil. Our results indicate that strain CAT-1 is a highly heterozygous diploid yeast strain, and the ~12-Mb genome of CAT-1, when compared with the reference S228c genome, contains ~36,000 homozygous and ~30,000 heterozygous single nucleotide polymorphisms, exhibiting an uneven distribution among chromosomes due to large genomic regions of loss of heterozygosity (LOH). In total, 58 % of the 6,652 predicted protein-coding genes of the CAT-1 genome constitute different alleles when compared with the genes present in the reference S288c genome. The CAT-1 genome contains a reduced number of transposable elements, as well as several gene deletions and duplications, especially at telomeric regions, some correlated with several of the physiological characteristics of this industrial fuel-ethanol strain. Phylogenetic analyses revealed that some genes were likely associated with traits important for bioethanol production. Identifying and characterizing the allelic variations controlling traits relevant to industrial fermentation should provide the basis for a forward genetics approach for developing better fermenting yeast strains.

Six RNA viruses and forty-one hosts: viral small RNAs and modulation of small RNA repertoires in vertebrate and invertebrate systems.

We have used multiplexed high-throughput sequencing to characterize changes in small RNA populations that occur during viral infection in animal cells. Small RNA-based mechanisms such as RNA interference (RNAi) have been shown in plant and invertebrate systems to play a key role in host responses to viral infection. Although homologs of the key RNAi effector pathways are present in mammalian cells, and can launch an RNAi-mediated degradation of experimentally targeted mRNAs, any role for such responses in mammalian host-virus interactions remains to be characterized. Six different viruses were examined in 41 experimentally susceptible and resistant host systems. We identified virus-derived small RNAs (vsRNAs) from all six viruses, with total abundance varying from "vanishingly rare" (less than 0.1% of cellular small RNA) to highly abundant (comparable to abundant micro-RNAs "miRNAs"). In addition to the appearance of vsRNAs during infection, we saw a number of specific changes in host miRNA profiles. For several infection models investigated in more detail, the RNAi and Interferon pathways modulated the abundance of vsRNAs. We also found evidence for populations of vsRNAs that exist as duplexed siRNAs with zero to three nucleotide 3' overhangs. Using populations of cells carrying a Hepatitis C replicon, we observed strand-selective loading of siRNAs onto Argonaute complexes. These experiments define vsRNAs as one possible component of the interplay between animal viruses and their hosts.

Non-housekeeping genes expressed in human trabecular meshwork cell cultures.

PURPOSE: To identify non-housekeeping genes definitively expressed in the human trabecular meshwork (TM). METHODS: Microarray gene expression data on TM cultured cells from four studies were compared. Genes that were queried in at least three studies and assessed to be expressed in at least three studies were considered definitively expressed genes of the human TM. Housekeeping genes were removed from this set of genes. The non-housekeeping TM gene profile was analyzed for pathway enrichment and microRNA targeting, using bioinformatics tools. The results were compared with results of previous non-array based studies. RESULTS: Nine hundred and sixty-two genes were identified as non-housekeeping TM expressed genes. Analysis of these by Kyoto Encyclopedia of Genes and Genomes led to identification of two enriched biologic pathways that achieved a highly significant Bonferroni p-value (p≤0.01): focal adhesion and extracellular matrix (ECM)-receptor interaction. Many of the genes were previously implicated in TM-related functions and the TM-associated disease glaucoma; however, some are novel. MicroRNAs known to be expressed in the trabecular meshwork were predicted to target some of the genes. Ten genes identified here, ALDH1A1 (aldehyde dehydrogenase 1 family, member A1), CDH11 (cadherin 11, type 2, OB-cadherin), CXCR7 (chemokine (C-X-C motif) receptor 7), CHI3L1 (chitinase 3-like 1), FGF2 (fibroblast growth factor 2), GNG11 (guanine nucleotide binding protein [G protein], gamma 11), IGFBP5 (insulin-like growth factor binding protein 5), PTPRM (protein tyrosine phosphatase, receptor type, M), RGS5 (regulator of G-protein signaling 5), and TUSC3 (tumor suppressor candidate 3), were also reported as TM expressed genes in three earlier non-microarray based studies. CONCLUSIONS: A transcriptome consisting of 962 non-housekeeping genes definitively expressed in the human TM was identified. Multiple genes and microRNAs are proposed for further study for a better understanding of TM physiology.

pfSNP: An integrated potentially functional SNP resource that facilitates hypotheses generation through knowledge syntheses.

Currently, >14,000,000 single nucleotide polymorphisms (SNPs) are reported. Identifying phenotype-affecting SNPs among these many SNPs pose significant challenges. Although several Web resources are available that can inform about the functionality of SNPs, these resources are mainly annotation databases and are not very comprehensive. In this article, we present a comprehensive, well-annotated, integrated pfSNP (potentially functional SNPs) Web resource (http://pfs.nus.edu.sg/), which is aimed to facilitate better hypothesis generation through knowledge syntheses mediated by better data integration and a user-friendly Web interface. pfSNP integrates >40 different algorithms/resources to interrogate >14,000,000 SNPs from the dbSNP database for SNPs of potential functional significance based on previous published reports, inferred potential functionality from genetic approaches as well as predicted potential functionality from sequence motifs. Its query interface has the user-friendly "auto-complete, prompt-as-you-type" feature and is highly customizable, facilitating different combination of queries using Boolean-logic. Additionally, to facilitate better understanding of the results and aid in hypotheses generation, gene/pathway-level information with text clouds highlighting enriched tissues/pathways as well as detailed-related information are also provided on the results page. Hence, the pfSNP resource will be of great interest to scientists focusing on association studies as well as those interested to experimentally address the functionality of SNPs.

Loss of function mutations in the gene encoding latent transforming growth factor beta binding protein 2, LTBP2, cause primary congenital glaucoma.

Glaucoma is a heterogeneous group of optic neuropathies that manifests by optic nerve head cupping or degeneration of the optic nerve, resulting in a specific pattern of visual field loss. Glaucoma leads to blindness if left untreated, and is considered the second leading cause of blindness worldwide. The subgroup primary congenital glaucoma (PCG) is characterized by an anatomical defect in the trabecular meshwork, and age at onset in the neonatal or infantile period. It is the most severe form of glaucoma. CYP1B1 was the first gene genetically linked to PCG, and CYP1B1 mutations are the cause of disease in 20-100% of patients in different populations. Here, we report that LTBP2 encoding latent transforming growth factor beta binding protein 2 is a PCG causing gene, confirming results recently reported. A disease-associated locus on chromosome 14 was identified by performing whole genome autozygosity mapping in Iranian PCG families using high density single nucleotide polymorphism chips, and two disease-segregating loss of function mutations in LTBP2, p.Ser472fsX3 and p.Tyr1793fsX55, were observed in two families while sequencing candidate genes in the locus. The p.Tyr1793fsX55 mutation affects an amino acid close to the C-terminal of the encoded protein. Subsequently, LTBP2 expression was shown in human eyes, including the trabecular meshwork and ciliary processes that are thought to be relevant to the etiology of PCG.

Genome-wide linkage analysis of a Parkinsonian-pyramidal syndrome pedigree by 500 K SNP arrays.

Robust SNP genotyping technologies and data analysis programs have encouraged researchers in recent years to use SNPs for linkage studies. Platforms used to date have been 10 K chip arrays, but the possible value of interrogating SNPs at higher densities has been considered. Here, we present a genome-wide linkage analysis by means of a 500 K SNP platform. The analysis was done on a large pedigree affected with Parkinsonian-pyramidal syndrome (PPS), and the results showed linkage to chromosome 22. Sequencing of candidate genes revealed a disease-associated homozygous variation (R378G) in FBXO7. FBXO7 codes for a member of the F-box family of proteins, all of which may have a role in the ubiquitin-proteosome protein-degradation pathway. This pathway has been implicated in various neurodegenerative diseases, and identification of FBXO7 as the causative gene of PPS is expected to shed new light on its role. The performance of the array was assessed and systematic analysis of effects of SNP density reduction was performed with the real experimental data. Our results suggest that linkage in our pedigree may have been missed had we used chips containing less than 100,000 SNPs across the genome.

Effect of PITX2 knockdown on transcriptome of primary human trabecular meshwork cell cultures.

PURPOSE: To identify genes whose expressions in primary human trabecular meshwork (TM) cell cultures are affected by the transcription factor pituitary homeobox 2 (PITX2) and to identify genes that may have roles in glaucoma. Known glaucoma causing genes account for disease in a small fraction of patients, and we aimed at identification of other genes that may have subtle and accumulative effects not easily identifiable by a genetic approach. METHODS: Expression profiles derived using microarrays were compared between TM control cells and cells treated with PITX2 siRNAs using three protocols so as to minimize false positive and negative results. The first protocol was based on the commonly used B statistic. The second and third protocols were based on fold change in expression. The second protocol used a threshold of at least 2 fold change in expression, whereas the third protocol used ranking in fold change without setting a threshold. The likelihood of a selected gene being a true positive was considered to correlate with the number of protocols by which it was selected. By considering all genes that were selected by at least one protocol, the likelihood of false negatives was expected to decrease. Effects on a subset of selected genes were verified by real time PCR, western blots, and immunocytochemistry. Effects on ALDH1A1, were further pursued because its protein product, aldehyde dehydrogenase 1 family, member A1, has roles in oxidative stress and because oxidative stress is known to be relevant to the etiology of glaucoma. RESULTS: The expression level of 41 genes was assessed by to be possibly affected by PITX2 knockdown. Twenty one genes were down-regulated and twenty were upregulated. The expression of five genes was assessed to be altered by all three analysis protocols. The five genes were DIRAS3 (DIRAS family, GTP-binding RAS-like 3), CXCL6 (chemokine (C-X-C motif) ligand 6), SAMD5 (sterile alpha motif domain containing 5), CBFB (core-binding factor, beta subunit), and MEIS2 (meis homeobox 2). Real time PCR experiments verified results on a subset of genes tested. Notably, the results were also confirmed in two independent TMs. Effects on CXCL6 and ALDH1A1 were also confirmed by western blots, and effects on ALDH1A1 were further shown by immunocytochemistry. Data consistent with PITX2 involvement in ALDH1A1 mediated response to oxidative stress were presented. CONCLUSIONS: Bioinformatics tools revealed that the genes identified affect functions and pathways relevant to glaucoma. Involvement of PITX2 in expression of some of the genes and in some of the pathways is being reported here for the first time. As many of the genes identified have not been studied vis-à-vis glaucoma, we feel they introduce new candidates for understanding this devastating disease.

PRKN, DJ-1, and PINK1 screening identifies novel splice site mutation in PRKN and two novel DJ-1 mutations.

We present results of mutation screening of PRKN gene in 93 Iranian Parkinson's disease (PD) patients with average age at onset (AAO) of 42.2 years. The gene was screened by direct sequencing and by a semi-quantitative PCR protocol for detection of sequence rearrangements. Heterozygous rearrangements were tested by reverse transcription-polymerase chain reaction (RT-PCR). Nine different PRKN mutations were found. One of these, IVS9+1G>A, affects splicing and is novel. Two mutated PRKN alleles were observed in each of 6 patients whose average AAO was 25.7 years. Only 1 patient carried a single mutated allele and his AAO was 41 years. Among patients with AAO of <30 years, 31.3% had two mutated alleles, while only 2.6% with AAO of >30 years carried a PRKN mutation. Analysis of PRKN by RT-PCR led to identification of a novel exon expressed in leukocytes of control and PD individuals. The alternatively spliced transcript if translated would code a protein without a RING Finger 2 domain. Its functional relevance remains to be shown. DJ-I and PINK1 were also screened. Two novel DJ-1 mutations, c.91-2A>G affecting splicing and c.319G>C causing Ala107Pro, were observed among patients with AAO of <31 years, suggesting that PD in a high fraction (>12%) of this group of Iranian patients may be due to mutations in DJ-1. Mutations in PINK1 were not observed. Our results complement previous findings on LRRK2 mutations among Iranian PD patients.