The present disclosure relates to compositions and methods for the diagnosis of schizophrenia. In particular, the instant disclosure is directed to identification of novel copy number variants of sequences associated with the VIPR2 gene, including certain micro-duplications and triplications, and correlation of these copy number variants with schizophrenia.
Major progress has been made in understanding the genetic architecture of schizophrenia and other neuropsychiatric disorders. Studies have established the role of rare copy number variants (CNVs) in the etiology of schizophrenia22-27. Substantial risk for schizophrenia is conferred by large (>500 kb) CNVs at several loci, including microdeletions at 1q21.123, 3q2928, 15q13.323 and 22q11.229 and microduplication at 16p11.230. All of the loci described above are hotspots for genomic rearrangement, where local segmental duplication architecture promotes frequent and nearly identical rearrangements by non-allelic homologous recombination (NAHR). Because of the high structural mutation rates at these loci, the strong phenotypic effects of the causal variants, and the excellent power of most array platforms to detect such large CNVs, these genomic hotspots were the first to be detected in studies of CNVs in schizophrenia. However, these CNVs collectively account for a small fraction (2-4%) of cases26 and do not explain most of the increased CNV burden in schizophrenia reported in earlier studies22, 23.
Since most of the genome lacks the duplication architecture of the NAHR hotspots described above and because a variety of mutational mechanisms can give rise to structural rearrangements, causal variants in other regions of the genome may consist of CNVs that are individually rarer and smaller (less than 500 KB) than those arising at NAHR hotspots. For example, microdeletions of the gene Neurexin-1 (NRXN1) are highly enriched in autism and schizophrenia31-33 and consist of overlapping deletions with non-recurrent breakpoints. NRXN1 deletions are not flanked by segmental duplications, and may occur by different mutational mechanisms such as non-homologous end joining (NHEJ) or DNA replication-mediated rearrangement
In light of the fact that the aforementioned CNVs do not explain most of the increased CNV burden in schizophrenia, it would be desirable to identify new compositions and methods for facilitating the diagnosis of schizophrenia, particularly those that are associated with the identification of novel CNVs, such as those that are found outside of NAHR hotspots.
The present disclosure relates to compositions and methods for the diagnosis of schizophrenia. In particular, the instant disclosure is directed to identification of novel copy number variants of sequences associated with the VIPR2 gene (VIPR2 CNV).
In certain embodiments, the VIPR2 CNV comprises a duplication or triplication of nucleic acid of chromosome 7.
In certain embodiments, the VIPR2 CNV comprises nucleic acid in region 7q36 of human chromosome 7.
In certain embodiments, the VIPR2 CNV comprises a region or segment of nucleic acid that is within 89 kb of the transcriptional start site of the VIPR2 gene.
In certain embodiments, the VIPR2 CNV comprises a region or segment of nucleic acid that is located, at least in part, in region 158,448,321-158,810,016; 158,731,401-158,810,016 or 158,448,321-158,605,936 of human chromosome 7, as described by the human genome build NCBI36/hg18 (produced by the International Human Genome Sequencing Consortium).
In certain embodiments, the VIPR2 CNV comprises a region or segment of exon 3 and/or exon 4 of the VIPR2 gene.
In certain embodiments, the VIPR2 CNV is associated with a second disorder or condition, for example, a psychiatric condition, disorder or phenotype. In certain embodiments, the VIPR2 CNV is associated with a pediatric neurodevelopment disorder. In certain embodiments, the VIPR2 CNV is associated with autism.
The present disclosure also relates to compositions and methods for the diagnosis of schizophrenia in a subject by detecting an increase in expression of a VIPR2 gene in a sample from the subject compared to VIPR2 expression in a sample from a control subject who does not have schizophrenia.
The present disclosure also relates to compositions and methods for the diagnosis of schizophrenia in a subject by detecting an increase in the level of cyclic-AMP in response to an agonist of VIPR2 such as, but not limited to, vasoactive intestinal peptide, or BAY55-9837, in a sample from the subject compared to the level of cyclic-AMP in a sample from a control subject who does not have schizophrenia. As one non-limiting example, the sample may comprise lymphoid or lymphoblastoid cells, which are then treated with the VIPR2 agonist prior to measurement of cAMP levels.
In certain embodiments, a VIPR2 CNV or VIPR2 expression product, can be detected through the use of, for example, the polymerase chain reaction (PCR), quantitative PCR, nucleic acid sequencing, nucleic acid microarray analysis or immunological detection.
In certain embodiments, the present disclosure provides for methods of treating a subject that has a VIPR2 CNV, an increased level of VIPR2 gene expression compared to a non-schizophrenic control, or an increased level of c-AMP in response to VIPR2 agonist compared to a non-schizophrenic control, wherein an agent is administered to the subject in an amount effective to decrease the level of VIPR2 in a sample from the subject.
In certain embodiments the agent inhibits the function of VIPR2 protein or reduces the level of functional VIPR2 protein present in a sample from the subject.
In certain embodiments, the agent is a VIPR2 protein antagonist or inhibitor.
In certain embodiments, the agent is an antisense molecule, RNAi molecule or siRNA molecule. In certain embodiments, the antisense, RNAi or siRNA molecule is complementary to a segment or region of a VIPR2 mRNA transcript. In certain embodiments, the antisense, RNAi or siRNA molecule hybridizes to and inhibits or reduces translation of VIPR2 mRNA.
In certain embodiments, the present disclosure provides for a kit for detecting at least one VIPR2 CNV, or VIPR2 expression product, wherein the kit comprises a plurality of oligonucleotide primers, each of which is capable of specifically hybridizing to genomic DNA associated with a VIPR2 CNV, or VIPR2 expression product, for example, genomic VIPR2 nucleic acid or VIPR2 mRNA.
In certain embodiments, the present disclosure provides for a kit for detecting VIPR2 expression, or cyclic-AMP level, wherein the kit comprises at least one antibody specific for a VIPR2 expression product, for example, VIPR2 protein, or cyclic-AMP, and a means of detecting the antibody when it is bound to the VIPR2 expression product or cyclic-AMP.
FIG. 1A-1-1H Detection and validation of microduplications and triplications of 7q36.3. (a) Map of CNVs detected in the primary and secondary cohorts from the UCSC genome browser. (b) Plots of probe intensity ratios for 16 CNVs detected in the primary and MGS datasets. All are cases, with the exception of two controls which are indicated with an asterisk (*). Regions with estimated copy numbers of 2, 3 and 4 are highlighted in gray, blue and green respectively. Locations of four Sequenom validation assays are shown (dashed lines). (c-f) CNV genotypes were confirmed by MeZOD cluster plots of probe intensity ratios of the proximal and distal regions and in the primary dataset (c and d respectively) and secondary dataset (e and f respectively). Absolute copy numbers were confirmed for duplications and triplications of the proximal (g) and distal (h) regions by Sequenom MASSarray genotyping.
FIGS. 5-1-5-3 Map of CNVs Detected in the 7q36.3 Region after Relaxing Stringent Filters on CNV Confidence Score and CNV Length. No additional CNVs were detected after eliminating the filter on the CNV confidence score. Additional small CNVs were detected in cases and controls after relaxing filters on CNV length (shown in green). The process of identifying of target regions in the primary cohort and testing association in the secondary cohort was repeated using the low-stringency CNV call set below. Target regions, as re-defined in the low-stringency call set, are shown in black. Association of target regions in the secondary cohort was statistically significant (permutation P-value=1.5*10−5, Peak OR=17.1 [3.3, Infinity]).
FIGS. 6-1-6-2 Fine Mapping of CNVs in MGS Subjects Using the NimbleGen HD2 Array. Fine mapping of 7q36.3 duplications was performed by re-scanning several individuals from the MGS study using the NimbleGen HD2 platform. NimbleGen results (right panel) were presented in
FIGS. 8-1-8-7 Evaluating sensitivity of CNV detection in MGS cases and MGS controls based on concordance of segmentation calls with genotyping calls for 7 large copy number polymorphsisms (CNPs). For seven large (>100 Kb) CNPs reported in HapMap phase 3, clusters of median Z-scores were assigned genotypes as illustrated. The cluster that overlapped with a median Z-score of zero was assigned a genotype of “normal”, and distinct clusters where all Z-scores were greater than or less than 0 were assigned genotypes of “gain” or “loss”, respectively. Genotypes obtained in this manner were then compared with segmentation calls, and sensitivity was defined as the average fraction of genotyped gains or losses that were detected by segmentation per individual. By this approach, a reduced sensitivity to detect CNVs in controls was not observed (see Table 7).
The specification further incorporates by reference the Sequence Listing submitted via EFS on Jul. 3, 2013. Pursuant to 37 C.F.R. §1.52(e)(5), the Sequence Listing text file, identified as 0700504364ConSeqlist.txt, is 5,081 bytes and was created on Jul. 3, 2013. The Sequence Listing, electronically filed herewith, does not extend beyond the scope of the specification and thus does not contain new matter.
The present disclosure relates to compositions and methods for the diagnosis of schizophrenia. In particular, the instant application is directed to identification of novel copy number variants (CNV) of sequences associated with the VIPR2 gene, including certain micro-duplications and triplications, and correlation of these copy number variants with schizophrenia.
The present application is based at least in part on the identification of an association of CNVs of chromosome 7q36.3 with schizophrenia. For example, a genome-wide analysis of CNV in a primary cohort of 802 patients and 742 controls identified 114 genomic regions where multiple overlapping CNVs were detected exclusively in cases. These 114 regions were interrogated in a second series of 7,488 cases and 6,689 controls. Statistically significant associations were detected for copy number gains at 7q36.3 (P=4.0×10-5, OR=16.14 [3.06, ∞]). Microduplications with variable breakpoints occurred within a 362 kb region of 7q36.3, and were detected in 29 of 8,290 (0.35%) patients versus two of 7,431 (0.03%) controls in the combined sample (P=5.7×10−7, odds ratio (OR)=14.1 [3.5, 123.9]). All duplications at 7q36.3 overlapped or were located within 89 Kb upstream of the vasoactive intestinal peptide receptor VIPR2, and four patients carried partial triplications of this gene. VIPR2 transcription and cyclic-AMP signaling were significantly increased in cultured lymphocytes from patients with microduplications of 7q36.3.
For clarity and not by way of limitation, this detailed description is divided into the following sub-portions:
(i) VIPR2 CNVs;
(ii) Methods of Diagnosing Schizophrenia;
(iii) Methods of Treating Schizophrenia; and
(iv) Kits.
The present disclosure provides methods, compositions, and kits for diagnosing schizophrenia. Methods for detecting VIPR2 CNV are useful for detecting and/or diagnosing those who have or are at risk of developing schizophrenia. These methods and compositions will be useful in the diagnosis and treatment of persons with schizophrenia.
The VIPR2 gene encodes a Vasoactive Intestinal Peptide (VIP) Receptor 2 (VPAC2), a G protein-coupled receptor that is expressed in variety of tissues including, for example, brain, suprachiasmatic nucleus, hippocampus, amygdala, and hypothalamus. In certain embodiments, VIPR2 is a human VIPR2 gene, for example, a human VIPR2 gene described by GenBank accession number NM—003382. In certain embodiments, the human VIPR2 gene encodes an amino acid sequence described by GenBank accession number NP—003373.
In certain embodiments, VPAC2 binds VIP. In certain embodiments, binding of VPAC2 to VIP activates cyclic-AMP signaling and protein kinase A (PKA). In certain embodiments, VPAC2 functions to regulate synaptic transmission in the hippocampus. In certain embodiments, VPAC2 functions to promote proliferation of neural progenitor cells in the dentate gyms. In certain embodiments, activation of VPAC2 by VIP modulates learning and memory. In certain embodiments, VPAC2 also modulates circadian oscillations in the suprachiasmatic nucleus.
In certain embodiments, a “VIPR2 CNV” refers to a copy number variant (CNV) of a region or segment of nucleic acid of chromosome 7. In certain embodiments, the chromosome 7 is a human chromosome 7. In certain embodiments, the CNV comprises nucleic acid that is located on human chromosome 7 at position 7q36.3. In certain embodiments, the VIPR2 CNV comprises a subtelomeric region or segment of chromosome 7.
In certain embodiments, the VIPR2 CNV comprises a duplication of a region or segment of nucleic acid of chromosome 7.
In certain embodiments, the VIPR2 CNV comprises a triplication of a region or segment of nucleic acid of chromosome 7.
In certain embodiments, the VIPR2 CNV comprises a region or segment of an exon or intron of a VIPR2 gene.
In certain embodiments, the VIPR2 CNV comprises a region or segment of exon 3 of a VIPR2 gene.
In certain embodiments, the VIPR2 CNV comprises a region or segment of exon 4 of a VIPR2 gene.
In certain embodiments, the VIPR2 CNV comprises a region or segment of nucleic acid that is between about 1 and about 100 kb; or between about 1 and about 95 kb; or between about 1 and about 90 kb; or between about 1 and about 85 kb; or between about 1 and about 80 kb; or between about 1 and about 75 kb; or between about 1 and about 70 kb; or between about 1 and about 65 kb; or between about 1 and about 60 kb; or between about 1 and about 55 kb; or between about 1 and about 50 kb; or between about 1 and about 45 kb; or between about 1 and about 40 kb; or between about 1 and about 35 kb; or between about 1 and about 30 kb; or between about 1 and about 25 kb; or between about 1 and about 20 kb; or between about 1 and about 15 kb; or between about 1 and about 10 kb; or between about 1 and about 5 kb from the transcriptional start site of a VIPR2 gene.
In certain embodiments, the VIPR2 CNV comprises a region or segment of nucleic acid that is within 89 kb of the transcriptional start site of a VIPR2 gene.
In certain embodiments, the VIPR2 CNV comprises a region or segment of nucleic acid that is located in region 158,448,321-158,810,016 of human chromosome 7, as described by the human genome build NCBI36/hg18 (produced by the International Human Genome Sequencing Consortium).
In certain embodiments, the VIPR2 CNV comprises a region or segment of nucleic acid that is located in region 158,731,401-158,810,016 of human chromosome 7, as described by the human genome build NCBI36/hg18.
In certain embodiments, the VIPR2 CNV comprises a region or segment of nucleic acid that is located in region 158,448,321-158,605,936 of human chromosome 7, as described by the human genome build NCBI36/hg18.
In certain embodiments, the VIPR2 CNV comprises a region or segment of nucleic acid that is located in region 152,249,238-158,820,241; or 158,448,322-158,651,373; or 158,703,311-158,810,016; or 158,400,168-158,782,022; or 158,699,364-158,820,241; or 158,322,402-158,652,547; or 158,322,880-158,650,683; or 158,307,494-158,650,321; or 158,265,451-158,652,547; or 158,171,999-158,826,119; or 158,589,435-158,820,241; or 158,618,998-158,780,576; or 158,275,536-158,779,706; or 158,699,364-158,820,241; or 157,571,900-158,677,051; or 158,096,725-158,502,254; or 158,121,784-158,588,486; or 158,276,605-158,819,766; or 158,384,773-158,512,459; or 158,569,154-158,819,536; or 158,664,993-158,819,536; or 158,690,919-158,819,536; or 158,690,919-158,819,536; or 158,690,919-158,819,536; or 158,712,866-158,819,536; or 158,712,866-158,819,536; or 158,712,866-158,819,536; or 158,719,441-158,819,536; or 158,719,441-158,819,536; or 158,543,969-158,614,169; or 158,540,185-158,621,330; or 158,396,852-158,456,333; or 158,562,524-158,614,169; or 158,550,023-158,614,169; or 158,654,446-158,726,336 of human chromosome 7, as described by the human genome build NCBI36/hg18.
According to the present disclosure, VIPR2 CNVs are present in individuals with schizophrenia, or at risk of developing schizophrenia or schizoaffective disorder. Detection of these VIPR2 CNVs can be used to diagnose schizophrenia or schizoaffective disorder in a subject and can be used in conjunction with other criteria, such as those set forth in the Diagnostic and Statistical Manual of Mental Disorders IV (DSMIV) to support a clinical psychiatric diagnosis of schizophrenia or schizoaffective disorder. Detection of a VIPR2 CNV is also helpful to identify any adverse health effects associated therewith, and thus detection can be useful for finding and treatment of schizophrenia. While the present disclosure is exemplified in humans, its extension to other species including mammals is contemplated. Assays such as RT-PCR, PCR, qPCR, DNA and RNA sequencing, microarray analysis and any other genome-based analyses known in the art, along with any immunoassays known in the art, may be used to detect a VIPR2 CNV in a sample. In addition, such analyses may be qualitative or quantitative.
According to the present disclosure, a “subject” or “patient” is a human or non-human animal. Although the animal subject is preferably a human, the concepts, compounds and compositions of the disclosure have application in veterinary medicine as well, e.g., for the treatment of domesticated species, farm animal species, and wild animals or zoological garden animals.
In humans, the VIPR2 CNVs disclosed herein may be detected individually or in combination to provide a diagnostic evaluation of schizophrenia. Other VIPR2 CNVs from other species may prove useful, alone or in combination, for similar purposes.
In certain embodiments, the present disclosure provides for methods of diagnosing schizophrenia in a subject comprising detecting an increased level of VIPR2 expression in a sample from the subject compared to the level of VIPR2 expression in a sample from a control subject that does not have schizophrenia. In certain embodiments, the level of expression of VIPR2 is the level of transcription of VIPR2 in the samples, for example, as determined by the level of VIPR2 mRNA. In certain embodiments, the level of expression of VIPR2 is the level of VIPR2 protein detected in the samples, for example, the level of Vasoactive Intestinal Peptide Receptor 2 (VPAC2) protein detected in the samples.
In certain embodiments, the present disclosure provides for methods of diagnosing schizophrenia in a subject by detecting an increased level of cyclic-AMP activation, signaling or accumulation in a sample from the subject when the sample is contacted with a VPAC2 agonist, for example, VIP or BAY55-9837, compared to the level of cyclic-AMP activation, signaling or accumulation in a sample from a control subject that is contacted with a VPAC2 agonist, wherein the control subject does not have schizophrenia. As one non-limiting example, the sample of the subject may comprise lymphoid or lymphoblastic cells from the subject, which are treated with VPAC2 agonist, and then the consequent change in cAMP level may be measured.
In certain embodiments, the samples described herein can be derived from any tissues, cells and/or cells in biological fluids from, for example, a mammal or human to be tested.
In certain embodiments, a VIPR2 CNV is associated with a second disorder or condition, for example, a psychiatric condition, disorder or phenotype. In certain embodiments, the VIPR2 CNV is associated with a pediatric neurodevelopment disorder. In certain embodiments, the VIPR2 CNV is associated with autism As used herein, the terms “associated with a second disorder” mean that the VIPR2 CNV and a second disorder (i.e., a disorder other than schizophrenia), exhibit a degree of linkage, wherein the VIPR2 CNV and the second disorder are present together in an individual at a higher frequency than if their occurrences were independent of each other.
The present disclosure encompasses nucleic acid segments that are complementary, or essentially complementary, to the nucleic acid regions or segments of chromosome 7 described herein, for example, the VIPR2 gene and its chromosomal region, for example 89 kb upstream of the transcriptional start site. Nucleic acid sequences that are “complementary” are those that are capable of base-pairing according to the standard Watson-Crick complementary rules. As used herein, the term “complementary sequences” means nucleic acid sequences that are substantially complementary, as defined as being capable of hybridizing to a specified nucleic acid segment, under relatively stringent conditions as known in the art. Such sequences may encode the entire VIPR2 protein product encompassed herein or functional or non-functional fragments thereof.
A nucleic acid may be contained in a host cell, in some cases, capable of expressing the product of that nucleic acid. In addition to diagnostic and therapeutic considerations, cells expressing nucleic acids of the present disclosure may prove useful in the context of screening for agents that induce, repress, inhibit, augment, interfere with, block, abrogate, stimulate or enhance the expression, distribution, turnover, or detectability of VIPR2 CNVs.
Hybridizing segments may be relatively short nucleic acids, often termed oligonucleotides. Sequences of at least 10 bases long, for example, sequences of at least 17 or at least 22 bases long, should occur only once in the human genome and, therefore, suffice to specify a unique target sequence. Although shorter oligomers are easier to make and increase in vivo accessibility, numerous other factors are involved in determining the specificity of hybridization. Both binding affinity and sequence specificity of an oligonucleotide to its complementary target increases with increasing length. It is contemplated that exemplary oligonucleotides of any number from 8 to 100 or more base pairs will be used, although others are contemplated. Longer polynucleotides are contemplated as well. Such oligonucleotides will find use, for example, as probes in Southern and Northern blots and as primers in amplification reactions.
Suitable hybridization conditions will be well known to those of skill in the art. Accordingly, the nucleotide sequences of the disclosure may be used for their ability to selectively form duplex molecules with complementary stretches of DNA or RNA fragments. Depending on the application envisioned, one will desire to employ varying conditions of hybridization to achieve varying degrees of selectivity of probe towards target sequence. For applications requiring high selectivity, one will typically desire to employ relatively stringent conditions to form the hybrids, e.g., one will select relatively low salt and/or high temperature conditions, such as provided by about 0.02 M to about 0.15 M NaCl at temperatures of about 50° C. to about 70° C. Such selective conditions tolerate little, if any, mismatch between the probe and the template or target strand, and would be particularly suitable for isolating protein-encoding DNA segments. Detection of DNA segments via hybridization is well-known to those of skill in the art, and the teachings of U.S. Pat. Nos. 4,965,188 and 5,176,995 (each specifically incorporated herein by reference) are exemplary of the methods of hybridization analyses.
In certain applications, lower stringency conditions are required. Under these conditions, hybridization may occur even though the sequences of probe and target strand are not perfectly complementary, with mismatches at one or more positions. Conditions may be rendered less stringent by increasing salt concentration and decreasing temperature. For example, a medium stringency condition could be provided by about 0.1 to 0.25 M NaCl at temperatures of about 37° C. to about 55° C., while a low stringency condition could be provided by about 0.15 M to about 0.9 M salt, at temperatures ranging from about 20° C. to about 55° C. Thus, hybridization conditions can be readily manipulated, and thus will generally be a method of choice depending on the desired results. Cross-hybridizing species can thereby be readily identified as positively hybridizing signals with respect to control hybridizations. In any case, it is generally appreciated that conditions can be rendered more stringent by the addition of increasing amounts of formamide, which serves to destabilize the hybrid duplex in the same manner as increased temperature.
In certain embodiments, it will be advantageous to employ nucleic acid sequences of the present disclosure in combination with an appropriate means, such as a label, for determining hybridization. A wide variety of appropriate indicator means are known in the art, including fluorescent, radioactive, enzymatic or other ligands, such as avidin/biotin, which are capable of giving a detectable signal. In certain embodiments, one will likely desire to employ a fluorescent label or an enzyme tag, such as urease, alkaline phosphatase, luciferase, or peroxidase, instead of radioactive or other environmental undesirable reagents. In the case of enzyme tags, colorimetric indicator substrates are known that can be employed to provide a means visible to the human eye or spectrophotometrically, to identify specific hybridization with complementary nucleic acid-containing samples.
In general, it is envisioned that the hybridization probes described herein will be useful both as reagents in solution hybridization as well as in embodiments employing a solid phase. In embodiments involving a solid phase, the test DNA (or RNA) is adsorbed or otherwise affixed to a selected matrix or surface. This fixed, single-stranded nucleic acid is then subjected to specific hybridization with selected probes under desired conditions. The selected conditions will depend on the particular circumstances based on the particular criteria required (depending, for example, on the G+C content, type of target nucleic acid, source of nucleic acid, size of hybridization probe, etc.). Following washing of the hybridized surface so as to remove nonspecifically bound probe molecules, specific hybridization is detected, or even quantitated, by means of the label.
Probes and primers of the present disclosure are useful for PCR, qPCR, nucleic acid sequencing, microarray analysis, site-directed, and site-specific mutagenesis. Site-specific mutagenesis is a technique useful in the preparation of individual peptides, or biologically functional equivalent proteins or peptides, through specific mutagenesis of the underlying DNA. The technique further provides a ready ability to prepare and test sequence variants, incorporating one or more of the foregoing considerations, by introducing one or more nucleotide sequence changes into the DNA. Site-specific mutagenesis allows the production of mutants through the use of specific oligonucleotide sequences which encode the DNA sequence of the desired mutation, as well as a sufficient number of adjacent nucleotides, to provide a primer sequence of sufficient size and sequence complexity to form a stable duplex on both sides of the deletion junction being traversed. Typically, a primer of about 17 to 25 nucleotides in length is preferred, with about 5 to 10 residues on both sides of the junction of the sequence being altered.
In certain embodiments of the disclosure, the presence of a VIPR2 CNV or VIPR2 expression product in a sample can be determined by using nucleic acid microarrays, or gene chip technology (see, e.g., U.S. Pat. No. 7,455,975). As used herein, a “microarray” is an array of distinct polynucleotides, oligonucleotides, polypeptides, peptides, or antibodies affixed to a substrate, such as paper, nylon, or other type of membrane; filter; chip; glass slide; silicone or any other type of suitable support. The use of microarray technology for detecting VIPR2 CNV is further described in the embodiments disclosed in the Examples of the present application.
In certain embodiments, the microarray technology involves the positioning of highly condensed and ordered arrays of nucleic acid probes, for example, DNA oligonucleotides, on a substrate, for example, a glass slide or nylon membrane. Each oligonucleotide may comprise a nucleotide sequence that is complementary to a portion of, for example, a VIPR2 CNV or VIPR2 expression product, wherein the oligonucleotide can be placed, for example, on a single glass slide or nylon membrane. The resulting microarrays can then be used to screen for the presence of a VIPR2 CNV or VIPR2 expression product expressed in a sample to be screened.
In certain embodiments, a nucleic acid microarray may be utilized by preparing labeled nucleic acid from a sample to be screened, and hybridizing such labeled nucleic acid with the array. In addition, labeled nucleic acid of a designated control sequences may be prepared (or in the event that the array is sold as part of a kit, could be supplied to the user). Radioactive, colorimetric, chemiluminescent or fluorescent tags may be used for labeling of nucleic acid sequences from the sample and for the control. Numerous techniques for scanning arrays, detecting fluorescent, chemiluminescent, or colorimetric output, are known in the art and may be used for detecting hybridization of a nucleic acid from a test sample to the microarray. For example, a high-throughput fluorescent microarray scanning system (ScanArray®, PerkinElmer Life And Analytical Sciences, Inc., Waltham, Mass., USA), or a colorimetric microarray scanner (ArrayIt® SpotWare™, TeleChem International, Inc., Sunnyvale, Calif., USA) can be used. Additional microarray systems that can be used according to the methods of the present disclosure include the NimbleGen platform (Roche NimbleGen, Inc., Madison, Wis.), GeneChip® Human Mapping Array Sets (Affymetrix, Inc., Santa Clara, Calif.) and Genome-Wide Human SNP Arrays (Affymetrix, Inc., Santa Clara, Calif.).
The present disclosure contemplates the preparation of one or more specialized microarrays (e.g., oligonucleotide microarrays or cDNA microarrays) comprising one or more polynucleotides encoding one or more VIPR2 CNV, VIPR2 nucleic acid sequence or complementary sequences, or fragments thereof. In accordance with this aspect of the disclosure, the oligonucleotide sequences or cDNA sequences include any of the disclosed VIPR2 CNV or VIPR2 polynucleotides or fragments or combinations thereof, and are contained on a microarray, e.g., a oligonucleotide microarray or cDNA microarray in association with, or introduced onto, any supporting materials, such as glass slides, nylon membrane filters, glass or polymer beads, or other types of suitable substrate material.
Methods for producing and using DNA microarrays are well known in the art (see, e.g. Rampal J. B., 2001, DNA Arrays: Methods and Protocols, Humana Press, Inc., Totowa, N.J.; Schena M., 2000, Microarray Biochip Technology, Eaton Publishing; and Schena M., 2002, Microarray Analysis, John Wiley & Sons). For example, to determine gene expression using microarray technology, polynucleotides, e.g., RNA, DNA, or cDNA, are isolated from a biological sample, e.g., cells expressing a VIPR2. The isolated nucleic acid is detectably labeled, e.g., by fluorescent, enzyme, or chemiluminescent label, and applied to a microarray, e.g., one or more nucleic acid microarrays provided by this disclosure which comprises, for example, oligonucleotides complimentary to the labeled cellular derived nucleic acid applied to the microarray. The array is then washed to remove unbound material and visualized by staining or fluorescence, or other means known in the art depending on the type of label utilized.
For the purpose of example, and not limitation, microarrays of the disclosure may be prepared by amplifying VIPR2 CNVs and/or VIPR2 cDNAs or fragments thereof by PCR or RT-PCR, and arraying the PCR products from a microtiter plate onto silyated microscope slides using high-speed robotics. Printed arrays may be incubated in a humid chamber to allow rehydration of the array elements and rinsed, in, for example, 0.2% SDS, water, and sodium borohydride solutions.
Further, a commercially available microarray and/or commercially available software may be used to detect VIPR2 CNV according to the invention (see, for instance, the working example below). Non-limiting examples of microarrays which may be used to detect a VIPR2 CNV include microarrays and associated software marketed by Affymetrix, such as Genome-wide Human SNP Array 6 and the associated Genotyping Console (e.g., version 4.1.1, prior or subsequent versions), or Illumina, such as Whole Genome Genotyping and Copy Number Variation Analysis, such as the Omni Family of Microarrays, and associated software.
In certain embodiments, the present disclosure provides methods for the detection of a VIPR2 CNV or VIPR2 expression product comprising the use of Fluorescence in situ Hybridization (FISH). The term “in situ hybridization” generally refers to hybridization of a nucleic acid probe to a nucleic acid target that is part of a cytological or histological preparation. Typically, FISH methods involve the following steps: (a) fixing the tissue or other biological material under investigation to a support (e.g., glass slide or wall of a micro titer well), (b) treatment of the tissue or material to increase accessibility of FISH probe to target nucleic acid, (c) contacting the tissue or material containing the target nucleic acid with probes to form specific hybridization complexes, (d) post hybridization washes of the complexes to selectively remove probes that are not specifically hybridized to the target, and (e) detection of probes that have formed hybridization complexes with target nucleic acid molecules. Such methods are described in a number of sources, including: Gall and Pardue, (1981) Methods of Enzymology 21:470-480; Henderson, (1982) International Review of Cytology, 76:1-46; and Angerer, et al., (1985) in Genetic Engineering: Principles and Methods (Setlow and Hollaender, Eds.) vol. 7, pp. 43-65, Plenum Press, New York.
In certain embodiments of the present disclosure, the presence of VIPR2 CNVs can be detected by polymerase chain reaction (PCR).
In certain embodiments, the VIPR2 CNV is amplified from a genomic DNA sample of a subject. Any one or more of the VIPR2 CNVs disclosed herein can be amplified through PCR by using at least one set of primers for each VIPR2 CNV. Following PCR amplification, the PCR amplification products can be sequenced using standard techniques known in the art, and the sequence can be compared to the sequence of a control sample, or to the corresponding wild type sequence.
In certain embodiments, the presence of the VIPR2 CNV can be determined by detecting a difference in size between the subject sample PCR product and the control sample PCR product. Such differences in size can be determined, for example, by gel electrophoresis or any other method known in the art for detecting the size of a PCR product.
According to the present disclosure, the expression level of a VIPR2 gene can be detected through the use of quantitative polymerase chain reaction (qPCR), or quantitative reverse transcriptase-polymerase chain reaction (qRT-PCR), which utilizes competitive techniques employing an internal homologous control that differs in size from the target, for example, by a small insertion or deletion. Non-competitive and kinetic quantitative PCR or RT-PCR may also be used. Experiments may combine real-time, kinetic PCR or RT-PCR detection together with an internal homologous control that can be simultaneously detected alongside the target sequences. In certain embodiments, real time quantitative PCR may provide the capability of measuring the level of VIPR2 gene product amplified through PCR. In certain embodiments, quantitative PCR may require only a nominal amount of a sample to perform such experiments.
Quantitative amplification is based on the monitoring of a signal (e.g., fluorescence of a probe) representing copies of a template in cycles of an amplification (e.g., PCR) reaction. In the initial cycles of the PCR, a very low signal is observed because the quantity of the amplification product formed does not support a measurable signal output from the assay. After the initial cycles, as the amount of formed amplification product increases, the signal intensity increases to a measurable level and reaches a plateau in later cycles when the PCR enters into a non-logarithmic phase. Through a plot of the signal intensity versus the cycle number, the specific cycle at which a measurable signal is obtained from the PCR reaction can be deduced and used to back-calculate the quantity of the target before the start of the PCR. The number of the specific cycles that is determined by this method is typically referred to as the cycle threshold (Ct). Exemplary methods are described in, e.g., U.S. Pat. Nos. 6,180,349; 6,033,854; and 5,972,602.
In certain embodiments, a method for detection of amplification products is, for example, the 5′-3′ exonuclease activity during PCR reaction (also referred to as the TaqMan™ assay) (see, e.g., U.S. Pat. Nos. 5,210,015 and 5,487,972). This assay detects the accumulation of a specific PCR product by hybridization and cleavage of a doubly labeled fluorogenic probe (the “TaqMan™” probe) during the amplification reaction. The fluorogenic probe consists of an oligonucleotide labeled with both a fluorescent reporter dye and a quencher dye. During PCR, this probe is cleaved by the 5′-exonuclease activity of DNA polymerase if it hybridizes to the segment being amplified. Cleavage of the probe generates an increase in the fluorescence intensity of the reporter dye.
In certain embodiments, detection of amplification products may utilize, by way of example, and not by way of limitation, energy transfer according to the “beacon probe” method described by Tyagi and Kramer (Nature Biotech. 1996; 14:303-309.), and disclosed by U.S. Pat. Nos. 5,119,801 and 5,312,728. This method employs oligonucleotide hybridization probes that can form hairpin structures. On one end of the hybridization probe (either the 5′ or 3′ end), there is a donor fluorophore, and on the other end, an acceptor moiety. In the Tyagi and Kramer method, the acceptor moiety is a quencher, wherein the acceptor absorbs energy released by the donor, but then does not itself fluoresce. Thus, when the beacon is in the open conformation, the fluorescence of the donor fluorophore is detectable, whereas when the beacon is in hairpin (closed) conformation, the fluorescence of the donor fluorophore is quenched. When employed in PCR, the molecular beacon probe, which hybridizes to one of the strands of the PCR product, is in the open conformation and fluorescence is detected, while those that remain unhybridized will not fluoresce (Tyagi and Kramer, 1996). As a result, the amount of fluorescence will increase as the amount of PCR product increases, and thus may be used as a measure of the progress of the PCR. Those of skill in the art will recognize that other methods of quantitative amplification are also available.
Various other techniques for performing quantitative amplification of nucleic acids are also known. For example, some methodologies employ one or more probe oligonucleotides that are structured such that a change in fluorescence is generated when the oligonucleotide(s) is hybridized to a target nucleic acid. For example, one such method involves is a dual fluorophore approach that exploits fluorescence resonance energy transfer (FRET), e.g., LightCycler™ hybridization probes, where two oligo probes anneal to the amplification product. The oligonucleotides are designed to hybridize in a head-to-tail orientation with the fluorophores separated at a distance that is compatible with efficient energy transfer. Other examples of labeled oligonucleotides that are structured to emit a signal when bound to a nucleic acid or incorporated into an extension product include: Scorpions™ probes (e.g., Whitcombe et al., Nature Biotechnology 1999; 17:804-807, and U.S. Pat. No. 6,326,145), Sunrise™ (or Amplifluor™) probes (e.g., Nazarenko et al., Nuc. Acids Res. 1997; 25:2516-2521, and U.S. Pat. No. 6,117,635), and probes that form a secondary structure that results in reduced signal without a quencher and that emits increased signal when hybridized to a target (e.g., Lux Probes™).
Quantitation of a specific amplified product at the end of an amplification reaction (i.e., end-point PCR) can be employed to quantify the sequences in the final amplified population that match the sequence of DNA which remained undigested by a restriction enzyme. The end-point PCR analysis may be employed under conditions in which the reaction can be analyzed before the reactant nears depletion for a quantitative comparison. Most typically this is done through a comparison of reaction products following a limited number of cycles. For example, a reaction is allowed to cycle 10 times, 15 times, 20 times or 30 times. The quantities of end point PCR products can be compared to each other and an analysis of sequences from the differential enzyme treatments of the DNA sample can be made.
In certain embodiments, the presence of a VIPR2 CNV in a subject or control sample can be determined through sequencing (i.e. determining the nucleotide order of a given DNA or RNA fragment) of a genomic DNA product present in the sample. Any sequencing methods known in the art may be used to determine the nucleotide order of the VIPR2 CNV DNA or RNA.
For example, and not by way of limitation, chain terminator sequencing (i.e. Sanger sequencing) may be used to sequence the VIPR2 CNV, wherein extension of a polynucleotide is initiated at a specific site on the template VIPR2 CNV nucleic acid (e.g., DNA) by using a short oligonucleotide “primer” complementary to the template at that region. The classical chain-termination method requires a single-stranded DNA template, a DNA primer, a DNA polymerase, radioactively or fluorescently labeled nucleotides, and modified nucleotides that terminate DNA strand elongation (e.g., dideoxynucleotides). The DNA sample may be divided into four separate sequencing reactions, containing all four of the standard deoxynucleotides (dATP, dGTP, dCTP and dTTP) and the DNA polymerase. One of the four dideoxynucleotides (ddATP, ddGTP, ddCTP, or ddTTP) are added to each of the four reactions, which are the chain-terminating nucleotides, lacking a 3′-OH group required for the formation of a phosphodiester bond between two nucleotides, thus terminating DNA strand extension and resulting in various DNA fragments of varying length.
Newly synthesized and labeled DNA fragments are heat denatured, and separated by size by, for example, gel electrophoresis, with each of the four reactions run in one of four individual lanes of the gel (lanes A, T, G, C). The DNA bands may be visualized by autoradiography or UV light, and the DNA sequence can be directly read off the X-ray film or gel image.
In one embodiment, the primer is labeled (e.g., a fluorescent or radioactive label). In other embodiments, the chain-terminator nucleotides are labeled, for example, in ‘dye terminator sequencing’. In dye terminator sequencing, complete sequencing may be performed in a single reaction, wherein each of the di-deoxynucleotide chain-terminators (e.g., ddATP, ddGTP, ddCTP, and ddTTP) are labeled with a separate fluorescent dye which fluoresces at a different wavelength. The sequence of the template may be determined by separating the synthesized polynucleotide by size and determining the order of the dye signals exhibited by the reaction products.
In certain embodiments, sequencing may be performed according to the “pyrosequencing” method as described in, for example, Ronaghi et al., Analytical Biochemistry 1996; 242(1):84-9; Ronaghi et al, Science 1998; 281:363-365; and Nyrén, Methods Mol Biology 2007; 373: 1-14. Pyrosequencing is a nucleic acid (e.g., DNA) sequencing technique that relies on detection of pyrophosphate release upon nucleotide incorporation rather than chain termination with dideoxynucleotides. Thus, detection of the nucleotide order of the polynucleotide synthesized in the synthesis reaction may be determined in real time as the polynucleotide is extended.
The sequencing of a nucleic acid sample (i.e. determining the nucleotide order of a given DNA or RNA fragment) is not limited to any one technique. The present disclosure contemplates the use of any sequencing technique known in the art and, for example, new sequencing techniques arising in the future of the sequencing art.
Immunologic Detection of VIPR2 Protein and c-AMP
In certain embodiments, the present disclosure entails the use of antibodies in the immunologic detection of VIPR2 protein or c-AMP present in a subject or control sample. According to the present disclosure, such immunological detection can be used to detect and quantify the amount of VIPR2 protein or c-AMP present in the subject or control samples. Various useful immunodetection methods have been described in the scientific literature, such as, e.g., Nakamura et al. Handbook of Experimental Immunology (4th Ed.), Weir, E., Herzenberg, L. A., Blackwell, C., Herzenberg, L. (Eds.), Vol. 1, Chapter 27, Blackwell Scientific Publ., Oxford, 1987. Immunoassays, in their most simple and direct sense, are binding assays. Certain immunoassays include, but are not limited to, enzyme linked immunosorbent assays (ELISAs), Western blots and radioimmunoassays (RIA). Immunohistochemical detection using tissue sections also is particularly useful. However, it will be readily appreciated that detection is not limited to such techniques. For example, Western blotting, dot blotting, FACS analyses, and the like also may be used in connection with the present disclosure. Other assays include immunoprecipitation of labeled ligands and immunocytochemistry, both in vitro and in vivo.
In certain embodiments, the immunological methods of the disclosure may detect the total VIPR2 protein or c-AMP level present in a subject and control samples.
In general, immunobinding methods include obtaining a sample suspected of containing a protein, peptide or antigen, and contacting the sample with an antibody or protein or peptide in accordance with the present disclosure, as the case may be, under conditions effective to allow the formation of immunocomplexes. Preferred samples, according to the present disclosure, include, but are not limited to, fluids, such as plasma, serum, cerebrospinal fluid, sputum, saliva, breast milk, tears, bile, semen, vaginal secretion, amniotic fluid, urine or stool sample, as well extracts of cells such as leukocytes, bone marrow cells, buccal cells, fibroblasts and tissue biopsies.
Contacting a biological sample with the protein, peptide or antibody under conditions effective and for a period of time sufficient to allow the formation of immune complexes (primary immune complexes) generally comprises adding the composition, for example an antibody, to the sample and incubating the mixture for a period of time long enough for the antibodies to form immune complexes with the VIPR2 protein or c-AMP. After this time, the VIPR2 protein- or c-AMP-antibody mixture will be washed to remove any non-specifically bound antibody species, allowing only those antibodies specifically bound within the primary immune complexes to be detected.
In general, the detection of immunocomplex formation is well known in the art and may be achieved through the application of numerous approaches. These methods are generally based upon the detection of a label or marker, such as any radioactive, fluorescent, biological or enzymatic tags or labels of standard use in the art. U.S. patents concerning the use of such labels include U.S. Pat. Nos. 3,817,837; 3,850,752; 3,939,350; 3,996,345; 4,277,437; 4,275,149 and 4,366,241, each of which are incorporated herein by reference. A secondary binding ligand, such as a second antibody or a biotin/avidin ligand binding arrangement, as is known in the art, may also be used to detect the antibody-VIPR2 protein or antibody-c-AMP immunocomplex.
In certain embodiments, the primary immune complexes may be detected by means of a second binding ligand that has binding affinity for the VIPR2 protein or c-AMP or for the VIPR2 protein- or c-AMP-specific first antibody. In these cases, the second binding ligand may be linked to a detectable label. The second binding ligand is itself often an antibody, which may thus be termed a “secondary” antibody. The primary immune complexes are contacted with the labeled, secondary binding ligand, or antibody, under conditions effective and for a period of time sufficient to allow the formation of secondary immune complexes. The secondary immune complexes are then generally washed to remove any non-specifically bound labeled secondary antibodies or ligands, and the remaining label in the secondary immune complexes is then detected.
Further methods include the detection of primary immune complexes by a two step approach. A second binding ligand, such as an antibody that has binding affinity for the VIPR2 protein or c-AMP is used to form secondary immune complexes, as described above. The second binding ligand contains an enzyme capable of processing a substrate to a detectable product and, hence, amplifying signal over time. After washing, the secondary immune complexes are contacted with substrate, permitting detection.
As a part of the practice of the present disclosure, the principles of an enzyme-linked immunosorbent assay (ELISA) may be used (Engvall and Perlmann, Immunochem., 1971; 8:871-873; Engvall, Methods Enzymol., 1980; 70 (A):419-39; Engvall, Lancet, 1976; 2(8000):1410; Engvall, Med. Biol., 1977; 55(4): 193-200; Gripenberg et al., Scand J. Immunol., 1978; 7(2):151-7). Such a method can be used for the detection and quantification of VIPR2 protein or c-AMP. ELISA allows for substances to be passively adsorbed to solid supports such as plastic to enable facile handling under laboratory conditions. For a comprehensive treatise on ELISA the skilled artisan is referred to “ELISA; Theory and Practice” (Crowther, Methods in Molecular Biology, 1995; 42:1-218, incorporated herein by reference). In general, in an ELISA, an antigen (or antibody specific for an antigen) is affixed to a surface, and then a specific antibody (or antigen) is washed over the surface so that it can bind to the antigen (or antibody). The antigen-antibody complex may then be detected, for example, by the conversion of a substrate to a detectable signal by an enzyme linked to the antibody, or through the use of labeled secondary or tertiary antibodies specific for the antigen-antibody complex.
The sensitivity of ELISA methods is dependent on the turnover of the enzyme used and the ease of detection of the product of the enzyme reaction. Enhancement of the sensitivity of these assay systems can be achieved by the use of fluorescent and radioactive substrates for the enzymes.
In certain embodiments, the disclosure comprises a “sandwich” ELISA, where anti-VIPR2 protein or ant-c-AMP antibodies are immobilized onto a selected surface, such as a well in a polystyrene microtiter plate or a dipstick. Then, a test composition suspected of containing VIPR2 protein or c-AMP, e.g., a clinical sample, is contacted with the surface. After binding and washing to remove non-specifically bound immunocomplexes, the bound antigen may be detected by a second antibody to the anti-VIPR2 protein or ant-c-AMP antibodies.
In certain embodiments, polypeptides from the sample are immobilized onto a surface and then contacted with the anti-VIPR2 protein or ant-c-AMP antibodies. After binding and washing to remove non-specifically bound immune complexes, the bound antibody is detected. Where the initial antibodies are linked to a detectable label, the primary immune complexes may be detected directly. Alternatively, the immune complexes may be detected using a second antibody that has binding affinity for the first antibody, with the second antibody being linked to a detectable label.
In certain embodiments, an ELISA in which the VIPR2 protein or c-AMP are immobilized utilizes antibody competition for detection. In this ELISA, labeled antibodies are added to the wells, allowed to bind to the VIPR2 protein or c-AMP, and detected by means of their label. The amount of VIPR2 protein or c-AMP in a sample is determined by mixing the sample with the labeled antibodies before or during incubation with coated wells. The presence of VIPR2 protein or c-AMP in the sample acts to reduce the amount of antibody available for binding to the well, and thus reduces the ultimate signal.
In another embodiment of the present disclosure, a subject that has a VIPR2 CNV, an increased level of VIPR2 expression compared to a non-schizophrenic control, or an increased level of c-AMP compared to a non-schizophrenic control is a candidate for VIPR2 therapy, wherein an agent is administered in an amount effective to decrease the level of VIPR2 in a sample from the subject.
In certain embodiments the agent inhibits the function of VIPR2 protein or reduces the level of functional VIPR2 protein.
The agent can be administered, for example, systemically (e.g. by intravenous injection, oral administration, inhalation, etc.), by intra-arterial, intramuscular, intradermal, transdermal, subcutaneous, oral, intraperitoneal, intraventricular, or intrathecal administration, or may be administered by any other means known in the art.
In certain embodiments, the agent is a VIPR2 protein antagonist or inhibitor, for example, as described in Moreno, et al., 2000, Peptides 21:1543-1549; and Chu et al., 2009, Mol. Pharmacol. 77:95-101.
In certain embodiments, the VIPR2 protein antagonist or inhibitor is administered in an amount effective to reduce or inhibit the ability of VIPR2 protein to bind to VIP.
In certain embodiments, the VIPR2 protein antagonist or inhibitor is administered in an amount effective to reduce or inhibit the ability of VIPR2 protein to activate cyclic-AMP signaling, for example, cyclic-AMP accumulation, or protein kinase A (PKA) activation.
In certain embodiments, the VIPR2 protein antagonist or inhibitor is administered in an amount effective to reduce or inhibit the ability of VIPR2 protein to regulate synaptic transmission in the hippocampus.
In certain embodiments, the VIPR2 protein antagonist or inhibitor is administered in an amount effective to reduce or inhibit the ability of VIPR2 protein to promote proliferation of neural progenitor cells, for example, in the dentate gyrus.
In certain embodiments, the VIPR2 protein antagonist or inhibitor is administered in an amount effective to reduce or inhibit the ability of VIPR2 protein to modulate circadian oscillations in, for example, the suprachiasmatic nucleus.
In certain embodiments, the VIPR2 protein antagonist or inhibitor is an antibody that binds to the VIPR2 protein. An antibody can be a polyclonal or a monoclonal antibody composition. Such antibodies may also include but are not limited to chimeric, human, humanized, single chain, Fab fragments, and a Fab expression library. Means for preparing and characterizing antibodies are well known in the art and can be readily prepared through use of well-known techniques, such as those exemplified in U.S. Pat. No. 4,196,265, incorporated herein by reference (see also, e.g., Kohler and Milstein, Nature, 1975; 256:495-497; and Kohler and Milstein, Eur. J. Immunol., 1976; 6:511-519; Goding, Monoclonal Antibodies: Principles and Practice, 2d ed., Academic Press, Orlando, Fla., pp 60-61, 71-74, 1986; and Harlow and Lane, Antibodies, a Laboratory Manual, Cold Spring Harbor Laboratory, pp 139-281, 1988). Specifically, techniques developed for the production of chimeric and humanized antibodies have been described by Neuberger, et al., Nature, 1984; 312:604-608; Brüggemann et al., Proc Natl Acad Sci USA, 1989; 86(17): 6709-6713; and Takeda et al., Nature, 1985; 314:452-454. Various techniques have been described for the production of single chain antibodies including those in U.S. Pat. Nos. 5,476,786; 5,132,405; 4,946,778 each of which is incorporated herein by reference, and these techniques can be adapted to produce for example, VIPR2 protein-specific single chain antibodies. An additional embodiment of the disclosure may utilize the techniques described for the construction of Fab expression libraries of Huse et al., Science, 1989; 246:1275-1281, to allow rapid and easy identification of monoclonal Fab fragments with the desired specificity for a VIPR2 protein.
In certain embodiments, the agent is an antisense molecule, RNAi molecule or siRNA molecule. In certain embodiments, the antisense, RNAi or siRNA molecule is complementary to a segment or region of a VIPR2 mRNA transcript. In certain embodiments, the antisense, RNAi or siRNA molecule hybridizes to and inhibits or reduces translation of VIPR2 mRNA. In certain embodiments, the antisense, RNAi or siRNA molecule hybridizes to VIPR2 mRNA and increases degradation of the VIPR2 mRNA.
In further embodiments, the present application provides kits, such as an immunological kit, for use in detecting a VIPR2 CNV, VIPR2 nucleic acid or protein, and/or c-AMP in a biological sample. Such kits will generally comprise one or more oligonucleotides and/or antibodies that have specificity for VIPR2 CNV, VIPR2 nucleic acid or protein, and/or c-AMP.
In certain embodiments, a kit for detection of a VIPR2 CNV or VIPR2 nucleic acid will comprise, in suitable container means, one or more control VIPR2 CNV or control VIPR2 samples, and one or more oligonucleotide that specifically hybridizes to a DNA region associated with a VIPR2 CNV or VIPR2, as set forth above, for use in PCR, RT-PCR, qPCR, qRT-PCR, microarray analysis or nucleic acid sequencing. The kit may also comprise one or more polymerase, reverse transcriptase, and nucleotide bases, wherein the nucleotide bases may be further detectably labeled.
In certain embodiments, the oligonucleotide primers are immobilized on a solid surface or support, for example, on a nucleic acid microarray, wherein the position of each oligonucleotide primer bound to the solid surface or support is known and identifiable.
In certain embodiments, the immunodetection kits will comprise, in suitable container means, one or more control VIPR2 protein or control c-AMP sample, and one or more antibodies that bind to VIPR2 protein or c-AMP, and antibodies that bind to other antibodies via Fc portions.
In certain embodiments, the VIPR2 protein or c-AMP or primary anti-VIPR2 protein antibody or anti-c-AMP antibody may be provided bound to a solid support, such as a column matrix or well of a microtitre plate. Alternatively, the support may be provided as a separate element of the kit.
The immunodetection reagents of the kit may include detectable labels that are associated with, or linked to, the given antibody or antigen itself. Detectable labels that are associated with or attached to a secondary binding ligand are also contemplated. Such detectable labels include, for example, chemiluminescent or fluorescent molecules (e.g., rhodamine, fluorescein, green fluorescent protein, luciferase, Cy3, Cy5, or ROX), radiolabels (e.g., 3H, 35S, 32P, 14C, 131I)) or enzymes (e.g., alkaline phosphatase, horseradish peroxidase).
The kits may further comprise suitable standards of predetermined amounts, including both oligonucleotides, antibodies and VIPR2 protein or c-AMP. These may be used to prepare a standard curve for a detection assay.
The kits of the disclosure, regardless of type, will generally comprise one or more containers into which the biological agents are placed and, preferably, suitably aliquoted. The components of the kits can be packaged either in aqueous media or in lyophilized form.
The container means of the kits will generally include at least one vial, test tube, flask, bottle, or even syringe or other container means, into which the antibody or antigen may be placed, and preferably, suitably aliquoted. Where a second or third binding ligand or additional component is provided, the kit will also generally contain a second, third or other additional container into which this ligand or component may be placed.
The kits of the present disclosure will also typically include a means for containing the nucleic acids, VIPR2 CNV control samples, VIPR2 protein samples, c-AMP samples, or antibodies and any other reagent containers in close confinement for commercial sale. Such containers may include injection or blow-molded plastic containers into which the desired vials are retained.
To identify novel schizophrenia genes, copy number variation genome-wide was investigated using an approach that detects enrichment of multiple overlapping rare variants. An unbiased, genome-wide rare CNV analysis was performed using a staged approach. Regions of interest were defined in a primary sample of 802 patients and 742 controls, and statistical significance was estimated in a larger sample of 7,488 patients and 6,689 controls. The primary cohort consisted of unrelated patients derived from family-based studies conducted by investigators at the University of Washington, McLean Hospital, Columbia University, Trinity College Dublin, NYU, and Harvard Medical School (see Table 2); all samples were analyzed by array comparative genomic hybridization (CGH) using the NimbleGen HD2 platform. The secondary cohort consisted of Affymetrix SNP Array 6.0 data from the MGS study of schizophrenia34, publicly available data from the International Schizophrenia Consortium (ISC), genotyped using Affymetrix 6.0 and 5.0 platforms23, and Affymetrix 6.0 data on an independent set of controls from the Bipolar Genome Study (BiGS)35, (see Table 2).
The Primary Cohort.
The initial discovery data set was composed of 1,761 unrelated subjects analyzed on the NimbleGen HD2 Array-CGH platform. The unfiltered sample consisted of 913 patients and 848 controls ascertained at ten sites (see Table 2). Ascertainment of these samples for family-based studies is described in previous publications1-8. Microarray hybridizations using the NimbleGen HD2 platform were performed at the service laboratory of Roche NimbleGen according to the manufacturer's specifications. Samples were filtered from the dataset based on data quality measures described in Sections 6.4 and 6.5. Five duplicate samples from Trinity College Dublin also present in the ISC dataset were removed from the primary dataset. The final discovery data set consisted of 906 unrelated patients and 742 controls (see Table 2).
The Secondary Cohort.
The secondary cohort consisted of Affymetrix Genome-Wide Human SNP Array 6.0 data from the Genetic Association Information Network (GAIN) supported case-control studies “Genome-Wide Association Study of Schizophrenia” (phs000021.v2.p1) and the “Molecular Genetics of Schizophrenia—nonGAIN Sample” (phs000167.v1.p1). Both of these studies are subsets of the larger “Molecular Genetics of Schizophrenia Genome-wide Association” study, so the datasets were combined and the merged cohort is referred to as “MGS”. In total, the combined MGS sample consisted of 4,195 unrelated cases and 3,804 controls. A further 429 controls assayed on the Affymetrix 6.0 by the Bipolar Genome Study (BiGS) were joined with the MGS cohort to increase the number of controls. Genotyping of the MGS and BiGS cohorts was performed at the Broad Institute Center for Genotyping and Analysis (http://www.broad.mit.node/306).
Several NIMH controls overlapped between primary and secondary cohorts. Therefore, se eliminated overlapping samples from the dataset by retaining the experiment with the greatest probe density. Sample and experiment redundancy were minimized in the secondary cohort by prioritizing studies in order of GAIN-MGS>nonGAIN-MGS>BiGS. In total, 4,097 cases and 3,508 controls assayed on the Affymetrix 6.0 platform were included in stage two of the analysis. CNV data from the International Schizophrenia Consortium (ISC)8 was further combined with the processed data to increase the final size of the secondary cohort to 7,488 cases and 6,689 controls. Because sample identifiers were not available for CNVs in the ISC dataset, each CNV in the ISC call set was treated as coming from an independent sample. This approximation is reasonable given that the mean rare CNV burden in this dataset was less than one CNV per sample8.
With the exception of the published CNV calls from the ISC, all data were processed and analyzed centrally as follows. Microarray intensity data were normalized, and CNV calls were generated using an analysis package developed for the present study called C-score. All CNV call sets were filtered in a consistent fashion. In order to minimize the differential sensitivity of the various array platforms to detect CNVs, the analysis was limited to CNVs≧100 Kb. This size range is readily detectable by all platforms used in this study. The same criteria have been previously applied to filter CNVs across studies23. Last, sensitivity to detect large (>100 Kb) copy number polymorphisms (CNPs) was evaluated at several locations in the genome, as described herein. Overall sensitivity to detect CNVs was comparable in cases and controls in both cohorts.
With the exception of data collected by ISC, all processing of microarray data was performed at Cold Spring Harbor Laboratory. Different methods were used for dual color intensity data from NimbleGen HD2 and single color intensity data from the Affymetrix Genome-Wide Human SNP 6.0 genotyping array.
NimbleGen HD2.
NimbleGen HD2 dual color intensity data were normalized in a two step process: (1) a spatial normalization of probes was performed to adjust for regional differences in intensities across the surface of the array, and (2) the Cy5 and Cy3 intensities were adjusted to a fitting curve by invariant set normalization.
Spatial normalization was performed using an R module provided by Roche-NimbleGen Inc. Invariant set normalization of intensity data involves selection of a probe set with minimal variability between the ranked test and reference autosomal probe intensities as described in Li et al10. The test intensities of this invariant autosomal probe set are then adjusted to the reference distribution. Based on these adjustments, a fitting curve is established to which all other intensities are shifted, preserving the variability in the data. The intensities of X and Y chromosomes were then extrapolated to the fitting curve. The process is repeated while exchanging the test and reference to simulate a dye swap experiment. The log 2 ratio for probe i is then estimated using the geometric mean of normalized and raw intensity data of test (Tst) and reference (Rfn) as follows:
Affymetrix Genome-Wide Human SNP Array 6.0.
A two-step normalization method was developed to process Affymetrix SNP Array single color intensity data (Affymetrix 500K, Affymetrix 5.0 & Affymetrix 6.0). First, all arrays are normalized by invariant set normalization to a single reference array. To obtain the single reference array, the median-most array was selected based on autosomal intensities; then the correlation matrix was built for the adjacent 100 experiments (50 lower and 50 higher than the median-most). The single array most highly correlated with all others was selected as the common reference for Invariant Set Normalization based on the sum of Pearson's correlation coefficients. The reference array itself was not normalized. Second, for each experiment the ratio of intensities in comparison to a “virtual reference genome” (VRG) was calculated. The VRG for each “test” is a composite of multiple individuals from the study population and selected to minimize the overall variance in probe ratios. The VRG consisted of the median of the 10 most highly-correlated experiments in the dataset based on the Pearson correlation of 5% of probes on the array. This defined probe set was comprised of every 20th autosomal SNP probe and the nearest adjacent copy number (CN) probe. Separate virtual reference genomes were computed for males and females in order to have a sex-matched VRG.
GC Correction of Log 2 Ratios.
The final step of data processing involved the correction of genomic waves effects in log 2 ratios due to regional correlations with GC content based on the fitted linear regression model proposed by Diskin et al11.
In this study, a two segmentation algorithms was used to discover copy number variants in the GC-corrected NimbleGen HD2 and Affymetrix 6.0 ratio data, HMMSeg12 and Genome Alteration Detection Analysis (GADA)13 (http://biron.usc.edu/˜piquereg/GADA). CNVs detected by both algorithms were used for downstream data processing and analysis. CNVs detected by only one algorithm were excluded. In addition, CNVs of the same type (i.e. deletion or duplication) that were separated by ≦3 probes were merged into one contiguous segment. The proximal and distal boundaries of overlapping and adjacent CNVs were defined by the minimal chromosomal start position and the maximal chromosomal end position of the CNVs. All CNV coordinates are based on the human genome build NCBI36/hg18.
Data filtering criteria were designed to minimize the false discovery rates of rare CNVs in the primary and secondary analysis data sets. In this section the filtering of CNVs and experiments in the NimbleGen HD2 and Affymetrix data sets are described. Further details describing the filtering criteria of rare CNVs are outlined below in Section 6.6.
NimbleGen HD2.
Filtering of NimbleGen HD2 experiments from analysis was based on the experiment quality and filtered CNV properties. Experiments with autosomal median absolute deviations (aMAD)>0.23, with >400 filtered autosomal CNVs and with >30 Mb in combined CNV length (˜1% of the human genome) were removed. In addition, any sample with aneuploidy of autosomal chromosomes was excluded as were experiments with conflicting reported gender and empirically derived X and Y median probe ratios.
CNVs that overlapped regions of the genome prone to somatic cell rearrangements were removed from the primary data set. In particular, CNVs intersecting or overlapping T-cell receptor regions (chr7:38,245,705-38,365,141, chr7:141,647,285-142,221,100, chr9:33,608,462-33,652,656, chr14:21,159,896-22,090,937) and abParts (chr2:88,937,989-89,411,302, chr2:88,966,183-89,377,035, chr2:89,589,457-89,897,555, chr14:105,065,301-106,352,275, chr22:20,715,572-21,595,082) were excluded. CNVs with median probe ratios (seg.median) between 0.83 and 1.15 were also removed as were CNVs less than 100 kb or greater than 10 Mb.
Affymetrix 6.0.
Similar CNV and experiment filtering criteria used to process the primary NimbleGen HD2 dataset were implemented for analysis of CNVs in the secondary data set, taking into account the non-uniformity of Affymetrix 6.0 probe distribution across the human genome. Probable somatic T-cell receptor and abPart rearrangements were removed and CNVs were filtered out based on segment median thresholds. Noisy experiments with Median Absolute Deviation (MAD)>0.2 and chromosomal aneuploidies were eliminated from further analysis.
Within each cohort rare CNV frequencies were determined in the combined set of cases and controls (within each ethnic group separately). CNV frequencies were estimated based on 50% reciprocal overlap between CNV calls of the same type. CNVs with frequency >1% were removed.
Rare CNV calls which passed the 1% frequency filter were further filtered based on the confidence score (CS). Thresholds for CS were adjusted for various size classes of CNVs and for each platform separately as shown in Table 7. As the CS score of a CNV call, the P-value derived from the outlier detection genotyping method, median Z-score outlier detection (MeZOD)14 was used. Appropriate CS thresholds for each size class and platform were determined by examining patterns of Mendelian inheritance in a set of mother-father-child trios (of confirmed parentage), and CS was adjusted to achieve a 5% rate of Mendelian inconsistency for rare CNVs (i.e. 5% of CNVs called in the child were not inherited from a parent). This process maintained an error rate of 5% across a range of CNV sizes.
Regional CNV burdens associated with schizophrenia were analyzed using a two-stage approach. First, using the primary cohort, regions of interest were defined as all genomic segments containing CNVs in at least two cases and in no controls. This discovery step yielded 114 genomic regions of interest. In the secondary cohort of 7,488 patients and 6,689 controls, the association of these regions with schizophrenia were assessed (Table 3). All CNVs overlapping each of the 114 regions of interest were collected, and CNV breakpoints falling within each region were used to partition the region into a series of non-overlapping segments or bins (see
Based on the CNV counts within a segment, association was quantified using the Exact Conditional test, with ancestry and study as covariates. The segment with the lowest one-sided p-value was the peak of association within a ROI. Because segments in different ROIs are driven by the underlying genetic architecture, their numbers and sizes varied widely. Furthermore, numbers of CNVs in nearby segments are highly correlated. To address these issues, a permutation-based p-value correction scheme was applied, where the observed one-sided p-value of the association peak is compared to the distribution of minimal one-sided p-values of any segment within the ROI, computed based on data with case/control labels shuffled at the sample level. This empirical p-value is reported as the p-value for the ROI. In the experiments, p-value was estimated by running 100,000 permutations and p-values of the 8 nominally significant regions were further refined by running 200,000 permutations. For each ROI both the region and peak ORs estimated by the use of Exact Conditional test are reported. In order to avoid infinite OR estimates which occur when the number of controls across all strata is equal to zero, in those situations 0.5 was added to all cells in the contingency table (Haldane correction). In Table 1 and Table 3 Haldane-corrected ORs and the corresponding corrected peak p-values have been marked with an asterisk (*). Within the permutation-based procedure no correction of the counts has been performed, because a mixture of corrected and uncorrected minimal p-values which may occur within a ROI would interfere with the accurate estimation of the null distribution.
Of the 114 regions detected in the first step, 88 had events in the secondary cohort. Four regions had statistically significant associations in the secondary sample after Bonferroni correction (α=0.05/114=4.4×10−4). Table 1 lists the four regions with significant p-values under this criterion and an additional four loci with nominally significant p-values (4.4×10−4−0.05) in the secondary cohort. Regions with the significant associations were loss of copy number at 22q11.2 (P<5×10−6, OR=14.2), gain at 7q36.3 (P=4.0×10−5, OR=16.4), gain at 16p11.2 (P=1.0×10−4, OR=16.1) and loss at 15q13.3 (P=1.6×10−4, OR=14.9). No significant heterogeneity was observed for these genomic regions across the three studies (Table 1).
15q13.3, 16p11.2 and 22q11.2 are well-documented loci conferring increased risk for schizophrenia23, 24, 30. All are hotspots for non-allelic homologous recombination (NAHR), and all alleles contributing to the association consist of large deletions with similar breakpoints. By contrast, microduplications at 7q36.3 have not been previously implicated in neuropsychiatric disorders. The 7q36.3 region harbored CNVs that overlapped but differed in size and breakpoint positions (
Sequenom MassArray. Absolute copy number (ACN) detection method by Sequenom MassArray combines real time competitive PCR (rPCR) with MassEXTEND procedures and matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF) (http://www.sequenom.com). The CNV assay involves spiking genomic DNA template with genomic DNA from a single reference chimpanzee (Pan troglodytes) that was purchased from Southwest National Primate Research Center, (San Antonio, Tex., USA). The genotype assay targets a specific single nucleotide difference between human and chimpanzee, and relative copy number is determined based the ratio of peak areas of the human and chimp alleles.
CNV Assay Design: CNV assays were designed for two segments (proximal and distal) of the 7q36.3 region as shown in
Sequence differences between human and chimpanzee were identified within CNV target regions by aligning the respective genomic DNA using UCSC Blat (http://genome.ucsc.edu). Nucleotide differences between human and chimpanzee sequences were then categorized based on the position of difference in the alignment, alleles in human and chimpanzee, and direction of alignment. All aligned human-chimpanzee loci were processed to identify the location of any variant bases (single nucleotide polymorphisms and insertions/deletions) within a given distance to the specified human-chimpanzee loci. All known human and chimpanzee SNPs and indels were excluded in CNV assay design. Single base nucleotide extension (SBE) assays were then designed to target the non-conserved nucleotides using Sequenom Assay Design v3.1 software. The primers used for Sequenom assays are listed in Table 7.
PCR: PCR was performed on Biorad thermocyclers in a 5 μl reaction using 15 ng of genomic human DNA (hsDNA), 15 ng of P. troglodytes DNA (ptDNA), 0.5 pmol of each primer in multiplex, 250 μM dNTPs, 0.5 μl of 10×20 mM MgCl PCR Buffer and 0.5 units of Taq polymerase (Roche, USA). The PCR cycling parameters were 95° C. for 15 min, followed by 45 cycles of 94° C. for 30 s, 56° C. for 30 s and 72° C. for 1 min, followed by 72° C. for 3 min. Following PCR, 5 μl PCR product was treated with 0.17 μl of 10×SAP (shrimp alkaline phosphate) buffer and 0.30 μl of SAP enzyme and incubated at 37° C. for 40 min followed by heat denaturation of enzyme at 85° C. for 5 min.
Single Base Primer Extension: Primer extension was performed on MJ thermocyclers in 9 μl reactions using the 7 μl SAP-treated PCR product, 7-14 pmol of each primer in multiplex, 0.20 μl of iplex buffer, 0.20 μl of iplex termination mix and 0.05 μl of iplex enzyme (Sequenom, San Diego, Calif., USA). The cycling parameters were 94° C. for 30 s, followed by 40 cycles of 94° C. for 5 s and a nested 5 cycles of 52° C. for 5 s and 80° C. for 5 s, followed by 72° C. for 3 min.
Quantification of absolute copy number: The absolute copy number for target regions was determined using the method described by Williams et al15. All copy number measurements from each sample were normalized against copy number measurements from the assays targeting the invariant (normal diploid copy) regions of the genome in the same sample to control for any sample to sample (i.e. well to well) loading variation. Data were analyzed using Splus 8.0.
As outlined in detail below, sensitivity and specificity of CNV calls in the 7q36 region were examined to determine the possibility of a spurious association. No additional duplications >100 kb were detected after reducing the stringency of the CNV filtering criteria. Second, a more sensitive targeted CNV calling algorithm, median Z-score Outlier Detection (MeZOD)30, was applied to obtain CNV genotype calls for the proximal and distal regions of 7q36.3. All 16 copy number gains, nine overlapping the proximal region (
Unexpectedly, manual examination of probe ratios in
Inheritance of the duplication at 7q36.3 could be evaluated in three families (
Fine Mapping of 7q36.3 Duplications in MGS Subjects Using the Nimblegen HD2 Array.
Duplications of 7q36.3 detected by the Affymetrix 6.0 array in the MGS sample were mapped at finer resolution using the Nimblegen HD2 platform. (
Based on these results (
Fluorescence In Situ Hybridization (FISH).
Cytogenetic confirmation was obtained for 2 samples with duplication of VIPR2 (05C43079 and 05C48386). Subtelomeric probes were obtained for chromosome 7p and 7q (Abbott Molecular). The genomic coordinates of the TelVysion 7q probe were Chr7:158551912-158650925 (hg18), which overlaps exons 1-4 of VIPR2. Hybridizations were performed according to manufacturer's protocols.
In images from fluorescence in situ hybridization (FISH) using a subtelomeric 7q probe, three distinct red signals can be seen in interphase cells, confirming a heterozygous duplication. In metaphase cells, all signals localize to the subtelomeric region of 7q, confirming that the duplications lie adjacent to each other in the 7q36.3 region (
Sensitivity of CNV calls in the 7q36.3 region were examined to determine the possibility of a spurious association. A spurious association can arise in CNV data, for instance, if there is comparatively reduced sensitivity to detect CNVs in the control sample compared with the case sample. Sensitivity of the segmentation-based CNV calling methods was evaluated by first examining CNVs within the 7q36.3 region after relaxing the stringent filtering criteria, and then by comparing segmentation calls in cases and controls to CNV genotype calls made in the 7q36.3 region and elsewhere in the genome using targeted genotyping algorithms with enhanced sensitivity. These efforts were focused only on the microarray datasets that were processed for this study. Published CNV call sets generated by the ISC study of schizophrenia have undergone evaluation in earlier studies8.
CNV Segmentation Calls in the 7q36.3 Region Examined Using a Lower Sensitivity Threshold.
After eliminating the primary filtering criterion, the CNV confidence score (CS), no additional CNVs were detected in the 7q36.3 region; therefore, filtering based on confidence score does not account for the observed differences in CNV frequencies. Second, rare CNVs in the 7q36.3 region were examined after relaxing the minimum size to 5 probes and 16 probes in the Nimblegen and Affymetrix data respectively. Additional small CNVs were detected that did not negatively influence the observed association (
Concordance Between Segmentation Calls at 706.3 Region and CNV Calls by a More Sensitive Outlier-Detection Based CNV Detection Method (MeZOD).
As a second approach to evaluating the sensitivity of the disclosed methods, a more sensitive outlier-detection based method, Median Z-score Outlier Detection (MeZOD), was applied to the analysis of the two CNV target regions at 7q36.3. Outlier detection-based methods14, 16 provide greater sensitivity for CNVs within defined target regions. If a significant number of CNVs overlapping the target region were undetected by the segmentation algorithms or filtered out of the call set, such variants can often be detected by MeZOD, leading to discordance between the segmentation calls and the MeZOD genotypes. When the distribution of z-scores for each target region was examined, punctuate clusters of data points were observed that are distinguishable from the overall distribution (displayed as cluster plots in
Sensitivity to Detect CNPs is Comparable in Cases and Controls.
Sensitivity at other loci throughout the genome was evaluated by examining sensitivity to detect a set of validated common CNPs (>100 Kb in size) characterized as part of HapMap phase 317. Sensitivity to detect common CNPs by the segmentation algorithms is equivalent to that of rare CNVs (because parameters of HMM-based segmentation algorithms are not adjusted based on prior knowledge of common CNVs). However, targeted genotyping methods have much greater sensitivity and accuracy. Therefore concordance between the segmentation calls and a predefined set of common CNP genotypes was used as a measure of the segmentor sensitivity,
First, all common CNPs >100 Kb were selected from HapMap 3 (http://hapmap.ncbi.nlm.nih.gov/downloads/cnv_data/hm3_cnv_submission.txt), and these markers were compared with a list of common CNP markers that were previously assembled for the NimbleGen HD2 platform. All CNPs that were present in both sets (total of 7) were selected. The seven CNPs were then genotyped by manual examination of clusters as illustrated in
Examining the Impact of Differential Sensitivity of Platforms Used in Different Studies.
The primary dataset consisted of NimbleGen 2.1M array data. The secondary dataset consisted of the MGS study (Affymetrix 6.0) and the ISC study of schizophrenia (which included data from Affymetrix 5.0 and 6.0 platforms). Combining data from multiple studies has the potential to create spurious associations. For instance, if data from two studies are combined and the study with the most sensitive platform had a large number of cases and a small number controls, a spurious enrichment of CNVs in cases could arise (if factors such as platform and study are not controlled for). There was no such a disparity in any of the 3 studies comprising the primary and secondary cohorts, which consisted of roughly equal numbers of cases and controls (Table 2). In the conditional Fisher's exact-based permutation test for association, platform and study were included as factors. Second, heterogeneity in the association across all 3 studies was tested for in the primary and secondary samples; no significant heterogeneity was observed across the 3 populations when tested using the Breslow-Day Tarone test (P=0.19−0.86). Therefore, evidence that the observed association can be attributed to differential sensitivity of platforms used across studies was not found.
Levels of VIPR2 mRNA were measured in lymphoblastoid cell lines obtained from the MGS study collection at the NIMH genetics repository (http://www.nimhgenetics.org), including patients with duplications and triplications of 7q36.3 and control individuals with normal diploid copy number of this region. The samples and genotypes used in these experiments were as follows: controls: 04C34751, 04C28567, 05C47472, 04C34070; duplication of VIPR2: 05C44574, 05C48386, 05C48694; Duplication of VIPR2 with a nested triplication of exons 3 and 4: 00C2204, 03C23250, 05C43079, 03C23091; duplication of upstream subtelomeric region: 05C51123, 05C46770. CNV genotypes were confirmed in cell lines by Sequenom assay (see Section 6.8 and
Cell lines from two subjects with each genotype were analyzed. From each cell line, a series of two RNA preparations were made in order to provide two biological replicates of each subject. Each biological replicate was in turn analyzed with triplicate qPCR reactions. Results presented for each subject represent the mean and standard error of 6 replicates. Standard error and P-values for pooled results were computed across individuals (
Tissue Culture.
EBV-transformed lymphoblastoid cells from MGS patients and controls were obtained from the NIMH genetics initiative repository. Cell lines were cultured in Gibco's RPMI-1640 supplemented with 15% FBS and 1× Penn-Strep at 37° C. and 5% CO2.
Sample Preparation.
Total RNA was prepared using Qiagen's RNEasy Plus kit, and cDNA was prepared using Quanta Qscript cDNA mix using 1 ug RNA for each 20 μl reaction.
Quantitative PCR (qPCR).
Measurement of VIPR2 mRNA abundance was carried out on mRNA from lymphoblastoid lines using the VIPR2 TaqMan® assay (Applied Biosystems Inc.). Assays were carried out by the Center for Aids Research (CFAR) core facility at UCSD (http://cfar.ucsd.edu/core-facilityes/genomics-core/protocols-and-resources-1/quantitative-real-time-rt-per-sample-preparation). cDNA from 500 ng of total RNA starting material was used for each qPCR reaction. Samples were done in triplicates using the VIPR2 and GAPDH primers provided by the manufacturer. qPCR cycling conditions were as follows: 50° C. for 2 min, 95° C. for 10 min, then 45 cycles of 95° C. for 15 sec, 60° C. for 1 min. Ct values obtained were normalized with GAPDH Ct values. VIPR2 mRNA was detectable in all cell lines analyzed in this study (
Low but measurable levels of VIPR2 expression was observed in all lymphoblastoid cell lines in this study (the range of Ct values was 33-36). This is consistent with the low VIPR2 transcript abundance in the published transcriptome data from the cell lines of 60 Nigerian subjects by Pickrell et al (21 reads/968 million) 18. This is also consistent with the very low level of VPAC2 receptor signaling observed in control cell lines.
In cell lines from patients carrying duplications of terminal 7q, an approximately 3-fold increase in VIPR2 expression relative to controls was observed (
Cyclic-AMP signaling has been implicated in schizophrenia51, 52. It was determined whether increases in VIPR2 transcription and VPAC2-mediated cAMP signaling was a consequence of the microduplications at 7q36.3. VIPR2 mRNA and cAMP accumulation in response to VIP and a VPAC2-selective agonist (BAY 55-9837) in lymphoblastoid cell lines from eight MGS study subjects was assessed: two with subtelomeric duplications, three with duplications of VIPR2, four with partial triplications, and four controls with normal copy number of the region.
Cyclic AMP accumulation was measured in lymphoblastoid cell lines (0.5 million per ml) pre-incubated for 20 minutes with the cyclic nucleotide phosphodiesterase inhibitor isobutylmethylxanthine (IBMX, 200 μM), before the addition of the stimulatory agonists forskolin (10 μM)+/−VIP, [100 nM], BAY 55-9837 (100 nM) or prostaglandin E2 [PGE2, 1 μM] for 10 min. Reactions were terminated by pelleting the cells, aspiration of the medium and addition of 100 μl of cold 7.5% (wt/vol) trichloroacetic acid (TCA). Cyclic AMP content in TCA extracts was determined by radioimmunoassay and normalized to the number of cells per well. Data are expressed as cAMP accumulation in response to the GPCR agonists relative to the response to non-GPCR agonist forskolin (10 μM) and IBMX (200 mM) alone. Results presented for each subject represents the mean and standard error of at least ten replicates. Standard error and P-values for pooled results were computed across individuals (
VIPR2 transcripts were present at detectable levels in all cell lines. VIPR2 mRNA levels were significantly increased in duplication carriers compared to controls (
In cell lines from patients carrying duplications of terminal 7q, a significant increase in cyclic-AMP accumulation in response to two different VPAC2 agonists relative to controls was observed (
In four MGS patients, complex duplication/triplication were observed consisting of a duplication of the entire gene and another duplication of exons 3 and 4 nested within it. The exact structure of this haplotype is not known. If the nested duplication of exons 3 and 4 lies adjacent and in a direct orientation, then this rearrangement may result in a mutant transcript with tandem duplication of exons 3 and 4 which would result in a shift of the open reading frame. However, if the nested duplication lies outside of the gene or in an inverted orientation, then nested duplication would not have an obvious impact on either of the two copies of VIPR2 (thus, the duplication/triplication might be functionally equivalent to a tandem duplication).
Using primers designed for exons 2 and 5, VIPR2 cDNA products were sequenced from the cell lines listed in
Both patients were from a family triad sample collected from the US5. All patients in the sample participated in a formal diagnostic interview conducted by an experienced, licensed clinician specially trained in the use of the Diagnostic Interview for Genetic Studies (DIGS). Diagnostic criteria were applied as described in the Diagnostic and Statistical Manual 4 (DSM-IV). All diagnoses were confirmed by a senior licensed clinical psychologist.
Subject 02-0016.
Subject 02-0016 is male, age 44, of Norwegian descent. His diagnosis is schizoaffective disorder, depressed type, with onset at the age of 21. Family history: mother suffers from depression. Patient has never been married. He completed 4 years of college. He currently resides in a halfway house and attends a day treatment program. He is unemployed and receives disability. At age 16, he reported first feeling paranoid and having trouble in school. At age 21, he reported experiencing his first psychotic break and was hospitalized for 4 months. His symptoms primarily entailed paranoia. He also reported being depressed during his hospitalization. At age 22, he was hospitalized for a second time for suicidal thoughts, depression and feelings of hopelessness. At age 26, he experienced a relapse with prominent paranoid delusions. The patient had a history of hyperactivity prior to the age of 16. Developmentally, he reported maturation lag (i.e., growing up more slowly than other children). He also reported motor and vocal tics.
Subject 02-0135.
Subject 02-0135 is female, age 50, of Irish/German/English/Scottish descent. Her diagnosis is schizoaffective disorder, bipolar type, with onset at the age of 22. Family history: None known. Patient has never been married and her education level includes some college course work. She currently resides in a semi-independent group home and has been unemployed due to her psychiatric problems. Her primary psychiatric symptoms are paranoia and irritability. She began having significant psychiatric problems when she was 21 while she was away from home at college. She became aggressive and loud and started showing poor judgment. She has been hospitalized numerous times (at least 8) and her functioning has deteriorated steadily over the years. She has experienced at least 4 manic episodes, each lasting several months, with the first occurring when she was 25. The patient exhibited pressured speech, irritable mood, and grandiose ideation during each manic episode. She also reported obsessions with patterns but did not meet diagnostic criteria for OCD.
Developmentally, the patient has a history of delayed speech (did not speak words until the age of 3) and delayed puberty (did not menstruate until college). She lacked close friendships during childhood.
Variable expressivity is often characteristic of pathogenic CNVs30. The spectrum of psychiatric phenotypes associated with 7q36.3 duplications was evaluated by screening for these events in individuals with bipolar disorder or autism spectrum disorder (ASD). Microarray data was available for 2,777 patients from the Bipolar Genome Study (BiGS), for 996 ASD patients from the Autism Genome Project Consortium (AGP), and from unpublished analyses of 114 patients with ASD using the NimbleGen HD2 platform. Microduplications of 7q36.3 (>100 Kb in size) were detected in three of 1,110 (0.27%) of patients with ASD; compared with the controls described above, P=0.018. Microduplications at 7q36.3 were detected in two of 2,777 (0.07%) patients with bipolar disorder; compared with the controls, P=0.23. These results provide preliminary evidence that the clinical phenotypes associated with 7q36.3 duplications may include pediatric neurodevelopmental disorders such as autism, but do not include bipolar disorder. Also worthy of note, larger chromosomal abnormalities involving 7q have been reported in association with neurodevelopmental disorders, including deletions of 7q36-7qter36, 37 and duplications of 7q35-7qter.38; a 550 Kb duplication of 7qter (of unknown clinical relevance) has been reported in a patient with neurofibromatosis39.
These genetic data implicate the gene VIPR2. All variants contributing to the observed association at 7q36.3 overlap with this gene or lie within the gene-less subtelomeric region <89 kb from the transcriptional start site of VIPR2. In addition, partial triplications of VIPR2 involving only exons 3 and 4 were detected in four patients and zero controls. VIPR2 encodes the Vasoactive Intestinal Peptide (VIP) Receptor VPAC2, a G protein-coupled receptor that is expressed in variety of tissues including, in the brain, the suprachiasmatic nucleus, hippocampus, amygdala, and hypothalamus40. VPAC2 binds VIP41, activates cyclic-AMP signaling and PKA, regulates synaptic transmission in the hippocampus42-44, and promotes the proliferation of neural progenitor cells in the dentate gyrus45. Genetic studies in mouse have established that VIP signaling plays a role in learning and memory46. VPAC2 also plays a role in sustaining normal circadian oscillations in the SCN47, 48, and VIPR2-null49 and VIPR2-overexpression50 mice exhibit abnormal rhythms of rest and activity.
The observed expression patterns suggest that a variety of different genomic duplications can influence the transcription of VIPR2. Without being bound by theory, it is thought that duplications of 7q36 impact the regulation of VIPR2. Tandem duplication of regulatory sequences, for instance, could affect expression of the gene. Alternatively, the subtelomeric location of VIPR2 could be relevant to the mechanism. Intrinsic regulation of telomere structure and function often impacts the transcriptional regulation of adjacent genes, a phenomenon known as Telomere Position Effect (TPE)53, 54. If VIPR2 is under such epigenetic regulation, any large tandem duplication of the subtelomeric region could potentially cause the gene to escape repression.
The present study implicates VIPR2 as a genetic susceptibility factor for schizophrenia. It has been shown that increased copy number of VIPR2 or the adjacent region is highly enriched in patients with schizophrenia. Moreover, disease mutations lead to increased VIPR2 transcription and VIP-induced cyclic-AMP signaling in patient cells. These results implicate VIP signaling as a molecular mechanism underlying schizophrenia. In light of the emerging roles of VIPR2 in the brain, the results described herein support the contention that the pathogenesis of schizophrenia, in some patients, involves the dysregulation of cellular processes such as adult neurogenesis and synaptic transmission and of the corresponding cognitive processes of learning and memory. Furthermore, in light of the brain expression patterns of VIPR240, the present results support the involvement of certain brain regions, such as hippocampus, amygdala and suprachiasmatic nucleus.
The link between VIPR2 duplications and schizophrenia may have significant implications for the development of molecular diagnostics and treatments for this disorder. Genetic testing for duplications of the 7q36 region could enable the early detection of a subtype of patients characterized by overexpression of VIPR2. Significant potential also exists for the development of therapeutics targeting this receptor. For instance, a selective antagonist of VPAC2R could have therapeutic potential in patients who carry duplications of the VIPR2 region. Peptide derivatives55, 56 and small molecules57 have been identified that are selective VPAC2 inhibitors, and these pharmacological studies offer potential leads in the development of new drugs.
All patents, patent applications, publications, product descriptions and protocols, cited in this specification are hereby incorporated by reference in their entirety. In case of a conflict in terminology, the present disclosure controls.
While it will be apparent that the disclosure herein described is well calculated to achieve the benefits and advantages set forth above, the present disclosure is not to be limited in scope by the specific embodiments described herein. It will be appreciated that the disclosure is susceptible to modification, variation and change without departing from the spirit thereof.
This application is a continuation of International Application No. PCT/US12/020683 filed Jan. 9, 2012, which claims the benefit of and priority to U.S. Provisional Application Ser. No. 61/430,883, filed Jan. 7, 2011, both of which are hereby incorporated by reference in their entireties, and to each of which priority is claimed.
This invention was made with government support under grant number NIH MH061399 awarded by the National Institutes of Health. The government has certain rights in the invention.
Number | Date | Country | |
---|---|---|---|
61430883 | Jan 2011 | US |
Number | Date | Country | |
---|---|---|---|
Parent | PCT/US2012/020683 | Jan 2012 | US |
Child | 13935248 | US |