Your browser doesn't support javascript.
loading
Show: 20 | 50 | 100
Results 1 - 20 de 106
Filter
1.
J Virol ; : e0008324, 2024 Jul 12.
Article in English | MEDLINE | ID: mdl-38995042

ABSTRACT

Mosquitoes can transmit several pathogenic viruses to humans, but their natural viral community is also composed of a myriad of other viruses such as insect-specific viruses (ISVs) and those that infect symbiotic microorganisms. Besides a growing number of studies investigating the mosquito virome, the majority are focused on few urban species, and relatively little is known about the virome of sylvatic mosquitoes, particularly in high biodiverse biomes such as the Brazilian biomes. Here, we characterized the RNA virome of 10 sylvatic mosquito species from Atlantic forest remains at a sylvatic-urban interface in Northeast Brazil employing a metatranscriptomic approach. A total of 16 viral families were detected. The phylogenetic reconstructions of 14 viral families revealed that the majority of the sequences are putative ISVs. The phylogenetic positioning and, in most cases, the association with a high RNA-dependent RNA polymerase amino acid divergence from other known viruses suggests that the viruses characterized here represent at least 34 new viral species. Therefore, the sylvatic mosquito viral community is predominantly composed of highly divergent viruses highlighting the limited knowledge we still have about the natural virome of mosquitoes in general. Moreover, we found that none of the viruses recovered were shared between the species investigated, and only one showed high identity to a virus detected in a mosquito sampled in Peru, South America. These findings add further in-depth understanding about the interactions and coevolution between mosquitoes and viruses in natural environments. IMPORTANCE: Mosquitoes are medically important insects as they transmit pathogenic viruses to humans and animals during blood feeding. However, their natural microbiota is also composed of a diverse set of viruses that cause no harm to the insect and other hosts, such as insect-specific viruses. In this study, we characterized the RNA virome of sylvatic mosquitoes from Northeast Brazil using unbiased metatranscriptomic sequencing and in-depth bioinformatic approaches. Our analysis revealed that these mosquitoes species harbor a diverse set of highly divergent viruses, and the majority comprises new viral species. Our findings revealed many new virus lineages characterized for the first time broadening our understanding about the natural interaction between mosquitoes and viruses. Finally, it also provided several complete genomes that warrant further assessment for mosquito and vertebrate host pathogenicity and their potential interference with pathogenic arboviruses.

2.
An Acad Bras Cienc ; 96(2): e20230894, 2024.
Article in English | MEDLINE | ID: mdl-38922277

ABSTRACT

The need for the identification of risk factors associated to COVID-19 disease severity remains urgent. Patients' care and resource allocation can be potentially different and are defined based on the current classification of disease severity. This classification is based on the analysis of clinical parameters and routine blood tests, which are not standardized across the globe. Some laboratory test alterations have been associated to COVID-19 severity, although these data are conflicting partly due to the different methodologies used across different studies. This study aimed to construct and validate a disease severity prediction model using machine learning (ML). Seventy-two patients admitted to a Brazilian hospital and diagnosed with COVID-19 through RT-PCR and/or ELISA, and with varying degrees of disease severity, were included in the study. Their electronic medical records and the results from daily blood tests were used to develop a ML model to predict disease severity. Using the above data set, a combination of five laboratorial biomarkers was identified as accurate predictors of COVID-19 severe disease with a ROC-AUC of 0.80 ​±â€‹ 0.13. Those biomarkers included prothrombin activity, ferritin, serum iron, ATTP and monocytes. The application of the devised ML model may help rationalize clinical decision and care.


Subject(s)
Biomarkers , COVID-19 , Machine Learning , SARS-CoV-2 , Severity of Illness Index , Humans , COVID-19/blood , COVID-19/diagnosis , Female , Male , Biomarkers/blood , Middle Aged , Prognosis , Adult , Ferritins/blood , Aged , Brazil , Hematologic Tests/methods , ROC Curve , Risk Factors
3.
NAR Genom Bioinform ; 6(2): lqae056, 2024 Jun.
Article in English | MEDLINE | ID: mdl-38800829

ABSTRACT

