Pairwise Population Research Articles

Phylogeography Reveals Routes of Colonization of the Bark Beetle Dendroctonus approximatus Dietz in Mexico

We used mitochondrial DNA (mtDNA) sequence data and allele frequencies at eight microsatellite loci to examine the population genetic structure, estimate the divergence times of distinct lineages, and infer patterns associated with host colonization in populations of the bark beetle Dendroctonus approximatus in Mexico. Two haplotype groups were identified using mtDNA sequences in 71 individuals from 15 populations. The first group was distributed in the Sierra Madre Occidental (SMOc, Western Mexico), with some populations in the Faja Volcánica Transmexicana (Central Mexico), and the second was found in the Sierra Madre Oriental (SMOr, Eastern Mexico), with populations in the Sierra Madre del Sur (Southern Mexico). The estimated split between groups occurred in the late Pleistocene, around 0.195 Mya. Microsatellite allele frequencies revealed high genetic differentiation between pairwise populations, and genetic differentiation values indicated a genetic structure of isolation by distance. Both mtDNA sequence data and microsatellite allele frequencies indicated that D. approximatus had two independent colonization routes in Mexico, one through the SMOc and another along the SMOr. The widespread geographic distribution of D. approximatus in Mexico follows a model of population range expansion of two haplotype groups in which gene flow is restricted by the geographic separation between hosts imposed by physical barriers between populations.

Corrigendum

We wish to bring attention to the following correction to Marsden et al. (2012). The legend for Fig. 3 contained an error with regards to the description of which genetic marker corresponds to the different symbols in the graph. Specifically, the printed legend stated that DLA-DRB1 is represented by white symbols (not black), and microsatellite loci by black symbols (not white). These discrepancies do not influence any of the results or conclusions of our work. A corrected version for the Fig. 3 legend is shown here. Pairwise population genetic distance (Dest) plotted against geographical distance (km) for DLA-DRB1 (black dots) and microsatellite (white dots) loci.

Tracking the Emergence of a New Breed Using 49,034 SNP in Sheep

Domestic animals are unique in that they have been organised into managed populations called breeds. The strength of genetic divergence between breeds may vary dependent on the age of the breed, the scenario under which it emerged and the strength of reproductive isolation it has from other breeds. In this study, we investigated the Gulf Coast Native breed of sheep to determine if it contains lines of animals that are sufficiently divergent to be considered separate breeds. Allele sharing and principal component analysis (PCA) using nearly 50,000 SNP loci revealed a clear genetic division that corresponded with membership of either the Florida or Louisiana Native lines. Subsequent analysis aimed to determine if the strength of the divergence exceeded that found between recognised breed pairs. Genotypes from 14 breeds sampled from Europe and Asia were used to obtain estimates of pair-wise population divergence measured as F ST. The divergence separating the Florida and Louisiana Native (F ST = 6.2%) was approximately 50% higher than the average divergence separating breeds developed within the same region of Europe (F ST = 4.2%). This strongly indicated that the two Gulf Coast Native lines are sufficiently different to be considered separate breeds. PCA using small SNP sets successfully distinguished between the Florida and Louisiana Native animals, suggesting that allele frequency differences have accumulated across the genome. This is consistent with a population history involving geographic separation and genetic drift. Suggestive evidence was detected for divergence at the poll locus on sheep chromosome 10; however drift at neutral markers has been the largest contributor to the genetic separation observed. These results document the emergence of populations that can be considered separate breeds, with practical consequences for bio-conservation priorities, animal registration and the establishment of separate breed societies.

Adaptation shapes spike train correlations: theory and application to tinnitus

weakened adaptation reported in neurons from animals with behavioral evidence of tinnitus both increases firing rate and spike train correlations. Our work shows that the biophysical correlates of tinnitus are consistent with the spike train correlates of tinnitus, and that this link is not obvious from past theories linking firing rates and spike train correlations because of the presence of SFA. These results reveal that the biophysical mechanisms giving rise to hyperactivity in tinnitus and other neural pathologies have specific implications for the pairwise statistics and population coding of the affected neurons.

Genetic Variation in Wild and Cultured Populations of the Sea Squirt Halocynthia roretzi Inferred from Microsatellite DNA Analysis

