Genomic subtyping of liver cancers with prognostic application.

BACKGROUND: Cancer subtyping has mainly relied on pathological and molecular means. Massively parallel sequencing-enabled subtyping requires genomic markers to be developed based on global features rather than individual mutations for effective implementation. METHODS: In the present study, the whole genome sequences (WGS) of 110 liver cancers of Japanese patients published with different pathologies were analyzed with respect to their single nucleotide variations (SNVs) comprising both gain-of-heterozygosity (GOH) and loss-of-heterozygosity (LOH) mutations, the signatures of combined GOH and LOH mutations, along with recurrent copy number variations (CNVs). RESULTS: The results, obtained based on the WGS sequences as well as the Exome subset within the WGSs that covered ~ 2.0% of the WGS and the AluScan-subset within the WGSs that were amplifiable by Alu element-consensus primers and covered ~ 2.1% of the WGS, indicated that the WGS samples could be employed with the mutational parameters of SNV load, LOH%, the Signature α%, and survival-associated recurrent CNVs (srCNVs) as genomic markers for subtyping to stratify liver cancer patients prognostically into the long and short survival subgroups. The usage of the AluScan-subset data, which could be implemented with sub-micrograms of DNA samples and vastly reduced sequencing analysis task, outperformed the usage of WGS data when LOH% was employed as stratifying criterion. CONCLUSIONS: Thus genomic subtyping performed with novel genomic markers identified in this study was effective in predicting patient-survival duration, with cohorts of hepatocellular carcinomas alone and those including intrahepatic cholangiocarcinomas. Such relatively heterogeneity-insensitive genomic subtyping merits further studies with a broader spectrum of cancers.

Noninvasive prenatal paternity testing by means of SNP-based targeted sequencing.

OBJECTIVE: To develop a method for noninvasive prenatal paternity testing based on targeted sequencing of single nucleotide polymorphisms (SNPs). METHOD: SNPs were selected based on population genetics data. Target-SNPs in cell-free DNA extracted from maternal blood (maternal cfDNA) were analyzed by targeted sequencing wherein target enrichment was based on multiplex amplification using QIAseq Targeted DNA Panels with Unique Molecular Identifiers. Fetal SNP genotypes were called using a novel bioinformatics algorithm, and the combined paternity indices (CPIs) and resultant paternity probabilities were calculated. RESULTS: Fetal SNP genotypes obtained from targeted sequencing of maternal cfDNA were 100% concordant with those from amniotic fluid-derived fetal genomic DNA. From an initial panel of 356 target-SNPs, an average of 148 were included in paternity calculations in 15 family trio cases, generating paternity probabilities of greater than 99.9999%. All paternity results were confirmed by short-tandem-repeat analysis. The high specificity of the methodology was validated by successful paternity discrimination between biological fathers and their siblings and by large separations between the CPIs calculated for the biological fathers and those for 60 unrelated men. CONCLUSION: The novel method is highly effective, with substantial improvements over similar approaches in terms of reduced number of target-SNPs, increased accuracy, and reduced costs.

Forward and reverse mutations in stages of cancer development.

BACKGROUND: Massive occurrences of interstitial loss of heterozygosity (LOH) likely resulting from gene conversions were found by us in different cancers as a type of single-nucleotide variations (SNVs), comparable in abundance to the commonly investigated gain of heterozygosity (GOH) type of SNVs, raising the question of the relationships between these two opposing types of cancer mutations. METHODS: In the present study, SNVs in 12 tetra sample and 17 trio sample sets from four cancer types along with copy number variations (CNVs) were analyzed by AluScan sequencing, comparing tumor with white blood cells as well as tissues vicinal to the tumor. Four published "nontumor"-tumor metastasis trios and 246 pan-cancer pairs analyzed by whole-genome sequencing (WGS) and 67 trios by whole-exome sequencing (WES) were also examined. RESULTS: Widespread GOHs enriched with CG-to-TG changes and associated with nearby CNVs and LOHs enriched with TG-to-CG changes were observed. Occurrences of GOH were 1.9-fold higher than LOH in "nontumor" tissues more than 2 cm away from the tumors, and a majority of these GOHs and LOHs were reversed in "paratumor" tissues within 2 cm of the tumors, forming forward-reverse mutation cycles where the revertant LOHs displayed strong lineage effects that pointed to a sequential instead of parallel development from "nontumor" to "paratumor" and onto tumor cells, which was also supported by the relative frequencies of 26 distinct classes of CNVs between these three types of cell populations. CONCLUSIONS: These findings suggest that developing cancer cells undergo sequential changes that enable the "nontumor" cells to acquire a wide range of forward mutations including ones that are essential for oncogenicity, followed by revertant mutations in the "paratumor" cells to avoid growth retardation by excessive mutation load. Such utilization of forward-reverse mutation cycles as an adaptive mechanism was also observed in cultured HeLa cells upon successive replatings. An understanding of forward-reverse mutation cycles in cancer development could provide a genomic basis for improved early diagnosis, staging, and treatment of cancers.

