Molecular systematics and biogeography of Crawfurdia, Metagentiana and Tripterospermum (Gentianaceae) based on nuclear ribosomal and plastid DNA sequences

Background and Aims The systematic position of the genus Metagentiana and its phylogenetic relationships with Crawfurdia, Gentiana and Tripterospermum have not been explicitly addressed. These four genera belong to one of two subtribes (Gentianinae) of Gentianeae. The aim of this paper is to examine the systematic position of Crawfurdia, Metagentiana and Tripterospermum and to clarify their phylogenetic affinities more clearly using ITS and trnL intron sequences.Methods Nucleotide sequences from the internal transcribed spacers (ITS) of nuclear ribosomal DNA and the plastid DNA trnL (UAA) intron were analysed phylogenetically. Ten of fourteen Metagentiana species were sampled, together with 40 species of other genera in the subtribe Gentianinae.Key Results The data support several previously published conclusions relating to the separation of Metagentiana from Gentiana and its closer relationships to Crawfurdia and Tripterospermum based on studies of gross morphology, floral anatomy, chromosomes, palynology, embryology and previous molecular data. The molecular clock hypothesis for the tested sequences in subtribe Gentianinae was not supported by the data (P < 0.05), so the clock-independent non-parametric rate smoothing method was used to estimate divergence time. This indicates that the separation of Crawfurdia, Metagentiana and Tripterospermum from Gentiana occurred about 11.4-21.4 Mya (million years ago), and the current species of these three genera diverged at times ranging from 0.4 to 6.2 Mya.Conclusions The molecular analyses revealed that Crawfurdia, Metagentiana and Tripterospermum do not merit status as three separate genera, because sampled species of Crawfurdia and Tripterospermum are embedded within Metagentiana. The speciation and rapid radiation of these three genera is likely to have occurred in western China as a result of upthrust of the Himalayas during the late Miocene and the Pleistocene.

In vivo characterization of the circadian clock output gene usp2

Life on earth is subject to the repeated change between day and night periods. All organisms that undergo these alterations have to anticipate consequently the adaptation of their physiology and possess an endogenous periodicity of about 24 hours called circadian rhythm from the Latin circa (about) and diem (day). At the molecular level, virtually all cells of an organism possess a molecular clock which drives rhythmic gene expression and output functions. Besides altered rhythmicity in constant conditions, impaired clock function causes pathophysiological conditions such as diabetes or hypertension. These data unveil a part of the mechanisms underlying the well-described epidemiology of shift work and highlight the function of clock-driven regulatory mechanisms. The post-translational modification of proteins by the ubiquitin polypeptide is a central mechanism to regulate their stability and activity and is capital for clock function. Similarly to the majority of biological processes, it is reversible. Deubiquitylation is carried out by a wide variety of about ninety deubiquitylating enzymes and their function remains poorly understood, especially in vivo. This class of proteolytic enzymes is parted into five families including the Ubiquitin-Specific Proteases (USP), which is the most important with about sixty members. Among them, the Ubiquitin-Specific Protease 2 (Usp2) gene encodes two protein isoforms, USP2-45 and USP2-69. The first is ubiquitously expressed under the control of the circadian clock and displays all features of core clock genes or its closest outputs effectors. Additionally, Usp2-45 was also found to be induced by the mineralocorticoid hormone aldosterone and thought to participate in Na+ reabsorption and blood pressure regulation by Epithelial Na+ Channel ENaC in the kidneys. During my thesis, I aimed to characterize the role of Usp2 in vivo with respect to these two areas, by taking advantage of a total constitutive knockout mouse model. In the first project I aimed to validate the role of USP2-45 in Na+ homeostasis and blood pressure regulation by the kidneys. I found no significant alterations of diurnal Na+ homeostasis and blood pressure in these mice, indicating that Usp2 does not play a substantial role in this process. In urine analyses, we found that our Usp2-KO mice are actually hypercalciuric. In a second project, I aimed to understand the causes of this phenotype. I found that the observed hypercalciuria results essentially from intestinal hyperabsorption. These data reveal a new role for Usp2 as an output effector of the circadian clock in dietary Ca2+ metabolism in the intestine.

A molecular phylogeny of the stingless bee genus Melipona (Hymenoptera: Apidae)

