Genetic Origins and Sex-Biased Admixture of the Huis

The Hui people are unique among Chinese ethnic minorities in that they speak the same language as Han Chinese (HAN) but practice Islam. However, as the second-largest minority group in China numbering well over 10 million, the Huis are under-represented in both global and regional genomic studies. Here, we present the first whole-genome sequencing effort of 234 Hui individuals (NXH) aged over 60 who have been living in Ningxia, where the Huis are mostly concentrated. NXH are genetically more similar to East Asian than to any other global populations. In particular, the genetic differentiation between NXH and HAN (F_ST = 0.0015) is only slightly larger than that between northern and southern HAN (F_ST = 0.0010), largely attributed to the western ancestry in NXH (∼10%). Highly differentiated functional variants between NXH and HAN were identified in genes associated with skin pigmentation (e.g., SLC24A5), facial morphology (e.g., EDAR), and lipid metabolism (e.g., ABCG8). The Huis are also distinct from other Muslim groups such as the Uyghurs (F_ST = 0.0187), especially, NXH derived much less western ancestry (∼10%) compared with the Uyghurs (∼50%). Modeling admixture history indicated that NXH experienced an episode of two-wave admixture. An ancient admixture occurred ∼1,025 years ago, reflecting the intensive west–east contacts during the late Tang Dynasty, and the Five Dynasties and Ten Kingdoms period. A recent admixture occurred ∼500 years ago, corresponding to the Ming Dynasty. Notably, we identified considerable sex-biased admixture, that is, excess of western males and eastern females contributing to the NXH gene pool. The origins and the genomic diversity of the Hui people imply the complex history of contacts between western and eastern Eurasians.     

Intra- and interpopulation variability of southwestern Kamchatka sockeye salmon Oncorhynchus nerka inferred from the data on single nucleotide polymorphism

The variability of 45 single nucleotide polymorphism (SNP) loci was studied in nine samples of the sockeye salmon Oncorhynchus nerka from the rivers of southwestern Kamchatka. The Wahlund effect, gametic disequilibrium at some loci, and a decrease in interpopulation genetic diversity indices observed in samples from the Bolshaya River outlet can be attributed to the samples’ heterogeneity. Partitioning of the mixed samples using some biological characteristics of the individuals led to a noticeable decrease in the frequency of these phenomena. It was demonstrated that the allelic diversity between the populations within the river accounted for the larger part of genetic variation, as compared to the differentiation between the basins. The SNP loci responsible for intra- and interpopulation differentiation of sockeye salmon from the rivers of southwestern Kamchatka were identified. Some recommendations for field population genetic studies of Asian sockeye salmon were formulated.     

Distinctly different sex ratios in African and European populations of Drosophila melanogaster inferred from chromosomewide single nucleotide polymorphism data

It has been hypothesized that the ratio of X-linked to autosomal sequence diversity is influenced by unequal sex ratios in Drosophila melanogaster populations. We conducted a genome scan of single nucleotide polymorphism (SNP) of 378 autosomal loci in a derived European population and of a subset of 53 loci in an ancestral African population. On the basis of these data and our already available X-linked data, we used a coalescent-based maximum-likelihood method to estimate sex ratios and demographic histories simultaneously for both populations. We confirm our previous findings that the African population experienced a population size expansion while the European population suffered a population size bottleneck. Our analysis also indicates that the female population size in Africa is larger than or equal to the male population size. In contrast, the European population shows a huge excess of males. This unequal sex ratio and the bottleneck alone, however, cannot account for the overly strong decrease of X-linked diversity in the European population (compared to the reduction on the autosome). The patterns of the frequency spectrum and the levels of linkage disequilibrium observed in Europe suggest that, in addition, positive selection must have acted in the derived population.     

EST-derived single nucleotide polymorphism markers for assembling genetic and physical maps of the barley genome

In a panel of seven genotypes, 437 expressed sequence tag (EST)-derived DNA fragments were sequenced. Single nucleotide polymorphisms (SNPs) that were polymorphic between the parents of three mapping populations were mapped by heteroduplex analysis and a genome-wide consensus map comprising 216 EST-derived SNPs and 4 InDel (insertion/deletion) markers was constructed. The average frequency of SNPs amounted to 1/130 bp and 1/107.8 bp for a set of randomly selected and a set of mapped ESTs, respectively. The calculated nucleotide diversities (π) ranged from 0 to 40.0 × 10⁻³ (average 3.1 × 10⁻³) and 0.52 × 10⁻³ to 39.51 × 10^–3 (average 4.37 × 10⁻³) for random and mapped ESTs, respectively. The polymorphism information content value for mapped SNPs ranged from 0.24 to 0.50 with an average of 0.34. As expected, combination of SNPs present in an amplicon (haplotype) exhibited a higher information content ranging from 0.24 to 0.85 with an average of 0.50. Cleaved amplified polymorphic sequence assays (including InDels) were designed for a total of 87 (39.5%) SNP markers. The high abundance of SNPs in the barley genome provides avenues for the systematic development of saturated genetic maps and their integration with physical maps. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users. Both R. Kota and R.K. Varshney contributed equally to this work.     

