Multiple Deeply Divergent Denisovan Ancestries in Papuans.

Genome sequences are known for two archaic hominins-Neanderthals and Denisovans-which interbred with anatomically modern humans as they dispersed out of Africa. We identified high-confidence archaic haplotypes in 161 new genomes spanning 14 island groups in Island Southeast Asia and New Guinea and found large stretches of DNA that are inconsistent with a single introgressing Denisovan origin. Instead, modern Papuans carry hundreds of gene variants from two deeply divergent Denisovan lineages that separated over 350 thousand years ago. Spatial and temporal structure among these lineages suggest that introgression from one of these Denisovan groups predominantly took place east of the Wallace line and continued until near the end of the Pleistocene. A third Denisovan lineage occurs in modern East Asians. This regional mosaic suggests considerable complexity in archaic contact, with modern humans interbreeding with multiple Denisovan groups that were geographically isolated from each other over deep evolutionary time.

Physiological and Genetic Adaptations to Diving in Sea Nomads.

Understanding the physiology and genetics of human hypoxia tolerance has important medical implications, but this phenomenon has thus far only been investigated in high-altitude human populations. Another system, yet to be explored, is humans who engage in breath-hold diving. The indigenous Bajau people ("Sea Nomads") of Southeast Asia live a subsistence lifestyle based on breath-hold diving and are renowned for their extraordinary breath-holding abilities. However, it is unknown whether this has a genetic basis. Using a comparative genomic study, we show that natural selection on genetic variants in the PDE10A gene have increased spleen size in the Bajau, providing them with a larger reservoir of oxygenated red blood cells. We also find evidence of strong selection specific to the Bajau on BDKRB2, a gene affecting the human diving reflex. Thus, the Bajau, and possibly other diving populations, provide a new opportunity to study human adaptation to hypoxia tolerance. VIDEO ABSTRACT.

Climate velocities and species tracking in global mountain regions.

Mountain ranges contain high concentrations of endemic species and are indispensable refugia for lowland species that are facing anthropogenic climate change1,2. Forecasting biodiversity redistribution hinges on assessing whether species can track shifting isotherms as the climate warms3,4. However, a global analysis of the velocities of isotherm shifts along elevation gradients is hindered by the scarcity of weather stations in mountainous regions5. Here we address this issue by mapping the lapse rate of temperature (LRT) across mountain regions globally, both by using satellite data (SLRT) and by using the laws of thermodynamics to account for water vapour6 (that is, the moist adiabatic lapse rate (MALRT)). By dividing the rate of surface warming from 1971 to 2020 by either the SLRT or the MALRT, we provide maps of vertical isotherm shift velocities. We identify 17 mountain regions with exceptionally high vertical isotherm shift velocities (greater than 11.67 m per year for the SLRT; greater than 8.25 m per year for the MALRT), predominantly in dry areas but also in wet regions with shallow lapse rates; for example, northern Sumatra, the Brazilian highlands and southern Africa. By linking these velocities to the velocities of species range shifts, we report instances of close tracking in mountains with lower climate velocities. However, many species lag behind, suggesting that range shift dynamics would persist even if we managed to curb climate-change trajectories. Our findings are key for devising global conservation strategies, particularly in the 17 high-velocity mountain regions that we have identified.

Rainforest transformation reallocates energy from green to brown food webs.

