Clustering exact matches of pairwise sequence alignments by weighted linear regression

<p>Abstract</p> <p>Background</p> <p>At intermediate stages of genome assembly projects, when a number of contigs have been generated and their validity needs to be verified, it is desirable to align these contigs to a reference genome when it is available. The interest is not to analyze a detailed alignment between a contig and the reference genome at the base level, but rather to have a rough estimate of where the contig aligns to the reference genome, specifically, by identifying the starting and ending positions of such a region. This information is very useful in ordering the contigs, facilitating post-assembly analysis such as gap closure and resolving repeats. There exist programs, such as BLAST and MUMmer, that can quickly align and identify high similarity segments between two sequences, which, when seen in a dot plot, tend to agglomerate along a diagonal but can also be disrupted by gaps or shifted away from the main diagonal due to mismatches between the contig and the reference. It is a tedious and practically impossible task to visually inspect the dot plot to identify the regions covered by a large number of contigs from sequence assembly projects. A forced global alignment between a contig and the reference is not only time consuming but often meaningless.</p> <p>Results</p> <p>We have developed an algorithm that uses the coordinates of all the exact matches or high similarity local alignments, clusters them with respect to the main diagonal in the dot plot using a weighted linear regression technique, and identifies the starting and ending coordinates of the region of interest.</p> <p>Conclusion</p> <p>This algorithm complements existing pairwise sequence alignment packages by replacing the time-consuming seed extension phase with a weighted linear regression for the alignment seeds. It was experimentally shown that the gain in execution time can be outstanding without compromising the accuracy. This method should be of great utility to sequence assembly and genome comparison projects.</p

ERK1/2 signaling dominates over RhoA signaling in regulating early changes in RNA expression induced by endothelin-1 in neonatal rat cardiomyocytes

Cardiomyocyte hypertrophy is associated with changes in gene expression. Extracellular signal-regulated kinases 1/2 (ERK1/2) and RhoA [activated by hypertrophic agonists (e.g. endothelin-1)] regulate gene expression and are implicated in the response, but their relative significance in regulating the cardiomyocyte transcriptome is unknown. Our aim was to establish the significance of ERK1/2 and/or RhoA in the early cardiomyocyte transcriptomic response to endothelin-1.Cardiomyocytes were exposed to endothelin-1 (1 h) with/without PD184352 (to inhibit ERK1/2) or C3 transferase (C3T, to inhibit RhoA). RNA expression was analyzed using microarrays and qPCR. ERK1/2 signaling positively regulated approximately 65% of the early gene expression response to ET-1 with a small (approximately 2%) negative effect, whereas RhoA signaling positively regulated approximately 10% of the early gene expression response to ET-1 with a greater (approximately 14%) negative contribution. Of RNAs non-responsive to endothelin-1, 66 or 448 were regulated by PD184352 or C3T, respectively, indicating that RhoA had a more significant effect on baseline RNA expression. mRNAs upregulated by endothelin-1 encoded a number of receptor ligands (e.g. Ereg, Areg, Hbegf) and transcription factors (e.g. Abra/Srf) that potentially propagate the response.ERK1/2 dominates over RhoA in the early transcriptomic response to endothelin-1. RhoA plays a major role in maintaining baseline RNA expression but, with upregulation of Abra/Srf by endothelin-1, RhoA may regulate changes in RNA expression over longer times. Our data identify ERK1/2 as a more significant node than RhoA in regulating the early stages of cardiomyocyte hypertrophy

Warped Riemannian metrics for location-scale models

The present paper shows that warped Riemannian metrics, a class of Riemannian metrics which play a prominent role in Riemannian geometry, are also of fundamental importance in information geometry. Precisely, the paper features a new theorem, which states that the Rao-Fisher information metric of any location-scale model, defined on a Riemannian manifold, is a warped Riemannian metric, whenever this model is invariant under the action of some Lie group. This theorem is a valuable tool in finding the expression of the Rao-Fisher information metric of location-scale models defined on high-dimensional Riemannian manifolds. Indeed, a warped Riemannian metric is fully determined by only two functions of a single variable, irrespective of the dimension of the underlying Riemannian manifold. Starting from this theorem, several original contributions are made. The expression of the Rao-Fisher information metric of the Riemannian Gaussian model is provided, for the first time in the literature. A generalised definition of the Mahalanobis distance is introduced, which is applicable to any location-scale model defined on a Riemannian manifold. The solution of the geodesic equation is obtained, for any Rao-Fisher information metric defined in terms of warped Riemannian metrics. Finally, using a mixture of analytical and numerical computations, it is shown that the parameter space of the von Mises-Fisher model of $n$-dimensional directional data, when equipped with its Rao-Fisher information metric, becomes a Hadamard manifold, a simply-connected complete Riemannian manifold of negative sectional curvature, for $n = 2,\ldots,8$. Hopefully, in upcoming work, this will be proved for any value of $n$.Comment: first version, before submissio

