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.     

The mutational dynamics of the SARS-CoV-2 virus in serial passages in vitro

Since its outbreak in 2019, Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) keeps surprising the medical community by evolving diverse immune escape mutations in a rapid and effective manner. To gain deeper insight into mutation frequency and dynamics, we isolated ten ancestral strains of SARS-CoV-2 and performed consecutive serial incubation in ten replications in a suitable and common cell line and subsequently analysed them using RT-qPCR and whole genome sequencing. Along those lines we hoped to gain fundamental insights into the evolutionary capacity of SARS-CoV-2 in vitro. Our results identified a series of adaptive genetic changes, ranging from unique convergent substitutional mutations and hitherto undescribed insertions. The region coding for spike proved to be a mutational hotspot, evolving a number of mutational changes including the already known substitutions at positions S:484 and S:501. We discussed the evolution of all specific adaptations as well as possible reasons for the seemingly inhomogeneous potential of SARS-CoV-2 in the adaptation to cell culture. The combination of serial passage in vitro with whole genome sequencing uncovers the immense mutational potential of some SARS-CoV-2 strains. The observed genetic changes of SARS-CoV-2 in vitro could not be explained solely by selectively neutral mutations but possibly resulted from the action of directional selection accumulating favourable genetic changes in the evolving variants, along the path of increasing potency of the strain. Competition among a high number of quasi-species in the SARS-CoV-2 in vitro population gene pool may reinforce directional selection and boost the speed of evolutionary change.     

Revisiting the phylogeography and demography of European badgers (Meles meles) based on broad sampling,multiple markers and simulations

Although the phylogeography of European mammals has been extensively investigated since the 1990s, many studies were limited in terms of sampling distribution, the number of molecular markers used and the analytical techniques employed, frequently leading to incomplete postglacial recolonisation scenarios. The broad-scale genetic structure of the European badger (Meles meles) is of interest as it may result from historic restriction to glacial refugia and/or recent anthropogenic impact. However, previous studies were based mostly on samples from western Europe, making it difficult to draw robust conclusions about the location of refugia, patterns of postglacial expansion and recent demography. In the present study, continent-wide sampling and analyses with multiple markers provided evidence for two glacial refugia (Iberia and southeast Europe) that contributed to the genetic variation observed in badgers in Europe today. Approximate Bayesian computation provided support for a colonisation of Scandinavia from both Iberian and southeastern refugia. In the whole of Europe, we observed a decline in genetic diversity with increasing latitude, suggesting that the reduced diversity in the peripheral populations resulted from a postglacial expansion processes. Although MSVAR v.1.3 also provided evidence for recent genetic bottlenecks in some of these peripheral populations, the simulations performed to estimate the method''s power to correctly infer the past demography of our empirical populations suggested that the timing and severity of bottlenecks could not be established with certainty. We urge caution against trying to relate demographic declines inferred using MSVAR with particular historic or climatological events.     

Cyclooxygenase-2 contributes to elevated renin in the early postnatal period in rats

We asked whether cyclooxygenase (COX) activity controls the renin-angiotensin system in the postnatal period. During kidney development, renin peaked at postnatal days 0-1 at the mRNA, tissue protein [renal renin concentration (RRC)], and plasma renin concentration (PRC) levels and was widely expressed along preglomerular vessels. PRC and renin mRNA expression was elevated until weaning in the 4th postnatal week compared with adult rats. Renocortical COX-2 was restricted to Tamm-Horsfall protein-positive cells in the thick ascending limb of Henle's loop, and cortical COX-2 mRNA and protein expression were elevated along with PRC in the 2nd and 3rd postnatal weeks. In contrast, cortical COX-1 expression was constant, but medullary COX-1 expression increased eightfold from the 1st to 4th postnatal week. A COX-2-selective blocker, parecoxib, and a nonselective blocker, indomethacin, given in a period with COX-2 induction from postnatal day 6 to day 12, markedly decreased PRC, but not renin mRNA or RRC. Inhibition of angiotensin AT(1) receptors by candesartan from postnatal day 1 to day 5 increased COX-2 mRNA (2.5-fold), protein, and distribution, renin mRNA (7-fold) and PRC (20- to 70-fold), but had no influence on COX-1 mRNA. Thus, due to very low levels of expression, COX-2 is unlikely to be responsible for the birth peak of renin, but COX-2 activity supports renin secretion later in the suckling period. ANG II negatively feeds back on renocortical COX-2 expression in the 1st postnatal days with high activity of the renin system. We suggest that suckling in the rat is correlated to an enhanced, COX-2-mediated, secretory activity of renin-producing juxtaglomerular cells.     

