Computational Modelling of Metastasis Development in Renal Cell Carcinoma

The biology of the metastatic colonization process remains a poorly understood phenomenon. To improve our knowledge of its dynamics, we conducted a modelling study based on multi-modal data from an orthotopic murine experimental system of metastatic renal cell carcinoma. The standard theory of metastatic colonization usually assumes that secondary tumours, once established at a distant site, grow independently from each other and from the primary tumour. Using a mathematical model that translates this assumption into equations, we challenged this theory against our data that included: 1) dynamics of primary tumour cells in the kidney and metastatic cells in the lungs, retrieved by green fluorescent protein tracking, and 2) magnetic resonance images (MRI) informing on the number and size of macroscopic lesions. Critically, when calibrated on the growth of the primary tumour and total metastatic burden, the predicted theoretical size distributions were not in agreement with the MRI observations. Moreover, tumour expansion only based on proliferation was not able to explain the volume increase of the metastatic lesions. These findings strongly suggested rejection of the standard theory, demonstrating that the time development of the size distribution of metastases could not be explained by independent growth of metastatic foci. This led us to investigate the effect of spatial interactions between merging metastatic tumours on the dynamics of the global metastatic burden. We derived a mathematical model of spatial tumour growth, confronted it with experimental data of single metastatic tumour growth, and used it to provide insights on the dynamics of multiple tumours growing in close vicinity. Together, our results have implications for theories of the metastatic process and suggest that global dynamics of metastasis development is dependent on spatial interactions between metastatic lesions.     

Retroviral Oligonucleotide Distributions Correlate with Biased Nucleotide Compositions of Retrovirus Sequences,Suggesting a Duplicative Stepwise Molecular Evolution

A computer-assisted analysis was made of 24 complete nucleotide sequences selected from the vertebrate retroviruses to represent the ten viral groups. The conclusions of this analysis extend and strengthen the previously made hypothesis on the Moloney murine leukemia virus: The evolution of the nucleotide sequence appears to have occurred mainly through at least three overlapping levels of duplication: (1) The distributions of overrepresented (3–6)-mers are consistent with the universal rule of a trend toward TG/CT excess and with the persistence of a certain degree of symmetry between the two strands of DNA. This suggests one or several original tandemly repeated sequences and some inverted duplications. (2) The existence of two general core consensuses at the level of these (3–6)-mers supports the hypothesis of a common evolutionary origin of vertebrate retroviruses. Consensuses more specific to certain sequences are compatible with phylogenetic trees established independently. The consensuses could correspond to intermediary evolutionary stages. (3) Most of the (3–6)-mers with a significantly higher than average frequency appear to be internally repeated (with monomeric or oligomeric internal iterations) and seem to be at least partly the cause of the bias observed by other researchers at the level of retroviral nucleotide composition. They suggest a third evolutionary stage by slippage-like stepwise local duplications. Received: 3 January 1996 / Accepted: 27 March 1996     

Accurate genome relative abundance estimation for closely related species in a metagenomic sample

Background

Metagenomics has a great potential to discover previously unattainable information about microbial communities. An important prerequisite for such discoveries is to accurately estimate the composition of microbial communities. Most of prevalent homology-based approaches utilize solely the results of an alignment tool such as BLAST, limiting their estimation accuracy to high ranks of the taxonomy tree.

Results

We developed a new homology-based approach called Taxonomic Analysis by Elimination and Correction (TAEC), which utilizes the similarity in the genomic sequence in addition to the result of an alignment tool. The proposed method is comprehensively tested on various simulated benchmark datasets of diverse complexity of microbial structure. Compared with other available methods designed for estimating taxonomic composition at a relatively low taxonomic rank, TAEC demonstrates greater accuracy in quantification of genomes in a given microbial sample. We also applied TAEC on two real metagenomic datasets, oral cavity dataset and Crohn’s disease dataset. Our results, while agreeing with previous findings at higher ranks of the taxonomy tree, provide accurate estimation of taxonomic compositions at the species/strain level, narrowing down which species/strains need more attention in the study of oral cavity and the Crohn’s disease.

Conclusions

By taking account of the similarity in the genomic sequence TAEC outperforms other available tools in estimating taxonomic composition at a very low rank, especially when closely related species/strains exist in a metagenomic sample.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2105-15-242) contains supplementary material, which is available to authorized users.     