ViralFlow v1.0 is a computational workflow developed for viral genomic surveillance. Several key changes turned ViralFlow into a general-purpose reference-based genome assembler for all viruses with an available reference genome. New virus-agnostic modules were implemented to further study nucleotide and amino acid mutations. ViralFlow v1.0 runs on a broad range of computational infrastructures, from laptop computers to high-performance computing (HPC) environments, and generates standard and well-formatted outputs suited for both public health reporting and scientific problem-solving. ViralFlow v1.0 is available at: https://viralflow.github.io/index-en.html.

4.
medRxiv ; 2024 May 17.
Article in English | MEDLINE | ID: mdl-38798319

ABSTRACT

Dengue virus (DENV) is currently causing epidemics of unprecedented scope in endemic settings and expanding to new geographical areas. It is therefore critical to track this virus using genomic surveillance. However, the complex patterns of viral genomic diversity make it challenging to use the existing genotype classification system. Here we propose adding two sub-genotypic levels of virus classification, named major and minor lineages. These lineages have high thresholds for phylogenetic distance and clade size, rendering them stable between phylogenetic studies. We present an assignment tool to show that the proposed lineages are useful for regional, national and sub-national discussions of relevant DENV diversity. Moreover, the proposed lineages are robust to classification using partial genome sequences. We provide a standardized neutral descriptor of DENV diversity with which we can identify and track lineages of potential epidemiological and/or clinical importance. Information about our lineage system, including methods to assign lineages to sequence data and propose new lineages, can be found at: dengue-lineages.org.

5.
Virus Evol ; 10(1): veae031, 2024.
Article in English | MEDLINE | ID: mdl-38756986

ABSTRACT

The highly pathogenic avian influenza viruses of clade 2.3.4.4b have caused unprecedented deaths in South American wild birds, poultry, and marine mammals. In September 2023, pinnipeds and seabirds appeared dead on the Uruguayan Atlantic coast. Sixteen influenza virus strains were characterized by real-time reverse transcription PCR and genome sequencing in samples from sea lions (Otaria flavescens), fur seals (Arctocephalus australis), and terns (Sterna hirundinacea). Phylogenetic and ancestral reconstruction analysis showed that these strains have pinnipeds most likely as the ancestral host, representing a recent introduction of clade 2.3.4.4b in Uruguay. The Uruguayan and closely related strains from Peru (sea lions) and Chile (sea lions and a human case) carry mammalian adaptative residues 591K and 701N in the viral polymerase basic protein 2 (PB2). Our findings suggest that clade 2.3.4.4b strains in South America may have spread from mammals to mammals and seabirds, revealing a new transmission route.

6.
medRxiv ; 2024 Apr 18.
Article in English | MEDLINE | ID: mdl-38699354

ABSTRACT

During the ongoing western equine encephalitis virus (WEEV) outbreak in South America, we described three fatal cases in horses from Rio Grande do Sul, Brazil. We sequenced WEEV strains and identified a novel lineage causing these cases. Continued surveillance and horse immunization are needed to mitigate the WEEV burden.

7.
J Mol Biol ; 436(11): 168577, 2024 Jun 01.
Article in English | MEDLINE | ID: mdl-38642883

ABSTRACT

The Red Queen Hypothesis (RQH), derived from Lewis Carroll's "Through the Looking-Glass", postulates that organisms must continually adapt in response to each other to maintain relative fitness. Within the context of host-pathogen interactions, the RQH implies an evolutionary arms race, wherein viruses evolve to exploit hosts and hosts evolve to resist viral invasion. This study delves into the dynamics of the RQH in the context of virus-cell interactions, specifically focusing on virus receptors and cell receptors. We observed multiple virus-host systems and noted patterns of co-evolution. As viruses evolved receptor-binding proteins to effectively engage with cell receptors, cells countered by altering their receptor genes. This ongoing mutual adaptation cycle has influenced the molecular intricacies of receptor-ligand interactions. Our data supports the RQH as a driving force behind the diversification and specialization of both viral and host cell receptors. Understanding this co-evolutionary dance offers insights into the unpredictability of emerging viral diseases and potential therapeutic interventions. Future research is crucial to dissect the nuanced molecular changes and the broader ecological consequences of this ever-evolving battle. Here, we combine phylogenetic inferences, structural modeling, and molecular dynamics analyses to describe the epidemiological characteristics of major Brazilian DENV strains that circulated from 1990 to 2022 from a combined perspective, thus providing us with a more detailed picture on the dynamics of such interactions over time.


