Elucidating the cancer-specific genetic alteration spectrum of glioblastoma derived cell lines from whole exome and RNA sequencing

Cell lines derived from tumor tissues have been used as a valuable system to study gene regulation and cancer development. Comprehensive characterization of the genetic background of cell lines could provide clues on novel genes responsible for carcinogenesis and help in choosing cell lines for particular studies. Here, we have carried out whole exome and RNA sequencing of commonly used glioblastoma (GBM) cell lines (U87, T98G, LN229, U343, U373 and LN18) to unearth single nucleotide variations (SNVs), indels, differential gene expression, gene fusions and RNA editing events. We obtained an average of 41,071 SNVs out of which 1,594 (3.88%) were potentially cancer-specific. The cell lines showed frequent SNVs and indels in some of the genes that are known to be altered in GBM-EGFR, TP53, PTEN, SPTA1 and NF1. Chromatin modifying genes-ATRX, MLL3, MLL4, SETD2 and SRCAP also showed alterations. While no cell line carried IDH1 mutations, five cell lines showed hTERT promoter activating mutations with a concomitant increase in hTERT transcript levels. Five significant gene fusions were found of which NUP93-CYB5B was validated. An average of 18,949 RNA editing events was also obtained. Thus we have generated a comprehensive catalogue of genetic alterations for six GBM cell lines.

Elucidating the cancer-specific genetic alteration spectrum of glioblastoma derived cell lines from whole exome and RNA sequencing

Cell lines derived from tumor tissues have been used as a valuable system to study gene regulation and cancer development. Comprehensive characterization of the genetic background of cell lines could provide clues on novel genes responsible for carcinogenesis and help in choosing cell lines for particular studies. Here, we have carried out whole exome and RNA sequencing of commonly used glioblastoma (GBM) cell lines (U87, T98G, LN229, U343, U373 and LN18) to unearth single nucleotide variations (SNVs), indels, differential gene expression, gene fusions and RNA editing events. We obtained an average of 41,071 SNVs out of which 1,594 (3.88%) were potentially cancer-specific. The cell lines showed frequent SNVs and indels in some of the genes that are known to be altered in GBM-EGFR, TP53, PTEN, SPTA1 and NF1. Chromatin modifying genes-ATRX, MLL3, MLL4, SETD2 and SRCAP also showed alterations. While no cell line carried IDH1 mutations, five cell lines showed hTERT promoter activating mutations with a concomitant increase in hTERT transcript levels. Five significant gene fusions were found of which NUP93-CYB5B was validated. An average of 18,949 RNA editing events was also obtained. Thus we have generated a comprehensive catalogue of genetic alterations for six GBM cell lines.

Formation of linear plasmid multimers promoted by the phage lambda Red-system in lon mutants of Escherichia coli

We report here the formation of plasmid linear multimers promoted by the Red-system of phage lambda using a multicopy plasmid comprised of lambda red alpha and red beta genes, under the control of the lambda cI857 repressor. Our observations have revealed that the multimerization of plasmid DNA is dependent on the red beta and recA genes, suggesting a concerted role for these functions in the formation of plasmid multimers. The formation of multimers occurred in a recBCD+ sbcB+ xthA+ lon genetic background at a higher frequency than in the isogenic lon+ host cells. The multimers comprised tandem repeats of monomer plasmid DNA. Treatment of purified plasmid DNA with exonuclease III revealed the presence of free double-chain ends in the molecules. Determination of the size of multimeric DNA, by pulse field gel electrophoresis, revealed that the bulk of the DNA was in the range 50-240 kb, representing approximately 5-24 unit lengths of monomeric plasmid DNA. We provide a conceptual framework for Red-system-promoted formation and enhanced accumulation of plasmid linear multimers in lon mutants of E. coli.

Involvement of the Global Regulator H-NS in the Survival of Escherichia coli in Stationary Phase

