Genome-based local dynamics of canine rabies virus epidemiology, transmission, and evolution in Davao City, Philippines, 2018-2019.

Rabies is a fatal zoonotic and neglected tropical disease caused by the rabies virus (RABV) and is associated with neuronal dysfunction and death, with dogs as the predominant carrier. The Philippines plans to eradicate rabies by 2022, but this is challenged with sub-optimal coverage of vaccination programs coupled with sustained transmission chains, making it unable to eradicate the disease. We investigated the dynamics of canine rabies in the highly urbanized Davao City of the Philippines and its neighboring localities by assessing genetic relationships, transmission patterns, selection pressure, and recombination events using the whole genome sequence of 49 RABV cases from June 2018 to May 2019, majority of which (46%) were from the district of Talomo, Davao City. Although phylogeographic clustering was observed, local variants also exhibited genetic sub-lineages. Phylogenetic and spatial transmission analysis provided evidence for intra- and inter-city transmission predominantly through the Talomo district of Davao City. Around 84% of the cases were owned dogs, but the genetic similiarity of RABVs from stray and owned dogs further alluded to the role of the former as transmission vectors. The high rate of improper vaccination among the affected dogs (80%) was also a likely contributor to transmission. The RABV population under Investigation is generally under strong purifying selection with no evidence of vaccine evasion due to the genetic homogeneity of viruses from vaccinated and improperly vaccinated dogs. However, some homologous recombination (HR) events were identified along the G and L genes, also predominantly associated with viruses from Talomo. The complementary findings on epidemiology, transmission, and recombination for Talomo suggest that high incidence areas can be seeds for virus dispersal and evolution. We recommend further Investigations on the possibility of HR in future large-scale genome studies. Finally, districts associated with these phenomena can be targeted for evidence-based local strategies that can help break RABV transmission chains and prevent emergence of novel strains in Davao City.

Filipino DNA variation at 12 X-chromosome short tandem repeat markers.

Demands for solving complex kinship scenarios where only distant relatives are available for testing have risen in the past years. In these instances, other genetic markers such as X-chromosome short tandem repeat (X-STR) markers are employed to supplement autosomal and Y-chromosomal STR DNA typing. However, prior to use, the degree of STR polymorphism in the population requires evaluation through generation of an allele or haplotype frequency population database. This population database is also used for statistical evaluation of DNA typing results. Here, we report X-STR data from 143 unrelated Filipino male individuals who were genotyped via conventional polymerase chain reaction-capillary electrophoresis (PCR-CE) using the 12 X-STR loci included in the Investigator® Argus X-12 kit (Qiagen) and via massively parallel sequencing (MPS) of seven X-STR loci included in the ForenSeq™ DNA Signature Prep kit of the MiSeq® FGx™ Forensic Genomics System (Illumina). Allele calls between PCR-CE and MPS systems were consistent (100% concordance) across seven overlapping X-STRs. Allele and haplotype frequencies and other parameters of forensic interest were calculated based on length (PCR-CE, 12 X-STRs) and sequence (MPS, seven X-STRs) variations observed in the population. Results of our study indicate that the 12 X-STRs in the PCR-CE system are highly informative for the Filipino population. MPS of seven X-STR loci identified 73 X-STR alleles compared with 55 X-STR alleles that were identified solely by length via PCR-CE. Of the 73 sequence-based alleles observed, six alleles have not been reported in the literature. The population data presented here may serve as a reference Philippine frequency database of X-STRs for forensic casework applications.

Early Austronesians: into and out of Taiwan.

A Taiwan origin for the expansion of the Austronesian languages and their speakers is well supported by linguistic and archaeological evidence. However, human genetic evidence is more controversial. Until now, there had been no ancient skeletal evidence of a potential Austronesian-speaking ancestor prior to the Taiwan Neolithic ~6,000 years ago, and genetic studies have largely ignored the role of genetic diversity within Taiwan as well as the origins of Formosans. We address these issues via analysis of a complete mitochondrial DNA genome sequence of an ~8,000-year-old skeleton from Liang Island (located between China and Taiwan) and 550 mtDNA genome sequences from 8 aboriginal (highland) Formosan and 4 other Taiwanese groups. We show that the Liangdao Man mtDNA sequence is closest to Formosans, provides a link to southern China, and has the most ancestral haplogroup E sequence found among extant Austronesian speakers. Bayesian phylogenetic analysis allows us to reconstruct a history of early Austronesians arriving in Taiwan in the north ~6,000 years ago, spreading rapidly to the south, and leaving Taiwan ~4,000 years ago to spread throughout Island Southeast Asia, Madagascar, and Oceania.

Complete mtDNA genomes of Filipino ethnolinguistic groups: a melting pot of recent and ancient lineages in the Asia-Pacific region.

