Straightforward inference of ancestry and admixture proportions through ancestry-informative insertion deletion multiplexing

Ancestry-informative markers (AIMs) show high allele frequency divergence between different ancestral or geographically distant populations. These genetic markers are especially useful in inferring the likely ancestral origin of an individual or estimating the apportionment of ancestry components in admixed individuals or populations. The study of AIMs is of great interest in clinical genetics research, particularly to detect and correct for population substructure effects in case-control association studies, but also in population and forensic genetics studies. This work presents a set of 46 ancestry-informative insertion deletion polymorphisms selected to efficiently measure population admixture proportions of four different origins (African, European, East Asian and Native American). All markers are analyzed in short fragments (under 230 basepairs) through a single PCR followed by capillary electrophoresis (CE) allowing a very simple one tube PCR-to-CE approach. HGDP-CEPH diversity panel samples from the four groups, together with Oceanians, were genotyped to evaluate the efficiency of the assay in clustering populations from different continental origins and to establish reference databases. In addition, other populations from diverse geographic origins were tested using the HGDP-CEPH samples as reference data. The results revealed that the AIM-INDEL set developed is highly efficient at inferring the ancestry of individuals and provides good estimates of ancestry proportions at the population level. In conclusion, we have optimized the multiplexed genotyping of 46 AIM-INDELs in a simple and informative assay, enabling a more straightforward alternative to the commonly available AIM-SNP typing methods dependent on complex, multi-step protocols or implementation of large-scale genotyping technologies.     

山东汉、回族男性人群父系遗传结构研究

本研究基于75个Y-SNP位点和23个Y-STR基因座对山东汉、回族男性人群进行研究,旨在揭示两个人群的父系遗传结构,为法医学应用及群体遗传学等提供基础数据。研究基于微测序技术检测187份山东汉族和130份山东回族样本,获取75个Y-SNP位点分型;采用PowerPlex^®Y23试剂盒检测23个Y-STR基因座;采用直接计数法统计等位基因频率、单倍型频率及单倍群频率,根据公式D=n(1-∑pi²)/(n-1)计算基因多样性、单倍型多样性以及单倍群多样性;根据Median-joining方法,使用NETWORK 5.0和NETWORK Publisher构建并展示网络图。研究结果显示,单倍群O-M175、C-M130、N-M231、Q-M242为山东汉族男性人群主要的Y单倍群,单倍群O-M175、J-M304、R-M207、C-M130、N-M231为山东回族男性人群最主要的单倍群;23个Y-STR基因座在山东汉族男性样本中检出187种单倍型,单倍型多样性为1.0000,在山东回族中检出121种单倍型,单倍型多样性为0.9988;网络图显示同一Y单倍群的样本相对独立地聚集在一起,山东汉族与回族人群之间存在共享单倍群,同时也存在一些特异性单倍群,如单倍群J-M304、R-M207均以山东回族为主,单倍群Q-M242则以山东汉族为主。山东汉族和回族男性人群的主要单倍群均为单倍群O-M175;单倍群J-M304、R-M207在山东回族中的高频分布,单倍群Q-M242则在山东汉族中高频分布。研究表明山东回族人群中保留有一定比例的欧亚西部和中东特有的Y染色体类型。     

Sexual selection and genetic colour polymorphisms in animals

Genetic colour polymorphisms are widespread across animals and often subjected to complex selection regimes. Traditionally, colour morphs were used as simple visual markers to measure allele frequency changes in nature, selection, population divergence and speciation. With advances in sequencing technology and analysis methods, several model systems are emerging where the molecular targets of selection are being described. Here, we discuss recent studies on the genetics of sexually selected colour polymorphisms, aiming at (i) reviewing the evidence of sexual selection on colour polymorphisms, (ii) highlighting the genetic architecture, molecular and developmental basis underlying phenotypic colour diversification and (iii) discuss how the maintenance of such polymorphisms might be facilitated or constrained by these. Studies of the genetic architecture of colour polymorphism point towards the importance of tight clustering of colour loci with other trait loci, such as in the case of inversions and supergene structures. Other interesting findings include linkage between colour loci and mate preferences or sex determination, and the role of introgression and regulatory variation in fuelling polymorphisms. We highlight that more studies are needed that explicitly integrate fitness consequences of sexual selection on colour with the underlying molecular targets of colour to gain insights into the evolutionary consequences of sexual selection on polymorphism maintenance.     

