Oligonucleotide properties determination and primer designing: a critical examination of predictions

Motivation: Precise prediction of melting temperature (Tm), secondary structures and design of oligonucleotides determine the efficiency and success of experimentation in molecular biology. Availability of a plethora of software and the users unawareness about their limitations compromises the accuracy and reliability of the predictions. Results: Comparative analysis of 56 modules was done for Tm prediction using a large set of oligonucleotide sequences spanning the whole range of GC-content and length. Allawi module of the calculator ‘MELTING’, Nearest Neighbor (NN) of oligo calculator (McLab), NN of Tm Calculation for Oligos (Biomath Calculator, Promega) and HYTHER provided the most precise Tm predictions. A model has also been proposed to calculate the optimum annealing temperature integrating the already reported formulations. Secondary structure predictions of oligonucleotides reveal a large number of structures in contrast to the experimental observations. Of the 11 primer designing tools evaluated, Primer 3 and WebPrimer performed the best for the AT-rich templates, Exon Primer for AT = GC templates, and Primer Design Assistant, Primer3 and Primer Quest for GC-rich templates. This study provides optimal choice for application to the user, increasing the success of a variety of experimentations, especially those that have high-throughput and complex assay designs

Gene prioritization in Type 2 Diabetes using domain interactions and network analysis

<p>Abstract</p> <p>Background</p> <p>Identification of disease genes for Type 2 Diabetes (T2D) by traditional methods has yielded limited success. Based on our previous observation that T2D may result from disturbed protein-protein interactions affected through disrupting modular domain interactions, here we have designed an approach to rank the candidates in the T2D linked genomic regions as plausible disease genes.</p> <p>Results</p> <p>Our approach integrates Weight value (Wv) method followed by prioritization using clustering coefficients derived from domain interaction network. Wv for each candidate is calculated based on the assumption that disease genes might be functionally related, mainly facilitated by interactions among domains of the interacting proteins. The benchmarking using a test dataset comprising of both known T2D genes and non-T2D genes revealed that Wv method had a sensitivity and specificity of 0.74 and 0.96 respectively with 9 fold enrichment. The candidate genes having a Wv > 0.5 were called High Weight Elements (HWEs). Further, we ranked HWEs by using the network property-the clustering coefficient (C_i). Each HWE with a C_i< 0.015 was prioritized as plausible disease candidates (HWEc) as previous studies indicate that disease genes tend to avoid dense clustering (with an average C_iof 0.015). This method further prioritized the identified disease genes with a sensitivity of 0.32 and a specificity of 0.98 and enriched the candidate list by 6.8 fold. Thus, from the dataset of 4052 positional candidates the method ranked 435 to be most likely disease candidates. The gene ontology sharing for the candidates showed higher representation of metabolic and signaling processes. The approach also captured genes with unknown functions which were characterized by network motif analysis.</p> <p>Conclusions</p> <p>Prioritization of positional candidates is essential for cost-effective and an expedited discovery of disease genes. Here, we demonstrate a novel approach for disease candidate prioritization from numerous loci linked to T2D.</p

A CEP215-HSET complex links centrosomes with spindle poles and drives centrosome clustering in cancer.

Numerical centrosome aberrations underlie certain developmental abnormalities and may promote cancer. A cell maintains normal centrosome numbers by coupling centrosome duplication with segregation, which is achieved through sustained association of each centrosome with a mitotic spindle pole. Although the microcephaly- and primordial dwarfism-linked centrosomal protein CEP215 has been implicated in this process, the molecular mechanism responsible remains unclear. Here, using proteomic profiling, we identify the minus end-directed microtubule motor protein HSET as a direct binding partner of CEP215. Targeted deletion of the HSET-binding domain of CEP215 in vertebrate cells causes centrosome detachment and results in HSET depletion at centrosomes, a phenotype also observed in CEP215-deficient patient-derived cells. Moreover, in cancer cells with centrosome amplification, the CEP215-HSET complex promotes the clustering of extra centrosomes into pseudo-bipolar spindles, thereby ensuring viable cell division. Therefore, stabilization of the centrosome-spindle pole interface by the CEP215-HSET complex could promote survival of cancer cells containing supernumerary centrosomes.S.C. is supported by UK Medical Research Council (MC_U105185859). This work was made possible by funding from Cancer Research UK (C14303/A17197). We acknowledge the support of the University of Cambridge and Hutchison Whampoa Ltd.This is the final version of the article. It first appeared from Nature Publishing Group via https://doi.org/10.1038/ncomms1100