Mechanisms of peroxisome proliferator-induced DNA hypomethylation in rat liver

Genomic hypomethylation is a consistent finding in both human and animal tumors and mounting experimental evidence suggests a key role for epigenetic events in tumorigenesis. Furthermore, it has been suggested that early changes in DNA methylation and histone modifications may serve as sensitive predictive markers in animal testing for carcinogenic potency of environmental agents. Alterations in metabolism of methyl donors, disturbances in activity and/or expression of DNA methyltransferases, and presence of DNA single-strand breaks could contribute to the loss of cytosine methylation during carcinogenesis; however, the precise mechanisms of genomic hypomethylation induced by chemical carcinogens remain largely unknown. This study examined the mechanism of DNA hypomethylation during hepatocarcinogenesis induced by peroxisome proliferators WY-14,643 (4-chloro-6-(2,3-xylidino)-pyrimidynylthioacetic acid) and DEHP (di-(2-ethylhexyl)phthalate), agents acting through non-genotoxic mode of action. In the liver of male Fisher 344 rats exposed to WY-14,643 (0.1% (w/w), 5 months), the level of genomic hypomethylation increased by approximately 2-fold, as compared to age-matched controls, while in the DEHP group (1.2% (w/w), 5 months) DNA methylation did not change. Global DNA hypomethylation in livers from WY-14,643 group was accompanied by the accumulation of DNA single-strand breaks, increased cell proliferation, and diminished expression of DNA methyltransferase 1, while the metabolism of methyl donors was not affected. In contrast, none of these parameters changed significantly in rats fed DEHP. Since WY-14,643 is much more potent carcinogen than DEHP, we conclude that the extent of loss of DNA methylation may be related to the carcinogenic potential of the chemical agent, and that accumulation of DNA single-strand breaks coupled to the increase in cell proliferation and altered DNA methyltransferase expression may explain genomic hypomethylation during peroxisome proliferator-induced carcinogenesis.     

Two Rubisco activase isoforms may play different roles in photosynthetic heat acclimation in the rice plant

Studies on some plant species have shown that increasing the growth temperature gradually or pretreating with high temperature can lead to obvious photosynthetic acclimation to high temperature. To test whether this acclimation arises from heat adaptation of ribulose 1,5‐bisphosphate carboxylase/oxygenase (Rubisco, EC 4.1.1.39) activation mediated by Rubisco activase (RCA), gene expression of RCA large isoform (RCA_L) and RCA small isoform (RCA_S) in rice was determined using a 4‐day heat stress treatment [40/30°C (day/night)] followed by a 3‐day recovery under control conditions [30/22°C (day/night)]. The heat stress significantly induced the expression of RCA_L as determined by both mRNA and protein levels. Correlative analysis indicated that RCA_S protein content was extremely significantly related to Rubisco initial activity and net photosynthetic rate (Pn) under both heat stress and normal conditions. Immunoblot analysis of the Rubisco–RCA complex revealed that the ratio of RCA_L to Rubisco increased markedly in heat‐acclimated rice leaves. Furthermore, transgenic rice plants expressing enhanced amounts of RCA_L exhibited higher thermotolerance in Pn and Rubisco initial activity and grew better at high temperature than wild‐type (WT) plants and transgenic rice plants expressing enhanced amounts of RCA_S. Under normal conditions, the transgenic rice plants expressing enhanced amounts of RCA_S showed higher Pn and produced more biomass than transgenic rice plants expressing enhanced amounts of RCA_L and wild‐type plants. Together, these suggest that the heat‐induced RCA_L may play an important role in photosynthetic acclimation to moderate heat stress in vivo, while RCA_S plays a major role in maintaining Rubisco initial activity under normal conditions.     

Comparative study of the glycan specificities of cell-bound human tandem-repeat-type galectin-4, -8 and -9

