Great work @dyno I prefer you contributing instead of joking around

Most of the peeps here got's no charisma man. 

NP

Maybe as slaves by pirates 

The geographic distribution of complete mtDNA sequences affiliated with the Malagasy motif (B4a1a1b haplogroup), its precursor haplogroup B4a1a1 (the Polynesian motif), and related subclades (supplementary table S1, Supplementary Material online) shows a clear phylogeographic distribution pattern (fig. 1 ). On the one hand, the Polynesian motif and all of its subclades, with the exception of the Malagasy motif subclade (B4a1a1b), are only found east of Wallace’s line (with minor exceptions). On the other hand, the Malagasy motif is only detected in the WIO (Madagascar, and now Somalia and Yemen).

FIG. 1.

Open in new tabDownload slide

—Schematic representations of the B4a1a1 phylogeny based on complete mtDNA sequences. Subclades are represented by triangles. Subclades are colored according to their geographic origin, as shown on the map. (A) B4a1a1 tree. (B) B4a1a1b tree. (C) Geographic distribution of the subclades. KA, thousand years ago.

 

 

 

In Island Southeast Asia, the Banjar, who have been shown to be linguistically (Adelaar 1989, 2017) and genetically (Brucato et al. 2017) the descendants of the ancestors of the Asian background found in Madagascar and the Comoros, carry no maternal lineage affiliated with the Malagasy motif, and only one individual (out of 99) carries a maternal lineage affiliated to the Polynesian motif (supplementary table S3, Supplementary Material online). The other Banjar mitochondrial haplogroups are similar to the observed haplogroup diversity in other populations in Indonesia (Tumonggor et al. 2013; Kusuma et al. 2015).

The phylogenetic tree based on the B4a1a1b high-quality mitogenomes available (fig. 1 and supplementary fig. S1, Supplementary Material online) shows the absence of geographic structuring among these sequences, together with their very low diversity. The two sequences from Yemen and Madagascar are identical to the basal B4a1a1b lineage and the other sequences carry only one additional mutation each, mostly from mutational hotspots in the hypervariable noncoding region (16291T, 16357C, 16368C, and 16222T; Lott et al. 2013). This low diversity is in agreement with the recent coalescence age estimated for haplogroup B4a1a1b of 1,500–1,800 years BP (table 1 and supplementary table S4, Supplementary Material online), a narrower window than previously estimated (1,200–16,800 years BP; Razafindrazaka et al. 2010). This new age estimate predates, but is consistent with, the earliest Austronesian settlement of the East African offshore islands (Madagascar and Comoros), which is dated to 800–1,200 years BP from recent archaeological (Crowther et al. 2016) and genetic studies (Brucato et al. 2016, 2018).

Genome-Wide Analysis

To determine whether genome-wide genetic data can shed further light on the history of individuals carrying the Malagasy motif, the genetic ancestries present in the genome-wide data set were decomposed with ADMIXTURE v.1.3 (fig. 2). The lowest cross-validation values were obtained with 26 ancestries (supplementary fig. S2, Supplementary Material online). For clarity, using a reduced admixture data set including representative populations from each region and the population of interest, K = 4 also displays the main African, European, South Asia, and Southeast Asian ancestries (supplementary fig. S3, Supplementary Material online). We observed that the Somali (S25) and two Yemeni (Y270 and Y115) individuals have relatively similar ancestry patterns to their local populations, including a dual Asian ancestry reflected by a South Asian (India) and Southeast Asian (Indonesia) components (fig. 2 and supplementary fig. S3, Supplementary Material online). Moreover, Yemenese and Somali do not have significantly more Asian ancestry (0.87% and 1.18%, respectively) than paired more inland populations (supplementary fig. S4, Supplementary Material online). The two Yemeni individuals have one major component (pinky red gradient: 65–80%) shared with other Middle East populations. Other minor ancestries are shared with populations from the horn of Africa (gray gradient: 13%), sub-Saharan Bantu speakers (green gradient: 7–10%), and Europeans (blue gradient 2–5%). The Asian component (yellow gradient) represents just a small component: from 1% in Y115—similar to those in Yemeni populations, up to 10% in Y270. The Somali individual has one major component (green gradient: 60%) shared with Bantu speakers, and a smaller component shared with populations from the horn of Africa (gray gradient: 25%), the Middle East (pink red gradient: 5%), and very minor components from Asia (yellow gradient <1%).

FIG. 2.

Open in new tabDownload slide

—ADMIXTURE analysis (K = 26) for 15 population groups and the two Yemeni and one Somali individuals carrying the Malagasy motif B4a1a1b maternal lineage. Each colored line represents a sampled population whose genetic background can be decomposed into 26 genetic components.

 

 

 