Terrestrial animal biodiversity is increasingly being lost because of land-use change1,2. However, functional and energetic consequences aboveground and belowground and across trophic levels in megadiverse tropical ecosystems remain largely unknown. To fill this gap, we assessed changes in energy fluxes across 'green' aboveground (canopy arthropods and birds) and 'brown' belowground (soil arthropods and earthworms) animal food webs in tropical rainforests and plantations in Sumatra, Indonesia. Our results showed that most of the energy in rainforests is channelled to the belowground animal food web. Oil palm and rubber plantations had similar or, in the case of rubber agroforest, higher total animal energy fluxes compared to rainforest but the key energetic nodes were distinctly different: in rainforest more than 90% of the total animal energy flux was channelled by arthropods in soil and canopy, whereas in plantations more than 50% of the energy was allocated to annelids (earthworms). Land-use change led to a consistent decline in multitrophic energy flux aboveground, whereas belowground food webs responded with reduced energy flux to higher trophic levels, down to -90%, and with shifts from slow (fungal) to fast (bacterial) energy channels and from faeces production towards consumption of soil organic matter. This coincides with previously reported soil carbon stock depletion3. Here we show that well-documented animal biodiversity declines with tropical land-use change4-6 are associated with vast energetic and functional restructuring in food webs across aboveground and belowground ecosystem compartments.

Net greenhouse gas balance of fibre wood plantation on peat in Indonesia.

Tropical peatlands cycle and store large amounts of carbon in their soil and biomass1-5. Climate and land-use change alters greenhouse gas (GHG) fluxes of tropical peatlands, but the magnitude of these changes remains highly uncertain6-19. Here we measure net ecosystem exchanges of carbon dioxide, methane and soil nitrous oxide fluxes between October 2016 and May 2022 from Acacia crassicarpa plantation, degraded forest and intact forest within the same peat landscape, representing land-cover-change trajectories in Sumatra, Indonesia. This allows us to present a full plantation rotation GHG flux balance in a fibre wood plantation on peatland. We find that the Acacia plantation has lower GHG emissions than the degraded site with a similar average groundwater level (GWL), despite more intensive land use. The GHG emissions from the Acacia plantation over a full plantation rotation (35.2 ± 4.7 tCO2-eq ha-1 year-1, average ± standard deviation) were around two times higher than those from the intact forest (20.3 ± 3.7 tCO2-eq ha-1 year-1), but only half of the current Intergovernmental Panel on Climate Change (IPCC) Tier 1 emission factor (EF)20 for this land use. Our results can help to reduce the uncertainty in GHG emissions estimates, provide an estimate of the impact of land-use change on tropical peat and develop science-based peatland management practices as nature-based climate solutions.

Genome of a middle Holocene hunter-gatherer from Wallacea.

Much remains unknown about the population history of early modern humans in southeast Asia, where the archaeological record is sparse and the tropical climate is inimical to the preservation of ancient human DNA1. So far, only two low-coverage pre-Neolithic human genomes have been sequenced from this region. Both are from mainland Hòabìnhian hunter-gatherer sites: Pha Faen in Laos, dated to 7939-7751 calibrated years before present (yr cal BP; present taken as AD 1950), and Gua Cha in Malaysia (4.4-4.2 kyr cal BP)1. Here we report, to our knowledge, the first ancient human genome from Wallacea, the oceanic island zone between the Sunda Shelf (comprising mainland southeast Asia and the continental islands of western Indonesia) and Pleistocene Sahul (Australia-New Guinea). We extracted DNA from the petrous bone of a young female hunter-gatherer buried 7.3-7.2 kyr cal BP at the limestone cave of Leang Panninge2 in South Sulawesi, Indonesia. Genetic analyses show that this pre-Neolithic forager, who is associated with the 'Toalean' technocomplex3,4, shares most genetic drift and morphological similarities with present-day Papuan and Indigenous Australian groups, yet represents a previously unknown divergent human lineage that branched off around the time of the split between these populations approximately 37,000 years ago5. We also describe Denisovan and deep Asian-related ancestries in the Leang Panninge genome, and infer their large-scale displacement from the region today.

Reflecting on Indonesia's young academy movement.