Adhesion/growth-regulatory galectins (gals) exert their functionality by the cis/trans-cross-linking of distinct glycans after initial one-point binding. In order to define the specificity of ensuing association events leading to cross-linking, we recently established a cell-based assay using fluorescent glycoconjugates as flow cytometry probes and tested it on two human gals (gal-1 and -3). Here we present a systematic study of tandem-repeat-type gal-4, -8 and -9 loaded on Raji cells resulting in the following key insights: (i) all three gals bound to oligolactosamines; (ii) binding to ligands with Galβ1-3GlcNAc or Galβ1-3GalNAc as basic motifs was commonly better than that to canonical Galβ1-4GlcNAc; (iii) all three gals bound to 3'-O-sulfated and 3'-sialylated disaccharides mentioned above better than that to parental neutral forms and (iv) histo-blood group ABH antigens were the highest affinity ligands in both the cell and the solid-phase assay. Fine specificity differences were revealed as follows: (i) gal-8 and -9, but not gal-4, bound to disaccharide Galβ1-3GlcNAc; (ii) increase in binding due to negatively charged substituents was marked only in the case of gal-4 and (iii) gal-4 and -8 bound preferably to histo-blood group A glycans, whereas gal-9 targeted B-type glycans. Experiments with single carbohydrate recognition domains (CRDs) of gal-4 showed that the C-CRD preferably bound to ABH glycans, whereas the N-CRD associated with oligolactosamines. In summary, the comparative analysis disclosed the characteristic profiles of glycan reactivity for the accessible CRD of cell-bound gals. These results indicate the distinct sets of functionality for these three members of the same subgroup of human gals.     

From Mesolithic hunters to Iron Age herders: a unique record of woodland use from eastern central Europe (Czech Republic)

Vegetation History and Archaeobotany - In a continuous, perfectly stratified sedimentary sequence which was discovered under a large sandstone overhang in northern Bohemia, Czech Republic, we...     

Structure of a Two-Domain N-Terminal Fragment of Ribosomal Protein L10 from Methanococcus jannaschii Reveals a Specific Piece of the Archaeal Ribosomal Stalk

Ribosomal stalk is involved in the formation of the so-called “GTPase-associated site” and plays a key role in the interaction of ribosome with translation factors and in the control of translation accuracy. The stalk is formed by two or three copies of the L7/L12 dimer bound to the C-terminal tail of protein L10. The N-terminal domain of L10 binds to a segment of domain II of 23S rRNA near the binding site for ribosomal protein L11. The structure of bacterial L10 in complex with three L7/L12 N-terminal dimers has been determined in the isolated state, and the structure of the first third of archaeal L10 bound to domain II of 23S rRNA has been solved within the Haloarcula marismortui 50S ribosomal subunit. A close structural similarity between the RNA-binding domain of archaeal L10 and the RNA-binding domain of bacterial L10 has been demonstrated. In this work, a long RNA-binding N-terminal fragment of L10 from Methanococcus jannaschii has been isolated and crystallized. The crystal structure of this fragment (which encompasses two-thirds of the protein) has been solved at 1.6 Å resolution. The model presented shows the structure of the RNA-binding domain and the structure of the adjacent domain that exist in archaeal L10 and eukaryotic P0 proteins only. Furthermore, our model incorporated into the structure of the H. marismortui 50S ribosomal subunit allows clarification of the structure of the archaeal ribosomal stalk base.     

LncRNA882 regulates leukemia inhibitory factor (LIF) by sponging miR-15b in the endometrial epithelium cells of dairy goat

Despite the fact that long noncoding RNAs (lncRNAs) play roles in almost all biological processes, little is known about their biological function in the endometrium during the formation of endometrial receptivity. In this study, a comprehensive analysis of lncRNAs in goat endometrial tissues on Day 5 (prereceptive endometrium, PE) and Day 15 (receptive endometrium, RE) of pregnancy was performed by using RNA-Seq. As a result, 668 differentially expressed lncRNAs (DELs) were found between the PE and RE. Further study showed that lncRNA882, regulated by estrogen (E2) and progestin (P4), could act as competing endogenous RNAs (ceRNAs) for miR-15b, which inhibited the expression of transforming growth factor-b-activated kinase 1 binding protein 3 (TAB3) and then indirectly regulated the level of leukemia inhibitory factor (LIF). This was helpful for the formation of endometrial receptivity in dairy goats. In conclusion, we elucidated the endometrium lncRNA profiles of PE and RE in dairy goats; lncRNA882 acted as a ceRNA for miR-15b and then indirectly regulated the level of LIF in goat endometrial epithelium cells. Thus, this study helped us to better understand the molecular regulation of endometrial receptivity in dairy goats.