Lidar reveals pre-Hispanic low-density urbanism in the Bolivian Amazon.

Archaeological remains of agrarian-based, low-density urbananism1-3 have been reported to exist beneath the tropical forests of Southeast Asia, Sri Lanka and Central America4-6. However, beyond some large interconnected settlements in southern Amazonia7-9, there has been no such evidence for pre-Hispanic Amazonia. Here we present lidar data of sites belonging to the Casarabe culture (around AD 500 to AD 1400)10-13 in the Llanos de Mojos savannah-forest mosaic, southwest Amazonia, revealing the presence of two remarkably large sites (147 ha and 315 ha) in a dense four-tiered settlement system. The Casarabe culture area, as far as known today, spans approximately 4,500 km2, with one of the large settlement sites controlling an area of approximately 500 km2. The civic-ceremonial architecture of these large settlement sites includes stepped platforms, on top of which lie U-shaped structures, rectangular platform mounds and conical pyramids (which are up to 22 m tall). The large settlement sites are surrounded by ranked concentric polygonal banks and represent central nodes that are connected to lower-ranked sites by straight, raised causeways that stretch over several kilometres. Massive water-management infrastructure, composed of canals and reservoirs, complete the settlement system in an anthropogenically modified landscape. Our results indicate that the Casarabe-culture settlement pattern represents a type of tropical low-density urbanism that has not previously been described in Amazonia.

Early Holocene crop cultivation and landscape modification in Amazonia.

The onset of plant cultivation is one of the most important cultural transitions in human history1-4. Southwestern Amazonia has previously been proposed as an early centre of plant domestication, on the basis of molecular markers that show genetic similarities between domesticated plants and wild relatives4-6. However, the nature of the early human occupation of southwestern Amazonia, and the history of plant cultivation in this region, are poorly understood. Here we document the cultivation of squash (Cucurbita sp.) at about 10,250 calibrated years before present (cal. yr BP), manioc (Manihot sp.) at about 10,350 cal. yr BP and maize (Zea mays) at about 6,850 cal. yr BP, in the Llanos de Moxos (Bolivia). We show that, starting at around 10,850 cal. yr BP, inhabitants of this region began to create a landscape that ultimately comprised approximately 4,700 artificial forest islands within a treeless, seasonally flooded savannah. Our results confirm that the Llanos de Moxos is a hotspot for early plant cultivation and demonstrate that-ever since their arrival in Amazonia-humans have markedly altered the landscape, with lasting repercussions for habitat heterogeneity and species conservation.

Diverse mathematical knowledge among indigenous Amazonians.

We investigate number and arithmetic learning among a Bolivian indigenous people, the Tsimane', for whom formal schooling is comparatively recent in history and variable in both extent and consistency. We first present a large-scale meta-analysis on child number development involving over 800 Tsimane' children. The results emphasize the impact of formal schooling: Children are only found to be full counters when they have attended school, suggesting the importance of cultural support for early mathematics. We then test especially remote Tsimane' communities and document the development of specialized arithmetical knowledge in the absence of direct formal education. Specifically, we describe individuals who succeed on arithmetic problems involving the number five-which has a distinct role in the local economy-even though they do not succeed on some lower numbers. Some of these participants can perform multiplication with fives at greater accuracy than addition by one. These results highlight the importance of cultural factors in early mathematics and suggest that psychological theories of number where quantities are derived from lower numbers via repeated addition (e.g., a successor function) are unlikely to explain the diversity of human mathematical ability.

Natural selection of immune and metabolic genes associated with health in two lowland Bolivian populations.

A growing body of work has addressed human adaptations to diverse environments using genomic data, but few studies have connected putatively selected alleles to phenotypes, much less among underrepresented populations such as Amerindians. Studies of natural selection and genotype-phenotype relationships in underrepresented populations hold potential to uncover previously undescribed loci underlying evolutionarily and biomedically relevant traits. Here, we worked with the Tsimane and the Moseten, two Amerindian populations inhabiting the Bolivian lowlands. We focused most intensively on the Tsimane, because long-term anthropological work with this group has shown that they have a high burden of both macro and microparasites, as well as minimal cardiometabolic disease or dementia. We therefore generated genome-wide genotype data for Tsimane individuals to study natural selection, and paired this with blood mRNA-seq as well as cardiometabolic and immune biomarker data generated from a larger sample that included both populations. In the Tsimane, we identified 21 regions that are candidates for selective sweeps, as well as 5 immune traits that show evidence for polygenic selection (e.g., C-reactive protein levels and the response to coronaviruses). Genes overlapping candidate regions were strongly enriched for known involvement in immune-related traits, such as abundance of lymphocytes and eosinophils. Importantly, we were also able to draw on extensive phenotype information for the Tsimane and Moseten and link five regions (containing PSD4, MUC21 and MUC22, TOX2, ANXA6, and ABCA1) with biomarkers of immune and metabolic function. Together, our work highlights the utility of pairing evolutionary analyses with anthropological and biomedical data to gain insight into the genetic basis of health-related traits.

