Worldwide spread of Dengue virus type 1

BACKGROUND: DENV-1 is one of the four viral serotypes that causes Dengue, the most common mosquito-borne viral disease of humans. The prevalence of these viruses has grown in recent decades and is now present in more than 100 countries. Limited studies document the spread of DENV-1 over the world despite its importance for human health. METHODOLOGY/PRINCIPAL FINDINGS: We used representative DENV-1 envelope gene sequences to unravel the dynamics of viral diffusion under a Bayesian phylogeographic approach. Data included strains from 45 distinct geographic locations isolated from 1944 to 2009. The estimated mean rate of nucleotide substitution was 6.56 × 10⁻⁴ substitutions/site/year. The larger genotypes (I, IV and V) had a distinctive phylogenetic structure and since 1990 they experienced effective population size oscillations. Thailand and Indonesia represented the main sources of strains for neighboring countries. Besides, Asia broadcast lineages into the Americas and the Pacific region that diverged in isolation. Also, a transmission network analysis revealed the pivotal role of Indochina in the global diffusion of DENV-1 and of the Caribbean in the diffusion over the Americas. CONCLUSIONS/SIGNIFICANCE: The study summarizes the spatiotemporal DENV-1 worldwide spread that may help disease control.

Bicistronic DNA vaccines simultaneously encoding HIV, HSV and HPV antigens promote CD8⁺ T cell responses and protective immunity

Millions of people worldwide are currently infected with human papillomavirus (HPV), herpes simplex virus (HSV) or human immunodeficiency virus (HIV). For this enormous contingent of people, the search for preventive and therapeutic immunological approaches represents a hope for the eradication of latent infection and/or virus-associated cancer. To date, attempts to develop vaccines against these viruses have been mainly based on a monovalent concept, in which one or more antigens of a virus are incorporated into a vaccine formulation. In the present report, we designed and tested an immunization strategy based on DNA vaccines that simultaneously encode antigens for HIV, HSV and HPV. With this purpose in mind, we tested two bicistronic DNA vaccines (pIRES I and pIRES II) that encode the HPV-16 oncoprotein E7 and the HIV protein p24 both genetically fused to the HSV-1 gD envelope protein. Mice i.m. immunized with the DNA vaccines mounted antigen-specific CD8⁺ T cell responses, including in vivo cytotoxic responses, against the three antigens. Under experimental conditions, the vaccines conferred protective immunity against challenges with a vaccinia virus expressing the HIV-derived protein Gag, an HSV-1 virus strain and implantation of tumor cells expressing the HPV-16 oncoproteins. Altogether, our results show that the concept of a trivalent HIV, HSV, and HPV vaccine capable to induce CD8⁺ T cell-dependent responses is feasible and may aid in the development of preventive and/or therapeutic approaches for the control of diseases associated with these viruses.

The dengue virus non-structural 1 protein: risks and benefits

The dengue virus (DENV) non-structural 1 (NS1) protein plays a critical role in viral RNA replication and has a central position in DENV pathogenesis. DENV NS1 is a glycoprotein expressed in infected mammalian cells as soluble monomers that dimerize in the lumen of the endoplasmic reticulum; NS1 is subsequently transported to the cell surface, where it remains membrane associated or is secreted into the extracellular milieu as a hexameric complex. During the last three decades, the DENV NS1 protein has also been intensively investigated as a potential target for vaccines and antiviral drugs. In addition, NS1 is the major diagnostic marker for dengue infection. This review highlights some important issues regarding the role of NS1 in DENV pathogenesis and its biotechnological applications, both as a target for the development of safe and effective vaccines and antiviral drugs and as a tool for the generation of accurate diagnostic methods

Meccanismi evolutivi dei parapoxvirus: caratterizzazione genomica di pseudocowpoxvirus e messa a punto di sistemi per lo studio delle ricombinazioni

