Extensive genetic admixture between Tai-Kadai-speaking people and their neighbours in the northeastern region of the Yungui Plateau inferred from genome-wide variations.

BACKGROUND: Yungui Plateau in Southwest China is characterized by multi-language and multi-ethnic communities and is one of the regions with the wealthiest ethnolinguistic, cultural and genetic diversity in East Asia. There are numerous Tai-Kadai (TK)-speaking populations, but their detailed evolutionary history and biological adaptations are still unclear. RESULTS: Here, we genotyped genome-wide SNP data of 77 unrelated TK-speaking Zhuang and Dong individuals from the Yungui Plateau and explored their detailed admixture history and adaptive features using clustering patterns, allele frequency differentiation and sharing haplotype patterns. TK-speaking Zhuang and Dong people in Guizhou are closely related to geographically close TK and Hmong-Mien (HM)-speaking populations. Besides, we identified that Guizhou TK-speaking people have a close genetic relationship with Austronesian (AN)-speaking Atayal and Paiwan people, which is supported by the common origin of the ancient Baiyue tribe. We additionally found subtle genetic differences among the newly studied TK people and previously reported Dais via the fine-scale genetic substructure analysis based on the shared haplotype chunks. Finally, we identified specific selection candidate signatures associated with several essential human immune systems and neurological disorders, which could provide evolutionary evidence for the allele frequency distribution pattern of genetic risk loci. CONCLUSIONS: Our comprehensive genetic characterization of TK people suggested the strong genetic affinity within TK groups and extensive gene flow with geographically close HM and Han people. We also provided genetic evidence that supported the common origin hypothesis of TK and AN people. The best-fitted admixture models further suggested that ancestral sources from northern millet farmers and southern inland and coastal people contributed to the formation of the gene pool of the Zhuang and Dong people.

Genome-wide phylogeography reveals cryptic speciation in the circumglobal planktonic calcifier Limacina bulimoides.

Little is known about when and how planktonic species arise and persist in the open ocean without apparent dispersal barriers. Pteropods are planktonic snails with thin shells susceptible to dissolution that are used as bio-indicators of ocean acidification. However, distinct evolutionary units respond to acidification differently, and defining species boundaries is therefore crucial for predicting the impact of changing ocean conditions. In this global population genomic study of the shelled pteropod Limacina bulimoides, we combined genetic (759,000 single nucleotide polymorphisms) and morphometric data from 161 individuals, revealing three major genetic lineages (FST = 0.29-0.41): an "Atlantic lineage" sampled across the Atlantic, an "Indo-Pacific lineage" sampled in the North Pacific and Indian Ocean, and a "Pacific lineage" sampled in the North and South Pacific. A time-calibrated phylogeny suggests that the lineages diverged about 1 million years ago, with estimated effective population size remaining high (~10 million) throughout Pleistocene glacial cycles. We do not observe any signatures of recent hybridization, even in areas of sympatry in the North Pacific. While the lineages are reproductively isolated, they are morphologically cryptic, with overlapping shell shape and shell colour distributions. Despite showing that the circumglobal L. bulimoides consists of multiple species with smaller ranges than initially thought, we found that these pteropods still possess high levels of genetic variability. Our study adds to the growing evidence that speciation is often overlooked in the open ocean, and suggests the presence of distinct biological species within many other currently defined circumglobal planktonic species.

Genome-wide SNPs in the spiny lobster Panulirus homarus reveal a hybrid origin for its subspecies.