We used nine microsatellite DNA markers to estimate genetic variation among wild and cultured populations of the sea squirt Halocynthia roretzi. The loci were polymorphic, with 6-32 alleles, and allelic richness ranged from 6.0 to 26.1 in each population. The wild and the cultured populations had similar mean heterozygosities (<TEX>$H_O$</TEX> and <TEX>$H_E$</TEX>), allele numbers, and allelic richness. One cultured population with softness syndrome had a lower mean in the observed heterozygosity (<TEX>$H_O$</TEX> = 0.57) and higher mean inbreeding coefficient (<TEX>$F_{IS}$</TEX> = 0.261) than any other populations. This suggests that the loss of genetic variation in the diseased population might be due to increased inbreeding. A neighbor-joining tree and pairwise population estimates of <TEX>$F_{ST}$</TEX> showed moderate genetic differentiation between the wild and the cultured populations. Additionally, the softness syndrome population was genetically divergent from wild populations, but it was genetically close to the cultured populations.

Pronounced genetic structure and low genetic diversity in European red-billed chough (Pyrrhocorax pyrrhocorax) populations

The red-billed chough (Pyrrhocorax pyrrhocorax) is of conservation concern in the British Isles and continental Europe, with historically declining populations and a highly fragmented distribution. We quantified the distribution of genetic variation within and among European populations to identify isolated populations that may need to be managed as demographically independent units, and assess whether individual populations are denuded of genetic diversity and so may show reduced viability. We genotyped 326 choughs from ten wild populations and 22 from one captive population at 16 nuclear microsatellite loci, and sequenced 34 individuals across three mitochondrial regions to quantify genetic structure, diversity and phylogeography. Microsatellite diversity was low (often <4 alleles per locus), but pairwise population differentiation was high (often D est > 0.1), with a signature of isolation-by-distance. Bayesian-inferred a posteriori genetic clusters coincided with a priori populations, supporting strong genetic structure. Microsatellites also allowed us to identify the probable origin of the captive choughs and one recently founded wild population. Mitochondrial DNA sequence diversity was low (π = 0.00103). Phylogeographic structure was consequently poorly resolved, but indicated that sampled continental-European populations are ancestral to British Isles populations, which comprised a single clade. Our data suggest that British Isles chough populations are relatively isolated with infrequent gene flow and relatively genetically depauperate, potentially requiring genetic management. These findings should be integrated into conservation management policy to ensure long-term viability of chough populations.

Population genetic structure and phylogeography of the ascidian, Halocynthia roretzi, along the coasts of Korea and Japan, inferred from mitochondrial DNA sequence analysis

Abstract The genetic structure and phylogeography of the solitary ascidian Halocynthia roretzi (sea pineapple) in Korea and Japan were studied based on mitochondrial DNA cytochrome oxidase subunit I ( COI ) gene nucleotide variations. The COI haplotypes found were divided into two clades in a statistical parsimony network. From the extensive haplotype diversity, limited nucleotide diversity, neutrality test, and mismatch distribution analysis, the species appears to have undergone rapid population expansion from an ancestral population with a small effective population size due to glacial bottlenecking in the late Pleistocene. Pairwise population F ST estimates inferred distinct genetic differentiation between the population from the northernmost east coast of Korea and the remaining pooled coastal water populations, while other populations generally exhibited low levels of differentiation. Analysis of molecular variance showed a weak but significant differentiation among three geographic population groups. Overall, the moderate population structure of H. roretzi may have been influenced by high levels of gene flow caused by hydrographic conditions such as sea currents.

Microsatellite genetic diversity and differentiation of native and introduced grass carp populations in three continents

Grass carp (Ctenopharyngodon idella), a freshwater species native to China, has been introduced to about 100 countries/regions and poses both biological and environmental challenges to the receiving ecosystems. In this study, we analyzed genetic variation in grass carp from three introduced river systems (Mississippi River Basin in US, Danube River in Hungary, and Tone River in Japan) as well as its native ranges (Yangtze, Pearl, and Amur Rivers) in China using 21 novel microsatellite loci. The allelic richness, observed heterozygosity, and within-population gene diversity were found to be lower in the introduced populations than in the native populations, presumably due to the small founder population size of the former. Significant genetic differentiation was found between all pairwise populations from different rivers. Both principal component analysis and Bayesian clustering analysis revealed obvious genetic distinction between the native and introduced populations. Interestingly, genetic bottlenecks were detected in the Hungarian and Japanese grass carp populations, but not in the North American population, suggesting that the Mississippi River Basin grass carp has experienced rapid population expansion with potential genetic diversification during the half-century since its introduction. Consequently, the combined forces of the founder effect, introduction history, and rapid population expansion help explaining the observed patterns of genetic diversity within and among both native and introduced populations of the grass carp.

