Review : alternative splicing (AS) of genes as an approach for generating protein complexity

Prior to the completion of the human genome project, the human genome was thought to have a greater number of genes as it seemed structurally and functionally more complex than other simpler organisms. This along with the belief of “one gene, one protein”, were demonstrated to be incorrect. The inequality in the ratio of gene to protein formation gave rise to the theory of alternative splicing (AS). AS is a mechanism by which one gene gives rise to multiple protein products. Numerous databases and online bioinformatic tools are available for the detection and analysis of AS. Bioinformatics provides an important approach to study mRNA and protein diversity by various tools such as expressed sequence tag (EST) sequences obtained from completely processed mRNA. Microarrays and deep sequencing approaches also aid in the detection of splicing events. Initially it was postulated that AS occurred only in about 5%; of all genes but was later found to be more abundant. Using bioinformatic approaches, the level of AS in human genes was found to be fairly high with 35-59%; of genes having at least one AS form. Our ability to determine and predict AS is important as disorders in splicing patterns may lead to abnormal splice variants resulting in genetic diseases. In addition, the diversity of proteins produced by AS poses a challenge for successful drug discovery and therefore a greater understanding of AS would be beneficial.

Circulating microRNAs involved in multiple sclerosis

Multiple sclerosis (MS) is an immune-mediated, demyelinating and neurodegenerative disease of the central nervous system. After traumatic brain injury, it is the leading cause of neurology disability in young adults. Considerable advances have been made in identifying genes involved in MS but the genetic and phenotypic complexity associated with this disease significantly hinders any progress. A novel class of small RNA molecules, microRNAs (miRNAs) has acquired much attention because they regulate the expression of up to 30% of protein-coding genes and may play a pivotal role in the development of many, if not all, complex diseases. Seven published studies investigated miRNAs from peripheral blood mononuclear cells, CD4+, CD8+ T cell, B lymphocytes, peripheral blood leukocytes, whole blood and brain astrocytes with MS risk. The absence of MS studies investigating plasma miRNA prompted the current investigation of identifying a circulating miRNA signature in MS. We conducted a microarray analysis of over 900 known miRNA transcripts from plasma samples collected from four MS individuals and four sex-aged and ethnicity matched healthy controls. We identified six plasma miRNA (miR-614, miR-572, miR-648, miR-1826, miR-422a and miR-22) that were significantly up-regulated and one plasma miRNA (miR-1979) that was significantly down-regulated in MS individuals. Both miR-422a and miR-22 have previously been implicated in MS. The present study is the first to show a circulating miRNA signature involved in MS that could serve as a potential prognostic and diagnostic biomarker for MS.

Modeling the cumulative genetic risk for multiple sclerosis from genome-wide association data

Background: Multiple sclerosis (MS) is the most common cause of chronic neurologic disability beginning in early to middle adult life. Results from recent genome-wide association studies (GWAS) have substantially lengthened the list of disease loci and provide convincing evidence supporting a multifactorial and polygenic model of inheritance. Nevertheless, the knowledge of MS genetics remains incomplete, with many risk alleles still to be revealed. Methods: We used a discovery GWAS dataset (8,844 samples, 2,124 cases and 6,720 controls) and a multi-step logistic regression protocol to identify novel genetic associations. The emerging genetic profile included 350 independent markers and was used to calculate and estimate the cumulative genetic risk in an independent validation dataset (3,606 samples). Analysis of covariance (ANCOVA) was implemented to compare clinical characteristics of individuals with various degrees of genetic risk. Gene ontology and pathway enrichment analysis was done using the DAVID functional annotation tool, the GO Tree Machine, and the Pathway-Express profiling tool. Results: In the discovery dataset, the median cumulative genetic risk (P-Hat) was 0.903 and 0.007 in the case and control groups, respectively, together with 79.9% classification sensitivity and 95.8% specificity. The identified profile shows a significant enrichment of genes involved in the immune response, cell adhesion, cell communication/ signaling, nervous system development, and neuronal signaling, including ionotropic glutamate receptors, which have been implicated in the pathological mechanism driving neurodegeneration. In the validation dataset, the median cumulative genetic risk was 0.59 and 0.32 in the case and control groups, respectively, with classification sensitivity 62.3% and specificity 75.9%. No differences in disease progression or T2-lesion volumes were observed among four levels of predicted genetic risk groups (high, medium, low, misclassified). On the other hand, a significant difference (F = 2.75, P = 0.04) was detected for age of disease onset between the affected misclassified as controls (mean = 36 years) and the other three groups (high, 33.5 years; medium, 33.4 years; low, 33.1 years). Conclusions: The results are consistent with the polygenic model of inheritance. The cumulative genetic risk established using currently available genome-wide association data provides important insights into disease heterogeneity and completeness of current knowledge in MS genetics.