Stingless bees (Meliponini) constitute a diverse group of highly eusocial insects that occur throughout tropical regions around the world. The meliponine genus Melipona is restricted to the New World tropics and has over 50 described species. Melipona, like Apis, possesses the remarkable ability to use representational communication to indicate the location of foraging patches. Although Melipona has been the subject of numerous behavioral, ecological, and genetic studies, the evolutionary history of this genus remains largely unexplored. Here, we implement a multigene phylogenetic approach based on nuclear, mitochondrial, and ribosomal loci, coupled with molecular clock methods, to elucidate the phylogenetic relationships and antiquity of subgenera and species of Melipona. Our phylogenetic analysis resolves the relationship among subgenera and tends to agree with morphology-based classification hypotheses. Our molecular clock analysis indicates that the genus Melipona shared a most recent common ancestor at least similar to 14-17 million years (My) ago. These results provide the groundwork for future comparative analyses aimed at understanding the evolution of complex communication mechanisms in eusocial Apidae. (C) 2010 Elsevier Inc. All rights reserved.

Bayesian coalescent analysis reveals a high rate of molecular evolution in GB virus C

GB virus C/hepatitis G (GBV-C) is an RNA virus of the family Flaviviridae. Despite replicating with an RNA-dependent RNA polymerase, some previous estimates of rates of evolutionary change in GBV-C suggest that it fixes mutations at the anomalously low rate of similar to 100(-7) nucleotide substitution per site, per year. However, these estimates were largely based on the assumption that GBV-C and its close relative GBV-A (New World monkey GB viruses) codiverged with their primate hosts over millions of years. Herein, we estimated the substitution rate of GBV-C using the largest set of dated GBV-C isolates compiled to date and a Bayesian coalescent approach that utilizes the year of sampling and so is independent of the assumption of codivergence. This revealed a rate of evolutionary change approximately four orders of magnitude higher than that estimated previously, in the range of 10(-2) to 10(-3) sub/site/year, and hence in line with those previously determined for RNA viruses in general and the Flaviviridae in particular. In addition, we tested the assumption of host-virus codivergence in GBV-A by performing a reconciliation analysis of host and virus phylogenies. Strikingly, we found no statistical evidence for host-virus codivergence in GBV-A, indicating that substitution rates in the GB viruses should not be estimated from host divergence times.

Molecular phylogeny and zoogeography of the freshwater crayfish genus Cherax Erichson (Decapoda: Parastacidae) in Australia

The evolutionary history and biogeography of freshwater-dependent taxa in Australia is of intrinsic interest given the present-day aridity of this continent. Cherax is the most widespread and one of the most species-rich of Australia's nine freshwater crayfish genera. The phylogenetic relationships amongst 19 of the 23 Australian Cherax were established from mitochondrial DNA sequences representing the 12S rRNA and 16S rRNA gene regions. The relationships among species support an initial east–west separation, followed by a north–south divergence in eastern Australia. Molecular clock estimations suggest that these divergences date back to the Miocene. The phylogenetic relationships support endemic speciation within geographical regions and indicate that long-distance dispersal has not led to recent speciation as previously hypothesized. This new evolutionary scenario is consistent with the climatic history of Australia and the evolutionary history of other similarly distributed freshwater-dependent organisms in Australia.

Caracterização molecular de cepas do Vírus dengue 1 isoladas no Brasil entre os anos de 1994 a 2008

A febre do dengue é uma das mais importantes arboviroses distribuída por todas as áreas tropicais do mundo. O vírus dengue (VDEN) é transmitido principalmente pela picada do mosquito Aedes aegypti infectado. A dispersão do vetor e o aumento do fluxo migratório entre países possibilitaram a ocorrência de grandes epidemias e manifestações clínicas severas, como febre hemorrágica do dengue (FHD) e Síndrome do choque do dengue (SCD). O objetivo deste trabalho foi realizar a caracterização molecular de isolados do VDEN sorotipo 1 (VDEN-1) no Brasil ao nível dos genes estruturais C/prM/M/E de 29 cepas isoladas durante epidemias ocorridas no Brasil no período de 1994 a 2008. A identidade nucleotídica entre as cepas de VDEN-1 do estudo em relação às outras isoladas no Brasil variou de 96,1% a 100%, enquanto o percentual de identidade de aminoácidos foi determinado entre 98,4% a 100%. As diferenças de aminoácidos entre as cepas do estudo, quando comparadas com a cepa FGA/89 (Guiana Francesa), mostraram a presença de importantes substituições não-sinônimas com mudança de caráter bioquímico, tais como os resíduos E297 (Met Tre) e E338 (Ser Leu), sendo necessário estudos para verificar se essas alterações podem ou não estar relacionadas à virulência. A análise filogenética para a proteína E, realizada por meio do método de máxima verossimilhança para as cepas do estudo e outras cepas selecionadas do banco de dados do Genbank, mostraram que as cepas de VDEN-1 isoladas no Brasil desde 1982 pertencem ao genótipo V, corroborando com os resultados publicados anteriormente. O tempo de divergência do VDEN-1, estimado através da hipótese de relógio molecular, mostrou que este vírus teve sua origem a partir de uma linhagem ancestral, há aproximadamente, 113 anos e observou-se ainda que as cepas de VDEN-1 circulantes no Brasil e as provenientes da África possuem um ancestral comum, sendo necessário estudos de filogeografia que mostrem as possíveis rotas de entrada do VDEN-1 no Brasil.