Gabrb2-knockout mice displayed schizophrenia-like and comorbid phenotypes with interneuron-astrocyte-microglia dysregulation.

Intronic polymorphisms of the GABAA receptor ß2 subunit gene (GABRB2) under adaptive evolution were associated with schizophrenia and reduced expression, especially of the long isoform which differs in electrophysiological properties from the short isoform. The present study was directed to examining the gene dosage effects of Gabrb2 in knockout mice of both heterozygous (HT) and homozygous (KO) genotypes with respect to possible schizophrenia-like and comorbid phenotypes. The KO mice, and HT mice to a lesser extent, were found to display prepulse inhibition (PPI) deficit, locomotor hyperactivity, stereotypy, sociability impairments, spatial-working and spatial-reference memory deficits, reduced depression and anxiety, and accelerated pentylenetetrazol (PTZ)-induced seizure. In addition, the KO mice were highly susceptible to audiogenic epilepsy. Some of the behavioral phenotypes showed evidence of imprinting, gender effect and amelioration by the antipsychotic risperidone, and the audiogenic epilepsy was inhibited by the antiepileptic diazepam. GABAergic parvalbumin (PV)-positive interneuron dystrophy, astrocyte dystrophy, and extensive microglia activation were observed in the frontotemporal corticolimbic regions, and reduction of newborn neurons was observed in the hippocampus by immunohistochemical staining. The neuroinflammation indicated by microglial activation was accompanied by elevated brain levels of oxidative stress marker malondialdehyde (MDA) and the pro-inflammatory cytokines tumor necrosis factor-alpha (TNF-α) and interleukin-6 (IL-6). These extensive schizophrenia-like and comorbid phenotypes brought about by Gabrb2 knockout, in conjunction with our previous findings on GABRB2 association with schizophrenia, support a pivotal role of GABRB2 in schizophrenia etiology.

Variation of global DNA methylation levels with age and in autistic children.

BACKGROUND: The change in epigenetic signatures, in particular DNA methylation, has been proposed as risk markers for various age-related diseases. However, the course of variation in methylation levels with age, the difference in methylation between genders, and methylation-disease association at the whole genome level is unclear. In the present study, genome-wide methylation levels in DNA extracted from peripheral blood for 2116 healthy Chinese in the 2-97 age range and 280 autistic trios were examined using the fluorescence polarization-based genome-wide DNA methylation quantification method developed by us. RESULTS: Genome-wide or global DNA methylation levels proceeded through multiple phases of variation with age, consisting of a steady increase from age 2 to 25 (r = 0.382) and another rise from age 41 to 55 to reach a peak level of ~80 % (r = 0.265), followed by a sharp decrease to ~40 % in the mid-1970s (age 56 to 75; r = -0.395) and leveling off thereafter. Significant gender effect in methylation levels was observed only for the 41-55 age group in which methylation in females was significantly higher than in males (p = 0.010). In addition, global methylation level was significantly higher in autistic children than in age-matched healthy children (p < 0.001). CONCLUSIONS: The multiphasic nature of changes in global methylation levels with age was delineated, and investigation into the factors underlying this profile will be essential to a proper understanding of the aging process. Furthermore, this first report of global hypermethylation in autistic children also illustrates the importance of age-matched controls in characterization of disease-associated variations in DNA methylation.

Feature co-localization landscape of the human genome.