Genetic Diversity of Pinus nigra Arn. Populations in Southern Spain and Northern Morocco Revealed By Inter-Simple Sequence Repeat Profiles

Eight Pinus nigra Arn. populations from Southern Spain and Northern Morocco were examined using inter-simple sequence repeat markers to characterize the genetic variability amongst populations. Pair-wise population genetic distance ranged from 0.031 to 0.283, with a mean of 0.150 between populations. The highest inter-population average distance was between PaCU from Cuenca and YeCA from Cazorla, while the lowest distance was between TaMO from Morocco and MA Sierra Mágina populations. Analysis of molecular variance (AMOVA) and Nei’s genetic diversity analyses revealed higher genetic variation within the same population than among different populations. Genetic differentiation (Gst) was 0.233. Cuenca showed the highest Nei’s genetic diversity followed by the Moroccan region, Sierra Mágina, and Cazorla region. However, clustering of populations was not in accordance with their geographical locations. Principal component analysis showed the presence of two major groups—Group 1 contained all populations from Cuenca while Group 2 contained populations from Cazorla, Sierra Mágina and Morocco—while Bayesian analysis revealed the presence of three clusters. The low genetic diversity observed in PaCU and YeCA is probably a consequence of inappropriate management since no estimation of genetic variability was performed before the silvicultural treatments. Data indicates that the inter-simple sequence repeat (ISSR) method is sufficiently informative and powerful to assess genetic variability among populations of P. nigra.

Genetic Delineation of Local Provenance in Persoonia longifolia: Implications for Seed Sourcing for Ecological Restoration

Restoration of diverse native plant communities typically requires the collection of large amounts of seed. Thus, practitioners often struggle to find adequate supplies near project sites and need to know from how far they can collect without compromising restoration success—how far does local provenance extend? We addressed this issue by assessing genetic variation within, and differentiation among, 12 potential seed source populations of Persoonia longifolia, a key component of the jarrah forest of Western Australia. An analysis of molecular variance of 66 polymorphic amplified fragment length polymorphism markers partitioned 92% of the total genetic variation within populations and 8% among populations, indicating relatively weak but statistically significant population genetic differentiation. Ordination of these genetic data showed marked west/east and north/south gradients. Pairwise population genetic dissimilarity was correlated with both geographic distance and environmental distance derived from five climate variables. However, partial Mantel tests showed that the relationship between genetic and geographic distance was not independent of environmental distance, suggesting a non‐neutral signature in these markers. Bayesian outlier analysis identified two markers, and spatial analysis method tests identified highly significant associations between these two markers and three environmental variables. Frequency differences at these markers across populations suggested the possibility of climatically adapted provenances. The global significance value from analyses of similarities for these two markers correlated to a general provenance distance of 47 km, in contrast to a threshold of 60 km for the complete dataset. Guidelines for seed sourcing that consider these population genetic data should lead to more effective ecological restoration with this species.

A forensic DNA profiling system for Northern European brown bears (Ursus arctos)

A set of 13 dinucleotide STR loci (G1A, G10B, G1D, G10L, MU05, MU09, MU10, MU15, MU23, MU26, MU50, MU51, MU59) were selected as candidate markers for a DNA forensic profiling system for Northern European brown bear (Ursus arctos). We present results from validation of the markers with respect to their sensitivity, species specificity and performance (precision, heterozygote balance and stutter ratios). All STRs were amplified with 0.6ng template input, and there were no false bear genotypes in the cross-species amplification tests. The validation experiments showed that stutter ratios and heterozygote balance was more pronounced than in the tetranucleotide loci used in human forensics. The elevated ratios of stutter and heterozygote balance at the loci validated indicate that these dinucleotide STRs are not well suited for interpretation of individual genotypes in mixtures. Based on the results from the experimental validations we discuss the challenges related to genotyping dinucleotide STRs in single source samples. Sequence studies of common alleles showed that, in general, the size variation of alleles corresponded with the variation in number of repeats. The samples characterized by sequence analysis may serve as standard DNA samples for inter laboratory calibration. A total of 479 individuals from eight Northern European brown bear populations were analyzed in the 13 candidate STRs. Locus MU26 was excluded as a putative forensic marker after revealing large deviations from expected heterozygosity likely to be caused by null-alleles at this locus. The remaining STRs did not reveal significant deviations from Hardy–Weinberg equilibrium expectations except for loci G10B and MU10 that showed significant deviations in one population each, respectively. There were 9 pairwise locus comparisons that showed significant deviation from linkage equilibrium in one or two out of the eight populations. Substantial genetic differentiation was detected in some of the pairwise population comparisons and the average estimate of population substructure (FST) was 0.09. The average estimate of inbreeding (FIS) was 0.005. Accounting for population substructure and inbreeding the total average probability of identity in each of the eight populations was lower than 1.1×10−9 and the total average probability of sibling identity was lower than 1.3×10−4. The magnitude of these measurements indicates that if applying these twelve STRs in a DNA profiling system this would provide individual specific evidence.