Molecular insight into the genesis of ranked caste populations of western India based upon polymorphisms across non-recombinant and recombinant regions in genome

Background

Large-scale trade and cultural contacts between coastal populations of western India and Western-Eurasians paved for extensive immigration and genesis of wide spectrum of admixed gene pool. To trace admixture and genesis of caste populations of western India, we have examined polymorphisms across non-recombining 20 Y-SNPs, 20 Y-STRs, 18 mtDNA diagnostic sites, HVS-1 plus HVS-2 regions; and recombining 15 highly polymorphic autosomal STRs in four predominant caste populations- upper-ranking Desasth-brahmin and Chitpavan-brahmin; a middle-ranking Kshtriya Maratha; and a lower-rank peasant Dhangar.

Results

The generated genomic data was compared with putative parental populations- Central Asians, West Asians and Europeans using AMOVA, PC plot, and admixture estimates. Overall, disparate uniparental ancestries, and l.1% GST value for biparental markers among four studied caste populations linked well with their exchequer demographic histories. Marathi-speaking ancient Desasth-brahmin shows substantial admixture from Central Asian males but Paleolithic maternal component support their Scytho-Dravidian origin. Chitpavanbrahmin demonstrates younger maternal component and substantial paternal gene flow from West Asia, thus giving credence to their recent Irano-Scythian ancestry from Mediterranean or Turkey, which correlated well with European-looking features of this caste. This also explains their untraceable ethno-history before 1000 years, brahminization event and later amalgamation by Maratha. The widespread Palaeolithic mtDNA haplogroups in Maratha and Dhangar highlight their shared Proto-Asian ancestries. Maratha males harboured Anatolianderived J2 lineage corroborating the blending of farming communities. Dhangar heterogeneity is ascribable to predominantly South-Asian males and West-Eurasian females.

Conclusions

The genomic data-sets of this study provide ample genomic evidences of diverse origins of four ranked castes and synchronization of caste stratification with asymmetrical gene flows from Indo-European migration during Upper Paleolithic, Neolithic, and later dates. However, subsequent gene flows among these castes living in geographical proximity, have diminished significant genetic differentiation as indicated by AMOVA and structure.     

The quest for Y-chromosomal markers - methodological strategies for mammalian non-model organisms

Tracing maternal and paternal lineages independently to explore breeding systems and dispersal strategies in natural populations has been high on the wish-list of evolutionary biologists. As males are the heterogametic sex in mammals, such sex-specific patterns can be indirectly observed when Y chromosome polymorphism is combined with mitochondrial sequence information. Over the past decade, Y-chromosomal markers applied to human populations have revealed remarkable differences in the demographic history and behaviour between the sexes. However, with a few exceptions, genetic data tracing the paternal line are lacking in most other mammalian species. This deficit can be attributed to the difficulty of developing Y-specific genetic markers in non-model organisms and the general low levels of polymorphisms observed on the Y chromosome. Here, we present an overview of the currently employed strategies for developing paternal markers in mammals. Moreover, we review the practical feasibility and requirements of various methodological strategies and highlight their future prospects when combined with new molecular techniques such as next generation sequencing.     

Disclosing the genetic structure of Brazil through analysis of male lineages with highly discriminating haplotypes

In a large variety of genetic studies, probabilistic inferences are made based on information available in population databases. The accuracy of the estimates based on population samples are highly dependent on the number of chromosomes being analyzed as well as the correct representation of the reference population. For frequency calculations the size of a database is especially critical for haploid markers, and for countries with complex admixture histories it is important to assess possible substructure effects that can influence the coverage of the database. Aiming to establish a representative Brazilian population database for haplotypes based on 23 Y chromosome STRs, more than 2,500 Y chromosomes belonging to Brazilian, European and African populations were analyzed. No matter the differences in the colonization history of the five geopolitical regions that currently exist in Brazil, for the Y chromosome haplotypes of the 23 studied Y-STRs, a lack of genetic heterogeneity was found, together with a predominance of European male lineages in all regions of the country. Therefore, if we do not consider the diverse Native American or Afro-descendent isolates, which are spread through the country, a single Y chromosome haplotype frequency database will adequately represent the urban populations in Brazil. In comparison to the most commonly studied group of 17 Y-STRs, the 23 markers included in this work allowed a high discrimination capacity between haplotypes from non-related individuals within a population and also increased the capacity to discriminate between paternal relatives. Nevertheless, the expected haplotype mutation rate is still not enough to distinguish the Y chromosome profiles of paternally related individuals. Indeed, even for rapidly mutating Y-STRs, a very large number of markers will be necessary to differentiate male lineages from paternal relatives.     