The Poxviruses are a family of double stranded DNA (dsDNA) viruses that cause disease in many species, both vertebrate and invertebrate. Their genomes range in size from 135 to 365 kbp and show conservation in both organization and content. In particular, the central genomic regions of the chordopoxvirus subfamily (those capable of infecting vertebrates) contain 88 genes which are present in all the virus species characterised to date and which mostly occur in the same order and orientation. In contrast, however, the terminal regions of the genomes frequently contain genes that are species or genera-specific and that are not essential for the growth of the virus in vitro but instead often encode factors with important roles in vivo including modulation of the host immune response to infection and determination of the host range of the virus. The Parapoxviruses (PPV), of which Orf virus is the prototypic species, represent a genus within the chordopoxvirus subfamily of Poxviridae and are characterised by their ability to infect ruminants and humans. The genus currently contains four recognised species of virus, bovine papular stomatitis virus (BPSV) and pseudocowpox virus (PCPV) both of which infect cattle, orf virus (OV) that infects sheep and goats, and parapoxvirus of red deer in New Zealand (PVNZ). The ORFV genome has been fully sequenced, as has that of BPSV, and is ~138 kb in length encoding ~132 genes. The vast majority of these genes allow the virus to replicate in the cytoplasm of the infected host cell and therefore encode proteins involved in replication, transcription and metabolism of nucleic acids. These genes are well conserved between all known genera of poxviruses. There is however another class of genes, located at either end of the linear dsDNA genome, that encode proteins which are non-essential for replication and generally dictate host range and virulence of the virus. The non-essential genes are often the most variable within and between species of virus and therefore are potentially useful for diagnostic purposes. Given their role in subverting the host-immune response to infection they are also targets for novel therapeutics. The function of only a relatively small number of these proteins has been elucidated and there are several genes whose function still remains obscure principally because there is little similarity between them and proteins of known function in current sequence databases. It is thought that by selectively removing some of the virulence genes, or at least neutralising the proteins in some way, current vaccines could be improved. The evolution of poxviruses has been proposed to be an adaptive process involving frequent events of gene gain and loss, such that the virus co-evolves with its specific host. Gene capture or horizontal gene transfer from the host to the virus is considered an important source of new viral genes including those likely to be involved in host range and those enabling the virus to interfere with the host immune response to infection. Given the low rate of nucleotide substitution, recombination can be seen as an essential evolutionary driving force although it is likely underestimated. Recombination in poxviruses is intimately linked to DNA replication with both viral and cellular proteins participate in this recombination-dependent replication. It has been shown, in other poxvirus genera, that recombination between isolates and perhaps even between species does occur, thereby providing another mechanism for the acquisition of new genes and for the rapid evolution of viruses. Such events may result in viruses that have a selective advantage over others, for example in re-infections (a characteristic of the PPV), or in viruses that are able to jump the species barrier and infect new hosts. Sequence data related to viral strains isolated from goats suggest that possible recombination events may have occurred between OV and PCPV (Ueda et al. 2003). The recombination events are frequent during poxvirus replication and comparative genomic analysis of several poxvirus species has revealed that recombinations occur frequently on the right terminal region. Intraspecific recombination can occur between strains of the same PPV species, but also interspecific recombination can happen depending on enough sequence similarity to enable recombination between distinct PPV species. The most important pre-requisite for a successful recombination is the coinfection of the individual host by different virus strains or species. Consequently, the following factors affecting the distribution of different viruses to shared target cells need to be considered: dose of inoculated virus, time interval between inoculation of the first and the second virus, distance between the marker mutations, genetic homology. At present there are no available data on the replication dynamics of PPV in permissive and non permissive hosts and reguarding co-infetions there are no information on the interference mechanisms occurring during the simultaneous replication of viruses of different species. This work has been carried out to set up permissive substrates allowing the replication of different PPV species, in particular keratinocytes monolayers and organotypic skin cultures. Furthermore a method to isolate and expand ovine skin stem cells was has been set up to indeep further aspects of viral cellular tropism during natural infection. The study produced important data to elucidate the replication dynamics of OV and PCPV virus in vitro as well as the mechanisms of interference that can arise during co-infection with different viral species. Moreover, the analysis carried on the genomic right terminal region of PCPV 1303/05 contributed to a better knowledge of the viral genes involved in host interaction and pathogenesis as well as to locate recombination breakpoints and genetic homologies between PPV species. Taken together these data filled several crucial gaps for the study of interspecific recombinations of PPVs which are thought to be important for a better understanding of the viral evolution and to improve the biosafety of antiviral therapy and PPV-based vectors.

Evaluation of glycoconjugate antigens as vaccine candidates against group A Streptococcus and human immunodeficiency virus infections