Chapare Hemorrhagic Fever and Virus Detection in Rodents in Bolivia in 2019.

BACKGROUND: In June 2019, the Bolivian Ministry of Health reported a cluster of cases of hemorrhagic fever that started in the municipality of Caranavi and expanded to La Paz. The cause of these cases was unknown. METHODS: We obtained samples for next-generation sequencing and virus isolation. Human and rodent specimens were tested by means of virus-specific real-time quantitative reverse-transcriptase-polymerase-chain-reaction assays, next-generation sequencing, and virus isolation. RESULTS: Nine cases of hemorrhagic fever were identified; four of the patients with this illness died. The etiologic agent was identified as Mammarenavirus Chapare mammarenavirus, or Chapare virus (CHAPV), which causes Chapare hemorrhagic fever (CHHF). Probable nosocomial transmission among health care workers was identified. Some patients with CHHF had neurologic manifestations, and those who survived had a prolonged recovery period. CHAPV RNA was detected in a variety of human body fluids (including blood; urine; nasopharyngeal, oropharyngeal, and bronchoalveolar-lavage fluid; conjunctiva; and semen) and in specimens obtained from captured small-eared pygmy rice rats (Oligoryzomys microtis). In survivors of CHHF, viral RNA was detected up to 170 days after symptom onset; CHAPV was isolated from a semen sample obtained 86 days after symptom onset. CONCLUSIONS: M. Chapare mammarenavirus was identified as the etiologic agent of CHHF. Both spillover from a zoonotic reservoir and possible person-to-person transmission were identified. This virus was detected in a rodent species, O. microtis. (Funded by the Bolivian Ministry of Health and others.).

Very large fluxes of methane measured above Bolivian seasonal wetlands.

Methane (CH4) mole fractions from the large semiseasonal Llanos de Moxos wetlands (∼70,000 km2) in northern Bolivia were measured by aircraft flights and ground sampling during early March 2019 (late wet season). Daily fluxes of CH4 determined from the measurements using box models and inverse modeling were between 168 (± 50) and 456 (± 145) mg CH4⋅m-2⋅d-1 for the areas overflown, very high compared with those of previous Amazon basin studies. If the seasonality of the CH4 emissions is comparable to other parts of the Amazon Basin, the region could contribute as much as 8% of annual Amazonian CH4 emissions.

Pediatric neurosurgical medulloblastoma outcomes in La Paz, Bolivia: How a Lower Middle-Income Country (LMIC) institution in South America compares to the United States.

BACKGROUND: How pediatric medulloblastoma patients fare in Lower Middle-Income Country (LMICs) in South America is not well understood. Correspondingly, the aim of this study was to summarize the pediatric neurosurgical experience of an institution in La Paz, and compare outcomes to that of a generalized High Income Country (HIC) United States (US) experience. METHODS: A retrospective review of all pediatric neurosurgical medulloblastoma patients at the Children's Hospital of La Paz, Bolivia (Hospital del Niño "Dr. Ovidio Aliaga Uria") between 2014 and 2023 was conducted and compared to a generalized US experience abstracted from the US National Cancer Database (NCDB) and National Inpatient Sample (NIS) databases. Categorical, continuous and survival data were statistically summarized and compared. RESULTS: A total of 24 pediatric medulloblastoma patients underwent neurosurgical treatment at the Hospital del Niño. In this La Paz cohort, there were 15 (63%) males and 9 (38%) females, with a mean age of 5.6 years old at diagnosis. The majority of patients underwent subtotal resection (STR, 79%), while the remaining patients underwent biopsy only. Ten (42%) patients expired during their hospitalization, and mean length of stay overall was 39 days. Only 8 (33%) patients received adjuvant treatment after surgery. Median overall survival from diagnosis in the La Paz cohort was 1.9 months. Compared to the US databases, the La Paz cohort experienced significantly more emergency room admissions for surgery, less gross total resection, more STR, more return to operating room for ventriculoperitoneal shunting, more bacteremia, more tracheostomy procedures, more percutaneous gastrostomy placements, longer lengths of stay, less adjuvant chemotherapy, less radiation therapy, shorter follow-up, and ultimately, significantly shorter overall survival (all P < 0.050). CONCLUSIONS: Pediatric neurosurgical medulloblastoma outcomes at the Children's Hospital of La Paz, Bolivia are significantly inferior to that of a generalized US experience. Future research is required to identify institution- and country-specific initiatives to improve discrepancies between institutions in LMICs in South America compared to HICs.