Subject(s)
Cell Adhesion Molecules , Dengue Virus , Evolution, Molecular , Host-Pathogen Interactions , Receptors, Cell Surface , Viral Envelope Proteins , Viral Envelope , Humans , Brazil , Cell Adhesion Molecules/metabolism , Cell Adhesion Molecules/genetics , Cell Adhesion Molecules/chemistry , Dengue/virology , Dengue Virus/genetics , Dengue Virus/metabolism , Host-Pathogen Interactions/genetics , Lectins, C-Type/metabolism , Lectins, C-Type/genetics , Lectins, C-Type/chemistry , Molecular Dynamics Simulation , Phylogeny , Protein Binding , Receptors, Cell Surface/metabolism , Receptors, Cell Surface/genetics , Receptors, Cell Surface/chemistry , Receptors, Virus/metabolism , Receptors, Virus/chemistry , Receptors, Virus/genetics , Viral Envelope/metabolism , Viral Envelope Proteins/genetics , Viral Envelope Proteins/metabolism , Viral Envelope Proteins/chemistry
8.
Microbiol Spectr ; 12(6): e0421823, 2024 Jun 04.
Article in English | MEDLINE | ID: mdl-38651879

ABSTRACT

SARS-CoV-2 virus emerged as a new threat to humans and spread around the world, leaving a large death toll. As of January 2023, Brazil is among the countries with the highest number of registered deaths. Nonpharmacological and pharmacological interventions have been heterogeneously implemented in the country, which, associated with large socioeconomic differences between the country regions, has led to distinct virus spread dynamics. Here, we investigate the spatiotemporal dispersion of SARS-CoV-2 lineages in the Pernambuco state (Northeast Brazil) throughout the distinct epidemiological scenarios that unfolded in the first 2 years of the pandemic. We generated a total of 1,389 new SARS-CoV-2 genomes from June 2020 to August 2021. This sampling captured the arrival, communitary transmission, and the circulation of the B1.1, B.1.1.28, and B.1.1.33 lineages; the emergence of the former variant of interest P.2; and the emergence and fast replacement of all previous variants by the more transmissible variant of concern P.1 (Gamma). Based on the incidence and lineage spread pattern, we observed an East-to-West to inner state pattern of transmission, which is in agreement with the transmission of more populous metropolitan areas to medium- and small-size country-side cities in the state. Such transmission patterns may be partially explained by the main routes of traffic across municipalities in the state. Our results highlight that the fine-grained intrastate analysis of lineages and incidence spread can provide actionable insights for planning future nonpharmacological intervention for air-borne transmissible human pathogens.IMPORTANCEDuring the COVID-19 pandemic, Brazil was one of the most affected countries, mainly due its continental-size, socioeconomic differences among regions, and heterogeneous implementation of intervention methods. In order to investigate SARS-CoV-2 dynamics in the state of Pernambuco, we conducted a spatiotemporal dispersion study, covering the period from June 2020 to August 2021, to comprehend the dynamics of viral transmission during the first 2 years of the pandemic. Throughout this study, we were able to track three significant epidemiological waves of transmission caused by B1.1, B.1.1.28, B.1.1.33, P.2, and P.1 lineages. These analyses provided valuable insights into the evolution of the epidemiological landscape, contributing to a deeper understanding of the dynamics of virus transmission during the early years of the pandemic in the state of Pernambuco.


Subject(s)
COVID-19 , SARS-CoV-2 , COVID-19/transmission , COVID-19/epidemiology , COVID-19/virology , Humans , Brazil/epidemiology , SARS-CoV-2/genetics , SARS-CoV-2/classification , Spatio-Temporal Analysis , Genome, Viral , Phylogeny , Pandemics
10.
Trop Med Infect Dis ; 9(3)2024 Feb 24.
Article in English | MEDLINE | ID: mdl-38535877

ABSTRACT