The Philippines is a strategic point in the Asia-Pacific region for the study of human diversity, history and origins, as it is a cross-road for human migrations and consequently exhibits enormous ethnolinguistic diversity. Following on a previous in-depth study of Y-chromosome variation, here we provide new insights into the maternal genetic history of Filipino ethnolinguistic groups by surveying complete mitochondrial DNA (mtDNA) genomes from a total of 14 groups (11 groups in this study and 3 groups previously published) including previously published mtDNA hypervariable segment (HVS) data from Filipino regional center groups. Comparison of HVS data indicate genetic differences between ethnolinguistic and regional center groups. The complete mtDNA genomes of 14 ethnolinguistic groups reveal genetic aspects consistent with the Y-chromosome, namely: diversity and heterogeneity of groups, no support for a simple dichotomy between Negrito and non-Negrito groups, and different genetic affinities with Asia-Pacific groups that are both ancient and recent. Although some mtDNA haplogroups can be associated with the Austronesian expansion, there are others that associate with South Asia, Near Oceania and Australia that are consistent with a southern migration route for ethnolinguistic group ancestors into the Asia-Pacific, with a timeline that overlaps with the initial colonization of the Asia-Pacific region, the initial colonization of the Philippines and a possible separate post-colonization migration into the Philippine archipelago.

Genome-wide data substantiate Holocene gene flow from India to Australia.

The Australian continent holds some of the earliest archaeological evidence for the expansion of modern humans out of Africa, with initial occupation at least 40,000 y ago. It is commonly assumed that Australia remained largely isolated following initial colonization, but the genetic history of Australians has not been explored in detail to address this issue. Here, we analyze large-scale genotyping data from aboriginal Australians, New Guineans, island Southeast Asians and Indians. We find an ancient association between Australia, New Guinea, and the Mamanwa (a Negrito group from the Philippines), with divergence times for these groups estimated at 36,000 y ago, and supporting the view that these populations represent the descendants of an early "southern route" migration out of Africa, whereas other populations in the region arrived later by a separate dispersal. We also detect a signal indicative of substantial gene flow between the Indian populations and Australia well before European contact, contrary to the prevailing view that there was no contact between Australia and the rest of the world. We estimate this gene flow to have occurred during the Holocene, 4,230 y ago. This is also approximately when changes in tool technology, food processing, and the dingo appear in the Australian archaeological record, suggesting that these may be related to the migration from India.

Bridging near and remote Oceania: mtDNA and NRY variation in the Solomon Islands.

Although genetic studies have contributed greatly to our understanding of the colonization of Near and Remote Oceania, important gaps still exist. One such gap is the Solomon Islands, which extend between Bougainville and Vanuatu, thereby bridging Near and Remote Oceania, and include both Austronesian-speaking and Papuan-speaking groups. Here, we describe patterns of mitochondrial DNA (mtDNA) and nonrecombining Y chromosome (NRY) variation in over 700 individuals from 18 populations in the Solomons, including 11 Austronesian-speaking groups, 3 Papuan-speaking groups, and 4 Polynesian Outliers (descended via back migration from Polynesia). We find evidence for ancient (pre-Lapita) colonization of the Solomons in old NRY paragroups as well as from M2-M353, which probably arose in the Solomons ∼9,200 years ago and is the most frequent NRY haplogroup there. There are no consistent genetic differences between Austronesian-speaking and Papuan-speaking groups, suggesting extensive genetic contact between them. Santa Cruz, which is located in Remote Oceania, shows unusually low frequencies of mtDNA and NRY haplogroups of recent Asian ancestry. This is in apparent contradiction with expectations based on archaeological and linguistic evidence for an early (∼3,200 years ago), direct colonization of Santa Cruz by Lapita people from the Bismarck Archipelago, via a migration that "leapfrogged" over the rest of the Solomons. Polynesian Outliers show dramatic island-specific founder events involving various NRY haplogroups. We also find that NRY, but not mtDNA, genetic distance is correlated with the geographic distance between Solomons groups and that historically attested spheres of cultural interaction are associated with the recent genetic structure of Solomons groups, as revealed by mtDNA HV1 sequence and Y-STR haplotype diversity. Our results fill an important lacuna in human genetic studies of Oceania and aid in understanding the colonization and genetic history of this region.

Denisova admixture and the first modern human dispersals into Southeast Asia and Oceania.

