Teosinte Pollen Drive guides maize diversification and domestication by RNAi.

Selfish genetic elements contribute to hybrid incompatibility and bias or 'drive' their own transmission1,2. Chromosomal drive typically functions in asymmetric female meiosis, whereas gene drive is normally post-meiotic and typically found in males. Here, using single-molecule and single-pollen genome sequencing, we describe Teosinte Pollen Drive, an instance of gene drive in hybrids between maize (Zea mays ssp. mays) and teosinte mexicana (Z. mays ssp. mexicana) that depends on RNA interference (RNAi). 22-nucleotide small RNAs from a non-coding RNA hairpin in mexicana depend on Dicer-like 2 (Dcl2) and target Teosinte Drive Responder 1 (Tdr1), which encodes a lipase required for pollen viability. Dcl2, Tdr1 and the hairpin are in tight pseudolinkage on chromosome 5, but only when transmitted through the male. Introgression of mexicana into early cultivated maize is thought to have been critical to its geographical dispersal throughout the Americas3, and a tightly linked inversion in mexicana spans a major domestication sweep in modern maize4. A survey of maize traditional varieties and sympatric populations of teosinte mexicana reveals correlated patterns of admixture among unlinked genes required for RNAi on at least four chromosomes that are also subject to gene drive in pollen from synthetic hybrids. Teosinte Pollen Drive probably had a major role in maize domestication and diversification, and offers an explanation for the widespread abundance of 'self' small RNAs in the germ lines of plants and animals.

Dengue virus infection modifies mosquito blood-feeding behavior to increase transmission to the host.

Mosquito blood-feeding behavior is a key determinant of the epidemiology of dengue viruses (DENV), the most-prevalent mosquito-borne viruses. However, despite its importance, how DENV infection influences mosquito blood-feeding and, consequently, transmission remains unclear. Here, we developed a high-resolution, video-based assay to observe the blood-feeding behavior of Aedes aegypti mosquitoes on mice. We then applied multivariate analysis on the high-throughput, unbiased data generated from the assay to ordinate behavioral parameters into complex behaviors. We showed that DENV infection increases mosquito attraction to the host and hinders its biting efficiency, the latter resulting in the infected mosquitoes biting more to reach similar blood repletion as uninfected mosquitoes. To examine how increased biting influences DENV transmission to the host, we established an in vivo transmission model with immuno-competent mice and demonstrated that successive short probes result in multiple transmissions. Finally, to determine how DENV-induced alterations of host-seeking and biting behaviors influence dengue epidemiology, we integrated the behavioral data within a mathematical model. We calculated that the number of infected hosts per infected mosquito, as determined by the reproduction rate, tripled when mosquito behavior was influenced by DENV infection. Taken together, this multidisciplinary study details how DENV infection modulates mosquito blood-feeding behavior to increase vector capacity, proportionally aggravating DENV epidemiology. By elucidating the contribution of mosquito behavioral alterations on DENV transmission to the host, these results will inform epidemiological modeling to tailor improved interventions against dengue.

De novo evolution of transmissible tumours in hydra.

While most cancers are not transmissible, there are rare cases where cancer cells can spread between individuals and even across species, leading to epidemics. Despite their significance, the origins of such cancers remain elusive due to late detection in host populations. Using Hydra oligactis, which exhibits spontaneous tumour development that in some strains became vertically transmitted, this study presents the first experimental observation of the evolution of a transmissible tumour. Specifically, we assessed the initial vertical transmission rate of spontaneous tumours and explored the potential for optimizing this rate through artificial selection. One of the hydra strains, which evolved transmissible tumours over five generations, was characterized by analysis of cell type and bacteriome, and assessment of life-history traits. Our findings indicate that tumour transmission can be immediate for some strains and can be enhanced by selection. The resulting tumours are characterized by overproliferation of large interstitial stem cells and are not associated with a specific bacteriome. Furthermore, despite only five generations of transmission, these tumours induced notable alterations in host life-history traits, hinting at a compensatory response. This work, therefore, makes the first contribution to understanding the conditions of transmissible cancer emergence and their short-term consequences for the host.

Approaches and methods to study wildlife cancer.