Reconstructing the genetic admixture history of Tai-Kadai and Sinitic people: Insights from genome-wide SNP data from South China

South China (SC) was a region with mixed rice–millet farming during the Middle Neolithic period and was also suggested to be the homeland of Tai-Kadai (TK)-speaking people. However, the formations of inland TK-speaking people and southwestern Hans are far from clear due to very few studies on this subject. Here, we reveal the spatiotemporally demographic history of SC by analyzing newly-generated genome-wide SNP data of 115 modern southwestern individuals and find that inland TK-speaking Dongs and Bouyeis have a close genomic affinity to coastal TK/Austronesian (AN)-speaking people and Neolithic Yangtze River basin (YZRB) farmers, while southwestern Hans and TK-speaking Gelaos possess a close genomic affinity to Neolithic Yellow River basin (YRB) farmers. Genetic differentiations are identified among TK people from SC and Southeast Asia, and between northern and southern inland Chinese TK people, in which the identified shared genetic ancestry between TK and AN people highlights a common origin of AN/TK groups. Conclusively, our findings indicate that millet farmers deriving from the YRB and rice farmers deriving from the YZRB substantially contribute to the present-day inland TK speakers and southwestern Hans via a two-way admixture scenario of bi-directional gene-flow events, which facilitates the formation of a modern two-way genetic admixture profile.     

Phylogeography and adaptation genetics of stickleback from the Haida Gwaii archipelago revealed using genome‐wide single nucleotide polymorphism genotyping

Threespine stickleback populations are model systems for studying adaptive evolution and the underlying genetics. In lakes on the Haida Gwaii archipelago (off western Canada), stickleback have undergone a remarkable local radiation and show phenotypic diversity matching that seen throughout the species distribution. To provide a historical context for this radiation, we surveyed genetic variation at >1000 single nucleotide polymorphism (SNP) loci in stickleback from over 100 populations. SNPs included markers evenly distributed throughout genome and candidate SNPs tagging adaptive genomic regions. Based on evenly distributed SNPs, the phylogeographic pattern differs substantially from the disjunct pattern previously observed between two highly divergent mtDNA lineages. The SNP tree instead shows extensive within watershed population clustering and different watersheds separated by short branches deep in the tree. These data are consistent with separate colonizations of most watersheds, despite underlying genetic connections between some independent drainages. This supports previous suppositions that morphological diversity observed between watersheds has been shaped independently, with populations exhibiting complete loss of lateral plates and giant size each occurring in several distinct clades. Throughout the archipelago, we see repeated selection of SNPs tagging candidate freshwater adaptive variants at several genomic regions differentiated between marine–freshwater populations on a global scale (e.g. EDA, Na/K ATPase). In estuarine sites, both marine and freshwater allelic variants were commonly detected. We also found typically marine alleles present in a few freshwater lakes, especially those with completely plated morphology. These results provide a general model for postglacial colonization of freshwater habitat by sticklebacks and illustrate the tremendous potential of genome‐wide SNP data sets hold for resolving patterns and processes underlying recent adaptive divergences.     

Detecting signatures of inter‐regional and inter‐specific hybridization among the Chinese rhesus macaque specific pathogen‐free (SPF) population using single nucleotide polymorphic (SNP) markers

Background While rates of gene flow between rhesus and longtail macaque populations near their hybrid zone in Indochina have been quantified elsewhere, this study demonstrates that the inter‐specific introgression is not limited to the Indochinese hybrid zone but is more geographically widespread. Methods Twelve rhesus and longtail macaque populations were analyzed using single nucleotide polymorphic (SNP) loci. Results There is evidence for inter‐specific admixture between Chinese rhesus and mainland longtails, with implications for genetic diversity both in the Chinese super‐SPF population at the California National Primate Research Center and in other primate facilities. Eastern Chinese rhesus appeared more highly derived than western Chinese rhesus, and allele sharing between longtails and Chinese rhesus was not random with regard to geographic distance, but no significant nuclear genetic differences between eastern and western Chinese rhesus were detected among the 245 genic SNPs assayed. Conclusion The implications of this inter‐specific admixture for the use of Chinese rhesus and mainland longtail in biomedical research should be considered.     