In the past three decades, there has been a rise in young academy movements in the Global North and South. Such movements, in at least Germany and the Netherlands, have been shown to be quite effective in connecting scientific work with society. Likewise, these movements share a common goal of developing interdisciplinary collaboration among young scientists, which contributes to the growth of a nation's-but also global-scientific endeavors. This paper focuses on the young academy movement in the fourth-largest country hosting the biggest Muslim population in the world, which is also the third-most populous democracy: Indonesia. We observe that there has been rising awareness among the young generation of scientists in Indonesia of the need to advocate for the use of sciences in responding to upcoming and current multidimensional crises. Science advocacy can be seen in their peer-based identification of Indonesia's future challenges, encompassing the fundamental areas for scientific inquiry, discovery, and intervention. We focus on the Indonesian Young Academy of Sciences (ALMI) and its network of young scientists. We describe ALMI's science communication practice, specifically SAINS45 and Science for Indonesia's Biodiversity, and how they have been useful for policymakers, media, and school engagements. The article closes with a reflection on future directions for the young academy movement in Indonesia and beyond.

Land in limbo: Nearly one third of Indonesia's cleared old-growth forests left idle.

Indonesia has experienced rapid primary forest loss, second only to Brazil in modern history. We examined the fates of Indonesian deforested areas, immediately after clearing and over time, to quantify deforestation drivers in Indonesia. Using time-series satellite data, we tracked degradation and clearing events in intact and degraded natural forests from 1991 to 2020, as well as land use trajectories after forest loss. While an estimated 7.8 Mha (SE = 0.4) of forest cleared during this period had been planted with oil palms by 2020, another 8.8 Mha (SE = 0.4) remained unused. Of the 28.4 Mha (SE = 0.7) deforested, over half were either initially left idle or experienced crop failure before a land use could be detected, and 44% remained unused for 5 y or more. A majority (54%) of these areas were cleared mechanically (not by escaped fires), and in cases where idle lands were eventually converted to productive uses, oil palm plantations were by far the most common outcome. The apparent deliberate creation of idle deforested land in Indonesia and subsequent conversion of idle areas to oil palm plantations indicates that speculation and land banking for palm oil substantially contribute to forest loss, although failed plantations could also contribute to this dynamic. We also found that in Sumatra, few lowland forests remained, suggesting that a lack of remaining forest appropriate for palm oil production, together with an extensive area of banked deforested land, may partially explain slowing forest loss in Indonesia in recent years.

Last appearance of Homo erectus at Ngandong, Java, 117,000-108,000 years ago.

Homo erectus is the founding early hominin species of Island Southeast Asia, and reached Java (Indonesia) more than 1.5 million years ago1,2. Twelve H. erectus calvaria (skull caps) and two tibiae (lower leg bones) were discovered from a bone bed located about 20 m above the Solo River at Ngandong (Central Java) between 1931 and 19333,4, and are of the youngest, most-advanced form of H. erectus5-8. Despite the importance of the Ngandong fossils, the relationship between the fossils, terrace fill and ages have been heavily debated9-14. Here, to resolve the age of the Ngandong evidence, we use Bayesian modelling of 52 radiometric age estimates to establish-to our knowledge-the first robust chronology at regional, valley and local scales. We used uranium-series dating of speleothems to constrain regional landscape evolution; luminescence, 40argon/39argon (40Ar/39Ar) and uranium-series dating to constrain the sequence of terrace evolution; and applied uranium-series and uranium series-electron-spin resonance (US-ESR) dating to non-human fossils to directly date our re-excavation of Ngandong5,15. We show that at least by 500 thousand years ago (ka) the Solo River was diverted into the Kendeng Hills, and that it formed the Solo terrace sequence between 316 and 31 ka and the Ngandong terrace between about 140 and 92 ka. Non-human fossils recovered during the re-excavation of Ngandong date to between 109 and 106 ka (uranium-series minimum)16 and 134 and 118 ka (US-ESR), with modelled ages of 117 to 108 thousand years (kyr) for the H. erectus bone bed, which accumulated during flood conditions3,17. These results negate the extreme ages that have been proposed for the site and solidify Ngandong as the last known occurrence of this long-lived species.

Coupling of Indo-Pacific climate variability over the last millennium.