This PhD thesis discusses the rationale for design and use of synthetic oligosaccharides for the development of glycoconjugate vaccines and the role of physicochemical methods in the characterization of these vaccines. The study concerns two infectious diseases that represent a serious problem for the national healthcare programs: human immunodeficiency virus (HIV) and Group A Streptococcus (GAS) infections. Both pathogens possess distinctive carbohydrate structures that have been described as suitable targets for the vaccine design. The Group A Streptococcus cell membrane polysaccharide (GAS-PS) is an attractive vaccine antigen candidate based on its conserved, constant expression pattern and the ability to confer immunoprotection in a relevant mouse model. Analysis of the immunogenic response within at-risk populations suggests an inverse correlation between high anti-GAS-PS antibody titres and GAS infection cases. Recent studies show that a chemically synthesized core polysaccharide-based antigen may represent an antigenic structural determinant of the large polysaccharide. Based on GAS-PS structural analysis, the study evaluates the potential to exploit a synthetic design approach to GAS vaccine development and compares the efficiency of synthetic antigens with the long isolated GAS polysaccharide. Synthetic GAS-PS structural analogues were specifically designed and generated to explore the impact of antigen length and terminal residue composition. For the HIV-1 glycoantigens, the dense glycan shield on the surface of the envelope protein gp120 was chosen as a target. This shield masks conserved protein epitopes and facilitates virus spread via binding to glycan receptors on susceptible host cells. The broadly neutralizing monoclonal antibody 2G12 binds a cluster of high-mannose oligosaccharides on the gp120 subunit of HIV-1 Env protein. This oligomannose epitope has been a subject to the synthetic vaccine development. The cluster nature of the 2G12 epitope suggested that multivalent antigen presentation was important to develop a carbohydrate based vaccine candidate. I describe the development of neoglycoconjugates displaying clustered HIV-1 related oligomannose carbohydrates and their immunogenic properties.

Reverse genetic studies of Benyvirus - Polymyxa betae molecular interaction: Role of the RNA4-encoded protein in virus transmission