A polymorphism in the HLA-DPB1 gene is associated with susceptibility to multiple sclerosis

We conducted an association study across the human leukocyte antigen (HLA) complex to identify loci associated with multiple sclerosis (MS). Comparing 1927 SNPs in 1618 MS cases and 3413 controls of European ancestry, we identified seven SNPs that were independently associated with MS conditional on the others (each ). All associations were significant in an independent replication cohort of 2212 cases and 2251 controls () and were highly significant in the combined dataset (). The associated SNPs included proxies for HLA-DRB1*15:01 and HLA-DRB1*03:01, and SNPs in moderate linkage disequilibrium (LD) with HLA-A*02:01, HLA-DRB1*04:01 and HLA-DRB1*13:03. We also found a strong association with rs9277535 in the class II gene HLA-DPB1 (discovery set , replication set , combined ). HLA-DPB1 is located centromeric of the more commonly typed class II genes HLA-DRB1, -DQA1 and -DQB1. It is separated from these genes by a recombination hotspot, and the association is not affected by conditioning on genotypes at DRB1, DQA1 and DQB1. Hence rs9277535 represents an independent MS-susceptibility locus of genome-wide significance. It is correlated with the HLA-DPB1*03:01 allele, which has been implicated previously in MS in smaller studies. Further genotyping in large datasets is required to confirm and resolve this association.

Association analysis of chromosome 1 migraine candidate genes

Migraine with aura (MA) is a subtype of typical migraine. Migraine with aura (MA) also encompasses a rare severe subtype Familial Hemiplegic Migraine (FHM) with several known genetic loci. The type 2 FHM (FHM-2) susceptibility locus maps to chromosome 1q23 and mutations in the ATP1A2 gene at this site have recently been implicated. We have previously provided evidence of linkage of typical migraine (predominantly MA) to microsatellite markers on chromosome 1, in the 1q31 and 1q23 regions. In this study, we have undertaken a large genomic investigation involving candidate genes that lie within the chromosome 1q23 and 1q31 regions using an association analysis approach. Methods We have genotyped a large population of case-controls (243 unrelated Caucasian migraineurs versus 243 controls) examining a set of 5 single nucleotide polymorphisms (SNPs) and the Fas Ligand dinucleotide repeat marker, located within the chromosome 1q23 and 1q31 regions. Results Several genes have been studied including membrane protein (ATP 1 subtype A4 and FasL), cytoplasmic glycoprotein (CASQ 1) genes and potassium (KCN J9 and KCN J10) and calcium (CACNA1E) channel genes in 243 migraineurs (including 85% MA and 15% of migraine without aura (MO)) and 243 matched controls. After correction for multiple testing, chi-square results showed non-significant P values (P > 0.008) across all SNPs (and a CA repeat) tested in these different genes, however results with the KCN J10 marker gave interesting results (P = 0.02) that may be worth exploring further in other populations. Conclusion These results do not show a significant role for the tested candidate gene variants and also do not support the hypothesis that a common chromosome 1 defective gene influences both FHM and the more common forms of migraine.

Gene expression studies in multiple sclerosis