Oxidative potential of particulate matter and its association to respiratory health endpoints in high-altitude cities in Bolivia.

Exposure to particulate matter (PM) pollution is a significant health risk, driving the search for innovative metrics that more accurately reflect the potential harm to human health. Among these, oxidative potential (OP) has emerged as a promising health-based metric, yet its application and relevance across different environments remain to be further explored. This study, set in two high-altitude Bolivian cities, aims to identify the most significant sources of PM-induced oxidation in the lungs and assess the utility of OP in assessing PM health impacts. Utilizing two distinct assays, OPDTT and OPDCFH, we measured the OP of PM samples, while also examining the associations between PM mass, OP, and black carbon (BC) concentrations with hospital visits for acute respiratory infections (ARI) and pneumonia over a range of exposure lags (0-2 weeks) using a Poisson regression model adjusted for meteorological conditions. The analysis also leveraged Positive Matrix Factorization (PMF) to link these health outcomes to specific PM sources, building on a prior source apportionment study utilizing the same dataset. Our findings highlight anthropogenic combustion, particularly from traffic and biomass burning, as the primary contributors to OP in these urban sites. Significant correlations were observed between both OPDTT and PM2.5 concentration exposure and ARI hospital visits, alongside a notable association with pneumonia cases and OPDTT levels. Furthermore, PMF analysis demonstrated a clear link between traffic-related pollution and increased hospital admissions for respiratory issues, affirming the health impact of these sources. These results underscore the potential of OPDTT as a valuable metric for assessing the health risks associated with acute PM exposure, showcasing its broader application in environmental health studies.

Genomic analysis of inbreeding level, kinship and breed relationships in Creole cattle from South America.

The conservation of animal genetic resources refers to measures taken to prevent the loss of genetic diversity in livestock populations, including the protection of breeds from extinction. Creole cattle populations have suffered a drastic reduction in recent decades owing to absorbent crosses or replacement with commercial breeds of European or Indian origin. Genetic characterization can serve as a source of information for conservation strategies to maintain genetic variation. The objective of this work was to evaluate the levels of inbreeding and kinship through the use of genomic information. A total of 903 DNAs from 13 cattle populations from Argentina, Bolivia and Uruguay were genotyped using an SNP panel of 48 K. Also, a dataset of 76 K SNPs from Peruvian Creole was included. Two inbreeding indices (FROH and Fhat2) and kinship relationships were calculated. In addition, effective population size (Ne), linkage disequilibrium, population composition and phylogenetic relationships were estimated. In Creole cattle, FROH ranged from 0.14 to 0.03, and Fhat2 was close to zero. The inferred Ne trends exhibited a decline toward the present for all populations, whereas Creole cattle presented a lower magnitude of Ne than foreign breeds. Cluster analysis clearly differentiated the taurine and Zebu components (K2) and showed that Bolivian Creole cattle presented Zebu gene introgression. Despite the population reduction, Creole populations did not present extreme values of consanguinity and kinship and maintain high levels of genetic diversity. The information obtained in this work may be useful for planning conservation programmes for these valuable local animal genetic resources.

Cultural adaptation to Bolivian Quechua and psychometric analysis of the Patient Health Questionnaire PHQ-9.