Long-term batch cultures of Escherichia coli grown in nutrient-rich medium accumulate mutations that provide a growth advantage in the stationary phase (GASP). We have examined the survivors of prolonged stationary phase to identify loci involved in conferring a growth advantage and show that a mutation in the hns gene causing reduced activity of the global regulator H-NS confers a GASP phenotype under specific conditions. The hns-66 allele bears a point mutation within the termination codon of the H-NS open reading frame, resulting in a longer protein that is partially functional. Although isolated from a long-term stationary-phase culture of the parent carrying the rpoS819 allele that results in reduced RpoS activity, the hns-66 survivor showed a growth disadvantage in the early stationary phase (24 to 48 h) when competed against the parent. The hns-66 mutant is also unstable and reverts at a high frequency in the early stationary phase by accumulating second-site suppressor mutations within the ssrA gene involved in targeting aberrant proteins for proteolysis. The mutant was more stable and showed a moderate growth advantage in combination with the rpoS819 allele when competed against a 21-day-old parent. These studies show that H-NS is a target for mutations conferring fitness gain that depends on the genetic background as well as on the stage of the stationary phase.

SUB1 Plays a Negative Role during Starvation Induced Sporulation Program in Saccharomyces cerevisiae

Saccharomyces cerevisiae Sub1 is involved in several cellular processes such as, transcription initiation, elongation, mRNA processing and DNA repair. It has also been reported to provide cellular resistance during conditions of oxidative DNA damage and osmotic stress. Here, we report a novel role of SUB1 during starvation stress-induced sporulation, which leads to meiosis and spore formation in diploid yeast cells. Deletion of SUB1 gene significantly increased sporulation efficiency as compared to the wild-type cells in S288c genetic background. Whereas, the sporulation functions of the sub1(Y66A) missense mutant were similar to Sub1. SUB1 transcript and protein levels are downregulated during sporulation, in highly synchronized and sporulation proficient wild-type SK1 cells. The changes in Sub1 levels during sporulation cascade correlate with the induction of middle sporulation gene expression. Deletion of SUB1 increased middle sporulation gene transcript levels with no effect on their induction kinetics. In wild-type cells, Sub1 associates with chromatin at these loci in a temporal pattern that correlates with their enhanced gene expression seen in sub1. cells. We show that SUB1 genetically interacts with HOS2, which led us to speculate that Sub1 might function with Set3 repressor complex during sporulation. Positive Cofactor 4, human homolog of Sub1, complemented the sub1. sporulation phenotype, suggesting conservation of function. Taken together, our results suggest that SUB1 acts as a negative regulator of sporulation.

A linear programming approach for estimating the structure of a sparse linear genetic network from transcript profiling data

Background: A genetic network can be represented as a directed graph in which a node corresponds to a gene and a directed edge specifies the direction of influence of one gene on another. The reconstruction of such networks from transcript profiling data remains an important yet challenging endeavor. A transcript profile specifies the abundances of many genes in a biological sample of interest. Prevailing strategies for learning the structure of a genetic network from high-dimensional transcript profiling data assume sparsity and linearity. Many methods consider relatively small directed graphs, inferring graphs with up to a few hundred nodes. This work examines large undirected graphs representations of genetic networks, graphs with many thousands of nodes where an undirected edge between two nodes does not indicate the direction of influence, and the problem of estimating the structure of such a sparse linear genetic network (SLGN) from transcript profiling data. Results: The structure learning task is cast as a sparse linear regression problem which is then posed as a LASSO (l1-constrained fitting) problem and solved finally by formulating a Linear Program (LP). A bound on the Generalization Error of this approach is given in terms of the Leave-One-Out Error. The accuracy and utility of LP-SLGNs is assessed quantitatively and qualitatively using simulated and real data. The Dialogue for Reverse Engineering Assessments and Methods (DREAM) initiative provides gold standard data sets and evaluation metrics that enable and facilitate the comparison of algorithms for deducing the structure of networks. The structures of LP-SLGNs estimated from the INSILICO1, INSILICO2 and INSILICO3 simulated DREAM2 data sets are comparable to those proposed by the first and/or second ranked teams in the DREAM2 competition. The structures of LP-SLGNs estimated from two published Saccharomyces cerevisae cell cycle transcript profiling data sets capture known regulatory associations. In each S. cerevisiae LP-SLGN, the number of nodes with a particular degree follows an approximate power law suggesting that its degree distributions is similar to that observed in real-world networks. Inspection of these LP-SLGNs suggests biological hypotheses amenable to experimental verification. Conclusion: A statistically robust and computationally efficient LP-based method for estimating the topology of a large sparse undirected graph from high-dimensional data yields representations of genetic networks that are biologically plausible and useful abstractions of the structures of real genetic networks. Analysis of the statistical and topological properties of learned LP-SLGNs may have practical value; for example, genes with high random walk betweenness, a measure of the centrality of a node in a graph, are good candidates for intervention studies and hence integrated computational – experimental investigations designed to infer more realistic and sophisticated probabilistic directed graphical model representations of genetic networks. The LP-based solutions of the sparse linear regression problem described here may provide a method for learning the structure of transcription factor networks from transcript profiling and transcription factor binding motif data.