Although feature co-localizations could serve as useful guide-posts to genome architecture, a comprehensive and quantitative feature co-localization map of the human genome has been lacking. Herein we show that, in contrast to the conventional bipartite division of genomic sequences into genic and inter-genic regions, pairwise co-localizations of forty-two genomic features in the twenty-two autosomes based on 50-kb to 2,000-kb sequence windows indicate a tripartite zonal architecture comprising Genic zones enriched with gene-related features and Alu-elements; Proximal zones enriched with MIR- and L2-elements, transcription-factor-binding-sites (TFBSs), and conserved-indels (CIDs); and Distal zones enriched with L1-elements. Co-localizations between single-nucleotide-polymorphisms (SNPs) and copy-number-variations (CNVs) reveal a fraction of sequence windows displaying steeply enhanced levels of SNPs, CNVs and recombination rates that point to active adaptive evolution in such pathways as immune response, sensory perceptions, and cognition. The strongest positive co-localization observed between TFBSs and CIDs suggests a regulatory role of CIDs in cooperation with TFBSs. The positive co-localizations of cancer somatic CNVs (CNVT) with all Proximal zone and most Genic zone features, in contrast to the distinctly more restricted co-localizations exhibited by germline CNVs (CNVG), reveal disparate distributions of CNVTs and CNVGs indicative of dissimilarity in their underlying mechanisms.

GABRB2 Haplotype Association with Heroin Dependence in Chinese Population.

Substance dependence is a frequently observed comorbid disorder in schizophrenia, but little is known about genetic factors possibly shared between the two psychotic disorders. GABRB2, a schizophrenia candidate gene coding for GABAA receptor ß2 subunit, is examined for possible association with heroin dependence in Han Chinese population. Four single nucleotide polymorphisms (SNPs) in GABRB2, namely rs6556547 (S1), rs1816071 (S3), rs18016072 (S5), and rs187269 (S29), previously associated with schizophrenia, were examined for their association with heroin dependence. Two additional SNPs, rs10051667 (S31) and rs967771 (S32), previously associated with alcohol dependence and bipolar disorder respectively, were also analyzed. The six SNPs were genotyped by direct sequencing of PCR amplicons of target regions for 564 heroin dependent individuals and 498 controls of Han Chinese origin. Interestingly, it was found that recombination between the haplotypes of all-derived-allele (H1; OR = 1.00) and all-ancestral-allele (H2; OR = 0.74) at S5-S29 junction generated two recombinants H3 (OR = 8.51) and H4 (OR = 5.58), both conferring high susceptibility to heroin dependence. Additional recombination between H2 and H3 haplotypes at S1-S3 junction resulted in a risk-conferring haplotype H5 (OR = 1.94x109). In contrast, recombination between H1 and H2 haplotypes at S3-S5 junction rescued the risk-conferring effect of recombination at S5-S29 junction, giving rise to the protective haplotype H6 (OR = 0.68). Risk-conferring effects of S1-S3 and S5-S29 crossovers and protective effects of S3-S5 crossover were seen in both pure heroin dependent and multiple substance dependence subgroups. In conclusion, significant association was found with haplotypes of the S1-S29 segment in GABRB2 for heroin dependence in Han Chinese population. Local recombination was an important determining factor for switching haplotypes between risk-conferring and protective statuses. The present study provide evidence for the schizophrenia candidate gene GABRB2 to play a role in heroin dependence, but replication of these findings is required.

SAMSVM: A tool for misalignment filtration of SAM-format sequences with support vector machine.

Sequence alignment/map (SAM) formatted sequences [Li H, Handsaker B, Wysoker A et al., Bioinformatics 25(16):2078-2079, 2009.] have taken on a main role in bioinformatics since the development of massive parallel sequencing. However, because misalignment of sequences poses a significant problem in analysis of sequencing data that could lead to false positives in variant calling, the exclusion of misaligned reads is a necessity in analysis. In this regard, the multiple features of SAM-formatted sequences can be treated as vectors in a multi-dimension space to allow the application of a support vector machine (SVM). Applying the LIBSVM tools developed by Chang and Lin [Chang C-C, Lin C-J, ACM Trans Intell Syst Technol 2:1-27, 2011.] as a simple interface for support vector classification, the SAMSVM package has been developed in this study to enable misalignment filtration of SAM-formatted sequences. Cross-validation between two simulated datasets processed with SAMSVM yielded accuracies that ranged from 0.89 to 0.97 with F-scores ranging from 0.77 to 0.94 in 14 groups characterized by different mutation rates from 0.001 to 0.1, indicating that the model built using SAMSVM was accurate in misalignment detection. Application of SAMSVM to actual sequencing data resulted in filtration of misaligned reads and correction of variant calling.