OBJECTIVE: Cultural adaptation of the Patient Health Questionnaire-PHQ-9 to Bolivian Quechua and analysis of the internal structure validity, reliability, and measurement invariance by sociodemographic variables. METHODS: The PHQ-9 was translated and back-translated (English-Quechua-English) to optimise translation. For the cultural adaptation, experts, and people from the target population (e.g., in focus groups) verified the suitability of the translated PHQ-9. For the psychometric analysis, we performed a Confirmatory Factor Analysis (CFA) to evaluate internal validity, calculated α and ω indices to assess reliability, and performed a Multiple Indicator, Multiple Cause (MIMIC) model for evaluating measurement invariance by sex, age, marital status, educational level and residence. We used standard goodness-of-fit indices to interpret both CFA results. RESULTS: The experts and focus groups improved the translated PHQ-9, making it clear and culturally equivalent. For the psychometric analysis, we included data from 397 participants, from which 73.3% were female, 33.0% were 18-30 years old, 56.7% reported primary school studies, 63.2% were single, and 62.0% resided in urban areas. In the CFA, the single-factor model showed adequate fit (Comparative Fit Index = 0.983; Tucker-Lewis Index = 0.977; Standardized Root Mean Squared Residual = 0.046; Root Mean Squared Error of Approximation = 0.069), while the reliability was optimal (α = 0.869-0.877; ω = 0.874-0.885). The invariance was confirmed across all sociodemographic variables (Change in Comparative Fit Index (delta) or Root Mean Square Error of Approximation (delta) < 0.01). CONCLUSIONS: The PHQ-9 adapted to Bolivian Quechua offers a valid, reliable and invariant unidimensional measurement across groups by sex, age, marital status, educational level and residence.

Global and country-level estimates of human population at high altitude.

Estimates of the global population of humans living at high altitude vary widely, and such data at the country level are unavailable. Herein, we use a geographic information system (GIS)-based approach to quantify human population at 500-m elevation intervals for each country. Based on georeferenced data for population (LandScan Global 2019) and elevation (Global Multiresolution Terrain Elevation Data), 500.3 million humans live at ≥1,500 m, 81.6 million at ≥2,500 m, and 14.4 million at ≥3,500 m. Ethiopia has the largest absolute population at ≥1,500 m and ≥2,500 m, while China has the greatest at ≥3,500 m. Lesotho has the greatest percentage of its population above 1,500 m, while Bolivia has the greatest at ≥2,500 m and ≥3,500 m. High altitude presents a myriad of environmental stresses that provoke physiological responses and adaptation, and consequently impact disease prevalence and severity. While the majority of high-altitude physiology research is based upon lowlanders from western, educated, industrialized, rich, and democratic countries ascending to high altitude, the global population distribution of high-altitude residents encourages an increased emphasis on understanding high-altitude physiology, adaptation, epidemiology, and public health in the ∼500 million permanent high-altitude residents.

Quinoa, potatoes, and llamas fueled emergent social complexity in the Lake Titicaca Basin of the Andes.

The Lake Titicaca basin was one of the major centers for cultural development in the ancient world. This lacustrine environment is unique in the high, dry Andean altiplano, and its aquatic and terrestrial resources are thought to have contributed to the florescence of complex societies in this region. Nevertheless, it remains unclear to what extent local aquatic resources, particularly fish, and the introduced crop, maize, which can be grown in regions along the lakeshores, contributed to facilitating sustained food production and population growth, which underpinned increasing social political complexity starting in the Formative Period (1400 BCE to 500 CE) and culminating with the Tiwanaku state (500 to 1100 CE). Here, we present direct dietary evidence from stable isotope analysis of human skeletal remains spanning over two millennia, together with faunal and floral reference materials, to reconstruct foodways and ecological interactions in southern Lake Titicaca over time. Bulk stable isotope analysis, coupled with compound-specific amino acid stable isotope analysis, allows better discrimination between resources consumed across aquatic and terrestrial environments. Together, this evidence demonstrates that human diets predominantly relied on C3 plants, particularly quinoa and tubers, along with terrestrial animals, notably domestic camelids. Surprisingly, fish were not a significant source of animal protein, but a slight increase in C4 plant consumption verifies the increasing importance of maize in the Middle Horizon. These results underscore the primary role of local terrestrial food resources in securing a nutritious diet that allowed for sustained population growth, even in the face of documented climate and political change across these periods.

Divergent Genotype of Hepatitis A Virus in Alpacas, Bolivia, 2019.

Hepatitis A virus (HAV) is a common human pathogen found exclusively in primates. In a molecular and serologic study of 64 alpacas in Bolivia, we detected RNA of distinct HAV in ≈9% of animals and HAV antibodies in ≈64%. Complete-genome analysis suggests a long association of HAV with alpacas.

Rio Negro Virus Infection, Bolivia, 2021.

In May 2021, an agricultural worker originally from Trementinal, Argentina, sought treatment for febrile illness in Tarija, Bolivia, where he resided at the time of illness onset. The patient tested negative for hantavirus RNA, but next-generation sequencing of a serum sample yielded a complete genome for Rio Negro virus.