Epidemiologia molecular do Vírus Oropouche (Bunyaviridae: Orthobunyavirus) na Amazônia brasileira

O Virus Oropouche (VORO; Bunyaviridae, Orthobunyavirus) é um dos mais importantes arbovírus que infecta humanos na Amazônia brasileira, e é causador da febre do Oropouche. Entre 1961 e 2009, um grande número de epidemias foi registrado em diferentes centros urbanos dos Estados Brasileiros do Acre, Amapá, Amazonas, Maranhão, Pará, Rondônia e Tocantins, e também no Panamá, Peru e Trinidad & Tobago. Este trabalho teve por objetivo desenvolver um estudo retrospectivo dos aspectos epidemiológicos e moleculares do VORO enfatizando sua distribuição, a dinâmica das epidemias ocorridas no período, bem como a dispersão de diferentes genótipos na América Latina e no Brasil como contribuição à epidemiologia molecular do VORO. Para tanto 66 isolamentos do VORO pertencentes ao acervo do Instituto Evandro Chagas foram propagados em camundongos e em cultura de células VERO, seguida da extração do RNA viral e obtenção do cDNA por RTPCR; os amplicons foram purificados e submetidos ao sequenciamento nucleotídico para análises moleculares e evolução, incluindo o rearranjo genético, estudo de relógio molecular e análise de dispersão viral. Foi demonstrada a presença de quatro linhagens distintas do VORO na Amazônia brasileira (genótipos I, II, III e IV), sendo os genótipos I e II, respectivamente os mais frequentemente encontrados em áreas da Amazônia ocidental e oriental. Esses e o genótipo III estão constantemente evoluindo, mediante o mecanismo “boom and boost” que resulta na emergência seguida de substituição das sublinhagens (subgenótipos) circulantes por outras mais recentes. O genótipo III do VORO, previamente encontrado somente no Panamá, foi descrito na Amazônia e Sudeste do Brasil. Os dados obtidos pela análise filogenética comparativa das topologias para os segmentos PRNA e MRNA sugerem que o VORO utiliza o rearranjo genético como mecanismo de geração de biodiversidade viral, sendo o genótipo I o mais estável e o II o mais instável e, portanto, mais sensível às pressões evolutivas; foi reconhecido um novo genótipo do VORO neste estudo em amostras isoladas em Manaus no ano de 1980, que foi denominado de genótipo IV. O estudo do relógio molecular mostrou que a emergência do VORO se deu no Estado do Pará provavelmente há 223 anos e daí ao longo dos anos se dispersou pela PanAmazônia bem como para o Caribe, sendo que o genótipo I foi o que originou os demais genótipos do VORO.

Molecular and phylogenetic analyses of human Parvovirus B19 isolated from Brazilian patients with sickle cell disease and beta-thalassemia major and healthy blood donors