It has recently been shown that ancestors of New Guineans and Bougainville Islanders have inherited a proportion of their ancestry from Denisovans, an archaic hominin group from Siberia. However, only a sparse sampling of populations from Southeast Asia and Oceania were analyzed. Here, we quantify Denisova admixture in 33 additional populations from Asia and Oceania. Aboriginal Australians, Near Oceanians, Polynesians, Fijians, east Indonesians, and Mamanwa (a "Negrito" group from the Philippines) have all inherited genetic material from Denisovans, but mainland East Asians, western Indonesians, Jehai (a Negrito group from Malaysia), and Onge (a Negrito group from the Andaman Islands) have not. These results indicate that Denisova gene flow occurred into the common ancestors of New Guineans, Australians, and Mamanwa but not into the ancestors of the Jehai and Onge and suggest that relatives of present-day East Asians were not in Southeast Asia when the Denisova gene flow occurred. Our finding that descendants of the earliest inhabitants of Southeast Asia do not all harbor Denisova admixture is inconsistent with a history in which the Denisova interbreeding occurred in mainland Asia and then spread over Southeast Asia, leading to all its earliest modern human inhabitants. Instead, the data can be most parsimoniously explained if the Denisova gene flow occurred in Southeast Asia itself. Thus, archaic Denisovans must have lived over an extraordinarily broad geographic and ecological range, from Siberia to tropical Asia.

The Y-chromosome landscape of the Philippines: extensive heterogeneity and varying genetic affinities of Negrito and non-Negrito groups.

The Philippines exhibits a rich diversity of people, languages, and culture, including so-called 'Negrito' groups that have for long fascinated anthropologists, yet little is known about their genetic diversity. We report here, a survey of Y-chromosome variation in 390 individuals from 16 Filipino ethnolinguistic groups, including six Negrito groups, from across the archipelago. We find extreme diversity in the Y-chromosome lineages of Filipino groups with heterogeneity seen in both Negrito and non-Negrito groups, which does not support a simple dichotomy of Filipino groups as Negrito vs non-Negrito. Filipino non-recombining region of the human Y chromosome lineages reflect a chronology that extends from after the initial colonization of the Asia-Pacific region, to the time frame of the Austronesian expansion. Filipino groups appear to have diverse genetic affinities with different populations in the Asia-Pacific region. In particular, some Negrito groups are associated with indigenous Australians, with a potential time for the association ranging from the initial colonization of the region to more recent (after colonization) times. Overall, our results indicate extensive heterogeneity contributing to a complex genetic history for Filipino groups, with varying roles for migrations from outside the Philippines, genetic drift, and admixture among neighboring groups.

Y-STR DNA analysis of 154 female child sexual assault cases in the Philippines.

The laboratory evaluated 154 sexual assault cases from four Child Protection Units in the Philippines involving female child victims aged from 2 years to 18 years old. All child victims sought medical attention within 72 h after sexual contact. In 130 cases, the child victim knew the alleged offender and identified them during the interview with the social worker. Penile ejaculation was reported by 68 child victims with varying reports of washing after contact. Overall, 84 child victims admitted having wiped their genitalia prior to the collection of biological samples for DNA testing. Laboratory personnel examined vaginal smears in only 109 cases using a light microscope and reported 23 samples to be positive for sperm cells. Using the PowerPlex® short tandem repeat of the Y chromosome (Y-STR) DNA multiplex system, male DNA was detected in vaginal swab samples from 63 child victims. In 39 cases, positive amplification at 11 Y-STR DNA markers consistent with a single male DNA profile was observed. Twenty-eight of these full single Y-STR DNA profiles were found to be unique when searched in worldwide Y-STR DNA population databases (~40,000 haplotypes), eight haplotypes matching Filipinos and/or Asian haplotypes and one Y-STR DNA profile only matching European, Caucasian, and Latin American haplotypes. Y-STR DNA profiles generated will be compared with reference DNA profiles of alleged offenders once reference samples are submitted to the laboratory.

The human genetic history of East Asia: weaving a complex tapestry.

East Asia encompasses a wide variety of environments, peoples, cultures and languages. Although this review focuses on East Asia, no geographic region can be considered in isolation in terms of human population history, and migrations to and from East Asia have had a major impact. Here, we review the following topics: the initial colonization of East Asia, the direction of migrations between southeast Asia and northern Asia, the genetic relationships of East Asian hunter-gatherers and the genetic impact of various social practices on East Asian populations. By necessity we focus on insights derived from mitochondrial DNA and/or Y-chromosome data; ongoing and future studies of genome-wide SNP or multi-locus re-sequencing data, combined with the use of simulation, model-based methods to infer demographic parameters, will undoubtedly provide additional insights into the population history of East Asia.

Mapping human genetic diversity in Asia.

Asia harbors substantial cultural and linguistic diversity, but the geographic structure of genetic variation across the continent remains enigmatic. Here we report a large-scale survey of autosomal variation from a broad geographic sample of Asian human populations. Our results show that genetic ancestry is strongly correlated with linguistic affiliations as well as geography. Most populations show relatedness within ethnic/linguistic groups, despite prevalent gene flow among populations. More than 90% of East Asian (EA) haplotypes could be found in either Southeast Asian (SEA) or Central-South Asian (CSA) populations and show clinal structure with haplotype diversity decreasing from south to north. Furthermore, 50% of EA haplotypes were found in SEA only and 5% were found in CSA only, indicating that SEA was a major geographic source of EA populations.