Insights into Iberian population origins through the construction of highly informative Y-chromosome haplotypes using biallelic markers,STRs, and the MSY1 minisatellite

To investigate the diversity of Y chromosomes in the Iberian Peninsula and the North African population of Maghreb, we constructed superhaplotypes on the basis of 10 biallelic markers, 7 microsatellites, and 1 minisatellite located in the nonrecombining portion of the human Y chromosome. The analysis of extremely high MSY1 variability was performed by reducing the MVR-codes to modular structures. Y-STRs and MSY1 data provide information about the relationship between closely related populations such as those of Iberia. Analysis of biallelic markers allowed us to identify 7 of 12 haplogroups defined by those polymorphisms. The haplogroup background showed clear differences between Iberian populations and the North African one. The use of differently mutating Y-chromosome markers allowed us to infer different population events at different time scales: the Paleolithic background of the Iberian Peninsula, the Neolithic fingerprint on Y-chromosome lineages, and the Iron Age influence in the populations of Iberia. Implications of our results for the highly debated origin of Basques are also discussed.     

Automation of Molecular-Based Analyses: A Primer on Massively Parallel Sequencing

Recent advances in genetics have been enabled by new genetic sequencing techniques called massively parallel sequencing (MPS) or next-generation sequencing. Through the ability to sequence in parallel hundreds of thousands to millions of DNA fragments, the cost and time required for sequencing has dramatically decreased. There are a number of different MPS platforms currently available and being used in Australia. Although they differ in the underlying technology involved, their overall processes are very similar: DNA fragmentation, adaptor ligation, immobilisation, amplification, sequencing reaction and data analysis. MPS is being used in research, translational and increasingly now also in clinical settings. Common applications include sequencing of whole genomes, whole exomes or targeted genes for disease-causing gene discovery, genetic diagnosis and targeted cancer therapy. Even though the revolution that is occurring with MPS is exciting due to its increasing use, improving and emerging technologies and new applications, significant challenges still exist. Particularly challenging issues are the bioinformatics required for data analysis, interpretation of results and the ethical dilemma of ‘incidental findings’.     

Signature of recent historical events in the European Y-chromosomal STR haplotype distribution

Previous studies of human Y-chromosomal single-nucleotide polymorphisms (Y-SNPs) established a link between the extant Y-SNP haplogroup distribution and the prehistoric demography of Europe. By contrast, our analysis of seven rapidly evolving Y-chromosomal short tandem repeat loci (Y-STRs) in over 12,700 samples from 91 different locations in Europe reveals a signature of more recent historic events, not previously detected by other genetic markers. Cluster analysis based upon molecular variance yields two clearly identifiable sub-clusters of Western and Eastern European Y-STR haplotypes, and a diverse transition zone in central Europe, where haplotype spectra change more rapidly with longitude than with latitude. This and other observed patterns of Y-STR similarity may plausibly be related to particular historical incidents, including, for example, the expansion of the Franconian and Ottoman Empires. We conclude that Y-STRs may be capable of resolving male genealogies to an unparalleled degree and could therefore provide a useful means to study local population structure and recent demographic history.L. Roewer and P.J.P. Croucher contributed equally to this paper.This revised version was published online in February 2005 with corrections to the addresses of the authors of the Forensic Y-Chromosome Research Group.     

Massively Parallel Sequencing of Patients with Intellectual Disability,Congenital Anomalies and/or Autism Spectrum Disorders with a Targeted Gene Panel