Interactions between the Nse3 and Nse4 Components of the SMC5-6 Complex Identify Evolutionarily Conserved Interactions between MAGE and EID Families

The SMC5-6 protein complex is involved in the cellular response to DNA damage. It is composed of 6-8 polypeptides, of which Nse1, Nse3 and Nse4 form a tight sub-complex. MAGEG1, the mammalian ortholog of Nse3, is the founding member of the MAGE (melanoma-associated antigen) protein family and Nse4 is related to the EID (E1A-like inhibitor of differentiation) family of transcriptional repressors.Using site-directed mutagenesis, protein-protein interaction analyses and molecular modelling, we have identified a conserved hydrophobic surface on the C-terminal domain of Nse3 that interacts with Nse4 and identified residues in its N-terminal domain that are essential for interaction with Nse1. We show that these interactions are conserved in the human orthologs. Furthermore, interaction of MAGEG1, the mammalian ortholog of Nse3, with NSE4b, one of the mammalian orthologs of Nse4, results in transcriptional co-activation of the nuclear receptor, steroidogenic factor 1 (SF1). In an examination of the evolutionary conservation of the Nse3-Nse4 interactions, we find that several MAGE proteins can interact with at least one of the NSE4/EID proteins.We have found that, despite the evolutionary diversification of the MAGE family, the characteristic hydrophobic surface shared by all MAGE proteins from yeast to humans mediates its binding to NSE4/EID proteins. Our work provides new insights into the interactions, evolution and functions of the enigmatic MAGE proteins

Diabetes causes marked inhibition of mitochondrial metabolism in pancreatic β-cells

Diabetes is a global health problem caused primarily by the inability of pancreatic β-cells to secrete adequate levels of insulin. The molecular mechanisms underlying the progressive failure of β-cells to respond to glucose in type-2 diabetes remain unresolved. Using a combination of transcriptomics and proteomics, we find significant dysregulation of major metabolic pathways in islets of diabetic βV59M mice, a non-obese, eulipidaemic diabetes model. Multiple genes/proteins involved in glycolysis/gluconeogenesis are upregulated, whereas those involved in oxidative phosphorylation are downregulated. In isolated islets, glucose-induced increases in NADH and ATP are impaired and both oxidative and glycolytic glucose metabolism are reduced. INS-1 β-cells cultured chronically at high glucose show similar changes in protein expression and reduced glucose-stimulated oxygen consumption: targeted metabolomics reveals impaired metabolism. These data indicate hyperglycaemia induces metabolic changes in β-cells that markedly reduce mitochondrial metabolism and ATP synthesis. We propose this underlies the progressive failure of β-cells in diabetes.Peer reviewe

Clinicopathologic characteristics and prognostic factors of ovarian fibrosarcoma: the results of a multi-center retrospective study