The last few years have seen a surge of interest from field ecologists and evolutionary biologists to study neoplasia and cancer in wildlife. This contributes to the One Health Approach, which investigates health issues at the intersection of people, wild and domestic animals, together with their changing environments. Nonetheless, the emerging field of wildlife cancer is currently constrained by methodological limitations in detecting cancer using non-invasive sampling. In addition, the suspected differential susceptibility and resistance of species to cancer often make the choice of a unique model species difficult for field biologists. Here, we provide an overview of the importance of pursuing the study of cancer in non-model organisms and we review the currently available methods to detect, measure and quantify cancer in the wild, as well as the methodological limitations to be overcome to develop novel approaches inspired by diagnostic techniques used in human medicine. The methodology we propose here will help understand and hopefully fight this major disease by generating general knowledge about cancer, variation in its rates, tumour-suppressor mechanisms across species as well as its link to life history and physiological characters. Moreover, this is expected to provide key information about cancer in wildlife, which is a top priority due to the accelerated anthropogenic change in the past decades that might favour cancer progression in wild populations.

Sociodemographic Variables Can Guide Prioritized Testing Strategies for Epidemic Control in Resource-Limited Contexts.

BACKGROUND: Targeted surveillance allows public health authorities to implement testing and isolation strategies when diagnostic resources are limited, and can be implemented via the consideration of social network topologies. However, it remains unclear how to implement such surveillance and control when network data are unavailable. METHODS: We evaluated the ability of sociodemographic proxies of degree centrality to guide prioritized testing of infected individuals compared to known degree centrality. Proxies were estimated via readily available sociodemographic variables (age, gender, marital status, educational attainment, household size). We simulated severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) epidemics via a susceptible-exposed-infected-recovered individual-based model on 2 contact networks from rural Madagascar to test applicability of these findings to low-resource contexts. RESULTS: Targeted testing using sociodemographic proxies performed similarly to targeted testing using known degree centralities. At low testing capacity, using proxies reduced infection burden by 22%-33% while using 20% fewer tests, compared to random testing. By comparison, using known degree centrality reduced the infection burden by 31%-44% while using 26%-29% fewer tests. CONCLUSIONS: We demonstrate that incorporating social network information into epidemic control strategies is an effective countermeasure to low testing capacity and can be implemented via sociodemographic proxies when social network data are unavailable.

SARS-CoV-2 BW lineage, a fast-growing Omicron variant from southeast Mexico bearing relevant escape mutations.

PURPOSE: The swift expansion of the BW.1 SARS-CoV-2 variant coincided with a rapid increase of COVID-19 cases occurring in Southeast Mexico in October, 2022, which marked the start of Mexico's sixth epidemiological wave. In Yucatan, up to 92% (58 of 73) of weekly sequenced genomes between epidemiological week 42 and 47 were identified as either BW.1 or its descendant, BW.1.1 in the region, during the last trimester of 2022. In the current study, a comprehensive genomic comparison was carried out to characterize the evolutionary history of the BW lineage, identifying its origins and its most important mutations. METHODS: An alignment of all the genomes of the BW lineage and its parental BA.5.6.2 variant was carried out to identify their mutations. A phylogenetic and ancestral sequence reconstruction analysis with geographical inference, as well as a longitudinal analysis of point mutations, were performed to trace back their origin and contrast them with key RBD mutations in variant BQ.1, one of the fastest-growing lineages to date. RESULTS: Our ancestral reconstruction analysis portrayed Mexico as the most probable origin of the BW.1 and BW.1.1 variants. Two synonymous substitutions, T7666C and C14599T, support their Mexican origin, whereas other two mutations are specific to BW.1: S:N460K and ORF1a:V627I. Two additional substitutions and a deletion are found in its descending subvariant, BW.1.1. Mutations found in the receptor binding domain, S:K444T, S:L452R, S:N460K, and S:F486V in BW.1 have been reported to be relevant for immune escape and are also key mutations in the BQ.1 lineage. CONCLUSIONS: BW.1 appears to have arisen in the Yucatan Peninsula in Southeast Mexico sometime around July 2022 during the fifth COVID-19 wave. Its rapid growth may be in part explained by the relevant escape mutations also found in BQ.1.

Transmissible cancers in mammals and bivalves: How many examples are there?: Predictions indicate widespread occurrence.