A Haplotype at the UCP1 Gene Locus Contributes to Genetic Risk for Type 2 Diabetes in Asian Indians (CURES-72)

Background: The gene encoding for uncoupling protein-1 (UCP1) is considered to be a candidate gene for type 2 diabetes because of its role in thermogenesis and energy expenditure. The objective of the study was to examine whether genetic variations in the UCP1 gene are associated with type 2 diabetes and its related traits in Asian Indians. Methods: The study subjects, 810 type 2 diabetic subjects and 990 normal glucose tolerant (NGT) subjects, were chosen from the Chennai Urban Rural Epidemiological Study (CURES), an ongoing population-based study in southern India. The polymorphisms were genotyped using the polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) method. Linkage disequilibrium (LD) was estimated from the estimates of haplotypic frequencies. Results: The three polymorphisms, namely -3826A -> G, an A -> C transition in the 5'-untranslated region (UTR) and Met229Leu, were not associated with type 2 diabetes. However, the frequency of the A-C-Met (-3826A -> G-5'UTR A -> C-Met229Leu) haplotype was significantly higher among the type 2 diabetic subjects (2.67%) compared with the NGT subjects (1.45%, P < 0.01). The odds ratio for type 2 diabetes for the individuals carrying the haplotype A-C-Met was 1.82 (95% confidence interval, 1.29-2.78, P = 0.009). Conclusions: The haplotype, A-C-Met, in the UCP1 gene is significantly associated with the increased genetic risk for developing type 2 diabetes in Asian Indians.

No-reference image quality assessment using modified extreme learning machine classifier

In this paper, we present a machine learning approach to measure the visual quality of JPEG-coded images. The features for predicting the perceived image quality are extracted by considering key human visual sensitivity (HVS) factors such as edge amplitude, edge length, background activity and background luminance. Image quality assessment involves estimating the functional relationship between HVS features and subjective test scores. The quality of the compressed images are obtained without referring to their original images ('No Reference' metric). Here, the problem of quality estimation is transformed to a classification problem and solved using extreme learning machine (ELM) algorithm. In ELM, the input weights and the bias values are randomly chosen and the output weights are analytically calculated. The generalization performance of the ELM algorithm for classification problems with imbalance in the number of samples per quality class depends critically on the input weights and the bias values. Hence, we propose two schemes, namely the k-fold selection scheme (KS-ELM) and the real-coded genetic algorithm (RCGA-ELM) to select the input weights and the bias values such that the generalization performance of the classifier is a maximum. Results indicate that the proposed schemes significantly improve the performance of ELM classifier under imbalance condition for image quality assessment. The experimental results prove that the estimated visual quality of the proposed RCGA-ELM emulates the mean opinion score very well. The experimental results are compared with the existing JPEG no-reference image quality metric and full-reference structural similarity image quality metric.

Timing the Emergence of Resistance to Anti-HIV Drugs with Large Genetic Barriers

New antiretroviral drugs that offer large genetic barriers to resistance, such as the recently approved inhibitors of HIV-1 protease, tipranavir and darunavir, present promising weapons to avert the failure of current therapies for HIV infection. Optimal treatment strategies with the new drugs, however, are yet to be established. A key limitation is the poor understanding of the process by which HIV surmounts large genetic barriers to resistance. Extant models of HIV dynamics are predicated on the predominance of deterministic forces underlying the emergence of resistant genomes. In contrast, stochastic forces may dominate, especially when the genetic barrier is large, and delay the emergence of resistant genomes. We develop a mathematical model of HIV dynamics under the influence of an antiretroviral drug to predict the waiting time for the emergence of genomes that carry the requisite mutations to overcome the genetic barrier of the drug. We apply our model to describe the development of resistance to tipranavir in in vitro serial passage experiments. Model predictions of the times of emergence of different mutant genomes with increasing resistance to tipranavir are in quantitative agreement with experiments, indicating that our model captures the dynamics of the development of resistance to antiretroviral drugs accurately. Further, model predictions provide insights into the influence of underlying evolutionary processes such as recombination on the development of resistance, and suggest guidelines for drug design: drugs that offer large genetic barriers to resistance with resistance sites tightly localized on the viral genome and exhibiting positive epistatic interactions maximally inhibit the emergence of resistant genomes.