BACKGROUND: Evolutionary divergence and speciation often occur at a slower rate in the marine realm due to the higher potential for long-distance reproductive interaction through larval dispersal. One common evolutionary pattern in the Indo-Pacific, is divergence of populations and species at the peripheries of widely-distributed organisms. However, the evolutionary and demographic histories of such divergence are yet to be well understood. Here we address these issues by coupling genome-wide SNP data with mitochondrial DNA sequences to test the patterns of genetic divergence and possible secondary contact among geographically distant populations of the highly valuable spiny lobster Panulirus homarus species complex, distributed widely through the Indo-Pacific, from South Africa to the Marquesas Islands. RESULT: After stringent filtering, 2020 SNPs were used for population genetic and demographic analyses, revealing strong regional structure (FST = 0.148, P < 0001), superficially in accordance with previous analyses. However, detailed demographic analyses supported a much more complex evolutionary history of these populations, including a hybrid origin of a North-West Indian Ocean (NWIO) population, which has previously been discriminated morphologically, but not genetically. The best-supported demographic models suggested that the current genetic relationships among populations were due to a complex series of past divergences followed by asymmetric migration in more recent times. CONCLUSION: Overall, this study suggests that alternating periods of marine divergence and gene flow have driven the current genetic patterns observed in this lobster and may help explain the observed wider patterns of marine species diversity in the Indo-Pacific.

Historical Introgression from Wild Relatives Enhanced Climatic Adaptation and Resistance to Pneumonia in Sheep.

How animals, particularly livestock, adapt to various climates and environments over short evolutionary time is of fundamental biological interest. Further, understanding the genetic mechanisms of adaptation in indigenous livestock populations is important for designing appropriate breeding programs to cope with the impacts of changing climate. Here, we conducted a comprehensive genomic analysis of diversity, interspecies introgression, and climate-mediated selective signatures in a global sample of sheep and their wild relatives. By examining 600K and 50K genome-wide single nucleotide polymorphism data from 3,447 samples representing 111 domestic sheep populations and 403 samples from all their seven wild relatives (argali, Asiatic mouflon, European mouflon, urial, snow sheep, bighorn, and thinhorn sheep), coupled with 88 whole-genome sequences, we detected clear signals of common introgression from wild relatives into sympatric domestic populations, thereby increasing their genomic diversities. The introgressions provided beneficial genetic variants in native populations, which were significantly associated with local climatic adaptation. We observed common introgression signals of alleles in olfactory-related genes (e.g., ADCY3 and TRPV1) and the PADI gene family including in particular PADI2, which is associated with antibacterial innate immunity. Further analyses of whole-genome sequences showed that the introgressed alleles in a specific region of PADI2 (chr2: 248,302,667-248,306,614) correlate with resistance to pneumonia. We conclude that wild introgression enhanced climatic adaptation and resistance to pneumonia in sheep. This has enabled them to adapt to varying climatic and environmental conditions after domestication.

Genomic analysis of a reef-building coral, Acropora digitifera, reveals complex population structure and a migration network in the Nansei Islands, Japan.

Understanding the structure and connectivity of coral populations is fundamental for developing marine conservation policies, especially in patchy environments such as archipelagos. The Nansei Islands, extending more than 1000 km in southwestern Japan, are characterized by high levels of biodiversity and endemism, supported by coral reefs, which make this region ideal for assessing genetic attributes of coral populations. In this study, we conducted population genomic analyses based on genome-wide, single-nucleotide polymorphisms (SNPs) of Acropora digitifera, a common species in the Nansei Islands. By merging newly obtained genome resequencing data with previously published data, we identified more than 4 million genome-wide SNPs in 303 colonies collected at 22 locations, with sequencing coverage ranging from 3.91× to 27.41×. While population structure analyses revealed genetic similarities between the southernmost and northernmost locations, separated by >1000 km, several subpopulations in intermediate locations suggested limited genetic admixture, indicating conflicting migration tendencies in the Nansei Islands. Although migration networks revealed a general tendency of northward migration along the Kuroshio Current, a substantial amount of southward migration was also detected, indicating important contributions of minor ocean currents to coral larval dispersal. Moreover, heterogeneity in the transition of effective population sizes among locations suggests different histories for individual subpopulations. The unexpected complexity of both past and present population dynamics in the Nansei Islands implies that heterogeneity of ocean currents and local environments, past and present, have influenced the population structure of this species, and similar unexpected population complexities may be expected for other marine species with similar reproductive modes.

