Genome-wide Single Nucleotide Polymorphism Analysis Research Articles

Respiratory infection- and asthma-prone, low vaccine responder children demonstrate distinct mononuclear cell DNA methylation pathways

BackgroundInfants with frequent viral and bacterial respiratory infections exhibit compromised immunity to routine immunizations. They are also more likely to develop chronic respiratory diseases in later childhood. This study investigated the feasibility of epigenetic profiling to reveal endotype-specific molecular pathways with potential for early identification and immuno-modulation. Peripheral blood mononuclear cells from respiratory infection allergy/asthma-prone (IAP) infants and non-infection allergy/asthma prone (NIAP) were retrospectively selected for genome-wide DNA methylation and single nucleotide polymorphism analysis. The IAP infants were enriched for the low vaccine responsiveness (LVR) phenotype (Fisher's exact p-value = 0.02).ResultsAn endotype signature of 813 differentially methylated regions (DMRs) comprising 238 lead CpG associations (FDR < 0.05) emerged, implicating pathways related to asthma, mucin production, antigen presentation and inflammasome activation. Allelic variation explained only a minor portion of this signature. Stimulation of mononuclear cells with monophosphoryl lipid A (MPL), a TLR agonist, partially reversed this signature at a subset of CpGs, suggesting the potential for epigenetic remodeling.ConclusionsThis proof-of-concept study establishes a foundation for precision endotyping of IAP children and highlights the potential for immune modulation strategies using adjuvants for future investigation.

Large-scale genome-wide SNP analysis reveals the rugged (and ragged) landscape of global ancestry, phylogeny, and demographic history in chicken breeds.

The worldwide chicken gene pool encompasses a remarkable, but shrinking, number of divergently selected breeds of diverse origin. This study was a large-scale genome-wide analysis of the landscape of the complex molecular architecture, genetic variability, and detailed structure among 49 populations. These populations represent a significant sample of the world's chicken breeds from Europe (Russia, Czech Republic, France, Spain, UK, etc.), Asia (China), North America (USA), and Oceania (Australia). Based on the results of breed genotyping using the Illumina 60K single nucleotide polymorphism (SNP) chip, a bioinformatic analysis was carried out. This included the calculation of heterozygosity/homozygosity statistics, inbreeding coefficients, and effective population size. It also included assessment of linkage disequilibrium and construction of phylogenetic trees. Using multidimensional scaling, principal component analysis, and ADMIXTURE-assisted global ancestry analysis, we explored the genetic structure of populations and subpopulations in each breed. An overall 49-population phylogeny analysis was also performed, and a refined evolutionary model of chicken breed formation was proposed, which included egg, meat, dual-purpose types, and ambiguous breeds. Such a large-scale survey of genetic resources in poultry farming using modern genomic methods is of great interest both from the viewpoint of a general understanding of the genetics of the domestic chicken and for the further development of genomic technologies and approaches in poultry breeding. In general, whole genome SNP genotyping of promising chicken breeds from the worldwide gene pool will promote the further development of modern genomic science as applied to poultry.

Negative Selection on a SOD1 Mutation Limits Canine Degenerative Myelopathy While Avoiding Inbreeding.

Several hundred disease-causing mutations are currently known in domestic dogs. Breeding management is therefore required to minimize their spread. Recently, genetic methods such as direct-to-consumer testing have gained popularity; however, their effects on dog populations are unclear. Here, we aimed to evaluate the influence of genetic testing on the frequency of mutations responsible for canine degenerative myelopathy and assess the changes in the genetic structure of a Pembroke Welsh corgi population from Japan. Genetic testing of 5,512 dogs for the causative mutation in superoxide dismutase 1 (SOD1) (c.118G>A (p.E40K)) uncovered a recent decrease in frequency, plummeting from 14.5% (95/657) in 2019 to 2.9% (24/820) in 2022. Weir and Cockerham population differentiation (FST) based on genome-wide single-nucleotide polymorphism (SNP) of 117 selected dogs detected the SNP with the highest FST located in the intron of SOD1 adjacent to the c.118G>A mutation, supporting a selection signature on SOD1. Further genome-wide SNP analyses revealed no obvious changes in inbreeding levels and genetic diversity between the 2019 and 2022 populations. Our study highlights that genetic testing can help inform improved mating choices in breeding programs to reduce the frequency of risk variants and avoid inbreeding. This combined strategy could decrease the genetic risk of canine degenerative myelopathy, a fatal disease, within only a few years.

Host-Shift Speciation Proceeded with Gene Flow in Algae Covering Shells.

AbstractHost shifts represent the advancement of a novel niche and often lead to speciation in symbionts. However, its mechanisms are not well understood. Here, we focused on the alga Pseudocladophora conchopheria growing on the shells of intertidal snails. Previous surveys have shown that the alga has host specificity-only attaching to the shell of Lunella correensis-but we discovered that the alga attaches to the shells of multiple sympatric snails. A genome-wide single-nucleotide polymorphism analysis (MIG-seq) was performed to determine whether host-associated speciation occurred in the algae. As a result, there was no gene flow or limited gene flow among the algae from different hosts, and some algae were genetically differentiated among hosts. In addition, the demographic estimate revealed that speciation with gene flow occurred between the algae from different hosts. Therefore, these results support the idea that host-shift speciation gradually proceeded with gene flow in the algae, providing insight into the early evolution of host shifts.

Genome-wide SNP analyses reveal a substantial gene flow and isolated-genetic structure of sea cucumber Holothuria leucospilota populations in Western Central Pacific

The black sea cucumber Holothuria leucospilota is a commercially exploited species in the Western and Central Pacific region including Australia and Vietnam, where it has been overfished. There is interest to develop captive breeding, sea ranching and restocking programs for H. leucospilota. However, the wide distribution of this species could lead to genetic structure among site, which needs to be considered for responsible genetic management for example, restocking. This is the first study to investigate population genetics in H. leucospilota, between Australia and Vietnam, using single nucleotide polymorphisms (SNPs). The results indicated that there was high gene flow (Nm > 4) and genetic similarities among three Vietnamese populations, which inferred a high panmictic stock structure in Vietnam. Less gene flow (Nm > 1) was indicated between Vietnam and Mooloolaba (Eastern Australia) populations, and there was a significant genetic differentiation among any Vietnamese populations and Mooloolaba population (FST from 0.034 to 0.043; P < 0.05). Surprisingly, there was a major genetic difference (FST > 0.2) and a restricted gene flow (Nm <1) between Darwin (Northern Australia) population and other sea cucumber populations studied (Vietnam and Eastern Australia), which revealed a genetic break/divergence in Darwin. Based on these findings, we can propose that the Vietnamese sea cucumber stocks be considered as a single genetic stock with concordant implications for genetic management. We consider whether the genetic differences in the Darwin stock arise from bottlenecking and are not reflective of adaptive differences. In summary, the results from this study provide basic genetic information on H. leucospilota populations, that would help to design restocking management plans, genetic conservation initiatives and sea ranching programs for H. leucospilota.

Taxonomic study of Scytosiphon (Phaeophyceae) from temperate coasts of Argentina.

Scytosiphon is a common intertidal genus widely distributed on temperate coasts worldwide. Recently, eight species have been delimited with molecular tools. Although S. lomentaria is the only species that predominates in the macroalgal literature of the Southwestern Atlantic Ocean (SwAO), unpublished molecular data obtained for a population study of S. lomentaria revealed hidden species diversity of Scytosiphon among the individuals collected from four localities at the SwAO. The aim of this study was to revise the identity and phylogenetic relationships of Scytosiphon from temperate coasts of the SwAO using DNA data. Thalli were collected from the Argentinean coast between 39° S and 43° S, from which cox1 and rbcL gene sequences were obtained. Phylogenies and haplotype networks were inferred and morphology of gametophytes was studied. Four species were recognized, S. lomentaria, S. promiscuus, S. shibazakiorum, and one species that belongs to a complex of species known as "Scytosiphon Atlantic complex." This complex was known to occur only in the North Atlantic, however, the results found in this study revealed that it has an extended distribution range that includes the southern hemisphere, where its populations have high genetic diversity and unique haplotypes. The morphological differences among the four species were subtle; denoting that previous Scytosiphon records from the SwAO attributed to the renowned S. lomentaria could represent different species. In addition, sex ratio and genome-wide single nucleotide polymorphisms (SNPs) analyses were done for populations of S. promiscuus presumably introduced to the SwAO, and the results indicated that they included female-dominant parthenogenetic populations, which were probably introduced from Japan.

Germline variation in RASAL2 may predict survival in patients with RAS-activated colorectal cancer.

Therapeutic agents that specifically target patients with RAS mutant colorectal cancer (CRC) are needed. We sought potential drug targets by relating genome-wide association study and survival data in patients with advanced CRC profiled for mitogen-activated protein kinase (MAPK) pathway mutations. In total, 694 patients from the clinical trials COIN and COIN-B had MAPK-activated CRCs (assigned as KRAS, NRAS, or BRAF mutant). Genome-wide single nucleotide polymorphism (SNP), gene, and gene-set analyses were performed to identify determinants of survival. For rs12028023 in RAS protein activator-like 2 (RASAL2), we studied its effect by MAPK pathway activation status (by comparing to 760 patients without MAPK-activated CRCs), MAPK gene mutation status, surface area of the primary tumor (as a marker of proliferation), and expression on RASAL2. In MAGMA genome-wide analyses, RASAL2 was the most significant gene associated with survival (p=2.0× 10-5 ). Patients carrying the minor (A) allele in the lead SNP, rs12028023 in intron 1 of RASAL2, had a median increase in survival of 167 days as compared with patients carrying the major allele. rs12028023 was predictive for survival by MAPK-activation status (pZ-test = 2.1× 10-3 ). Furthermore, rs12028023 improved survival in patients with RAS mutant (hazard ratio [HR]=0.62, 95% confidence intervals [CI]=0.5-0.8, p= 3.4× 10-5 ) but not BRAF mutant (p= 0.87) CRCs. The rs12028023 A-allele was associated with reduced surface area of the primary tumor (Beta=-0.037, standard error [SE]=0.017, p= 3.2× 10-2 ) and reduced RASAL2 expression in cultured fibroblasts (p= 1.6× 10-11 ). Our data demonstrate a prognostic role for RASAL2 in patients with MAPK-activated CRCs, with potential as a therapeutic target.

Genetic Diversity and Population Structure of Fusarium oxysporum f. sp. conglutinans Race 1 in Northern China Samples.

Fusarium oxysporum f. sp. conglutinans (FOC), the causal agent of cabbage fusarium wilt, is a serious threat to cabbage production in northern China, and most Chinese FOC isolates were identified as FOC race 1 (FOC1). To better understand the genetic diversity of FOC1 in northern China, we collected FOC isolates from five provinces in northern China and identified them as FOC1 through pathogenicity and race test. To evaluate the genome-level diversity of FOC1, we performed a genome assembly for a FOC1 isolate (FoYQ-1) collected from Yanqing, Beijing, where cabbage fusarium wilt was first reported in China. Using resequencing data of FOC1 isolates, we conducted a genome-wide SNP (single nucleotide polymorphism) analysis to investigate the genetic diversity and population structure of FOC1 isolates in northern China. Our study indicated that Chinese FOC1 can be grouped into four populations and revealed that the genetic diversity of FOC1 were closely associated with geographical locations. Our study further suggests that genetic differentiation occurred when FOC1 spread to the northwest provinces from Beijing Province in China. The FOC1 genetic diversity based on whole-genome SNPs could deepen our understanding of FOC1 variation and provide clues for the control of cabbage fusarium wilt in China.

Genome-wide survey on three local horse populations with a focus on runs of homozygosity pattern.

Purosangue Orientale Siciliano, Sanfratellano and Siciliano represent the Sicilian equine genetic resource. This study aimed to investigate the genetic diversity, population structure and the pattern of autozygosity of Sicilian horse populations using genome‐wide single‐nucleotide polymorphism (SNP) data generated with the Illumina Equine SNP70 array. The genotyping data of 17 European and Middle East populations were also included in the study. The patterns of genetic differentiation, model‐based clustering and Neighbour‐Net showed the expected positioning of Sicilian populations within the wide analysed framework and the close connections between the Purosangue Orientale Siciliano and the Arab as well as between Sanfratellano, Siciliano and Maremmano. The highest expected heterozygosity (H e) and contemporary effective population size (cNe) were reported in Siciliano (H e = 0.323, cNe = 397), and the lowest were reported in Purosangue Orientale Siciliano (H e = 0.277, cNe = 10). The analysis of the runs of homozygosity and the relative derived inbreeding revealed high internal homogeneity in Purosangue Orientale Siciliano and Arab horses, intermediate values in Maremmano and Sanfratellano and high heterogeneity in the Siciliano population. The genome‐wide SNP analysis showed the selective pressure on Purosangue Orientale Siciliano towards traits related to endurance performance. Our results underline the importance of planning adequate conservation and exploitation programmes to reduce the level of inbreeding and, therefore, the loss of genetic diversity.

Enhanced heterozygosity from male meiotic chromosome chains is superseded by hybrid female asexuality in termites

Although males are a ubiquitous feature of animals, they have been lost repeatedly in diverse lineages. The tendency for obligate asexuality to evolve is thought to be reduced in animals whose males play a critical role beyond the contribution of gametes, for example, via care of offspring or provision of nuptial gifts. To our knowledge, the evolution of obligate asexuality in such species is unknown. In some species that undergo frequent inbreeding, males are hypothesized to play a key role in maintaining genetic heterozygosity through the possession of neo-sex chromosomes, although empirical evidence for this is lacking. Because inbreeding is a key feature of the life cycle of termites, we investigated the potential role of males in promoting heterozygosity within populations through karyotyping and genome-wide single-nucleotide polymorphism analyses of the drywood termite Glyptotermes nakajimai We showed that males possess up to 15 out of 17 of their chromosomes as sex-linked (sex and neo-sex) chromosomes and that they maintain significantly higher levels of heterozygosity than do females. Furthermore, we showed that two obligately asexual lineages of this species-representing the only known all-female termite populations-arose independently via intraspecific hybridization between sexual lineages with differing diploid chromosome numbers. Importantly, these asexual females have markedly higher heterozygosity than their conspecific males and appear to have replaced the sexual lineages in some populations. Our results indicate that asexuality has enabled females to supplant a key role of males.

Genome-wide SNP discovery, linkage mapping, and analysis of QTL for morpho-physiological traits in rice during vegetative stage under drought stress.

The online version contains supplementary material available at 10.1007/s12298-021-01095-y.

Genome-wide SNP analysis of three moose subspecies at the southern range limit in the contiguous United States

Genome-wide evaluations of genetic diversity and population structure are important for informing management and conservation of trailing-edge populations. North American moose (Alces alces) are declining along portions of the southern edge of their range due to disease, species interactions, and marginal habitat, all of which may be exacerbated by climate change. We employed a genotyping by sequencing (GBS) approach in an effort to collect baseline information on the genetic variation of moose inhabiting the species’ southern range periphery in the contiguous United States. We identified 1920 single nucleotide polymorphisms (SNPs) from 155 moose representing three subspecies from five states: A. a. americana (New Hampshire), A. a. andersoni (Minnesota), and A. a. shirasi (Idaho, Montana, and Wyoming). Molecular analyses supported three geographically isolated clusters, congruent with currently recognized subspecies. Additionally, while moderately low genetic diversity was observed, there was little evidence of inbreeding. Results also indicated > 20% shared ancestry proportions between A. a. shirasi samples from northern Montana and A. a. andersoni samples from Minnesota, indicating a putative hybrid zone warranting further investigation. GBS has proven to be a simple and effective method for genome-wide SNP discovery in moose and provides robust data for informing herd management and conservation priorities. With increasing disease, predation, and climate related pressure on range edge moose populations in the United States, the use of SNP data to identify gene flow between subspecies may prove a powerful tool for moose management and recovery, particularly if hybrid moose are more able to adapt.

The Potential Role of FREM1 and Its Isoform TILRR in HIV-1 Acquisition through Mediating Inflammation.

FREM1 (Fras-related extracellular matrix 1) and its splice variant TILRR (Toll-like interleukin-1 receptor regulator) have been identified as integral components of innate immune systems. The potential involvement of FREM1 in HIV-1 (human immunodeficiency virus 1) acquisition was suggested by a genome-wide SNP (single nucleotide polymorphism) analysis of HIV-1 resistant and susceptible sex workers enrolled in the Pumwani sex worker cohort (PSWC) in Nairobi, Kenya. The studies showed that the minor allele of a FREM1 SNP rs1552896 is highly enriched in the HIV-1 resistant female sex workers. Subsequent studies showed that FREM1 mRNA is highly expressed in tissues relevant to mucosal HIV-1 infection, including cervical epithelial tissues, and TILRR is a major modulator of many genes in the NF-κB signal transduction pathway. In this article, we review the role of FREM1 and TILRR in modulating inflammatory responses and inflammation, and how their influence on inflammatory responses of cervicovaginal tissue could enhance the risk of vaginal HIV-1 acquisition.

Genome-wide SNP analysis to assess the genetic population structure and diversity of Acrocomia species.

Acrocomia (Arecaceae) is a genus widely distributed in tropical and subtropical America that has been achieving economic interest due to the great potential of oil production of some of its species. In particular A. aculeata, due to its vocation to supply oil with the same productive capacity as the oil palm (Elaeis guineenses) even in areas with water deficit. Although eight species are recognized in the genus, the taxonomic classification based on morphology and geographic distribution is still controversial. Knowledge about the genetic diversity and population structure of the species is limited, which has limited the understanding of the genetic relationships and the orientation of management, conservation, and genetic improvement activities of species of the genus. In the present study, we analyzed the genomic diversity and population structure of Acrocomia genus, including 172 samples from seven species, with a focus on A. aculeata with 117 samples covering a wide geographical area of occurrence of the species, using Single Nucleotide Polymorphism (SNP) markers originated from Genotyping By Sequencing (GBS).The genetic structure of the Acrocomia species were partially congruent with the current taxonomic classification based on morphological characters, recovering the separation of the species A. aculeata, A. totai, A. crispa and A. intumescens as distinct taxonomic groups. However, the species A. media was attributed to the cluster of A. aculeata while A. hassleri and A. glauscescens were grouped together with A. totai. The species that showed the highest and lowest genetic diversity were A. totai and A. media, respectively. When analyzed separately, the species A. aculeata showed a strong genetic structure, forming two genetic groups, the first represented mainly by genotypes from Brazil and the second by accessions from Central and North American countries. Greater genetic diversity was found in Brazil when compared to the other countries. Our results on the genetic diversity of the genus are unprecedented, as is also establishes new insights on the genomic relationships between Acrocomia species. It is also the first study to provide a more global view of the genomic diversity of A. aculeata. We also highlight the applicability of genomic data as a reference for future studies on genetic diversity, taxonomy, evolution and phylogeny of the Acrocomia genus, as well as to support strategies for the conservation, exploration and breeding of Acrocomia species and in particular A. aculeata.

Mitochondrial DNA Analysis Clarifies Taxonomic Status of the Northernmost Snow Sheep (Ovis nivicola) Population.

Currently, the intraspecific taxonomy of snow sheep (Ovis nivicola) is controversial and needs to be specified using DNA molecular genetic markers. In our previous work using whole-genome single nucleotide polymorphism (SNP) analysis, we found that the population inhabiting Kharaulakh Ridge was genetically different from the other populations of Yakut subspecies to which it was usually referred. Here, our study was aimed at the clarification of taxonomic status of Kharaulakh snow sheep using mitochondrial cytochrome b gene. A total of 87 specimens from five different geographic locations of Yakut snow sheep as well as 20 specimens of other recognized subspecies were included in this study. We identified 19 haplotypes, two of which belonged to the population from Kharaulakh Ridge. Median-joining network and Bayesian tree analyses revealed that Kharaulakh population clustered separately from all the other Yakut snow sheep. The divergence time between Kharaulakh population and Yakut snow sheep was estimated as 0.48 ± 0.19 MYA. Thus, the study of the mtDNA cytb sequences confirmed the results of genome-wide SNP analysis. Taking into account the high degree of divergence of Kharaulakh snow sheep from other groups, identified by both nuclear and mitochondrial DNA markers, we propose to classify the Kharaulakh population as a separate subspecies.

A genome-wide association study identifying the SNPs predictive of rapid joint destruction in patients with rheumatoid arthritis.

Rheumatoid arthritis (RA) is a common chronic autoimmune disease leading to joint destruction. The aim of the present study was to identify the genomic factors predictive of susceptibility to joint destruction in patients with RA by performing a genome-wide association study of genetic variants, including single nucleotide polymorphisms (SNPs). The study sample included 228 patients with a diagnosis of RA in the past 5 years. Patients were classified into rapid (total Sharp score/years of RA, ≥50) and slow (total Sharp score/years of RA, <50) joint destruction groups for analysis. The association between the genome-wide SNP analysis and joint destruction was evaluated. The following SNPs were strongly associated with rapid radiographic joint destruction: rs2295926 (P<1x10-7), belonging to the N-acetylgalactosaminyltransferase 12 (GALNT12) gene and rs11958855 (P<1x10-6), belonging to the KCNN2 gene (associated with the potassium calcium-activated channel subfamily). The identification of genetic predictors of rapid joint destruction in RA (GALNT12 and KCNN2) may provide information regarding potential therapeutic targets, and this information may be used to assist in the management RA disease progression, thereby improving the functional outcomes for patients.

Rapid Spread of Mutant Alleles in Worldwide SARS-CoV-2 Strains Revealed by Genome-Wide Single Nucleotide Polymorphism and Variation Analysis.

The novel coronavirus (SARS-CoV-2) has become a pandemic and is threatening human health globally. Here, we report nine newly evolved SARS-CoV-2 single nucleotide polymorphism (SNP) alleles those underwent a rapid increase (seven cases) or decrease (two cases) in their frequency for 30–80% in the initial four months, which are further confirmed by intrahost single nucleotide variation analysis using raw sequence data including 8,217 samples. The nine SNPs are mostly (8/9) located in the coding region and are mainly (6/9) nonsynonymous substitutions. The nine SNPs show a complete linkage in SNP pairs and belong to three different linkage groups, named LG_1 to LG_3. Analyses in population genetics show signatures of adaptive selection toward the mutants in LG_1, but no signal of selection for LG_2. Population genetic analysis results on LG_3 show geological differentiation. Analyses on geographic COVID-19 cases and published clinical data provide evidence that the mutants in LG_1 and LG_3 benefit virus replication and those in LG_1 have a positive correlation with the disease severity in COVID-19-infected patients. The mutants in LG_2 show a bias toward mildness of the disease based on available public clinical data. Our findings may be instructive for epidemiological surveys and disease control of COVID-19 in the future.

Genome-Wide SNP Analysis Reveals the Population Structure and the Conservation Status of 23 Italian Chicken Breeds.

Simple SummaryTo assess the conservation status and the genomic variability of Italian chicken breeds, we carried out genome-wide analyses using single nucleotide polymorphism (SNP) molecular markers. These results increase knowledge and can aid with the development of conservation plans for local Italian chicken breeds to safeguard their genetic variability.The genomic variability of local Italian chicken breeds, which were monitored under a conservation plan, was studied using single nucleotide polymorphisms (SNPs) to understand their genetic diversity and population structure. A total of 582 samples from 23 local breeds and four commercial stocks were genotyped using the Affymetrix 600 K Chicken SNP Array. In general, the levels of genetic diversity, investigated through different approaches, were lowest in the local chicken breeds compared to those in the commercial stocks. The level of genomic inbreeding, based on runs of homozygosity (FROH), was markedly different among the breeds and ranged from 0.121 (Valdarnese) to 0.607 (Siciliana). In all breeds, short runs of homozygosity (ROH) (<4 Mb in length) were more frequent than long segments. The patterns of genetic differentiation, model-based clustering, and neighbor networks showed that most breeds formed non-overlapping clusters and were clearly separate populations, which indicated the presence of gene flow, especially among breeds that originated from the same geographical area. Four genomic regions were identified as hotspots of autozygosity (islands) among the breeds, where the candidate genes are involved in morphological traits, such as body weight and feed conversion ratio. We conclude that the investigated breeds have conserved authentic genetic patterns, and these results can improve conservation strategies; moreover, the conservation of local breeds may play an important role in the local economy as a source of high-quality products for consumers.

Mimicry diversification in Papilio dardanus via a genomic inversion in the regulatory region of engrailed-invected.

Polymorphic Batesian mimics exhibit multiple protective morphs that each mimic a different noxious model. Here, we study the genomic transitions leading to the evolution of different mimetic wing patterns in the polymorphic Mocker Swallowtail Papilio dardanus. We generated a draft genome (231 Mb over 30 chromosomes) and re-sequenced individuals of three morphs. Genome-wide single nucleotide polymorphism (SNP) analysis revealed elevated linkage disequilibrium and divergence between morphs in the regulatory region of engrailed, a developmental gene previously implicated in the mimicry switch. The diverged region exhibits a discrete chromosomal inversion (of 40 kb) relative to the ancestral orientation that is associated with the cenea morph, but not with the bottom-recessive hippocoonides morph or with non-mimetic allopatric populations. The functional role of this inversion in the expression of the novel phenotype is currently unknown, but by preventing recombination, it allows the stable inheritance of divergent alleles enabling geographic spread and local coexistence of multiple adaptive morphs.

Admixture and natural selection shaped genomes of an Austronesian-speaking population in the Solomon Islands

People in the Solomon Islands today are considered to have derived from Asian- and Papuan-related ancestors. Papuan-related ancestors colonized Near Oceania about 47,000 years ago, and Asian-related ancestors were Austronesian (AN)-speaking population, called Lapita, who migrated from Southeast Asia about 3,500 years ago. These two ancestral populations admixed in Near Oceania before the expansion of Lapita people into Remote Oceania. To understand the impact of the admixture on the adaptation of AN-speaking Melanesians in Near Oceania, we performed the genome-wide single nucleotide polymorphism (SNP) analysis of 21 individuals from Munda, the main town of the New Georgia Islands in the western Solomon Islands. Population samples from Munda were genetically similar to other Solomon Island population samples. The analysis of genetic contribution from the two different ancestries to the Munda genome revealed significantly higher proportions of Asian- and Papuan-related ancestries in the region containing the annexin A1 (ANXA1) gene (Asian component > 82.6%) and in the human leukocyte antigen (HLA) class II region (Papuan component > 85.4%), respectively. These regions were suspected to have undergone natural selection since the time of admixture. Our results suggest that admixture had affected adaptation of AN-speaking Melanesians in the Solomon Islands.