Human Parvovirus B19 (B19V) is a recognized cause of life-threatening conditions among patients with hemoglobinopathies. This study investigates B19V infection in patients with sickle cell disease and beta-thalassemia using different experimental approaches. A total of 183 individuals (144 with sickle cell disease and 39 with beta-thalassemia major) and 100 healthy blood donors were examined for B19V using anti-B19V IgG enzyme immunoassay, quantitative PCR, DNA sequencing, and phylogenetic analysis. Viremia was documented in 18.6% of patients and 1% of donors, and was generally characterized by low viral load (VL); however, acute infections were also observed. Anti-B19V IgG was detected in 65.9% of patients with sickle cell disease and in 60% of donors, whereas the patients with thalassemia exhibited relatively low seroreactivity. The seroprevalence varied among the different age groups. In patients, it progressively increased with age, whereas in donors it reached a plateau. Based on partial NS1 fragments, all isolates detected were classified as subgenotype 1A with a tendency to elicit genetically complex infections. Interestingly, quasispecies occurred in the plasma of not only patients but also donors with even higher heterogeneity. The partial NS1 sequence examined did not exhibit positive selection. Quantitation of B19V with a conservative probe is a technically and practically useful approach. The extensive spread of B19V subgenotype 1A in patients and donors and its recent introduction into the countryside of the Sao Paulo State, Brazil were demonstrated; however, it is difficult to establish a relationship between viral sequences and the clinical outcomes of the infection. J. Med. Virol. 84:16521665, 2012. (c) 2012 Wiley Periodicals, Inc.

Molecular phylogeny and historical zoogeography of four genera of Eastern Atlantic skates (Rajiformes, Rajidae)

The aim of this study was to reconstruct a solid phylogeny of four genera of the Rajidae family (Chondrichthyans: Batoidea) using a concatenated alignment of mtDNA genes. Then use the resultant tree to estimate divergence time between taxa based on molecular clock and fossil calibration and conduct biogeographic analysis. The intent was to prove that the actual distribution of species of Eastern Atlantic and Mediterranean skates is due to a series of vicariant events. The species considered belongs to two different tribe: Rajini (Raja and Dipturus) and Amblyrajini (Leucoraja and Rajella). The choice of this genera is due to their high presence in the area of interest and to the richness of endemic species. The results show that despite the ancient origin of Rajidae (97 MYA), the Eastern Atlantic and Mediterranean faunas originated more recently, during Middle Miocene-Late Pliocene, after the closure of connection between these areas and the Indo-Pacific ocean (15 MYA). The endemic species of the Mediterranean (Raja asterias, R. radula, R. polystigma and Leucoraja melitensis) originated after the Messinian salinity crisis (7-5 MYA), when the recolonization of the basin occurred, and are still maintained in allopatric distribution by the presence of biogeographic barriers. Moreover from 4 to 2.6 MYA we can observe the formation of sister species for Raja, Leucoraja and Rajella, one of which has a Northern distribution, and the other has a Southern distribution (R. clavata vs R. straeleni, L. wallacei vs L. naevus, R. fyllae vs R. caudaspinosa and R. kukujevi vs R. leopardus + R. barnardi). The Quaternary and present oceanographic discontinuities that occur along the western African continental shelf (e.g., Cape Blanc and the Angola–Benguela Front) might contribute to the maintenance of low or null levels of gene flow between these closely related siblings species. Also sympatric speciation must be invoked to explain the evolution of skates, for example for the division between R. leopardus and R. barnardi. The speciation processes followed a south-to-north pathways for Dipturus and a north-to-south pathways for Raja, Leucoraja and Rajella underling that the evolution of the genera occurred independently. In the end, it is conceivable that the evolutionary pathways of the tribes followed the costal line during the gondwana fragmentation. The results demonstrate that the evolution of this family is characterized by a series of parallel and independent speciation events, strictly correlated to the tectonic movement of continental masses and paleogeographic and paleoclimatic events and so can be explained by a panbiogeographical (vicariance) model.

Bacillus anthracis: Molecular taxonomy, population genetics, phylogeny and patho-evolution

Bacillus anthracis, the etiological agent of anthrax, manifests a particular bimodal lifestyle. This bacterial species alternates between short replication phases of 20-40 generations that strictly require infection of the host, normally causing death, interrupted by relatively long, mostly dormant phases as spores in the environment. Hence, the B. anthracis genome is highly homogeneous. This feature and the fact that strains from nearly all parts of the world have been analysed for canonical single nucleotide polymorphisms (canSNPs) and variable number tandem repeats (VNTRs) has allowed the development of molecular epidemiological and molecular clock models to estimate the age of major diversifications in the evolution of B. anthracis and to trace the global spread of this pathogen, which was mostly promoted by movement of domestic cattle with settlers and by international trade of contaminated animal products. From a taxonomic and phylogenetic point of view, B. anthracis is a member of the Bacillus cereus group. The differentiation of B. anthracis from B. cereus sensu strict, solely based on chromosomal markers, is difficult. However, differences in pathogenicity clearly differentiate B. anthracis from B. cereus and are marked by the strict presence of virulence genes located on the two virulence plasmids pXO1 and pXO2, which both are required by the bacterium to cause anthrax. Conversely, anthrax-like symptoms can also be caused by organisms with chromosomal features that are more closely related to B. cereus, but which carry these virulence genes on two plasmids that largely resemble the B. anthracis virulence plasmids. (C) 2011 Elsevier B.V. All rights reserved.