Multiple sclerosis (MS) is a serious neurological disorder affecting young Caucasian individuals, usually with an age of onset at 18 to 40 years old. Females account for approximately 60x of MS cases and the manifestation and course of the disease is highly variable from patient to patient. The disorder is characterised by the development of plaques within the central nervous system (CNS). Many gene expression studies have been undertaken to look at the specific patterns of gene transcript levels in MS. Human tissues and experimental mice were used in these gene-profiling studies and a very valuable and interesting set of data has resulted from these various expression studies. In general, genes showing variable expression include mainly immunological and inflammatory genes, stress and antioxidant genes, as well as metabolic and central nervous system markers. Of particular interest are a number of genes localised to susceptible loci previously shown to be in linkage with MS. However due to the clinical complexity of the disease, the heterogeneity of the tissues used in expression studies, as well as the variable DNA chips/membranes used for the gene profiling, it is difficult to interpret the available information. Although this information is essential for the understanding of the pathogenesis of MS, it is difficult to decipher and define the gene pathways involved in the disorder. Experiments in gene expression profiling in MS have been numerous and lists of candidates are now available for analysis. Researchers have investigated gene expression in peripheral mononuclear white blood cells (PBMCs), in MS animal models Experimental Allergic Encephalomyelitis (EAE) and post mortem MS brain tissues. This review will focus on the results of these studies.

Quantitative and qualitative changes in gene expression patterns characterize the activity of plaques in multiple sclerosis

Multiple sclerosis (MS) is a complex autoimmune disorder of the CNS with both genetic and environmental contributing factors. Clinical symptoms are broadly characterized by initial onset, and progressive debilitating neurological impairment. In this study, RNA from MS chronic active and MS acute lesions was extracted, and compared with patient matched normal white matter by fluorescent cDNA microarray hybridization analysis. This resulted in the identification of 139 genes that were differentially regulated in MS plaque tissue compared to normal tissue. Of these, 69 genes showed a common pattern of expression in the chronic active and acute plaque tissues investigated (Pvalue<0.0001, ρ=0.73, by Spearman's ρ analysis); while 70 transcripts were uniquely differentially expressed (≥1.5-fold) in either acute or chronic active tissues. These results included known markers of MS such as the myelin basic protein (MBP) and glutathione S-transferase (GST) M1, nerve growth factors, such as nerve injury-induced protein 1 (NINJ1), X-ray and excision DNA repair factors (XRCC9 and ERCC5) and X-linked genes such as the ribosomal protein, RPS4X. Primers were then designed for seven array-selected genes, including transferrin (TF), superoxide dismutase 1 (SOD1), glutathione peroxidase 1 (GPX1), GSTP1, crystallin, alpha-B (CRYAB), phosphomannomutase 1 (PMM1) and tubulin β-5 (TBB5), and real time quantitative (Q)-PCR analysis was performed. The results of comparative Q-PCR analysis correlated significantly with those obtained by array analysis (r=0.75, Pvalue<0.01, by Pearson's bivariate correlation). Both chronic active and acute plaques shared the majority of factors identified suggesting that quantitative, rather than gross qualitative differences in gene expression pattern may define the progression from acute to chronic active plaques in MS.

Association and linkage analyses of restriction fragment length polymorphisms for the human renin and antithrombin III genes in essential hypertension

