Gene structure-based splice variant deconvolution using a microarray platform

MOTIVATION: Alternative splicing allows a single gene to generate multiple mRNAs, which can be translated into functionally and structurally diverse proteins. One gene can have multiple variants coexisting at different concentrations. Estimating the relative abundance of each variant is important for the study of underlying biological function. Microarrays are standard tools that measure gene expression. But most design and analysis has not accounted for splice variants. Thus splice variant-specific chip designs and analysis algorithms are needed for accurate gene expression profiling. RESULTS: Inspired by Li and Wong (2001), we developed a gene structure-based algorithm to determine the relative abundance of known splice variants. Probe intensities are modeled across multiple experiments using gene structures as constraints. Model parameters are obtained through a maximum likelihood estimation (MLE) process/framework. The algorithm produces the relative concentration of each variant, as well as an affinity term associated with each probe. Validation of the algorithm is performed by a set of controlled spike experiments as well as endogenous tissue samples using a human splice variant array.     

Genome-wide Epigenetic Data Facilitate Understanding of Disease Susceptibility Association Studies

Complex traits such as susceptibility to diseases are determined in part by variants at multiple genetic loci. Genome-wide association studies can identify these loci, but most phenotype-associated variants lie distal to protein-coding regions and are likely involved in regulating gene expression. Understanding how these genetic variants affect complex traits depends on the ability to predict and test the function of the genomic elements harboring them. Community efforts such as the ENCODE Project provide a wealth of data about epigenetic features associated with gene regulation. These data enable the prediction of testable functions for many phenotype-associated variants.     

Two-dimensional RFLP analyses reveal megabase-sized clusters of rRNA gene variants in Arabidopsis thaliana, suggesting local spreading of variants as the mode for gene homogenization during concerted evolution

Eukaryotic genes encoding the precursor of 18S, 5.8S and 25S ribosomal RNA (rRNA genes or rDNA) are virtually identical within a species, yet they evolve rapidly between species, a phenomenon known as concerted evolution. The mechanisms by which sequence homogenization and fixation of new rRNA gene variants occurs within a genome are not clear. In diploid Arabidopsis thaliana , approximately 1500 rRNA genes are tandemly arrayed at two nucleolus organizer regions, one on chromosome 2 ( NOR2 ), the other on chromosome 4 ( NOR4 ). This paper shows that NOR2 and NOR4 are similar in size, each spanning approximately 3.5–4.0 Mbp. Using two-dimensional mapping techniques involving a combination of pulsed-field and conventional gel electrophoresis, the distributions of four distinct rRNA gene variants at NOR2 and NOR4 have been determined. rRNA genes at NOR4 are homogeneous with respect to a Hin dIII site occurring once per gene. In contrast, fewer than 10% of the rRNA genes at NOR2 are Hin dIII-bearing variants. A single intergenic spacer length is found among rRNA genes at NOR2 but three classes of spacer length variants are present at NOR4 . The NOR4 variants are not intermingled with one another; instead, they are highly clustered over distances as large as 1.5 Mbp. These data suggest that in the concerted evolution of rRNA genes, homogenization is a consequence of local spreading of new rRNA gene variants.     

Mutation analysis of oxisterol-binding-protein gene in patients with age-related macular degeneration

Allelic variants of several genes are increasingly recognized as susceptibility factors in age-related macular degeneration (AMD). Because of its metabolic characteristics the macula is sensitive to oxidative damage, and supplementation with antioxidants has been shown to be effective in slowing the progression of disease in AMD patients. The oxisterol-binding-protein (OSBP2) gene is expressed mainly in the retinal pigmented epithelium underlying the macular region. Its product specifically binds and transports oxisterols, the cytotoxic effects of which may be involved in macular damage. The aim of this study was to search for allelic variants of OSBP2 gene, as well as to evaluate several risk factors in 24 patients with AMD; 17 with nonexudative (NE) and 7 with neovascular (NV) form. Total cholesterol was elevated in 66% of the patients, high-density lipoprotein (HDL) cholesterol was reduced in 12%; vitamin A or vitamin E deficiency was not observed. OSBP2 gene analysis was performed in AMD patients and in 110 control subjects by single-stranded conformational polymorphism (SSCP) analysis followed by direct sequencing. Six allelic variants were detected: 2 nonpolymorphic unique exonic variants in 2 AMD subjects and 4 polymorphic variants (2 exonic and 2 intronic). These data indicate a possible role of OSBP2 gene in the pathogenesis of oxidative damage to the macula induced by oxysterols in AMD patients.     

