Gene-gene interaction and gene-pathway analysis of genome-wide association for schizophrenia

Schizophrenia (SZ) is a complex disorder with high heritability and variable phenotypes that has limited success in finding causal genes associated with the disease development. Pathway-based analysis is an effective approach in investigating the molecular mechanism of susceptible genes associated with complex diseases. The etiology of complex diseases could be a network of genetic factors and within the genes, interaction may occur. In this work we argue that some genes might be of small effect that by itself are neither sufficient nor necessary to cause the disease however, their effect may induce slight changes to the gene expression or affect the protein function, therefore, analyzing the gene-gene interaction mechanism within the disease pathway would play crucial role in dissecting the genetic architecture of complex diseases, making the pathway-based analysis a complementary approach to GWAS technique. ^ In this study, we implemented three novel linkage disequilibrium based statistics, the linear combination, the quadratic, and the decorrelation test statistics, to investigate the interaction between linked and unlinked genes in two independent case-control GWAS datasets for SZ including participants of European (EA) and African (AA) ancestries. The EA population included 1,173 cases and 1,378 controls with 729,454 genotyped SNPs, while the AA population included 219 cases and 288 controls with 845,814 genotyped SNPs. We identified 17,186 interacting gene-sets at significant level in EA dataset, and 12,691 gene-sets in AA dataset using the gene-gene interaction method. We also identified 18,846 genes in EA dataset and 19,431 genes in AA dataset that were in the disease pathways. However, few genes were reported of significant association to SZ. ^ Our research determined the pathways characteristics for schizophrenia through the gene-gene interaction and gene-pathway based approaches. Our findings suggest insightful inferences of our methods in studying the molecular mechanisms of common complex diseases.^

Patenting genes and gene sequences: Analysis of public health implications and the law

At issue is whether or not isolated DNA is patent eligible under the U.S. Patent Law and the implications of that determination on public health. The U.S. Patent and Trademark Office has issued patents on DNA since the 1980s, and scientists and researchers have proceeded under that milieu since that time. Today, genetic research and testing related to the human breast cancer genes BRCA1 and BRCA2 is conducted within the framework of seven patents that were issued to Myriad Genetics and the University of Utah Research Foundation between 1997 and 2000. In 2009, suit was filed on behalf of multiple researchers, professional associations and others to invalidate fifteen of the claims underlying those patents. The Court of Appeals for the Federal Circuit, which hears patent cases, has invalidated claims for analyzing and comparing isolated DNA but has upheld claims to isolated DNA. The specific issue of whether isolated DNA is patent eligible is now before the Supreme Court, which is expected to decide the case by year's end. In this work, a systematic review was performed to determine the effects of DNA patents on various stakeholders and, ultimately, on public health; and to provide a legal analysis of the patent eligibility of isolated DNA and the likely outcome of the Supreme Court's decision. ^ A literature review was conducted to: first, identify principle stakeholders with an interest in patent eligibility of the isolated DNA sequences BRCA1 and BRCA2; and second, determine the effect of the case on those stakeholders. Published reports that addressed gene patents, the Myriad litigation, and implications of gene patents on stakeholders were included. Next, an in-depth legal analysis of the patent eligibility of isolated DNA and methods for analyzing it was performed pursuant to accepted methods of legal research and analysis based on legal briefs, federal law and jurisprudence, scholarly works and standard practice legal analysis. ^ Biotechnology, biomedical and clinical research, access to health care, and personalized medicine were identified as the principle stakeholders and interests herein. Many experts believe that the patent eligibility of isolated DNA will not greatly affect the biotechnology industry insofar as genetic testing is concerned; unlike for therapeutics, genetic testing does not require tremendous resources or lead time. The actual impact on biomedical researchers is uncertain, with greater impact expected for researchers whose work is intended for commercial purposes (versus basic science). The impact on access to health care has been surprisingly difficult to assess; while invalidating gene patents might be expected to decrease the cost of genetic testing and improve access to more laboratories and physicians' offices that provide the test, a 2010 study on the actual impact was inconclusive. As for personalized medicine, many experts believe that the availability of personalized medicine is ultimately a public policy issue for Congress, not the courts. ^ Based on the legal analysis performed in this work, this writer believes the Supreme Court is likely to invalidate patents on isolated DNA whose sequences are found in nature, because these gene sequences are a basic tool of scientific and technologic work and patents on isolated DNA would unduly inhibit their future use. Patents on complementary DNA (cDNA) are expected to stand, however, based on the human intervention required to craft cDNA and the product's distinction from the DNA found in nature. ^ In the end, the solution as to how to address gene patents may lie not in jurisprudence but in a fundamental change in business practices to provide expanded licenses to better address the interests of the several stakeholders. ^