Essential hypertension is a highly hereditable disorder in which genetic influences predominate over environmental factors. The molecular genetic profiles which predispose to essential hypertension are not known. In rats with genetic hypertension, there is some recent evidence pointing to linkage of renin gene alleles with blood pressure. The genes for renin and antithrombin III belong to a conserved synteny group which, in humans, spans the q21.3-32.3 region of chromosome I and, in rats, is linkage group X on chromosome 13. The present study examined the association of particular human renin gene (REN) and antithrombin III gene (AT3) polymorphisms with essential hypertension by comparing the frequency of specific alleles for each of these genes in 50 hypertensive offspring of hypertensive parents and 91 normotensive offspring of normotensive parents. In addition, linkage relationships were examined in hypertensive pedigrees with multiple affected individuals. Alleles of a REN HindIII restriction fragment length polymorphism (RFLP) were detected using a genomic clone, λHR5, to probe Southern blots of HindIII-cut leucocyte DNA, and those for an AT3 Pstl RFLP were detected by phATIII 113 complementary DNA probe. The frequencies of each REN allele in the hypertensive group were 0.76 and 0.24 compared with 0.74 and 0.26 in the normotensive group. For AT3, hypertensive allele frequencies were 0.49 and 0.51 compared with normotensive values of 0.54 and 0.46. These differences were not significant by χ2 analysis (P > 0.2). Linkage analysis of a family (data from 16 family members, 10 of whom were hypertensive), informative for both markers, without an age-of-onset correction, and assuming dominant inheritance of hypertension, complete penetrance and a disease frequency of 20%, did not indicate linkage of REN with hypertension, but gave a positive, although not significant, logarithm of the odds for linkage score of 0.784 at a recombination fraction of 0 for AT3 linkage to hypertension. In conclusion, the present study could find no evidence for an association of a REN HindIII RFLP with essential hypertension or for a linkage of the locus defined by this RFLP in a family segregating for hypertension. In the case of an AT3 Pstl RFLP, although association analysis was negative, linkage analysis suggested possible involvement (odds of 6:1 in favour) of a gene located near the 1q23 locus with hypertension in one informative family.

A transcription factor map as revealed by a genome-wide gene expression analysis of whole-blood mRNA transcriptome in multiple sclerosis

Background Several lines of evidence suggests that transcription factors are involved in the pathogenesis of Multiple Sclerosis (MS) but a complete mapping the whole network has been elusive. One of the reasons is that there are several clinical subtypes of MS and transcription factors which may be involved in one subtype may not be in others. We investigated the possibility that this network could be mapped using microarray technologies and modern bioinformatics methods on a dataset from whole blood in 99 untreated MS patients (36 Relapse Remitting MS, 43 Primary Progressive MS, and 20 Secondary Progressive MS) and 45 age-matched healthy controls, Methodology/Principal Findings We have used two different analytical methodologies: a differential expression analysis and a differential co-expression analysis, which have converged on a significant number of regulatory motifs that seem to be statistically overrepresented in genes which are either differentially expressed (or differentially co-expressed) in cases and controls (e.g. V$KROX_Q6, p-value < 3.31E-6; V$CREBP1_Q2, p-value < 9.93E-6, V$YY1_02, p-value < 1.65E-5). Conclusions/significance: Our analysis uncovered a network of transcription factors that potentially dysregulate several genes in MS or one or more of its disease subtypes. Analysing the published literature we have found that these transcription factors are involved in the early T-lymphocyte specification and commitment as well as in oligodendrocytes dedifferentiation and development. The most significant transcription factors motifs were for the Early Growth response EGR/KROX family, ATF2, YY1 (Yin and Yang 1), E2F-1/DP-1 and E2F-4/DP-2 heterodimers, SOX5, and CREB and ATF families.

Progesterone activates multiple innate immune pathways inChlamydia trachomatis-Infected endocervical cells

Problem Susceptibility to Chlamydia trachomatis infection is increased by oral con- traceptives and modulated by sex hormones. We therefore sought to determine the effects of female sex hormones on the innate immune response to C. trachomatis infection. Method of study ECC-1 endometrial cells, pre-treated with oestradiol or progesterone, were infected with C. trachomatis and the host transcriptome analysed by Illumina Sentrix HumanRef-8 microarray. Primary endocervical epithe- lial cells, prepared at either the proliferative or secretory phase of the menstrual cycle, were infected with C. trachomatis and cytokine gene expression determined by quantitative RT-PCR analysis. Results Chlamydia trachomatis yield from progesterone-primed ECC-1 cells was significantly reduced compared with oestradiol-treated cells. Genes upregulated in progesterone-treated and Chlamydia-infected cells only included multiple CC and CXC chemokines, IL-17C, IL-29, IL-32, TNF-a, DEFB4B, LCN2, S100A7-9, ITGAM, NOD2, JAK1, IL-6ST, type I and II interferon receptors, numerous interferon-stimulated genes and STAT6. CXCL10, CXCL11, CX3CL1 and IL-17C, which were also upregu- lated in infected secretory-stage primary cells, and there was a trend towards higher levels of immune mediators in infected secretory-phase compared with proliferative-phase cells. Conclusion Progesterone treatment primes multiple innate immune pathways in hormone-responsive epithelial cells that could potentially increase resis- tance to chlamydial infection.