Journal Staff

MicroRNAs (miRNAs) play a key role in regulating mRNA expression, and individual miRNAs have been proposed as diagnostic and therapeutic candidates. The identification of such candidates is complicated by the involvement of multiple miRNAs and mRNAs as well as unknown disease topology of the miRNAs. Here, we investigated if disease-associated miRNAs regulate modules of disease-associated mRNAs, if those miRNAs act complementarily or synergistically, and if single or combinations of miRNAs can be targeted to alter module functions. We first analyzed publicly available miRNA and mRNA expression data for five different diseases. Integrated target prediction and network-based analysis showed that the miRNAs regulated modules of disease-relevant genes. Most of the miRNAs acted complementarily to regulate multiple mRNAs. To functionally test these findings, we repeated the analysis using our own miRNA and mRNA expression data from CD4+ T cells from patients with seasonal allergic rhinitis. This is a good model of complex diseases because of its well-defined phenotype and pathogenesis. Combined computational and functional studies confirmed that miRNAs mainly acted complementarily and that a combination of two complementary miRNAs, miR-223 and miR-139-3p, could be targeted to alter disease-relevant module functions, namely, the release of type 2 helper T-cell (Th2) cytokines. Taken together, our findings indicate that miRNAs act complementarily to regulate modules of disease-related mRNAs and can be targeted to alter disease-relevant functions

Molecular determinants underlying functional innovations of TBP and their impact on transcription initiation.

Funder: RCUK | Medical Research Council (MRC); doi: https://doi.org/10.13039/501100000265TATA-box binding protein (TBP) is required for every single transcription event in archaea and eukaryotes. It binds DNA and harbors two repeats with an internal structural symmetry that show sequence asymmetry. At various times in evolution, TBP has acquired multiple interaction partners and different organisms have evolved TBP paralogs with additional protein regions. Together, these observations raise questions of what molecular determinants (i.e. key residues) led to the ability of TBP to acquire new interactions, resulting in an increasingly complex transcriptional system in eukaryotes. We present a comprehensive study of the evolutionary history of TBP and its interaction partners across all domains of life, including viruses. Our analysis reveals the molecular determinants and suggests a unified and multi-stage evolutionary model for the functional innovations of TBP. These findings highlight how concerted chemical changes on a conserved structural scaffold allow for the emergence of complexity in a fundamental biological process

Variable glutamine-rich repeats modulate transcription factor activity

Excessive expansions of glutamine (Q)-rich repeats in various human proteins are known to result in severe neurodegenerative disorders such as Huntington's disease and several ataxias. However, the physiological role of these repeats and the consequences of more moderate repeat variation remain unknown. Here, we demonstrate that Q-rich domains are highly enriched in eukaryotic transcription factors where they act as functional modulators. Incremental changes in the number of repeats in the yeast transcriptional regulator Ssn6 (Cyc8) result in systematic, repeat-length-dependent variation in expression of target genes that result in direct phenotypic changes. The function of Ssn6 increases with its repeat number until a certain threshold where further expansion leads to aggregation. Quantitative proteomic analysis reveals that the Ssn6 repeats affect its solubility and interactions with Tup1 and other regulators. Thus, Q-rich repeats are dynamic functional domains that modulate a regulator's innate function, with the inherent risk of pathogenic repeat expansions

Cis-regulation of microRNA expression by scaffold/matrix-attachment regions