Estimation of Genetic Correlation via Linkage Disequilibrium Score Regression and Genomic Restricted Maximum Likelihood.

Genetic correlation is a key population parameter that describes the shared genetic architecture of complex traits and diseases. It can be estimated by current state-of-art methods, i.e., linkage disequilibrium score regression (LDSC) and genomic restricted maximum likelihood (GREML). The massively reduced computing burden of LDSC compared to GREML makes it an attractive tool, although the accuracy (i.e., magnitude of standard errors) of LDSC estimates has not been thoroughly studied. In simulation, we show that the accuracy of GREML is generally higher than that of LDSC. When there is genetic heterogeneity between the actual sample and reference data from which LD scores are estimated, the accuracy of LDSC decreases further. In real data analyses estimating the genetic correlation between schizophrenia (SCZ) and body mass index, we show that GREML estimates based on ∼150,000 individuals give a higher accuracy than LDSC estimates based on ∼400,000 individuals (from combined meta-data). A GREML genomic partitioning analysis reveals that the genetic correlation between SCZ and height is significantly negative for regulatory regions, which whole genome or LDSC approach has less power to detect. We conclude that LDSC estimates should be carefully interpreted as there can be uncertainty about homogeneity among combined meta-datasets. We suggest that any interesting findings from massive LDSC analysis for a large number of complex traits should be followed up, where possible, with more detailed analyses with GREML methods, even if sample sizes are lesser.

The wall lizards of the Balkan peninsula: Tackling questions at the interface of phylogenomics and population genomics.

Wall lizards of the genus Podarcis (Sauria, Lacertidae) are the predominant reptile group in southern Europe, including 24 recognized species. Mitochondrial DNA data have shown that, with the exception of P. muralis, the Podarcis species distributed in the Balkan peninsula form a species group that is further sub-divided into two subgroups: the one of "P. tauricus" consisting of P. tauricus, P. milensis, P. gaigeae, and P. melisellensis, and the other of "P. erhardii" comprising P. erhardii, P. levendis, P. cretensis, and P. peloponnesiacus. In an attempt to explore the Balkan Podarcis phylogenomic relationships, assess the levels of genetic structure and to re-evaluate the number of extant species, we employed phylogenomic and admixture approaches on ddRADseq (double digested Restriction site Associated DNA sequencing) genomic data. With this efficient Next Generation Sequencing approach, we were able to obtain a large number of genomic loci randomly distributed throughout the genome and use them to resolve the previously obscure phylogenetic relationships among the different Podarcis species distributed in the Balkans. The obtained phylogenomic relationships support the monophyly of both aforementioned subgroups and revealed several divergent lineages within each subgroup, stressing the need for taxonomic re-evaluation of Podarcis' species in Balkans. The phylogenomic trees and the species delimitation analyses confirmed all recently recognized species (P. levendis, P. cretensis, and P. ionicus) and showed the presence of at least two more species, one in P. erhardii and the other in P. peloponnesiacus.

Phylogenomic resolution of the Ceratitis FARQ complex (Diptera: Tephritidae).

The Ceratitis FARQ complex (formerly FAR complex) includes four frugivorous tephritids, Ceratitis fasciventris, C. anonae, C. rosa and C. quilicii, the latter two causing important agricultural losses in Africa. Although FARQ species can be identified on the basis of subtle morphological differences, they cannot be resolved as monophyletic when trying phylogenetic tree reconstructions based on mitochondrial or nuclear gene fragments except for microsatellites. In this study, we used mitogenome and genome-wide SNPs to investigate the phylogenetic relationship within the complex as well as between all four Ceratitis subgenera. The analysis of 13 species supported the monophyly of the Ceratitis subgenera Ceratitis, Ceratalaspis, Pardalaspis, and recovered Pterandrus as paraphyletic but could not properly resolve species within the FARQ complex. Conversely, gene and species tree reconstructions based on 785,484 genome-wide SNPs could consistently resolve the FARQ taxa and provide insights into their phylogenetic relationships. Gene flow was detected by TreeMix analysis from C. quilicii to C. fasciventris, suggesting the existence of introgression events in the FARQ complex. Our results suggest that genome-wide SNPs represent a suitable tool for the molecular diagnosis of FARQ species and could possibly be used to develop rapid diagnostic methods or to trace the origins of intercepted samples.