The Indian Ocean Dipole (IOD) affects climate and rainfall across the world, and most severely in nations surrounding the Indian Ocean1-4. The frequency and intensity of positive IOD events increased during the twentieth century5 and may continue to intensify in a warming world6. However, confidence in predictions of future IOD change is limited by known biases in IOD models7 and the lack of information on natural IOD variability before anthropogenic climate change. Here we use precisely dated and highly resolved coral records from the eastern equatorial Indian Ocean, where the signature of IOD variability is strong and unambiguous, to produce a semi-continuous reconstruction of IOD variability that covers five centuries of the last millennium. Our reconstruction demonstrates that extreme positive IOD events were rare before 1960. However, the most extreme event on record (1997) is not unprecedented, because at least one event that was approximately 27 to 42 per cent larger occurred naturally during the seventeenth century. We further show that a persistent, tight coupling existed between the variability of the IOD and the El Niño/Southern Oscillation during the last millennium. Indo-Pacific coupling was characterized by weak interannual variability before approximately 1590, which probably altered teleconnection patterns, and by anomalously strong variability during the seventeenth century, which was associated with societal upheaval in tropical Asia. A tendency towards clustering of positive IOD events is evident in our reconstruction, which-together with the identification of extreme IOD variability and persistent tropical Indo-Pacific climate coupling-may have implications for improving seasonal and decadal predictions and managing the climate risks of future IOD variability.

Genetic architecture of gene regulation in Indonesian populations identifies QTLs associated with global and local ancestries.

Lack of diversity in human genomics limits our understanding of the genetic underpinnings of complex traits, hinders precision medicine, and contributes to health disparities. To map genetic effects on gene regulation in the underrepresented Indonesian population, we have integrated genotype, gene expression, and CpG methylation data from 115 participants across three island populations that capture the major sources of genomic diversity in the region. In a comparison with European datasets, we identify eQTLs shared between Indonesia and Europe as well as population-specific eQTLs that exhibit differences in allele frequencies and/or overall expression levels between populations. By combining local ancestry and archaic introgression inference with eQTLs and methylQTLs, we identify regulatory loci driven by modern Papuan ancestry as well as introgressed Denisovan and Neanderthal variation. GWAS colocalization connects QTLs detected here to hematological traits, and further comparison with European datasets reflects the poor overall transferability of GWAS statistics across diverse populations. Our findings illustrate how population-specific genetic architecture, local ancestry, and archaic introgression drive variation in gene regulation across genetically distinct and in admixed populations and highlight the need for performing association studies on non-European populations.

Prevalence of adolescent mental disorders in Kenya, Indonesia, and Viet Nam measured by the National Adolescent Mental Health Surveys (NAMHS): a multi-national cross-sectional study.