Beet necrotic yellow vein virus (BNYVV), the leading infectious agent that affects sugar beet, is included within viruses transmitted through the soil from plasmodiophorid as Polymyxa betae. BNYVV is the causal agent of Rhizomania, which induces abnormal rootlet proliferation and is widespread in the sugar beet growing areas in Europe, Asia and America; for review see (Peltier et al., 2008). In this latter continent, Beet soil-borne mosaic virus (BSBMV) has been identified (Lee et al., 2001) and belongs to the benyvirus genus together with BNYVV, both vectored by P. betae. BSBMV is widely distributed only in the United States and it has not been reported yet in others countries. It was first identified in Texas as a sugar beet virus morphologically similar but serologically distinct to BNYVV. Subsequent sequence analysis of BSBMV RNAs evidenced similar genomic organization to that of BNYVV but sufficient molecular differences to distinct BSBMV and BNYVV in two different species (Rush et al., 2003). Benyviruses field isolates usually consist of four RNA species but some BNYVV isolates contain a fifth RNA. RNAs -1 contains a single long ORF encoding polypeptide that shares amino acid homology with known viral RNA-dependent RNA polymerases (RdRp) and helicases. RNAs -2 contains six ORFs: capsid protein (CP), one readthrough protein, triple gene block proteins (TGB) that are required for cell-to-cell virus movement and the sixth 14 kDa ORF is a post-translation gene silencing suppressor. RNAs -3 is involved on disease symptoms and is essential for virus systemic movement. BSBMV RNA-3 can be trans-replicated, trans-encapsidated by the BNYVV helper strain (RNA-1 and -2) (Ratti et al., 2009). BNYVV RNA-4 encoded one 31 kDa protein and is essential for vector interactions and virus transmission by P. betae (Rahim et al., 2007). BNYVV RNA-5 encoded 26 kDa protein that improve virus infections and accumulation in the hosts. We are interest on BSBMV effect on Rhizomania studies using powerful tools as full-length infectious cDNA clones. B-type full-length infectious cDNA clones are available (Quillet et al., 1989) as well as A/P-type RNA-3, -4 and -5 from BNYVV (unpublished). A-type BNYVV full-length clones are also available, but RNA-1 cDNA clone still need to be modified. During the PhD program, we start production of BSBMV full-length cDNA clones and we investigate molecular interactions between plant and Benyviruses exploiting biological, epidemiological and molecular similarities/divergences between BSBMV and BNYVV. During my PhD researchrs we obtained full length infectious cDNA clones of BSBMV RNA-1 and -2 and we demonstrate that they transcripts are replicated and packaged in planta and able to substitute BNYVV RNA-1 or RNA-2 in a chimeric viral progeny (BSBMV RNA-1 + BNYVV RNA-2 or BNYVV RNA-1 + BSBMV RNA-2). During BSBMV full-length cDNA clones production, unexpected 1,730 nts long form of BSBMV RNA-4 has been detected from sugar beet roots grown on BSBMV infected soil. Sequence analysis of the new BSBMV RNA-4 form revealed high identity (~100%) with published version of BSBMV RNA-4 sequence (NC_003508) between nucleotides 1-608 and 1,138-1,730, however the new form shows 528 additionally nucleotides between positions 608-1,138 (FJ424610). Two putative ORFs has been identified, the first one (nucleotides 383 to 1,234), encode a protein with predicted mass of 32 kDa (p32) and the second one (nucleotides 885 to 1,244) express an expected product of 13 kDa (p13). As for BSBMV RNA-3 (Ratti et al., 2009), full-length BSBMV RNA-4 cDNA clone permitted to obtain infectious transcripts that BNYVV viral machinery (Stras12) is able to replicate and to encapsidate in planta. Moreover, we demonstrated that BSBMV RNA-4 can substitute BNYVV RNA-4 for an efficient transmission through the vector P. betae in Beta vulgaris plants, demonstrating a very high correlation between BNYVV and BSBMV. At the same time, using BNYVV helper strain, we studied BSBMV RNA-4’s protein expression in planta. We associated a local necrotic lesions phenotype to the p32 protein expression onto mechanically inoculated C. quinoa. Flag or GFP-tagged sequences of p32 and p13 have been expressed in viral context, using Rep3 replicons, based on BNYVV RNA-3. Western blot analyses of local lesions contents, using FLAG-specific antibody, revealed a high molecular weight protein, which suggest either a strong interaction of BSBMV RNA4’s protein with host protein(s) or post translational modifications. GFP-fusion sequences permitted the subcellular localization of BSBMV RNA4’s proteins. Moreover we demonstrated the absence of self-activation domains on p32 by yeast two hybrid system approaches. We also confirmed that p32 protein is essential for virus transmission by P. betae using BNYVV helper strain and BNYVV RNA-3 and we investigated its role by the use of different deleted forms of p32 protein. Serial mechanical inoculation of wild-type BSBMV on C. quinoa plants were performed every 7 days. Deleted form of BSBMV RNA-4 (1298 bp) appeared after 14 passages and its sequence analysis shows deletion of 433 nucleotides between positions 611 and 1044 of RNA-4 new form. We demonstrated that this deleted form can’t support transmission by P. betae using BNYVV helper strain and BNYVV RNA-3, moreover we confirmed our hypothesis that BSBMV RNA-4 described by Lee et al. (2001) is a deleted form. Interesting after 21 passages we identifed one chimeric form of BSBMV RNA-4 and BSBMV RNA-3 (1146 bp). Two putative ORFs has been identified on its sequence, the first one (nucleotides 383 to 562), encode a protein with predicted mass of 7 kDa (p7), corresponding to the N-terminal of p32 protein encoded by BSBMV RNA-4; the second one (nucleotides 562 to 789) express an expected product of 9 kDa (p9) corresponding to the C-terminal of p29 encoded by BSBMV RNA-3. Results obtained by our research in this topic opened new research lines that our laboratories will develop in a closely future. In particular BSBMV p32 and its mutated forms will be used to identify factors, as host or vector protein(s), involved in the virus transmission through P. betae. The new results could allow selection or production of sugar beet plants able to prevent virus transmission then able to reduce viral inoculum in the soil.

Epidemiologia e caratterizzazione molecolare di virus emergenti responsabili di Zoonosi: Calcivirus e Virus dell'Epatite E