Reconstruction of the Austronesian Diaspora in the Era of Genomics.

The Austronesian Diaspora is a 5,000-year account of how a small group of Taiwanese farmers expanded to occupy territories reaching halfway around the world. Reconstructing their detailed history has spawned many academic contests across many disciplines. An outline orthodox version has eventually emerged but still leaves many unanswered questions. The remarkable power of whole-genome technology has now been applied to people across the entire region. This review gives an account of this era of genetic investigation and discusses its many achievements, including revelation in detail of many unexpected patterns of population movement and the significance of this information for medical genetics.

Genome-wide SNP discovery and genetic diversity evaluation of Liangshan cattle in China.

Liangshan cattle are a very small indigenous breed with adult weight of less than 300 Kg and have been mainly distributed in the Liangshan Yi Autonomous Prefecture of Southwestern Sichuan, China. Due to its long-term adaptation to local environments, Liangshan cattle is a valuable genetic resource and should be paid with more attentions. However, the genetic diversity of Liangshan cattle have not been specifically investigated yet, which would be required when designing the appropriate conservation and utilization programs. In this study, we successfully employed the restriction-site-associated DNA sequencing (RAD-seq) approach to explore a total of 84,854 genome-wide and high-confidence SNPs of Liangshan cattle. All these SNPs were evenly distributed through all chromosomes with an average of 98 SNPs per 1-Mb region. The nucleotide diversity, expected heterozygosity, polymorphism information content of Liangshan cattle were 0.227, 0.223 and 0.183, respectively. Furthermore, there was no obvious difference on the genetic diversity among the three studied geographical populations. In conclusion, we provided a list of SNPs that could be used in the follow-up studies for Liangshan cattle and revealed a relatively high genetic variation in this gene pool.

Using spatial genetics to quantify mosquito dispersal for control programs.

BACKGROUND: Hundreds of millions of people get a mosquito-borne disease every year and nearly one million die. Transmission of these infections is primarily tackled through the control of mosquito vectors. The accurate quantification of mosquito dispersal is critical for the design and optimization of vector control programs, yet the measurement of dispersal using traditional mark-release-recapture (MRR) methods is logistically challenging and often unrepresentative of an insect's true behavior. Using Aedes aegypti (a major arboviral vector) as a model and two study sites in Singapore, we show how mosquito dispersal can be characterized by the spatial analyses of genetic relatedness among individuals sampled over a short time span without interruption of their natural behaviors. RESULTS: Using simple oviposition traps, we captured adult female Ae. aegypti across high-rise apartment blocks and genotyped them using genome-wide SNP markers. We developed a methodology that produces a dispersal kernel for distance which results from one generation of successful breeding (effective dispersal), using the distance separating full siblings and 2nd- and 3rd-degree relatives (close kin). The estimated dispersal distance kernel was exponential (Laplacian), with a mean dispersal distance (and dispersal kernel spread σ) of 45.2 m (95% CI 39.7-51.3 m), and 10% probability of a dispersal > 100 m (95% CI 92-117 m). Our genetically derived estimates matched the parametrized dispersal kernels from previous MRR experiments. If few close kin are captured, a conventional genetic isolation-by-distance analysis can be used, as it can produce σ estimates congruent with the close-kin method if effective population density is accurately estimated. Genetic patch size, estimated by spatial autocorrelation analysis, reflects the spatial extent of the dispersal kernel "tail" that influences, for example, the critical radii of release zones and the speed of Wolbachia spread in mosquito replacement programs. CONCLUSIONS: We demonstrate that spatial genetics can provide a robust characterization of mosquito dispersal. With the decreasing cost of next-generation sequencing, the production of spatial genetic data is increasingly accessible. Given the challenges of conventional MRR methods, and the importance of quantified dispersal in operational vector control decisions, we recommend genetic-based dispersal characterization as the more desirable means of parameterization.

