Stochastic noise and synchronisation during Dictyostelium aggregation make cAMP oscillations robust

The molecular network, which underlies the oscillations in the concentration of adenosine 3′, 5′-cyclic monophosphate (cAMP) during the aggregation phase of starvation-induced development in Dictyostelium discoideum, achieves remarkable levels of robust performance in the face of environmental variations and cellular heterogeneity. However, the reasons for this robustness remain poorly understood. Tools and concepts from the field of control engineering provide powerful methods for uncovering the mechanisms underlying the robustness of these types of biological systems. Using such methods, two important factors contributing to the robustness of cAMP oscillations in Dictyostelium are revealed. First, stochastic fluctuations in the molecular interactions of the intracellular network, arising from random or directional noise and biological sources, play an important role in preserving stable oscillations in the face of variations in the kinetics of the network. Second, synchronisation of the aggregating cells through the diffusion of extracellular cAMP appears to be a key factor in ensuring robustness to cell-to-cell variations of the oscillatory waves of cAMP observed in Dictyostelium cell cultures. The conclusions have important general implications for the robustness of oscillating biomolecular networks (whether seen at organism, cell, or intracellular levels and including circadian clocks or Ca2+ oscillations, etc.), and suggest that such analysis can be conducted more reliably by using models including stochastic simulations, even in the case where molecular concentrations are very high

Elucidating the mechanisms of cooperative calcium-calmodulin interactions: a structural systems biology approach

BACKGROUND: Calmodulin is an important multifunctional molecule that regulates the activities of a large number of proteins in the cell. Calcium binding induces conformational transitions in calmodulin that make it specifically active to particular target proteins. The precise mechanisms underlying calcium binding to calmodulin are still, however, quite poorly understood. RESULTS: In this study, we adopt a structural systems biology approach and develop a mathematical model to investigate various types of cooperative calcium-calmodulin interactions. We compare the predictions of our analysis with physiological dose-response curves taken from the literature, in order to provide a quantitative comparison of the effects of different mechanisms of cooperativity on calcium-calmodulin interactions. The results of our analysis reduce the gap between current understanding of intracellular calmodulin function at the structural level and physiological calcium-dependent calmodulin target activation experiments. CONCLUSION: Our model predicts that the specificity and selectivity of CaM target regulation is likely to be due to the following factors: variations in the target-specific Ca2+ dissociation and cooperatively effected dissociation constants, and variations in the number of Ca2+ ions required to bind CaM for target activation

Crosstalk between G-protein and Ca2+ pathways switches intracellular cAMP levels

Cyclic adenosine monophosphate and cyclic guanosine monophosphate are universal intracellular messengers whose concentrations are regulated by molecular networks comprised of different isoforms of the synthases adenylate cyclase or guanylate cyclase and the phosphodiesterases which degrade these compounds. In this paper, we employ a systems biology approach to develop mathematical models of these networks that, for the first time, take into account the different biochemical properties of the isoforms involved. To investigate the mechanisms underlying the joint regulation of cAMP and cGMP, we apply our models to analyse the regulation of cilia beat frequency in Paramecium by Ca(2+). Based on our analysis of these models, we propose that the diversity of isoform combinations that occurs in living cells provides an explanation for the huge variety of intracellular processes that are dependent on these networks. The inclusion of both G-protein receptor and Ca(2+)-dependent regulation of AC in our models allows us to propose a new explanation for the switching properties of G-protein subunits involved in nucleotide regulation. Analysis of the models suggests that, depending on whether the G-protein subunit is bound to AC, Ca(2+) can either activate or inhibit AC in a concentration-dependent manner. The resulting analysis provides an explanation for previous experimental results that showed that alterations in Ca(2+) concentrations can either increase or decrease cilia beat frequency over particular Ca(2+) concentration ranges

Computational modelling suggests dynamic interactions between Ca2+, IP3 and G protein-coupled modules are key to robust Dictyostelium aggregation