Glioma Association and Balancing Selection of ZFPM2.

ZFPM2, encoding a zinc finger protein and abundantly expressed in the brain, uterus and smooth muscles, plays important roles in cardiac and gonadal development. Abnormal expression of ZFPM2 in ovarian tumors and neuroblastoma has been reported but hitherto its genetic association with cancer and effects on gliomas have not been studied. In the present study, the hexamer insertion-deletion polymorphism rs71305152, located within a large haplotype block spanning intron 1 to intron 3 of ZFPM2, was genotyped in Chinese cohorts of glioma (n = 350), non-glioma cancer (n = 354) and healthy control (n = 463) by direct sequencing and length polymorphism in gel electrophoresis, and ZFPM2 expression in glioma tissues (n = 69) of different grades was quantified by real-time RT-PCR. Moreover, potential natural selection pressure acting on the gene was investigated. Disease-association analysis showed that the overall genotype of rs71305152 was significantly associated with gliomas (P = 0.016), and the heterozygous genotype compared to the combined homozygous genotypes was less frequent in gliomas than in controls (P = 0.005) or non-glioma cancers (P = 0.020). ZFPM2 mRNA expression was negatively correlated with the grades of gliomas (P = 0.002), with higher expression levels in the low-grade gliomas. In the astrocytoma subtype, higher ZFPM2 expression was also correlated with the rs71305152 heterozygous genotype (P = 0.028). In addition, summary statistics tests gave highly positive values, demonstrating that the gene is under the influence of balancing selection. These findings suggest that ZFPM2 is a glioma susceptibility gene, its genotype and expression showing associations with incidence and severity, respectively. Moreover, the balancing selection acting on ZFPM2 may be related to the important roles it has to play in multiple organ development or associated disease etiology.

Massive interstitial copy-neutral loss-of-heterozygosity as evidence for cancer being a disease of the DNA-damage response.

BACKGROUND: The presence of loss-of-heterozygosity (LOH) mutations in cancer cell genomes is commonly encountered. Moreover, the occurrences of LOHs in tumor suppressor genes play important roles in oncogenesis. However, because the causative mechanisms underlying LOH mutations in cancer cells yet remain to be elucidated, enquiry into the nature of these mechanisms based on a comprehensive examination of the characteristics of LOHs in multiple types of cancers has become a necessity. METHODS: We performed next-generation sequencing on inter-Alu sequences of five different types of solid tumors and acute myeloid leukemias, employing the AluScan platform which entailed amplification of such sequences using multiple PCR primers based on the consensus sequences of Alu elements; as well as the whole genome sequences of a lung-to-liver metastatic cancer and a primary liver cancer. Paired-end sequencing reads were aligned to the reference human genome to identify major and minor alleles so that the partition of LOH products between homozygous-major vs. homozygous-minor alleles could be determined at single-base resolution. Strict filtering conditions were employed to avoid false positives. Measurements of LOH occurrences in copy number variation (CNV)-neutral regions were obtained through removal of CNV-associated LOHs. RESULTS: We found: (a) average occurrence of copy-neutral LOHs amounting to 6.9% of heterologous loci in the various cancers; (b) the mainly interstitial nature of the LOHs; and (c) preference for formation of homozygous-major over homozygous-minor, and transitional over transversional, LOHs. CONCLUSIONS: The characteristics of the cancer LOHs, observed in both AluScan and whole genome sequencings, point to the formation of LOHs through repair of double-strand breaks by interhomolog recombination, or gene conversion, as the consequence of a defective DNA-damage response, leading to a unified mechanism for generating the mutations required for oncogenesis as well as the progression of cancer cells.

Copy number variation analysis based on AluScan sequences.