Dissecting the Pre-Columbian Genomic Ancestry of Native Americans along the Andes-Amazonia Divide.

Extensive European and African admixture coupled with loss of Amerindian lineages makes the reconstruction of pre-Columbian history of Native Americans based on present-day genomes extremely challenging. Still open questions remain about the dispersals that occurred throughout the continent after the initial peopling from the Beringia, especially concerning the number and dynamics of diffusions into South America. Indeed, if environmental and historical factors contributed to shape distinct gene pools in the Andes and Amazonia, the origins of this East-West genetic structure and the extension of further interactions between populations residing along this divide are still not well understood. To this end, we generated new high-resolution genome-wide data for 229 individuals representative of one Central and ten South Amerindian ethnic groups from Mexico, Peru, Bolivia, and Argentina. Low levels of European and African admixture in the sampled individuals allowed the application of fine-scale haplotype-based methods and demographic modeling approaches. These analyses revealed highly specific Native American genetic ancestries and great intragroup homogeneity, along with limited traces of gene flow mainly from the Andes into Peruvian Amazonians. Substantial amount of genetic drift differentially experienced by the considered populations underlined distinct patterns of recent inbreeding or prolonged isolation. Overall, our results support the hypothesis that all non-Andean South Americans are compatible with descending from a common lineage, while we found low support for common Mesoamerican ancestors of both Andeans and other South American groups. These findings suggest extensive back-migrations into Central America from non-Andean sources or conceal distinct peopling events into the Southern Continent.

ddRAD analyses reveal a credible phylogenetic relationship of the four main genera of Bambusa-Dendrocalamus-Gigantochloa complex (Poaceae: Bambusoideae).

The Bambusa-Dendrocalamus-Gigantochloa complex (BDG complex) is the most diversified and phylogenetically recalcitrant group of the paleotropical woody bamboos. Species of this complex occur in tropical and subtropical Asia and most of them are of great economic, cultural and ecological value. The lack of resolution achieved through the analyses of previous molecular datasets has long confounded its phylogenetic estimation and generic delimitation. Here, we adopted a ddRAD-seq strategy to investigate phylogenetic relationships of the four main genera (Bambusa, Dendrocalamus, Gigantochloa, and Melocalamus) in the BDG complex. A total of 102 species were sampled, and SNP data were generated. Both MP and ML analyses of the ddRAD-seq data resulted in a well-resolved topology with Gigantochloa and Melocalamus confirmed as monophyletic, and Melocalamus resolved as sister to the rest of the complex. Bambusa and Dendrocalamus were both resolved as paraphyletic. The phylogenetic relationships were mostly supported by morphological evidence including characters of the branch complement, rachilla, lodicules, filaments and stigma. We also generated and assembled complete plastid genomes of 48 representative species. There were conflicts between the plastome and the ddRAD topologies. Our study demonstrated that RAD-seq can be used to reconstruct evolutionary history of lineages such as the bamboos where ancient hybridization and polyploidy play a significant role. The four genera of the BDG complex have a complex evolutionary history which is likely a product of ancient introgression events.

Chad Genetic Diversity Reveals an African History Marked by Multiple Holocene Eurasian Migrations.