Aedes aegypti and Culex quinquefasciatus mosquitoes are vectors of different arboviruses that cause a large burden of disease in humans worldwide. A key step towards reducing the impact of arboviruses on humans can be achieved through integrated mosquito surveillance and control approaches. We carried out an integrated approach of mosquito surveillance and control actions to reduce populations of these insects along with a viral surveillance in a neighborhood of Recife (Northeastern Brazil) with high mosquito densities and arbovirus transmission. The actions were carried out in 40 different houses in the Nova Descoberta neighborhood. The area was divided into two groups, the control group using tools to monitor the mosquito density (1 OVT; 1 Double BR-ovt; monthly capture of alates) and the experimental group with control actions using surveillance tools in an intensified way (2 OVTs; 2 Double BR-ovts; fortnightly capture of alates; toxic baits). We evaluated the study's impact on the mosquito density via the Egg Density (ED) and Adult Density (AD) over a period of 12 cycles of 28 days each. The collected adult mosquitoes were processed via RT-qPCR for DENV, CHIKV and ZIKV and, subsequently, the Minimum Infection Rate (MIR) was calculated. After 12 cycles, we observed a 91% and 99% reduction in Aedes ED and AD in the monitored properties, as well as a 76% reduction in the AD of Cx. quinquefasciatus in the same properties. Moreover, we detected circulating arboviruses (DENV and ZIKV) in 19.52% of captured adult mosquitoes. We show that enhancing entomological surveillance tools can aid in the early detection of possible risk areas based on vector mosquito population numbers. Additionally, the detection of important arboviruses such as ZIKV and DENV raises awareness and allows for a better selection of risk areas and silent virus spread. It offers supplementary information for guiding emergency mosquito control measures in the target area. The goal is to minimize human-vector interactions and, subsequently, to lower the risk of transmitting circulating arboviruses.

11.
Genet Mol Biol ; 46(3 Suppl 1): e20230138, 2024.
Article in English | MEDLINE | ID: mdl-38373163

ABSTRACT

Transposable elements are repetitive and mobile DNA segments that can be found in virtually all organisms investigated to date. Their complex structure and variable nature are particularly challenging from the genomic annotation point of view. Many softwares have been developed to automate and facilitate TEs annotation at the genomic level, but they are highly heterogeneous regarding documentation, usability and methods. In this review, we revisited the existing software for TE genomic annotation, concentrating on the most often used ones, the methodologies they apply, and usability. Building on the state of the art of TE annotation software we propose best practices and highlight the strengths and weaknesses from the available solutions.

12.
Emerg Infect Dis ; 30(2): 310-320, 2024 Feb.
Article in English | MEDLINE | ID: mdl-38270216

ABSTRACT

We generated 238 Zika virus (ZIKV) genomes from 135 persons in Brazil who had samples collected over 1 year to evaluate virus persistence. Phylogenetic inference clustered the genomes together with previously reported ZIKV strains from northern Brazil, showing that ZIKV has been remained relatively stable over time. Temporal phylogenetic analysis revealed limited within-host diversity among most ZIKV-persistent infected associated samples. However, we detected unusual virus temporal diversity from >5 persons, uncovering the existence of divergent genomes within the same patient. All those patients showed an increase in neutralizing antibody levels, followed by a decline at the convalescent phase of ZIKV infection. Of interest, in 3 of those patients, titers of neutralizing antibodies increased again after 6 months of ZIKV infection, concomitantly with real-time reverse transcription PCR re-positivity, supporting ZIKV reinfection events. Altogether, our findings provide evidence for the existence of ZIKV reinfection events.


Subject(s)
Zika Virus Infection , Zika Virus , Humans , Zika Virus/genetics , Zika Virus Infection/epidemiology , Antibody Formation , Brazil/epidemiology , Phylogeny , Reinfection , Antibodies, Neutralizing
13.
Gene ; 888: 147785, 2023 Dec 20.
Article in English | MEDLINE | ID: mdl-37689222

ABSTRACT