Il virus dell’Epatite E (HEV) e i calicivirus (norovirus e sapovirus) causano rispettivamente epatite acuta e gastroenterite. Questi virus sono considerati agenti eziologici emergenti rappresentando un problema di sanità pubblica e di sicurezza alimentare. Per HEV, è ormai confermata la trasmissione zoonotica, e il suino è considerato il principale serbatoio asintomatico. Norovirus e sapovirus infettano sia i bambini che gli adulti. Sebbene questi virus siano stati identificati anche negli animali, la possibile trasmissione zoonotica non è stata dimostrata in modo conclusivo. Il lavoro sperimentale condotto durante il Dottorato di Ricerca è stato focalizzato sullo studio degli aspetti biologici ed epidemiologici dell’infezioni causate da HEV e da calicivirus. Per la prima volta in Italia, i risultati ottenuti hanno dimostrato la presenza del virus HEV nei fegati di suini in fase di macellazione ed hanno confermato, attraverso la ricerca di anticorpi, un’elevata esposizione degli animali al virus. Inoltre, mediante la produzione di antigeni e reattivi immunologici, sono stati messi a punto test diagnostici per la ricerca di anticorpi contro HEV nel suino e nei cinghiali. Il lavoro svolto per la ricerca di calicivirus nel suino e nel bovino ha dimostrato la circolazione dei sapovirus in popolazioni di suini asintomatici e la presenza di norovirus nei vitelli affetti da diarrea acuta.Sono stati inoltre sviluppati reattivi immunologici, utilizzando proteine del capside di norovirus umano e bovino espresse con il sistema ricombinante baculovirus. Questi hanno permesso di evidenziare la presenza di anticorpi contro norovirus umano e bovino, in sieri di veterinari professionalmente esposti. Inoltre, sono stati utilizzati per sviluppare metodi per la concentrazione dei virus da matrici a bassa concentrazione.Infine, le VLP sono state utilizzate per valutare l’attivazione del sistema immunitario umano ex vivo. I risultati hanno dimostrato che le VLP di NoV stimolano il sistema immunitario attivando risposte di tipo Th1 e Th2 .

Development of new therapeutic approaches in the treatments of Inflammatory Bowel Diseases: functional food and nutraceuticals vs synthetic peptides based vaccines

Inflammatory Bowel Diseases (IBD) are intestinal chronic relapsing diseases which ethiopathogenesis remains uncertain. Several group have attempted to study the role of factors involved such as genetic susceptibility, environmental factors such as smoke, diet, sex, immunological factors as well as the microbioma. None of the treatments available satisfy several criteria at the same time such as safety, long-term remission, histopatological healing, and specificity. We used two different approaches for the development of new therapeutic treatment for Inflammatory Bowel Disease. The first is focused on the understanding of the potential role of functional food and nutraceuticals nutrients in the treatment of IBD. To do so, we investigated the role of Curcuma longa in the treatment of chemical induced colitis in mice model. Since Curcma Longa has been investigated for its antinflammatory role related to the TNFα pathway as well investigators have reported few cases of patients with ulcerative colites treated with this herbs, we harbored the hypothesis of a role of Curcuma Longa in the treatment f IBD as well as we decided to assess its role in intestinal motility. The second part is based on an immunological approach to develop new drugs to induce suppression in Crohn’s disease or to induce mucosa immunity such as in colonrectal tumor. The main idea behind this approach is that we could manipulate relevant cell-cell interactions using synthetic peptides. We demonstrated the role of the unique interaction between molecules expressed on intestinal epithelial cells such as CD1d and CEACAM5 and on CD8+ T cells. In normal condition this interaction has a role for the expansion of the suppressor CD8+ T cells. Here, we characterized this interaction, we defined which are the epitope involved in the binding and we attempted to develop synthetic peptides from the N domain of CEACAM5 in order to manipulate it.

Functional characterization of modified vaccinia virus Ankara-encoded anti-apoptotic proteins

Candidate vaccines based on the highly attenuated orthopoxvirus strain MVA are tested against various infectious and cancer diseases and, more profound, vaccines based on wildtype and recombinant viruses have been found safe and immunogenic in clinical trials. Compared to conventional vaccine strains, MVA lacks many functional genes for potentially important regulators of virus-host interactions. However, some gene functions responsible for counteraction of cellular antiviral pathways are still conserved in the genome of MVA and the inhibition of apoptosis seems to be one important mechanism, the virus is still able to interact with.rnrnVaccinia viruses encode several proteins which prevent the induction of virus-induced apoptosis. The vaccinia virus anti-apoptotic protein F1 was shown to counteract the activation of the mitochondrial pathway of apoptosis in a highly effective manner. Another vaccinia virus protein, N1, like F1 shows structural and functional similarity to members of the cellular anti-apoptotic bcl-2 family and was also shown to inhibit apoptosis. The vaccinia virus early protein E3 inhibits programmed cell death by binding to and sequestration of dsRNA molecules, normally inducing cellular antiviral pathways also driving the induction of apoptosis. All three anti-apoptotic genes were functionally analyzed during this work.rn