Understanding human genetic diversity in Africa is important for interpreting the evolution of all humans, yet vast regions in Africa, such as Chad, remain genetically poorly investigated. Here, we use genotype data from 480 samples from Chad, the Near East, and southern Europe, as well as whole-genome sequencing from 19 of them, to show that many populations today derive their genomes from ancient African-Eurasian admixtures. We found evidence of early Eurasian backflow to Africa in people speaking the unclassified isolate Laal language in southern Chad and estimate from linkage-disequilibrium decay that this occurred 4,750-7,200 years ago. It brought to Africa a Y chromosome lineage (R1b-V88) whose closest relatives are widespread in present-day Eurasia; we estimate from sequence data that the Chad R1b-V88 Y chromosomes coalesced 5,700-7,300 years ago. This migration could thus have originated among Near Eastern farmers during the African Humid Period. We also found that the previously documented Eurasian backflow into Africa, which occurred ∼3,000 years ago and was thought to be mostly limited to East Africa, had a more westward impact affecting populations in northern Chad, such as the Toubou, who have 20%-30% Eurasian ancestry today. We observed a decline in heterozygosity in admixed Africans and found that the Eurasian admixture can bias inferences on their coalescent history and confound genetic signals from adaptation and archaic introgression.

The genetic legacy of the Yaghnobis: A witness of an ancient Eurasian ancestry in the historically reshuffled central Asian gene pool.

OBJECTIVES: The Yaghnobis are an ethno-linguistic minority historically settled along the Yaghnob River in the Upper-Zarafshan Valley in Tajikistan. They speak a language of Old Sogdian origin, which is the only present-day witness of the Lingua Franca used along the Silk Road in Late Antiquity. The aim of this study was to reconstruct the genetic history of this community in order to shed light on its isolation and genetic ancestry within the Euro-Asiatic context. MATERIALS AND METHODS: A total of 100 DNA samples were collected in the Yaghnob and Matcha Valleys during several expeditions and their mitochondrial, Y-chromosome and autosomal genome-wide variation were compared with that from a large set of modern and ancient Euro-Asiatic samples. RESULTS: Findings from uniparental markers highlighted the long-term isolation of the Yaghnobis. Mitochondrial DNA ancestry traced an ancient link with Middle Eastern populations, whereas Y-chromosome legacy showed more tight relationships with Central Asians. Admixture, outgroup-f3, and D-statistics computed on autosomal variation corroborated Y-chromosome evidence, pointing respectively to low Anatolian Neolithic and high Steppe ancestry proportions in Yaghnobis, and to their closer affinity with Tajiks than to Iranians. DISCUSSION: Although the Yaghnobis do not show evident signs of recent admixture, they could be considered a modern proxy for the source of gene flow for many Central Asian and Middle Eastern groups. Accordingly, they seem to retain a peculiar genomic ancestry probably ascribable to an ancient gene pool originally wide spread across a vast area and subsequently reshuffled by distinct demographic events occurred in Middle East and Central Asia.

Genomic diversity dynamics in conserved chicken populations are revealed by genome-wide SNPs.

BACKGROUND: Maintaining maximum genetic diversity and preserving breed viability in conserved populations necessitates the rigorous evaluation of conservation schemes. Three chicken breeds (Baier Yellow Chicken (BEC), Beijing You Chicken (BYC) and Langshan Chicken (LSC)) are currently in conservation programs in China. Changes in genetic diversity were measured by heterozygosity, genomic inbreeding coefficients, and autozygosity, using estimates derived from runs of homozygosity (ROH) that were identified using SNPs. RESULTS: Ninety DNA samples were collected from three generations for each breed. In the most recent generation, the highest genetic diversity was observed in LSC, followed by BEC and BYC. Inbreeding coefficients based on ROH for the three breeds declined slightly between the first and middle generations, and then rapidly increased. No inbreeding coefficients exceeded 0.1. Population structure assessments using neighbor-joining tree analysis, principal components analysis, and STRUCTURE analysis indicated that no genetic differentiation existed within breeds. LD decay and ROH at different cut-off lengths were used to identify traces left by recent or ancient inbreeding. Few long ROH were identified, indicating that inbreeding has been largely avoided with the current conservation strategy. The observed losses in genetic diversity and occurrences of inbreeding might be consequences of sub-optimal effective population sizes. CONCLUSIONS: The conserved Chinese chicken populations have high genomic diversity under the current conservation program (R: F). Furthermore, this study highlights the need to monitor dynamic changes in genetic diversity in conserved breeds over successive generations. Our research provides new insights into genetic diversity dynamics in conserved populations, and lays a solid foundation for improving conservation schemes.