BACKGROUND: AluScan combines inter-Alu PCR using multiple Alu-based primers with opposite orientations and next-generation sequencing to capture a huge number of Alu-proximal genomic sequences for investigation. Its requirement of only sub-microgram quantities of DNA facilitates the examination of large numbers of samples. However, the special features of AluScan data rendered difficult the calling of copy number variation (CNV) directly using the calling algorithms designed for whole genome sequencing (WGS) or exome sequencing. RESULTS: In this study, an AluScanCNV package has been assembled for efficient CNV calling from AluScan sequencing data employing a Geary-Hinkley transformation (GHT) of read-depth ratios between either paired test-control samples, or between test samples and a reference template constructed from reference samples, to call the localized CNVs, followed by use of a GISTIC-like algorithm to identify recurrent CNVs and circular binary segmentation (CBS) to reveal large extended CNVs. To evaluate the utility of CNVs called from AluScan data, the AluScans from 23 non-cancer and 38 cancer genomes were analyzed in this study. The glioma samples analyzed yielded the familiar extended copy-number losses on chromosomes 1p and 9. Also, the recurrent somatic CNVs identified from liver cancer samples were similar to those reported for liver cancer WGS with respect to a striking enrichment of copy-number gains in chromosomes 1q and 8q. When localized or recurrent CNV-features capable of distinguishing between liver and non-liver cancer samples were selected by correlation-based machine learning, a highly accurate separation of the liver and non-liver cancer classes was attained. CONCLUSIONS: The results obtained from non-cancer and cancerous tissues indicated that the AluScanCNV package can be employed to call localized, recurrent and extended CNVs from AluScan sequences. Moreover, both the localized and recurrent CNVs identified by this method could be subjected to machine-learning selection to yield distinguishing CNV-features that were capable of separating between liver cancers and other types of cancers. Since the method is applicable to any human DNA sample with or without the availability of a paired control, it can also be employed to analyze the constitutional CNVs of individuals.

Application of Machine Learning to Development of Copy Number Variation-based Prediction of Cancer Risk.

In the present study, recurrent copy number variations (CNVs) from non-tumor blood cell DNAs of Caucasian non-cancer subjects and glioma, myeloma, and colorectal cancer-patients, and Korean non-cancer subjects and hepatocellular carcinoma, gastric cancer, and colorectal cancer patients, were found to reveal for each of the two ethnic cohorts highly significant differences between cancer patients and controls with respect to the number of CN-losses and size-distribution of CN-gains, suggesting the existence of recurrent constitutional CNV-features useful for prediction of predisposition to cancer. Upon identification by machine learning, such CNV-features could extensively discriminate between cancer-patient and control DNAs. When the CNV-features selected from a learning-group of Caucasian or Korean mixed DNAs consisting of both cancer-patient and control DNAs were employed to make predictions on the cancer predisposition of an unseen test group of mixed DNAs, the average prediction accuracy was 93.6% for the Caucasian cohort and 86.5% for the Korean cohort.

Social cognitive role of schizophrenia candidate gene GABRB2.

The occurrence of positive selection in schizophrenia-associated GABRB2 suggests a broader impact of the gene product on population fitness. The present study considered the possibility of cognition-related GABRB2 involvement by examining the association of GABRB2 with psychosis and altruism, respectively representing psychiatric and psychological facets of social cognition. Four single nucleotide polymorphisms (SNPs) were genotyped for quantitative trait analyses and population-based association studies. Psychosis was measured by either the Positive and Negative Syndrome Scale (PANSS) or antipsychotics dosage, and altruism was based on a self-report altruism scale. The minor alleles of SNPs rs6556547, rs1816071 and rs187269 in GABRB2 were correlated with high PANSS score for positive symptoms in a Han Chinese schizophrenic cohort, whereas those of rs1816071 and rs1816072 were associated with high antipsychotics dosage in a US Caucasian schizophrenic cohort. Moreover, strongly significant GABRB2-disease associations were found among schizophrenics with severe psychosis based on high PANSS positive score, but no significant association was observed for schizophrenics with only mild psychosis. Interestingly, in addition to association with psychosis in schizophrenics, rs187269 was also associated with altruism in healthy Han Chinese. Furthermore, parallel to correlation with severe psychosis, its minor allele was correlated with high altruism scores. These findings revealed that GABRB2 is associated with psychosis, the core symptom and an endophenotype of schizophrenia. Importantly, the association was found across the breadth of the psychiatric (psychosis) to psychological (altruism) spectrum of social cognition suggesting GABRB2 involvement in human cognition.