BACKGROUND: Mental disorders are the leading global cause of health burden among adolescents. However, prevalence data for mental disorders among adolescents in low-income and middle-income countries are scarce with often limited generalisability. This study aimed to generate nationally representative prevalence estimates for mental disorders in adolescents in Kenya, Indonesia, and Viet Nam. METHODS: As part of the National Adolescent Mental Health Surveys (NAMHS), a multinational cross-sectional study, nationally representative household surveys were conducted in Kenya, Indonesia, and Viet Nam between March and December, 2021. Adolescents aged 10-17 years and their primary caregiver were interviewed from households selected randomly according to sampling frames specifically designed to elicit nationally representative results. Six mental disorders (social phobia, generalised anxiety disorder, major depressive disorder, post-traumatic stress disorder, conduct disorder, and attention-deficit hyperactivity disorder) were assessed with the Diagnostic Interview Schedule for Children, Version 5. Suicidal behaviours and self-harm in the past 12 months were also assessed. Prevalence in the past 12 months and past 4 weeks was calculated for each mental disorder and collectively for any mental disorder (ie, of the six mental disorders assessed). Prevalence of suicidal behaviours (ie, ideation, planning, and attempt) and self-harm in the past 12 months was calculated, along with adjusted odds ratios (aORs) to show the association with prevalence of any mental disorder in the past 12 months. Inverse probability weighting was applied to generate national estimates with corresponding 95% CIs. FINDINGS: Final samples consisted of 5155 households (ie, adolescent and primary caregiver pairs) from Kenya, 5664 households from Indonesia, and 5996 households from Viet Nam. In Kenya, 2416 (46·9%) adolescents were male and 2739 (53·1%) were female; in Indonesia, 2803 (49·5%) adolescents were male and 2861 (50·5%) were female; and in Viet Nam, 3151 (52·5%) were male and 2845 (47·4%) were female. Prevalence of any mental disorder in the past 12 months was 12·1% (95% CI 10·9-13·5) in Kenya, 5·5% (4·3-6·9) in Indonesia, and 3·3% (2·7-4·1) in Viet Nam. Prevalence in the past 4 weeks was 9·4% (8·3-10·6) in Kenya, 4·4% (3·4-5·6) in Indonesia, and 2·7% (2·2-3·3) in Viet Nam. The prevalence of suicidal behaviours in the past 12 months was low in all three countries, with suicide ideation ranging from 1·4% in Indonesia (1·0-2·0) and Viet Nam (1·0-1·9) to 4·6% (3·9-5·3) in Kenya, suicide planning ranging from 0·4% in Indonesia (0·3-0·8) and Viet Nam (0·2-0·6) to 2·4% (1·9-2·9) in Kenya, and suicide attempts ranging from 0·2% in Indonesia (0·1-0·4) and Viet Nam (0·1-0·3) to 1·0% (0·7-1·4) in Kenya. The prevalence of self-harm in the past 12 months was also low in all three countries, ranging from 0·9% (0·6-1·3) in Indonesia to 1·2% (0·9-1·7) in Kenya. However, the prevalence of suicidal behaviours and self-harm in the past 12 months was significantly higher among those with any mental disorder in the past 12 months than those without (eg, aORs for suicidal ideation ranged from 7·1 [3·1-15·9] in Indonesia to 14·7 [7·5-28·6] in Viet Nam). INTERPRETATION: NAMHS provides the first national adolescent mental disorders prevalence estimates for Kenya, Indonesia, and Viet Nam. These data can inform mental health and broader health policies in low-income and middle-income countries. FUNDING: The University of Queensland in America (TUQIA) through support from Pivotal Ventures, a Melinda French Gates company.

Earliest hunting scene in prehistoric art.

Humans seem to have an adaptive predisposition for inventing, telling and consuming stories1. Prehistoric cave art provides the most direct insight that we have into the earliest storytelling2-5, in the form of narrative compositions or 'scenes'2,5 that feature clear figurative depictions of sets of figures in spatial proximity to each other, and from which one can infer actions taking place among the figures5. The Upper Palaeolithic cave art of Europe hosts the oldest previously known images of humans and animals interacting in recognizable scenes2,5, and of therianthropes6,7-abstract beings that combine qualities of both people and animals, and which arguably communicated narrative fiction of some kind (folklore, religious myths, spiritual beliefs and so on). In this record of creative expression (spanning from about 40 thousand years ago (ka) until the beginning of the Holocene epoch at around 10 ka), scenes in cave art are generally rare and chronologically late (dating to about 21-14 ka)7, and clear representations of therianthropes are uncommon6-the oldest such image is a carved figurine from Germany of a human with a feline head (dated to about 40-39 ka)8. Here we describe an elaborate rock art panel from the limestone cave of Leang Bulu' Sipong 4 (Sulawesi, Indonesia) that portrays several figures that appear to represent therianthropes hunting wild pigs and dwarf bovids; this painting has been dated to at least 43.9 ka on the basis of uranium-series analysis of overlying speleothems. This hunting scene is-to our knowledge-currently the oldest pictorial record of storytelling and the earliest figurative artwork in the world.

Boat encounter with the 2019 Java bioluminescent milky sea: Views from on-deck confirm satellite detection.