Transmissible cancers are elusive and understudied parasitic life forms caused by malignant clonal cells (nine lineages are known so far). They emerge by completing sequential steps that include breaking cell cooperation, evade anti-cancer defences and shedding cells to infect new hosts. Transmissible cancers impair host fitness, and their importance as selective force is likely largely underestimated. It is, therefore, crucial to determine how common they might be in the wild. Here, we draw a parallel between the steps required for a transmissible cancer to emerge and the steps required for an intelligent civilisation to emerge in the Milky Way using a modified Drake equation. Using numerical analyses, we estimate the potential number of extant marine and bivalve species in which transmissible cancers might exist. Our results suggest that transmissible cancers are more common than expected, and that new lineages can be found by screening a large number of species.

Transmissible cancer and the evolution of sex.

The origin and subsequent maintenance of sex and recombination are among the most elusive and controversial problems in evolutionary biology. Here, we propose a novel hypothesis, suggesting that sexual reproduction not only evolved to reduce the negative effects of the accumulation of deleterious mutations and processes associated with pathogen and/or parasite resistance but also to prevent invasion by transmissible selfish neoplastic cheater cells, henceforth referred to as transmissible cancer cells. Sexual reproduction permits systematic change of the multicellular organism's genotype and hence an enhanced detection of transmissible cancer cells by immune system. Given the omnipresence of oncogenic processes in multicellular organisms, together with the fact that transmissible cancer cells can have dramatic effects on their host fitness, our scenario suggests that the benefits of sex and concomitant recombination will be large and permanent, explaining why sexual reproduction is, despite its costs, the dominant mode of reproduction among eukaryotes.

Can postfertile life stages evolve as an anticancer mechanism?

Why a postfertile stage has evolved in females of some species has puzzled evolutionary biologists for over 50 years. We propose that existing adaptive explanations have underestimated in their formulation an important parameter operating both at the specific and the individual levels: the balance between cancer risks and cancer defenses. During their life, most multicellular organisms naturally accumulate oncogenic processes in their body. In parallel, reproduction, notably the pregnancy process in mammals, exacerbates the progression of existing tumors in females. When, for various ecological or evolutionary reasons, anticancer defenses are too weak, given cancer risk, older females could not pursue their reproduction without triggering fatal metastatic cancers, nor even maintain a normal reproductive physiology if the latter also promotes the growth of existing oncogenic processes, e.g., hormone-dependent malignancies. At least until stronger anticancer defenses are selected for in these species, females could achieve higher inclusive fitness by ceasing their reproduction and/or going through menopause (assuming that these traits are easier to select than anticancer defenses), thereby limiting the risk of premature death due to metastatic cancers. Because relatively few species experience such an evolutionary mismatch between anticancer defenses and cancer risks, the evolution of prolonged life after reproduction could also be a rare, potentially transient, anticancer adaptation in the animal kingdom.

A mechanistic and data-driven reconstruction of the time-varying reproduction number: Application to the COVID-19 epidemic.

The effective reproduction number Reff is a critical epidemiological parameter that characterizes the transmissibility of a pathogen. However, this parameter is difficult to estimate in the presence of silent transmission and/or significant temporal variation in case reporting. This variation can occur due to the lack of timely or appropriate testing, public health interventions and/or changes in human behavior during an epidemic. This is exactly the situation we are confronted with during this COVID-19 pandemic. In this work, we propose to estimate Reff for the SARS-CoV-2 (the etiological agent of the COVID-19), based on a model of its propagation considering a time-varying transmission rate. This rate is modeled by a Brownian diffusion process embedded in a stochastic model. The model is then fitted by Bayesian inference (particle Markov Chain Monte Carlo method) using multiple well-documented hospital datasets from several regions in France and in Ireland. This mechanistic modeling framework enables us to reconstruct the temporal evolution of the transmission rate of the COVID-19 based only on the available data. Except for the specific model structure, it is non-specifically assumed that the transmission rate follows a basic stochastic process constrained by the observations. This approach allows us to follow both the course of the COVID-19 epidemic and the temporal evolution of its Reff(t). Besides, it allows to assess and to interpret the evolution of transmission with respect to the mitigation strategies implemented to control the epidemic waves in France and in Ireland. We can thus estimate a reduction of more than 80% for the first wave in all the studied regions but a smaller reduction for the second wave when the epidemic was less active, around 45% in France but just 20% in Ireland. For the third wave in Ireland the reduction was again significant (>70%).

Identification of factors associated with residual malaria transmission using school-based serological surveys in settings pursuing elimination.