Developmental delay and/or intellectual disability (DD/ID) affects 1–3% of all children. At least half of these are thought to have a genetic etiology. Recent studies have shown that massively parallel sequencing (MPS) using a targeted gene panel is particularly suited for diagnostic testing for genetically heterogeneous conditions. We report on our experiences with using massively parallel sequencing of a targeted gene panel of 355 genes for investigating the genetic etiology of eight patients with a wide range of phenotypes including DD/ID, congenital anomalies and/or autism spectrum disorder. Targeted sequence enrichment was performed using the Agilent SureSelect Target Enrichment Kit and sequenced on the Illumina HiSeq2000 using paired-end reads. For all eight patients, 81–84% of the targeted regions achieved read depths of at least 20×, with average read depths overlapping targets ranging from 322× to 798×. Causative variants were successfully identified in two of the eight patients: a nonsense mutation in the ATRX gene and a canonical splice site mutation in the L1CAM gene. In a third patient, a canonical splice site variant in the USP9X gene could likely explain all or some of her clinical phenotypes. These results confirm the value of targeted MPS for investigating DD/ID in children for diagnostic purposes. However, targeted gene MPS was less likely to provide a genetic diagnosis for children whose phenotype includes autism.     

Genetic admixture, relatedness, and structure patterns among Mexican populations revealed by the Y-chromosome

Y-linked markers are suitable loci to analyze genetic diversity of human populations, offering knowledge of medical, forensic, and anthropological interest. In a population sample of 206 Mestizo males from western Mexico, we analyzed two binary loci (M3 and YAP) and six Y-STRs, adding to the analysis data of Mexican Mestizos and Amerindians, and relevant worldwide populations. The paternal ancestry estimated in western Mexican-Mestizos was mainly European (60-64%), followed by Amerindian (25-21%), and African ( approximately 15%). Significant genetic heterogeneity was established between Mestizos from western (Jalisco State) and northern Mexico (Chihuahua State) compared with Mexicans from the center of the Mexican Republic (Mexico City), this attributable to higher European ancestry in western and northern than in central and southeast populations, where higher Amerindian ancestry was inferred. This genetic structure has important implications for medical and forensic purposes. Two different Pre-Hispanic evolutionary processes were evident. In Mesoamerican region, populations presented higher migration rate (N(m) = 24.76), promoting genetic homogeneity. Conversely, isolated groups from the mountains and canyons of the Western and Northern Sierra Madre (Huichols and Tarahumaras, respectively) presented a lower migration rate (N(m) = 10.27) and stronger genetic differentiation processes (founder effect and/or genetic drift), constituting a Pre-Hispanic population substructure. Additionally, Tarahumaras presented a higher frequency of Y-chromosomes without Q3 that was explained by paternal European admixture (15%) and, more interestingly, by a distinctive Native-American ancestry. In Purepechas, a special admixture process involving preferential integration of non-Purepecha women in their communities could explain contrary genetic evidences (autosomal vs. Y-chromosome) for this tribe.     

Spiked GBS: a unified,open platform for single marker genotyping and whole-genome profiling

Background

In plant breeding, there are two primary applications for DNA markers in selection: 1) selection of known genes using a single marker assay (marker-assisted selection; MAS); and 2) whole-genome profiling and prediction (genomic selection; GS). Typically, marker platforms have addressed only one of these objectives.

Results

We have developed spiked genotyping-by-sequencing (sGBS), which combines targeted amplicon sequencing with reduced representation genotyping-by-sequencing. To minimize the cost of targeted assays, we utilize a small percent of sequencing capacity available in runs of GBS libraries to “spike” amplified targets of a priori alleles tagged with a different set of unique barcodes. This open platform allows multiple, single-target loci to be assayed while simultaneously generating a whole-genome profile. This dual-genotyping approach allows different sets of samples to be evaluated for single markers or whole genome-profiling. Here, we report the application of sGBS on a winter wheat panel that was screened for converted KASP markers and newly-designed markers targeting known polymorphisms in the leaf rust resistance gene Lr34.

Conclusions

The flexibility and low-cost of sGBS will enable a range of applications across genetics research. Specifically in breeding applications, the sGBS approach will allow breeders to obtain a whole-genome profile of important individuals while simultaneously targeting specific genes for a range of selection strategies across the breeding program.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1404-9) contains supplementary material, which is available to authorized users.     

Will genetics really revolutionize the drug discovery process?

Genetics has played only a modest role in drug discovery, but new technologies will radically change this. Whole genome sequencing will identify new drug discovery targets, and emerging methods for the determination of gene function will increase the ability to select robust targets. Detection of single nucleotide polymorphisms and common polymorphisms will enhance the investigation of polygenic diseases and the use of genetics in drug development. Oligonucleotide arraying technologies will allow analysis of gene expression patterns in novel ways.     