Whole-genome single nucleotide polymorphism and mating compatibility studies reveal the presence of distinct species in sub-Saharan Africa Bemisia tabaci whiteflies

In sub-Saharan Africa cassava growing areas, two members of the Bemisia tabaci species complex termed sub-Saharan Africa 1 (SSA1) and SSA2 have been reported as the prevalent whiteflies associated with the spread of viruses that cause cassava mosaic disease (CMD) and cassava brown streak disease (CBSD) pandemics. At the peak of CMD pandemic in the late 1990s, SSA2 was the prevalent whitefly, although its numbers have diminished over the last two decades with the resurgence of SSA1 whiteflies. Three SSA1 subgroups (SG1 to SG3) are the predominant whiteflies in East Africa and vary in distribution and biological properties. Mating compatibility between SSA1 subgroups and SSA2 whiteflies was reported as the possible driver for the resurgence of SSA1 whiteflies. In this study, a combination of both phylogenomic methods and reciprocal crossing experiments were applied to determine species status of SSA1 subgroups and SSA2 whitefly populations. Phylogenomic analyses conducted with 26 548 205 bp whole genome single nucleotide polymorphisms (SNPs) and the full mitogenomes clustered SSA1 subgroups together and separate from SSA2 species. Mating incompatibility between SSA1 subgroups and SSA2 further demonstrated their distinctiveness from each other. Phylogenomic analyses conducted with SNPs and mitogenomes also revealed different genetic relationships among SSA1 subgroups. The former clustered SSA1-SG1 and SSA1-SG2 together but separate from SSA1-SG3, while the latter clustered SSA1-SG2 and SSA1-SG3 together but separate from SSA1-SG1. Mating compatibility was observed between SSA1-SG1 and SSA1-SG2, while incompatibility occurred between SSA1-SG1 and SSA1-SG3, and SSA1-SG2 and SSA1-SG3. Mating results among SSA1 subgroups were coherent with phylogenomics results based on SNPs but not the full mitogenomes. Furthermore, this study revealed that the secondary endosymbiont—Wolbachia—did not mediate reproductive success in the crossing assays carried out. Overall, using genome wide SNPs together with reciprocal crossings assays, this study established accurate genetic relationships among cassava-colonizing populations, illustrating that SSA1 and SSA2 are distinct species while at least two species occur within SSA1 species.     

Relationship between single nucleotide polymorphism of interleukin-18 and susceptibility to pulmonary tuberculosis in the Chinese Han population

Interleukin-18 (IL-18) is a multi-functional cytokine capable of inducing either Th1 or Th2 polarization depending on the immunologic milieu. IL-18 may influence the host response to Mycobacterium tuberculosis (M.tb) infection. To investigate the relationship between single nucleotide polymorphisms of the IL-18 and susceptibility to pulmonary tuberculosis in the Chinese Han population, the IL-18 gene was sequenced to detect polymorphisms and to examine the genotype frequencies in 300 patients and 702 healthy controls. DNA sequencing revealed three IL-18 variants: rs1946518, rs5744247, and rs549908. It also revealed that allele A of rs1946518 confers a 1.47-fold increased risk of developing tuberculosis (TB) (P = 0.0001, OR [95%CI] = 1.47 [1.21-1.78]), and that the C allele of rs5744247 confers a 0.77-fold decreased risk of disease (P = 0.01, R [95%CI] = 0.77 [0.632-0.937]). The genotypes rs1946518, rs5744247 and rs549908 were found to be significantly associated with TB. Estimation of the frequencies of haplotypes revealed a potential risk haplotype AGA (P = 0.01, OR [95%CI] = 1.41 [1.15-1.72]) and a protective haplotype CCA (P = 0.01, OR [95%CI] = 0.70 [0.57-0.85]) for TB. The present findings suggest that polymorphisms in the IL-18 gene may affect susceptibility to TB and increase the risk of developing the disease in the Chinese Han population.     

Linked selection,ancient polymorphism,and ecological adaptation shape the genomic landscape of divergence in Quercus dentata