BACKGROUND: Targeted research on residual malaria transmission is important to improve strategies in settings pursuing elimination, where transmission reductions prove challenging. This study aimed to detect and characterize spatial heterogeneity and factors associated with Plasmodium falciparum infections and exposure, P. falciparum apical membrane antigen 1 (PfAMA1) antibody (Ab) response, in the Central Highlands of Madagascar (CHL). METHODS: From May to July 2014, a cross-sectional school-based survey was carried out in 182 fokontany (villages) within 7 health districts of the CHL. Rapid diagnostic tests (RDTs) and a bead-based immunoassay including PfAMA1 antigen biomarker were used to estimate malaria prevalence and seroprevalence, respectively. Local Moran's I index was used to detect spatial "hotspots". Remotely sensed environmental data-temperature, vegetation indices, land covers, and elevation-were used in multivariable mixed-effects logistic regression models to characterize factors associated with malaria infection and cumulative exposure. RESULTS: Among 6,293 school-children ages 2-14 years surveyed, RDT prevalence was low at 0.8% (95% CI 0.6-1.1%), while PfAMA1 Ab seroprevalence was 7.0% (95% CI 6.4-7.7%). Hotspots of PfAMA1 Ab seroprevalence were observed in two districts (Ankazobe and Mandoto). Seroprevalence increased for children living > 5 km from a health centre (adjusted odds ratio (OR) = 1.6, 95% CI 1.2-2.2), and for those experiencing a fever episode in the previous 2 weeks (OR 1.7, 95% CI 1.2-2.4), but decreased at higher elevation (for each 100-m increase, OR = 0.7, 95% CI 0.6-0.8). A clear age pattern was observed whereby children 9-10 years old had an OR of 1.8 (95% CI 1.2-2.4), children 11-12 years an OR of 3.7 (95% CI 2.8-5.0), and children 13-14 years an OR of 5.7 (95% CI 4.0-8.0) for seropositivity, compared with younger children (2-8 years). CONCLUSION: The use of serology in this study provided a better understanding of malaria hotspots and associated factors, revealing a pattern of higher transmission linked to geographical barriers in health care access. The integration of antibody-assays into existing surveillance activities could improve exposure assessment, and may help to monitor the effectiveness of malaria control efforts and adapt elimination interventions.

Retrospective analysis of SARS-CoV-2 omicron invasion over delta in French regions in 2021-22: a status-based multi-variant model.

BACKGROUND: SARS-CoV-2 is a rapidly spreading disease affecting human life and the economy on a global scale. The disease has caused so far more then 5.5 million deaths. The omicron outbreak that emerged in Botswana in the south of Africa spread around the globe at further increased rates, and caused unprecedented SARS-CoV-2 infection incidences in several countries. At the start of December 2021 the first omicron cases were reported in France. METHODS: In this paper we investigate the spreading potential of this novel variant relatively to the delta variant that was also in circulation in France at that time. Using a dynamic multi-variant model accounting for cross-immunity through a status-based approach, we analyze screening data reported by Santé Publique France over 13 metropolitan French regions between 1st of December 2021 and the 30th of January 2022. During the investigated period, the delta variant was replaced by omicron in all metropolitan regions in approximately three weeks. The analysis conducted retrospectively allows us to consider the whole replacement time window and compare regions with different times of omicron introduction and baseline levels of variants' transmission potential. As large uncertainties regarding cross-immunity among variants persist, uncertainty analyses were carried out to assess its impact on our estimations. RESULTS: Assuming that 80% of the population was immunized against delta, a cross delta/omicron cross-immunity of 25% and an omicron generation time of 3.5 days, the relative strength of omicron to delta, expressed as the ratio of their respective reproduction rates, [Formula: see text], was found to range between 1.51 and 1.86 across regions. Uncertainty analysis on epidemiological parameters led to [Formula: see text] ranging from 1.57 to 2.34 on average over the metropolitan French regions, weighted by population size. CONCLUSIONS: Upon introduction, omicron spread rapidly through the French territory and showed a high fitness relative to delta. We documented considerable geographical heterogeneities on the spreading dynamics. The historical reconstruction of variant emergence dynamics provide valuable ground knowledge to face future variant emergence events.

Prioritizing COVID-19 vaccination efforts and dose allocation within Madagascar.