Discovery of single-nucleotide polymorphisms (SNPs) in the uncharacterized genome of the ascomycete Ophiognomonia clavigignenti-juglandacearum from 454 sequence data

The benefits from recent improvement in sequencing technologies, such as the Roche GS FLX (454) pyrosequencing, may be even more valuable in non-model organisms, such as many plant pathogenic fungi of economic importance. One application of this new sequencing technology is the rapid generation of genomic information to identify putative single-nucleotide polymorphisms (SNPs) to be used for population genetic, evolutionary, and phylogeographic studies on non-model organisms. The focus of this research was to sequence, assemble, discover and validate SNPs in a fungal genome using 454 pyrosequencing when no reference sequence is available. Genomic DNA from eight isolates of Ophiognomonia clavigignenti-juglandacearum was pooled in one region of a four-region sequencing run on a Roche 454 GS FLX. This yielded 71 million total bases comprising 217,000 reads, 80% of which collapsed into 16,125,754 bases in 30,339 contigs upon assembly. By aligning reads from multiple isolates, we detected 298 SNPs using Roche's GS Mapper. With no reference sequence available, however, it was difficult to distinguish true polymorphisms from sequencing error. Eagleview software was used to manually examine each contig that contained one or more putative SNPs, enabling us to discard all but 45 of the original 298 putative SNPs. Of those 45 SNPs, 13 were validated using standard Sanger sequencing. This research provides a valuable genetic resource for research into the genus Ophiognomonia, demonstrates a framework for the rapid and cost-effective discovery of SNP markers in non-model organisms and should prove especially useful in the case of asexual or clonal fungi with limited genetic variability.     

Detection of DNA sequence polymorphisms by enzymatic amplification and direct genomic sequencing.

The discovery of RFLPs and their utilization as genetic markers has revolutionized research in human molecular genetics. However, only a fraction of the DNA sequence polymorphisms in the human genome affect the length of a restriction fragment and hence result in an RFLP. Polymorphisms that are not detected as RFLPs are typically passed over in the screening process though they represent a potentially important source of informative genetic markers. We have used a rapid method for the detection of naturally occurring DNA sequence variations that is based on enzymatic amplification and direct sequencing of genomic DNA. This approach can detect essentially all useful sequence variations within the region screened. We demonstrate the feasibility of the technique by applying it to the human retinoblastoma susceptibility locus. We screened 3,712 bp of genomic DNA from each of nine individuals and found four DNA sequence polymorphisms. At least one of these DNA sequence polymorphisms was informative in each of three families with hereditary retinoblastoma that were not informative with any of the known RFLPs at this locus. We believe that direct sequencing is a reasonable alternative to other methods of screening for DNA sequence polymorphisms and that it represents a step forward for obtaining informative markers at well-characterized loci that have been minimally informative in the past.     

A Hu sheep genome with the first ovine Y chromosome reveal introgression history after sheep domestication

The Y chromosome plays key roles in male fertility and reflects the evolutionary history of paternal lineages. Here, we present a de novo genome assembly of the Hu sheep with the first draft assembly of ovine Y chromosome(o MSY), using nanopore sequencing and Hi-C technologies. The o MSY that we generated spans 10.6 Mb from which 775 Y-SNPs were identified by applying a large panel of whole genome sequences from worldwide sheep and wild Iranian mouflons. Three major paternal lineages(HY1a, HY1b and HY2) were defined across domestic sheep, of which HY2 was newly detected. Surprisingly, HY2 forms a monophyletic clade with the Iranian mouflons and is highly divergent from both HY1a and HY1b. Demographic analysis of Y chromosomes, mitochondrial and nuclear genomes confirmed that HY2 and the maternal counterpart of lineage C represented a distinct wild mouflon population in Iran that diverge from the direct ancestor of domestic sheep, the wild mouflons in Southeastern Anatolia. Our results suggest that wild Iranian mouflons had introgressed into domestic sheep and thereby introduced this Iranian mouflon specific lineage carrying HY2 to both East Asian and Africa sheep populations.     

Ancestral informative marker selection and population structure visualization using sparse Laplacian eigenfunctions