Under conditions of starvation, Dictyostelium cells begin a programme of development during which they aggregate to form a multicellular structure by chemotaxis, guided by propagating waves of cyclic AMP that are relayed robustly from cell to cell. In this paper, we develop and analyse a new model for the intracellular and extracellular cAMP dependent processes that regulate Dictyostelium migration. The model allows, for the first time, a quantitative analysis of the dynamic interactions between calcium, IP(3) and G protein-dependent modules that are shown to be key to the generation of robust cAMP oscillations in Dictyostelium cells. The model provides a mechanistic explanation for the transient increase in cytosolic free Ca(2+) concentration seen in recent experiments with the application of the calmodulin inhibitor calmidazolium (R24571) to Dictyostelium cells, and also allows elucidation of the effects of varying both the conductivity of stretch-activated channels and the concentration of external phosphodiesterase on the oscillatory regime of an individual cell. A rigorous analysis of the robustness of the new model shows that interactions between the different modules significantly reduce the sensitivity of the resulting cAMP oscillations to variations in the kinetics of different Dictyostelium cells, an essential requirement for the generation of the spatially and temporally synchronised chemoattractant cAMP waves that guide Dictyostelium aggregation

The nature and organization of satellite DNAs in Petunia hybrida, related, and ancestral genomes

IntroductionThe garden petunia, Petunia hybrida (Solanaceae) is a fertile, diploid, annual hybrid species (2n=14) originating from P. axillaris and P. inflata 200 years ago. To understand the recent evolution of the P. hybrida genome, we examined tandemly repeated or satellite sequences using bioinformatic and molecular cytogenetic analysis.MethodsRaw reads from available genomic assemblies and survey sequences of P. axillaris N (PaxiN), P. inflata S6, (PinfS6), P. hybrida (PhybR27) and the here sequenced P. parodii S7 (PparS7) were used for graph and k-mer based cluster analysis of TAREAN and RepeatExplorer. Analysis of repeat specific monomer lengths and sequence heterogeneity of the major tandem repeat families with more than 0.01% genome proportion were complemented by fluorescent in situ hybridization (FISH) using consensus sequences as probes to chromosomes of all four species.ResultsSeven repeat families, PSAT1, PSAT3, PSAT4, PSAT5 PSAT6, PSAT7 and PSAT8, shared high consensus sequence similarity and organisation between the four genomes. Additionally, many degenerate copies were present. FISH in P. hybrida and in the three wild petunias confirmed the bioinformatics data and gave corresponding signals on all or some chromosomes. PSAT1 is located at the ends of all chromosomes except the 45S rDNA bearing short arms of chromosomes II and III, and we classify it as a telomere associated sequence (TAS). It is the most abundant satellite repeat with over 300,000 copies, 0.2% of the genomes. PSAT3 and the variant PSAT7 are located adjacent to the centromere or mid-arm of one to three chromosome pairs. PSAT5 has a strong signal at the end of the short arm of chromosome III in P. axillaris and P.inflata, while in P. hybrida additional interstitial sites were present. PSAT6 is located at the centromeres of chromosomes II and III. PSAT4 and PSAT8 were found with only short arrays.DiscussionThese results demonstrate that (i) repeat families occupy distinct niches within chromosomes, (ii) they differ in the copy number, cluster organization and homogenization events, and that (iii) the recent genome hybridization in breeding P. hybrida preserved the chromosomal position of repeats but affected the copy number of repetitive DNA

Integration of Banana Streak Badnavirus into theMusaGenome: Molecular and Cytogenetic Evidence

AbstractBreeding and tissue culture of certain cultivars of bananas (Musa) have led to high levels of banana streak badnavirus (BSV) infection in progeny from symptomless parents. BSV DNA hybridized to genomic DNA of one such parent, Obino l'Ewai, suggesting integration of viral sequences. Sequencing of clones of Obino l'Ewai genomic DNA revealed an interface between BSV andMusasequences and a complex BSV integrant.In situhybridization revealed two different BSV sequence locations in Obino l'Ewai chromosomes and a complex arrangement of BSV andMusasequences was shown by probing stretched DNA fibers. This is the first report of integrated sequences that possibly lead to a plant pararetrovirus episomal infection by a mechanism differing markedly from animal retroviral systems

Participation of Multifunctional RNA in Replication, Recombination and Regulation of Endogenous Plant Pararetroviruses (EPRVs)