The admixture scenario from f3 statistics (supplementary tables S5–S7, Supplementary Material online) shows that the two Yemeni (Y115 and Y270) and the Somali (S25) samples cannot be modeled as any kind of simple mixture, using any Comoros, Madagascar, or Banjar groups as one population source and one of their respective parental or neighboring populations as the other source. To explore their admixture further, the admixture history and Asian component of these individuals was subjected to PCAdmix analysis (table 2). This indicated that two individuals (Y115 and S25) have significantly more Malagasy fragments than are observed in their respective parental populations (Z-score = 1,572 and 5,556, respectively), thus suggesting a likely Malagasy origin of their Asian ancestry. The other Yemeni individual (Y270) carries more Indonesian fragments than observed in its parental population (i.e., Asian and African Bantu fragments are not associated) (Z-score = 8,368), suggesting an origin of its Asian component from a population with limited African Bantu admixture, and thus excluding the Malagasy as a source (in average, Malagasy populations analyzed so far have around 60% African Bantu ancestry and 40% Asian ancestry; Pierron et al. 2014, 2017; Brucato et al. 2016, 2018).

The date of inheritance of the Asian component found in the Yemeni and Somali Malagasy motif carriers can be estimated using an exponential decay function to explore the decline in total genome-wide Asian ancestry. Figure 3 shows that the two individuals with <1% Asian SNPs (Y115 and S25) likely had an Asian ancestor before 5–7 generations ago (lower bound prior to 250–190 years BP) and probably much older. As the Asian ancestry in these individuals is similar to other individuals in the local population who do not carry the Malagasy motif, the inheritance is possibly far earlier. In contrast, the individual with 10% Asian SNPs (Y270) likely had an Asian ancestor 3–4 generations ago (160–130 years BP), and because his total Asian ancestry is much higher than others in his local population (which has similar Asian ancestry to neighboring populations, supplementary fig. S4, Supplementary Material online), the inheritance date is probably relatively close to this date estimate. There may therefore be multiple reasons why the Malagasy motif was inherited by these three individuals.

FIG. 3.

Open in new tabDownload slide

—Decay curves showing expected autosomal Asian ancestry at each generation of intermarriage to partners with no Asian ancestry. Decay rates for individuals with Indonesian (green) and Malagasy (blue) origins. Horizontal lines mark observed Asian ancestry for the Somali and Yemeni populations (orange), Indonesian descendant Yemeni 270 (green), and Malagasy descendants Somali 25 and Yemeni 115 (blue).

 

 

 

Surprisingly, these results suggest that despite the sharing of a clear Asian component, from their maternal lineage and autosomal DNA, the geographic origin and timing of admixture of this Asian component may be different among the three Yemeni and Somali individuals, revealing a more complex dynamic of Austronesian gene flow in the Western Indian Ocean region.

Discussion

Our results bring new information on the Austronesian dispersal westward across the Indian Ocean. We show that the key maternal lineage B4a1a1b, the so-called Malagasy motif, has a coalescence age (1,500–1,800 years BP) either predating or very early during the period of Austronesian arrivals to the East African offshore islands of Madagascar and the Comoros (table 1 and supplementary table S4, Supplementary Material online and fig. 1). This age estimate gives some temporal support for the emergence of the Malagasy motif in Island Southeast Asia (also in agreement with the fact that demographic expansion predates the geographic expansion). Although it remains undetected in the Banjar population, leaving its population source in ISEA unknown (supplementary table S3, Supplementary Material online), its recent contribution detected in this study (nineteenth century) from a population without African admixture implies that it should probably still exist in some location in ISEA.

The absence of B4a1a1b in the Banjar is surprising as it represents the most frequent Asian maternal lineage in Madagascar (22%, Razafindrazaka et al. 2010; Cox et al. 2012), and considering the cultural, linguistic and genetic links between the Banjar and Madagascar/Comoros populations (Adelaar 1989, 2017; Beaujard 2012a, 2012b; Brucato et al. 2016, 2018), this lineage was expected to be present in the Asian parental population (the Banjar). Its absence may suggests that the main Austronesian maternal lineage (B4a1a1b) found in the WIO may not have originated from the Banjar population. However, an origin of the Malagasy motif in situ in Madagascar or the Comoros after the arrival of the Polynesian motif carriers is also unlikely (Razafindrazaka et al. 2010), as it requires that 1) two mutations (1,473 and 3,423) arose in Madagascar in the last 1,300–1,000 years (Brucato et al. 2018), 2) diffused into population groups across the entire island, and 3) the Malagasy motif precursor (the Polynesian motif), later disappeared from these populations.