Identification of a small panel of population structure informative markers can reduce genotyping cost and is useful in various applications, such as ancestry inference in association mapping, forensics and evolutionary theory in population genetics. Traditional methods to ascertain ancestral informative markers usually require the prior knowledge of individual ancestry and have difficulty for admixed populations. Recently Principal Components Analysis (PCA) has been employed with success to select SNPs which are highly correlated with top significant principal components (PCs) without use of individual ancestral information. The approach is also applicable to admixed populations. Here we propose a novel approach based on our recent result on summarizing population structure by graph laplacian eigenfunctions, which differs from PCA in that it is geometric and robust to outliers. Our approach also takes advantage of the priori sparseness of informative markers in the genome. Through simulation of a ring population and the real global population sample HGDP of 650K SNPs genotyped in 940 unrelated individuals, we validate the proposed algorithm at selecting most informative markers, a small fraction of which can recover the similar underlying population structure efficiently. Employing a standard Support Vector Machine (SVM) to predict individuals' continental memberships on HGDP dataset of seven continents, we demonstrate that the selected SNPs by our method are more informative but less redundant than those selected by PCA. Our algorithm is a promising tool in genome-wide association studies and population genetics, facilitating the selection of structure informative markers, efficient detection of population substructure and ancestral inference.     

牦牛分子遗传多样性研究进展

遗传多样性研究可有效地揭示牦牛的遗传变异, 是牦牛群体遗传学研究的主要内容之一。自20世纪70年代以来, 人们已对牦牛的体形外貌特征、染色体核型(带型)、生理生化特性和DNA序列变异等进行了较为深入地研究。随着分子遗传学和DNA测序技术的迅猛发展, 近年来的研究主要集中在牦牛的分子遗传多样性。文章对近15年来牦牛mtDNA和核基因组分子标记及侯选基因多样性的研究现状进行了综述, 对前景进行展望, 以期为牦牛群体基因组学等研究提供依据。     

Next generation deep sequencing and vaccine design: today and tomorrow

Next generation sequencing (NGS) technologies have redefined the modus operandi in both human and microbial genetics research, allowing the unprecedented generation of very large sequencing datasets on a short time scale and at affordable costs. Vaccine development research is rapidly taking full advantage of the advent of NGS. This review provides a concise summary of the current applications of NGS in relation to research seeking to develop vaccines for human infectious diseases, incorporating studies of both the pathogen and the host. We focus on rapidly mutating viral pathogens, which are major targets in current vaccine research. NGS is unraveling the complex dynamics of viral evolution and host responses against these viruses, thus contributing substantially to the likelihood of successful vaccine development.     

Skin pigmentation,biogeographical ancestry and admixture mapping

Ancestry informative markers (AIMs) are genetic loci showing alleles with large frequency differences between populations. AIMs can be used to estimate biogeographical ancestry at the level of the population, subgroup (e.g. cases and controls) and individual. Ancestry estimates at both the subgroup and individual level can be directly instructive regarding the genetics of the phenotypes that differ qualitatively or in frequency between populations. These estimates can provide a compelling foundation for the use of admixture mapping (AM) methods to identify the genes underlying these traits. We present details of a panel of 34 AIMs and demonstrate how such studies can proceed, by using skin pigmentation as a model phenotype. We have genotyped these markers in two population samples with primarily African ancestry, viz. African Americans from Washington D.C. and an African Caribbean sample from Britain, and in a sample of European Americans from Pennsylvania. In the two African population samples, we observed significant correlations between estimates of individual ancestry and skin pigmentation as measured by reflectometry (R(2)=0.21, P<0.0001 for the African-American sample and R(2)=0.16, P<0.0001 for the British African-Caribbean sample). These correlations confirm the validity of the ancestry estimates and also indicate the high level of population structure related to admixture, a level that characterizes these populations and that is detectable by using other tests to identify genetic structure. We have also applied two methods of admixture mapping to test for the effects of three candidate genes (TYR, OCA2, MC1R) on pigmentation. We show that TYR and OCA2 have measurable effects on skin pigmentation differences between the west African and west European parental populations. This work indicates that it is possible to estimate the individual ancestry of a person based on DNA analysis with a reasonable number of well-defined genetic markers. The implications and applications of ancestry estimates in biomedical research are discussed.