Epigenetic regulation on GABRB2 isoforms expression: developmental variations and disruptions in psychotic disorders.

INTRODUCTION: To improve the understanding of psychotic abnormalities and their non-Mendelian inheritance patterns, the epigenetic regulation of the psychotic disorder-associated GABRB2, gene for the type A γ-aminobutyric acid receptor ß(2)-subunit, was investigated. METHODS: Expression of GABRB2, and the epigenetic regulatory enzymes histone deacetylases (HDACs) and DNA methyltransferases (DNMTs) in mouse and postmortem human brains was analyzed using real-time PCR. RESULTS: Results showed that expression of GABRB2 isoforms significantly increased over time in both mouse and human, especially for the long splicing isoform. In the brains of non-psychiatric controls (CON), a significant positive correlation of GABRB2 expression with age was observed in individuals with MM genotypes of the single nucleotide polymorphisms (SNPs) rs187269 and rs1816072. This was reversed to a significant negative correlation in schizophrenics (SCZ). A similar reversal was also displayed by bipolar disorder (BPD) patients. In parallel, a significant co-variation of HDAC1 with GABRB2 expression observed in CON remained significant in BPD but not in SCZ; comparably, a significant co-variation of HDAC2 with GABRB2 expression observed in CON became non-significant in both SCZ and BPD. Moreover, co-variations of DNMT1 and DNMT3B with GABRB2, not observable in CON, became significant in BPD. CONCLUSION: These findings demonstrated that GABRB2 expression was under epigenetic regulation that varied with development, genotype and disease status, and these regulatory mechanisms were observably disrupted in SCZ and BPD. This study provided insight into the complex inheritance patterns of psychiatric disorders, and pointed to the involvement of epigenetic dysregulation in the disease process of major psychotic disorders.

AluScan: a method for genome-wide scanning of sequence and structure variations in the human genome.

BACKGROUND: To complement next-generation sequencing technologies, there is a pressing need for efficient pre-sequencing capture methods with reduced costs and DNA requirement. The Alu family of short interspersed nucleotide elements is the most abundant type of transposable elements in the human genome and a recognized source of genome instability. With over one million Alu elements distributed throughout the genome, they are well positioned to facilitate genome-wide sequence amplification and capture of regions likely to harbor genetic variation hotspots of biological relevance. RESULTS: Here we report on the use of inter-Alu PCR with an enhanced range of amplicons in conjunction with next-generation sequencing to generate an Alu-anchored scan, or 'AluScan', of DNA sequences between Alu transposons, where Alu consensus sequence-based 'H-type' PCR primers that elongate outward from the head of an Alu element are combined with 'T-type' primers elongating from the poly-A containing tail to achieve huge amplicon range. To illustrate the method, glioma DNA was compared with white blood cell control DNA of the same patient by means of AluScan. The over 10 Mb sequences obtained, derived from more than 8,000 genes spread over all the chromosomes, revealed a highly reproducible capture of genomic sequences enriched in genic sequences and cancer candidate gene regions. Requiring only sub-micrograms of sample DNA, the power of AluScan as a discovery tool for genetic variations was demonstrated by the identification of 357 instances of loss of heterozygosity, 341 somatic indels, 274 somatic SNVs, and seven potential somatic SNV hotspots between control and glioma DNA. CONCLUSIONS: AluScan, implemented with just a small number of H-type and T-type inter-Alu PCR primers, provides an effective capture of a diversity of genome-wide sequences for analysis. The method, by enabling an examination of gene-enriched regions containing exons, introns, and intergenic sequences with modest capture and sequencing costs, computation workload and DNA sample requirement is particularly well suited for accelerating the discovery of somatic mutations, as well as analysis of disease-predisposing germline polymorphisms, by making possible the comparative genome-wide scanning of DNA sequences from large human cohorts.

Effects of flavone 6-substitutions on GABAA receptors efficacy.