Sorveglianza dell’infezione da virus dell’epatite E (HEV) in Italia: from farm to table

L'epatite E è una malattia umana con caratteristiche di epatite acuta, causata da un ssRNA virus (HEV). Nel 1997, HEV è stato identificato per la prima volta nei suini (SwHEV). In seguito, diverse evidenze, tra cui la vicinanza genetica tra ceppi umani e suini, suggerirono la trasmissione zoonotica del virus. Nella presente tesi, l’identificazione di SwHEV è stata condotta mediante ricerca di porzioni di genoma virale attraverso RT-PCR. Dal 2011 al 2013, sono stati analizzati 343 campioni fecali (da 19 allevamenti) e 70 bili (da 2 macelli) prelevati da altrettanti suini, in diverse Regioni italiane. E’ stato inoltre condotto uno studio retrospettivo su 78 feci (da 3 allevamenti) raccolte nel 2000. Il virus è stato identificato nel 24,5% e 19,2% delle feci raccolte rispettivamente nel 2011-2013 e nel 2000. Nessuna bile è risultata positiva. Mediante sequenziamento del genoma intero di uno dei virus identificati, è stata condotta l’analisi filogenetica per valutarne il grado di correlazione con alti ceppi suini e umani. La presenza di HEV è stata valutata lungo la filiera di produzione suina, dal macello al punto vendita. Trentaquattro campioni di feci, fegato e muscolo sono stati raccolti in un macello da altrettanti suini sani (età:6-7 mesi). Quattordici feci e 2 fegati, sono risultati positivi per HEV. Sono state prelevate 129 salsicce sia allo stabilimento di trasformazione sia alla vendita, ma nessuna è risultata positiva. La presenza di HEV è stata valutata anche nelle salsicce di fegato, fresche e secche, acquistate presso una macelleria. Il genoma virale è stato rilevato nel 22,2% delle salsicce fresche e nel 4,3 % di quelle secche ma la vitalità del virus non è stata dimostrata. In conclusione, lo studio condotto ha confermato l’ampia circolazione di HEV nei suini e la possibile contaminazione dei prodotti carnei derivati, confermando la necessità di una continua sorveglianza.

Molecular variability of Meningococcal antigens in carriage and disease strains and in other Neisseria species

This PhD thesis is focused on the study of the molecular variability of some specific proteins, part of the outer membrane of the pathogen Neisseria meningitidis, and described as protective antigens and important virulence factors. These antigens have been employed as components of the vaccine developed by Novartis Vaccines against N. meningitidis of serogroup B, and their variability in the meningococcal population is a key aspect when the effect of the vaccine is evaluated. The PhD project has led to complete three major studies described in three different manuscritps, of which two have been published and the third is in preparation. The thesis is structured in three main chapters, each of them dedicated to the three studies. The first, described in Chapter 1, is specifically dedicated to the analysis of the molecular conservation of meningococcal antigens in the genomes of all species classified in the genus Neisseria (Conservation of Meningococcal Antigens in the Genus Neisseria. A. Muzzi et al.. 2013. mBio 4 (3)). The second study, described in Chapter 2, focuses on the analysis of the presence and conservation of the antigens in a panel of bacterial isolates obtained from cases of the disease and from healthy individuals, and collected in the same year and in the same geographical area (Conservation of fHbp, NadA, and NHBA in carrier and pathogenic isolates of Neisseria meningitidis collected in the Czech Republic in 1993. A. Muzzi et al.. Manuscript in preparation). Finally, Chapter 3 describes the molecular features of the antigens in a panel of bacterial isolates collected over a period of 50 years, and representatives of the epidemiological history of meningococcal disease in the Netherlands (An Analysis of the Sequence Variability of Meningococcal fHbp, NadA and NHBA over a 50-Year Period in the Netherlands. S. Bambini et al.. 2013. PloS one e65043).