MedVir: 3D visual interface applied to gene profile analysis

The use of data mining techniques for the gene profile discovery of diseases, such as cancer, is becoming usual in many researches. These techniques do not usually analyze the relationships between genes in depth, depending on the different variety of manifestations of the disease (related to patients). This kind of analysis takes a considerable amount of time and is not always the focus of the research. However, it is crucial in order to generate personalized treatments to fight the disease. Thus, this research focuses on finding a mechanism for gene profile analysis to be used by the medical and biologist experts. Results: In this research, the MedVir framework is proposed. It is an intuitive mechanism based on the visualization of medical data such as gene profiles, patients, clinical data, etc. MedVir, which is based on an Evolutionary Optimization technique, is a Dimensionality Reduction (DR) approach that presents the data in a three dimensional space. Furthermore, thanks to Virtual Reality technology, MedVir allows the expert to interact with the data in order to tailor it to the experience and knowledge of the expert.

The TIGR Gene Indices: analysis of gene transcript sequences in highly sampled eukaryotic species

While genome sequencing projects are advancing rapidly, EST sequencing and analysis remains a primary research tool for the identification and categorization of gene sequences in a wide variety of species and an important resource for annotation of genomic sequence. The TIGR Gene Indices (http://www.tigr.org/tdb/tgi.shtml) are a collection of species-specific databases that use a highly refined protocol to analyze EST sequences in an attempt to identify the genes represented by that data and to provide additional information regarding those genes. Gene Indices are constructed by first clustering, then assembling EST and annotated gene sequences from GenBank for the targeted species. This process produces a set of unique, high-fidelity virtual transcripts, or Tentative Consensus (TC) sequences. The TC sequences can be used to provide putative genes with functional annotation, to link the transcripts to mapping and genomic sequence data, to provide links between orthologous and paralogous genes and as a resource for comparative sequence analysis.

Evidence from disruption of the lmcpb gene array of Leishmania mexicana that cysteine proteinases are virulence factors.

The mammalian form of the protozoan parasite Leishmania mexicana contains high activity of a cysteine proteinase (LmCPb) encoded on a tandem array of 19 genes (lmcpb). Homozygous null mutants for lmcpb have been produced by targeted gene disruption. All life-cycle stages of the mutant can be cultured in vitro, demonstrating that the gene is not essential for growth or differentiation of the parasite. However, the mutant exhibits a marked phenotype affecting virulence-- its infectivity to macrophages is reduced by 80%. The mutants are as efficient as wild-type parasites in invading macrophages but they only survive in a small proportion of the cells. However, those parasites that successfully infect these macrophages grow normally. Despite their reduced virulence, the mutants are still able to produce subcutaneous lesions in mice, albeit at a slower rate than wild-type parasites. The product of a single copy of lmcpb re-expressed in the null mutant was enzymatically active and restored infectivity toward macrophages to wild-type levels. Double null mutants created for lmcpb and lmcpa (another cathepsin L-like cysteine proteinase) have a similar phenotype to the lmcpb null mutant, showing that LmCPa does not compensate for the loss of LmCPb.

Human plectin: organization of the gene, sequence analysis, and chromosome localization (8q24).

Plectin, a 500-kDa intermediate filament binding protein, has been proposed to provide mechanical strength to cells and tissues by acting as a cross-linking element of the cytoskeleton. To set the basis for future studies on gene regulation, tissue-specific expression, and pathological conditions involving this protein, we have cloned the human plectin gene, determined its coding sequence, and established its genomic organization. The coding sequence contains 32 exons that extend over 32 kb of the human genome. Most of the introns reside within a region encoding the globular N-terminal domain of the molecule, whereas the entire central rod domain and the entire C-terminal globular domain were found to be encoded by single exons of remarkable length, >3 kb and >6 kb, respectively. Overall, the organization of the human plectin gene was strikingly similar to that of human bullous pemphigoid antigen 1 (BPAG1), confirming that both proteins belong to the same gene family. Comparison of the deduced protein sequences for human and rat plectin revealed that they were 93% identical. By using fluorescence in situ hybridization, we have mapped the plectin gene to the long arm of chromosome 8 within the telomeric region. This gene locus (8q24) has previously been implicated in the human blistering skin disease epidermolysis bullosa simplex Ogna. Detailed knowledge of the structure of the plectin gene and its chromosome localization will aid in the elucidation of whether this or any other pathological conditions are linked to alterations in the plectin gene.

Using the R package Shiny to create web applications that facilitate quantitative gene expression analysis

Thesis (Master's)--University of Washington, 2016-06

Comparative gene expression analysis of NKT cell subpopulations

Natural killer T (NKT) cells are a lymphocyte lineage, which has diverse immune regulatory activities in many disease settings. Most previous studies have investigated the functions of this family of cells as a single entity, but more recent evidence highlights the distinct functional and phenotypic properties of NKT cell subpopulations. It is likely that the diverse functions of NKT cells are regulated and coordinated by these different NKT subsets. Little is known about how NKT subsets differ in their interactions with the host. We have undertaken the first microarray analysis comparing the gene expression profiles of activated human NKT cell subpopulations, including CD8(+) NKT cells, which have often been overlooked. We describe the significant gene expression differences among NKT cell subpopulations and some of the molecules likely to confer their distinct functional roles. Several genes not associated previously with NKT cells were shown to be expressed differentially in specific NKT cell subpopulations. Our findings provide new insights into the NKT cell family, which may direct further research toward better manipulation of NKT cells for therapeutic applications.

Genetic variation of temperature-regulated curd induction in cauliflower: elucidation of floral transition by genome-wide association mapping and gene expression analysis

Cauliflower (Brassica oleracea var. botrytis) is a vernalization-responsive crop. High ambient temperatures delay harvest time. The elucidation of the genetic regulation of floral transition is highly interesting for a precise harvest scheduling and to ensure stable market supply. This study aims at genetic dissection of temperature-dependent curd induction in cauliflower by genome-wide association studies and gene expression analysis. To assess temperature dependent curd induction, two greenhouse trials under distinct temperature regimes were conducted on a diversity panel consisting of 111 cauliflower commercial parent lines, genotyped with 14,385 SNPs. Broad phenotypic variation and high heritability (0.93) were observed for temperature-related curd induction within the cauliflower population. GWA mapping identified a total of 18 QTL localized on chromosomes O1, O2, O3, O4, O6, O8, and O9 for curding time under two distinct temperature regimes. Among those, several QTL are localized within regions of promising candidate flowering genes. Inferring population structure and genetic relatedness among the diversity set assigned three main genetic clusters. Linkage disequilibrium (LD) patterns estimated global LD extent of r(2) = 0.06 and a maximum physical distance of 400 kb for genetic linkage. Transcriptional profiling of flowering genes FLOWERING LOCUS C (BoFLC) and VERNALIZATION 2 (BoVRN2) was performed, showing increased expression levels of BoVRN2 in genotypes with faster curding. However, functional relevance of BoVRN2 and BoFLC2 could not consistently be supported, which probably suggests to act facultative and/or might evidence for BoVRN2/BoFLC-independent mechanisms in temperature regulated floral transition in cauliflower. Genetic insights in temperature-regulated curd induction can underpin genetically informed phenology models and benefit molecular breeding strategies toward the development of thermo-tolerant cultivars.

S100A1 deficiency results in prolonged ventricular repolarization in response to sympathetic activation.

S100A1 is a Ca(2+)-binding protein and predominantly expressed in the heart. We have generated a mouse line of S100A1 deficiency by gene trap mutagenesis to investigate the impact of S100A1 ablation on heart function. Electrocardiogram recordings revealed that after beta-adrenergic stimulation S100A1-deficient mice had prolonged QT, QTc and ST intervals and intraventricular conduction disturbances reminiscent of 2 : 1 bundle branch block. In order to identify genes affected by the loss of S100A1, we profiled the mutant and wild type cardiac transcriptomes by gene array analysis. The expression of several genes functioning to the electrical activity of the heart were found to be significantly altered. Although the default prediction would be that mRNA and protein levels are highly correlated, comprehensive immunoblot analyses of salient up- or down-regulated candidate genes of any cellular network revealed no significant changes on protein level. Taken together, we found that S100A1 deficiency results in cardiac repolarization delay and alternating ventricular conduction defects in response to sympathetic activation accompanied by a significantly different transcriptional regulation.

Platelets release novel thiol isomerase enzymes which are recruited to the cell surface following activation

The thiol isomerase enzymes protein disulphide isomerase (PDI) and endoplasmic reticulum protein 5 (ERp5) are released by resting and activated platelets. These re-associate with the cell surface where they modulate a range of platelet responses including adhesion, secretion and aggregation. Recent studies suggest the existence of yet uncharacterised platelet thiol isomerase proteins. This study aimed to identify which other thiol isomerase enzymes are present in human platelets. Through the use of immunoblotting, flow cytometry, cell-surface biotinylation and gene array analysis, we report the presence of five additional thiol isomerases in human and mouse platelets and megakaryocytes, namely; ERp57, ERp72, ERp44, ERp29 and TMX3. ERp72, ERp57, ERp44 and ERp29 are released by platelets and relocate to the cell surface following platelet activation. The transmembrane thiol isomerase TMX3 was also detected on the platelet surface but does not increase following activation. Extracellular PDI is also implicated in the regulation of coagulation by the modulation of tissue factor activity. ERp57 was identified within platelet-derived microparticle fractions, suggesting that ERp57 may also be involved in the regulation of coagulation as well as platelet function. These data collectively implicate the expanding family of platelet-surface thiol isomerases in the regulation of haemostasis.

Basement membrane remodeling in skeletal muscles of patients with limb ischemia involves regulation of matrix metalloproteinases and tissue inhibitor of matrix metalloproteinases

BACKGROUND/AIM: Because the pericapillary basement membrane in skeletal muscles of patients with chronic critical limb ischemia (CLI) is thickened, we determined the expression patterns of genes involved in collagen metabolism, using samples from 9 CLI patients, 4 patients with acute limb ischemia and 4 healthy controls. METHODS: Gene array analysis, quantitative RT-PCR and semiquantitative grading of immunohistochemical reactivity were performed to determine mRNA/cDNA and protein concentrations. RESULTS: In CLI patients compared to controls, cDNA levels of matrix metalloproteinase (MMP)-9 and MMP-19 were higher, collagen type IV chains A1 and A2, tissue inhibitor of matrix metalloproteinase (TIMP)-1 and TIMP-2 were similar and MMP-2 were lower. On the protein level, MMP-2, MMP-9, MMP-19 and TIMP-1 were more abundantly expressed. In skeletal muscles from patients with acute limb ischemia, cDNA and protein levels of MMP-9, MMP-19, collagen type IV chains, TIMP-1 and TIMP-2 were high. MMP-2 was elevated at the protein but decreased on the cDNA level. CONCLUSION: Expression of basement membrane components in skeletal muscles of CLI and acute limb ischemia patients is altered, possibly contributing to the pathogenesis of peripheral arterial disease.