Flavones have been studied for their activities via benzodiazepine site on the type-A γ-aminobutyric acid (GABA(A)) receptors, for which knowledge on structure-efficacy relationships has been rather limited in comparison to that on structure-affinity relationships. The present study focused on flavone 6-substitution, implied in previous studies being relevant to efficacy. Structure analogs, each varying only at position 6, were compared, including 6-fluoroflavone, 6-chloroflavone, 6-bromoflavone, and 2'-hydroxyflavone analyzed in the present study, as well as 6,2'-dihydroxyflavone reported earlier. Radio-ligand binding assays, whole-cell patch-clamp, and mouse behavioral experiments were performed. In consistent with a previous report, the present whole-cell patch-clamp and animal behavior experiments demonstrated 6-bromoflavone to be a positive modulator at GABA(A) receptors acting through flumazenil-sensitive high-affinity benzodiazepine site. In contrast, the other two 6-haloflavones were both neutralizing modulators. In vitro electrophysiological and in vivo animal experiments showed that 2'-hydroxyflavone was a neutralizing modulator, different in efficacy from its structural analog, 6,2'-dihydroxyflavone, a negative modulator of GABA(A) receptors. The fact that flavone analogs differing only at position 6 showed drastically different pharmacological properties clearly points to 6-substitution being an important determinant of efficacy. The results suggest that a large width of the first atom on the 6-substituent favors a high binding affinity of the 6-substituted flavone, whereas a large overall volume of the 6-substituent favors positive modulator activity, which could be modified by, e.g., 2'-hydroxyl substitution. These findings have contributed to the understanding of quantitative structure-efficacy relationships for flavones acting at GABA(A) receptors, and hence facilitation of flavone-based drug development.

A recombination hotspot in a schizophrenia-associated region of GABRB2.

BACKGROUND: Schizophrenia is a major disorder with complex genetic mechanisms. Earlier, population genetic studies revealed the occurrence of strong positive selection in the GABRB2 gene encoding the beta(2) subunit of GABA(A) receptors, within a segment of 3,551 bp harboring twenty-nine single nucleotide polymorphisms (SNPs) and containing schizophrenia-associated SNPs and haplotypes. METHODOLOGY/PRINCIPAL FINDINGS: In the present study, the possible occurrence of recombination in this 'S1-S29' segment was assessed. The occurrence of hotspot recombination was indicated by high resolution recombination rate estimation, haplotype diversity, abundance of rare haplotypes, recurrent mutations and torsos in haplotype networks, and experimental haplotyping of somatic and sperm DNA. The sub-segment distribution of relative recombination strength, measured by the ratio of haplotype diversity (H(d)) over mutation rate (theta), was indicative of a human specific Alu-Yi6 insertion serving as a central recombining sequence facilitating homologous recombination. Local anomalous DNA conformation attributable to the Alu-Yi6 element, as suggested by enhanced DNase I sensitivity and obstruction to DNA sequencing, could be a contributing factor of the increased sequence diversity. Linkage disequilibrium (LD) analysis yielded prominent low LD points that supported ongoing recombination. LD contrast revealed significant dissimilarity between control and schizophrenic cohorts. Among the large array of inferred haplotypes, H26 and H73 were identified to be protective, and H19 and H81 risk-conferring, toward the development of schizophrenia. CONCLUSIONS/SIGNIFICANCE: The co-occurrence of hotspot recombination and positive selection in the S1-S29 segment of GABRB2 has provided a plausible contribution to the molecular genetics mechanisms for schizophrenia. The present findings therefore suggest that genome regions characterized by the co-occurrence of positive selection and hotspot recombination, two interacting factors both affecting genetic diversity, merit close scrutiny with respect to the etiology of common complex disorders.

GABRB2 in schizophrenia and bipolar disorder: disease association, gene expression and clinical correlations.