Genetic diversity and population structure of Dactylorhiza hatagirea (Orchidaceae) in cold desert Ladakh region of India.

Random amplified polymorphic DNA (RAPD) analysis was used to characterize the genetic diversity and population genetic structure within and among nine natural populations ofDactylorhiza hatagirea, a critically endangered or rare terrestrial medicinal orchid in cold desert of Ladakh region. Out of the 177 bands generated from twenty random primers, 174 were polymorphic. The genetic diversity of D. hatagirea that was revealed by observed number of alleles (Na), expected number of alleles (Ne), Nei’s diversity index (H), Shannon’s diversity index (I), amplificated loci, polymorphic loci and the percentage of polymorphic loci (PPL). Pair-wise population genetic distances ranged from 0.05 to 0.23. Analysis of molecular variance (AMOVA) revealed that 57% RAPD of variability was partitioned among population. The Principal coordinate’s analysis (PCoA) and UPGMA supported the grouping of all 96 accessions of nine populations into four cluster groups. However, a moderate genetic differentiation among population was detected based on different measures (Nei’s genetic diversity analysis: Gst = 0.2538; AMOVA analysis: Fst = 0.254). Key words: Conservation biology, endangered orchid, genetic diversity, genetic structure.

Development and Characterization of 18 Novel EST-SSRs from the Western Flower Thrips, Frankliniella occidentalis (Pergande)

The western flower thrips, Frankliniella occidentalis (Pergande), is an invasive species and the most economically important pest within the insect order Thysanoptera. For a better understanding of the genetic makeup and migration patterns of F. occidentalis throughout the world, we characterized 18 novel polymorphic EST-derived microsatellites. The mutational mechanism of these EST-SSRs was also investigated to facilitate the selection of appropriate combinations of markers for population genetic studies. Genetic diversity of these novel markers was assessed in 96 individuals from three populations in China (Harbin, Dali, and Guiyang). The results showed that all these 18 loci were highly polymorphic; the number of alleles ranged from 2 to 15, with an average of 5.50 alleles per locus. The observed (HO) and expected (HE) heterozygosities ranged from 0.072 to 0.707 and 0.089 to 0.851, respectively. Furthermore, only two locus/population combinations (WFT144 in Dali and WFT50 in Guiyang) significantly deviated from Hardy–Weinberg equilibrium (HWE). Pairwise FST analysis showed a low but significant differentiation (0.026 < FST < 0.032) among all three pairwise population comparisons. Sequence analysis of alleles per locus revealed a complex mutational pattern of these EST-SSRs. Thus, these EST-SSRs are useful markers but greater attention should be paid to the mutational characteristics of these microsatellites when they are used in population genetic studies.

Genetic characterization and gene flow in different geographical-distance neighbouring natural populations of wild soybean (Glycine soja Sieb. &amp; Zucc.) and implications for protection from GM soybeans

To design appropriate strategies for ex situ and in situ conservation of wild soybean and safeguard the biosafety of the wild soybean gene pool when genetically modified soybeans are grown, it is important to understand its genetic characteristics, and to quantify gene flow and kinship within and between neighboring populations. We analyzed 9 pairs of neighboring populations of wild soybeans using 20 pairs of nuclear SSR markers. Results showed that Chinese wild soybean natural populations had outcrossing rates of 0–3.5% and that most populations contained many kinship families. The kinship families could be attributed to the accumulation of outcrossed offspring within populations during the history of population colonization. Wild soybean is very sensitive to environmental selection, which results in genetic differentiation of populations, and the emergence of specific alleles. We used an index τ to explain why genetic differences would exist between the pairwise populations; the interpopulation genetic differentiation chiefly consisted in the differences of allele frequencies over the genome. We found long-distance dispersal (1.5 km) of wild soybean seeds in a land ecosystem. There was close correlation between genetic and geographical distance among natural populations of Chinese wild soybean. Within a distance of 50 km, there was greater gene flow when the distance between populations was shorter. These findings have implications for ex situ and in situ conservation in an ecogeographical region, and also for protection of the gene pool from contamination by GM soybeans in wild soybean species.