microRNAs (miRNAs) spatio-temporally modulate gene expression; however, very little is known about the regulation of their expression. Here, we hypothesized that the well-known cis-regulatory elements of gene expression, scaffold/matrix-attachment regions (MARs) could modulate miRNA expression. Accordingly, we found MARs to be enriched in the upstream regions of miRNA genes. To determine their role in cell type-specific expression of miRNAs, we examined four individual miRNAs (let-7b, miR-17, miR-93 and miR-221) and the miR-17–92 cluster, known to be overexpressed in neuroblastoma. Our results show that MARs indeed define the cell-specific expression of these miRNAs by tethering the chromatin to nuclear matrix. This is brought about by cell type-specific binding of HMG I/Y protein to MARs that then promotes the local acetylation of histones, serving as boundary elements for gene activation. The binding, chromatin tethering and gene activation by HMG I/Y was not observed in fibroblast control cells but were restricted to neuroblastoma cells. This study implies that the association of MAR binding proteins to MARs could dictate the tissue/context specific regulation of miRNA genes by serving as a boundary element signaling the transcriptional activation

Common Variants in CRP and LEPR Influence High Sensitivity C-Reactive Protein Levels in North Indians

BACKGROUND: High sensitivity C-reactive protein (hsCRP) levels are shown to be influenced by genetic variants in Europeans; however, little is explored in Indian population. METHODS: Herein, we comprehensively evaluated association of all previously reported genetic determinants of hsCRP levels, including 18 cis (proximal to CRP gene) and 73 trans-acting (distal to CRP gene) variants in 4,200 North Indians of Indo-European ethnicity. First, we evaluated association of 91 variants from 12 candidate loci with hsCRP levels in 2,115 North Indians (1,042 non-diabetic subjects and 1,073 patients with type 2 diabetes). Then, cis and trans-acting variants contributing maximally to hsCRP level variation were further replicated in an independent 2,085 North Indians (1,047 patients with type 2 diabetes and 1,038 non-diabetic subjects). RESULTS: We found association of 12 variants from CRP, LEPR, IL1A, IL6, and IL6R with hsCRP levels in non-diabetic subjects. However, only rs3093059-CRP [β = 0.33, P = 9.6×10⁻⁵] and the haplotype harboring rs3093059 risk allele [β = 0.32 µg/mL, P = 1.4×10⁻⁴/P(perm) = 9.0×10⁻⁴] retained significance after correcting for multiple testing. The cis-acting variant rs3093059-CRP had maximum contribution to the variance in hsCRP levels (1.14%). Among, trans-acting variants, rs1892534-LEPR was observed to contribute maximally to hsCRP level variance (0.59%). Associations of rs3093059-CRP and rs1892534-LEPR were confirmed by replication and attained higher significance after meta-analysis [β(meta) = 0.26/0.22; P(meta) = 4.3×10⁻⁷/7.4×10⁻³ and β(meta) = -0.15/-0.12; P(meta) = 2.0×10⁻⁶/1.6×10⁻⁶ for rs3093059 and rs1892534, respectively in non-diabetic subjects and all subjects taken together]. CONCLUSION: In conclusion, we identified rs3093059 in CRP and rs1892534 in LEPR as major cis and trans-acting contributor respectively, to the variance in hsCRP levels in North Indian population

Genome-wide association study for type 2 diabetes in Indians identifies a new susceptibility locus at 2q21.

Indians undergoing socioeconomic and lifestyle transitions will be maximally affected by epidemic of type 2 diabetes (T2D). We conducted a two-stage genome-wide association study of T2D in 12,535 Indians, a less explored but high-risk group. We identified a new type 2 diabetes-associated locus at 2q21, with the lead signal being rs6723108 (odds ratio 1.31; P = 3.32 × 10⁻⁹). Imputation analysis refined the signal to rs998451 (odds ratio 1.56; P = 6.3 × 10⁻¹²) within TMEM163 that encodes a probable vesicular transporter in nerve terminals. TMEM163 variants also showed association with decreased fasting plasma insulin and homeostatic model assessment of insulin resistance, indicating a plausible effect through impaired insulin secretion. The 2q21 region also harbors RAB3GAP1 and ACMSD; those are involved in neurologic disorders. Forty-nine of 56 previously reported signals showed consistency in direction with similar effect sizes in Indians and previous studies, and 25 of them were also associated (P < 0.05). Known loci and the newly identified 2q21 locus altogether explained 7.65% variance in the risk of T2D in Indians. Our study suggests that common susceptibility variants for T2D are largely the same across populations, but also reveals a population-specific locus and provides further insights into genetic architecture and etiology of T2D