"Milky seas" are massive swaths of uniformly and steadily glowing ocean seen at night. The phenomenon is thought to be caused by luminous bacteria, but details of milky sea composition, structure, cause, and implications in nature remain largely uncertain. Between late July and early September 2019, specialized low-light satellite sensors detected a possible bioluminescent milky sea south of Java, Indonesia, spanning >100,000 km2. Upon learning of these findings, crew members of the yacht Ganesha reached out to confirm and share details of their personal encounter with this same event. Here, we document Ganesha's experience as recalled by the crew, compare their course to satellite data, and assess their photography of this milky sea.

Carbon emissions reductions from Indonesia's moratorium on forest concessions are cost-effective yet contribute little to Paris pledges.

International initiatives for reducing carbon emissions from deforestation and forest degradation (REDD+) could make critical, cost-effective contributions to tropical countries' nationally determined contributions (NDCs). Norway, a key donor of such initiatives, had a REDD+ partnership with Indonesia, offering results-based payments in exchange for emissions reductions calculated against a historical baseline. Central to this partnership was an area-based moratorium on new oil palm, timber, and logging concessions in primary and peatland forests. We evaluate the effectiveness of the moratorium between 2011 and 2018 by applying a matched triple difference strategy to a unique panel dataset. Treated dryland forest inside moratorium areas retained, at most, an average of 0.65% higher forest cover compared to untreated dryland forest outside the moratorium. By contrast, carbon-rich peatland forest was unaffected by the moratorium. Cumulative avoided dryland deforestation from 2011 until 2018 translates into 67.8 million to 86.9 million tons of emissions reductions, implying an effective carbon price below Norway's US$5 per ton price. Based on Norway's price, our estimated cumulative emissions reductions are equivalent to a payment of US$339 million to US$434.5 million. Annually, our estimates suggest a 3 to 4% contribution to Indonesia's NDC commitment of a 29% emissions reduction by 2030. Despite the Indonesia-Norway partnership ending in 2021, reducing emissions from deforestation remains critical for meeting this commitment. Future area-based REDD+ initiatives could build on the moratorium's outcomes by reforming its incentives and institutional arrangements, particularly in peatland forest areas.

Vaccine Take of RV3-BB Rotavirus Vaccine Observed in Indonesian Infants Regardless of HBGA Status.

BACKGROUND: Histo-blood group antigen (HBGA) status may affect vaccine efficacy due to rotavirus strains binding to HBGAs in a P genotype-dependent manner. This study aimed to determine if HBGA status affected vaccine take of the G3P[6] neonatal vaccine RV3-BB. METHODS: DNA was extracted from stool samples collected in a subset (n = 164) of the RV3-BB phase IIb trial in Indonesian infants. FUT2 and FUT3 genes were amplified and sequenced, with any single-nucleotide polymorphisms analyzed to infer Lewis and secretor status. Measures of positive cumulative vaccine take were defined as serum immune response (immunoglobulin A or serum-neutralizing antibody) and/or stool excretion of RV3-BB virus. Participants were stratified by HBGA status and measures of vaccine take. RESULTS: In 147 of 164 participants, Lewis and secretor phenotype were determined. Positive vaccine take was recorded for 144 (97.9%) of 147 participants with the combined phenotype determined. Cumulative vaccine take was not significantly associated with secretor status (relative risk, 1.00 [95% CI, .94-1.06]; P = .97) or Lewis phenotype (relative risk, 1.03 [95% CI, .94-1.14]; P = .33), nor was a difference observed when analyzed by each component of vaccine take. CONCLUSIONS: The RV3-BB vaccine produced positive cumulative vaccine take, irrespective of HBGA status in Indonesian infants.

Molecular characterization of recombinant LSDV isolates from 2022 outbreak in Indonesia through phylogenetic networks and whole-genome SNP-based analysis.