BACKGROUND: While mass COVID-19 vaccination programs are underway in high-income countries, limited availability of doses has resulted in few vaccines administered in low and middle income countries (LMICs). The COVID-19 Vaccines Global Access (COVAX) is a WHO-led initiative to promote vaccine access equity to LMICs and is providing many of the doses available in these settings. However, initial doses are limited and countries, such as Madagascar, need to develop prioritization schemes to maximize the benefits of vaccination with very limited supplies. There is some consensus that dose deployment should initially target health care workers, and those who are more vulnerable including older individuals. However, questions of geographic deployment remain, in particular associated with limits around vaccine access and delivery capacity in underserved communities, for example in rural areas that may also include substantial proportions of the population. METHODS: To address these questions, we developed a mathematical model of SARS-CoV-2 transmission dynamics and simulated various vaccination allocation strategies for Madagascar. Simulated strategies were based on a number of possible geographical prioritization schemes, testing sensitivity to initial susceptibility in the population, and evaluating the potential of tests for previous infection. RESULTS: Using cumulative deaths due to COVID-19 as the main outcome of interest, our results indicate that distributing the number of vaccine doses according to the number of elderly living in the region or according to the population size results in a greater reduction of mortality compared to distributing doses based on the reported number of cases and deaths. The benefits of vaccination strategies are diminished if the burden (and thus accumulated immunity) has been greatest in the most populous regions, but the overall strategy ranking remains comparable. If rapid tests for prior immunity may be swiftly and effectively delivered, there is potential for considerable gain in mortality averted, but considering delivery limitations modulates this. CONCLUSION: At a subnational scale, our results support the strategy adopted by the COVAX initiative at a global scale.

Spatiotemporal patterns of dengue and Zika incidence during the 2015-2018 outbreak of Zika in Mexico.

OBJECTIVE: Evaluate spatially and temporally simultaneous presence of clusters of dengue and Zika clinical cases and their relationship with expected dengue transmission risk. MATERIALS AND METHODS: A classification of dengue risk transmission was carried out for whole country, and spatial autocorrelation analyses to identify clusters of confirmed clinical cases of dengue and Zika from 2015 to 2018 was conducted using Moran's Index statistics. RESULTS: Clusters of both diseases were identified in dengue-high risk munici-palities at the beginning of the outbreak, but, at the end of the outbreak, Zika clusters occurred in dengue low-risk mu-nicipalities. CONCLUSION: This study identified Zika clusters in low-risk dengue areas suggesting participation of several factors that favor virus introduction and dissemination, such as differences in entomological and control interventions, and the possibility of cross-immunity in the population.

A preoperative nomogram for predicting long-term survival after resection of large hepatocellular carcinoma (>10 cm).

BACKGROUND: It has previously been demonstrated that a fraction of patients with hepatocellular carcinoma (HCC) > 10 cm can benefit from liver resection. However, there is still a lack of effective decision-making tools to inform intervention in these patients. METHODS: We analysed a comprehensive set of clinical data from 234 patients who underwent liver resection for HCC >10 cm at the National Cancer Institute of Peru between 1990 and 2015, monitored their survival, and constructed a nomogram to predict the surgical outcome based on preoperative variables. RESULTS: We identified cirrhosis, multifocality, macroscopic vascular invasion, and spontaneous tumour rupture as independent predictors of survival and integrated them into a nomogram model. The nomogram's ability to forecast survival at 1, 3, and 5 years was subsequently confirmed with high concordance using an internal validation. Through applying this nomogram, we stratified three groups of patients with different survival probabilities. CONCLUSION: We constructed a preoperative nomogram to predict long-term survival in patients with HCC >10 cm. This nomogram is useful in determining whether a patient with large HCC might truly benefit from liver resection, which is paramount in low- and middle-income countries where HCC is often diagnosed at advanced stages.

Is adaptive therapy natural?

Research suggests that progression-free survival can be prolonged by integrating evolutionary principles into clinical cancer treatment protocols. The goal is to prevent or slow the proliferation of resistant malignant cell populations. The logic behind this therapy relies on ecological and evolutionary processes. These same processes would be available to natural selection in decreasing the probability of an organism's death due to cancer. We propose that organisms' anticancer adaptions include not only ones for preventing cancer but also ones for directing and retarding the evolution of life-threatening cancer cells. We term this last strategy natural adaptive therapy (NAT). The body's NAT might include a lower than otherwise possible immune response. A restrained immune response might forego maximum short-term kill rates. Restraint would forestall immune-resistant cancer cells and produce long-term durable control of the cancer population. Here, we define, develop, and explore the possibility of NAT. The discovery of NAT mechanisms could identify new strategies in tumor prevention and treatments. Furthermore, we discuss the potential risks of immunotherapies that force the immune system to ramp up the short-term kill rates of malignant cancer cells in a manner that undermines the body's NAT and accelerates the evolution of immune resistance.