Genetic Structure in Wild Populations of Ayu Plecoglossus altivelis in Korea and Japan

We investigated the genetic structure of Korean and Japanese ayu Plecoglossus altivelis populations by examining 669 individuals from 14 populations using three microsatellite loci. Genetic variation did not differ significantly among the populations examined in terms of allelic number and heterozygosity. Korean populations were genetically close to each other, implying that persistent gene flow has occurred in these populations. This suggests that eastern populations in Korea form a single large population and all of the Korean populations are distinct from the Japanese populations. Pairwise population <TEX>$F_{ST}$</TEX> estimates, principal component analyses, and a neighbor-joining tree showed that genetic separation between the southern and pooled eastern coast populations was probably influenced by restricted gene flow. Hierarchical analysis of molecular variance (AMOVA) revealed a weak but significant genetic structure among three ayu groups (eastern and southern coasts of Korea and the Japan coast), and no genetic variation within groups. The estimated genetic population structure and potential applications of microsatellite markers may aid in the proper management of ayu populations.

Oceanic dispersal in a sedentary reef shark (Triaenodon obesus): genetic evidence for extensive connectivity without a pelagic larval stage

AbstractAim Most reef fishes are site‐attached, but can maintain a broad distribution through their highly dispersive larval stage. The whitetip reef shark (Triaenodon obesus) is site‐attached, yet maintains the largest Indo‐Pacific distribution of any reef shark while lacking the larval stage of bony (teleost) fishes. Here we use mitochondrial DNA (mtDNA) sequence data to evaluate the enigma of the sedentary reef shark that maintains a distribution across two‐thirds of the planet.Location Tropical Pacific and Indian Oceans.Methods We analysed 1025 base pairs of the mtDNA control region in 310 individuals from 25 locations across the Indian and Pacific Oceans. Phylogeographic and population genetic analyses were used to reveal the dispersal and recent evolutionary history of the species.Results We resolved 15 mtDNA control region haplotypes, but two comprised 87% of the specimens and were detected at nearly every location. Similar to other sharks, genetic diversity was low (h = 0.550 ± 0.0254 and π = 0.00213 ± 0.00131). Spatial analyses of genetic variation demonstrated strong isolation across the Indo‐Pacific Barrier and between western and central Pacific locations. Pairwise ΦST comparisons indicated high connectivity among archipelagos of the central Pacific but isolation across short distances of contiguous habitat (Great Barrier Reef) and intermittent habitat (Hawaiian Archipelago). In the eastern Pacific only a single haplotype (the most common one in the central Pacific) was observed, indicating recent dispersal (or colonization) across the East Pacific Barrier.Main conclusions The shallow haplotype network indicates recent expansion of modern populations within the last half million years from a common ancestor. Based on the distribution of mtDNA diversity, this began with an Indo‐West Pacific centre of origin, with subsequent dispersal to the Central Pacific and East Pacific. Genetic differences between Indian and Pacific Ocean populations are consistent with Pleistocene closures of the Indo‐Pacific Barrier associated with glacial cycles. Pairwise population comparisons reveal weak but significant isolation by distance, and notably do not indicate the high coastal connectivity observed in other shark species. The finding of population structure among semi‐contiguous habitats, but population connectivity among archipelagos, may indicate a previously unsuspected oceanic dispersal behaviour in whitetip reef sharks.

Shallow Population Genetic Structures of Thread-sail Filefish (Stephanolepis cirrhifer) Populations from Korean Coastal Waters.

Genetic diversities, population genetic structures and demographic histories of the thread-sail filefish Stephanolepis cirrhifer were investigated by nucleotide sequencing of 336 base pairs of the mitochondrial DNA (mtDNA) control region in 111 individuals collected from six populations in Korean coastal waters. A total of 70 haplotypes were defined by 58 variable nucleotide sites. The neighbor-joining tree of the 70 haplotypes was shallow and did not provide evidence of geographical associations. Expansion of S. cirrhifer populations began approximate 51,000 to 102,000 years before present, correlating with the period of sea level rise since the late Pleistocene glacial maximum. High levels of haplotype diversities (0.974±0.029 to 1.000±0.076) and nucleotide diversities (0.014 to 0.019), and low levels of genetic differentiation among populations inferred from pairwise population FST values (−0.007 to 0.107), support an expansion of the S. cirrhifer population. Hierarchical analysis of molecular variance (AMOVA) revealed weak but significant genetic structures among three groups (FCT = 0.028, p<0.05), and no genetic variation within groups (0.53%; FSC = 0.005, p = 0.23). These results may help establish appropriate fishery management strategies for stocks of S. cirrhifer and related species.