MOLECULAR AND GENOMIC BASED INSIGHTS INTO THE EVOLUTION OF ENTEROCOCCUS FAECIUM FROM COMMENSAL TO HOSPITAL-ADAPTED PATHOGEN

The basis for the recent transition of Enterococcus faecium from a primarily commensal organism to one of the leading causes of hospital-acquired infections in the United States is not yet understood. To address this, the first part of my project assessed isolates from early outbreaks in the USA and South America using sequence analysis, colony hybridizations, and minimal inhibitory concentrations (MICs) which showed clinical isolates possess virulence and antibiotic resistance determinants that are less abundant or lacking in community isolates. I also revealed that the level of ampicillin resistance increased over time in clinical strains. By sequencing the pbp5 gene, I demonstrated an ~5% difference in the pbp5 gene between strains with MICs <4ug/ml and those with MICs >4µg/ml, but no specific sequence changes correlated with increases in MICs within the latter group. A 3-10% nucleotide difference was also seen in three other genes analyzed, which suggested the existence of two distinct subpopulations of E. faecium. This led to the second part of my project analyzing concatenated core gene sequences, SNPs, the 16S rRNA, and phylogenetics of 21 E. faecium genomes confirming two distinct clades; a community-associated (CA) clade and hospital-associated (HA) clade. Molecular clock calculations indicate that these two clades likely diverged ~ 300,000 to > 1 million years ago, long before the modern antibiotic era. Genomic analysis also showed that, in addition to core genomic differences, HA E. faecium harbor specific accessory genetic elements that may confer selection advantages over CA E. faecium. The third part of my project discovered 6 E. faecium genes with the newly identified “WxL” domain. My analyses, using RT-PCR, western blots, patient sera, whole-cell ELISA, and immunogold electron microscopy, indicated that E. faecium WxL genes exist in operons, encode bacterial cell surface localized proteins, that WxL proteins are antigenic in humans, and are more exposed on the surface of clinical isolates versus community isolates (even though they are ubiquitous in both clades). ELISAs and BIAcore analyses also showed that proteins encoded by these operons bind several different host extracellular matrix proteins, as well as to each other, suggesting a novel cell-surface complex. In summary, my studies provide new insights into the evolution of E. faecium by showing that there are two distantly related clades; one being more successful in the hospital setting. My studies also identified operons encoding WxL proteins whose characteristics could also contribute to colonization and virulence within this species.

Sociality and the rate of molecular evolution

The molecular clock does not tick at a uniform rate in all taxa but maybe influenced by species characteristics. Eusocial species (those with reproductive division of labor) have been predicted to have faster rates of molecular evolution than their nonsocial relatives because of greatly reduced effective population size; if most individuals in a population are nonreproductive and only one or few queens produce all the offspring, then eusocial animals could have much lower effective population sizes than their solitary relatives, which should increase the rate of substitution of nearly neutral mutations. An earlier study reported faster rates in eusocial honeybees and vespid wasps but failed to correct for phylogenetic nonindependence or to distinguish between potential causes of rate variation. Because sociality has evolved independently in many different lineages, it is possible to conduct a more wide-ranging study to test the generality of the relationship. We have conducted a comparative analysis of 25 phylogenetically independent pairs of social lineages and their nonsocial relatives, including bees, wasps, ants, termites, shrimps, and mole rats, using a range of available DNA sequences (mitochondrial and nuclear DNA coding for proteins and RNAs, and nontranslated sequences). By including a wide range of social taxa, we were able to test whether there is a general influence of sociality on rates of molecular evolution and to test specific predictions of the hypothesis: (1) that social species have faster rates because they have reduced effective population sizes; (2) that mitochondrial genes would show a greater effect of sociality than nuclear genes; and (3) that rates of molecular evolution should be correlated with the degree of sociality. We find no consistent pattern in rates of molecular evolution between social and nonsocial lineages and no evidence that mitochondrial genes show faster rates in social taxa. However, we show that the most highly eusocial Hymenoptera do have faster rates than their nonsocial relatives. We also find that social parasites (that utilize the workers from related species to produce their own offspring) have faster rates than their social relatives, which is consistent with an effect of lower effective population size on rate of molecular evolution. Our results illustrate the importance of allowing for phylogenetic nonindependence when conducting investigations of determinants of variation in rate of molecular evolution.