Transposable elements (TEs) are mobile repetitive DNA sequences that can transfer horizontally between species. Due to their mutagenic characteristics, TEs are associated with different evolutionary events, including chromosomal rearrangements that are abundant in the beetle Euchroma gigantea. In order to understand more in depth the impact of TEs on the genomic evolution of E. gigantea, we characterized the E. gigantea mobilome and evaluated the horizontal transfer of Tc1-Mariner elements. Genomic sequencing data was generated on the Illumina Hiseq plataform, from a specimen (Northeast lineage) collected in Recife, Pernambuco - Brazil. The TEs were characterized by two independent approaches based on the clustering and assembly of highly repetitive sequences, the RepeatExplorer and dnaPipeTE. The sequences obtained were further characterized using ORFfinder and CD-Search, to obtain the TEs' potential coding proteins and verify the presence and integrity of known TE domains. Evidence for horizontal transfer was evaluated by nucleotide and protein genetic distance between TEs from E. gigantea and other species and phylogenetic incongruences detected between TEs and hosts phylogenetic trees. The mobilome of E. gigantea represents about 21 to 26% of its genome. This mobilome is composed of TEs from 31 superfamilies, belonging to different classes and most known orders of TEs. Several types of TEs with intact domains were observed with emphasis on Tc1-Mariner suggesting the presence of potentially autonomous elements. This superfamily also stands out for having the greatest abundance and diversity, with TEs being classified into four families. When compared to TEs deposited in databases, Mariner TEs stood out as having the highest nucleotide identity (above 90%) with TEs from phylogenetically distant species, such as ants and bees. Altogether these results suggest that E. gigantea Mariner TEs underwent multiple horizontal transfer events to other insect species.


Subject(s)
Coleoptera , Humans , Animals , Coleoptera/genetics , Phylogeny , DNA Transposable Elements/genetics , Genomics , Nucleotides , Evolution, Molecular
14.
Mem Inst Oswaldo Cruz ; 118: e220259, 2023.
Article in English | MEDLINE | ID: mdl-37531506

ABSTRACT

BACKGROUND: Chikungunya is a mosquito-borne virus that has been causing large outbreaks in the Americas since 2014. In Brazil, Asian-Caribbean (AC) and East-Central-South-African (ECSA) genotypes have been detected and lead to large outbreaks in several Brazilian states. In Rio Grande do Sul (RS), the southernmost state of Brazil, the first cases were reported in 2016. OBJECTIVES AND METHODS: We employed genome sequencing and epidemiological investigation to characterise the Chikungunya fever (CHIKF) burden in RS between 2017-2021. FINDINGS: We detected an increasing CHIKF burden linked to travel associated introductions and communitary transmission of distinct lineages of the ECSA genotype during this period. MAIN CONCLUSIONS: Until 2020, CHIKV introductions were most travel associated and transmission was limited. Then, in 2021, the largest outbreak occurred in the state associated with the introduction of a new ECSA lineage. CHIKV outbreaks are likely to occur in the near future due to abundant competent vectors and a susceptible population, exposing more than 11 million inhabitants to an increasing infection risk.


Subject(s)
Chikungunya Fever , Chikungunya virus , Animals , Humans , Chikungunya virus/genetics , Brazil/epidemiology , Travel , Phylogeny , Mosquito Vectors , Disease Outbreaks , Genotype
15.
Heliyon ; 9(8): e18994, 2023 Aug.
Article in English | MEDLINE | ID: mdl-37600420

ABSTRACT

In late 2021, a new variant of SARS-CoV-2 called Omicron emerged, replacing Delta worldwide. Although it has been associated with a lower risk of hospitalization and severe forms of COVID-19, there is little evidence of its relationship with specific symptoms and viral load. The aim of this study was to verify the relationship between Delta and Omicron variants of concern, viral load, and the occurrence of symptoms in individuals with COVID-19. Nasopharyngeal swab samples were collected and sequenced from patients with COVID-19 from the Northeast Region of Brazil between August 2021 and March 2022. The results showed a gradual replacement of the Delta variant by the Omicron variant during the study period. A total of 316 samples (157 Delta and 159 Omicron) were included. There was a higher prevalence of symptoms in Delta-infected individuals, such as coryza, olfactory and taste disturbances, headache, and myalgia. There was no association between viral load and the variants analyzed. The results reported here contribute to the understanding of the symptoms associated with the Delta and Omicron variants in individuals affected by COVID-19.

18.
Virol J ; 20(1): 83, 2023 05 02.
Article in English | MEDLINE | ID: mdl-37131237

ABSTRACT

Hepatitis E virus (HEV) circulation in humans and swine has been extensively studied in South America over the last two decades. Nevertheless, only 2.1% of reported HEV strains are available as complete genome sequences. Therefore, many clinical, epidemiological, and evolutionary aspects of circulating HEV in the continent still need to be clarified. Here, we conducted a retrospective evolutionary analysis of one human case and six swine HEV strains previously reported in northeastern, southern, and southeastern Brazil. We obtained two complete and four nearly complete genomic sequences. Evolutionary analysis comparing the whole genomic and capsid gene sequences revealed high genetic variability. This included the circulation of at least one unrecognized unique South American subtype. Our results corroborate that sequencing the whole capsid gene could be used as an alternative for HEV subtype assignment in the absence of complete genomic sequences. Moreover, our results substantiate the evidence for zoonotic transmission by comparing a larger genomic fragment recovered from the sample of the autochthonous human hepatitis E case. Further studies should continuously investigate HEV genetic diversity and zoonotic transmission of HEV in South America.