Genetic Diversity and Population Structure of Ethiopian Chickpea (Cicer arietinum L.) Germplasm Accessions from Different Geographical Origins as Revealed by Microsatellite Markers

Genetic diversity and population structure were studied in 155 chickpea (Cicer arietinum L.) entries using 33 polymorphic microsatellite (SSR) markers. Molecular analysis of variance showed variations of 73% within and 27% among populations. Introduced genotypes exhibited highest polymorphism (70.27%) than the landraces (36–57%). Collections from Shewa, Harerge, W. Gojam and S. Gonder regions also showed higher polymorphism (50–57%) than the rest of the local accessions (36–45%). Analyses of pairwise population Nei genetic distance and PhiPT coefficients, expected heterozygosity (He) and unbiased expected heterozygosity (UHe), Shannon’s information index (I) and percent polymorphism (% P) showed existence of high genetic variation between geographical regions. Accessions from adjoining geographical regions mostly showed more genetic similarities than those from origins far isolated apart. This could be associated with the ease and likelihood of inter-regional gene flow and seed movement particularly during times of drought. The 155 entries were grouped into five clusters following analysis of population structure. The first cluster (C1) constituted accessions from Arsi; the second (C2) from Gojam and Gonder; the third (C3) from Harerge and E. and N. Shewa; and the fourth (C4) from W. Shewa, Tigray, and Wello regions. The fifth cluster (C5) was entirely consisted of improved genotypes. Improved genotypes of both Kabuli and Desi types distinctly fell into cluster five (C5) regardless of their difference in seed types. The result has firmly established that introduction of genetic materials from exotic sources has broadened the genetic base of the national chickpea breeding program. Further implications of the findings as regards to chickpea germplasm management and its utilization in breeding program are also discussed.

Matching genetics with oceanography: directional gene flow in a Mediterranean fish species

Genetic connectivity and geographic fragmentation are two opposing mechanisms determining the population structure of species. While the first homogenizes the genetic background across populations the second one allows their differentiation. Therefore, knowledge of processes affecting dispersal of marine organisms is crucial to understand their genetic distribution patterns and for the effective management of their populations. In this study, we use genetic analyses of eleven microsatellites in combination with oceanographic satellite and dispersal simulation data to determine distribution patterns for Serranus cabrilla, a ubiquitous demersal broadcast spawner, in the Mediterranean Sea. Pairwise population F(ST) values ranged between -0.003 and 0.135. Two genetically distinct clusters were identified, with a clear division located between the oceanographic discontinuities at the Ibiza Channel (IC) and the Almeria-Oran Front (AOF), revealing an admixed population in between. The Balearic Front (BF) also appeared to dictate population structure. Directional gene flow on the Spanish coast was observed as S.cabrilla dispersed from west to east over the AOF, from north to south on the IC and from south of the IC towards the Balearic Islands. Correlations between genetic and oceanographic data were highly significant. Seasonal changes in current patterns and the relationship between ocean circulation patterns and spawning season may also play an important role in population structure around oceanographic fronts.

The structure of cross-cultural musical diversity

Human cultural traits, such as languages, musics, rituals and material objects, vary widely across cultures. However, the majority of comparative analyses of human cultural diversity focus on between-culture variation without consideration for within-culture variation. In contrast, biological approaches to genetic diversity, such as the analysis of molecular variance (AMOVA) framework, partition genetic diversity into both within- and between-population components. We attempt here for the first time to quantify both components of cultural diversity by applying the AMOVA model to music. By employing this approach with 421 traditional songs from 16 Austronesian-speaking populations, we show that the vast majority of musical variability is due to differences within populations rather than differences between. This demonstrates a striking parallel to the structure of genetic diversity in humans. A neighbour-net analysis of pairwise population musical divergence shows a large amount of reticulation, indicating the pervasive occurrence of borrowing and/or convergent evolution of musical features across populations.