A timescale for diatom evolution based on four molecular markers: reassessment of ghost lineages and major steps defining diatom evolution.

ABSTRACT. – Phylogenies and molecular clocks of the diatoms have largely been inferred from SSU rDNA sequences. A new phylogeny of diatoms was estimated using four gene markers SSU and LSU rDNA rbcL and psbA (total 4352 bp) with 42 diatom species. The four gene trees analysed with a maximum likelihood (ML) and Baysian (BI) analysis recovered a monophyletic origin of the new diatom classes with high bootstrap support, which has been controversial with single gene markers using single outgroups and alignments that do not take secondary structure of the SSU gene into account. The divergence time of the classes were calculated from a ML tree in the MultliDiv Time program using a Bayesian estimation allowing for simultaneous constraints from the fossil record and varying rates of molecular evolution of different branches in the phylogenetic tree. These divergence times are generally in agreement with those proposed by other clocks using single genes with the exception that the pennates appear much earlier and suggest a longer Cretaceous fossil record that has yet to be sampled. Ghost lineages (i.e. the discrepancy between first appearance (FA) and molecular clock age of origin from an extant taxon) were revealed in the pennate lineage, whereas those ghost lineages in the centric lineages previously reported by others are reviewed and referred to earlier literature.

A timescale for diatom evolution based on four molecular markers: reassessment of ghost lineages and major steps defining diatom evolution.

ABSTRACT. – Phylogenies and molecular clocks of the diatoms have largely been inferred from SSU rDNA sequences. A new phylogeny of diatoms was estimated using four gene markers SSU and LSU rDNA rbcL and psbA (total 4352 bp) with 42 diatom species. The four gene trees analysed with a maximum likelihood (ML) and Baysian (BI) analysis recovered a monophyletic origin of the new diatom classes with high bootstrap support, which has been controversial with single gene markers using single outgroups and alignments that do not take secondary structure of the SSU gene into account. The divergence time of the classes were calculated from a ML tree in the MultliDiv Time program using a Bayesian estimation allowing for simultaneous constraints from the fossil record and varying rates of molecular evolution of different branches in the phylogenetic tree. These divergence times are generally in agreement with those proposed by other clocks using single genes with the exception that the pennates appear much earlier and suggest a longer Cretaceous fossil record that has yet to be sampled. Ghost lineages (i.e. the discrepancy between first appearance (FA) and molecular clock age of origin from an extant taxon) were revealed in the pennate lineage, whereas those ghost lineages in the centric lineages previously reported by others are reviewed and referred to earlier literature.

Generalized binomial Tau-leap method for biochemical kinetics incorporating both delay and intrinsic noise

The delay stochastic simulation algorithm (DSSA) by Barrio et al. [Plos Comput. Biol.2, 117–E (2006)] was developed to simulate delayed processes in cell biology in the presence of intrinsic noise, that is, when there are small-to-moderate numbers of certain key molecules present in a chemical reaction system. These delayed processes can faithfully represent complex interactions and mechanisms that imply a number of spatiotemporal processes often not explicitly modeled such as transcription and translation, basic in the modeling of cell signaling pathways. However, for systems with widely varying reaction rate constants or large numbers of molecules, the simulation time steps of both the stochastic simulation algorithm (SSA) and the DSSA can become very small causing considerable computational overheads. In order to overcome the limit of small step sizes, various τ-leap strategies have been suggested for improving computational performance of the SSA. In this paper, we present a binomial τ- DSSA method that extends the τ-leap idea to the delay setting and avoids drawing insufficient numbers of reactions, a common shortcoming of existing binomial τ-leap methods that becomes evident when dealing with complex chemical interactions. The resulting inaccuracies are most evident in the delayed case, even when considering reaction products as potential reactants within the same time step in which they are produced. Moreover, we extend the framework to account for multicellular systems with different degrees of intercellular communication. We apply these ideas to two important genetic regulatory models, namely, the hes1 gene, implicated as a molecular clock, and a Her1/Her 7 model for coupled oscillating cells.