Somatic variation of H-2Kk expression and structure in a T-cell lymphoma: instability, stabilization, high production and structural mutation

In the heterozygous T lymphoma line LDHB, variants which have lost the expression of individual H-2 class I genes are spontaneously generated in vitro at a frequency of 10(-1)-10(-2). A cell line (HK13) in which the class I gene Kk is stably expressed (frequency of loss variants less than 10(-4) was selected from LDHB cells by fluorescence activated cell sorting. Further selection of HK13 cells for high Kk expression led to the isolation of the HK22 line which expresses twice as much Kk as HK13. From HK13 and HK22 cells, spontaneous structural variants of Kk having lost individual serological determinants of the Kk wild-type molecule, were isolated by fluorescence activated cell sorting. Such variants occur at a frequency of 10(-6)-10(-7) per cell per generation. The analysis of these variants indicates that they carry mutations in the Kk structural gene and that HK13 cells express a single Kk gene which appears to be duplicated in HK22 cells. We did not find evidence for the generation of variants expressing 'alien' class I products in the LDHB cell line. The instability of class I gene expression in LDHB cells and the transition to stable expression may represent steps of T-cell differentiation in the thymus.     

Selection and characterization of cells resistant to diphtheria toxin and pseudomonas exotoxin A: presumptive translational mutants.

Two classes of diphtheria toxin-resistant variants were selected from Chinese hamster ovary (CHO-K1) cells: permeability variants, in which uptake of toxin was impaired, and a new class of cytoplasmic variants, which were cross-resistant to Pseudomonas exotoxin. EF-2 prepared from the cytoplasmic variants was resistant to ADP-ribosylation by either toxin. The evidence presented suggests that these are translational variants possessing a mutationally altered EF-2 gene product. These studies also confirmed that Pseudomonas toxin ADP-ribosylates EF-2 in toxin-sensitive intact cells, as well as in cell-free systems.     

Gene Flow,Recombination, and Selection in Cyanobacteria: Population Structure of Geographically Related Planktothrix Freshwater Strains

Several Planktothrix strains, each producing a distinct oligopeptide profile, have been shown to coexist within Lake Steinsfjorden (Norway). Using nonribosomal peptide synthetase (NRPS) genes as markers, it has been shown that the Planktothrix community comprises distinct genetic variants displaying differences in bloom dynamics, suggesting a Planktothrix subpopulation structure. Here, we investigate the Planktothrix variants inhabiting four lakes in southeast of Norway utilizing both NRPS and non-NRPS genes. Phylogenetic analyses showed similar topologies for both NRPS and non-NRPS genes, and the lakes appear to have similar structuring of Planktothrix genetic variants. The structure of distinct variants was also supported by very low genetic diversity within variants compared to the between-variant diversity. Incongruent topologies and split decomposition revealed recombination events between Planktothrix variants. In several strains the gene variants seem to be a result of recombination. Both NRPS and non-NRPS genes are dominated by purifying selection; however, sites subjected to positive selection were also detected. The presence of similar and well-separated Planktothrix variants with low internal genetic diversity indicates gene flow within Planktothrix populations. Further, the low genetic diversity found between lakes (similar range as within lakes) indicates gene flow also between Planktothrix populations and suggests recent, or recurrent, dispersals. Our data also indicate that recombination has resulted in new genetic variants. Stability within variants and the development of new variants are likely to be influenced by selection patterns and within-variant homologous recombination.     

Predicting phenotypic severity of uncertain gene variants in the RET proto-oncogene

Although reported gene variants in the RET oncogene have been directly associated with multiple endocrine neoplasia type 2 and hereditary medullary thyroid carcinoma, other mutations are classified as variants of uncertain significance (VUS) until the associated clinical phenotype is made clear. Currently, some 46 non-synonymous VUS entries exist in curated archives. In the absence of a gold standard method for predicting phenotype outcomes, this follow up study applies feature selected amino acid physical and chemical properties feeding a Bayes classifier to predict disease association of uncertain gene variants into categories of benign and pathogenic. Algorithm performance and VUS predictions were compared to established phylogenetic based mutation prediction algorithms. Curated outcomes and unpublished RET gene variants with known disease association were used to benchmark predictor performance. Reliable classification of RET uncertain gene variants will augment current clinical information of RET mutations and assist in improving prediction algorithms as knowledge increases.     

Mining the LIPG allelic spectrum reveals the contribution of rare and common regulatory variants to HDL cholesterol

Genome-wide association studies (GWAS) have successfully identified loci associated with quantitative traits, such as blood lipids. Deep resequencing studies are being utilized to catalogue the allelic spectrum at GWAS loci. The goal of these studies is to identify causative variants and missing heritability, including heritability due to low frequency and rare alleles with large phenotypic impact. Whereas rare variant efforts have primarily focused on nonsynonymous coding variants, we hypothesized that noncoding variants in these loci are also functionally important. Using the HDL-C gene LIPG as an example, we explored the effect of regulatory variants identified through resequencing of subjects at HDL-C extremes on gene expression, protein levels, and phenotype. Resequencing a portion of the LIPG promoter and 5' UTR in human subjects with extreme HDL-C, we identified several rare variants in individuals from both extremes. Luciferase reporter assays were used to measure the effect of these rare variants on LIPG expression. Variants conferring opposing effects on gene expression were enriched in opposite extremes of the phenotypic distribution. Minor alleles of a common regulatory haplotype and noncoding GWAS SNPs were associated with reduced plasma levels of the LIPG gene product endothelial lipase (EL), consistent with its role in HDL-C catabolism. Additionally, we found that a common nonfunctional coding variant associated with HDL-C (rs2000813) is in linkage disequilibrium with a 5' UTR variant (rs34474737) that decreases LIPG promoter activity. We attribute the gene regulatory role of rs34474737 to the observed association of the coding variant with plasma EL levels and HDL-C. Taken together, the findings show that both rare and common noncoding regulatory variants are important contributors to the allelic spectrum in complex trait loci.     

Genetic heterogeneity of Borrelia afzelii in the natural focus of the Middle Urals

As shown by sequencing the spacer rrf (5S)--rrl (23S) in 72 isolates of B. afzelii (one of the causative agents of Ixodes tick borne Borrelia infections) and the chromosomal gene coding protein P66 in 22 isolates, that in the natural focus located in the Middle Urals two different genetic subgroups (VS461 and NT28) of this genospecies simultaneously circulate. These subgroups are represented by 5 gene variants (rrf) 5S--(rrl) 23S and 5 allelic variants in gene p66. The latter, similarly to spacer gene variants, are not linked with a definite host and occur in different rrf--rrl variants of the infective agent. At the same time the definite species of vectors and carriers may be the host of several different B. afzelii variants, both in the spacer and in the gene coding protein P66, which maintains the genetic heterogeneity of B. afzelii population in the natural focus.     

Denaturing high-performance liquid chromatography in the detection of ABCA1 gene mutations in familial HDL deficiency

Mutations in the ABCA1 gene are the cause of familial high density lipoprotein deficiency (FHD). Because these mutations are spread over the entire gene, their detection requires the sequencing of all 50 exons. The aim of this study was to validate denaturing high-performance liquid chromatography (DHPLC) in mutation detection as an alternative to systematic sequencing. Exons of the ABCA1 gene were amplified using primers employed for sequencing. Temperatures for DHPLC were deducted from a software and empirically defined for each amplicon. To assess DHPLC reliability, we tested 30 sequence variants found in FHD patients and controls. Combined DHPLC and sequencing was applied to the genotyping of new FHD patients. Most of the amplicons required from two to five temperature conditions to obtain partially denatured DNA over the entire amplicon length. Twenty-nine of the variants found by sequencing were detected by DHPLC (97% sensitivity). The detection of the last variant (in exon 40) required different primers and amplification conditions. DHPLC and sequencing analysis of new FHD patients revealed that all amplicons showing a heteroduplex DHPLC profile contained sequence variants. No variants were detected in amplicons with a homoduplex profile. DHPLC is a sensitive and reliable method for the detection of ABCA1 gene mutations.     

Identification of regulatory variants of carboxylesterase 1 (CES1): A proof-of-concept study for the application of the Allele-Specific Protein Expression (ASPE) assay in identifying cis-acting regulatory genetic polymorphisms

It is challenging to study regulatory genetic variants as gene expression is affected by both genetic polymorphisms and non-genetic regulators. The mRNA allele-specific expression (ASE) assay has been increasingly used for the study of cis-acting regulatory variants because cis-acting variants affect gene expression in an allele-specific manner. However, poor correlations between mRNA and protein expressions were observed for many genes, highlighting the importance of studying gene expression regulation at the protein level. In the present study, we conducted a proof-of-concept study to utilize a recently developed allele-specific protein expression (ASPE) assay to identify the cis-acting regulatory variants of CES1 using a large set of human liver samples. The CES1 gene encodes for carboxylesterase 1 (CES1), the most abundant hepatic hydrolase in humans. Two cis-acting regulatory variants were found to be significantly associated with CES1 ASPE, CES1 protein expression, and its catalytic activity on enalapril hydrolysis in human livers. Compared to conventional gene expression-based approaches, ASPE demonstrated an improved statistical power to detect regulatory variants with small effect sizes since allelic protein expression ratios are less prone to the influence of non-genetic regulators (e.g., diseases and inducers). This study suggests that the ASPE approach is a powerful tool for identifying cis-regulatory variants.     

Introduction to Deep Sequencing and Its Application to Drug Addiction Research with a Focus on Rare Variants

Through linkage analysis, candidate gene approach, and genome-wide association studies (GWAS), many genetic susceptibility factors for substance dependence have been discovered such as the alcohol dehydrogenase gene (ALDH2) for alcohol dependence (AD) and nicotinic acetylcholine receptor (nAChR) subunit variants on chromosomes 8 and 15 for nicotine dependence (ND). However, these confirmed genetic factors contribute only a small portion of the heritability responsible for each addiction. Among many potential factors, rare variants in those identified and unidentified susceptibility genes are supposed to contribute greatly to the missing heritability. Several studies focusing on rare variants have been conducted by taking advantage of next-generation sequencing technologies, which revealed that some rare variants of nAChR subunits are associated with ND in both genetic and functional studies. However, these studies investigated variants for only a small number of genes and need to be expanded to broad regions/genes in a larger population. This review presents an update on recently developed methods for rare-variant identification and association analysis and on studies focused on rare-variant discovery and function related to addictions.     

Identification of novel variants of metadherin in breast cancer

Metadherin (MTDH, also known as AEG-1, and Lyric) has been demonstrated to play a potential role in several significant aspects of tumor progression. It has been reported that overexpression of MTDH is associated with progression of disease and poorer prognosis in breast cancer. However, there are no studies to date assessing variants of the MTDH gene and their potential relationship with breast cancer susceptibility. Thus, we investigated all variants of the MTDH gene and explored the association of the variants with breast cancer development. Our cohort consisted of full-length gene sequencing of 108 breast cancer cases and 100 healthy controls; variants were detected in 11 breast cancer cases and 13 controls. Among the variants detected, 9 novel variants were discovered and 2 were found to be associated with the susceptibility of breast cancer. However, additional studies need to be conducted in larger sample sizes to validate these findings and to further investigate whether these variants are prognostic in breast cancer patients.     

Rare and common regulatory variation in population-scale sequenced human genomes

Population-scale genome sequencing allows the characterization of functional effects of a broad spectrum of genetic variants underlying human phenotypic variation. Here, we investigate the influence of rare and common genetic variants on gene expression patterns, using variants identified from sequencing data from the 1000 genomes project in an African and European population sample and gene expression data from lymphoblastoid cell lines. We detect comparable numbers of expression quantitative trait loci (eQTLs) when compared to genotypes obtained from HapMap 3, but as many as 80% of the top expression quantitative trait variants (eQTVs) discovered from 1000 genomes data are novel. The properties of the newly discovered variants suggest that mapping common causal regulatory variants is challenging even with full resequencing data; however, we observe significant enrichment of regulatory effects in splice-site and nonsense variants. Using RNA sequencing data, we show that 46.2% of nonsynonymous variants are differentially expressed in at least one individual in our sample, creating widespread potential for interactions between functional protein-coding and regulatory variants. We also use allele-specific expression to identify putative rare causal regulatory variants. Furthermore, we demonstrate that outlier expression values can be due to rare variant effects, and we approximate the number of such effects harboured in an individual by effect size. Our results demonstrate that integration of genomic and RNA sequencing analyses allows for the joint assessment of genome sequence and genome function.