Identification of exhumed remains of fire tragedy victims using conventional methods and autosomal/Y-chromosomal short tandem repeat DNA profiling.

In a fire tragedy in Manila in December 1998, one of the worst tragic incidents which resulted in the reported death of 23 children, identity could not be established initially resulting in the burial of still unidentified bodies. Underscoring the importance of identifying each of the human remains, the bodies were exhumed 3 months after the tragedy. We describe here our work, which was the first national case handled by local laboratories wherein conventional and molecular-based techniques were successfully applied in forensic identification. The study reports analysis of DNA obtained from skeletal remains exposed to conditions of burning, burial, and exhumation. DNA typing methods using autosomal and Y-chromosomal short tandem repeat (Y-STR) markers reinforced postmortem examinations using conventional identification techniques. The strategy resulted in the identification of 18 out of the 21 human remains analyzed, overcoming challenges encountered due to the absence of established procedures for the recovery of mass disaster remains. There was incomplete antemortem information to match the postmortem data obtained from the remains of 3 female child victims. Two victims were readily identified due to the availability of antemortem tissues. In the absence of this biologic material, parentage testing was performed using reference blood samples collected from parents and relatives. Data on patrilineal lineage based on common Y-STR haplotypes augmented autosomal DNA typing, particularly in deficiency cases.

Y-STR analysis for detection and objective confirmation of child sexual abuse.

We evaluated 26 child sexual assault cases for the incorporation of Y-STR screening in the routine detection and objective confirmation of sexual contact between the child victim and the perpetrator. Various samples, e.g. vaginal or anal swabs from patients aged 2-17 years old (25 females, 1 male), were collected 6-72 h after the incident. Due to the limited amounts of DNA in these samples, total DNA was extracted using a one-step procedure and screened with autosomal STRs to detect signs of a victim-assailant DNA mixture and with Y-STRs for assailant DNA. Autosomal STRs failed to give signs of victim-assailant DNA mixtures while Y-STRs were detected in 24 of the 26 cases corresponding to a success rate of 92.3%. With the possible presence of both male sperm and/or male epithelial cells in forensic evidence, Y-STR DNA markers were detected regardless of external ejaculation, microscopic detection of sperm and with post-coital intervals of up to 72 h. While only partial profiles were generated owing to low quantities of male DNA present, Y-STR screening results can serve as objective evidence of sexual contact in child sexual abuse cases involving victims who do not have any previous sexual history. This type of evidence can corroborate child victim testimony and spare the child victim from further trauma caused by prolonged forensic investigations and court proceedings. Alternatively, Y-STR screening can provide objective proof of non-involvement of an accused with the victim.

Asian online Y-STR Haplotype Reference Database.

For several years Y-chromosomal microsatellites (short tandem repeats, STRs) have been well established in forensic practice. In this context, the genetic characteristics of the Y chromosome (i.e. its paternal inheritance and lack of recombination) render STRs particularly powerful. However, genetic differences between male populations appear to be larger for Y-STRs than for autosomal STRs, a fact that is most likely due to the higher sensitivity of Y-chromosomal lineages to genetic drift (Forensic Sci Int 118 (2001) 153). The assessment of probabilities for matches between haplotyped male persons or traces/persons requires the typing of a large number of haplotypes in the appropriate reference populations. The haplotype data of a large number of European as well as South and North American populations have been collected and are continuously published online (Y-STR Haplotype Reference Database--YHRD; http://www.ystr.org). The most recent multicentric effort has led to the establishment of an Asian YHRD (http://www.ystr.org/asia) which has been available since January 2002. All databases are maintained and curated at the Institute of Legal Medicine, Humboldt-University, Berlin and will soon be fused to a global repository including populations from all continents.

Identification of two fire victims by comparative nuclear DNA typing of skeletal remains and stored umbilical tissues.

We describe here our collaborative efforts in identifying 2 fatalities of a fire disaster by using a variety of identification techniques. Postmortem findings in both cases were reinforced using Short Tandem Repeat (STR) DNA technology to establish with a high degree of certainty the identities of 2 child victims. STR markers used in the present study include HUMAMEL, HUMCSFIPO, HUMTHO1, HUMvWA, HUMFES/FPS, HUMF13A01, HUMFOLP23, D8S3O6, HUMFGA, and HUMTPOX. Unambiguous identification was made possible through matching DNA profiles generated from skeletal remains with those from umbilical tissues. These tissues were kept by their mothers in accordance with a Philippine tradition and were submitted for DNA analysis. Of the DNA profiles generated from exhumed bone samples of 21 child victims, comparison with the genetic profiles of children A and B obtained from umbilical tissues showed consistent DNA matches with remains 1756 and 1758, respectively.