Recent Historical Migrations Have Shaped the Gene Pool of Arabs and Berbers in North Africa.

North Africa is characterized by its diverse cultural and linguistic groups and its genetic heterogeneity. Genomic data has shown an amalgam of components mixed since pre-Holocean times. Though no differences have been found in uniparental and classical markers between Berbers and Arabs, the two main ethnic groups in the region, the scanty genomic data available have highlighted the singularity of Berbers. We characterize the genetic heterogeneity of North African groups, focusing on the putative differences of Berbers and Arabs, and estimate migration dates. We analyze genome-wide autosomal data in five Berber and six Arab groups, and compare them to Middle Easterns, sub-Saharans, and Europeans. Haplotype-based methods show a lack of correlation between geographical and genetic populations, and a high degree of genetic heterogeneity, without strong differences between Berbers and Arabs. Berbers enclose genetically diverse groups, from isolated endogamous groups with high autochthonous component frequencies, large homozygosity runs and low effective population sizes, to admixed groups with high frequencies of sub-Saharan and Middle Eastern components. Admixture time estimates show a complex pattern of recent historical migrations, with a peak around the 7th century C.E. coincident with the Arabization of the region; sub-Saharan migrations since the 1st century B.C. in agreement with Roman slave trade; and a strong migration in the 17th century C.E., coincident with a huge impact of the trans-Atlantic and trans-Saharan trade of sub-Saharan slaves in the Modern Era. The genetic complexity found should be taken into account when selecting reference groups in population genetics and biomedical studies.

Genetic relatedness of indigenous ethnic groups in northern Borneo to neighboring populations from Southeast Asia, as inferred from genome-wide SNP data.

The region of northern Borneo is home to the current state of Sabah, Malaysia. It is located closest to the southern Philippine islands and may have served as a viaduct for ancient human migration onto or off of Borneo Island. In this study, five indigenous ethnic groups from Sabah were subjected to genome-wide SNP genotyping. These individuals represent the "North Borneo"-speaking group of the great Austronesian family. They have traditionally resided in the inland region of Sabah. The dataset was merged with public datasets, and the genetic relatedness of these groups to neighboring populations from the islands of Southeast Asia, mainland Southeast Asia and southern China was inferred. Genetic structure analysis revealed that these groups formed a genetic cluster that was independent of the clusters of neighboring populations. Additionally, these groups exhibited near-absolute proportions of a genetic component that is also common among Austronesians from Taiwan and the Philippines. They showed no genetic admixture with Austro-Melanesian populations. Furthermore, phylogenetic analysis showed that they are closely related to non-Austro-Melansian Filipinos as well as to Taiwan natives but are distantly related to populations from mainland Southeast Asia. Relatively lower heterozygosity and higher pairwise genetic differentiation index (FST ) values than those of nearby populations indicate that these groups might have experienced genetic drift in the past, resulting in their differentiation from other Austronesians. Subsequent formal testing suggested that these populations have received no gene flow from neighboring populations. Taken together, these results imply that the indigenous ethnic groups of northern Borneo shared a common ancestor with Taiwan natives and non-Austro-Melanesian Filipinos and then isolated themselves on the inland of Sabah. This isolation presumably led to no admixture with other populations, and these individuals therefore underwent strong genetic differentiation. This report contributes to addressing the paucity of genetic data on representatives from this strategic region of ancient human migration event(s).