Multiple repeats of a promoter segment causes transcription factor autoregulation in red apples

Mutations in the genes encoding for either the biosynthetic or transcriptional regulation of the anthocyanin pathway have been linked to color phenotypes. Generally, this is a loss of function resulting in a reduction or a change in the distribution of anthocyanin. Here, we describe a rearrangement in the upstream regulatory region of the gene encoding an apple (Malus x domestica) anthocyanin-regulating transcription factor, MYB10. We show that this modification is responsible for increasing the level of anthocyanin throughout the plant to produce a striking phenotype that includes red foliage and red fruit flesh. This rearrangement is a series of multiple repeats, forming a minisatellite-like structure that comprises five direct tandem repeats of a 23-bp sequence. This MYB10 rearrangement is present in all the red foliage apple varieties and species tested but in none of the white fleshed varieties. Transient assays demonstrated that the 23-bp sequence motif is a target of the MYB10 protein itself, and the number of repeat units correlates with an increase in transactivation by MYB10 protein. We show that the repeat motif is capable of binding MYB10 protein in electrophoretic mobility shift assays. Taken together, these results indicate that an allelic rearrangement in the promoter of MYB10 has generated an autoregulatory locus, and this autoregulation is sufficient to account for the increase in MYB10 transcript levels and subsequent ectopic accumulation of anthocyanins throughout the plant.

Mapping eQTLs in the Norfolk Island genetic isolate identifies candidate genes for CVD risk traits

Cardiovascular disease (CVD) affects millions of people worldwide and is influenced by numerous factors, including lifestyle and genetics. Expression quantitative trait loci (eQTLs) influence gene expression and are good candidates for CVD risk. Founder-effect pedigrees can provide additional power to map genes associated with disease risk. Therefore, we identified eQTLs in the genetic isolate of Norfolk Island (NI) and tested for associations between these and CVD risk factors. We measured genome-wide transcript levels of blood lymphocytes in 330 individuals and used pedigree-based heritability analysis to identify heritable transcripts. eQTLs were identified by genome-wide association testing of these transcripts. Testing for association between CVD risk factors (i.e., blood lipids, blood pressure, and body fat indices) and eQTLs revealed 1,712 heritable transcripts (p < 0.05) with heritability values ranging from 0.18 to 0.84. From these, we identified 200 cis-acting and 70 trans-acting eQTLs (p < 1.84 × 10(-7)) An eQTL-centric analysis of CVD risk traits revealed multiple associations, including 12 previously associated with CVD-related traits. Trait versus eQTL regression modeling identified four CVD risk candidates (NAAA, PAPSS1, NME1, and PRDX1), all of which have known biological roles in disease. In addition, we implicated several genes previously associated with CVD risk traits, including MTHFR and FN3KRP. We have successfully identified a panel of eQTLs in the NI pedigree and used this to implicate several genes in CVD risk. Future studies are required for further assessing the functional importance of these eQTLs and whether the findings here also relate to outbred populations.

Genome-wide analysis in a murine Dnmt1 knockdown model identifies epigenetically silenced genes in primary human pituitary tumors