<p>Abstract</p> <p>Background</p> <p>Ovarian fibrosarcomas are very rare tumors, and therefore, few case studies have evaluated the prognostic factors of this disease. To our knowledge, this study represents the largest study to evaluate the clinical and pathologic factors associated with ovarian fibrosarcoma patients.</p> <p>Methods</p> <p>Thirty-one cases of ovarian fibrosarcoma were retrospectively reviewed, which included medical records for eight patients, and 23 published case reports from 1995 through 2009. Patient treatment regimens included total hysterectomy with bilateral adnexectomy and an omentectomy (BAO) (n = 9), oophorectomy (OR) (n = 8), chemotherapy (CT) (n = 1), BAO followed by chemotherapy (BAO+CT) (n = 11), BAO followed by radiotherapy (BAO+RT) (n = 1), and oophorectomy followed by radiotherapy (OR + RT) (n = 1).</p> <p>Results</p> <p>The patients of this cohort were staged according to the guidelines of the Federation of Gynecology and Obstetrics (FIGO), with 15, 6, 9, and 1 stage I-IV cases identified, respectively. Mitotic count values were also evaluated from 10 high-power fields (HPFs), and 3 cases had an average mitotic count < 4, 18 cases were between 4 and 10, and 10 cases had an average mitotic count value ≥ 10. The Ki-67 (MIB-1) proliferation index values were grouped according to values that as follows: < 10% (n = 5), between 10% and 50% (n = 9), and ≥ 50% (n = 5). Positive expression of vimentin (100%, 22/22) and negative expression of CD117 (0%, 5/5) were also detected. Moreover, expression of smooth muscle actin (2/18), desmin (1/13), epithelial membrane antigen (0/11), S-100 (1/19), CD99 (0/6), CD34 (1/5), α-inhibin (7/15), estrogen receptor (1/6), and progesterone receptor (1/6) were reported for subsets of the cases examined. After a median follow-up period of 14 months (range, 2-120), the 2-year overall survival rates (OS) and disease-free survival (DFS) rates for all patients were 55.9% and 45.4%, respectively. Cox proportional hazard regression analysis of survival showed that FIGO stage (<it>P </it>= 0.007) and treatment (<it>P </it>= 0.008) were predictive of poor prognosis. Furthermore, patients with stage I tumors that received BAO+CT were associated with a better prognosis.</p> <p>Conclusions</p> <p>Mitotic activity, and cells positive for Ki-67 were identified as important factors in the diagnosis of ovarian fibrosarcoma. Furthermore, FIGO stage and treatment modalities have the potential to be prognostic factors of survival, with BAO followed by adjuvant chemotherapy associated with an improved treatment outcome.</p

The five dimensions of B cell tolerance

B cell tolerance has been generally understood to be an acquired property of the immune system that governs antibody specificity in ways that avoid auto‐toxicity. As useful as this understanding has proved, it fails to fully explain the existence of auto‐reactive specificities in healthy individuals and contribution these may have to health. Mechanisms underlying B cell tolerance are considered to select a clonal repertoire that generates a collection of antibodies that do not bind self, ie tolerance operates more or less in three dimensions that largely spare autologous cells and antigens. Yet, most B lymphocytes in humans and probably in other vertebrates are auto‐reactive and absence of these auto‐reactive B cells is associated with disease. We suggest that auto‐reactivity can be embodied by extending the concept of tolerance by two further dimensions, one of time and circumstance and one that allows healthy cells to actively resist injury. In this novel concept, macromolecular recognition by the B cell receptor leading to deletion, anergy, receptor editing or B cell activation is extended by taking account of the time of development of normal immune responses (4th dimension) and the accommodation (or tolerance) of normal cells to bound antibody, activation of complement, and interaction with inflammatory cells (fifth dimension). We discuss how these dimensions contribute to understanding B cell biology in health or disease.Peer Reviewedhttps://deepblue.lib.umich.edu/bitstream/2027.42/153034/1/imr12813.pdfhttps://deepblue.lib.umich.edu/bitstream/2027.42/153034/2/imr12813_am.pd

A novel class of microRNA-recognition elements that function only within open reading frames.

MicroRNAs (miRNAs) are well known to target 3' untranslated regions (3' UTRs) in mRNAs, thereby silencing gene expression at the post-transcriptional level. Multiple reports have also indicated the ability of miRNAs to target protein-coding sequences (CDS); however, miRNAs have been generally believed to function through similar mechanisms regardless of the locations of their sites of action. Here, we report a class of miRNA-recognition elements (MREs) that function exclusively in CDS regions. Through functional and mechanistic characterization of these 'unusual' MREs, we demonstrate that CDS-targeted miRNAs require extensive base-pairing at the 3' side rather than the 5' seed; cause gene silencing in an Argonaute-dependent but GW182-independent manner; and repress translation by inducing transient ribosome stalling instead of mRNA destabilization. These findings reveal distinct mechanisms and functional consequences of miRNAs that target CDS versus the 3' UTR and suggest that CDS-targeted miRNAs may use a translational quality-control-related mechanism to regulate translation in mammalian cells