The SCZ (schizophrenia)-associated GABA(A) receptor (gamma-aminobutyric acid type A receptor) beta(2) subunit gene GABRB2 was recently associated with BPD (bipolar disorder). Although weaker than its association with SCZ, significant association of GABRB2 with BPD was found in both German and Chinese, especially for the haplotypes rs1816071-rs187269 and rs1816072-rs187269 for which the M-M variants showed higher frequency in disease than the control. Significant genotype-dependent reduction in GABRB2 expression was shown for BPD, but to a lesser extent than that for SCZ. Temporal effects on GABRB2 expression were observed. Moreover, for the homozygous major genotypes of rs1816071, rs1816072 and rs187269, expression increased with time in CON but decreased in SCZ and BPD. The genotypes of these three SNPs (single nucleotide polymorphisms) were further correlated with antipsychotics dosage in SCZ cohorts. The findings highlight the importance of GABRB2 in neuropsychiatric disease aetiology, with respect to haplotype association, as well as reduction of and temporal effects on gene expression in both SCZ and BPD, but to a lesser extent in the latter, supporting the suggestion that functional psychosis can be conceptualized as a continuous spectrum of clinical phenotypes rather than as distinct categories.

Alternative-splicing in the exon-10 region of GABA(A) receptor beta(2) subunit gene: relationships between novel isoforms and psychotic disorders.

BACKGROUND: Non-coding single nucleotide polymorphisms (SNPs) in GABRB2, the gene for beta(2)-subunit of gamma-aminobutyric acid type A (GABA(A)) receptor, have been associated with schizophrenia (SCZ) and quantitatively correlated to mRNA expression and alternative splicing. METHODS AND FINDINGS: Expression of the Exon 10 region of GABRB2 from minigene constructs revealed this region to be an "alternative splicing hotspot" that readily gave rise to differently spliced isoforms depending on intron sequences. This led to a search in human brain cDNA libraries, and the discovery of two novel isoforms, beta(2S1) and beta(2S2), bearing variations in the neighborhood of Exon-10. Quantitative real-time PCR analysis of postmortem brain samples showed increased beta(2S1) expression and decreased beta(2S2) expression in both SCZ and bipolar disorder (BPD) compared to controls. Disease-control differences were significantly correlated with SNP rs187269 in BPD males for both beta(2S1) and beta(2S2) expressions, and significantly correlated with SNPs rs2546620 and rs187269 in SCZ males for beta(2S2) expression. Moreover, site-directed mutagenesis indicated that Thr(365), a potential phosphorylation site in Exon-10, played a key role in determining the time profile of the ATP-dependent electrophysiological current run-down. CONCLUSION: This study therefore provided experimental evidence for the importance of non-coding sequences in the Exon-10 region in GABRB2 with respect to beta(2)-subunit splicing diversity and the etiologies of SCZ and BPD.

GABA A receptor subtype selectivity underlying selective anxiolytic effect of baicalin.

Baicalin, a naturally occurring flavonoid, was previously reported to induce anxiolytic-like effect devoid of sedation and myorelaxation in mice, acting through type A gamma-aminobutyric acid (GABA(A)) receptor benzodiazepine (BZ) site. The present study further expanded the behavioral pharmacology profile of baicalin and subtype selectivity was explored as a possible mechanism underlying its in vivo effects on mice. Baicalin was characterized using convulsion, memory, and motor function related animal tests; and its selectivity towards recombinant GABA(A) receptor subtypes expressed in HEK 293T cells was determined by radioligand binding assay and electrophysiological studies. In the picrotoxin-induced seizure, step-through passive avoidance and rotarod tests, the anticonvulsant, amnesic and motor incoordination effects commonly associated with classical BZs were not observed when baicalin was administered at effective anxiolytic doses, demonstrating a separation of the anticonvulsant, amnesic and motor incoordination effects from the anxiolytic-like effect. Although baicalin exhibited higher binding affinity for the alpha1-containing GABA(A) subtype compared with alpha2-, alpha3-, and alpha5-containing subtypes, this was not statistically significant. In contrast to the classical BZ diazepam, baicalin showed significant preference for alpha2- and alpha3-containing subtypes compared to alpha1- and alpha5-containing subtypes in whole-cell patch clamp studies (P < 0.01). Its subtype selectivity suggested that baicalin exerted its in vivo anxiolytic-like effect mainly through the alpha2- and alpha3-containing subtypes. Therefore, the present study revealed an underlying mechanism for the selective anxiolytic profile of baicalin, suggesting alpha2- and alpha3-containing subtypes were important drug targets for flavonoid-based anxiolytics.