Structural Basis of Drug Resistance by Genetic Variants of HIV Type 1 Clade C Protease from India

Using computer modeling of three-dimensional structures and structural information available on the crystal structures of HIV-1 protease, we investigated the structural effects of mutations, in treatment-naive and treatment-exposed individuals from India and postulated mechanisms of resistance in clade C variants. A large number of models (14) have been generated by computational mutation of the available crystal structures of drug bound proteases. Localized energy minimization was carried out in and around the sites of mutation in order to optimize the geometry of interactions present. Most of the mutations result in structural differences at the flap that favors the semiopen state of the enzyme. Some of the mutations were also found to confer resistance by affecting the geometry of the active site. The E35D mutation affects the flap structure in clade B strains and E35N and E35K mutation, seen in our modeled strains, have a more profound effect. Common polymorphisms at positions 36 and 63 in clade C also affected flap structure. Apart from a few other residues Gln-58, Asn-83, Asn-88, and Gln-92 and their interactions are important for the transition from the closed to the open state. Development of protease inhibitors by structure-based design requires investigation of mechanisms operative for clade C to improve the efficacy of therapy.

Contribution from normal and starburst galaxies to the extragalactic gamma-ray background (EGRB)

The extragalactic diffuse emission at gamma-ray energies has interesting cosmological implications since these photons suffer little or no attenuation during their propagation from the site of origin. The emission could originate from either truly diffuse processes or from unresolved point sources such as AGNs, normal galaxies and starburst galaxies. Here, we examine the unresolved point source origin of the extragalactic gamma-ray background emission from normal galaxies and starburst galaxies. gamma-ray emission from normal galaxies is primarily coming from cosmic-ray interactions with interstellar matter and radiation (similar to 90%) along with a small contribution from discrete point sources (similar to 10%). Starburst galaxies are expected to have enhanced supernovae activity which leads to higher cosmic-ray densities, making starburst galaxies sufficiently luminous at gamma-ray energies to be detected by the current gamma-ray mission(Fermi Gamma-ray Space Telescope).

Genetic analysis of foot-and-mouth disease virus serotype A of Indian origin and detection of positive selection and recombination in leader protease- and capsid-coding regions

The leader protease (L-pro) and capsid-coding sequences (P1) constitute approximately 3 kb of the foot-and-mouth disease virus (FMDV). We studied the phylogenetic relationship of 46 FMDV serotype A isolates of Indian origin collected during the period 1968-2005 and also eight vaccine strains using the neighbour-joining tree and Bayesian tree methods. The viruses were categorized under three major groups - Asian, Euro-South American and European. The Indian isolates formed a distinct genetic group among the Asian isolates. The Indian isolates were further classified into different genetic subgroups (<5% divergence). Post-1995 isolates were divided into two subgroups while a few isolates which originated in the year 2005 from Andhra Pradesh formed a separate group. These isolates were closely related to the isolates of the 1970s. The FMDV isolates seem to undergo reverse mutation or onvergent evolution wherein sequences identical to the ancestors are present in the isolates in circulation. The eight vaccine strains included in the study were not related to each other and belonged to different genetic groups. Recombination was detected in the L-pro region in one isolate (A IND 20/82) and in the VP1 coding 1D region in another isolate (A RAJ 21/96). Positive selection was identified at aa positions 23 in the L-pro (P<0.05; 0.046*) and at aa 171 in the capsid protein VP1 (P<0.01; 0.003**).

Inbreeding and the incidence of childhood genetic disorders in Karnataka, South India

Consanguineous marriages are strongly favoured among the populations of South India. In a study conducted on 407 infants and children, a total of 35 genetic diseases was diagnosed in 63 persons: 44 with single gene defects, 12 with polygenic disorders, and seven with Down's syndrome. The coefficient of inbreeding of the total study group, F = 0.0414, was significantly higher than that previously calculated for the general population, F = 0.0271, and autosomal recessive disorders formed the largest single disease category diagnosed. The results suggest that long term inbreeding may not have resulted in appreciable elimination of recessive lethals and sub-lethals from the gene pool.