Multiple evolutionary forces contribute to heterogeneous genomic landscapes; however, disentangling their relative contributions is challenging. We sampled nine populations across the distribution of Quercus dentata, a dominant forest tree in East Asia, and used whole-genome sequencing data to investigate mechanisms underlying divergence. We identified two genetic groups (north and south) that diverged ~1.84 million years ago, consistent with the uplift of the Qinling Mountains during the Pleistocene. The north group experienced a bottleneck during the middle–late Pleistocene and expanded from multiple refugia. The south group experienced a more severe bottleneck and showed high population differentiation, probably due to long-term isolation and habitat fragmentation. We detected genomic islands with elevated relative differentiation (F_ST) scattered across the genome. Among these, 65.9% showed reduced absolute divergence (d_XY) consistent with linked selection, while the remaining (34.1%) showed elevated d_XY suggestive of divergent sorting of ancient polymorphisms. The recombination rate in genomic islands was lower than background, suggesting the importance of genome structure in shaping the genomic landscape. We detected 108 single nucleotide polymorphisms significantly associated with environmental factors, 12 of which clustered in a region of ~500 kb. This region showed multiple signals of positive selection in the north group, including the enrichment of XP-extended haplotype homozygosity scores, an elevated population branch statistic, and an excess of high-frequency derived alleles. In addition, we found that linkage disequilibrium was low and derived haplotypes declined rapidly in this region, indicating selection on standing variation. Our results clarify the evolutionary processes driving genomic divergence in Q. dentata.     

Phylogeography and Pleistocene refugia of the adder (Vipera berus) as inferred from mitochondrial DNA sequence data

In order to contribute to the debate about southern glacial refugia used by temperate species and more northern refugia used by boreal or cold-temperate species, we examined the phylogeography of a widespread snake species (Vipera berus) inhabiting Europe up to the Arctic Circle. The analysis of the mitochondrial DNA (mtDNA) sequence variation in 1043 bp of the cytochrome b gene and in 918 bp of the noncoding control region was performed with phylogenetic approaches. Our results suggest that both the duplicated control region and cytochrome b evolve at a similar rate in this species. Phylogenetic analysis showed that V. berus is divided into three major mitochondrial lineages, probably resulting from an Italian, a Balkan and a Northern (from France to Russia) refugial area in Eastern Europe, near the Carpathian Mountains. In addition, the Northern clade presents an important substructure, suggesting two sequential colonization events in Europe. First, the continent was colonized from the three main refugial areas mentioned above during the Lower-Mid Pleistocene. Second, recolonization of most of Europe most likely originated from several refugia located outside of the Mediterranean peninsulas (Carpathian region, east of the Carpathians, France and possibly Hungary) during the Mid-Late Pleistocene, while populations within the Italian and Balkan Peninsulas fluctuated only slightly in distribution range, with larger lowland populations during glacial times and with refugial mountain populations during interglacials, as in the present time. The phylogeographical structure revealed in our study suggests complex recolonization dynamics of the European continent by V. berus, characterized by latitudinal as well as altitudinal range shifts, driven by both climatic changes and competition with related species.     

Tool for genomic selection and breeding to evolutionary adaptation: Development of a 100K single nucleotide polymorphism array for the honey bee

High‐throughput high‐density genotyping arrays continue to be a fast, accurate, and cost‐effective method for genotyping thousands of polymorphisms in high numbers of individuals. Here, we have developed a new high‐density SNP genotyping array (103,270 SNPs) for honey bees, one of the most ecologically and economically important pollinators worldwide. SNPs were detected by conducting whole‐genome resequencing of 61 honey bee drones (haploid males) from throughout Europe. Selection of SNPs for the chip was done in multiple steps using several criteria. The majority of SNPs were selected based on their location within known candidate regions or genes underlying a range of honey bee traits, including hygienic behavior against pathogens, foraging, and subspecies. Additionally, markers from a GWAS of hygienic behavior against the major honey bee parasite Varroa destructor were brought over. The chip also includes SNPs associated with each of three major breeding objectives—honey yield, gentleness, and Varroa resistance. We validated the chip and make recommendations for its use by determining error rates in repeat genotypings, examining the genotyping performance of different tissues, and by testing how well different sample types represent the queen's genotype. The latter is a key test because it is highly beneficial to be able to determine the queen's genotype by nonlethal means. The array is now publicly available and we suggest it will be a useful tool in genomic selection and honey bee breeding, as well as for GWAS of different traits, and for population genomic, adaptation, and conservation questions.     

Genetic variation and population genetic structure of Laodelphax striatellus via genome‐wide single nucleotide polymorphisms from specific‐locus amplified fragment‐sequencing