Was the COVID-19 pandemic avoidable? A call for a "solution-oriented" approach in pathogen evolutionary ecology to prevent future outbreaks.

Concerns about the prospect of a global pandemic have been triggered many times during the last two decades. These have been realised through the current COVID-19 pandemic, due to a new coronavirus SARS-CoV2, which has impacted almost every country on Earth. Here, we show how considering the pandemic through the lenses of the evolutionary ecology of pathogens can help better understand the root causes and devise solutions to prevent the emergence of future pandemics. We call for better integration of these approaches into transdisciplinary research and invite scientists working on the evolutionary ecology of pathogens to contribute to a more "solution-oriented" agenda with practical applications, emulating similar movements in the field of economics in recent decades.

Exposure to Ebola Virus and Risk for Infection with Malaria Parasites, Rural Gabon.

An association between malaria and risk for death among patients with Ebola virus disease has suggested within-host interactions between Plasmodium falciparum parasites and Ebola virus. To determine whether such an interaction might also influence the probability of acquiring either infection, we used a large snapshot surveillance study from rural Gabon to test if past exposure to Ebola virus is associated with current infection with Plasmodium spp. during nonepidemic conditions. We found a strong positive association, on population and individual levels, between seropositivity for antibodies against Ebola virus and the presence of Plasmodium parasites in the blood. According to a multiple regression model accounting for other key variables, antibodies against Ebola virus emerged as the strongest individual-level risk factor for acquiring malaria. Our results suggest that within-host interactions between malaria parasites and Ebola virus may underlie epidemiologic associations.

The evolution of resistance and tolerance as cancer defences.

Although there is a plethora of cancer associated-factors that can ultimately culminate in death (cachexia, organ impairment, metastases, opportunistic infections, etc.), the focal element of every terminal malignancy is the failure of our natural defences to control unlimited cell proliferation. The reasons why our defences apparently lack efficiency is a complex question, potentially indicating that, under Darwinian terms, solutions other than preventing cancer progression are also important contributors. In analogy with host-parasite systems, we propose to call this latter option 'tolerance' to cancer. Here, we argue that the ubiquity of oncogenic processes among metazoans is at least partially attributable to both the limitations of resistance mechanisms and to the evolution of tolerance to cancer. Deciphering the ecological contexts of alternative responses to the cancer burden is not a semantic question, but rather a focal point in understanding the evolutionary ecology of host-tumour relationships, the evolution of our defences, as well as why and when certain cancers are likely to be detrimental for survival.

The Kendrick modelling platform: language abstractions and tools for epidemiology.

BACKGROUND: Mathematical and computational models are widely used to study the transmission, pathogenicity, and propagation of infectious diseases. Unfortunately, complex mathematical models are difficult to define, reuse and reproduce because they are composed of several concerns that are intertwined. The problem is even worse for computational models because the epidemiological concerns are also intertwined with low-level implementation details that are not easily accessible to non-computing scientists. Our goal is to make compartmental epidemiological models easier to define, reuse and reproduce by facilitating implementation of different simulation approaches with only very little programming knowledge. RESULTS: We achieve our goal through the definition of a domain-specific language (DSL), Kendrick, that relies on a very general mathematical definition of epidemiological concerns as stochastic automata that are combined using tensor-algebra operators. A very large class of epidemiological concerns, including multi-species, spatial concerns, control policies, sex or age structures, are supported and can be defined independently of each other and combined into models to be simulated by different methods. Implementing models does not require sophisticated programming skills any more. The various concerns involved within a model can be changed independently of the others as well as reused within other models. They are not plagued by low-level implementation details. CONCLUSIONS: Kendrick is one of the few DSLs for epidemiological modelling that does not burden its users with implementation details or required sophisticated programming skills. It is also currently the only language for epidemiology modelling that supports modularity through clear separation of concerns hence fostering reproducibility and reuse of models and simulations. Future work includes extending Kendrick to support non-compartmental models and improving its interoperability with existing complementary tools.