This “missing” lineage in Southeast Borneo 1) may have been lost in Indonesia due to genetic drift in small island populations or 2) may remain undetected due to the lower sampling coverage in Borneo compared with other studied regions. As suggested previously (Razafindrazaka et al. 2010; Cox et al. 2012; Kusuma et al. 2015), the motif may have emerged in eastern Indonesia—where the Polynesian motif, its immediate precursor, is geographically restricted—and later brought into the Banjar gene pool. This may have occurred during long distance trading activities and admixture among ISEA regional populations favored by the emergence Hindu Malay Kingdoms, such as Śrīvijaya (sixth to thirteenth centuries), which developed the Banjarmasin trading post in Southeast Borneo (Ras 1968; Beaujard 2012a; Brucato et al. 2017). This eastern Indonesian input into the Southeast Borneo Banjar finds some support in the fact that the Banjar are one of the rare populations in western Indonesia, along with the sea nomad Bajo from Borneo, to carry the Polynesian motif (supplementary table S3, Supplementary Material online, Tumonggor et al. 2013; Kusuma et al. 2015, 2017), demonstrating their previous connection with eastern Indonesians.

In the Indian Ocean rim west of Island Southeast Asia (Mainland Southeast Asia and the Indian subcontinent), we note the absence of the Polynesian motif or any of its subclades (fig. 1). This suggests that the Austronesian population that brought the Polynesian Motif subclade B4a1a1b to the WIO (e.g., to Madagascar) did not leave detectable genetic traces on the northern Indian Ocean rim, in agreement with an absence of Indonesian gene flow to this region as inferred from recent genome-wide data (Brucato et al. 2017). It is still unclear whether this is due to Indonesian populations restricting their interactions in these regions to trading and/or cultural activities, or because Indonesian traders used a more direct route across the Indian Ocean to cover the 7,500 km between Indonesia and Madagascar and the Comoros, which is possible based on ocean current and monsoon weather patterns (Fitzpatrick and Callaghan 2008). Both scenarios merit further study.

For the first time, we observe the Malagasy motif outside of Madagascar, in one individual from East Africa (Somalia) and two from the South Arabian Peninsula (Yemen). To date, a long-term Austronesian presence in the Western Indian Ocean region has only been supported on the island territories of Madagascar and the Comoros and is very tenuous on the African continent and Arabian Peninsula, from genetic (Brucato et al. 2017, 2018), historical (Beaujard 2012a, 2012b), and archaeological (Crowther et al. 2016) perspectives. Our results suggest clear, but limited impact, of Austronesian genetic input into these East African and South Arabian regions (0.02%, 3 out of 14,461 individuals from mtDNA Austronesian key marker, up to around 1% when considering the entire Asian autosomal signal), but leaves open the question of their origin.

When considering the Asian component in the autosomal DNA, beyond the fact that no more than 1% of Asian substrate is detected in East Africa/Arabia (fig. 2 and supplementary fig. S3, Supplementary Material online), the three new Malagasy motif carriers displayed a more complex pattern. For two individuals (Y115 and S25), an inheritance from Malagasy individuals before the late eighteenth century and potentially much earlier is likely. For the other individual (Y270), more recent gene flow in the late nineteenth century from a different population source (likely with primary Asian ancestry) seems to be the most parsimonious explanation (figs. 2 and3 and table 2).

This complex pattern of Asian genetic input suggests that these three individuals obtained their Asian genetic ancestry from two different processes or events, via both a primary and secondary Asian input to the East African coast and South Arabian Peninsula. This involved at least two different population sources: an African–Asian admixed population (e.g., potentially from Madagascar) for the oldest admixture event; and a population without African admixture as the source of the most recent admixture event. These two events likely took place during the last millennium with the intensification of the Indian Ocean trading network among Africa, the Middle East and Asia, leading to the exchange of ideas, goods, and people (Beaujard 2012a, 2012b; Brucato et al. 2017), the arrival of the Austronesians in the Western Indian Ocean region around the eighth to thirteenth centuries AD, and the development of the Swahili corridor (from southern Somalia in the north to the Comoros archipelago, Madagascar, and Central Mozambique in the south; Horton 1987; Brucato et al. 2018).