Pararetroviruses, taxon Caulimoviridae, are typical of retroelements with reverse transcriptase and share a common origin with retroviruses and LTR retrotransposons, presumably dating back 1.6 billion years and illustrating the transition from an RNA to a DNA world. After transcription of the viral genome in the host nucleus, viral DNA synthesis occurs in the cytoplasm on the generated terminally redundant RNA including inter- and intra-molecule recombination steps rather than relying on nuclear DNA replication. RNA recombination events between an ancestral genomic retroelement with exogenous RNA viruses were seminal in pararetrovirus evolution resulting in horizontal transmission and episomal replication. Instead of active integration, pararetroviruses use the host DNA repair machinery to prevail in genomes of angiosperms, gymnosperms and ferns. Pararetrovirus integration – leading to Endogenous ParaRetroViruses, EPRVs – by illegitimate recombination can happen if their sequences instead of homologous host genomic sequences on the sister chromatid (during mitosis) or homologous chromosome (during meiosis) are used as template. Multiple layers of RNA interference exist regulating episomal and chromosomal forms of the pararetrovirus. Pararetroviruses have evolved suppressors against this plant defense in the arms race during co-evolution which can result in deregulation of plant genes. Small RNAs serve as signaling molecules for Transcriptional and Post-Transcriptional Gene Silencing (TGS, PTGS) pathways. Different populations of small RNAs comprising 21–24 nt and 18–30 nt in length have been reported for Citrus, Fritillaria, Musa, Petunia, Solanum and Beta. Recombination and RNA interference are driving forces for evolution and regulation of EPRVs

Reticulate evolution in Panicum (Poaceae): the origin of tetraploid broomcorn millet, P. miliaceum.

Panicum miliaceum (broomcorn millet) is a tetraploid cereal, which was among the first domesticated crops, but is now a minor crop despite its high water use efficiency. The ancestors of this species have not been determined; we aimed to identify likely candidates within the genus, where phylogenies are poorly resolved. Nuclear and chloroplast DNA sequences from P. miliaceum and a range of diploid and tetraploid relatives were used to develop phylogenies of the diploid and tetraploid species. Chromosomal in situ hybridization with genomic DNA as a probe was used to characterize the genomes in the tetraploid P. miliaceum and a tetraploid accession of P. repens. In situ hybridization showed that half the chromosomes of P. miliaceum hybridized more strongly with labelled genomic DNA from P. capillare, and half with labelled DNA from P. repens. Genomic DNA probes differentiated two sets of 18 chromosomes in the tetraploid P. repens. Our phylogenetic data support the allotetraploid origin of P. miliaceum, with the maternal ancestor being P. capillare (or a close relative) and the other genome being shared with P. repens. Our P. repens accession was also an allotetraploid with two dissimilar but closely related genomes, the maternal genome being similar to P. sumatrense. Further collection of Panicum species, particularly from the Old World, is required. It is important to identify why the water-efficient P. miliaceum is now of minimal importance in agriculture, and it may be valuable to exploit the diversity in this species and its ancestors

Genome assembly of Musa beccarii shows extensive chromosomal rearrangements and genome expansion during evolution of Musaceae genomes

Background Musa beccarii (Musaceae) is a banana species native to Borneo, sometimes grown as an ornamental plant. The basic chromosome number of Musa species is x = 7, 10, or 11; however, M. beccarii has a basic chromosome number of x = 9 (2n = 2x = 18), which is the same basic chromosome number of species in the sister genera Ensete and Musella. Musa beccarii is in the section Callimusa, which is sister to the section Musa. We generated a high-quality chromosome-scale genome assembly of M. beccarii to better understand the evolution and diversity of genomes within the family Musaceae. Findings The M. beccarii genome was assembled by long-read and Hi-C sequencing, and genes were annotated using both long Iso-seq and short RNA-seq reads. The size of M. beccarii was the largest among all known Musaceae assemblies (∼570 Mbp) due to the expansion of transposable elements and increased 45S ribosomal DNA sites. By synteny analysis, we detected extensive genome-wide chromosome fusions and fissions between M. beccarii and the other Musa and Ensete species, far beyond those expected from differences in chromosome number. Within Musaceae, M. beccarii showed a reduced number of terpenoid synthase genes, which are related to chemical defense, and enrichment in lipid metabolism genes linked to the physical defense of the cell wall. Furthermore, type III polyketide synthase was the most abundant biosynthetic gene cluster (BGC) in M. beccarii. BGCs were not conserved in Musaceae genomes. Conclusions The genome assembly of M. beccarii is the first chromosome-scale genome assembly in the Callimusa section in Musa, which provides an important genetic resource that aids our understanding of the evolution of Musaceae genomes and enhances our knowledge of the pangenome