Genetic portrait of the Amazonian communities of Peru and Bolivia: The legacy of the Takanan-speaking people.

During the colonial period in South America, many autochthonous populations were affected by relocation by European missionary reductions and other factors that impacted and reconfigured their genetic makeup. Presently, the descendants of some "reduced" and other isolated groups are distributed in the Amazonian areas of Peru, Bolivia, and Brazil, and among them, speakers of Takanan and Panoan languages. Based on linguistics, these peoples should be closely related, but so far no DNA comparison studies have been conducted to corroborate a genetic relationship. To clarify these questions, we used a set of 15 short tandem repeats of the non-recombining part of the Y-chromosome (Y-STRs) and mitochondrial DNA (mtDNA) control region sequence data. Paternal line comparisons showed the Takanan-speaking peoples from Peru and Bolivia descended from recent common ancestors; one group was related to Arawakan, Jivaroan, and Cocama and the other to Panoan speakers, consistent with linguistics. Also, a genetic affinity for maternal lines was observed between some Takanan speakers and individuals who spoke different Amazonian languages. Our results supported a shared ancestry of Takanan, Panoan, Cocama, and Jivaroan-speaking communities who appeared to be related to each other and came likely from an early Arawak expansion in the western Amazonia of South America.

Concepts Are Restructured During Language Contact: The Birth of Blue and Other Color Concepts in Tsimane'-Spanish Bilinguals.