Although the Asian inheritance for two individuals (Y115 and S25) may reflect the development of trading posts by Arab merchants in the East African region, and deportation of African slaves northward to Arabia and South Asia (Hellenthal et al. 2014; Blench 2014; Brucato et al. 2017), we cannot reject an older origin from the early Austronesian period. The Asian ancestry of the other individual (Y270) represents a different pattern with a likely origin in Island Southeast Asia, which might be explained by the emergence of several sultanates in the fifteenth to sixteenth centuries AD in East Africa and in ISEA (Beaujard 2012b). Diaspora communities spread around the Indian Ocean rim may have played an important role to develop trading activities and also to maintain contact with their native land (Beaujard 2012b). The Hadrami community from Hadramawt in Southeast Yemen—to which the Yemeni individual (Y270) belongs—emerge as one plausible cause to reconcile the different geographic and temporal origins of the Asian inheritance found in the three Yemeni and Somali individuals. The Hadrami diaspora, established in the Comoros at the end of the first millennium (Beaujard 2012b) in South Asia (Gujarat in India) and in Island Southeast Asia (Malacca in Malaysia) from the beginning of the second millennium, was involved in trading activities between Africa, Middle East, and Asia. Despite strict intra community marriage rules, admixture with local population was frequent, as the men could choose their bride from among the local African and Asian populations (Boxberger 2002). In the nineteenth century, some economically successful Hadramis (i.e., from Singapore and Hyderabad, and other places of the diaspora) returned to their native land in the Yemen bringing back their admixed family (Manger 2010).

The pattern emerging is that trading activities and diaspora communities may have favored small scale genetic admixture through direct contact between distant regions and populations from the Arabian Peninsula, Africa, and Island Southeast Asia, and that these involved Indonesian populations carrying this peculiar maternal lineage, the Malagasy motif B4a1a1b.

Conclusions

We gathered comparative genome-wide data from previously published studies of populations from Africa, Madagascar, the Middle East, ISEA, Southeast Asia, South Asia, East Asia, and Europe (supplementary table S2, Supplementary Material online). Before analysis, data quality controls were performed using PLINK v.1.9 (Chang et al. 2015): 1) to avoid close relatives, relatedness was measured between all pairs of individuals within each population using an identity-by-descent estimation with an upper threshold of 0.25 (second-degree relatives); 2) SNPs that failed the Hardy–Weinberg exact test (P < 10⁻⁶) in each group were excluded; and 3) samples with an overall call rate <0.99 and individual SNPs with missing rates >0.05 across all samples within each population were excluded. Two data sets were subsequently constituted. For frequency-based analyses (ADMIXTURE, f3-statistics), our data set included 3,480 individuals from 193 populations genotyped for 171,728 SNPs, obtained after pruning for variants in high linkage disequilibrium (LD) with Plink v.1.9 (r² > 0.5; 50 SNP sliding windows). For haplotype-based analyses (PCAdmix), the data set included 411,442 SNPs for 524 individuals, representing 5 metapopulations (Southeast Africa, East Africa, Middle East, South Asia, and ISEA), the Yemeni and Somali populations, and the two Yemeni individuals and the Somali individual of interest. The five metapopulations are composed of 50 individuals with Southeast African Bantu ancestry (randomly selected from Kenyan Luhya, South African Bantu, and Swahili groups), East African (horn of Africa) ancestry (randomly selected from Oromo, Ethiopian from Somalia), Middle Eastern ancestry (randomly selected from Oman and Saudi Arabia), South Asia (randomly selected from Gujarat Brahmin), and Indonesian ancestry (randomly selected from Indonesian Banjar, Samihim, and Malay). Genotypes were then phased with SHAPEIT v.2 (Delaneau et al. 2011) using the 1000 Genomes Project phased data (Delaneau et al. 2014) as a reference panel and the HapMap phase II genetic map.

To address specific questions regarding the ancestries of these three individuals, we performed ADMIXTURE analyses (Alexander et al. 2009), with default settings, for components K = 2 to K = 30. Ten iterations with randomized seeds were run and compiled with CLUMPAK v.1 (Kopelman et al. 2015). The minimum average cross-validation value was used to define the most informative K component (here, K = 26).

To identify potential admixture events, three-population (f3) statistics (Pickrell and Pritchard 2012) were computed using each new sample as recipient and two population sources from the genome-wide data set, with one of the sources representing the Asian input (Comoros, Madagascar, or Banjar) and the other source representing one of their respective parental or neighboring populations. The population trios yielding a Z-score smaller than −2 were considered significantly admixed.

Local ancestry analysis in the three individuals was performed with PCAdmix v.1.050 with the five parental metapopulations (Southeast African Bantu, East African, Middle Eastern, South Asia, and Indonesian). The phased data were screened with LD information so that the probability of common ancestry of each haplotype with each “parental” metapopulation could be defined. The Viterbi algorithm was then used to identify local ancestries of all haplotypes in the two Yemeni and one Somali individuals, as well as in the Yemeni and Somali populations (excluding the three individuals of interest). This allows estimation of the proportion of Malagasy ancestry (defined as Southeast African Bantu and Indonesian associated fragments, based on Malagasy genetic diversity history, Brucato et al. 2016) and Indonesian ancestry (Indonesian fragments isolated from Southeast African fragments) in these three individuals compared with their respective parental populations. A Z-score was computed to estimate the deviation of Malagasy/Asian proportions in the three individuals of interest in comparison to the rest of their population.

Results