The small brown planthopper (SBPH), Laodelphax striatellus, is one of the most destructive agricultural pests that causes serious economic loss in the main rice‐producing areas of China. To clarify issues such as the genetic differentiation, gene flow and population genetic structure of SBPH populations, we investigated the genetic diversity, genetic structure and phylogeography of 27 SBPH populations at 23 sampling sites from three climatic zones of China using specific‐locus amplified fragment sequencing (SLAF‐seq) for large‐scale single nucleotide polymorphism (SNP) detection. In total, 115.95 M reads, 56,355 polymorphic specific‐locus amplified fragments were developed, and 32,556 reliable single nucleotides (SNPs) were detected. The results indicated that the genotypes of many polymorphism sites had low heterozygosity in every population. Overall, the pairwise F_ST values between the populations varied from 0.056 to 0.092; it suggested the lack of strong differentiation among three climatic zone groups, respectively. It also suggested a strong level of gene flow (Nm) among populations in different climate zones, ranging from 28.864 to 35.197. Phylogenetic analyses, principal component analysis (PCA), Bayesian clustering method and AMOVA revealed that there was no evidence genetic clustering in three main bioclimatic zones. Neutrality testing provided strong evidence for a recent rapid expansion without any recent genetic bottleneck in these populations. Accordingly, the results of the present study should be beneficial for SBPH management and provide insight into the genetics of SBPH.     

Association between Glu298Asp/677C-T single nucleotide polymorphism in the eNOS/MTHRF gene and blood stasis syndrome of ischemic stroke

Blood stasis syndrome of ischemic stroke (BSS-IS) is a common clinical phenotype that may be affected by certain mutagenic environmental factors or chemotherapeutic drugs; however, the role of susceptibility genes remains unclear. Previous studies have shown that ischemic stroke (IS) was closely associated with the Glu298Asp polymorphism in the eNOS gene and the 677C-T (Ala→Val) polymorphism in methylenetetrahydrofolate reductase (MTHRF) gene. Therefore, these two single nucleotide polymorphisms (SNPs) were selected to detect their associations with BSS-IS in this study. A SNP chip was employed to screen the SNP variation between both groups, and the results were verified using denaturing high-performance liquid chromatography (DHPLC) and restriction fragment length polymorphism (RFLP). The results confirmed that the TT genotype of Glu298Asp in the eNOS gene may be one of the risk factors associated with BSS-IS, while the genotype of 677C-T (Ala→Val) in the MTHRF gene may not be relevant to BSS-IS.     

突尼斯两种麦瘿蚊的遗传变异和相关性分析

小麦黑森瘿蚊Mayetiola destructor Say和大麦茎干瘿蚊M. hordei Kieffer是在突尼斯每年都可导致谷物重大损失的两个植食性姊妹种。通常认为为害小麦的瘿蚊是小麦黑森瘿蚊,但是不同谷类物种（小麦或大麦）与麦瘿蚊的这两个种（ destructor 或 hordei ）之间寄主关系并不很严格。提出有效的害虫管理方案首先要求对瘿蚊基因型进行精确分析。本研究应用随机扩增多态DNA（RAPD）技术,结合交配分析和线粒体DNA分型技术,对位于突尼斯北部的一个为害小麦的麦瘿蚊种群的遗传变异程度和分类关联性进行了评估。基于RAPD结果的系统发育分析表明,所研究的种群具有较大的遗传变异范围,这可能源于被分析的小麦样品有被两种瘿蚊共同侵害的复杂背景。虽然交配分析表明有少数不能成功产卵(2/14),但是基于细胞色素b基因限制性酶切分析显示全部样品的线粒体分型均属M. destructor。本文结果进一步证明以RAPD可变性作为分类推断依据不可靠,还为突尼斯M. destructor和M. hordei属于异域分布的观点提供了补充证据。     

Genomic footprints of an old affair: Single nucleotide polymorphism data reveal historical hybridization and the subsequent evolution of reproductive barriers in two recently diverged grasshoppers with partly overlapping distributions

Secondary contact in close relatives can result in hybridization and the admixture of previously isolated gene pools. However, after an initial period of hybridization, reproductive isolation can evolve through different processes and lead to the interruption of gene flow and the completion of the speciation process. Omocestus minutissimus and O. uhagonii are two closely related grasshoppers with partially overlapping distributions in the Central System mountains of the Iberian Peninsula. To analyse spatial patterns of historical and/or contemporary hybridization between these two taxa and understand how species boundaries are maintained in the region of secondary contact, we sampled sympatric and allopatric populations of the two species and obtained genome‐wide single nucleotide polymorphism data using a restriction site‐associated DNA sequencing approach. We used Bayesian clustering analyses to test the hypothesis of contemporary hybridization in sympatric populations and employed a suite of phylogenomic approaches and a coalescent‐based simulation framework to evaluate alternative hypothetical scenarios of interspecific gene flow. Our analyses rejected the hypothesis of contemporary hybridization but revealed past introgression in the area where the distributions of the two species overlap. Overall, these results point to a scenario of historical gene flow after secondary contact followed by the evolution of reproductive isolation that currently prevents hybridization among sympatric populations.