DNA methylation at promoter CpG islands (CGI) is an epigenetic modification associated with inappropriate gene silencing in multiple tumor types. In the absence of a human pituitary tumor cell line, small interfering RNA-mediated knockdown of the maintenance methyltransferase DNA methyltransferase (cytosine 5)-1 (Dnmt1) was used in the murine pituitary adenoma cell line AtT-20. Sustained knockdown induced reexpression of the fully methylated and normally imprinted gene neuronatin (Nnat) in a time-dependent manner. Combined bisulfite restriction analysis (COBRA) revealed that reexpression of Nnat was associated with partial CGI demethylation, which was also observed at the H19 differentially methylated region. Subsequent genome-wide microarray analysis identified 91 genes that were significantly differentially expressed in Dnmt1 knockdown cells (10% false discovery rate). The analysis showed that genes associated with the induction of apoptosis, signal transduction, and developmental processes were significantly overrepresented in this list (P < 0.05). Following validation by reverse transcription-PCR and detection of inappropriate CGI methylation by COBRA, four genes (ICAM1, NNAT, RUNX1, and S100A10) were analyzed in primary human pituitary tumors, each displaying significantly reduced mRNA levels relative to normal pituitary (P < 0.05). For two of these genes, NNAT and S100A10, decreased expression was associated with increased promoter CGI methylation. Induced expression of Nnat in stable transfected AtT-20 cells inhibited cell proliferation. To our knowledge, this is the first report of array-based "epigenetic unmasking" in combination with Dnmt1 knockdown and reveals the potential of this strategy toward identifying genes silenced by epigenetic mechanisms across species boundaries.

Signaling pathway genes for blood pressure, folate and cholesterol levels among hypertensives : an epistasis analysis

Irregular atrial pressure, defective folate and cholesterol metabolism contribute to the pathogenesis of hypertension. However, little is known about the combined roles of the methylenetetrahydrofolate reductase (MTHFR), apolipoprotein-E (ApoE) and angiotensin-converting enzyme (ACE) genes, which are involved in metabolism and homeostasis. The objective of this study is to investigate the association of the MTHFR 677 C>T and 1298A>C, ACE insertion–deletion (I/D) and ApoE genetic polymorphisms with hypertension and to further explore the epistasis interactions that are involved in these mechanisms. A total of 594 subjects, including 348 normotensive and 246 hypertensive ischemic stroke subjects were recruited. The MTHFR 677 C>T and 1298A>C, ACE I/D and ApoEpolymorphisms were genotyped and the epistasis interaction were analyzed. The MTHFR 677 C>T and ApoE polymorphisms demonstrated significant associations with susceptibility to hypertension in multiple logistic regression models, multifactor dimensionality reduction and a classification and regression tree. In addition, the logistic regression model demonstrated that significant interactions between the ApoE E3E3, E2E4, E2E2 and MTHFR 677 C>T polymorphisms existed. In conclusion, the results of this epistasis study indicated significant association between the ApoE and MTHFR polymorphisms and hypertension.

Association between single-nucleotide polymorphisms in growth factor genes and quality of life in men with prostate cancer and the general population

Purpose Improved survival for men with prostate cancer has led to increased attention to factors influencing quality of life (QOL). As protein levels of vascular endothelial growth factor (VEGF) and insulin-like growth factor 1 (IGF-1) have been reported to be associated with QOL in people with cancer, we sought to identify whether single-nucleotide polymorphisms (SNPs) of these genes were associated with QOL in men with prostate cancer. Methods Multiple linear regression of two data sets (including approximately 750 men newly diagnosed with prostate cancer and 550 men from the general population) was used to investigate SNPs of VEGF and IGF-1 (10 SNPs in total) for associations with QOL (measured by the SF-36v2 health survey). Results Men with prostate cancer who carried the minor ‘T’ allele for IGF-1 SNP rs35767 had higher mean Role-Physical scale scores (≥0.3 SD) compared to non-carriers (p < 0.05). While this association was not identified in men from the general population, one IGF-1 SNP rs7965399 was associated with higher mean Bodily Pain scale scores in men from the general population that was not found in men with prostate cancer. Men from the general population who carried the rare ‘C’ allele had higher mean Bodily Pain scale scores (≥0.3 SD) than non-carriers (p < 0.05). Conclusions Through identifying SNPs that are associated with QOL in men with prostate cancer and men from the general population, this study adds to the mapping of complex interrelationships that influence QOL and suggests a role for IGF-I in physical QOL outcomes. Future research may identify biomarkers associated with increased risk of poor QOL that could assist in the provision of pre-emptive support for those identified at risk.