Words and the concepts they represent vary across languages. Here we ask if mother-tongue concepts are altered by learning a second language. What happens when speakers of Tsimane', a language with few consensus color terms, learn Bolivian Spanish, a language with more terms? Three possibilities arise: Concepts in Tsimane' may remain unaffected, or they may be remapped, either by Tsimane' terms taking on new meanings or by borrowing Bolivian-Spanish terms. We found that adult bilingual speakers (n = 30) remapped Tsimane' concepts without importing Bolivian-Spanish terms into Tsimane'. All Tsimane' terms become more precise; for example, concepts of monolingual shandyes and yushñus (~green or blue, used synonymously by Tsimane' monolinguals; n = 71) come to reflect the Bolivian-Spanish distinction of verde (~green) and azul (~blue). These results show that learning a second language can change the concepts in the first language.

Spiny but photogenic: Amateur sightings complement herbarium specimens to reveal the bioregions of cacti.

PREMISE: Cacti are characteristic elements of the Neotropical flora and of major interest for biogeographic, evolutionary, and ecological studies. We tested global biogeographic boundaries for Neotropical Cactaceae using specimen-based occurrences, coupled with data from visual observations, as a means to tackle the known collection biases in the family. METHODS: Species richness and record density were assessed for preserved specimens and human observations, and a bioregional scheme tailored to Cactaceae was produced using the interactive web application Infomap Bioregions, based on data from 261,272 point records cleaned through automated and manual steps. RESULTS: We found that areas in Mexico and southwestern USA, in eastern Brazil, and along the Andean region have the greatest density of records and the highest species richness. Human observations complement information from preserved specimens substantially, especially along the Andes. We propose 24 cactus bioregions, among which the most species-rich are northern Mexico/southwestern USA, central Mexico, southern central Mexico, Central America, Mexican Pacific coast, central and southern Andes, northwestern Mexico/extreme southwestern USA, southwestern Bolivia, northeastern Brazil, and Mexico/Baja California. CONCLUSIONS: The bioregionalization proposed shows biogeographic boundaries specific to cacti and can thereby aid further evolutionary, biogeographic, and ecological studies by providing a validated framework for further analyses. This classification builds upon, and is distinctive from, other expert-derived regionalization schemes for other taxa. Our results showcase how observation data, including citizen-science records, can complement traditional specimen-based data for biogeographic research, particularly for taxa with specific specimen collection and preservation challenges and those that are threatened or internationally protected.

When the flood passes, does health return? A short panel examining water and food insecurity, nutrition, and disease after an extreme flood in lowland Bolivia.

OBJECTIVES: Flooding is the most frequent extreme-weather disaster and disproportionately burdens marginalized populations. This article examines how food and water insecurity, blood pressure (BP), nutritional status, and diarrheal and respiratory illnesses changed during the 2 months following a historic flood in lowland Bolivia. METHODS: Drawing on longitudinal data from Tsimane' forager-horticulturalist (n = 118 household heads; n = 129 children) directly after a historic 2014 flood and ~2 months later, we use fixed effects linear regression and random effects logistic regression models to test changes in the markers of well-being and health over the recovery process. RESULTS: Results demonstrated that water insecurity scores decreased significantly 2 month's postflood, while food insecurity scores remained high. Adults' systolic and diastolic BP significantly declined 2 months after the flood's conclusion. Adults experienced losses in measures of adiposity (BMI, sum of four skinfolds, waist circumference). Children gained weight and BMI-for-age Z-scores indicating buffering of children by adults from food stress that mainly occurred in the community closer to the main market town with greater access to food aid. Odds of diarrhea showed a nonsignificant decline, while cough increased significantly for both children and adults 2 months postflood. CONCLUSIONS: Water insecurity and BP improved during the recovery process, while high levels of food insecurity persisted, and nutritional stress and respiratory illness worsened. Not all indicators of well-being and health recover at the same rate after historic flooding events. Planning for multiphase recovery is critical to improve health of marginalized populations after flooding.

Community and provider perceptions and experiences of cervical cancer screening in Rural Bolivia: a qualitative study.

BACKGROUND: Despite efforts to increase cervical cancer screening access in rural Bolivia, uptake remains low. Bolivia has one of the highest cervical cancer mortality rates in the Americas. As it redoubles efforts to deliver Universal Health Care, the Bolivian government needs information on the factors constraining cervical cancer screening access and utilization, especially in rural areas. METHODS: Our qualitative study explored cervical cancer screening barriers and described community and provider perceptions and experiences of care. Bolivian and US researchers analyzed data collected from eight focus groups with male and female community members (n = 80) and interviews with healthcare providers (n = 6) in four purposively selected rural communities in Hernando Siles, Bolivia. Deductive and inductive codes were used to thematically analyze data using MaxQDA software. RESULTS: Four themes emerged from the data: lack of knowledge/misconceptions, health system inadequacy, lack of confidence in providers, and opportunities for improvement. Both men and women displayed misconceptions about the causes of cervical cancer, its consequences, the recommended screening frequency, and the means of accessing care. Providers noted community members' lack of knowledge and low risk-perception as utilization barriers but also highlighted poor health service quality and inconsistent health education as factors. Poor healthcare quality was a significant barrier; this included poor patient-provider communication, lack of transportation to screening facilities, and severe delays in receiving test results. Providers also noted problems with provider training and physical space for screening. Community members reported low confidence in nurses to perform screening, preferring doctors and specialists. They also expressed discomfort in having male healthcare providers conduct screening. Suggestions for improvements included more intensive cervical cancer outreach to rural areas and having specialists train lower-level providers to perform screening. CONCLUSIONS: Our findings suggest that poor healthcare quality has affected screening uptake in addition to physical barriers to care. They indicate a need for initiatives to reduce reporting time for Pap test results, the incorporation of community-based HPV self-sampling into screening protocols, and the implementation of programs to improve community confidence in providers' ability to perform screening.

An epidemiological evaluation of COVID-19 in La paz, Bolivia.

INTRODUCTION: The Plurinational State of Bolivia (Bolivia) has experienced four major waves of coronavirus disease 2019 (COVID-19) so far. Although the ministry of health has been tracking morbidity and mortality through each wave, epidemiology of COVID-19 in Bolivia is not well defined, despite a need for more accurate measurement of the number of cases and deaths to allow for forecasting of the pandemic. This study examined prevalence of COVID-19 at community level, determinants of its occurrence and vaccine effectiveness. METHODS: We conducted a cross-sectional study in La Paz city on 2,775 individuals between March 2020 and February 2022. A structured questionnaire was used to collect data on COVID-19 morbidity, mortality and vaccination status. RESULTS: Of the 2,775 participants, 1,586 (57.1%) were infected with COVID-19, and 187 (6.7%) were suspected cases. The mortality rate was 2.9%. Sinopharm, Johnson & Johnson, Gamaleya, Pfizer-BioNtech, Moderna and AstraZeneka vaccines are in use, and all vaccines have demonstrated effectiveness in reducing the risk of onset. Risk for mortality was significantly lower in the vaccinated group with an odds ratio of 0.037 (95％ confidential interval: 0.01-0.10, p-value: <0.001). CONCLUSIONS: Actual prevalence of COVID-19 in La Paz (the prevalence rate: 63.8%, including suspected case) was higher than that reported by the Ministry of Health and Sports in Bolivia (7.5%). In addition, vaccination has contributed significantly to the control of the COVID-19 epidemic in Bolivia. We believe that our report will be useful for COVID-19 prevention strategies in Bolivia for the future.