Lumpy skin disease (LSD) is a transboundary viral disease of cattle and water buffaloes caused by the LSD virus, leading to high morbidity, low mortality, and a significant economic impact. Initially endemic to Africa only, LSD has spread to the Middle East, Europe, and Asia in the past decade. The most effective control strategy for LSD is the vaccination of cattle with live-attenuated LSDV vaccines. Consequently, the emergence of two groups of LSDV strains in Asian countries, one closely related to the ancient Kenyan LSDV isolates and the second made of recombinant viruses with a backbone of Neethling-vaccine and field isolates, emphasized the need for constant molecular surveillance. This current study investigated the first outbreak of LSD in Indonesia in 2022. Molecular characterization of the isolate circulating in the country based on selected LSDV-marker genes: RPO30, GPCR, EEV glycoprotein gene, and B22R, as well as whole genome analysis using several analytical tools, indicated the Indonesia LSDV isolate as a recombinant of LSDV_Neethling_vaccine_LW_1959 and LSDV_NI-2490. The analysis clustered the Indonesia_LSDV with the previously reported LSDV recombinants circulating in East and Southeast Asia, but different from the recombinant viruses in Russia and the field isolates in South-Asian countries. Additionally, this study has demonstrated alternative accurate ways of LSDV whole genome analysis and clustering of isolates, including the recombinants, instead of whole-genome phylogenetic tree analysis. These data will strengthen our understanding of the pathogens' origin, the extent of their spread, and determination of suitable control measures required.

Highly Pathogenic Avian Influenza A(H5N1) Virus Clade 2.3.4.4b in Domestic Ducks, Indonesia, 2022.

Highly pathogenic avian influenza A(H5N1) clade 2.3.4.4b viruses were isolated from domestic ducks in South Kalimantan, Indonesia, during April 2022. The viruses were genetically similar to those detected in East Asia during 2021-2022. Molecular surveillance of wild birds is needed to detect potential pandemic threats from avian influenza virus.

Efficacy of Wolbachia-Infected Mosquito Deployments for the Control of Dengue.

BACKGROUND: Aedes aegypti mosquitoes infected with the wMel strain of Wolbachia pipientis are less susceptible than wild-type A. aegypti to dengue virus infection. METHODS: We conducted a cluster-randomized trial involving releases of wMel-infected A. aegypti mosquitoes for the control of dengue in Yogyakarta, Indonesia. We randomly assigned 12 geographic clusters to receive deployments of wMel-infected A. aegypti (intervention clusters) and 12 clusters to receive no deployments (control clusters). All clusters practiced local mosquito-control measures as usual. A test-negative design was used to assess the efficacy of the intervention. Patients with acute undifferentiated fever who presented to local primary care clinics and were 3 to 45 years of age were recruited. Laboratory testing was used to identify participants who had virologically confirmed dengue (VCD) and those who were test-negative controls. The primary end point was symptomatic VCD of any severity caused by any dengue virus serotype. RESULTS: After successful introgression of wMel into the intervention clusters, 8144 participants were enrolled; 3721 lived in intervention clusters, and 4423 lived in control clusters. In the intention-to-treat analysis, VCD occurred in 67 of 2905 participants (2.3%) in the intervention clusters and in 318 of 3401 (9.4%) in the control clusters (aggregate odds ratio for VCD, 0.23; 95% confidence interval [CI], 0.15 to 0.35; P = 0.004). The protective efficacy of the intervention was 77.1% (95% CI, 65.3 to 84.9) and was similar against the four dengue virus serotypes. The incidence of hospitalization for VCD was lower among participants who lived in intervention clusters (13 of 2905 participants [0.4%]) than among those who lived in control clusters (102 of 3401 [3.0%]) (protective efficacy, 86.2%; 95% CI, 66.2 to 94.3). CONCLUSIONS: Introgression of wMel into A. aegypti populations was effective in reducing the incidence of symptomatic dengue and resulted in fewer hospitalizations for dengue among the participants. (Funded by the Tahija Foundation and others; AWED ClinicalTrials.gov number, NCT03055585; Indonesia Registry number, INA-A7OB6TW.).