Subject(s)
Hepatitis E virus , Swine , Humans , Animals , Hepatitis E virus/genetics , Brazil/epidemiology , Retrospective Studies , Sequence Analysis, DNA , Genotype , Phylogeny
19.
Sci Rep ; 13(1): 7306, 2023 05 05.
Article in English | MEDLINE | ID: mdl-37147348

ABSTRACT

The rapid spread of the SARS-CoV-2 Variant of Concern (VOC) Gamma in Amazonas during early 2021 fueled a second large COVID-19 epidemic wave and raised concern about the potential role of reinfections. Very few cases of reinfection associated with the VOC Gamma have been reported to date, and their potential impact on clinical, immunological, and virological parameters remains largely unexplored. Here we describe 25 cases of SARS-CoV-2 reinfection in Brazil. SARS-CoV-2 genomic analysis confirmed that individuals were primo-infected with distinct viral lineages between March and December 2020 (B.1.1, B.1.1.28, B.1.1.33, B.1.195, and P.2) and reinfected with the VOC Gamma between 3 to 12 months after primo-infection. We found a similar mean cycle threshold (Ct) value and limited intra-host viral diversity in both primo-infection and reinfection samples. Sera of 14 patients tested 10-75 days after reinfection displayed detectable neutralizing antibodies (NAb) titers against SARS-CoV-2 variants that circulated before (B.1.*), during (Gamma), and after (Delta and Omicron) the second epidemic wave in Brazil. All individuals had milder or no symptoms after reinfection, and none required hospitalization. These findings demonstrate that individuals reinfected with the VOC Gamma may display relatively high RNA viral loads at the upper respiratory tract after reinfection, thus contributing to onward viral transmissions. Despite this, our study points to a low overall risk of severe Gamma reinfections, supporting that the abrupt increase in hospital admissions and deaths observed in Amazonas and other Brazilian states during the Gamma wave was mostly driven by primary infections. Our findings also indicate that most individuals analyzed developed a high anti-SARS-CoV-2 NAb response after reinfection that may provide some protection against reinfection or disease by different SARS-CoV-2 variants.


Subject(s)
COVID-19 , SARS-CoV-2 , Humans , SARS-CoV-2/genetics , Brazil/epidemiology , COVID-19/epidemiology , Antibody Diversity , Gamma Rays , Reinfection , Patient Acuity
20.
Genet Mol Biol ; 46(1): e20220218, 2023.
Article in English | MEDLINE | ID: mdl-37036390

ABSTRACT

Endogenous virus elements (EVEs) are viral-derived sequences integrated into their host genomes. EVEs of the Jingchuvirales order were detected in a wide range of insect genomes covering several distantly related families. Moreover, Jingchuvirales-derived glycoproteins were recently associated by our group with the origin of a putative new retrovirus based on a glycoprotein captured by a mosquito retrotransposon. But, except for mosquitoes, there is a lack of a more detailed understanding of the endogenization mechanism, timing, and frequency per Jingchuvirales viral lineages. Here we screened Jingchuvirales glycoprotein-derived EVEs (Jg-EVEs) in eukaryotic genomes. We found six distinct endogenization events of Jg-EVEs, that belong to two out of five known Jingchuvirales families (Chuviridae and Natareviridae). For seven arthropod families bearing Jg-EVEs there is no register of bona fide circulating chuvirus infection. Hence, our results show that Jingchuvirales viruses infected or still infect these host families. Although we found abundant evidence of LTR-Gypsy retrotransposons fragments associated with the glycoprotein in Hymenoptera and other insect orders, our results show that the widespread distribution of Jingchuvirales glycoproteins in extant Arhtropods is a result of multiple ancient endogenization events and that these virus fossils are being vertically inherited in Arthropods genomes for millions of years.

SELECTION OF CITATIONS
SEARCH DETAIL
...