A record of the Plain Tiger,Danaus chrysippus,in Turkey

Summary

In this publication, a subdivision of the genus Annenica O. Boettger, 1877 sensu Auct. is proposed. One result of our study is the recognition of the genus Sprattia O. Boettger, 1883 for most of the South Anatolian species of the former genus Armenica. The genus Astrogena Szekeres, 1970, which was synonymised by H. NORDSIECK (1975: 98) with Armenica is redefined and used as a subgenus of Armenica for some species. Extensive recent collections contain new species and subspecies of this genus, which are described here.     

DNA oligonucleotides with A, T, G or C opposite an abasic site: structure and dynamics

Abasic sites are common DNA lesions resulting from spontaneous depurination and excision of damaged nucleobases by DNA repair enzymes. However, the influence of the local sequence context on the structure of the abasic site and ultimately, its recognition and repair, remains elusive. In the present study, duplex DNAs with three different bases (G, C or T) opposite an abasic site have been synthesized in the same sequence context (5′-CCA AAG₆ XA₈C CGG G-3′, where X denotes the abasic site) and characterized by 2D NMR spectroscopy. Studies on a duplex DNA with an A opposite the abasic site in the same sequence has recently been reported [Chen,J., Dupradeau,F.-Y., Case,D.A., Turner,C.J. and Stubbe,J. (2007) Nuclear magnetic resonance structural studies and molecular modeling of duplex DNA containing normal and 4′-oxidized abasic sites. Biochemistry, 46, 3096–3107]. Molecular modeling based on NMR-derived distance and dihedral angle restraints and molecular dynamics calculations have been applied to determine structural models and conformational flexibility of each duplex. The results indicate that all four duplexes adopt an overall B-form conformation with each unpaired base stacked between adjacent bases intrahelically. The conformation around the abasic site is more perturbed when the base opposite to the lesion is a pyrimidine (C or T) than a purine (G or A). In both the former cases, the neighboring base pairs (G6-C21 and A8-T19) are closer to each other than those in B-form DNA. Molecular dynamics simulations reveal that transient H-bond interactions between the unpaired pyrimidine (C20 or T20) and the base 3′ to the abasic site play an important role in perturbing the local conformation. These results provide structural insight into the dynamics of abasic sites that are intrinsically modulated by the bases opposite the abasic site.     

The FRP1 F-box gene has different functions in sexuality, pathogenicity and metabolism in three fungal pathogens

Plant‐pathogenic fungi employ a variety of infection strategies; as a result, fungi probably rely on different sets of proteins for successful infection. The F‐box protein Frp1, only present in filamentous fungi belonging to the Sordariomycetes, Leotiomycetes and Dothideomycetes, is required for nonsugar carbon catabolism and pathogenicity in the root‐infecting fungus Fusarium oxysporum. To assess the role of Frp1 in other plant‐pathogenic fungi, FRP1 deletion mutants were generated in Fusarium graminearum and Botrytis cinerea, and their phenotypes were analysed. Deletion of FgFRP1 in F. graminearum led to impaired infection of barley roots, but not of aerial plant parts. Deletion of BcFRP1 in B. cinerea did not show any effect on pathogenicity. Sexual reproduction, however, was impaired in both F. graminearum and B. cinerea FRP1 deletion mutants. The mutants of all three fungi displayed different phenotypes when grown on an array of carbon sources. The F. oxysporum and B. cinerea deletion mutants showed opposite growth phenotypes on sugar and nonsugar carbon sources. Replacement of FoFRP1 in F. oxysporum with the B. cinerea BcFRP1 resulted in the restoration of pathogenicity, but also in a switch from impaired growth on nonsugar carbon sources to impaired growth on sugar carbon sources. This effect could be ascribed in part to the B. cinerea BcFRP1 promoter sequence. In conclusion, the function of the F‐box protein Frp1, despite its high sequence conservation, is not conserved between different fungi, leading to differential requirements for pathogenicity and carbon source utilization.