Clonal Lineages of Erysipelothrix rhusiopathiae Responsible for Acute Swine Erysipelas in Japan Identified by Using Genome-Wide Single-Nucleotide Polymorphism Analysis.

Erysipelothrix rhusiopathiae causes swine erysipelas, an important infectious disease in the swine industry. In Japan, the incidence of acute swine erysipelas due to E. rhusiopathiae serovar 1a has recently increased markedly. To study the genetic relatedness of the strains from the recent cases, we analyzed 34 E. rhusiopathiae serovar 1a swine isolates collected between 1990 and 2011 and further investigated the possible association of the live Koganei 65-0.15 vaccine strain (serovar 1a) with the increase in cases. Pulsed-field gel electrophoresis analysis revealed no marked variation among the isolates; however, sequencing analysis of a hypervariable region in the surface-protective antigen A gene (spaA) revealed that the strains isolated after 2007 exhibited the same spaA genotype and could be differentiated from older strains. Phylogenetic analysis based on genome-wide single-nucleotide polymorphisms (SNPs) revealed that the Japanese strains examined were closely related, showing a relatively small number of SNPs among them. The strains were classified into four major lineages, with Koganei 65-0.15 (lineage III) being phylogenetically separated from the other three lineages. The strains isolated after 2007 and the two older strains constituted one major lineage (lineage IV) with a specific spaA genotype (M203/I257-SpaA), while the recent isolates were further divided into two geographic groups. The remaining older isolates belonged to either lineage I, with the I203/L257-SpaA type, or lineage II, with the I203/I257-SpaA type. These results indicate that the recent increased incidence of acute swine erysipelas in Japan is associated with two sublineages of lineage IV, which have independently evolved in two different geographic regions.IMPORTANCE Using large-scale whole-genome sequence data from Erysipelothrix rhusiopathiae isolates from a wide range of hosts and geographic origins, a recent study clarified the existence of three distinct clades (clades 1, 2, and 3) that are found across multiple continents and host species, representing both livestock and wildlife, and an "intermediate" clade between clade 2 and the dominant clade 3 within the species. In this study, we found that the E. rhusiopathiae Japanese strains examined exhibited remarkably low levels of genetic diversity and confirmed that all of the Japanese and Chinese swine isolates examined in this study belong to clonal lineages within the intermediate clade. We report that spaA genotyping of E. rhusiopathiae strains is a practical alternative to whole-genome sequencing analysis of the E. rhusiopathiae isolates from eastern Asian countries.

Genome-wide SNP analysis reveals population structure and demographic history of the ryukyu islanders in the southern part of the Japanese archipelago.

The Ryukyu Islands are located to the southwest of the Japanese archipelago. Archaeological evidence has revealed the existence of prehistoric cultural differentiation between the northern Ryukyu islands of Amami and Okinawa, and the southern Ryukyu islands of Miyako and Yaeyama. To examine a genetic subdivision in the Ryukyu Islands, we conducted genome-wide single nucleotide polymorphism typing of inhabitants from the Okinawa Islands, the Miyako Islands, and the Yaeyama Islands. Principal component and cluster analyses revealed genetic differentiation among the island groups, especially between Okinawa and Miyako. No genetic affinity was observed between aboriginal Taiwanese and any of the Ryukyu populations. The genetic differentiation observed between the inhabitants of the Okinawa Islands and the Miyako Islands is likely to have arisen due to genetic drift rather than admixture with people from neighboring regions. Based on the observed genetic differences, the divergence time between the inhabitants of Okinawa and Miyako islands was dated to the Holocene. These findings suggest that the Pleistocene inhabitants, whose bones have been found on the southern Ryukyu Islands, did not make a major genetic contribution, if any, to the present-day inhabitants of the southern Ryukyu Islands.